Mitochondrial genomes of three Bostrichiformia species and phylogenetic analysis of Polyphaga (Insecta, Coleoptera).

Here we determined mitogenomes of three Bostrichiformia species. These data were combined with 51 previously sequenced Polyphaga mitogenomes to explore the higher-level relationships within Polyphaga by using four different mitogenomic datasets and three tree inference approaches. Among Polyphaga mitogenomes we observed heterogeneity in nucleotide composition and evolutionary rates, which may have affected phylogenetic inferences across the different mitogenomic datasets. Elateriformia, Cucujiformia, and Scarabaeiformia were each inferred to be monophyletic by all analyses, as was Bostrichiformia by most analyses based on two datasets with low heterogeneity. The large series Staphyliniformia was never recovered as monophyletic in our analyses. The Bayesian tree using a degenerated nucleotide dataset (P123_Degen) and a site-heterogeneous mixture model in PhyloBayes was supported as the best Polyphaga phylogeny: (Scirtiformia, (Elateriformia, ((Bostrichiformia, Cucujiformia), (Scarabaeiformia + Staphyliniformia)))). For Cucujiformia, the largest series, we inferred a superfamily-level phylogeny: ((Cleroidea, Coccinelloidea), (Tenebrionoidea, (Cucujoidea + Curculionoidea + Chrysomeloidea))).

Mitochondrial phylogeny, divergence history and high-altitude adaptation of grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) inhabiting the Tibetan Plateau.

Grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) are the most important pests in alpine meadows of the Tibetan Plateau (TP) and have well adapted to high-altitude environments. To further understand the evolutionary history and their adaptation to the TP, we newly determined seven complete TP Gynaephora mitogenomes. Compared to single genes, whole mitogenomes provided the best phylogenetic signals and obtained robust results, supporting the monophyly of the TP Gynaephora species and a phylogeny of Arctiinae + (Aganainae + Lymantriinae). Incongruent phylogenetic signals were found among single mitochondrial genes, none of which recovered the same phylogeny as the whole mitogenome. We identified six best-performing single genes using Shimodaira-Hasegawa tests and found that the combinations of rrnS and either cox1 or cox3 generated the same phylogeny as the whole mitogenome, indicating the phylogenetic potential of these three genes for future evolutionary studies of Gynaephora. The TP Gynaephora species were estimated to radiate on the TP during the Pliocene and Quaternary, supporting an association of the diversification and speciation of the TP Gynaephora species with the TP uplifts and associated climate changes during this time. Selection analyses revealed accelerated evolutionary rates of the mitochondrial protein-coding genes in the TP Gynaephora species, suggesting that they accumulated more nonsynonymous substitutions that may benefit their adaptation to high altitudes. Furthermore, signals of positive selection were detected in nad5 of two Gynaephora species with the highest altitude-distributions, indicating that this gene may contribute to Gynaephora's adaptation to divergent altitudes. This study adds to the understanding of the TP Gynaephora evolutionary relationships and suggests a link between mitogenome evolution and ecological adaptation to high-altitude environments in grassland caterpillars.

Mitochondrial genome of Chrysochares punctatus (Coleoptera: Chrysomelidae: Eumolpinae) and phylogenetic analysis.

Here, we determined the nearly complete mitochondrial genome (mitogenome) of Chrysochares punctatus (Coleoptera: Chrysomelidae: Eumolpinae), an important insect pest on Apocynum venetum in Northwestern China. This mitogenome was 14,451 bp long, encoding 13 protein-coding genes (PCGs), 21 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes. The C. punctatus mitogenome presented an A + T content of 75.11%, with a positive AT-skew (0.064) and a negative GC-skew (-0.192). Ten PCGs started with a typical ATN codon, whereas the remaining three PCGs started with AAC (cox1) and TTG (nad1 and nad2). All tRNAs had a typical secondary cloverleaf structure, except for trnS1 which lacked the dihydrouridine arm. Bayesian phylogenetic analysis based on the nucleotide sequences of 13 PCGs recovered a phylogeny within Chrysomelidae: (((Chrysomelinae + Galerucinae), (((Eumolpinae, Lamprosomatinae), Cassidinae), Criocerinae)), Bruchinae).

Mitogenome evolution in ladybirds: Potential association with dietary adaptation.

Dietary shifts can alter the relative availability of different nutrients and are therefore associated with metabolic adaptation in animals. The Coccinellidae (ladybirds) exhibits three major types of feeding habits and provides a useful model to study the effects of dietary changes on the evolution of mitogenomes, which encode proteins directly involved in energy metabolism. Here, mitogenomes of three coccinellid species were newly sequenced. These data were combined with other ten previously sequenced coccinellid mitogenomes to explore the relationship between mitogenome evolution and diets. Our results indicate that mitogenomic data can be effectively used to resolve phylogenetic relationships of Coccinellidae. Strong codon usage bias in coccinellid mitogenomes was predominantly determined by nucleotide composition. The 13 mitochondrial protein-coding genes (PCGs) globally evolved under negative constraints, with some PCGs showing a stronger purifying selection. Six PCGs (nad3, nad4L, and nad5 from Complex I; cox1 and cox3 from Complex IV; and atp6 from Complex V) displayed signs of positive selection. Of these, adaptive changes in cox3 were potentially associated with metabolic differences resulting from dietary shifts in Coccinellidae. Our results provide insights into the adaptive evolution of coccinellid mitogenomes in response to both dietary shifts and other life history traits.

Adaptive evidence of mitochondrial genomes in Dolycoris baccarum (Hemiptera: Pentatomidae) to divergent altitude environments.

Given mitochondrion is the 'energy and oxygen usage factories', adaptive signatures of mitochondrial genes have been extensively investigated in vertebrates from different altitudes, but few studies focus on insects. Here, we sequenced the complete mitochondrial genome (mitogenome) of Dolycoris. baccarum living in the Tibetan Plateau (DBHC, ∼3200 m above sea level (asl)) and conducted a detailed comparative analysis with another D. baccarum mitogenome (DBQY) from relatively low altitude (∼1300 m asl). All the 37 mitochondrial genes were highly conserved and under purifying selection, except for two mitochondrial protein-coding genes (MPCGs) (atp6 and nad5) that showed positively selected signatures. We therefore further examined non-synonymous substitutions in atp6 and nad5, by sequencing more individuals from three populations with different altitudes. We found that these non-synonymous substitutions were polymorphic in these populations, likely due to relaxed selection constraints in different altitudes. Purifying selection in all mitochondrial genes may be due to their functional importance for the precision of ATP production usually. Length difference in mitochondrial control regions between DBHC and DBQY was also conversed at the population level, indicating that sequence size adjustments in control region may be associated with adaptation to divergent altitudes. Quantitatively real-time PCR analysis for 12 MPCGs showed that gene expression patterns had a significant change between the two populations, suggesting that expression levels of MPCGs could be modulated by divergent environmental pressures (e.g. oxygen content and ambient temperature). These results provided an important guide for further uncovering genetic mechanisms of ecological adaptation in insects.

The Complete Mitogenome of Pyrrhocoris tibialis (Hemiptera: Pyrrhocoridae) and Phylogenetic Implications.

We determined the complete mitogenome of Pyrrhocoris tibialis (Hemiptera: Heteroptera: Pyrrhocoridae) to better understand the diversity and phylogeny within Pentatomomorpha, which is the second largest infra-order of Heteroptera. Gene content, gene arrangement, nucleotide composition, codon usage, ribosomal RNA (rRNA) structures, and sequences of the mitochondrial transcription termination factor were well conserved in Pyrrhocoroidea. Different protein-coding genes have been subject to different evolutionary rates correlated with the G + C content. The size of control regions (CRs) was highly variable among mitogenomes of three sequenced Pyrrhocoroidea species, with the P. tibialis CR being the largest. All the transfer RNA genes found in Pyrrhocoroidea had the typical clover leaf secondary structure, except for trnS1 (AGN), which lacked the dihydrouridine arm and possessed an unusual anticodon stem (9 bp vs. the normal 5 bp). A total of three different phylogenetic relationships among the five super-families of Pentatomomorpha were obtained using three analytical methods (MrBayes and RAxML under site-homogeneous models and PhyloBayes under a site-heterogeneous CAT + GTR model) and two mitogenomic datasets (nucleotides and amino acids). The tree topology test using seven methods statistically supported a phylogeny of (Aradoidea + (Pentatomoidea + (Lygaeoidea + (Pyrrhocoroidea + Coreoidea)))) as the best topology, as recognized by both RAxML and MrBayes based on the two datasets.

A small set of differentially expressed genes was associated with two color morphs in natural populations of the pea aphid Acyrthosiphon pisum.

Color polymorphism is an ecologically important trait, which is related to local adaptation and ecological speciation. The pea aphid Acyrthosiphon pisum shows color polymorphism: the red and green color morphs where differences in ecological adaptation have been observed. Here, we measured genome-wide gene expression profiles of two color morphs in natural populations of A. pisum to explore the genetic basis of differentiated ecological adaptation. The results showed that only 32 genes were significantly differentially expressed between the two morphs, of which 18 had functional annotations. Among them, 13 genes were up-regulated [e.g. genes encoding protoheme IX farnesyltransferase (LOC100570971), carotene dehydrogenase (tor) and V-type proton ATPase subunit B (LOC100169462)] and 5 genes were down-regulated in the red morph (e.g. genes encoding transcription factors and heat shock proteins). To assess the functional importance of these differentially expressed genes (DEGs), we selected three highly expressed DEGs (LOC100169462, LOC100570971 and tor) with functional annotations and analyzed their expression levels in the red morph under three low temperatures (1 °C, 4 °C, and 8 °C) for 24 h. These three DEGs showed an interesting expression response to the cold acclimating conditions which resulted in an obvious phenotypic change of the red individuals to be greenish variants. This study suggests a link between gene expressions and body color polymorphisms in the pea aphid and provides important clues for further studying molecular mechanisms of ecological adaptation in aphids.

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.