Mitogenomics: digging deeper with complete mitochondrial genomes

Characterization of the complete mitochondrial genomes of two species of the genus Aphaena Guérin-Méneville (Hemiptera: Fulgoridae) and its phylogenetic implications

The complete mitochondrial genomes (mitogenomes) of Aphaena (Callidepsa) amabilis and Aphaena (Aphaena) discolor nigrotibiata were sequenced. The mitogenomes of these two species are 16,237 bp and 16,116 bp in length with an A + T content of 77.9% and 77.0%, respectively. Each contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (A + T-rich region). All PCGs initiate with the standard start codon of ATN and terminate with the complete stop codon of TAA or TAG except for atp6, where nad1 ends with an incomplete T codon. All tRNAs have the typical clover-leaf structure except for trnS1 and trnV which have a reduced DHU arm. Moreover, these two mitogenomes have trnL2, trnR and trnT with an unpaired base in the acceptor stem. The putative A + T-rich region includes multiple types of tandem repeat regions. These phylogenetic analyses are reconstructed based on 13 protein-coding genes of 25 auchenorrhynchan mitogenomes, with both maximum likelihood and Bayesian analyses yielding robust identical phylogenetic trees. These results support a monophyletic Auchenorrhyncha and the relationship (Pyrops + (Lycorma + Aphaena)) within Fulgoridae.

Complete mitochondrial genomes confirm the generic placement of the plateau vole, Neodon fuscus

The plateau vole, Neodon fuscus is endemic to China and is distributed mainly in Qinghai Province. It is of public health interest, as it is, a potential reservoir of Toxoplasma gondii and the intermediate host of Echinococcus multilocularis. However, genetic data of this species are lacking, and its name and taxonomy are still a controversy. In the present study, we determined the nucleotide sequence of the entire mitochondrial (mt) genome of N. fuscus and analyzed its evolutionary relationship. The mitogenome was 16328 bp in length and contained 13 protein-coding genes, 22 genes for transfer RNAs (tRNA), two ribosomal RNA genes and two major noncoding regions (O_L region and D-loop region). Most genes were located on the heavy strand. All tRNA genes had typical cloverleaf structures except for tRNA^{Ser (GCU)}. The mt genome of N. fuscus was rich in A+T (58.45%). Maximum likelihood (ML) and Bayesian methods yielded phylogenetic trees from 33 mt genomes of Arvicolinae, in which N. fuscus formed a sister group with Neodon irene and Neodon sikimensis to the exclusion of species of Microtus and other members of the Arvicolinae. Further phylogenetic analyses (ML only) based on the cytb gene sequences also demonstrated that N. fuscus had a close relationship with N. irene. The complete mitochondrial genome was successfully assembled and annotated, providing the necessary information for the phylogenetic analyses. Although the name Lasiopodomys fuscus was used in the book ‘Wilson & Reeder’s Mammal Species of the World’, we have confirmed here that its appropriate name is N. fuscus through an analysis of the evolutionary relationships.

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.

Mitochondrial phylogeny and comparative mitogenomics of closely related pine moth pests (Lepidoptera: Dendrolimus)

BACKGROUND

Pine moths, Dendrolimus spp. (Lasiocampidae), are serious economic pests of conifer forests. Six closely related species (Dendrolimus punctatus, D. tabulaeformis, D. spectabilis, D. superans, D. houi, and D. kikuchii) occur in China and cause serious damage to coniferophyte. The complete mito genomes of Dendrolimus genus are significant to resolve the phylogenetic relationship and provide theoretical support in pest control.

METHODS

The complete mitogenomes of three species (D. superans, D. houi, and D. kikuchii) were sequenced based on PCR-amplified with universal primers, which were used to amplify initial fragments. Phylogenetic analyses were carried out with 78 complete mitogenomes of lepidopteran species from 10 superfamilies.

RESULTS

The complete mitochondrial genomes of these three species were 15,417, 15,381, and 15,377 bp in length, separately. The phylogenetic analyses produced consistent results for six Dendrolimus species based on complete mitogenomes, two major clades were formed, one containing D. spectabilis clustered with D. punctatus + D. tabulaeformis, and D. superans as the sister group to this three-taxon clade, the other containing D. kikuchii and D. houi. Comparative analyses of the congeneric mitochondrial genomes were performed, which showed that non-coding regions were more variable than the A+T rich region. The mitochondrial nucleotide diversity was more variable when compared within than among genus, and the concatenated tRNA region was the most conserved and the nd6 genes was the most variable.

The complete mitochondrial genome of Aquila nipalensis and its phylogenetic position

Mitochondrial genome sequences are valuable resources for systematics and conservation biology studies. In this paper, we present the complete mitogenome of Aquila nipalensis which was 18,450 bp in length. The gene content and arrangement were typical for avian mtDNA. The overall A + T content of was 54.1%, and the AT skew was calculated as 0.12 for the complete mitogenome of A. nipalensis. The maximum-likelihood (ML) tree based on the concatenated 12 protein-coding genes (PCGs) revealed the basal phylogenetic position of A. nipalensis in Aquila.

The complete mitochondrial genome and phylogenetic analysis of Spondylus violaceus

Spondylus violaceus is a species of the bivalve family Spondylidae, commonly known as thorny oyster. Here, we report the complete mitochondrial genome of S. violaceus. The mitochondrial genome of S. violaceus is 30,160 bp in length, which contains 12 protein-coding genes, 2 ribosomal RNA genes and 19 transfer RNA genes. The overall AT content (58.31%) was higher than GC content (41.69%). This is the first record of complete mitochondrial genome in the family Spondylidae, providing the basis for future taxonomic and phylogenetic studies of this bivalve family.

The first mitochondrial genome of the family Epicopeiidae and higher-level phylogeny of Macroheterocera (Lepidoptera: Ditrysia)

The Macroheterocera clade contains most of the Lepidoptera species. However, extensive comparative and phylogenetic analyses of this group using complete mitochondrial genomes (mitogenome) are limited particularly in the context of increasing macroheteroceran mitogenomes reported to date. In this study, complete mitogenome of the Epicopeia hainesii is determined as the first Epicopeiidae species with mitogenome available. The whole mitogenome is circular with 15,395 bp long, and is highly biased toward A + T nucleotides (80.6%) in nucleotide composition. Comparative analyses show that gene content and arrangement of macroheteroceran mitogenomes are generally conservative and are typical of Lepidoptera but exceptions exist. In newly sequenced mitogenome, the motif "ATACTAA" is putatively located at the end of gene nad1, rather than in intergenic sequences between trnS2 and nad1 genes routinely observed in Lepidoptera. Interestingly, multiple phylogenetic analyses recover the six macroheteroceran superfamilies as Mimallonoidea + (Drepanoidea + ((Bombycoidea + Lasiocampoidea) + (Noctuoidea + Geometroidea))), providing supports for a large-scale transcriptomic study rather than various mitogenome- and multiple-gene-based investigations. In addition, our analyses consistently place the Epicopeiidae as sister group with Geometroidea, firstly demonstrating that this family is closer with Geometroidea than Drepanoidea based on mitogenome data.

Comparative analysis of the mitochondrial genomes of Orthonectida: insights into the evolution of an invertebrate parasite species

Among invertebrates, only a few groups still have uncertain phylogenetic position, Orthonectida, a small group of rare multi-cellular parasites of marine invertebrates, being one of them. Recent molecular and morphological findings suggest that orthonectids belong to Lophotrochozoa and are close to Annelida. Nevertheless, phylogenetic relationships between orthonectids and annelids are unclear, and the phylogeny within the group itself has never been studied. Sequencing of mitochondrial genomes is used here to clarify this issue. Complete mt genomes of the orthonectids Intoshia variabili and Rhopalura litoralis were characterized and compared with Intoshia linei mt genome. Our results show that Orthonectida mt genomes have undergone reduction and gene loss, and that they have complicated organization revealed in strand asymmetry in nucleotide composition, in some features of intergenic non-coding regions, tRNA duplication and folding. Moreover, all species of Orthonectida have a unique gene order with complicated rearrangement landscape. Significant differences in mitochondrial genomes in the three orthonectid species could be explained by the fact that their host species belong to different taxa (flat worms, nemertines and gastropods). Among the analyzed mt genomes of Orthonectida, I. linei possesses the closest gene order to the ancestral genome. All Orthonectida species are monophyletic, and in the phylogenetic tree are close to Pleistoannelida, and specifically, to Clitellata.

The first complete mitochondrial genome of marigold pest thrips, Neohydatothrips samayunkur (Sericothripinae) and comparative analysis

Complete mitogenomes from the order Thysanoptera are limited to representatives of the subfamily Thripinae. Therefore, in the present study, we sequenced the mitochondrial genome of Neohydatothrips samayunkur (15,295 bp), a member of subfamily Sericothripinae. The genome possesses the canonical 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs) as well as two putative control regions (CRs). The majority strand was 77.42% A + T content, and 22.58% G + C with weakly positive AT skew (0.04) and negative GC skew (-0.03). The majority of PCGs start with ATN codons as observed in other insect mitochondrial genomes. The GCG codon (Alanine) was not used in N. samayunkur. Most tRNAs have the typical cloverleaf secondary structure, however the DHU stem and loop were absent in trnV and trnS1, while the TΨC loop was absent in trnR and trnT. The two putative control regions (CR1 and CR2) show 99% sequence similarity indicated a possible duplication, and shared 57 bp repeats were identified. N. samayunkur showed extensive gene rearrangements, with 11 PCGs, 22 tRNAs, and two rRNAs translocated when compared to the ancestral insect. The gene trnL2 was separated from the 'trnL2-cox2' gene block, which is a conserved, ancestral gene order found in all previously sequenced thrips mitogenomes. Both maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees resulted in similar topologies. The phylogenetic position of N. samayunkur indicates that subfamily Sericothripinae is sister to subfamily Thripinae. More molecular data from different taxonomic groups is needed to understand thrips phylogeny and evolution.

Comparative analysis of the mitochondrial genomes of oriental spittlebug trible Cosmoscartini: insights into the relationships among closely related taxa

Background

Cosmoscartini (Hemiptera: Cercopoidea: Cercopidae) is a large and brightly colored Old World tropical tribe, currently containing over 310 phytophagous species (including some economically important pests of eucalyptus in China) in approximately 17 genera. However, very limited information of Cosmoscartini is available except for some scattered taxonomic studies. Even less is known about its phylogenetic relationship, especially among closely related genera or species. In this study, the detailed comparative genomic and phylogenetic analyses were performed on nine newly sequenced mitochondrial genomes (mitogenomes) of Cosmoscartini, with the purpose of exploring the taxonomic status of the previously defined genus Okiscarta and some closely related species within the genus Cosmoscarta.

Results

Mitogenomes of Cosmoscartini display similar genomic characters in terms of gene arrangement, nucleotide composition, codon usage and overlapping regions. However, there are also many differences in intergenic spacers, mismatches of tRNAs, and the control region. Additionally, the secondary structures of rRNAs within Cercopidae are inferred for the first time.

Based on comparative genomic (especially for the substitution pattern of tRNA secondary structure) and phylogenetic analyses, the representative species of Okiscarta uchidae possesses similar structures with other Cosmoscarta species and is placed consistently in Cosmoscarta. Although Cosmoscarta bimacula is difficult to be distinguished from Cosmoscarta bispecularis by traditional morphological methods, evidence from mitogenomes highly support the relationships of (C. bimacula + Cosmoscarta rubroscutellata) + (C. bispecularis + Cosmoscarta sp.).

Conclusions

This study presents mitogenomes of nine Cosmoscartini species and represents the first detailed comparative genomic and phylogenetic analyses within Cercopidae. It is indicated that knowledge of mitogenomes can be effectively used to resolve phylogenetic relationships at low taxonomic levels. Sequencing more mitogenomes at various taxonomic levels will also improve our understanding of mitogenomic evolution and phylogeny in Cercopidae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5365-7) contains supplementary material, which is available to authorized users.

Comparative analysis of seven mitochondrial genomes of Phymatostetha (Hemiptera: Cercopidae) and phylogenetic implications

In this study, we present seven mitochondrial genomes (mitogenomes) of Phymatostetha. Each mitogenome contains the entire set of 37 genes, which arranged in the same order as the putative ancestral pattern of insects. The nucleotide composition of Phymatostetha mitogenomes is biased toward A/T, with rRNAs and PCG12 (i.e. the first and second codon positions of PCGs) exhibit the highest and lowest A + T content, respectively. Relative synonymous codon usage of PCGs also show that degenerate codons are biased to use more A/T than G/C. All tRNAs exhibit typical clover-leaf structure, with the exception of trnS1. Additionally, unpaired nucleotides are detected in trnS1 anticodon stem and trnR acceptor stem. Phylogenetic relationships, based on the dataset of 13 PCGs, 22tRNAs, and two rRNAs, are analyzed using both the Bayesian and maximum likelihood methods. Our results clearly revealed the systematic status of Phymatostetha species and robustly supported the monophyly of this genus, in which Phymatostetha semele is sister to other Phymatostetha species. It was demonstrated that mitogenome was an effective molecular marker to adequately resolve phylogeny at low taxonomic levels.

The complete mitochondrial genome of a medicinal insect, Hydrillodes repugnalis (Lepidoptera: Noctuoidea: Erebidae), and related phylogenetic analysis

The complete mitochondrial genome (mitogenome) of an important medicinal insect Hydrillodes repugnalis (Lepidoptera: Noctuoidea) was sequenced and analyzed. The mitogenome is circular with 15,570 bp long, and shows typical gene content and arrangement. Nucleotide composition is highly biased toward A + T nucleotides (81.1%). All protein-coding genes (PCGs) initiate with canonical start codon ATN, except for cox1 being CGA. The typical stop codon TAA is used for most PCGs, while the nad4l uses the TAG, and cox1 and cox2 use incomplete termination codon T. All tRNAs have a typical clover-leaf structure, except for trnS1 (AGN) lacking the dihydrouridine arm. Comparative mitogenome analysis showed that the motif "ATGATAA" between atp8 and atp6, and the motif "ATACTAA" between trnS2 and nad1 are commonly present in noctuoid mitogenomes. In A + T-rich regions, the motif "ATAGA" and subsequent poly-T structure, the motif "ATTTA" and followed macrosatellite (AT)_n element and an "A"-rich 3' end upstream of the trnM gene can be recognized across noctuoid mitogenomes. Phylogenetic analyses showed that H. repugnalis is nested into the Erebidae clade, consistently being sister to the Aganainae. Within Noctuoidea, the (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))) was consistently recovered firstly based on multiple mitochondrial datasets.

Mitochondrial genome of Argopecten irradians reveals higher-level phylogenetic relationships in Anisomyaria

The complete mitochondrial genome of Argopecten irradians strain Zhongkehong was sequenced and annotated: it is 16,212 bp in length and contains twelve protein-coding genes (atp8 is absent, as in most species in Anisomyaria), two ribosomal RNA genes, and 21 transfer RNA genes (trnS is absent and there are two copies of trnF). The heavy strand has an overall A + T content of 57.3%; GC and AT skews are 0.249 and -0.262, respectively, indicating more Gs and more Ts than Cs and As. Phylogenetic analysis based on Bayesian Inference and Maximum Likelihood of the twelve protein-coding genes shows that A. irradians has close relationships with A. purpuratus and A. ventricosus; this indicated that A. irradians belongs to the Pectinidae family. The Pectinidae was sister to (Ostreidae + Mytilidae). This work provides general information on the evolution of cultured scallops.

Comparative mitochondrial genome analysis of Grammodes geometrica and other noctuid insects reveals conserved mitochondrial genome organization and phylogeny

The mitochondrial genome (mitogenome) plays an important role in revealing molecular evolution. In this study, the complete mitogenome of Grammodes geometrica (G. geometrica) (Lepidoptera: Erebidae) was sequenced and characterized. The nucleotide composition of the genome is highly A + T biased, accounting for 80.49%. Most protein-coding genes (PCGs) are initiated by ATN codons except for the cytochrome oxidase subunit 1 (cox1) gene, which was initiated by CGA. The order and orientation of genes with the order trnM-trnI-trnQ-nad2 is a typical rearrangement compared with those ancestral insects in which trnM is located between trnQ and nad2. Most tRNA genes were folded into the typical cloverleaf structure except for trnS1 (AGN). The A + T-rich region contains the conserved motif "ATAGA" followed by a 19 bp poly-T stretch, which was also observed in other Noctuoidea species. In addition, we reconstructed phylogenetic trees among the nucleotide alignments of five families of Noctuoidea species except the Oenosandridae. Finally, we achieved a well-supported tree, which showed that G. geometrica belongs to the Erebidae family. Moreover, the relationships at the family-level can be displayed as follows: (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))).

Complete mitochondrial genome of the first deep-sea spongicolid shrimp Spongiocaris panglao (Decapoda: Stenopodidea): Novel gene arrangement and the phylogenetic position and origin of Stenopodidea

Stenopodidea Claus, 1872 (Crustacea: Decapoda) is one of the major groups of decapods crustaceans. Hitherto, only one complete mitochondrial genome (mitogenome) from the family Stenopodidae is available for the infraorder Stenopodidea. Here, we determined the complete mitogenome of Spongiocaris panglao de Grave and Saito, 2016 using Illumina sequencing, representing the first species from the family Spongicolidae. The 15,909 bp genome is a circular molecule and consists of 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. Although the overall genome organization is typical for metazoans, the mitogenome of S. panglao shows some derived characters. A + T content of 77.42% in S. pamglao mitogenome is second-highest among the dacapods described to date. The trnR gene exhibit modified secondary structure with the TψC loop completely missing, which might be a putative autapomorphy of S. pamglao mitogenome. Compared with the shallow-water stenopodidean species S. hispidus, the control region of S. pamglao exhibits three characteristics: larger size, higher A + T content, and more tandem repeat sequences. The gene order exhibited difference from the ancestral mitogenome pattern of the Pancrustacea, with 5 tRNA genes rearrangement. The result from BI was agreed with most morphological characters and molecular evidences, revealing that Stenopodidea and Reptantia had the closest relationship, as the sister group of Caridea. Still, the alternative hypothesis supported from ML topology cannot be completely rejected based on the current data. Estimated times revealed that the two stenopodideans families Stenopodidae and Spongicolidae diverged from each other around 122 Mya. The divergence time of spongicolid shrimp is in good agreement with the origin of their hexactinellid hosts (78-144 Mya).

Information from the mitochondrial genomes of two egg parasitoids, Gonatocerus sp. and Telenomus sp., reveals a controversial phylogenetic relationship between Mymaridae and Scelionidae

The taxonomic status and phylogenetic affinities of Mymaridae and Scelionidae are controversial, based on similarities between these families in the characteristics of adults, larvae, and eggs. In this study, we sequenced the mitochondrial (mt) genomes of representatives from these two families and found that the derived secondary structure of tRNA-Arg was the same in each family due to the absence of the D-stem. The segment of "cox1 trnL₂cox2 trnK trnD atp8 atp6 cox3" in Gonatocerus sp. (Mymaridae) is conserved and distinct from those of four other species of Chalcidoidea but similar to that in Proctotrupoidea and Platygastroidea. However, phylogenetic analysis indicated that Gonatocerus sp. was sister group to other species of Chalcidoidea. Comparisons based on complete gene orders may be more useful in a phylogenetic and systematic context, as different branches may exhibit partially homoplastic gene orders.

Characterization of the Complete Mitochondrial Genome Sequences of Three Croakers (Perciformes, Sciaenidae) and Novel Insights into the Phylogenetics

The three croakers (Nibea coibor, Protonibea diacanthus and Argyrosomus amoyensis, Perciformes, Sciaenidae) are important commercial species inhabiting the Eastern Indian Ocean and Western Pacific. Molecular data employed in previous research on phylogenetic reconstruction have not been adequate and complete, and systematic and comprehensive phylogenetic relationships for these fish are unresolved. We sequenced the complete mitochondrial genomes of the three croakers using next-generation sequencing for the first time. We analyzed the composition and phylogenies between 19 species in the family Sciaenidae using the mitochondrial protein coding sequences of 204 species in the Series Eupercaria. We present the characterization of the complete mitochondrial genome sequences of the three croakers. Gene arrangement and distribution of the three croakers are canonically identical and consistent with other vertebrates. We found that the family Sciaenidae is an independent branch that is isolated from the order Perciformes and does not belong to any extant classification. Therefore, this family is expected to belong to a new classification at the order level and needs further analysis. The evolution of Sciaenidae has lagged far behind the Perciformes differentiation. This study presents a novel insight into the phylogenetics of the family Sciaenidae from the order Perciformes and facilitates additional studies on the evolution and phylogeny of Series Eupercaria.

Complete mitochondrial genome of the soybean leaffolder, Omiodes indicata (Lepidoptera: Pyraloidea: Crambidae), and phylogenetic analysis for Pyraloidea

Complete mitochondrial genome (mitogenome) of the Omiodes indicata was sequenced and characterized. The circular mitogenome is 15,367bp long, including 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an A+T-rich region. Nucleotide composition is highly biased toward A+T nucleotides (81.6%). All 13 PCGs initiate with canonical start codon (ATN), except for cox1 that initiates with CGA. All tRNAs have a typical clover-leaf structure, except for trnS1 (AGN) in which the base pairs of the dihydrouridine (DHU) arm are reduced. In O. indicata, the motifs "ATGATAA" and "ATACTAA" between atp8 and atp6, trnS2 and nad1, respectively, and the motifs "ATAG" and "ATTTA" in the A+T-rich region can be identified. Comparative phylogenetic analyses based on four datasets show that the dataset including all coding positions of 13 PCGs exhibit the highest informativeness in resolving higher phylogeny of Pyraloidea. Bayesian inference (BI) and maximum likelihood (ML) analyses yield generally well-supported phylogenetic relationships among the eleven pyraloid subfamilies involved. However, the relationships among the five subfamilies (Acentropinae, Crambinae, Glaphyriinae, Schoenobiinae and Scopariinae) in ML analysis are ambiguous, which might be resolved by ample sampling in future mitogenome-based phylogenetic studies.

Characterization of the complete mitogenome of Indian Mouse Deer, Moschiola indica (Artiodactyla: Tragulidae) and its evolutionary significance

The mitochondrial genome of Indian mouse deer (Moschiola indica) was sequenced, assembled and characterized for the first time using 22 pairs of polymerase chain reaction (PCR) primers. The mitogenome of M. indica which is 16,444 bp in size was found very similar to most vertebrates in organisation that harbours 13 protein-coding genes, 22 transfer RNA, 2 ribosomal RNA and 1A + T-rich region. Its comparison with over 52 mitogenomes of the order Artiodactyla, showed a conserved nature of gene organisation, codon usage, gene orientation and evolutionary rates of proteins except that M. indica possesses an extra copy of trnF. The complete mitogenome and protein-coding genes of M. indica were found to be highly A + T biased. Rate of protein evolution was highest in atp8 and lowest in cox3. Further, a higher purifying selection pressure was found to be acting on family Tragulidae compared to Bovidae and Cervidae. The phylogenetic analysis of M. indica placed the Tragulidae as sister-group of all other ruminants, similar to previous analyses. Moschiola forms the sister-group to the other two tragulid genera Tragulus (from Asia) and Hyemoschus (from Africa), which is unexpected as usually the Asian species are thought to form a monophyletic group.

The complete mitochondrial genome of Damora sagana and phylogenetic analyses of the family Nymphalidae

The monotypic genus Damora (Nymphalidae, Heliconiinae) contains a single species, Damora sagana, which is widely distributed across southern China. Herein, its complete mitogenome was sequenced to further understand lepidopteran mitogenome characteristics, reconstruct the nymphalid family phylogeny, and infer the subdivision of Heliconiinae species. The circular mitogenome was 15,151 bp long, abundant in A and T, and comprised of 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and one control region with a gene arrangement typical of lepidopteran mitogenomes. ATN codons initiated all PCGs, except cytochrome c oxidase subunit 1 (COX1), which was initiated by a CGA sequence as has been observed in other lepidopterans. Three PCGs (COX1, COX2 and ND4) employed a single T termination signal, whereas others had the typical complete termination codon (TAA). All tRNA genes were folded into the typical cloverleaf structure except for tRNA-Ser (AGN). The A+T-rich region included the conserved motif 'ATAGA' followed by a 17 bp poly-T stretch, which was also observed in tribe Argynnini mitogenomes. A phylogenetic tree was constructed via multiple methods using the 13 PCGs data of D. sagana and other available mitogenomes of nymphalid species. All three phylogenetic trees yielded the same topology. These results were consistent with those from previous studies of most major nymphalid groups, except those regarding tribe subdivision in certain subfamilies such as Argynnini + (Acraeini + Heliconiini) for Heliconiine. Furthermore, our analyses identified that the genus Cethosia was grouped with the genus Acraea composing the tribe Acraeini with strong support.

Deep-level phylogeny of Cicadomorpha inferred from mitochondrial genomes sequenced by NGS

Recent development and advancement of next-generation sequencing (NGS) technologies have enabled the determination of mitochondrial genome (mitogenome) at extremely efficiency. In this study, complete or partial mitogenomes for 19 cicadomorphan species and six fulgoroid species were reconstructed by using the method of high-throughput sequencing from pooled DNA samples. Annotation analyses showed that the mitogenomes obtained have the typical insect mitogenomic content and structure. Combined with the existing hemipteran mitogenomes, a series of datasets with all 37 mitochondrial genes (up to 14,381 nt total) under different coding schemes were compiled to test previous hypotheses of deep-level phylogeny of Cicadomorpha. Thirty-seven species representing Cicadomorpha constituted the ingroup. A taxon sampling with nine species from Fulgoroidea and six from Heteroptera comprised the outgroup. The phylogenetic reconstructions congruently recovered the monophyly of each superfamily within Cicadomorpha. Furthermore, the hypothesis (Membracoidea + (Cicadoidea + Cercopoidea)) was strongly supported under the heterogeneous CAT model.

Comparative mitochondrial genome analysis of Spilarctia subcarnea and other noctuid insects

This study was performed to better understand the phylogenetic relationships within the lepidopteran superfamily Noctuoidea. The mitochondrial genome (mitogenome) has been extensively used for studying phylogenetic relationships at different taxonomic levels. In this study, the complete mitogenome of Spilarctia subcarnea (Noctuoidea: Erebidae) was sequenced and annotated. The mitogenome is 15,441bp in length, containing 13 typical protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a noncoding control region (CR). The order and orientation of genes of S. subcarnea mitogenome with the order trnM-trnI-trnQ-nad2 is different from the ancestral insects in which trnM is located between trnQ and nad2 (trnI-trnQ-trnM-nad2). The phylogenetic relationships based on mitochondrial sequences using Bayesian inference and Maximum likelihood methods showed that S. subcarnea was closely related to Lemyra melli, supporting that S. subcarnea belongs to Erebidae. These analyses confirm that Lymantriidae should be included as subfamilies within Erebidae. The Erebidae was sister to (Nolidae+(Euteliidae+Noctuidae)); Notodontidae is sister to the other families of Noctuoidea in our study.

Characterization of the complete mitochondrial genome of Acanthoscelides obtectus (Coleoptera: Chrysomelidae: Bruchinae) with phylogenetic analysis

Acanthoscelides obtectus is a common species of the subfamily Bruchinae and a worldwide-distributed seed-feeding beetle. The complete mitochondrial genome of A. obtectus is 16,130 bp in length with an A + T content of 76.4%. It contains a positive AT skew and a negative GC skew. The mitogenome of A. obtectus contains 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes and a non-coding region (D-loop). All PCGs start with an ATN codon, and seven (ND3, ATP6, COIII, ND3, ND4L, ND6, and Cytb) of them terminate with TAA, while the remaining five (COI, COII, ND1, ND4, and ND5) terminate with a single T, ATP8 terminates with TGA. Except tRNA ^Ser , the secondary structures of 21 tRNAs that can be folded into a typical clover-leaf structure were identified. The secondary structures of lrRNA and srRNA were also predicted in this study. There are six domains with 48 helices in lrRNA and three domains with 32 helices in srRNA. The control region of A. obtectus is 1354 bp in size with the highest A + T content (83.5%) in a mitochondrial gene. Thirteen PCGs in 19 species have been used to infer their phylogenetic relationships. Our results show that A. obtectus belongs to the family Chrysomelidae (subfamily-Bruchinae). This is the first study on phylogenetic analyses involving the mitochondrial genes of A. obtectus and could provide basic data for future studies of mitochondrial genome diversities and the evolution of related insect lineages.

Complete mitochondrial genome of Clistocoeloma sinensis (Brachyura: Grapsoidea): Gene rearrangements and higher-level phylogeny of the Brachyura

Deciphering the animal mitochondrial genome (mitogenome) is very important to understand their molecular evolution and phylogenetic relationships. In this study, the complete mitogenome of Clistocoeloma sinensis was determined. The mitogenome of C. sinensis was 15,706 bp long, and its A+T content was 75.7%. The A+T skew of the mitogenome of C. sinensis was slightly negative (−0.020). All the transfer RNA genes had the typical cloverleaf structure, except for the trnS1 gene, which lacked a dihydroxyuridine arm. The two ribosomal RNA genes had 80.2% A+T content. The A+T-rich region spanned 684 bp. The gene order within the complete mitogenome of C. sinensis was identical to the pancrustacean ground pattern except for the translocation of trnH. Additionally, the gene order of trnI-trnQ-trnM in the pancrustacean ground pattern becomes trnQ-trnI-trnM in C. sinensis. Our phylogenetic analysis showed that C. sinensis and Sesarmops sinensis cluster together with high nodal support values, indicating that C. sinensis and S. sinensis have a sister group relationship. The results support that C. sinensis belongs to Grapsoidea, Sesarmidae. Our findings also indicate that Varunidae and Sesarmidae species share close relationships. Thus, mitogenomes are likely to be valuable tools for systematics in other groups of Crustacea.

The complete mitochondrial genome sequence of the Sichuan Digging Frog, Kaloula rugifera (Anura: Microhylidae) and its phylogenetic implications

The Sichuan Digging Frog (Kaloula rugifera) belongs to the family Dicroglossidae, which is endemic to northeastern Sichuan and southernmost Gansu provinces, in southwestern China. In this study, the complete mitochondrial genome of K. rugifera was sequenced. The mitogenome was 17,074bp in length, consisting of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a non-coding control region. As in other vertebrates, most mitochondrial genes are encoded on the heavy strand, except for ND6 and eight tRNA genes which are encoded on the light strand. The overall base composition of the K. rugifera is 30.32% A, 25.76% C, 29.72% T, and 14.20% G, which is consistent with the lowest frequency for G content in typical amphibian animals' mitochondrial genomes. The alignment of the Kaloula species control regions exhibited high genetic variability and rich A+T content. Besides, 3 types of tandem repeat units were also identified in the control region. Phylogenetic tree demonstrated that K. rugifera was clustered together with K. borealis and K. verrucosa and they had a close relationship with each other. The complete mitogenome of K. rugifera can provide an important data for the studies on phylogenetic relationship to further explore the taxonomic status of Kaloula species.

Genetic diversity and population structure of Cirrhinus mrigala revealed by mitochondrial ATPase 6 gene

The genetic diversity and population structure of Cirrhinus mrigala were studied using partial mitochondrial ATPase 6 gene. 155 individuals were collected from six distant rivers of peninsular India. A total of 11 haplotypes were found with haplotype diversity value of 0.46711 and nucleotide diversity 0.00154 indicating a low degree of genetic divergence among populations. The result of AMOVA analysis also revealed that 97.53% of the genetic variation contained within populations and 2.47% occurred among populations. An overall fixation index was found to be 0.02466 (p < .05). Pairwise Fst analysis showed significant genetic differentiation between the Kaveri and rest of the populations. The minimum spanning tree of the haplotype data demonstrated a common ancestral origin of the mrigal populations in peninsular rivers. From the overall results of the present study it may be inferred that mrigal population in peninsular India have low genetic differentiation. This information would be very essential for setting up of genetic improvement programmes, effective management, and conservation of natural genetic resources of mrigal.

Mitochondrial genomes of the hoverflies Episyrphus balteatus and Eupeodes corollae (Diptera: Syrphidae), with a phylogenetic analysis of Muscomorpha

The hoverflies Episyrphus balteatus and Eupeodes corollae (Diptera: Muscomorpha: Syrphidae) are important natural aphid predators. We obtained mitochondrial genome sequences from these two species using methods of PCR amplification and sequencing. The complete Episyrphus mitochondrial genome is 16,175 bp long while the incomplete one of Eupeodes is 15,326 bp long. All 37 typical mitochondrial genes are present in both species and arranged in ancestral positions and directions. The two mitochondrial genomes showed a biased A/T usage versus G/C. The cox1, cox2, cox3, cob and nad1 showed relatively low level of nucleotide diversity among protein-coding genes, while the trnM was the most conserved one without any nucleotide variation in stem regions within Muscomorpha. Phylogenetic relationships among the major lineages of Muscomorpha were reconstructed using a complete set of mitochondrial genes. Bayesian and maximum likelihood analyses generated congruent topologies. Our results supported the monophyly of five species within the Syrphidae (Syrphoidea). The Platypezoidea was sister to all other species of Muscomorpha in our phylogeny. Our study demonstrated the power of the complete mitochondrial gene set for phylogenetic analysis in Muscomorpha.

Genetic diversity of Ixodes ricinus (Ixodida: Ixodidae) ticks in sympatric and allopatric zones in Baltic countries

Ixodes ricinus (Linnaeus 1758) and Ixodes persulcatus (Schulze 1930) ticks are involved in the transmission of a wide variety of pathogens with considerable impact on human and animal health. The co-distribution zone of these two tick species is situated in the Baltic countries, which provides a special setting for the population studies. In the present study, genetic variability of I. ricinus ticks collected in allopatric and sympatric locations in the Baltic countries has been investigated using a sequence analysis of the mitochondrial DNA control region, 16S rRNA and cytb genes. There were 32 haplotypes (Hd: 0.8551) and 27 haplotypes (Hd:0.8213) of control region sequences from ticks in allopatric and sympatric zones detected, respectively. Out of 47 16S rRNA gene haplotypes, 32 haplotypes (Hd: 0.7213) were found in the allopatric zone and 27 (Hd:0.9572) in the sympatric zone. The Cytb gene was very conserved and monomorphic in ticks from the allopatric zone, whereas three unique haplotypes were observed in the sympatric zone. The higher number of unique haplotypes of the control region was detected in the allopatric zone. Median joining network and Fst analysis did not reveal a clear separation between ticks from the two zones.

Characterization of complete mitochondrial genome of fives tripe wrasse (Thalassoma quinquevittatum, Lay & Bennett, 1839) and phylogenetic analysis

To further supplement the genome-level features in related species, T. quinquevittatum complete mtDNA was firstly sequenced and de novo assembled by next-generation sequencing. The full-length mtDNA of T. quinquevittatum was a 16,896bp fragment, which was atypical of Labridae, with 2 ribosomal RNA (rRNA) genes, 13 protein-coding genes (PCGs), 23 transfer RNA (tRNA) genes, and a major non-coding control region (D-loop region). Additionally, the mtDNA of T. quinquevittatum exhibited characteristics of A (27.1%), T (29.3%), G (17.8%), and C (25.8%) with a high A+T content (56.4%). Furthermore, the analysis of the average Ka/Ks in the 13 PCGs of three Labridae species indicated a strong purifying selection within this group. Additionally, the phylogenetic analysis based on 13 concatenated PCGs nucleotide and amino acid datasets, showed high value support for the following sister clade among the four genera (T. quinquevittatum, Halichoeres trimaculatus, Halichoeres melanurus, Parajulis poecilepterus). The complete mtDNA of the T. quinquevittatum provided important information for the study in population genetics and evolutionary theory.

Next-Generation Sequencing of Two Mitochondrial Genomes from Family Pompilidae (Hymenoptera: Vespoidea) Reveal Novel Patterns of Gene Arrangement

Animal mitochondrial genomes have provided large and diverse datasets for evolutionary studies. Here, the first two representative mitochondrial genomes from the family Pompilidae (Hymenoptera: Vespoidea) were determined using next-generation sequencing. The sequenced region of these two mitochondrial genomes from the species Auplopus sp. and Agenioideus sp. was 16,746 bp long with an A + T content of 83.12% and 16,596 bp long with an A + T content of 78.64%, respectively. In both species, all of the 37 typical mitochondrial genes were determined. The secondary structure of tRNA genes and rRNA genes were predicted and compared with those of other insects. Atypical trnS1 using abnormal anticodons TCT and lacking D-stem pairings was identified. There were 49 helices belonging to six domains in rrnL and 30 helices belonging to three domains in rrns present. Compared with the ancestral organization, four and two tRNA genes were rearranged in mitochondrial genomes of Auplopus and Agenioideus, respectively. In both species, trnM was shuffled upstream of the trnI-trnQ-trnM cluster, and trnA was translocated from the cluster trnA-trnR-trnN-trnS1-trnE-trnF to the region between nad1 and trnL1, which is novel to the Vespoidea. In Auplopus, the tRNA cluster trnW-trnC-trnY was shuffled to trnW-trnY-trnC. Phylogenetic analysis within Vespoidea revealed that Pompilidae and Mutillidae formed a sister lineage, and then sistered Formicidae. The genomes presented in this study have enriched the knowledge base of molecular markers, which is valuable in respect to studies about the gene rearrangement mechanism, genomic evolutionary processes and phylogeny of Hymenoptera.

Application of RNA-seq for mitogenome reconstruction, and reconsideration of long-branch artifacts in Hemiptera phylogeny

Hemiptera make up the largest nonholometabolan insect assemblage. Despite previous efforts to elucidate phylogeny within this group, relationships among the major sub-lineages remain uncertain. In particular, mitochondrial genome (mitogenome) data are still sparse for many important hemipteran insect groups. Recent mitogenomic analyses of Hemiptera have usually included no more than 50 species, with conflicting hypotheses presented. Here, we determined the nearly complete nucleotide sequence of the mitogenome for the aphid species of Rhopalosiphum padi using RNA-seq plus gap filling. The 15,205 bp mitogenome included all mitochondrial genes except for trnF. The mitogenome organization and size for R. padi are similar to previously reported aphid species. In addition, the phylogenetic relationships for Hemiptera were examined using a mitogenomic dataset which included sequences from 103 ingroup species and 19 outgroup species. Our results showed that the seven species representing the Aleyrodidae exhibit extremely long branches, and always cluster with long-branched outgroups. This lead to the failure of recovering a monophyletic Hemiptera in most analyses. The data treatment of Degen-coding for protein-coding genes and the site-heterogeneous CAT model show improved suppression of the long-branch effect. Under these conditions, the Sternorrhyncha was often recovered as the most basal clade in Hemiptera.

Multiple Lines of Evidence from Mitochondrial Genomes Resolve Phylogenetic Relationships of Parasitic Wasps in Braconidae

The rapid increase in the number of mitochondrial genomes in public databases provides opportunities for insect phylogenetic studies; but it also provides challenges because of gene rearrangements and variable substitution rates among both lineages and sites. Typically, phylogenetic studies use mitochondrial sequence data but exclude other features of the mitochondrial genome from analyses. Here, we undertook large-scale sequencing of mitochondrial genomes from a worldwide collection of specimens belonging to Braconidae, one of the largest families of Metazoa. The strand-asymmetry of base composition in the mitochondrial genomes of braconids is reversed, providing evidence for monophyly of the Braconidae. We have reconstructed a backbone phylogeny of the major lineages of Braconidae from gene order of the mitochondrial genomes. Standard phylogenetic analyses of DNA sequences provided strong support for both Cyclostomes and Noncyclostomes. Four subfamily complexes, that is, helconoid, euphoroid, sigalphoid, and microgastroid, within the Noncyclostomes were reconstructed robustly, the first three of which formed a monophyletic group sister to the last one. Aphidiinae was recovered as a lineage sister to other groups of Cyclostomes, while the Ichneutinae was recovered as paraphyletic. Separate analyses of the subdivided groups showed congruent relationships, employing different matrices and methods, for the internal nodes of the Cyclostomes and the microgastroid complex of subfamilies. This research, using multiple lines of evidence from mitochondrial genomes, illustrates multiple uses of mitochondrial genomes for phylogenetic inference in Braconidae.

Four new mitochondrial genomes of the genus zosterops (aves: passeriformes: zosteropidae) from East Africa with a phylogenetic evaluation of the group

The white-eye birds of the genus Zosterops have been recognized for their high speciation rates in the past, but the relationships of the East African populations are not yet fully resolved. We sequenced and annotated mitogenomes of four populations currently assigned to three East African white-eye species, Zosterops senegalensis, Z. abyssinicus and Z. poliogaster. For Z. senegalensis specimens from two distant populations were sequenced; for the other taxa we used samples collected at one site. The mitogenomes ranged between 17,827 and 17,974 bp, in size similar to previously published mitogenomes analyzed for this genus from other geographic regions. The mitogenomes contain the classical set of 13 coding genes, two structural rRNA genes and 22 tRNA genes. We constructed a phylogeny using all complete mitogenomes currently available for the genus. The phylogeny supports an Asian or Oceanic origin of the genus Zosterops. The East African species represent a monophyletic clade, but the two specimens of Zosterops senegalensis from different regions do not group together, supporting previous hypotheses of cryptic species within the genus. The new genetic resources provided here may help to further explore the relationships and evolution of the genus.

Merging scleractinian genera: the overwhelming genetic similarity between solitary Desmophyllum and colonial Lophelia

Background

In recent years, several types of molecular markers and new microscale skeletal characters have shown potential as powerful tools for phylogenetic reconstructions and higher-level taxonomy of scleractinian corals. Nonetheless, discrimination of closely related taxa is still highly controversial in scleractinian coral research. Here we used newly sequenced complete mitochondrial genomes and 30 microsatellites to define the genetic divergence between two closely related azooxanthellate taxa of the family Caryophylliidae: solitary Desmophyllum dianthus and colonial Lophelia pertusa.

Results

In the mitochondrial control region, an astonishing 99.8 % of nucleotides between L. pertusa and D. dianthus were identical. Variability of the mitochondrial genomes of the two species is represented by only 12 non-synonymous out of 19 total nucleotide substitutions. Microsatellite sequence (37 loci) analysis of L. pertusa and D. dianthus showed genetic similarity is about 97 %. Our results also indicated that L. pertusa and D. dianthus show high skeletal plasticity in corallum shape and similarity in skeletal ontogeny, micromorphological (septal and wall granulations) and microstructural characters (arrangement of rapid accretion deposits, thickening deposits).

Conclusions

Molecularly and morphologically, the solitary Desmophyllum and the dendroid Lophelia appear to be significantly more similar to each other than other unambiguous coral genera analysed to date. This consequently leads to ascribe both taxa under the generic name Desmophyllum (priority by date of publication). Findings of this study demonstrate that coloniality may not be a robust taxonomic character in scleractinian corals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0654-8) contains supplementary material, which is available to authorized users.

Complete mitochondrial genomes of two gelechioids, Mesophleps albilinella and Dichomeris ustalella (Lepidoptera: Gelechiidae), with a description of gene rearrangement in Lepidoptera

We sequenced the entire mitochondrial genome (mitogenome) of two gelechioids, Mesophleps albilinella and Dichomeris ustalella, and compared their genome organization and sequence composition to those of available gelechioid mitogenomes for an enhanced understanding of Gelechioidea genomic characteristics. We compared all available lepidopteran mitogenome arrangements, including that of M. albilinella, which is unique in Gelechioidea, to comprehend the extensiveness and mechanisms of gene rearrangement in Lepidoptera. The genomes of M. albilinella and D. ustalella are 15,274 and 15,410 bp in size, respectively, with the typical sets of mitochondrial (mt) genes. The COI gene begins with CGA (arginine) in all sequenced gelechioids, including M. albilinella and D. ustalella, reinforcing the feature as a synapomorphic trait, at least in the Gelechioidea. Each 353- and 321-bp long A + T-rich region of M. albilinella and D. ustalella contains one (D. ustalella) or two (M. albilinella) tRNA-like structures. The M. albilinella mitogenome has a unique gene arrangement among the Gelechioidea: ARNESF (the underline signifies an inverted gene) at the ND3 and ND5 junction, as opposed to the ARNSEF that is found in ancestral insects. An extensive search of available lepidopteran mitogenomes, including that of M. albilinella, turned up six rearrangements that differ from those of ancestral insects. Most of the rearrangements can be explained by the tandem duplication-random loss model, but inversion, which requires recombination, is also found in two cases, including M. albilinella. Excluding the MIQ rearrangement at the A + T-rich region and ND2 junction, which is found in nearly all Ditrysia, most of the remaining rearrangements found in Lepidoptera appear to be independently derived in that they are automorphic at several taxonomic scales, although current mitogenomic data are limited, particularly for congeneric data.