The Evolution of Dark Matter in the Mitogenome of Seed Beetles

The first complete mitochondrial genome of a Lethrus species (Coleoptera, Geotrupidae) with phylogenetic implications

The flightless beetle genus Lethrus Scopoli, 1777 (Geotrupidae, Scarabaeoidea) has a large distribution area throughout Eurasia and is characterized by many species, especially in Middle Asia. Despite this diversity and the potential importance as models for speciation, Lethrus species are underrepresented in molecular databases. To fill this gap, we report the complete mitochondrial genome of Lethrusscoparius obtained using third-generation sequencing technology. The circular mitogenome is 24,944 bp long and has a structure characteristic of coleopterans. It contains 37 genes, including 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and an A+T-rich non-coding region, the control region between the 12S rRNA and tRNA-Ile (GAU). The phylogenetic analysis of the superfamily Scarabaeoidea placed L.scoparius in the monophyletic family Geotrupidae, which is related to the family Scarabaeidae. The assembled mitochondrial genome is a valuable new genomic resource in the genus Lethrus and contributes to a better understanding of the evolutionary history of the genus and the entire family Geotrupidae.

Toward the Development of the Trojan Female Technique in Pest Insects: Male-Specific Influence of Mitochondrial Haplotype on Reproductive Output in the Seed Beetle Acanthoscelides obtectus

Biocontrol techniques that impair reproductive capacity of insect pests provide opportunities to control the dynamics of their populations while minimizing collateral damage to non-target species and the environment. The Trojan Female Technique, or TFT, is a method of the trans-generational fertility-based population control through the release of females that carry mitochondrial DNA mutations that negatively affect male, but not female, reproductive output. TFT is based on the evolutionary hypothesis that, due to maternal inheritance of mitochondria, mutations which are beneficial or neutral in females but harmful in males can accumulate in the mitochondrial genome without selection acting against them. Although TFT has been theoretically substantiated, empirical work to date has focused only on Drosophila melanogaster populations, while the existence of male-biased mutations and the TFT approach in economically important pest species remain unexplored. Here, we examined the sex-specific effects of three distinct and naturally occurring mitochondrial haplotypes (MG1a, MG1d, and MG3b) on several reproductive and life history traits in the seed beetle Acanthoscelides obtectus. Our results revealed that males harboring the MG3b mitotype exhibited lower early fecundity and fertility, while there were no effects on females or longevity in either sex. Our experiments provide support for the existence of the mitochondrial variant that specifically impairs male reproductive output in pest insects. These results can be harnessed to further develop TFT as a novel form of biocontrol with broad applicability to economic pests and disease vector insects.

Comparative mitochondrial genomics of Thelebolaceae in Antarctica: insights into their extremophilic adaptations and evolutionary dynamics

Species of Antarctomyces and Thelebolus (Thelebolaceae), primarily found in Antarctic environments, exhibit psychrophilic adaptations, yet their mitochondrial genomes have not been extensively studied. Furthermore, few studies have compared the mitochondrial genomes of psychrophilic, psychrotrophic, and mesophilic fungi. After successful sequencing and assembly, this study annotated the mitochondrial genomes of Antarctomyces psychrotrophicus CPCC 401038 and Thelebolus microsporus CPCC 401041. We also performed a comparative analysis with the previously characterized mitochondrial genomes of psychrotrophic and mesophilic fungi. The analysis revealed that nad4L was the most conserved gene across the mitochondrial genomes, characterized by its synonymous and non-synonymous substitution rates (Ks and Ka), genetic distance, and GC content and skew within the protein-coding genes (PCGs). Additionally, the mitochondrial genomes of psychrophilic and psychrotrophic fungi showed a higher proportion of protein-coding regions and a lower GC content compared to those of mesophilic fungi, underscoring the genetic basis of cold adaptation. Phylogenetic analyses based on these mitochondrial genes also confirmed the phylogenetic relationships of Thelebolaceae in the class Leotiomycetes. These findings advance our understanding of the phylogenetic relationships and evolutionary dynamics within the family Thelebolaceae, highlighting how different environmental temperatures influence fungal mitochondrial genomic structure and adaptation.

New taxonomic insights for Brazilian Syrbatus Reitter (Coleoptera: Staphylinidae: Pselaphinae), including three new species and their mitochondrial genomes

Here we present a taxonomic treatment for the Brazilian species of Syrbatus (Reitter, 1882), including the description of three new species (Syrbatus moustache Asenjo & Valois sp. nov., Syrbatus obsidian Asenjo & Valois sp. nov. and Syrbatus superciliata Asenjo & Valois sp. nov.) from the Quadrilátero Ferrífero (Minas Gerais, Brazil). In addition, we designated lectotypes for the Brazilian species of species-group 2, Syrbatus centralis (Raffray, 1898), Syrbatus hetschkoi (Reitter, 1888), Syrbatus hiatusus (Reitter, 1888), Syrbatus transversalis (Raffray, 1898), and Syrbatus trinodulus (Schaufuss, 1887), besides recognizing the holotype for Syrbatus brevispinus (Reitter, 1882), Syrbatus bubalus (Raffray, 1898), and Syrbatus grouvellei (Raffray, 1898). The mitochondrial genomes (mitogenomes) of the three new species are presented, for which we present the phylogenetic placement among Staphylinidae with previously published data.

The complete mitochondrial genome of Ceresium sinicum ornaticolle Pic, 1907 (Coleoptera: Cerambycidae)

Ceresium sinicum ornaticolle Pic, 1907 (Coleoptera: Cerambycidae) is one of the main pests of pomegranate and citrus trees. In this study, we described the complete mitochondrial genome of C. sinicum ornaticolle. The total length of the mitochondrial genome was 15,817 bp, and the entire content of GC was 27.8%. The genome encoded 2 ribosomal RNA genes (rRNAs), 13 protein-coding genes (PCGs), and 22 transfer RNA genes (tRNAs). The phylogenetic tree showed that C. sinicum ornaticolle was clustered with Allotraeus orientalis and Zoodes fulguratus. These results will provide the genetic information for understanding the genetic evolution of C. sinicum ornaticolle and the insights to control cerambycid pests.

A chromosome-level assembly of the seed beetle Callosobruchus maculatus genome with annotation of its repetitive elements

Callosobruchus maculatus is a major agricultural pest of legume crops worldwide and an established model system in ecology and evolution. Yet, current molecular biological resources for this species are limited. Here, we employ Hi-C sequencing to generate a greatly improved genome assembly and we annotate its repetitive elements in a dedicated in-depth effort where we manually curate and classify the most abundant unclassified repeat subfamilies. We present a scaffolded chromosome-level assembly, which is 1.01 Gb in total length with 86% being contained within the 9 autosomes and the X chromosome. Repetitive sequences accounted for 70% of the total assembly. DNA transposons covered 18% of the genome, with the most abundant superfamily being Tc1-Mariner (9.75% of the genome). This new chromosome-level genome assembly of C. maculatus will enable future genetic and evolutionary studies not only of this important species but of beetles more generally.

Whole-genome sequence of the Cooley spruce gall adelgid, Adelges cooleyi (Hemiptera: Sternorrhyncha: Adelgidae)

The adelgids (Adelgidae) are a small family of sap-feeding insects, which, together with true aphids (Aphididae) and phylloxerans (Phylloxeridae), make up the infraorder Aphidomorpha. Some adelgid species are highly destructive to forest ecosystems such as Adelges tsugae, Adelges piceae, Adelges laricis, Pineus pini, and Pineus boerneri. Despite this, there are no high-quality genomic resources for adelgids, hindering advanced genomic analyses within Adelgidae and among Aphidomorpha. Here, we used PacBio continuous long-read and Illumina RNA-sequencing to construct a high-quality draft genome assembly for the Cooley spruce gall adelgid, Adelges cooleyi (Gillette), a gall-forming species endemic to North America. The assembled genome is 270.2 Mb in total size and has scaffold and contig N50 statistics of 14.87 and 7.18 Mb, respectively. There are 24,967 predicted coding sequences, and the assembly completeness is estimated at 98.1 and 99.6% with core BUSCO gene sets of Arthropoda and Hemiptera, respectively. Phylogenomic analysis using the A. cooleyi genome, 3 publicly available adelgid transcriptomes, 4 phylloxera transcriptomes, the Daktulosphaira vitifoliae (grape phylloxera) genome, 4 aphid genomes, and 2 outgroup coccoid genomes fully resolves adelgids and phylloxerans as sister taxa. The mitochondrial genome is 24 kb, among the largest in insects sampled to date, with 39.4% composed of noncoding regions. This genome assembly is currently the only genome-scale, annotated assembly for adelgids and will be a valuable resource for understanding the ecology and evolution of Aphidomorpha.

The complete mitochondrial genome and description of a new cryptic Brazilian species of Metopiellus Raffray (Coleoptera: Staphylinidae: Pselaphinae)

Metopiellus Raffray, 1908 is a genus of South American rove beetles typically found in tropical humid forests. Here we describe a new cryptic species from Eastern Amazon, in northern Brazil, Metopiellus crypticus Asenjo sp. nov., and its major morphologic diagnostic features, which were photographed and illustrated. In addition, we bring the complete mitochondrial genome sequence of M. crypticus sp. nov., and its position within the phylogenetic context of the family, including previously available mitogenomes of Staphylinidae species.

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

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

New Dielis species and structural dichotomy of the mitochondrial cox2 gene in Scoliidae wasps

Some mitochondrial protein-coding genes of protists and land plants have split over the course of evolution into complementary genes whose products can form heteromeric complexes that likely substitute for the undivided proteins. One of these genes, cox2, has also been found to have split in animals, specifically in Scoliidae wasps (Hymenoptera: Apocrita) of the genus Dielis (Campsomerini), while maintaining the conventional structure in related Scolia (Scoliini). Here, a hitherto unrecognized Nearctic species of Dielis, D. tejensis, is described based on its phenotype and mtDNA. The mitogenome of D. tejensis sp. nov. differs from that of the sympatric sibling species Dielis plumipes fossulana by the reduced size of the cox2-dividing insert, which, however, still constitutes the fifth part of the mtDNA; an enlarged nad2-trnW intergenic region; the presence of two trnK^ttt paralogues; and other features. Both species of Dielis have a unique insertion of a threonine in COXIIA, predicted to be involved in COXIIA-COXIIB docking, and substitutions of two hydrophobic residues with redox-active cysteines around the Cu_A centre in COXIIB. Importantly, the analysis of mtDNA from another Campsomerini genus, Megacampsomeris, shows that its cox2 gene is also split. The presented data highlight evolutionary processes taking place in hymenopteran mitogenomes that do not fall within the mainstream of animal mitochondrion evolution.

Complete mitochondrial genome of Episymploce splendens (Blattodea: Ectobiidae): A large intergenic spacer and lacking of two tRNA genes

The complete mitochondrial genome of Episymploce splendens, 15,802 bp in length, was determined and annotated in this study. The mito-genome included 13 PCGs, 20 tRNAs and 2 rRNAs. Unlike most typical mito-genomes with conservative gene arrangement and exceptional economic organization, E. splendens mito-genome has two tRNAs (tRNA-Gln and tRNA-Met) absence and a long intergenic spacer sequence (93 bp) between tRNA-Val and srRNA, showing the diversified features of insect mito-genomes. This is the first report of the tRNAs deletion in blattarian mito-genomes and we supported the duplication/random loss model as the origin mechanism of the long intergenic spacer. Two Numts, Numt-1 (557 bp) and Numt-2 (975 bp) transferred to the nucleus at about 14.15 Ma to 22.34 Ma, and 19.19 Ma to 24.06 Ma respectively, were found in E. splendens. They can be used as molecular fossils in insect phylogenetic relationship inference. Our study provided useful data for further studies on the evolution of insect mito-genome.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

The Mitochondrial Genomes of 18 New Pleurosticti (Coleoptera: Scarabaeidae) Exhibit a Novel trnQ-NCR-trnI-trnM Gene Rearrangement and Clarify Phylogenetic Relationships of Subfamilies within Scarabaeidae

The availability of next-generation sequencing (NGS) in recent years has facilitated a revolution in the availability of mitochondrial (mt) genome sequences. The mt genome is a powerful tool for comparative studies and resolving the phylogenetic relationships among insect lineages. The mt genomes of phytophagous scarabs of the subfamilies Cetoniinae and Dynastinae were under-represented in GenBank. Previous research found that the subfamily Rutelinae was recovered as a paraphyletic group because the few representatives of the subfamily Dynastinae clustered into Rutelinae, but the subfamily position of Dynastinae was still unclear. In the present study, we sequenced 18 mt genomes from Dynastinae and Cetoniinae using next-generation sequencing (NGS) to re-assess the phylogenetic relationships within Scarabaeidae. All sequenced mt genomes contained 37 sets of genes (13 protein-coding genes, 22 tRNAs, and two ribosomal RNAs), with one long control region, but the gene order was not the same between Cetoniinae and Dynastinae species. All mt genomes of Dynastinae species showed the same gene rearrangement of trnQ-NCR-trnI-trnM, whereas all mt genomes of Cetoniinae species showed the ancestral insect gene order of trnI-trnQ-trnM. Phylogenetic analyses (IQ-tree and MrBayes) were conducted using 13 protein-coding genes based on nucleotide and amino acid datasets. In the ML and BI trees, we recovered the monophyly of Rutelinae, Cetoniinae, Dynastinae, and Sericinae, and the non-monophyly of Melolonthinae. Cetoniinae was shown to be a sister clade to (Dynastinae + Rutelinae).

The First Eight Mitogenomes of Leaf-Mining Dactylispa Beetles (Coleoptera: Chrysomelidae: Cassidinae) Shed New Light on Subgenus Relationships

The taxonomic classification of Dactylispa, a large genus of leaf-mining beetles, is problematic because it is currently based on morphology alone. Here, the first eight mitochondrial genomes of Dactylispa species, which were used to construct the first molecular phylogenies of this genus, are reported. The lengths of the eight mitogenomes range from 17,189 bp to 20,363 bp. All of the mitochondrial genomes include 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and 1 A + T-rich region. According to the nonsynonymous/synonymous mutation ratio (Ka/Ks) of all PCGs, the highest and the lowest evolutionary rates were found for atp8 and cox1, respectively, which is a common phenomenon among animals. According to relative synonymous codon usage, UUA(L) has the highest frequency. With two Gonophorini species as the outgroup, mitogenome-based phylogenetic trees of the eight Dactylispa species were constructed using maximum likelihood (ML) and Bayesian inference (BI) methods based on the PCGs, tRNAs, and rRNAs. Two DNA-based phylogenomic inferences and one protein-based phylogenomic inference support the delimitation of the subgenera Dactylispa s. str. and Platypriella as proposed in the system of Chen et al. (1986). However, the subgenus Triplispa is not recovered as monophyletic. The placement of Triplispa species requires further verification and testing with more species. We also found that both adult body shape and host plant relationship might explain the subgeneric relationships among Dactylispa beetles to a certain degree.

The Evidential Statistics of Genetic Assembly: Bootstrapping a Reference Sequence

The reference sequences play an essential role in genome assembly, like type specimens in taxonomy. Those references are also samples obtained at some time and location with a specific method. How can we evaluate or discriminate uncertainties of the reference itself and assembly methods? Here we bootstrapped 50 random read data sets from a small circular genome of a Escherichia coli bacteriophage, phiX174, and tried to reconstruct the reference with 14 free assembly programs. Nine out of 14 assembly programs were capable of circular genome reconstruction. Unicycler correctly reconstructed the reference for 44 out of 50 data sets, but each reconstructed contig of the failed six data sets had minor defects. The other assembly software could reconstruct the reference with minor defects. The defect regions differed among the assembly programs, and the defect locations were far from randomly distributed in the reference genome. All contigs of Trinity included one, but Minia had two perfect copies other than an imperfect reference copy. The centroid of contigs for assembly programs except Unicycler differed from the reference with 75bases at most. Nonmetric multidimensional scaling (NMDS) plots of the centroids indicated that even the reference sequence was located slightly off from the estimated location of the true reference. We propose that the combination of bootstrapping a reference, making consensus contigs as centroids in an edit distance, and NMDS plotting will provide an evidential statistic way of genetic assembly for non-fragmented base sequences.

The complete mitochondrial genome sequence of Oryctes rhinoceros (Coleoptera: Scarabaeidae) based on long-read nanopore sequencing

BACKGROUND

The coconut rhinoceros beetle (CRB, Oryctes rhinoceros) is a severe and invasive pest of coconut and other palms throughout Asia and the Pacific. The biocontrol agent, Oryctes rhinoceros nudivirus (OrNV), has successfully suppressed O. rhinoceros populations for decades but new CRB invasions started appearing after 2007. A single-SNP variant within the mitochondrial cox1 gene is used to distinguish the recently-invading CRB-G lineage from other haplotypes, but the lack of mitogenome sequence for this species hinders further development of a molecular toolset for biosecurity and management programmes against CRB. Here we report the complete circular sequence and annotation for CRB mitogenome, generated to support such efforts.

METHODS

Sequencing data were generated using long-read Nanopore technology from genomic DNA isolated from a CRB-G female. The mitogenome was assembled with Flye v.2.5, using the short-read Illumina sequences to remove homopolymers with Pilon, and annotated with MITOS. Independently-generated transcriptome data were used to assess the O. rhinoceros mitogenome annotation and transcription. The aligned sequences of 13 protein-coding genes (PCGs) (with degenerate third codon position) from O. rhinoceros, 13 other Scarabaeidae taxa and two outgroup taxa were used for the phylogenetic reconstruction with the Maximum likelihood (ML) approach in IQ-TREE and Bayesian (BI) approach in MrBayes.

RESULTS

The complete circular mitogenome of O. rhinoceros is 20,898 bp in length, with a gene content canonical for insects (13 PCGs, two rRNA genes, and 22 tRNA genes), as well as one structural variation (rearrangement of trnQ and trnI) and a long control region (6,204 bp). Transcription was detected across all 37 genes, and interestingly, within three domains in the control region. ML and BI phylogenies had the same topology, correctly grouping O. rhinoceros with one other Dynastinae taxon, and recovering the previously reported relationship among lineages in the Scarabaeidae. In silico PCR-RFLP analysis recovered the correct fragment set that is diagnostic for the CRB-G haplogroup. These results validate the high-quality of the O. rhinoceros mitogenome sequence and annotation.

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.

An experimental test of temperature-dependent selection on mitochondrial haplotypes in Callosobruchus maculatus seed beetles

Mitochondrial DNA (mtDNA) consists of few but vital maternally inherited genes that interact closely with nuclear genes to produce cellular energy. How important mtDNA polymorphism is for adaptation is still unclear. The assumption in population genetic studies is often that segregating mtDNA variation is selectively neutral. This contrasts with empirical observations of mtDNA haplotypes affecting fitness-related traits and thermal sensitivity, and latitudinal clines in mtDNA haplotype frequencies. Here, we experimentally test whether ambient temperature affects selection on mtDNA variation, and whether such thermal effects are influenced by intergenomic epistasis due to interactions between mitochondrial and nuclear genes, using replicated experimental evolution in Callosobruchus maculatus seed beetle populations seeded with a mixture of different mtDNA haplotypes. We also test for sex-specific consequences of mtDNA evolution on reproductive success, given that mtDNA mutations can have sexually antagonistic fitness effects. Our results demonstrate natural selection on mtDNA haplotypes, with some support for thermal environment influencing mtDNA evolution through mitonuclear epistasis. The changes in male and female reproductive fitness were both aligned with changes in mtDNA haplotype frequencies, suggesting that natural selection on mtDNA is sexually concordant in stressful thermal environments. We discuss the implications of our findings for the evolution of mtDNA.

Characterization of the complete mitochondrial genome of Caryopemon giganteus Pic (Coleoptera: Chrysomelidae: Bruchinae)

We sequenced, assembled, and annotated the complete mitochondrial genome of the seed beetle Caryopemon giganteus, which represents the first report in the tribe Caryopemini from the subfamily Bruchinae of Chrysomelidae. The circular mitochondrial genome of the species contains 15,727 bases, 13 protein-coding genes (PCGs), 22 tRNA genes, 2 rRNA genes, and a non-coding region. The GC content of the genome is 25.3%, which is higher than any other reported mitochondrial genomes within Bruchinae. The 16S ribosomal RNA gene and the 12S ribosomal RNA gene are 1284 and 835 bp in length, respectively. 12 PCGs started with the typical ATN codon, except for ND1 initiated with TTG. Five PCGs have the typical stop codon of TAA or TGA, while the remainder PCGs are terminated with incomplete stop codons (TA or T). The phylogenetic analysis based on a combination of 13 genes of the mitochondrial genomes of six species of Bruchinae and 23 species from other 10 subfamilies of Chrysomelidae recovered a generally well resolved and strongly supported tree topology, which shows that C. giganteus has the basalmost position in Bruchinae.

Description of soybean aphid (Aphis glycines Matsumura) mitochondrial genome and comparative mitogenomics of Aphididae (Hemiptera: Sternorrhyncha)

The complete mitochondrial genome of the soybean aphid (Aphis glycines Matsumura), a major agricultural pest in the world, is described for the first time, which consists of 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, as well as a large repeat region between tRNA-Glu and tRNA-Phe, and an AT-rich control region. The 17,954 bp mtgenome is the largest reported from the family Aphididae, and its gene order follows the ancestral insect mtgenome except for the repeat region, which contains a 195 bp unit repeated 11.9 times, representing the highest reported repeats among the known aphid mtgenomes to date. A new molecular phylogeny of Aphidae is reconstructed based on all available aphid mtgenomes, and it is shown that the mtgenome data can robustly resolve relationships at the subfamily level, but do not have sufficient phylogenetic information to resolve deep relationships. A phylogeny-based comparative analysis of mtgenomes has been performed to investigate the evolution of the repeat region between tRNA-Glu and tRNA-Phe. So far, among aphids, 13 species are known to have this repeat region of variable lengths, and a phylogenetic analysis of the repeat region shows that a large proportion of the sequences are conserved across the phylogeny, suggesting that the repeat region evolved in the most recent common ancestor of Aphidinae and Eriosomatinae, and that it has gone through numerous episodes of lineage-specific losses and expansions. Combined together, this study provides novel insights into how the repeat regions have evolved within aphids.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

Mitochondrial genome of a leaf-mining beetle Prionispa champaka Maulik (Coleoptera: Chrysomelidae: Cassidinae)

Prionispa champaka is a leaf-mining species which feeds on Pollia spp. and widely distributes in southern China. The complete mitogenomic sequence of P. champaka (Chrysomelidae: Cassidinae) was obtained and annotated, with a length of 20,494 bp. It was longer than those of other Chrysomelid species (not including Bruchinae) because of its much longer non-coding sequences. Gene arrangement and content of P. champaka was identical to the most common type in insects, and it was also biased toward AT (accounting for 78.4%). Phylogenetic analysis based on mitochondrial PCGs indicated that P. champaka was closely clustered with 5 other Cassidinae species, supporting the traditional morphological classification within Cassidinae.