High-level phylogeny of the Coleoptera inferred with mitochondrial genome sequences

Two new complete mitochondrial genomes of Dorcus stag beetles (Coleoptera, Lucanidae)

The systematics of Dorcus MacLeay has been a long-standing debate. Mitochondrial genomes were widely used to deeply understand the phylogeny of problematic taxa in virtue of their genetic importance and comprehensiveness. To provide more useful genetic data for resolving the systematic disputation of Dorcus stag beetles. The complete mitochondrial genomes of Dorcus hopei and Dorcus seguyi were obtained using the next generation sequencing. Characteristics of the two genomes are explicated through comparing their genome organization and base composition, protein-coding genes and codon usage, intergenic spacers and non-coding region, transfer and ribosomal RNA genes and control region. Phylogenetic relationships were reconstructed using Maximum likelihood and Bayesian inference analyses based on the concatenated nucleotide sequences of 13 PCGs from 9 stag beetles and 3 scarab beetles. The complete mitogenomes of D. hopei and D. seguyi was 16,026 bp/17,955 bp long, respectively. A tandem repeat with the length of 940 bp was presented in the A+T-rich region in D. hopei. An unexpected non-coding region of 332 bp was located between nad2 and trnW in D. seguyi. The phylogenetic analyses robustly supported that D. hopei formed a branch with the generic type of D. parallelipipedus. Whereas D. seguyi was not covered in the branch of (D. hopei + D. parallelipipedus), but was sister to them. The results indicated that D. hopei should be a good member of Dorcus MacLeay. The taxonomic status of D. seguyi remained to be studied furtherly.

Compositional and mutational rate heterogeneity in mitochondrial genomes and its effect on the phylogenetic inferences of Cimicomorpha (Hemiptera: Heteroptera)

Background

Mitochondrial genome (mt-genome) data can potentially return artefactual relationships in the higher-level phylogenetic inference of insects due to the biases of accelerated substitution rates and compositional heterogeneity. Previous studies based on mt-genome data alone showed a paraphyly of Cimicomorpha (Insecta, Hemiptera) due to the positions of the families Tingidae and Reduviidae rather than the monophyly that was supported based on morphological characters, morphological and molecular combined data and large scale molecular datasets. Various strategies have been proposed to ameliorate the effects of potential mt-genome biases, including dense taxon sampling, removal of third codon positions or purine-pyrimidine coding and the use of site-heterogeneous models. In this study, we sequenced the mt-genomes of five additional Tingidae species and discussed the compositional and mutational rate heterogeneity in mt-genomes and its effect on the phylogenetic inferences of Cimicomorpha by implementing the bias-reduction strategies mentioned above.

Results

Heterogeneity in nucleotide composition and mutational biases were found in mt protein-coding genes, and the third codon exhibited high levels of saturation. Dense taxon sampling of Tingidae and Reduviidae and the other common strategies mentioned above were insufficient to recover the monophyly of the well-established group Cimicomorpha. When the sites with weak phylogenetic signals in the dataset were removed, the remaining dataset of mt-genomes can support the monophyly of Cimicomorpha; this support demonstrates that mt-genomes possess strong phylogenetic signals for the inference of higher-level phylogeny of this group. Comparison of the ratio of the removal of amino acids for each PCG showed that ATP8 has the highest ratio while CO1 has the lowest. This pattern is largely congruent with the evolutionary rate of 13 PCGs that ATP8 represents the highest evolutionary rate, whereas CO1 appears to be the lowest. Notably, the value of Ka/Ks ratios of all PCGs is less than 1, indicating that these genes are likely evolving under purifying selection.

Conclusions

Our results demonstrate that mt-genomes have sites with strong phylogenetic signals for the inference of higher-level phylogeny of Cimicomorpha. Consequently, bioinformatic approaches to removing sites with weak phylogenetic signals in mt-genome without relying on an a priori tree topology would greatly improve this field.

Electronic supplementary material

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

DNA barcoding evaluation and implications for phylogenetic relationships in ladybird beetles (Coleoptera: Coccinellidae)

Ladybird beetles (Coleoptera: Coccinellidae), with broad morphological diversity, wide geographic distribution and substantial agricultural significance, are a challenging group for taxonomists and phylogenetics. As a promising tool to identify and discover new species, DNA barcoding might offer significant potential for identification, taxonomy and phylogeny of ladybird beetles. In the present study, a total of 1364 COI (cytochrome C oxidase subunit I) sequences representing 128 species from 52 genera of ladybird beetles were screened for barcoding evaluation and phylogenetic analysis. Our results from the barcoding analysis revealed that COI displays a similar level of species identification efficiency (nearly 90%) either based on Kimura two-parameter (K2P) distances calculation or on simplified neighbour-joining (NJ) tree construction. The phylogenetic relationships within the family Coccinellidae was analyzed by Bayesian-inference (BI) method. The phylogenetic results confirmed the monophyly of the subfamilies Microweisinae and Coccinellinae sensu Ślipiński (2007), and suggested that the subfamilies Coccidulinae, Chilocorinae and Scymninae are paraphyletic. However, the phylogenetic relationships among different subfamilies are not clearly defined and thus remain to be thoroughly studied. Overall, our study confirmed the usefulness of DNA barcoding for coccinellid species identification and phylogenetic inference.

The complete mitochondrial genomes of hycleus cichorii and hycleus phaleratus (Coleoptera: Meloidae)

The Hycleus cichorii and Hycleus phaleratus are two species of medicinal meloids widely distributed in southwest of China. We sequenced an annotated the complete mitochondrial genomes of H. cichorii and H. phaleratus, and the mitogenomes are 15,847 and 16,004 bp in length, respectively. Every mitochondrial genome encodes 13 proteins, 2 ribosomal RNAs, 22 tRNAs, and a control region with the identical arrangement to other beetles. The preliminary phylogenetic analysis with mitochondrial genomes of nine meloid species further confirmed the status of these two species.

Molecular phylogeny of Systellognatha (Plecoptera: Arctoperlaria) inferred from mitochondrial genome sequences

The infraorder Systellognatha is the most species-rich clade in the insect order Plecoptera and includes six families in two superfamilies: Pteronarcyoidea (Pteronarcyidae, Peltoperlidae, and Styloperlidae) and Perloidea (Perlidae, Perlodidae, and Chloroperlidae). To resolve the debatable phylogeny of Systellognatha, we carried out the first mitochondrial phylogenetic analysis covering all the six families, including three newly sequenced mitogenomes from two families (Perlodidae and Peltoperlidae) and 15 published mitogenomes. The three newly reported mitogenomes share conserved mitogenomic features with other sequenced stoneflies. For phylogenetic analyses, we assembled five datasets with two inference methods to assess their influence on topology and nodal support within Systellognatha. The results indicated that inclusion of the third codon positions of PCGs, exclusion of rRNA genes, the use of nucleotide datasets and Bayesian inference could improve the phylogenetic reconstruction of Systellognatha. The monophyly of Perloidea was supported in the mitochondrial phylogeny, but Pteronarcyoidea was recovered as paraphyletic and remained controversial. In this mitochondrial phylogenetic study, the relationships within Systellognatha were recovered as (((Perlidae + (Perlodidae + Chloroperlidae)) + (Pteronarcyidae + Styloperlidae)) + Peltoperlidae).

Complete mitochondrial genome of the bamboo snout beetle, Cyrotrachelus buqueti (Coleoptera: Curculionidae)

The bamboo snout beetle Cyrotrachelus buqueti (Coleoptera: Curculionidae) is a destructive forest pest and distributed widely in Southeast Asia. The 15,035 bp complete mitochondrial genome of the species consists of 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 21 transfer RNA genes (tRNAs) and a control region (GenBank accession no. MG674390). The trnl gene was not found in the C. buqueti mitogenome. The gene order and the orientation of C. buqueti were similar to those found in other Coleoptera species. The nucleotide composition was significantly biased (A, G, C, and T was 38.18%, 10.10%, 16.16%, and 35.56%, respectively) with A + T contents of 73.74%. ATG, ATA, ATT, AAT, and TTG were initiation codons and TAA, TAG, and T were termination codons. All the 21 tRNAs displayed a typical cloverleaf secondary structure, except for trnS1 which lacked the dihydrouridine arm. Phylogenetic analysis was performed using 13 PCGs with 14 other beetles showed that C. buqueti is closely related to Eucryptorhynchus brandti, which agree with the traditional classification.

The first two mitochondrial genomes from Taeniopterygidae (Insecta: Plecoptera): Structural features and phylogenetic implications

The complete mitochondrial genomes (mitogenomes) of Taeniopteryx ugola and Doddsia occidentalis (Plecoptera: Taeniopterygidae) were firstly sequenced from the family Taeniopterygidae. The 15,353-bp long mitogenome of T. ugola and the 16,020-bp long mitogenome of D. occidentalis each contained 37 genes including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (CR). The mitochondrial gene arrangement of the two taeniopterygids and other stoneflies was identical with the putative ancestral mitogenome of Drosophila yakuba. Most PCGs used standard ATN start codons and TAN termination codons. Twenty-one of the 22 tRNAs in each mitogenome could fold into the cloverleaf secondary structures, while the dihydrouridine (DHU) arm of trnSer (AGN) was reduced or absent. Stem-loop (SL) structures, poly-T stretch, poly-[AT]_n stretch and tandem repeats were found in the CRs of the two mitogenomes. The phylogenetic analyses using Bayesian inference (BI) and maximum likelihood methods (ML) generated identical results, both supporting the monophyly of all stonefly families and the two infraorders, Systellognatha and Euholognatha. Taeniopterygidae was grouped with another two families from Euholognatha. The relationships within Plecoptera were recovered as (((Perlidae+Peltoperlidae)+((Pteronarcyidae+Chloroperlidae)+Styloperlidae))+((Capniidae+Taeniopterygidae)+Nemouridae))+Gripopterygidae.

Mitochondrial genomic variation and phylogenetic relationships of three groups in the genus Scaphoideus (Hemiptera: Cicadellidae: Deltocephalinae)

The widespread leafhopper genus Scaphoideus Uhler is the most diverse genus in Scaphoideini and includes some species that are serious pests and vectors of plant pathogens. Here the first Scaphoideus mitogenome sequences are provided for three species, S. maai, S. nigrivalveus and S. varius, representing three main species groups in the Oriental region based on color pattern. The lengths of these three mitogenomes were 15,188, 15,235 and 15,207 bp, respectively. Gene order of three mitogenomes is highly conserved and identical to that of the putative ancestral insect. All three mitogenomes exhibited similar AT nucleotide bias, AT-, GC-skews and codon usage. One large 101 bp intergenic spacer between trnY and cox1 was in S. varius. All 22 tRNA genes had typical cloverleaf secondary structures, except for trnS1 (AGN) which appears to lack the dihydrouridine arm. Genes atp8, nad6 and nad2 were highly variable while cox1 showed the lowest nucleotide diversity. Phylogenetic analyses of three concatenated nucleotide datasets using maximum likelihood and Bayesian methods, comprising all 13 mitogenomes currently available for Membracoidea plus mitogenomes for eight outgroup species representing other cicadomorphan superfamilies, yielded the same topology in which Scaphoideus species formed a monophyletic group within a larger clade comprising three other included Deltocephalinae.

Characterization of the complete mitochondrial genomes of Maiestas dorsalis and Japananus hyalinus (Hemiptera: Cicadellidae) and comparison with other Membracoidea

Only six mitochondrial genomes (mitogenomes) have been previously published for Cicadellidae, the largest family of Hemiptera. This study provides complete, annotated mitogenomes of two additional cicadellid, species Maiestas dorsalis and Japananus hyalinus, and the first comparative mitogenome analysis across the superfamily Membracoidea. The mitogenomes of both sequenced species are similar to those of other studied hemipteran mitogenomes in organization and the lengths are 15,352 and 15,364 bp with an A + T content of 78.7% and 76.6%, respectively. In M. dorsalis, all sequenced genes are arranged in the putative ancestral insect gene arrangement, while the tRNA cluster trnW-trnC-trnY is rearranged to trnY-trnW-trnC in J. hyalinus, the first reported gene rearrangement in Membracoidea. Phylogenetic analyses of the 11 available membracoid mitogenomes and outgroups representing the other two cicadomorphan superfamilies supported the monophyly of Membracoidea, and indicated that treehoppers are a derived lineage of leafhoppers. ML and BI analyses yielded topologies that were congruent except for relationships among included representatives of subfamily Deltocephalinae. Exclusion of third codon positions of PCGs improved some node support values in ML analyses.

Mitochondrial genomes of blister beetles (Coleoptera, Meloidae) and two large intergenic spacers in Hycleus genera

BACKGROUND

Insect mitochondrial genomes (mitogenomes) exhibit high diversity in some lineages. The gene rearrangement and large intergenic spacer (IGS) have been reported in several Coleopteran species, although very little is known about mitogenomes of Meloidae.

RESULTS

We determined complete or nearly complete mitogenomes of seven meloid species. The circular genomes encode 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and two ribosomal RNAs (rRNAs), and contain a control region, with gene arrangement identical to the ancestral type for insects. The evolutionary rates of all PCGs indicate that their evolution is based on purifying selection. The comparison of tRNA secondary structures indicates diverse substitution patterns in Meloidae. Remarkably, all mitogenomes of the three studied Hycleus species contain two large intergenic spacers (IGSs). IGS1 is located between trnW and trnC, including a 9 bp consensus motif. IGS2 is located between trnS2 (UCN) and nad1, containing discontinuous repeats of a pentanucleotide motif and two 18-bp repeat units in both ends. To date, IGS2 is found only in genera Hycleus across all published Coleopteran mitogenomes. The duplication/random loss model and slipped-strand mispairing are proposed as evolutionary mechanisms for the two IGSs (IGS1, IGS2). The phylogenetic analyses using MrBayes, RAxML, and PhyloBayes methods based on nucleotide and amino acid datasets of 13 PCGs from all published mitogenomes of Tenebrionoids, consistently recover the monophylies of Meloidae and Tenebrionidae. Within Meloidae, the genus Lytta clusters with Epicauta rather than with Mylabris. Although data collected thus far could not resolve the phylogenetic relationships within Meloidae, this study will assist in future mapping of the Meloidae phylogeny.

CONCLUSIONS

This study presents mitogenomes of seven meloid beetles. New mitogenomes retain the genomic architecture of the Coleopteran ancestor, but contain two IGSs in the three studied Hycleus species. Comparative analyses of two IGSs suggest that their evolutionary mechanisms are duplication/random loss model and slipped-strand mispairing.

Mitochondrial genome of Sitona callosus (Coleoptera: Curculionidae) and phylogenetic analysis within Entiminae

In this study, we sequenced and annotated the nearly complete mitochondrial genome (mitogenome) of Sitona callosus (Coleoptera: Curculionidae). This mitogenome was 14,333 bp long and encoded 13 protein-coding genes, 19 transfer RNA genes (tRNAs), and two ribosomal RNA unit genes. Gene rearrangements were presented in a tRNA cluster of six tRNAs between nad3 and nad5, i.e. the ancestral order ARNSEF was changed to be RNSAEF. All tRNAs had a typical secondary cloverleaf structure, except for trnS1 which lacked the dihydrouridine arm. The Bayesian phylogenetic tree of 11 Entiminae species based on the concatenated nucleotide sequences of 13 PCGs showed that S. callosus and S. lineatus formed a clade which was at the basal position in the Entiminae phylogeny.

Comparative mitogenomic analysis of mirid bugs (Hemiptera: Miridae) and evaluation of potential DNA barcoding markers

The family Miridae is one of the most species-rich families of insects. To better understand the diversity and evolution of mirids, we determined the mitogenome of Lygus pratenszs and re-sequenced the mitogenomes of four mirids (i.e., Apolygus lucorum, Adelphocoris suturalis, Ade. fasciaticollis and Ade. lineolatus). We performed a comparative analysis for 15 mitogenomic sequences representing 11 species of five genera within Miridae and evaluated the potential of these mitochondrial genes as molecular markers. Our results showed that the general mitogenomic features (gene content, gene arrangement, base composition and codon usage) were well conserved among these mirids. Four protein-coding genes (PCGs) (cox1, cox3, nad1 and nad3) had no length variability, where nad5 showed the largest size variation; no intraspecific length variation was found in PCGs. Two PCGs (nad4 and nad5) showed relatively high substitution rates at the nucleotide and amino acid levels, where cox1 had the lowest substitution rate. The Ka/Ks values for all PCGs were far lower than 1 (<0.59), but the Ka/Ks values of cox1-barcode sequences were always larger than 1 (1.34 -15.20), indicating that the 658 bp sequences of cox1 may be not the appropriate marker due to positive selection or selection relaxation. Phylogenetic analyses based on two concatenated mitogenomic datasets consistently supported the relationship of Nesidiocoris + (Trigonotylus + (Adelphocoris + (Apolygus + Lygus))), as revealed by nad4, nad5, rrnL and the combined 22 transfer RNA genes (tRNAs), respectively. Taken sequence length, substitution rate and phylogenetic signal together, the individual genes (nad4, nad5 and rrnL) and the combined 22 tRNAs could been used as potential molecular markers for Miridae at various taxonomic levels. Our results suggest that it is essential to evaluate and select suitable markers for different taxa groups when performing phylogenetic, population genetic and species identification studies.

The complete mitochondrial genome of the cigarette beetle, Lasioderma serricorne (Coleoptera: Anobiidae)

The cigarette beetle, Lasioderma serricorne (Fabricius), is an important pest of stored commodities and distributed widely in the world. Here, we report the complete mitochondrial genome of L. serricorne which was 15,958 bp and composed of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The gene order and orientation of L. serricorne were identical to those of other Coleopteran mitogenomes. ATG, ATA, ATT, ATC, TTG were initiation codons and TAA, TAG, T were termination codons. All 22 tRNA genes were predicted with a typical cloverleaf structure except for trnS1 (AGN). Phylogenetic analysis performed using 13 PCGs with 14 other beetles showed that L. serricorne is closely related to Stegobium paniceum, which agree with the conventional taxonomy.

The complete mitochondrial genome of Sympiezomias velatus (Coleoptera: Curculionidae)

Here, we determined the complete mitogenome sequence of Sympiezomias velatus (Coleoptera: Curculionidae: Entiminae). This mitogenome was 15,592 bp long with an A + T content of 74.1% and contains 13 protein-coding genes (PCGs), 21 transfer RNA genes (tRNAs), 2 ribosomal RNA unit genes and a large non-coding region (putative control region). The trnI gene was not found in the S. velatus mitogenome. The order and orientation of the mitochondrial genes were identical to the inferred ancestral arrangement of insects except for trnR which was changed to be adjacent the nad3 gene. All tRNAs had the typical cloverleaf structure, except for trnS1 which lacked the dihydrouridine arm. The Bayesian phylogenetic tree of 10 Entiminae species based on the concatenated nucleotide sequences of 13 PCGs strong supported a sister relationship of S. velatus and Barynotus obscures.

A transfer RNA gene rearrangement in the lepidopteran mitochondrial genome

Gene arrangements in the mitochondrial genomes (mitogenomes) of insects are conserved across the major lineages, but can be rearranged within derived groups and may provide valuable phylogenetic characters. In this study, we sequenced the entire mitogenome of Parasa consocia, a moth of the family Limacodidae (Lepidoptera: Zygaenoidea). Compared with other lepidopterans and ancestral insects, the P. consocia mitogenome features a transfer RNA gene arrangement novel among lepidopterans between the ND3 and ND5 genes: RANSEF (the underline signifies an inverted gene), which differs from the ARNSEF arrangement of ancestral insects. This rearrangement can be explained by the tandem duplication-random loss model. We inferred a phylogenetic hypothesis for the lepidopteran superfamily based on mitochondrial amino-acid sequences using the Bayesian-inference and maximum-likelihood methods. Our results showed that P. consocia belongs to the Zygaenoidea superfamily and supported the following phylogenetic relationship: Yponomeutoidea + (Tortricoidea + Zygaenoidea + (Papilionoidea + (Pyraloidea + (Noctuoidea + (Geometroidea + Bombycoidea)))))). Comparative analyses indicated that mitogenomes are a useful phylogenetic tool at the subfamily level within the order Lepidoptera. Our findings also suggest that mitogenomes are likely to represent a valuable tool for systematics in other groups of lepidopterans.

Transcriptome sequencing and phylogenetic analysis of four species of luminescent beetles

The evolution of bioluminescence has prompted scientific attention to illuminate phylogenetic relationships of luminescent beetles. However, genomic resources are virtually lacking in rhagophthalmids (Rhagophthalmidae) and their related firefly beetles lampyrids (Lampyridae). Here, we employed the Illumina Hiseq 2000 platform and sequenced the whole-body transcriptomes of the four luminescent beetles: one rhagophthalmid (Rhagophthalmus sp.) and three fireflies (Asymmetricata circumdata, Aquatica ficta, and Pyrocoelia pectoralis). We obtained 55.4, 43.4, 38.6, and 36.7 million clean reads for the four species, respectively. All reads were assembled into contigs from which unigenes were derived. All unigenes were annotated by publicly available databases, and a total of 4325 orthologous genes were identified. Using multiple phylogenetic approaches, our transcriptome data confirmed the distinctiveness of Rhagophthalmidae from Lampyridae, which was also supported by our mitogenome analysis using three newly determined mitogenome sequences and 12 previously published ones. Together, this study is the first report of whole transcriptome sequencing data in Rhagophthalmidae and Lampyridae species, representing a valuable genomic resource for studying the origin and evolution of some remarkable traits in these beetles such as bioluminescence. Moreover, our transcriptome and mitogenome data provide useful phylogenetic information that could be of importance in future studies of phylogenetic inference.

New Mitogenomes of Two Chinese Stag Beetles (Coleoptera, Lucanidae) and Their Implications for Systematics

Although conspicuous and well-studied, stag beetles have been slow to join the genomic era. In this study, mitochondrial genomes of two stag beetles, Sinodendron yunnanense and Prosopocoilus confucius, are sequenced for the first time. Both of their genomes consisted of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region. The mitogenome of S. yunnanense was 16,921 bp in length, and P. confucius was 16,951 bp. The location of the gene trnL(UUR), between the A + T-rich and control region in S. yunnanense, is the first observed in Lucanidae. In P. confucius, an unexpected noncoding region of 580 bp was discovered. Maximum likelihood and Bayesian inference on the 13 mitochondrial PCGs were used to infer the phylogenetic relationships among 12 representative stag beetles and three scarab beetles. The topology of the two phylogenetic trees was almost identical: S. yunnanense was recovered as the most basal Lucanid, and the genus Prosopocoilus was polyphyletic due to P. gracilis being recovered sister to the genera Dorcus and Hemisodorcus. The phylogenetic results, genetic distances and mitogenomic characteristics call into question the cohesion of the genus Prosopocoilus. The genetic resources and findings herein attempts to redress understudied systematics and mitogenomics of the stag beetles.

The complete mitochondrial genome of Stegobium paniceum (Coleoptera: Anobiidae)

The complete mitochondrial genome of Stegobium paniceum (Coleoptera: Anobiidae) is a circular DNA molecule of 15,271 bp (GenBank accession number XK819317), and its nucleotide composition is biased towards A + T nucleotides (78.32%). This genome comprises 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA (tRNA) genes, and an A + T-rich region. The gene order of S. paniceum was similar to those found in other known Coleoptera species. Sixteen reading frame overlaps and six intergenic regions were found in the mitochondrial genome of S. paniceum. All 22 tRNA genes have the typical cloverleaf secondary structure, with an exception for trnS₁ (AGN). The phylogenetic relationships based on neighbour-joining method revealed that S. paniceum is closely related to Apatides fortis, which is consistent with the traditional morphological classification.