Fourteen complete mitochondrial genomes of butterflies from the genus Lethe (Lepidoptera,Nymphalidae, Satyrinae) with mitogenome-based phylogenetic analysis

The mitochondrial genome (mitogenome) can help us understand the phylogenetic relationships within the genus Lethe and the subfamily Satyrinae. In this study, we sequenced the complete mitogenomes of 14 Lethe species, which range in size from 15,225 to 15,271 bp, with both 37 genes (13 PCGs, 22 tRNAs, 2 rRNAs) and a noncoding A + T-rich region. The gene arrangement and orientation is similar to typical mitogenomes of Lepidoptera. The Ka/Ks ratio shows that cox1 has the slowest evolutionary rate. The secondary structure of trnN lacks the Pseudouracil loop (TψC loop) in most Lethe species. The inferred phylogenetic analyses show that Lethe is a well-supported monophyletic group, and reveal 2 major clades within the genus Lethe, which is consistent with previous morphological classifications.     

The Complete Mitochondrial Genome of two Tetragnatha Spiders (Araneae: Tetragnathidae): Severe Truncation of tRNAs and Novel Gene Rearrangements in Araneae

Mitogenomes can provide information for phylogenetic analysis and evolutionary biology. The Araneae is one of the largest orders of Arachnida with great economic importance. In order to develop mitogenome data for this significant group, we determined the complete mitogenomes of two long jawed spiders Tetragnatha maxillosa and T. nitens and performed the comparative analysis with previously published spider mitogenomes. The circular mitogenomes are 14578 bp long with A+T content of 74.5% in T. maxillosa and 14639 bp long with A+T content of 74.3% in T. nitens, respectively. Both the mitogenomes contain a standard set of 37 genes and an A+T-rich region with the same gene orientation as the other spider mitogenomes, with the exception of the different gene order by the rearrangement of two tRNAs (trnW and trnG). Most of the tRNAs lose TΨC arm stems and have unpaired amino acid acceptor arms. As interesting features, both trnS^AGN and trnS^UCN lack the dihydrouracil (DHU) arm and long tandem repeat units are presented in the A+T-rich region of both the spider mitogenomes. The phylogenetic relationships of 23 spider mitogenomes based on the concatenated nucleotides sequences of 13 protein-coding genes indicated that the mitogenome sequences could be useful in resolving higher-level relationship of Araneae. The molecular information acquired from the results of this study should be very useful for future researches on mitogenomic evolution and genetic diversities in spiders.     

Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing

Little is known about the variations of nematode mitogenomes (mtDNA). Sequencing a complete mtDNA using a PCR approach remains a challenge due to frequent genome reorganizations and low sequence similarities between divergent nematode lineages. Here, a genome skimming approach based on HiSeq sequencing (shotgun) was used to assemble de novo the first complete mtDNA sequence of a root-knot nematode (Meloidogyne graminicola). An AT-rich genome (84.3%) of 20,030 bp was obtained with a mean sequencing depth superior to 300. Thirty-six genes were identified with a semi-automated approach. A comparison with a gene map of the M. javanica mitochondrial genome indicates that the gene order is conserved within this nematode lineage. However, deep genome rearrangements were observed when comparing with other species of the superfamily Hoplolaimoidea. Repeat elements of 111 bp and 94 bp were found in a long non-coding region of 7.5 kb, as similarly reported in M. javanica and M. hapla. This study points out the power of next generation sequencing to produce complete mitochondrial genomes, even without a reference sequence, and possibly opening new avenues for species/race identification, phylogenetics and population genetics of nematodes.     

Two novel mitogenomes of Dipodidae species and phylogeny of Rodentia inferred from the complete mitogenomes

Two novel mitogenomes of Eozapus setchuanus (KJ648495) and Sicista concolor (KJ648496) were reported and their total lengths were 16,630 bp and 16,493 bp, respectively. Both mitogenomes which were analogous to other rodent mitogenomes, contained 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a control region. Specifically, the ND2 gene of S. concolor has three amino acids lesser than that of two other Dipodidae species (E. setchuanus and Jaculus jaculus) due to a premature termination codon in the 3′ end. We detected a tandem repeat cluster of 221 bp and 274 bp in the control region of S. concolor and E. setchuanus, respectively. Along with phylogenetic relationship analysis, we speculated that the tandem repeats in control regions might be common in Dipodinae species. Our phylogenetic analysis using concatenated mitochondrial gene datasets suggested five suborder and 16 family monophyletic groups in 54 rodent taxa sampled and strongly supported a basal position of the squirrel-related clade (PP = 1; BP = 100). Dipodidae had a sister-group relationship with Muroidea, and Sicistinae was in the base of Dipodidae clade. The complete mitochondrial genomes showed high resolution in deep-level phylogenetic relationship reconstructions of Rodentia.     

Complete mitochondrial genomes of Conogethes punctiferalis and C. pinicolalis (Lepidoptera: Crambidae): Genomic comparison and phylogenetic inference in Pyraloidea

Conogethes punctiferalis Guenée, 1854 (Lepidoptera: Crambidae) is a serious polyphagous pest that attacks more than one hundred species of plants. Previously, C. punctiferalis was determined to be composed of two ecotypes; later, Conogethes pinicolalis, was described as a separate species. Due to the prolonged negligence of C. pinicolalis as an independent species, the genetic perspective of the two species is limited. Thus, in this study, 15,332 and 15,336 bp-long complete mitochondrial genome (mitogenome) of the two species were sequenced and compared to each other and to 54 available mitogenomes of Pyraloidea. The comparison of each protein-coding gene (PCG) and rRNA gene of the two congeneric species showed substantial sequence divergence, ranging from 3.13% (ATP8) to 8.3% (COIII), with an average of 5.92%. Phylogenetic analyses using concatenated sequences of 13 PCGs and 2 rRNAs (12,458 bp including gaps), both by maximum likelihood (ML) and Bayesian inference (BI) methods, consistently supported the monophyly of each family (Crambidae and Pyralidae) and subfamily, generally with the highest nodal supports. The subfamilial relationships of ((((Acentropinae + Schoenobiinae) + (Scopariinae + Crambinae)) + Evergestinae) + (Spilomelinae + Pyraustinae)) in Crambidae and the subfamilial relationships of ((((Pyralinae + Epipaschiinae) + Phycitinae) + Galleriilinae)) in Pyralidae were obtained in both analyses. However, nodal supports were substantially low in this study, mainly due to limited taxa.     

Characterization of the complete mitochondrial genome of Phodopus roborovskii (Rodentia: Cricetidae) and systematic implications for Cricetinae phylogenetics

Phodopus roborovskii (subfamily Cricetinae) is widely distributed in the northern arid regions of China. This study reports its complete mitochondrial genome (mitogenome) for the first time. The complete sequence was 16,273 bp long, including 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and 1 major noncoding region. The base composition and codon usage were described. The putative origin of replication for the light strand (O_L) of P. roborovskii was approximately 45 bp long and was highly conserved in the stem-loop and adjacent sequences, but the starting sequence of replication varied between genera among Rodentia. We analyzed the three domains of the D-loop region, and the results indicated that the central domain had higher G + C content and lower A + T content than two peripheral domains. Phylogenetic analyses indicated high resolution in four main divergent clades using mitogenomes data within Cricetidae. Within Cricetinae clade, P. roborovskii was at basal position which was in line with previous researches, and it shared a common ancestor with other extant hamsters. This work validated previous molecular and karyotype researches using mitogenomes data, and provided a set of useful data on phylogeny and molecular evolution in Cricetidae species.     

Comparative analysis of mitochondrial genomes of Nirvanini and Evacanthini (Hemiptera: Cicadellidae) reveals an explicit evolutionary relationship

Mitogenomes of five leafhopper species, Chudania hellerina and Concaveplana rufolineata in Nirvanini, Carinata rufipenna, Evacanthus danmainus and E. heimianus representing Evacanthini, were sequenced. The lengths of these five mitogenomes range from 15,044 (C. hellerina) to 15,680 bp (E. heimianus). All five mitogenomes exhibit similar base composition, gene size and codon usage of protein-coding genes. All 22 tRNA genes have typical cloverleaf secondary structures, except for trnS1 (AGN) which appears to lack the dihydrouridine arm. The two included Nirvanini species employ the anticodon TCT instead of the commonly used GCT in trnS1 (AGN). Genes nad2, atp8 and nad6 were highly variable while cox1 and cob showed the lowest nucleotide diversity. Phylogenetic analyses of two concatenated nucleotide datasets, incorporating the newly sequenced taxa and other available membracoid mitogenomes, recovered each included leafhopper subfamily as monophyletic with evacanthine tribes Nirvanini and Evacanthini forming monophyletic sister clades. A relationship among Evacanthinae, Cicadellinae and Typhlocybinae received moderate branch support.     

Complete mitochondrial genomes of three skippers in the tribe Aeromachini (Lepidoptera: Hesperiidae: Hesperiinae) and their phylogenetic implications

The mitochondrial genome is now widely used in the study of phylogenetics and molecular evolution due to its maternal inheritance, fast evolutionary rate, and highly conserved gene content. To explore the phylogenetic relationships of the tribe Aeromachini within the subfamily Hesperiinae at the mitochondrial genomic level, we sequenced and annotated the complete mitogenomes of 3 skippers: Ampittia virgata, Halpe nephele, and Onryza maga (new mitogenomes for 2 genera) with a total length of 15,333 bp, 15,291 bp, and 15,381 bp, respectively. The mitogenomes all contain 13 protein‐coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a noncoding A + T‐rich region and are consistent with other lepidopterans in gene order and type. In addition, we reconstructed the phylogenetic trees of Hesperiinae using maximum likelihood (ML) and Bayesian inference (BI) methods based on mitogenomic data. Results show that the tribe Aeromachini in this study robustly constitute a monophyletic group in the subfamily Hesperiinae, with the relationships Coeliadinae + (Euschemoninae + (Pyrginae + ((Eudaminae + Tagiadinae) + (Heteropterinae + ((Trapezitinae + Barcinae) + Hesperiinae))))). Moreover, our study supports the view that Apostictopterus fuliginosus and Barca bicolor should be placed out of the subfamily Hesperiinae.     

Evolutionary relatedness of mackerels of the genus Scomber based on complete mitochondrial genomes: Strong support to the recognition of Atlantic Scomber colias and Pacific Scomber japonicus as distinct species

Mackerels of the genus Scomber are commercially important species, but their taxonomic status is still controversial. Although previous phylogenetic data support the recognition of Atlantic Scomber colias and Pacific Scomber japonicus as separate species, it is only based on the analysis of partial mitochondrial and nuclear DNA sequences. In an attempt to shed light on this relevant issue, we have determined the complete mitochondrial DNA sequence of S. colias, S. japonicus, and Scomber australasicus. The total length of the mitogenomes was 16,568 bp for S. colias and 16,570 bp for both S. japonicus and S. australasicus. All mitogenomes had a gene content (13 protein-coding, 2 rRNAs, and 22 tRNAs) and organization similar to that observed in Scomber scombrus and most other vertebrates. The major noncoding region (control region) ranged between 865 and 866 bp in length and showed the typical conserved blocks. Phylogenetic analyses revealed a monophyletic origin of Scomber species with regard to other scombrid fish. The major finding of this study is that S. colias and S. japonicus were significantly grouped in distinct lineages within Scomber cluster, which phylogenetically constitutes evidence that they may be considered as separate species. Additionally, molecular data here presented provide a useful tool for evolutionary as well as population genetic studies.     

Complete mitochondrial genomes of three Oxya grasshoppers (Orthoptera) and their implications for phylogenetic reconstruction

Oxya is a genus of grasshoppers (Orthoptera: Acridoidea) attacking rice and other gramineous plants in Africa and Asia. In the present study, we characterized complete mitochondrial genomes (mitogenomes) of three species, Oxya japonica japonica (15,427 bp), Oxya hainanensis (15,443 bp) and Oxya agavisa robusta (15,552 bp) collected from China. The three mitogenomes contained a typical gene set of metazoan mitogenomes and shared the same gene order with other Acridid grasshoppers, including the rearrangement of tRNA^Asp and tRNA^Lys. Analyses of pairwise genetic distances showed that ATP8 was the least conserved gene, while COI the most conserved. To determine the position of Oxya grasshoppers in the phylogeny of Acrididae, we reconstructed phylogenetic trees among 64 species from across 11 subfamilies using nucleotide sequences of mitogenomes. While the tree confirms traditional classifications of Acrididae at major higher-levels, it suggests a few modifications for classifications at lower-levels.     

Evolutionary dynamics of chloroplast genomes in subfamily Aroideae (Araceae)

Aroideae is the largest and most diverse subfamily of the plant family Araceae. Despite its agricultural and horticultural importance, the genomic resources are sparse for this subfamily. Here, we report de novo assembled and fully annotated chloroplast genomes of 13 Aroideae species. The quadripartite chloroplast genomes (size range of 158,177–170,037 bp) are comprised of a large single copy (LSC; 75,594–94,702 bp), a small single copy (SSC; 12,903–23,981 bp) and a pair of inverted repeats (IRs; 25,266–34,840 bp). Notable gene rearrangements and IRs contraction / expansions were found for Anchomanes hookeri and Zantedeschia aethiopica. Codon usage, amino acid frequencies, oligonucleotide repeats, GC contents, and gene features revealed similarities among the 13 species. The number of oligonucleotide repeats was uncorrelated with genome size or phylogenetic position of the species. Phylogenetic analyses corroborated the monophyly of Aroideae but were unable to resolve the positions of Calla and Schismatoglottis.     

The complete mitochondrial genome of cricket Sclerogryllus punctatus (Orthoptera: Gryllidae) and phylogenetic analysis

The crickets of genus Sclerogryllus Gorochov, 1985 belongs to subfamily Sclerogryllinae of family Gryllidae. In this study, we report the first complete mitogenome sequences of the genus Sclerogryllus, and analyze the features of mitogenomes of S. punctatus. The mitogenome of S. punctatus was 15,438 bp and consisted of 37 genes, coding for 13 proteins, 2 ribosomal RNA (rRNA) and 22 transfer RNA (tRNA), and a control region. S. punctatus shares the arrangement of trnE-trnS-trnN with most mitogenomes of Grylloidea. Besides, the tRNAs possess the typical cloverleaf secondary structure except for the trnS1 (AGN) gene. The phylogenetic analysis using 13 protein-coding genes and 2 rRNA represents that genus Sclerogryllus is included in subfamily Gryllinae. Our results uncover the phylogenetic position of genus Sclerogryllus by mitogenome data within the family Gryllidae.     

Complete mitochondrial genome of the hemp borer,Grapholita delineana (Lepidoptera: Tortricidae): Gene variability and phylogeny among Grapholita

The mitochondrial genome (mitogenome) has been extensively used in phylogenetics and species-level evolutionary investigations. The lepidopteran family Tortricidae (leaf-roller moths), including the genus Grapholita, contains numerous species of economic importance, but for the majority of Grapholita species, their mitogenomes remain poorly studied. Here, we sequence and annotate the full mitogenome of Grapholita delineana, an important pest of hemp worldwide and compare it with the mitogenomes of two congeneric species available from GenBank. The G. delineana mitogenome is 15,599 bp long, including 37 typical mitochondrial genes and an A + T-rich region. Gene content, order and orientation are identical to other reported tortricid mitogenomes. Analyses of nucleotide diversity, Ka/Ks, genetic distance and number of variable sites together suggest that nad6 is the fastest-evolving gene among the mitochondrial PCGs of Grapholita. Our analyses indicate that Grapholita, as presently defined, is not monophyletic, confirming previous morphological and multiple-gene studies, using mitogenomic evidence. Our study provides information on comparative mitogenomics of Tortricidae especially Grapholita.     

Two novel cricetine mitogenomes: Insight into the mitogenomic characteristics and phylogeny in Cricetinae (Rodentia: Cricetidae)

Both Cricetus cricetus and Phodopus sungorus mitochondrial genomes (mitogenomes) were sequenced and elaborated for the first time in the present study. Their mitogenomes contained 37 genes and showed typical characteristics of the vertebrate mitogenome. Comparative analysis of 10 cricetine mitogenomes indicated that they shared similar characteristics with those of other cricetines in terms of genes arrangement, nucleotide composition, codon usage, tRNA structure, nucleotide skew and the origin of replication of light strand. Phylogenetic relationship of the subfamily Cricetinae was reconstructed using mitogenomes data with the methods of Bayesian Inference and Maximum Likelihood. Phylogenetic analysis indicated that Cricetulus kamensis was at basal position and phylogenetically distant from all other Cricetulus species but had a close relationship with the group of Phodopus, and supported that the genus Urocricetus deserved as a separate genus rank. The phylogenetic status of Tscherskia triton represented a separate clade corresponding to a diversified cricetine lineage (Cricetulus, Allocricetulus, and Cricetus).     

Phylogenomics and multigene phylogenies decipher two new cryptic marine algae from California,Gelidium gabrielsonii and G. kathyanniae (Gelidiales,Rhodophyta)

Molecular surveys are leading to the discovery of many new cryptic species of marine algae. This is particularly true for red algal intertidal species, which exhibit a high degree of morphological convergence. DNA sequencing of recent collections of Gelidium along the coast of California, USA, identified two morphologically similar entities that differed in DNA sequence from existing species. To characterize the two new species of Gelidium and to determine their evolutionary relationships to other known taxa, phylogenomic, multigene analyses, and morphological observations were performed. Three complete mitogenomes and five plastid genomes were deciphered, including those from the new species candidates and the type materials of two closely related congeners. The mitogenomes contained 45 genes and had similar lengths (24,963–24,964 bp). The plastid genomes contained 232 genes and were roughly similar in size (175,499–177,099 bp). The organellar genomes showed a high level of gene synteny. The two Gelidium species are diminutive, turf‐forming, and superficially resemble several long established species from the Pacific Ocean. The phylogenomic analysis, multigene phylogeny, and morphological evidence confirms the recognition and naming of two new species, describe herein as G. gabrielsonii and G. kathyanniae. On the basis of the monophyly of G. coulteri, G. gabrielsonii, G. galapagense, and G. kathyanniae, we suggest that this lineage likely evolved in California. Organellar genomes provide a powerful tool for discovering cryptic intertidal species and they continue to improve our understanding of the evolutionary biology of red algae and the systematics of the Gelidiales.     

The complete mitochondrial genome of Diaphorina citri (Hemiptera: Psyllidae) and phylogenetic analysis

The Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Psyllidae) is the most serious pest of citrus as the vector of Huanglongbing (HLB), the citrus greening disease. In this study, the complete mitochondrial genome (mitogenome) of D. citri has been sequenced and annotated, and a comparative analysis is provided with known Psylloidea species. The mitogenome of D. citri is a typical circular molecular of 15,038 bp in length with an A + T content of 74.56%, contains the typical 37 genes and the gene order is identical to the other Psylloidea mitogenomes. The nucleotide composition and codon usage of D. citri are similar to the four Psylloidea species. All protein-coding genes (PCGs) use standard initiation codons (TAN), stop with TAA and TAG except ND2 and ND5 which stop with incomplete termination codon T. All tRNAs have the typical clover-leaf structure, with the exception of trnS1 lacking the dihydrouridine (DHU) arm. The control region is located between rrnL and the trnI gene with the highest A + T content among the five Psylloidea species. Phylogenetic analysis is conducted based on the 13 PCGs or/and 2 rRNAs of 23 Sternorrhycha mitogenomes. Both maximum likelihood (ML) and Bayesian inference (BI) analysis suggest a clear relationship of Psylloidea, Aleyrodoidea and Aphidoidea within Sternorrhycha, which support the traditional morphological classification.     

Mitogenome evolution in Cephini (Hymenoptera: Cephidae): Evidence for parallel adaptive evolution

Two mitogenomes of remarkable size of the stem-boring sawflies, Trachelus tabidus (18539 bp) and T. iudaicus (20730 bp), were characterized and compared with previously known mitogenomes of Cephini. A rearrangement in the IQM gene cluster, an unusual elongation in rrnS gene and supposedly functional tandem repeat sequences in the A+T-rich region are synapomorphies of Cephini. Mitogenome evolution of the Cephini was investigated in a dataset of seven species representing all genera. The noticeably divergent mitogenomes of the Cephini both at nucleotide and amino acid level, and the variable rates of ω values, indicate the effects of different selective forces, rather than neutral evolution. The effect of positive selection is revealed by radical amino acid changes with ten physicochemical properties (P_α, B_r, pK′, h, E_l, H_nc, α_c, α_n, R_a and H_p) throughout each branch of the tree. The radical changes in the ND and ATP complexes are mostly related to chemical and energetic properties due to different energy requirements among species. The phenology of Cephini species and their host plants, together with habitat structure, rather than host plant preferences or phylogenetic constraint, are suggested to explain the role of adaptive evolution in shaping Cephini mitogenomes.     

Complete mitochondrial genome of the black‐tailed hornet,Vespa ducalis (Hymenoptera: Vespidae): Genomic comparisons in Vespoidea

We sequenced the complete mitochondrial genome (mitogenome) of the black‐tailed hornet, Vespa ducalis (Hymenoptera: Vespidae). The genome was 15,779‐bp long and contained typical sets of genes [13 protein‐coding genes (PCGs), 22 tRNAs, and 2 rRNAs]. The V. ducalis A + T‐rich region was 166‐bp long and was the shortest of all sequenced Vespoidea genomes, including Vespa. The genome was highly biased toward A/T nucleotides—80.1 % in the whole genome, 77.8 % in PCGs, 83.4–85.6 % in RNAs, and 92.8 % in the A + T‐rich region. These values are well within the typical range for genes and regions of Vespoidea mitogenomes. Start and stop codons in several Vespa species—including V. ducalis—were diversified, despite these species belonging to the same genus. In comparison with the ancestral mitogenomes, Vespa mitogenomes—including that of V. ducalis—showed substantial gene rearrangement; however, we detected no gene rearrangement among Vespa species. We conducted phylogenetic reconstruction based on concatenated sequences of 13 PCGs and two rRNAs (12,755 bp ) in available species of Vespoidea—21 species in six subfamilies in two families (Vespidae and Formicidae). The Bayesian inference and maximum likelihood (ML) methods revealed that each family formed strong monophyletic groups [Bayesian posterior probability (BPP) = 1; ML, 100 %]. Moreover, V. ducalis and V. mandarinia formed a strong sister group (BPP = 1; ML, 94 %).