Complete mitochondrial genome sequence of the yellow-spotted long-horned beetle Psacothea hilaris (Coleoptera: Cerambycidae) and phylogenetic analysis among coleopteran insects

The first mitogenomic phylogenetic framework of Dorcus sensu lato (Coleoptera: Lucanidae), with an emphasis on generic taxonomy in Eastern Asia

BACKGROUND

Dorcus stag beetles in broad sense are one of the most diverse group in Lucanidae and important saproxylic insects playing a crucial role in nutrient recycling and forest biomonitoring. However, the dazzling morphological differentiations have caused numerous systematic confusion within the big genus, especially the puzzlingly generic taxonomy. So far, there is lack of molecular phylogenetic study to address the chaotic situation. In this study, we undertook mitochondrial genome sequencing of 42 representative species including 18 newly-sequenced ones from Eastern Asia and reconstructed the phylogenetic framework of stag beetles in Dorcus sensu lato for the first time.

RESULTS

The mitogenome datasets of Dorcus species have indicated the variable mitogenomic lengths ranged from 15,785 to 19,813 bp. Each mitogenome contained 13 PCGs, 2 rRNAs, 22 tRNAs, and a control region, and all PCGs were under strong purifying selection (Ka/Ks < 1). Notably, we have identified the presence of a substantial intergenic spacer (IGS) between the trnAser (UCN) and NAD1 genes, with varying lengths ranging from 129 bp (in D. hansi) to 158 bp (in D. tityus). The mitogenomic phylogenetic analysis of 42 species showed that Eastern Asia Dorcus was monophyletic, and divided into eight clades with significant genetic distance. Four of them, Clade VIII, VII, VI and I are clustered by the representative species of Serrognathus Motschulsky, Kirchnerius Schenk, Falcicornis Séguy and Dorcus s.s. respectively, which supported their fully generic positions as the previous morphological study presented. The topology also showed the remaining clades were distinctly separated from the species of Dorcus sensu lato, which implied that each of them might demonstrate independent generic status. The Linnaeus nomenclatures were suggested as Eurydorcus Didier stat. res., Eurytrachellelus Didier stat. res., Hemisodorcus Thomson stat. res. and Velutinodorcus Maes stat. res. For Clade V, IV, III and II respectively.

CONCLUSION

This study recognized the monophyly of Dorcus stag beetles and provided a framework for the molecular phylogeny of this group for the first time. The newly generated mitogenomic data serves as a valuable resource for future investigations on lucanid beetles. The generic relationship would facilitate the systematics of Dorcus stag beetles and thus be useful for exploring their evolutionary, ecological, and conservation aspects.

The Phylogenetic Relationship of Lamiinae (Coleoptera: Cerambycidae) Using Mitochondrial Genomes

Lamiinae is the largest subfamily of the Cerambycidae (longhorn beetles), with approximately 21,863 described species. Previous phylogenetic studies of Lamiinae showed that this subfamily was monophyletic, but the relationship between the tribes of Lamiinae is still controversial. Partial molecular data and species morphological characteristics are not sufficient to resolve species phylogenetic studies perfectly. At the same time, the full mitochondrial genome contains more comprehensive genetic data. Benefiting from the development of next-generation sequencing (NGS), mitochondrial genomes can be easily acquired and used as reliable molecular markers to investigate phylogenetic relationships within Cerambycidae. Using NGS technology, we obtained 11 mitochondrial genome sequences of Lamiinae species. Based on this newly generated mitochondrial genome dataset matrix, we reconstructed the phylogeny of Lamiinae. The Bayesian Inference and Maximum Likelihood analyses strongly support the monophyly of four tribes (Lamiini, Batocerini, Mesosini, and Saperdini), whereas the tribe Acanthocinini was identified as paraphyletic. Other mitochondrial structural features were also observed: the start codon in the nad1 gene of all 11 mitochondrial genomes is TTG; 17-22 bp intergenic spacers (IGS) with a 'TACTA' motif were found between trnS2 and nad1. Moreover, two long IGS were found in Mesosa myops and Batocera sp. Tandem repeats were found in the IGS of Batocera sp.

Functional morphology of the mouthparts of longhorn beetle adult Psacothea hilaris (Coleoptera: Cerambycidae) and sensilla comparisons between the sexes

Psacothea hilaris is an important wood boring beetle that causes significant ecological and economic damage. The mouthparts of P. hilaris contain feeding and sensory structures that play important roles in many behaviors. The study of their functional morphologies provides insight into feeding and sensory mechanisms. The fine structures of the mouthparts and the sensilla of both sexes were observed by scanning electron microscopy, with special attention to quantitative comparisons. The general structures of the mouthparts are similar in males and females. However, the maxillary and labial palps of females are more well-developed than those of males. Six types of sensilla were found in both sexes: sensilla basiconca, sensilla trichodea, sensilla chaetica, sensilla digitiformia, sensilla coeloconica, and sensilla palmata. This was the first attempt to describe the sensilla on the epipharynx and ligula of Cerambycidae. There were differences in the number of sensilla of males and females, especially on the mandibles, maxillary and labial palps. However, not all types of sensilla have a greater number in females compared to males. The results provide basic information on the gustatory sensation mechanism of Cerambycidae.

Complete mitochondrial genome of Pseudoechthistatus hei (Coleoptera: Cerambycidae: Lamiinae) and its phylogenetic analysis

The little-known genus Pseudoechthistatus Pic, 1917 belongs to the subfamily Lamiinae of the family Cerambycidae. Adult Pseudoechthistatus hei Xie and W. Wang, 2019 specimens were collected from Bāijì Hill, Xīntángfáng Village, Wéixī County, Yúnnán Province, China. The complete mitochondrial genome (GenBank accession number: ON641973.1) of P. hei was sequenced, annotated, and characterized; it is a circular DNA molecule of 16,103 bp with a 75.71% AT content, and it comprised 13 protein-coding genes (PCG), 22 tRNA genes, two rRNA genes, and 1 control region. The PCGs initiated with the typical ATN (Met) start codons, and were terminated by typical TAN stop codons. The Bayesian Inference phylogenetic tree was first constructed using JTT + F + I + G4 model for P. hei, which showed that P. hei was closely related to Monochamus alternatus alternatus.

Complete mitochondrial genome of the longhorn date palm stem borer Jebusaea hammerschmidtii (Reiche, 1878)

The 15,619 bp mitochondrial genome of Jebusaea hammerschmidtii was assembled from short reads, annotated, and compared to the genomes of other longhorn beetles (Cerambycidae). Gene content was typical of animal mitochondrial genomes and contained 13 protein-coding, 22 tRNA, and 2 rRNA genes. Gene organization was identical to that of other longhorn beetles. Phylogenetic analysis placed J. hammerschmidtii within the subfamily Cerambycinae, and strongly supported the monophyly of the Cerambycinae, Lamiinae, and Prioninae subfamilies.

Complete Mitochondrial Genomes of Metcalfa pruinosa and Salurnis marginella (Hemiptera: Flatidae): Genomic Comparison and Phylogenetic Inference in Fulgoroidea

The complete mitochondrial genomes (mitogenomes) of two DNA barcode-defined haplotypes of Metcalfa pruinosa and one of Salurnis marginella (Hemiptera: Flatidae) were sequenced and compared to those of other Fulgoroidea species. Furthermore, the mitogenome sequences were used to reconstruct phylogenetic relationships among fulgoroid families. The three mitogenomes, including that of the available species of Flatidae, commonly possessed distinctive structures in the 1702-1836 bp A+T-rich region, such as two repeat regions at each end and a large centered nonrepeat region. All members of the superfamily Fulgoroidea, including the Flatidae, consistently possessed a motiflike sequence (TAGTA) at the ND1 and trnS₂ junction. The phylogenetic analyses consistently recovered the familial relationships of (((((Ricaniidae + Issidae) + Flatidae) + Fulgoridae) + Achilidae) + Derbidae) in the amino acid-based analysis, with the placement of Cixiidae and Delphacidae as the earliest-derived lineages of fulgoroid families, whereas the monophyly of Delphacidae was not congruent between tree-constructing algorithms.

Mitogenome Analysis of Four Lamiinae Species (Coleoptera: Cerambycidae) and Gene Expression Responses by Monochamus alternatus When Infected with the Parasitic Nematode, Bursaphelenchus mucronatus

We determined the mitochondrial gene sequence of Monochamus alternatus and three other mitogenomes of Lamiinae (Insect: Coleoptera: Cerambycidae) belonging to three genera (Aulaconotus, Apriona and Paraglenea) to enrich the mitochondrial genome database of Lamiinae and further explore the phylogenetic relationships within the subfamily. Phylogenetic trees of the Lamiinae were built using the Bayesian inference (BI) and maximum likelihood (ML) methods and the monophyly of Monochamus, Anoplophora, and Batocera genera was supported. Anoplophora chinensis, An. glabripennis and Aristobia reticulator were closely related, suggesting they may also be potential vectors for the transmission of the pine wood pathogenic nematode (Bursaphelenchus xylophilus) in addition to M. alternatus, a well-known vector of pine wilt disease. There is a special symbiotic relationship between M. alternatus and Bursaphelenchus xylophilus. As the native sympatric sibling species of B. xylophilus, B. mucronatus also has a specific relationship that is often overlooked. The analysis of mitochondrial gene expression aimed to explore the effect of B. mucronatus on the energy metabolism of the respiratory chain of M. alternatus adults. Using RT-qPCR, we determined and analyzed the expression of eight mitochondrial protein-coding genes (COI, COII, COIII, ND1, ND4, ND5, ATP6, and Cty b) between M. alternatus infected by B. mucronatus and M. alternatus without the nematode. Expression of all the eight mitochondrial genes were up-regulated, particularly the ND4 and ND5 gene, which were up-regulated by 4-5-fold (p < 0.01). Since longicorn beetles have immune responses to nematodes, we believe that their relationship should not be viewed as symbiotic, but classed as parasitic.

The complete mitochondrial genome of Spondylis buprestoides (Coleoptera: Cerambycidae)

Spondylis buprestoides is a major boring pest of Cerambycidae, which mainly damaging conifers and also can carry pine wood nematode, Bursaphelenchus xylophilus. In this study, the complete mitochondrial genome of S. buprestoides was determined by Illumina sequencing technology. The whole mitogenome was 15,837 bp in length with 20.05% GC content, which contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). Phylogenetic analysis showed that S. buprestoides was closely related to Lepturinae. The sequence data of S. buprestoides could provide useful genetic information for the studies on phylogenetic and evolutionary of Cerambycidae.

Mitochondrial genomes of twelve species of hyperdiverse Trigonopterus weevils

Mitochondrial genomes of twelve species of Trigonopterus weevils are presented, ten of them complete. We describe their gene order and molecular features and test their potential for reconstructing the phylogeny of this hyperdiverse genus comprising > 1,000 species. The complete mitochondrial genomes examined herein ranged from 16,501 bp to 21,007 bp in length, with an average AT content of 64.2% to 69.7%. Composition frequencies and skews were generally lower across species for atp6, cox1-3, and cob genes, while atp8 and genes coded on the minus strand showed much higher divergence at both nucleotide and amino acid levels. Most variation within genes was found at the codon level with high variation at third codon sites across species, and with lesser degree at the coding strand level. Two large non-coding regions were found, CR1 (between rrnS and trnI genes) and CR2 (between trnI and trnQ), but both with large variability in length; this peculiar structure of the non-coding region may be a derived character of Curculionoidea. The nad1 and cob genes exhibited an unusually high interspecific length variation of up to 24 bp near the 3' end. This pattern was probably caused by a single evolutionary event since both genes are only separated by trnS2 and length variation is extremely rare in mitochondrial protein coding genes. We inferred phylogenetic trees using protein coding gene sequences implementing both maximum likelihood and Bayesian approaches, each for both nucleotide and amino acid sequences. While some clades could be retrieved from all reconstructions with high confidence, there were also a number of differences and relatively low support for some basal nodes. The best partition scheme of the 13 protein coding sequences obtained by IQTREE suggested that phylogenetic signal is more accurate by splitting sequence variation at the codon site level as well as coding strand, rather than at the gene level. This result corroborated the different patterns found in Trigonopterus regarding to A+T frequencies and AT and GC skews that also greatly diverge at the codon site and coding strand levels.

The complete mitochondrial genome of Monochamus alternatus alternatus (Coleoptera: Cerambycidae)

Monochamus alternatus alternatus is the major vector of pinewood nematode, Bursaphelenchus xylophilus, in Asia. The length of the complete mitochondria genome of M. alternatus alternatus was 15,880 bp with 21% GC content, including 39.7% A, 12.3% C, 8.7% G and 39.3% T. There were 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and one AT-rich region. This study provides a useful genetic information for subsequent study of the differences between M. alternatus subspecies.

Comparative mitochondrial genome analysis of Dendrolimus houi (Lepidoptera: Lasiocampidae) and phylogenetic relationship among Lasiocampidae species

Dendrolimus houi is one of the most common caterpillars infesting Gymnosperm trees, and widely distributed in several countries in Southeast Asia, and exists soley or coexists with several congeners and some Lasiocampidae species in various forest habitats. However, natural hybrids occasionally occur among some closely related species in the same habitat, and host preference, extreme climate stress, and geographic isolation probably lead to their uncertain taxonomic consensus. The mitochondrial DNA (mtDNA) of D. houi was extracted and sequenced by using high-throughput technology, and the mitogenome composition and characteristics were compared and analyzed of these species, then the phylogenetic relationship was constructed using the maximum likelihood method (ML) and the Bayesian method (BI) based on their 13 protein-coding genes (PCGs) dataset, which were combined and made available to download which were combined and made available to download among global Lasiocampidae species data. Mitogenome of D. houi was 15,373 bp in length, with 37 genes, including 13 PCGs, 22 tRNA genes (tRNAs) and 2 rRNA genes (rRNAs). The positions and sequences of genes were consistent with those of most known Lasiocampidae species. The nucleotide composition was highly A+T biased, accounting for ~80% of the whole mitogenome. All start codons of PCGs belonged to typical start codons ATN except for COI which used CGA, and most stop codons ended with standard TAA or TAG, while COI, COII, ND4 ended with incomplete T. Only tRNASer (AGN) lacked DHU arm, while the remainder formed a typical "clover-shaped" secondary structure. For Lasiocampidae species, their complete mitochondrial genomes ranged from 15,281 to 15,570 bp in length, and all first genes started from trnM in the same direction. And base composition was biased toward A and T. Finally, both two methods (ML and BI) separately revealed that the same phylogenetic relationship of D. spp. as ((((D. punctatus + D. tabulaeformis) + D. spectabilis) + D. superans) + (D. kikuchii of Hunan population + D. houi) as in previous research, but results were different in that D. kikuchii from a Yunnan population was included, indicating that different geographical populations of insects have differentiated. And the phylogenetic relationship among Lasiocampidae species was ((((Dendrolimus) + Kunugia) + Euthrix) + Trabala). This provides a better theoretical basis for Lasiocampidae evolution and classification for future research directions.

Characterization of the Complete Mitochondrial Genome of Epicauta impressicornis (Coleoptera: Meloidae) and Its Phylogenetic Implications for the Infraorder Cucujiformia

The complete mitochondrial genome (mitogenome) of Epicauta impressicornis Pic (Coleoptera: Meloidae) was determined. The circular genome is 15,713-bp long, and encodes 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a control region (CR). The 13 PCGs start with the typical ATN codon and terminate with the typical stop codon TAA (ND2, ND4L, ND6, ATP6, ATP8, and CYTB), TAG (ND1 and ND3), and T- (COX1, COX2, COX3, ND4, and ND5). The two rRNA genes (rrn12S and rrn16S) are encoded on the minority strand. All tRNAs genes except trnS1 (AGN) are predicted to fold into the typical cloverleaf structure. The longest overlap (10 bp) is observed between ATP8 and ATP6. CR mainly harbors a conserved poly-T stretch (15 bp), a short repeat unit (17 bp), some universal microsatellite-like repeats, and a canonical poly-A tail. Phylogenetic analysis using Bayesian inferences and maximum likelihood based on nucleotide and corresponding amino acid sequences of the 13 PCGs showed that E. impressicornis is closely related to E. chinensis, this relationship is and supported within Cucujiformia belonging to Meloidae (Tenebrionoidea). Our results further confirmed the monophyly of Tenebrionoidea, Lymexyloidea, Curculionoidea, Chrysomeloidea, Cucujoidea, Coccinelloidea, and Cleroidea within Cucujiformia, and revealed the sister relationships of (Cleroidea + Coccinelloidea), (Lymexyloidea + Tenebrionoidea), and ((Chrysomeloidea + Cucujoidea) + Curculionoidea). We believe that the complete mitogenome of E. impressicornis will contribute to further studies on molecular bases for the classification and phylogeny of Meloidae or even Cucujiformia.

The complete mitochondrial genome of Annamanum lunulatum (Coleoptera: Lamiinae) and its phylogeny

The complete mitochondrial genome of the Annamanum lunulatum is 15,610 bp in length, which contains 13 protein-coding genes, 22 transfer RNAs, two ribosomal RNAs, and the A + T-rich region. The arrangement of genes is identical to all know longhorn beetles mitochondrial genomes. The overall AT content of the mitochondrial genome is 75.3%, whereas the AT content of A + T-rich region is 84.3%. In ML and BI phylogenetic analyses, A. lunulatum is a sister clade to Blepephaeus succinctor, and the monophyly of Lamiinae is supported.

The complete mitochondrial genomes of five longicorn beetles (Coleoptera: Cerambycidae) and phylogenetic relationships within Cerambycidae

Cerambycidae is one of the most diversified groups within Coleoptera and includes nearly 35,000 known species. The relationships at the subfamily level within Cerambycidae have not been convincingly demonstrated and the gene rearrangement of mitochondrial genomes in Cerambycidae remains unclear due to the low numbers of sequenced mitogenomes. In the present study, we determined five complete mitogenomes of Cerambycidae and investigated the phylogenetic relationship among the subfamilies of Cerambycidae based on mitogenomes. The mitogenomic arrangement of all five species was identical to the ancestral Cerambycidae type without gene rearrangement. Remarkably, however, two large intergenic spacers were detected in the mitogenome of Pterolophia sp. ZJY-2019. The origins of these intergenic spacers could be explained by the slipped-strand mispairing and duplication/random loss models. A conserved motif was found between trnS2 and nad1 gene, which was proposed to be a binding site of a transcription termination peptide. Also, tandem repeat units were identified in the A + T-rich region of all five mitogenomes. The monophyly of Lamiinae and Prioninae was strongly supported by both MrBayes and RAxML analyses based on nucleotide datasets, whereas the Cerambycinae and Lepturinae were recovered as non-monophyletic.

The mitochondrial genomes of two walnut pests, Gastrolina depressa depressa and G. depressa thoracica (Coleoptera: Chrysomelidae), and phylogenetic analyses

In this study, the mitochondrial genomes (mitogenomes) of two walnut leaf insect pests, Gastrolina depressa depressa and G. depressa thoracica, were sequenced by Sanger sequencing technology. The mitogenome of G. depressa thoracica was complete at 16,109 bp in length, while the mitogenome of G. depressa depressa (14,277 bp) was partial. The genomic analyses indicated that both mitogenomes have the typical gene content and arrangement. The formerly identified elements, 'TAGTA' between trnSer(UCN) and nad2, and 'ATGATAA' between atp8 and atp6, were more conserved than that between nad4L and nad4, which was 'ATGTTAA' in Coleoptera excluding Polyphaga. Phylogenetic analyses of the 13 protein-coding genes from 36 coleopteran species well supported a close affinity between the subfamily Chrysomelinae including G. depressa thoracica and G. depressa depressa and Galerucinae, as well as a sister relationship of ((Eumolpinae + Cryptocephalinae) + Cassidinae) within Chrysomelidae.

The Complete Mitochondrial Genome of the Longhorn Beetle Dorysthenes paradoxus (Coleoptera: Cerambycidae: Prionini) and the Implication for the Phylogenetic Relationships of the Cerambycidae Species

The longhorn beetle Dorysthenes paradoxus (Faldermann, 1833) (Coleoptera: Cerambycidae) is not only a serious agricultural pest but also a traditionally edible insect in China. However, no genetic information on this species has been acquired. In the present study, we report the mitochondrial genome (mitogenome) of Do. paradoxus, as the first complete mitogenome of Prioninae. The circular mitogenome of 15,922 bp encodes 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), and two ribosomal RNAs (rRNAs), and it contains an A+T-rich region. This mitogenome exhibits the lowest A+T content (71.13%) but harbors the largest AT skew (0.116) among the completely sequenced Cerambycidae species. Eleven of the 13 PCGs have a typical ATN start codon, whereas COI and ND1 are tentatively designated by AAT and TTG, respectively. Only 4 of the 13 PCGs harbor a complete termination codon, and the remaining 9 possess incomplete termination codons (T or TA). Apart from tRNASer(AGN), the other 21 tRNAs can fold into a typical clover-leaf secondary structures. The Do. paradoxus A+T-rich region contains two poly-T stretches and a tandem repeat that comprises two 47-bp-long copies. Both Bayesian inference and Maximum likelihood analyses confirmed the subfamily ranks of Cerambycidae ([Prioninae + Cerambycinae] + Lamiinae) and the close relationship between Philinae and Prioninae/Cerambycinae. However, the data did not support the monophyly of Prioninae and Cerambycinae. The mitogenome presented here provides basic genetic information for this economically important species.

Complete mitochondrial genome of Callipogon relictus Semenov (Coleoptera: Cerambycidae): a natural monument and endangered species in Korea

In the present study, we report the first complete mitochondrial DNA genome of the genus Callipogon based on C. relictus, a natural monument and endangered species in South Korea. The mitogenome is 15,742 base pairs with 13 protein coding genes (PCGs), two rRNAs, 22 tRNAs, and a 1033 bp long AT-rich region. The overall base composition was 67.3% AT and 32.7 GC. Among 13 PCGs, seven genes (Nad2, Atp8, Atp6, Nad4L, Nad6, Cob, Nad1) harbour the typical stop codon TAA or TAG, whereas remaining five genes terminate with T. Interestingly, Cox3 employs AGA as the termination codon.

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.

Revisiting Coleoptera a + T-rich region: structural conservation, phylogenetic and phylogeographic approaches in mitochondrial control region of bioluminescent Elateridae species (Coleoptera)

The control region (CR) or A + T-rich region in Coleoptera mt genome is poorly characterized, including the Elateroidea bioluminescent species. Here, we provided the first attempt to characterize and compare the structure and organization of the CR of different species within Elateridae. We also revisited some sequenced Coleoptera CR and observed consensus T-stretches, non-conserved sequences near the stem-loop and unusual inner tRNAs-like sequences. All these features are probably involved in the replication start of the mt genome. The phylogenetic relationships in Elateridae bioluminescent groups using partial sequence of CR showed the monophyly of Pyrearinus pumilus group and Pyrearinus as a polyphyletic genus, corroborating our previous results. The wider genetic variation obtained by CR analysis could separate two different lineages that occur within P. termitilluminans populations. In Elateridae, the CR exhibited high polymorphism within and between populations, which was also observed in other Coleoptera species, suggesting that the CR could be described as a suitable molecular marker to be applied in phylogenetic and phylogeographic studies.

The complete mitochondrial genome sequence from the longicorn beetle Obrium sp. (Coleoptera: Cerambycidae)

The complete mitochondrial genome (mitogenome) of Obrium sp. was sequenced, which was the first representation from the longicorn beetle genus Obrium. The whole mitogenome is a typical circular DNA molecule of 15 680 bp, and has similar gene content and gene order to previously reported beetle species. In the sampled subfamilies of Cerambycidae, both Lamiinae and Cerambycinae were recovered with high supports in phylogenetic analysis. And the newly sequenced Obrium sp. was the sister group to Massicus raddei in Lamiinae.

Complete mitochondrial genome of Cryptolestes pusillus (Coleoptera: Laemophloeidae)

In this study, the complete mitochondrial genome of Cryptolestes pusillus (GenBank accession number KT070713) was sequenced by long PCR and primer walking methods. The total length of mitochondrial DNA is 15 502 bp and contains 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a A + T-rich region. The base composition of the genome is A (39.04%), T (37.07%), C (23.4%), and G (14.6%). Except for COI and ATP8 with TCC and ATC as start codon, respectively, the remaining protein-coding genes initiated with the three orthodox start codons. Two complete stop codons (TAA and TAG) and two incomplete stop codons (COIII stop with T and ND5 stop with TA) were used in the protein-coding genes. The A + T-rich region is located between 12s rRNA and tRNA(Ile) with the length of 859 bp. The phylogenetic relationships of Coleoptera species were constructed based on the nucleotide sequences of 13 protein-coding genes of mitogenome using the neighbor-joining method. The molecular-based phylogenetic analysis supported the traditional morphological classification on relationships within Coleoptera species.

Phylogenetic analysis of the complete mitogenome sequence of the raspberry weevil, Aegorhinus superciliosus (Coleoptera: Curculionidae), supports monophyly of the tribe Aterpini

The superfamily Curculionoidea is one of the most diverse groups of insects in the world, including many species which are crop pests. Within this group, the native raspberry weevil, Aegorhinus superciliosus (Guérin, 1830), is an important pest in blueberry and raspberry fields in southern South America. Using a 454 sequencing approach, we sequenced and annotated the mitogenome of A. superciliosus, it, providing the first such information for any species in the tribe Aterpini, subfamily Cyclominae. The assembled mitogenome is a circular DNA molecule 15,121bp in length containing all 37 genes normally found in metazoans. Mitogenome organization and transcriptional orientation in A. superciliosus showed the same pattern that characterizes the suborder Polyphaga. Bayesian and Maximum Likelihood phylogenetic analyses supported the monophyly of the tribe Aterpini and the subfamily Cyclominae, recovering this clade in a sister group relationship with Entiminae and Hyperinae. The monophyly of these three subfamilies defines a critical transition to an ectophagous lifestyle in the larvae, from an ancestrally endophagous larval lifestyle in all other lineages. The sequenced mitogenome of A. superciliosus can provide basic data for future studies investigating population history, molecular systematics, stress ecophysiology and phylogeography.

The mitochondrial genome of the multicolored Asian lady beetle Harmonia axyridis (Pallas) and a phylogenetic analysis of the Polyphaga (Insecta: Coleoptera)

Here, we report the mitochondrial genome sequence of the multicolored Asian lady beetle Harmonia axyridis (Pallas, 1773) (Coleoptera: Coccinellidae) (GenBank accession No. KR108208). This is the first species with sequenced mitochondrial genome from the genus Harmonia. The current length with partitial A + T-rich region of this mitochondrial genome is 16,387 bp. All the typical genes were sequenced except the trnI and trnQ. As in most other sequenced mitochondrial genomes of Coleoptera, there is no re-arrangement in the sequenced region compared with the pupative ancestral arrangement of insects. All protein-coding genes start with ATN codons. Five, five and three protein-coding genes stop with termination codon TAA, TA and T, respectively. Phylogenetic analysis using Bayesian method based on the first and second codon positions of the protein-coding genes supported that the Scirtidae is a basal lineage of Polyphaga. The Harmonia and the Coccinella form a sister lineage. The monophyly of Staphyliniformia, Scarabaeiformia and Cucujiformia was supported. The Buprestidae was found to be a sister group to the Bostrichiformia.

First record of Rhoptrocentruspiceus Marshall (Hymenoptera, Braconidae, Doryctinae) as parasitoid of Psacotheahilarishilaris (Pascoe) (Coleoptera, Cerambycidae)

The species Rhoptrocentruspiceus Marshall (Hymenoptera: Braconidae) was reared from the larvae of the xylophagous beetle Psacotheahilarishilaris (Pascoe) (Coleoptera: Cerambycidae), an exotic pest of Ficus and Morus species native to eastern Asia. It was recorded in the north of Italy in September 2005. This discovery is the first report of this species as parasitoids of the yellow spotted longicorn beetle all over the world.

The complete mitochondrial genome of the longhorn beetle Xylotrechus grayii (Coleoptera: Cerambycidae)

We sequenced the complete mitochondrial genome of the longhorn beetle, Xylotrechus grayii. The total length of the X. grayii mitogenome was 15,540 bp with an A + T content of 75.29%, consisting of 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes and an A + T-rich region. All the genes were arranged in the same order as that of the ancestral insect. All PCGs started with a typical ATN codon except for cox1 and nad1, which used TTG as start codon. Ten out of 13 PCGs terminated with incomplete codons (TA or T). The A + T-rich region was 893 bp in length with an A + T content of 85.89 %.

Complete mitochondrial genome of yellow meal worm (Tenebrio molitor)

The yellow meal worm (Tenebrio molitor L.) is an important resource insect typically used as animal feed additive. It is also widely used for biological research. The first complete mitochondrial genome of T. molitor was determined for the first time by long PCR and conserved primer walking approaches. The results showed that the entire mitogenome of T. molitor was 15 785 bp long, with 72.35% A+T content [deposited in GenBank with accession number KF418153]. The gene order and orientation were the same as the most common type suggested as ancestral for insects. Two protein-coding genes used atypical start codons (CTA in ND2 and AAT in COX1), and the remaining 11 protein-coding genes started with a typical insect initiation codon ATN. All tRNAs showed standard clover-leaf structure, except for tRNA(Ser) (AGN), which lacked a dihydrouridine (DHU) arm. The newly added T. molitor mitogenome could provide information for future studies on yellow meal worm.

Complete mitochondrial genome of the Japanese pine sawyer, Monochamus alternatus (Coleoptera: Cerambycidae)

We determined the complete mitochondrial genome of the Japanese pine sawyer Monochamus alternatus Hope (Coleoptera: Cerambycidae), which is a major forest pest in Asia. The genome is 15,874 bp in length containing 37 typical animal mitochondrial genes and one non-coding A+T-rich region. Its gene content and order are typical of other coleopteran mitochondrial genomes described to date. All protein-coding genes (PCGs) are initiated by ATN codons. Eight PCGs use complete stop codons TAG or TAA, whereas other PCGs end with a single T. All tRNA genes show typical secondary cloverleaf structures except for tRNA(Ser(AGN)), which lacks the dihydrouridine (DHU) arm. The large non-coding A+T-rich region of 1249 bp contains a 14 bp-long poly-T stretch and two microsatellite-like (AT)(TA)7 and (TA)8 elements.

The complete mitochondrial genome of Taeniogonalos taihorina (Bischoff) (Hymenoptera: Trigonalyidae) reveals a novel gene rearrangement pattern in the Hymenoptera

The family Trigonalyidae is considered to be one of the most basal lineages in the suborder Apocrita of Hymenoptera. Here, we determine the first complete mitochondrial genome of the Trigonalyidae, from the species Taeniogonalos taihorina (Bischoff, 1914). This mitochondrial genome is 15,927bp long, with a high A+T-content of 84.60%. It contains all of the 37 typical animal mitochondrial genes and an A+T-rich region. The orders and directions of all genes are different from those of previously reported hymenopteran mitochondrial genomes. Eight tRNA genes, three protein-coding genes and the A+T-rich region were rearranged, with the dominant gene rearrangement events being translocation and local inversion. The arrangements of three tRNA clusters, trnY-trnM-trnI-trnQ, trnW-trnL2-trnC, and trnH-trnA-trnR-trnN-trnS-trnE-trnF, and the position of the cox1 gene, are novel to the Hymenoptera, even the insects. Six long intergenic spacers are present in the genome. The secondary structures of the RNA genes are normal, except for trnS2, in which the D-stem pairing is absent.

The complete mitochondrial genome of the longhorn beetle, Massicus raddei

We sequenced the complete mitochondrial genome of Massicus raddei, which is the first beetle sequenced in Cerambycinae to date. The complete mitochondrial genome is 15,585 bp in length with an A + T content of 71.82%, and contains 13 protein-coding genes, 2 rRNAs, 22 tRNAs and a control region. The gene order and orientation are similar to that of typical insect species. These data will provide useful molecular information for phylogenetic relationships among the suborders of Coleoptera. By using 13 protein-coding genes as phylogenetic markers, the results support that the suborder Archostemata is a sister group to the remaining beetles and the most primitive suborder in any case; the suborder Myxophaga is sister to the suborder Adephaga.

Complete mitochondrial genome of the Burmese giant earthworm, Tonoscolex birmanicus (Clitellata: Megascolecidae)

Until now the complete mitochondrial genome (mitogenome) sequences of only three species of clitellate have been available. We have determined the complete mitogenome sequences of the elusive Burmese giant earthworm Tonoscolex birmanicus (Clitellata: Megascolecidae), which is endemic to Myanmar. The 15,170-bp long genome contains the 37 genes typical of metazoan mitogenomes [13 protein-coding genes (PCG), 2 rRNA genes and 22 tRNA genes] and 1 major non-coding region. All of the 37 genes are transcribed from the same DNA strand. The arrangement of the T. birmanicus mitogenome is identical to that of two within-ordinal species Lumbricus terrestris and Perionyx excavates. All 13 PCGs start with the ATG. For the stop codon, only six PCGs end with the TAA, whereas the remaining ones ends with the incomplete stop codon, T. Genes overlap in a total of 14 bp in five locations, and harbor a total of 16 bp of intergenic spacer sequences in nine locations.

Complete coding region of the mitochondrial genome of Monochamus alternatus hope (Coleoptera: Cerambycidae)

The Japanese pine sawyer, Monochamus alternatus Hope, 1842, an important forest pest, mainly occurs in Far East. It is the main vector of pine wood nematode Bursaphelenchus xylophilus, which causes pine wilt disease. We determined the complete mitochondrial genome coding region of M. alternatus using long PCR and conserved primer walking. Our results show that the entire mitogenome coding region is 14,649 bp long, with 78.22% A+T content [deposited in GenBank (JX987292)]. Positions and arrangement of the 37 genes encoded by the coding region are identical to those of two other longhorn beetles (Psacothea hilaris and Anoplophora glabripennis) for which the complete gene content and arrangement are known. All protein-coding genes start with a typical initiation codon ATN in insects. All tRNAs show standard clover-leaf structure, except the tRNA(Ser) (AGN), which lacks dihydrouridine (DHU) arm. The most unusual feature found is the use of TCT as tRNA(Ser) (AGN) anticodon instead of GCT, which is used in most other arthropods. This provides further insights into the diversity and evolution of the Cerambycidae family of long-horned beetles.

Characterization of the mitochondrial genome of the diamondback moth Plutella xylostella (Lepidoptera: Plutellidae) and phylogenetic analysis of advanced moths and butterflies

Here we determined the mitochondrial genome sequence of a notorious pest, the diamondback moth Plutella xylostella (Lepidoptera: Yponomeutoidea: Plutellidae). The mitochondrial genome contains 37 typical animal mitochondrial genes and an A+T-rich region. The gene arrangement is identical to that of other ditrysian lepidopteran mitochondrial genomes, but different from the ancestral gene arrangement in the non-ditrysian Hepialidae of Lepidoptera. The start codon of the cox1 gene is CGA, which is dissimilar to its homologs in most other insects. In Lepidoptera, cox1 and cox2 have low nucleotide diversities, while the nad6, nad2, and nad3 genes are highly variable. Phylogenetic analyses uncovered the reciprocal monophyly of Ditrysia, Apoditrysia, Obtectomera, and Macrolepidoptera, and the placement of the Hesperiidae within Papilionoidea. Our analyses suggest that the complete mitochondrial genome sequences are a promising marker toward fully resolving the phylogenetic relationships within Lepidoptera.

Complete mitochondrial genome of the seven-spotted lady beetle, Coccinella septempunctata (Coleoptera: Coccinellidae)

The complete mitochondrial genome (mitogenome) of the seven-spotted lady beetle, Coccinella septempunctata (Coleoptera: Coccinellidae), which is one of the best known insects capable of predation, is described with an emphasis on the noteworthy composition of the A+T-rich region. The C. septempunctata genome consists of 2 rRNAs, 22 tRNAs, 13 protein-coding genes, and 1 control region, designated as the A+T-rich region in insects. Along with an unusually long A+T-rich region (4469 bp), the 18,965-bp long C. septempunctata mitogenome was the largest in Coleoptera. The A+T-rich region is composed of a 2214-bp long non-repeat region composed of 80.17% A/T nucleotides and a 2256-bp long repeat region composed of 65.71% A/T nucleotides. The repeat region harbors 32 identical 70-bp long tandem repeats plus one 15-bp long incomplete first repeat. These repeat sequences may possibly have been caused by slipped-strand mispairing and unequal crossing-over events during DNA replication.

The first mitochondrial genome for the wasp superfamily Platygastroidea: the egg parasitoid Trissolcus basalis

The nearly complete mitochondrial (mt) genome of an egg parasitoid, Trissolcus basalis (Wollaston), was sequenced using both 454 and Illumina next-generation sequencing technologies. A portion of the noncoding region remained unsequenced, possibly owing to the presence of repeats. The sequenced portion of the genome is 15,768 bp and has a high A+T content (84.2%), as is typical for hymenopteran mt genomes. A total of 36 of the 37 genes normally present in animal mt genomes were located. The one exception was trnR; a truncated version of this gene is present between trnS(1) and nd5, but it is unclear whether this gene fragment could code for the entire trnR gene. The mt gene arrangement of T. basalis is different from other Proctotrupomorpha mt genomes, with a number of trn genes in different positions. However, no shared derived gene rearrangements were identified in the present study. Bayesian analyses of mt genomes from 29 hymenopteran taxa and seven other orders of holometabolous insects support some uncontroversial evolutionary relationships, but indicate that much higher levels of taxonomic sampling are necessary for the resolution of family and superfamily relationships.

The mitochondrial genome of the ascalaphid owlfly Libelloides macaronius and comparative evolutionary mitochondriomics of neuropterid insects

Background

The insect order Neuroptera encompasses more than 5,700 described species. To date, only three neuropteran mitochondrial genomes have been fully and one partly sequenced. Current knowledge on neuropteran mitochondrial genomes is limited, and new data are strongly required. In the present work, the mitochondrial genome of the ascalaphid owlfly Libelloides macaronius is described and compared with the known neuropterid mitochondrial genomes: Megaloptera, Neuroptera and Raphidioptera. These analyses are further extended to other endopterygotan orders.

Results

The mitochondrial genome of L. macaronius is a circular molecule 15,890 bp long. It includes the entire set of 37 genes usually present in animal mitochondrial genomes. The gene order of this newly sequenced genome is unique among Neuroptera and differs from the ancestral type of insects in the translocation of trnC. The L. macaronius genome shows the lowest A+T content (74.50%) among known neuropterid genomes. Protein-coding genes possess the typical mitochondrial start codons, except for cox1, which has an unusual ACG. Comparisons among endopterygotan mitochondrial genomes showed that A+T content and AT/GC-skews exhibit a broad range of variation among 84 analyzed taxa. Comparative analyses showed that neuropterid mitochondrial protein-coding genes experienced complex evolutionary histories, involving features ranging from codon usage to rate of substitution, that make them potential markers for population genetics/phylogenetics studies at different taxonomic ranks. The 22 tRNAs show variable substitution patterns in Neuropterida, with higher sequence conservation in genes located on the α strand. Inferred secondary structures for neuropterid rrnS and rrnL genes largely agree with those known for other insects. For the first time, a model is provided for domain I of an insect rrnL. The control region in Neuropterida, as in other insects, is fast-evolving genomic region, characterized by AT-rich motifs.

Conclusions

The new genome shares many features with known neuropteran genomes but differs in its low A+T content. Comparative analysis of neuropterid mitochondrial genes showed that they experienced distinct evolutionary patterns. Both tRNA families and ribosomal RNAs show composite substitution pathways. The neuropterid mitochondrial genome is characterized by a complex evolutionary history.

Conserved PCR primer set designing for closely-related species to complete mitochondrial genome sequencing using a sliding window-based PSO algorithm

Background

Complete mitochondrial (mt) genome sequencing is becoming increasingly common for phylogenetic reconstruction and as a model for genome evolution. For long template sequencing, i.e., like the entire mtDNA, it is essential to design primers for Polymerase Chain Reaction (PCR) amplicons which are partly overlapping each other. The presented chromosome walking strategy provides the overlapping design to solve the problem for unreliable sequencing data at the 5′ end and provides the effective sequencing. However, current algorithms and tools are mostly focused on the primer design for a local region in the genomic sequence. Accordingly, it is still challenging to provide the primer sets for the entire mtDNA.

Methodology/Principal Findings

The purpose of this study is to develop an integrated primer design algorithm for entire mt genome in general, and for the common primer sets for closely-related species in particular. We introduce ClustalW to generate the multiple sequence alignment needed to find the conserved sequences in closely-related species. These conserved sequences are suitable for designing the common primers for the entire mtDNA. Using a heuristic algorithm particle swarm optimization (PSO), all the designed primers were computationally validated to fit the common primer design constraints, such as the melting temperature, primer length and GC content, PCR product length, secondary structure, specificity, and terminal limitation. The overlap requirement for PCR amplicons in the entire mtDNA is satisfied by defining the overlapping region with the sliding window technology. Finally, primer sets were designed within the overlapping region. The primer sets for the entire mtDNA sequences were successfully demonstrated in the example of two closely-related fish species. The pseudo code for the primer design algorithm is provided.

Conclusions/Significance

In conclusion, it can be said that our proposed sliding window-based PSO algorithm provides the necessary primer sets for the entire mt genome amplification and sequencing.

Why barcode? High-throughput multiplex sequencing of mitochondrial genomes for molecular systematics

Mitochondrial genome sequences are important markers for phylogenetics but taxon sampling remains sporadic because of the great effort and cost required to acquire full-length sequences. Here, we demonstrate a simple, cost-effective way to sequence the full complement of protein coding mitochondrial genes from pooled samples using the 454/Roche platform. Multiplexing was achieved without the need for expensive indexing tags ('barcodes'). The method was trialled with a set of long-range polymerase chain reaction (PCR) fragments from 30 species of Coleoptera (beetles) sequenced in a 1/16th sector of a sequencing plate. Long contigs were produced from the pooled sequences with sequencing depths ranging from ∼10 to 100× per contig. Species identity of individual contigs was established via three 'bait' sequences matching disparate parts of the mitochondrial genome obtained by conventional PCR and Sanger sequencing. This proved that assembly of contigs from the sequencing pool was correct. Our study produced sequences for 21 nearly complete and seven partial sets of protein coding mitochondrial genes. Combined with existing sequences for 25 taxa, an improved estimate of basal relationships in Coleoptera was obtained. The procedure could be employed routinely for mitochondrial genome sequencing at the species level, to provide improved species 'barcodes' that currently use the cox1 gene only.

Complete mitochondrial genome of brown marmorated stink bug Halyomorpha halys (Hemiptera: Pentatomidae), and phylogenetic relationships of hemipteran suborders

The newly sequenced complete mitochondrial genome of the brown marmorated stink bug, Halyomorpha halys (Stål) (Hemiptera: Pentatomidae), is a circular molecule of 16,518 bp with a total A+T content of 76.4% and two extensive repeat regions in A+T rich region. Nucleotide composition and codon usage of H. halys are about average when compared with values observed in 19 other published hemipteran mitochondrial genomes. Phylogenetic analyses using these 20 hemipteran mitochondrial genomes support the currently accepted hypothesis that suborders Heteroptera and Auchenorrhyncha form a monophyletic group. The mitochondrial gene arrangements of the 20 genomes=are also consistent with our results.