The phylogeny of Galerucinae (Coleoptera: Chrysomelidae) and the performance of mitochondrial genomes in phylogenetic inference compared to nuclear rRNA genes

Comparative Mitogenomics and Phylogenetic Implications for Nine Species of the Subfamily Meconematinae (Orthoptera: Tettigoniidae)

Currently, the subfamily Meconematinae encompasses 1029 species, but whole-mitochondrial-genome assemblies have only been made available for 13. In this study, the whole mitochondrial genomes (mitogenomes) of nine additional species in the subfamily Meconematinae were sequenced. The size ranged from 15,627 bp to 17,461 bp, indicating double-stranded circular structures. The length of the control region was the main cause of the difference in mitochondrial genome length among the nine species. All the mitogenomes including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and a control region (CR). The majority strand encoded 23 genes, and the minority strand encoded 14 genes. A phylogenetic analysis reaffirmed the monophyletic status of each subfamily, but the monophysitism of Xizicus, Xiphidiopsis and Phlugiolopsis was not supported.

Mitogenome-based phylogeny of mosquitoes (Diptera: Culicidae)

Mosquitoes are of great medical significance as vectors of many deadly diseases. Mitogenomes have been widely used in phylogenetic studies, but mitogenome knowledge within the family Culicidae is limited, and Culicidae phylogeny is far from resolved. In this study, we surveyed the mitogenomes of 149 Culicidae species, including 7 newly sequenced species. Comparative analysis of 149 mosquito mitogenomes shows gene composition and order to be identical to that of an ancestral insect, and the AT bias, length variation, and codon usage are all consistent with that of other reported Dipteran mitogenomes. Phylogenetic analyses based on the DNA sequences of the 13 protein-coding genes from the 149 species robustly support the monophyly of the subfamily Anophelinae and the tribes Aedini, Culicini, Mansoniini, Sabethini, and Toxorhynchitini. To resolve ambiguous relationships between clades within the subfamily Culicinae, we performed topological tests and show that Aedini is a sister to Culicini and that Uranotaeniini is a sister to (Mansoniini + (Toxorhynchitini + Sabethini)). In addition, we estimated divergence times using a Bayesian relaxation clock based on the sequence data and 3 fossil calibration points. The results show mosquitoes diverged during the Early Jurassic with massive Culicinae radiations during the Cretaceous, coincident with the emergence of angiosperms and the burst of mammals and birds. Overall, this study, which uses the largest number of Culicidae mitogenomes sequenced to date, comprehensively reveals the mitogenome characteristics and mitogenome-based phylogeny and divergence times of Culicidae, providing information for further studies on the mitogenome, phylogeny, evolution, and taxonomic revision of Culicidae.

Characterization of Four Complete Mitogenomes of Monolepta Species and Their Related Phylogenetic Implications

Monolepta is one of the diverse genera in the subfamily Galerucinae, including 708 species and 6 sub-species worldwide. To explore the information on the mitogenome characteristics and phylogeny of the section "Monoleptites", especially the genus Monolepta, we obtained the newly completed mitochondrial genomes (mitogenomes) of four Monolepta species using high-throughput sequencing technology. The lengths of these four new mitochondrial genomes are 16,672 bp, 16,965 bp, 16,012 bp, and 15,866 bp in size, respectively. All four mitochondrial genomes include 22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), and one control region, which is consistent with other Coleoptera. The results of the nonsynonymous with synonymous substitution rates showed that ND6 had the highest evolution rate, while COI displayed the lowest evolution rate. The substitution saturation of three datasets (13 PCGs_codon1, 13 PCGs_codon2, 13 PCGs_codon3) showed that there was no saturation across all datasets. Phylogenetic analyses based on three datasets (ND1, 15 genes of mitogenomes, and 13 PCGs_AA) were carried out using maximum likelihood (ML) and Bayesian inference (BI) methods. The results showed that mitogenomes had a greater capacity to resolve the main clades than the ND1 gene at the suprageneric and species levels. The section "Monoleptites" was proven to be a monophyletic group, while Monolepta was a non-monophyletic group. Based on ND1 data, the newly sequenced species whose antennal segment 2 was shorter than 3 were split into several clades, while, based on the mitogenomic dataset, the four newly sequenced species had close relationships with Paleosepharia. The species whose antennal segment 2 was as long as 3 were split into two clades, which indicated that the characteristic of "antennal segment 2 as long as 3" of the true "Monolepta" evolved multiple times in several subgroups. Therefore, to explore the relationships among the true Monolepta, the most important thing is to perform a thorough revision of Monolepta and related genera in the future.

Can Pheromones Contribute to Phylogenetic Hypotheses? A Case Study of Chrysomelidae

Pheromones mediate species-level communication in the search for mates, nesting, and feeding sites. Although the role of pheromones has long been discussed by various authors, their existence was not proven until the mid-twentieth century when the first sex pheromone was identified. From this finding, much has been speculated about whether this communication mechanism has acted as a regulatory agent in the process of speciation, competition, and sexual selection since it acts as an intraspecific barrier. Chrysomelidae is one of the major Phytophaga lineages, with approximately 40,000 species. Due to this immense diversity the internal relationships remain unstable when analyzed only with morphological data, consequently recent efforts have been directed to molecular analyses to establish clarity for the relationships and found their respective monophyly. Therefore, our goals are twofold 1) to synthesize the current literature on Chrysomelidae sex pheromones and 2) to test whether Chrysomelidae sex pheromones and their chemical structures could be used in phylogenetic analysis for the group. The results show that, although this is the first analysis in Chrysomelidae to use pheromones as a phylogenetic character, much can be observed in agreement with previous analyses, thus confirming that pheromones, when known in their entirety within lineages, can be used as characters in phylogenetic analyses, bringing elucidation to the relationships and evolution of organisms.

Characterization of the complete mitochondrial genome of the longhorn beetle, Batocerahorsfieldi (Coleoptera, Cerambycidae) and its phylogenetic analysis with suitable longhorn beetles

Mitochondrial genome analysis is an important tool for studying insect phylogenetics. The longhorn beetle, Batocerahorsfieldi, is a significant pest in timber, economic and protection forests. This study determined the mitochondrial genome of B.horsfieldi and compared it with the mitochondrial genomes of other Cerambycidae with the aim of exploring the phylogenetic status of the pest and the evolutionary relationships among some Cerambycidae subgroups. The complete mitochondrial genome of B.horsfieldi was sequenced by the Illumina HiSeq platform. The mitochondrial genome was aligned and compared with the existing mitochondrial genomes of Batoceralineolata and B.rubus in GenBank (MF521888, MW629558, OM161963, respectively). The secondary structure of transfer RNA (tRNA) was predicted using tRNAScan-SE server v.1.21 and MITOS WebSever. Thirteen protein-coding genes (PCGs) and two ribosomal RNA gene sequences of 21 longhorn beetles, including B.horsfieldi, plus two outgroups, Dryopsernesti (Dryopidae) and Heterocerusparallelus (Heteroceridae), were analyzed. The phylogenetic tree was constructed using maximum likelihood and Bayesian inference methods. In this study, we successfully obtained the complete mitochondrial genome of B.horsfieldi for the first time, which is 15 425 bp in length. It contains 37 genes and an A + T-rich region, arranged in the same order as the recognized ancestor of longhorn beetles. The genome of B.horsfieldi is composed of 33.12% A bases, 41.64% T bases, 12.08% C bases, and 13.16% G bases. The structure, nucleotide composition, and codon usage of the new mitochondrial genome are not significantly different from other longhorn mitochondrial genomes. Phylogenetic analyses revealed that Cerambycidae formed a highly supported single clade, and Vesperidae was either clustered with Cerambycidae or formed a separate clade. Interestingly, B.horsfieldi, B.rubus and B.lineolata were clustered with Monochamus and Anoplophora species in both analyses, with high node support. Additionally, the VesperidaeSpiniphilusspinicornis and Vesperussanzi and the 19 Cerambycidae species formed a sister clade in the Bayesian analysis. Our results have produced new complete mitogenomic data, which will provide information for future phylogenetic and taxonomic research, and provide a foundation for future relevant research.

Giant Sex Chromosomes in Omophoita Species (Coleoptera, Chrysomelidae): Structural and Evolutionary Relationships Revealed by Zoo-FISH and Comparative Genomic Hybridization (CGH)

The beetles of the subtribe Oedionychina (Chrysomelidae, Alticinae) are the only ones that have the atypical giant and achiasmatic sex chromosomes, which are substantially larger than the autosomes. Previous cytogenetic analyses suggest a large accumulation of repetitive DNA in the sex chromosomes. In this study, we examined the similarity of X and Y chromosomes in four Omophoita species and compared genomic differentiation to better understand the evolutionary process and the giant sex chromosomes origin. Intraspecific genomic comparation using male and female genomes of O. octoguttata and interespecific analyses using genomic DNA of O. octoguttata, O. sexnotata, O. magniguttis, and O. personata were performed. In addition, whole chromosome painting (WCP) experiments were performed with X and Y chromosome probes of O. octogutatta. CGH analysis revealed great genomic similarity between the sexes and a sex-specific region on the Y chromosome, and interspecific analysis revealed a genomic divergence between species. In contrast, WCP results revealed that the sex chromosomes of O. octoguttata have high intra- and interspecific similarity with the studied species. Our data support a common origin under the canonical evolution of the sex chromosomes in this group, as they have high genomic similarity between them.

CO1 barcodes resolve an asymmetric biphyletic clade for Diabrotica undecimpunctata subspecies and provide nucleotide variants for differentiation from related lineages using real-time PCR

Diabrotica undecimpunctata is a multivoltine polyphagous beetle species that has long been documented as a significant agricultural pest throughout its native range in North America. This beetle can vector bacterial and viral plant pathogens that result in major losses to crops such as cucumber and soybean. Many countries outside the Americas treat D. undecimpunctata as a species of quarantine importance, while in the USA only the subspecies D. u. duodecimnotata is subject to quarantine, to prevent introduction from Mexico. Identification of D. undecimpunctata on the basis of morphology alone can be complicated given the use of conflicting characters in the description of some subspecific taxa. To better understand relationships among D. undecimpunctata subspecies and other related species, we sequenced mitochondrial cytochrome oxidase 1 (CO1) and nuclear internal transcribed spacer 2 (ITS2) DNA from individuals in different subspecific taxa and across different parts of the species range using museum samples and interceptions. When our data were combined with publicly available Diabrotica data, no pattern of divergence consistent with the currently recognized subspecific designations was found. In addition, we compared phylogenetic patterns in CO1 data from the congener D. virgifera to demonstrate the utility of mitochondrial data in resolving subspecies. From the CO1 data, a diagnostic real-time PCR assay was developed that could successfully identify all haplotypes within the large D. undecimpunctata clade for use in surveys and identification at ports of entry. These findings underscore the need to resolve molecular and morphological datasets into cogent, lineage-based groupings. Such efforts will provide an evolutionary context for the study of agriculturally important attributes of Diabrotica such as host preferences, xenobiotic metabolism, and natural and anthropogenic patterns of dispersal.

Up and down from North to South: Latitudinal Distribution of Flea Beetle Genera in Continental Africa (Coleoptera, Chrysomelidae, Galerucinae, Alticini)

The distribution of global biodiversity can be investigated based on comprehensive datasets and many methods to process them. The taxonomic diversity of phytophagous insects is typically linked to plant diversity, which increases from temperate to tropical latitudes. In this paper, we explored the latitudinal distribution of the flea beetle genera (Coleoptera, Chrysomelidae, Galerucinae, Alticini) on the African continent. We divided the area into latitudinal belts and looked for possible correlations with the number and types of vegetational divisions, the area of each belt, and the bioclimatic variables. The number of flea beetle genera is related to the number and types of vegetation divisions rather than the area of each belt. Some bioclimatic variables are highly related to the number of genera, which is higher within those belts where climate factors limit the oscillation of temperature over the year and favor high precipitations, especially in the warmest months. These biotic and abiotic factors lead to a two-peak trend in the taxonomic richness of flea beetle genera from north to south. Genera endemic to restricted areas are linked to the presence of high mountain systems and increase the taxonomic richness of the belt they belong to.

Comparative Mitogenome Analyses Uncover Mitogenome Features and Phylogenetic Implications of the Parrotfishes (Perciformes: Scaridae)

In order to investigate the molecular evolution of mitogenomes among the family Scaridae, the complete mitogenome sequences of twelve parrotfish species were determined and compared with those of seven other parrotfish species. The comparative analysis revealed that the general features and organization of the mitogenome were similar among the 19 parrotfish species. The base composition was similar among the parrotfishes, with the exception of the genus Calotomus, which exhibited an unusual negative AT skew in the whole mitogenome. The PCGs showed similar codon usage, and all of them underwent a strong purifying selection. The gene rearrangement typical of the parrotfishes was detected, with the tRNAMet inserted between the tRNAIle and tRNAGln, and the tRNAGln was followed by a putative tRNAMet pseudogene. The parrotfish mitogenomes displayed conserved gene overlaps and secondary structure in most tRNA genes, while the non-coding intergenic spacers varied among species. Phylogenetic analysis based on the thirteen PCGs and two rRNAs strongly supported the hypothesis that the parrotfishes could be subdivided into two clades with distinct ecological adaptations. The early divergence of the sea grass and coral reef clades occurred in the late Oligocene, probably related to the expansion of sea grass habitat. Later diversification within the coral reef clade could be dated back to the Miocene, likely associated with the geomorphology alternation since the closing of the Tethys Ocean. This work provided fundamental molecular data that will be useful for species identification, conservation, and further studies on the evolution of parrotfishes.

Characterization of the Complete Mitochondrial Genome of a Flea Beetle Luperomorpha xanthodera (Coleoptera: Chrysomelidae: Galerucinae) and Phylogenetic Analysis

In this study, the mitochondrial genome of Luperomorpha xanthodera was assembled and annotated, which is a circular DNA molecule including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA genes (12S rRNA and 16S rRNA), and 1388 bp non-coding regions (A + T rich region), measuring 16,021 bp in length. The nucleotide composition of the mitochondrial genome is 41.3% adenine (A), 38.7% thymine (T), 8.4% guanine (G), and 11.6% cytosine (C). Most of the protein-coding genes presented a typical ATN start codon (ATA, ATT, ATC, ATG), except for ND1, which showed the start codon TTG. Three-quarters of the protein-coding genes showed the complete stop codon TAR (TAA, TAG), except the genes COI, COII, ND4, and ND5, which showed incomplete stop codons (T- or TA-). All the tRNA genes have the typical clover-leaf structure, except tRNASer1 (AGN), which has a missing dihydrouridine arm (DHU). The phylogenetic results determined by both maximum likelihood and Bayesian inference methods consistently supported the monophyly of the subfamily Galerucinae and revealed that the subtribe Luperina and genus Monolepta are polyphyletic groups. Meanwhile, the classification status of the genus Luperomorpha is controversial.

Comparative mitochondrial genomic analysis provides new insights into the evolution of the subfamily Lamiinae (Coleoptera: Cerambycidae)

The genus Monochamus within the subfamily Lamiinae is the main vector of Bursaphelenchus xylophilus, which causes pine wilt disease and induces substantial economic and ecological losses. Only three complete mitochondrial genomes of the genus Monochamus have been sequenced to date, and no comparative mitochondrial genomic studies of Lamiinae have been conducted. Here, the mitochondrial genomes of two Monochamus species, M. saltuarius and M. urussovi, were newly sequenced and annotated. The composition and order of genes in the mitochondrial genomes of Monochamus species are conserved. All transfer RNAs exhibit the typical clover-leaf secondary structure, with the exception of trnS1. Similar to other longhorn beetles, Lamiinae mitochondrial genomes have an A + T bias. All 13 protein-coding genes have experienced purifying selection, and tandem repeat sequences are abundant in the A + T-rich region. Phylogenetic analyses revealed congruent topologies among trees inferred from the five datasets, with the monophyly of Acanthocinini, Agapanthiini, Batocerini, Dorcaschematini, Pteropliini, and Saperdini receiving high support. The findings of this study enhance our understanding of mitochondrial genome evolution and will provide a basis for future studies of population genetics and phylogenetic investigations in this group.

Analysis of the Mitogenomes of Two Helotid Species Provides New Insights into the Phylogenetic Relationship of the Basal Cucujoidea (Insecta: Coleoptera)

Helotid beetles are commonly found in places where sap flows from tree trunks and in crevices in bark. The Helotidae family is a rare and primitive group of Cucujoidea. To date, no complete mitochondrial (mt) genome has been sequenced for this family. To better understand the characteristics of the mt genome and the evolution of Cucujoidea, we sequenced and annotated the complete mt genomes of Helota thoracica (Ritsema, 1895) and Helota yehi Lee, 2017 using next-generation sequencing. These are the first record of Helotidae mt genomes. The RNA secondary structures of both species were also predicted in this study. The mt genomes of H. thoracica and H. yehi are circular, with total lengths of 16,112 bp and 16,401 bp, respectively. After comparing the mt genomes of H. thoracica and H. yehi, we observed the gene arrangement, codon usage patterns, base content, and RNA secondary structures of both species to be similar, which has also been noted in other Coleoptera insects. The nucleotide sequence of the coding regions and the control region has small differences. The phylogenetic analysis indicated that Helotidae and Protocucujidae are sister groups and revealed the relationship between seven families; however, the validity of the two series (Erotylid series and Nitidulid series) as larger groups in the superfamily was not supported. The mt phylogenomic relationships have strong statistical support. Therefore, the division of Cucujoidea into series should be re-examined. Our results will provide a better understanding of the mt genome and phylogeny of Helotidae and Cucujoidea and will provide valuable molecular markers for further genetic studies.

Complete mitochondrial genome of Distenia punctulatoides (Coleoptera: Chrysomeloidea: Disteniinae) and its phylogenetic implications

The family Disteniidae is a moderately large and widely distributed lineage. Distenia punctulatoides belongs to the family Disteniidae from the cerambycoid assemblage. Here, we report the complete mitogenome of D. punctulatoides, which is 15,675 bp in length. It contains 37 genes and a noncoding control region, which are arranged in the same order as that of the putative ancestor of beetles. The total base composition of the new mitogenome is 40.2% for A, 17.1% for C, 10.0% for G, and 32.7% for T. The new mitogenomic organization, nucleotide composition, and codon usage do not differ significantly from other beetles. Using available complete mitogenomes, the high-level phylogeny of the family Disteniidae was explored. The phylogenetic analyses showed that Disteniidae were monophyletic, and the genus Distenia grouped with the genus Clytomelegena as sister groups. Combining the morphological and molecular data, Typodryas Thomson, 1864 is suggested to be a junior synonym of Distenia Lepeletier and Audinet-Serville, 1828.

Argopistes sexvittatus and Argopistes capensis (Chrysomelidae: Alticini): Mitogenomics and Phylogeny of Two Flea Beetles Affecting Olive Trees

The genus Argopistes (Chrysomelidae: Alticini) is the only group of flea beetles specialized in plant hosts in the family Oleaceae. In southern Africa, Argopistes are often found feeding on African Wild Olive (Olea europaea subsp. cuspidata) and European cultivated olive (O. e. subsp. europaea), and heavy infestations can be devastating to mature trees and compromise the development of young trees. Despite their negative agricultural impact, African Argopistes are an understudied group for which no genetic data were available. We assessed the species diversity of olive flea beetles in the Western Cape province of South Africa, the largest olive-producing region in sub-Saharan Africa, by collecting adult specimens on wild and cultivated olive trees between 2015 and 2017. Argopistes sexvittatus Bryant, 1922 (n = 289) dominated at all sampling sites, and Argopistes capensis Bryant, 1944 (n = 2) was found only once. Argopistes oleae Bryant, 1922, a third species previously reported in the region, was not found. The complete mitogenomes of one A. capensis and two A. sexvittatus (striped and black morphotypes) individuals were sequenced for phylogenetic reconstruction in the context of other 64 species. The two olive flea beetle species form a monophyletic clade with other Argopistes, supporting the hypothesis that the exclusive feeding habit on Oleaceae is an evolutionary adaptation in this genus.

First report of complete mitochondrial genome in the subfamily Alleculinae and mitochondrial genome-based phylogenetics in Tenebrionidae (Coleoptera: Tenebrionoidea)

Despite worldwide distribution and rich diversity, the knowledge of mitochondrial genome (mtgenome) characteristics within the family Tenebrionidae (Coleoptera) is still very limited, and phylogenetics remain unresolved for the family. In this study, the mtgenomes for 12 species are newly sequenced and annotated. Three of the species represent the first report of complete mtgenomes from the subfamily Alleculinae. Comparative analysis of 36 tenebrionid mtgenomes finds that gene composition and order are the same as a theoretical ancestral insect mtgenome, and AT bias, length variation, and codon usage are consistent with other reported beetle mtgenomes. Some intergenic overlap and gap sequences may contain phylogenetically informative information, whereas neither the conserved nor repeat sequences in the control region do. The subfamilies Lagriinae, Pimeliinae, Stenochiinae, and Alleculinae are found to be monophyletic, and the subfamilies Diaperinae and Tenebrioninae polyphyletic in our study. Furthermore, Lagriinae is sister to the rest of the subfamilies. At the tribal level, the tribes Lagriini in the subfamily Lagrrinae; Strongyliini in Stenochiinae; Cteniopodini in Alleculinae; and Triboliini, Opatrini, and Amarygmini in Tenebrioninae are monophyletic; while Diaperini in Diaperinae, and Tenebrionini in Tenebrioninae are polyphyletic.

Phylogenetic Implication of Large Intergenic Spacers: Insights from a Mitogenomic Comparison of Prosopocoilus Stag Beetles (Coleoptera: Lucanidae)

Simple Summary

Insect mitochondrial genomes (mitogenomes) show high diversity in some lineages. In the mitogenome of some Coleoptera species, a large intergenic spacer (IGS) has been identified. However, very little is known about mitogenomes of lucanid beetles. In this work, to provide further insight into the phylogenic relationships among species in lucanid beetles (genus Prosopocoilus), two Prosopocoilus species (Prosopocoilus castaneus and Prosopocoilus laterotarsus) were newly sequenced and comparatively analyzed. Significantly, the two newly sequenced Prosopocoilus species contained a large IGS located between trnI and trnQ. Our phylogenomic analyses showed that P. castaneus and P. laterotarsus were clustered in a clade with typical Prosopocoilus species (Prosopocoilus confucius, Prosopocoilus blanchardi, and Prosopocoilusastacoides). These results provide valuable data for the future study of the phylogenetic relationships in this genus.

Abstract

To explore the characteristics of mitogenomes and discuss the phylogenetic relationships within the genus Prosopocoilus, the mitogenomes of two species (P. castaneus and P. laterotarsus) were newly sequenced and comparatively analyzed. The arrangement of the mitogenome in these two lucanid beetles was the same as that in the inferred ancestral insect, and the nucleotide composition was highly biased towards A + T as in other lucanids. The evolutionary rates of 13 protein-coding genes (PCGs) suggested that their evolution was based on purifying selection. Notably, we found evidence of the presence of a large IGS between trnI and trnQ genes, whose length varied from 375 bp (in P. castaneus) to 158 bp (in P. laterotarsus). Within the large IGS region, a short sequence (TAAAA) was found to be unique among these two species, providing insights into phylogenomic reconstruction. Phylogenetic analyses were performed using the maximum likelihood (IQ-TREE) and Bayesian (PhyloBayes) methods based on 13 protein-coding genes (PCGs) in nucleotides and amino acids (AA) from published mitogenomes (n = 29). The genus Prosopocoilus was found to constitute a distinct clade with high nodal support. Overall, our findings suggested that analysis of the characteristics of the large IGS (presence or absence, size, and location) in mitogenomes of the genus Prosopocoilus may be informative for the phylogenetic and taxonomic analyses and for evaluation of the genus Prosopocoilus, despite the dense sampling materials needed.

Higher-level phylogeny of Chrysomelidae based on expanded sampling of mitogenomes

Chrysomelidae is one of the most diverse lineages of beetles. The classification and phylogeny of Chrysomelidae have been contentious. In this study, we obtained 16 new mitogenome sequences by using next-generation sequencing. Combined with the published mitogenomes, we inferred the phylogenetic relationships of Chrysomelidae. Different data recoding strategies and substitution models were applied to phylogenetic reconstruction. In the Maximum likelihood analyses under the homogeneous model, Dayhoff recoding allowed for the improved phylogenetic resolution due to the decreased level of heterogeneous sequence divergence. Bayesian inference under the heterogeneous model yielded generally well resolved subfamily relationships. The present mitogenome data strongly supported Chrysomelidae as a monophyletic group. Consistent with previous work, we found three major subfamily clades within Chrysomelidae. However, the pattern of the “sagrine” clade plus the “eumolpine” clade being sister to the “chrysomeline” clade contrasted with the prior study. The placement of the genus Syneta with regards to these three clades was ambiguous. Relationships recovered suggested several major chrysomelid lineages, including: (1) Bruchinae in the “sagrine” clade; (2) Donaciinae + Criocerinae; (3) Spilopyrinae + (Cassidinae + (Eumolpinae + (Lamprosomatinae + Cryptocephalinae))); (4) Chrysomelinae + (Alticinae + Galerucinae). Results also suggested the placement of Timarcha outside the major Chrysomelinae.

Multiple origins of moss-inhabiting flea beetles (Coleoptera: Chrysomelidae): molecular phylogeny, overview of genera and a new genus from Africa

Flea beetles (Alticini) are a diverse lineage of leaf beetles (Coleoptera: Chrysomelidae) with the majority of genera living on vascular plants. However, several genera are associated with mosses or leaf litter, as an alternative strategy, often accompanied by changes in morphology and biology. Here, we present the first comprehensive DNA-based phylogenetic analysis of Alticini including 13 moss-inhabiting genera. We show that moss-association is a common phenomenon in Alticini; based on our ancestral-state mapping, we reveal at least nine independent origins of moss-inhabiting, and four independent origins of leaf-litter-inhabiting life-strategy. Our findings illustrate that both cryptic strategies induce parallel evolution of similar ecomorphological adaptations, such as flight loss and a convex body shape. An overview of known moss-inhabiting flea beetle genera is also provided. We describe Benedictoides munclingeri gen. et sp. nov. from the cloud forests of Mount Cameroon, the first flea beetle genus clearly assigned to the moss-inhabiting strategy from tropical Africa. Additionally, our phylogeny provides new insights into the relationships among various flea beetle genera. We estimate six main clades grouping the majority of studied genera, and propose several new generic groups. Along with our phylogenetic results, we discuss moss-inhabiting flea beetle evolution and perspectives of their systematics.

Lack of Genetic Structure Among Populations of Striped Flea Beetle Phyllotreta striolata (Coleoptera: Chrysomelidae) Across Southern China

The striped flea beetle (SFB) Phyllotreta striolata (Fabricius) (Coleoptera: Chrysomelidae) is a major pest of cruciferous vegetables in southern China. The population diversity and genetic structure of SFB are unknown. Here, we assembled a draft genome for the SFB and characterized the distribution of microsatellites. Then, we developed 12 novel microsatellite markers across the genome. We used a segment of the cox1 gene and newly developed microsatellite markers to genotype the genetic diversity of SFB across southern China. There were 44 mitochondrial haplotypes in the SFB populations, with haplotype 2 as the most widespread. The population genetic differentiation was very low, indicated by FST-values (<0.05 except for Guangxi population with other populations based on cox1), high gene flow (4.10 and 44.88 of cox1 and microsatellite, respectively) and Principal Coordinate Analysis across all populations. Mantel test showed genetic distance in SFB was significantly associated with geographic distance based on microsatellites (R2 = 0.2373, P = 0.014) while result based on cox1 (R2 = 0.0365, P = 0.155) showed no significant difference. The phylogenetic analysis did not find any geographically related clades among all haplotypes. Analyses based on microsatellites showed a lack of population genetic structure among all populations. Our study provides a foundation for the future understanding of the ecology and evolution of SFB and its management.

Phylogenetic signal of sub-arctic beetle communities

Postglacial dispersal and colonization processes have shaped community patterns in sub-Arctic regions such as Churchill, Manitoba, and Canada. This study investigates evolutionary community structure within the beetle (Coleoptera) families of Churchill and tests whether biological traits have played a role in governing colonization patterns from refugial and southerly geographic regions. This study quantifies sub-Arctic beetle phylogenetic community structure for each family using the net relatedness index (NRI) and nearest taxon index (NTI), calculated using publicly available data from the Barcode of Life Data Systems (BOLD); compares patterns across families with different traits (habitat, diet) using standard statistical analysis (ANOVA) as well as phylogenetic generalized least squares (PGLS) using a family-level beetle phylogeny obtained from the literature; and compares community structure in Churchill with a region in southern Canada (Guelph, Ontario). These analyses were also repeated at a genus level. The dominant pattern detected in our study was that aquatic families were much better represented in Churchill compared to terrestrial families, when compared against richness sampled from across Canada and Alaska. Individually, most families showed significant phylogenetic clustering in Churchill, likely due to the strong environmental filtering present in Arctic environments. There was no significant difference in phylogenetic structure between Churchill and Guelph but with a trend toward stronger clustering in the North. Fungivores were significantly more overdispersed than other feeding modes, predators were significantly more clustered, and aquatic families showed significantly stronger clustering compared to terrestrial. This study contributes to our understanding of the traits and processes structuring insect biodiversity and macroecological trends in the sub-Arctic.

New data on Henneguya postexilis Minchew, 1977, a parasite of channel catfish Ictalurus punctatus, with notes on resolution of molecular markers for myxozoan phylogeny

Previous morphological and histological data are supplemented with molecular and ultrastructural data for a Henneguya sp. isolated from farm-raised channel catfish Ictalurus punctatus in Mississippi, USA. Myxospores were cryptic, encapsulated within a thin layer of epithelium in the gill lamellae with spore measurements consistent with the original description of Henneguya postexilis Minchew, 1977. Myxospores were 42.7-49.1 µm in total length with spore bodies 12.1-17.2 × 3.6-4.8 × 2.9-3 µm. Polar capsules were of unequal length, with the longer capsule being 4.4-6.7 × 1.1-1.6 µm and the shorter capsule being 4.4-6.4 × 1.1-1.6 µm. Polar tubules had 6-8 turns. Caudal processes were 25.7-38.1 µm in length. Spores were encapsulated in a thin layer of epithelium in the gill lamellae. Molecular data from the most commonly used markers for myxozoan identification and phylogeny, partial 18S small subunit ribosomal gene (SSU), partial 28S large subunit ribosomal gene (LSU), and elongation factor 2 (EF2) were generated for H. postexilis. Additionally, novel data for LSU and EF2 were generated for archived myxozoan specimens from farm-raised catfish (H. mississippiensis, H. ictaluri, H. exilis, H. adiposa, H. sutherlandi, H. bulbosus, Unicauda fimbrethilae), as well as archived specimens from wild fish (H. laseeae [from Pylodictis olivaris], Hennegoides flockae [from Aphredoderus sayanus], Myxobolus cloutmani [from Cycleptus elongatus]. These include the first EF2 sequence data for the genera Hennegoides and Unicauda. Phylogenetic analyses using these data placed H. postexilis in well supported clades with other ictalurid-infecting Henneguya species. Phylogenetic signal assessments on these analyses suggest that while SSU provided the greatest phylogenetic signal, LSU yielded comparable signal, supporting previous work implying this region may be underutilised in reconstructing myxobolid phylogenies.

Comparative and phylogenomic analyses of mitochondrial genomes in Coccinellidae (Coleoptera: Coccinelloidea)

The Coccinellidae are one of the most familiar beetle families, the ladybirds. Despite the great ecological and economic significance, the phylogenetic relationships of Coccinellidae remain poorly understood. One of the reasons is that the sequenced mitogenomes available for this family are very limited. We sequenced complete or nearly complete mitogenomes from seven species of the tribe Coccinellini with next-generation sequencing. All species have the same gene content and gene order as the putatively ancestral insect mitogenome. A large intergenic spacer region (> 890 bp) was found located between trnI and trnQ. The potential for using secondary structures of the large and small ribosomal subunits for phylogenetic reconstruction was predicted. The phylogenetic relationships were explored through comparative analyses across more than 30 coccinellid species. We performed phylogenetic analyses with both concatenation methods (Maximum Likelihood and Bayesian Inference) and multispecies coalescent method (ASTRAL). Phylogenetic results strongly supported the monophyly of Coccinellidae. Within Coccinellidae, the Epilachnini and the Coccinellini including Halyziini were monophyletic, while the Scymnini and Coccidulini were non-monophyletic.

Nine Mitochondrial Genomes of the Pyraloidea and Their Phylogenetic Implications (Lepidoptera)

Simple Summary

The Pyraloidea is a large superfamily of Lepidoptera in species composition. To date, the higher-level phylogenetic relationships in this group remain unresolved, and many taxa, with taxonomic positions historically established by morphological characters, need to be confirmed through sequencing of DNA, including mitochondrial genome sequences (mitogenomes). Here, we newly generated nine complete mitogenomes for Pyraloidea that shared identical gene content, and arrangements that are typical of Lepidoptera. The current phylogenetic results confirmed previous multilocus studies, indicating the effectiveness of mitogenomes for inference of Pyraloidea higher-level relationships. Unexpectedly, Orybina Snellen was robustly placed as basal to the remaining Pyralidae taxa, rather than nested in the Pyralinae of Pyralidae as morphologically defined and placed. Our results bring a greater understanding to Pyraloidea phylogeny, and highlight the necessity of sequencing more pyraloid taxa to reevaluate their phylogenetic positions.

Abstract

The Pyraloidea is one of the species-rich superfamilies of Lepidoptera and contains numerous economically important pest species that cause great loss in crop production. Here, we sequenced and annotated nine complete mitogenomes for Pyraloidea, and further performed various phylogenetic analyses, to improve our understanding of mitogenomic evolution and phylogeny of this superfamily. The nine mitogenomes were circular, double-stranded molecules, with the lengths ranging from 15,214 bp to 15,422 bp, which are comparable to other reported pyraloid mitogenomes in size. Gene content and arrangement were highly conserved and are typical of Lepidoptera. Based on the hitherto most extensive mitogenomic sampling, our various resulting trees showed generally congruent topologies among pyraloid subfamilies, which are almost in accordance with previous multilocus studies, indicating the suitability of mitogenomes in inferring high-level relationships of Pyraloidea. However, nodes linking subfamilies in the “non-PS clade” were not completely resolved in terms of unstable topologies or low supports, and future investigations are needed with increased taxon sampling and molecular data. Unexpectedly, Orybina Snellen, represented in a molecular phylogenetic investigation for the first time, was robustly placed as basal to the remaining Pyralidae taxa across our analyses, rather than nested in Pyralinae of Pyralidae as morphologically defined. This novel finding highlights the need to reevaluate Orybina monophyly and its phylogenetic position by incorporating additional molecular and morphological evidence.

The mitochondrial genomes of Tortricidae: nucleotide composition, gene variation and phylogenetic performance

BACKGROUND

Mitochondrial genomes (mitogenomes) have greatly improved our understanding of the backbone phylogeny of Lepidoptera, but few studies on comparative mitogenomics below the family level have been conducted. Here, we generated 13 mitogenomes of eight tortricid species, reannotated 27 previously reported mitogenomes, and systematically performed a comparative analysis of nucleotide composition, gene variation and phylogenetic performance.

RESULTS

The lengths of completely sequenced mitogenomes ranged from 15,440 bp to 15,778 bp, and the gene content and organization were conserved in Tortricidae and typical for Lepidoptera. Analyses of AT-skew and GC-skew, the effective number of codons and the codon bias index all show a base bias in Tortricidae, with little heterogeneity among the major tortricid groups. Variations in the divergence rates among 13 protein-coding genes of the same tortricid subgroup and of the same PCG among tortricid subgroups were detected. The secondary structures of 22 transfer RNA genes and two ribosomal RNA genes were predicted and comparatively illustrated, showing evolutionary heterogeneity among different RNAs or different regions of the same RNA. The phylogenetic uncertainty of Enarmoniini in Tortricidae was confirmed. The synonymy of Bactrini and Olethreutini was confirmed for the first time, with the representative Bactrini consistently nesting in the Olethreutini clade. Nad6 exhibits the highest phylogenetic informativeness from the root to the tip of the resulting tree, and the combination of the third coding positions of 13 protein-coding genes shows extremely high phylogenetic informativeness.

CONCLUSIONS

This study presents 13 mitogenomes of eight tortricid species and represents the first detailed comparative mitogenomics study of Tortricidae. The results further our understanding of the evolutionary architectures of tortricid mitogenomes and provide a basis for future studies of population genetics and phylogenetic investigations in this group.

Comparative Mitogenomics and Phylogenetic Analyses of Pentatomoidea (Hemiptera: Heteroptera)

Pentatomoidea is the largest superfamily of Pentatomomorpha; however, the phylogenetic relationships among pentatomoid families have been debated for a long time. In the present study, we gathered the mitogenomes of 55 species from eight common families (Acanthosomatidae, Cydnidae, Dinidoridae, Scutelleridae, Tessaratomidae, Plataspidae, Urostylididae and Pentatomidae), including 20 newly sequenced mitogenomes, and conducted comparative mitogenomic studies with an emphasis on the structures of non-coding regions. Heterogeneity in the base composition, and contrasting evolutionary rates were encountered among the mitogenomes in Pentatomoidea, especially in Urostylididae, which may lead to unstable phylogenetic topologies. When the family Urostylididae is excluded in taxa sampling or the third codon positions of protein coding genes are removed, phylogenetic analyses under site-homogenous models could provide more stable tree topologies. However, the relationships between families remained the same in all PhyloBayes analyses under the site-heterogeneous mixture model CAT + GTR with different datasets and were recovered as (Cydnidae + (((Tessaratomidae + Dinidoridae) + (Plataspidae + Scutelleridae)) + ((Acanthosomatidae + Urostylididae) + Pentatomidae)))). Our study showed that data optimizing strategies after heterogeneity assessments based on denser sampling and the use of site-heterogeneous mixture models are essential for further analysis of the phylogenetic relationships of Pentatomoidea.

Mitogenome of Coprophanaeus ensifer and phylogenetic analysis of the Scarabaeidae family (Coleoptera)

Several studies about the phylogenetic relationships of the Scarabaeinae subfamily (Coleoptera: Scarabaeidae) have been performed, but some phylogenetic uncertainties persist including the relationship and monophyly of different tribes and some genera. The aim of this study was to characterize the mitogenome of Coprophanaeus ensifer in order to establish its position within the Scarabaeidae family and to contribute to the resolution of some phylogenetic uncertainties. The mitogenome was sequenced on the Illumina HiSeq 4000, assembled using the Mitobim software and annotated in MITOS WebServer. The phylogenetic trees were reconstructed by Bayesian inference. The C. ensifer mitogenome is a molecule of 14,964 bp that contains the number and organization of the genes similar to those of most Coleoptera species. Phylogenetic reconstruction suggests monophyly of the tribe Phanaeini and supports the hypothesis that Coprini is a sister group of Phanaeini. The results also revealed the position of the tribe Oniticellini which is grouped with Onthophagini and Onitini. The geographic distribution of these species that form the most ancestral clade suggests with Scarabaeinae originated in Africa.

Mitochondrial genome of Monolepta hieroglyphica (Coleoptera: Chrysomeloidea: Chrysomelidae) and phylogenetic analysis

Monolepta hieroglyphica (Motschulsky, 1858) is a major pest of potato, maize, cotton and sorghum in China. In this study, we sequenced and analyzed the complete mitochondrial genome (mitogenome) of M. hieroglyphica. This mitogenome was 15,761 bp long and encoded 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs) and two ribosomal RNA unit genes (rRNAs). Gene order was conserved and identical to most other previously sequenced Galerucinae Most PCGs of M. hieroglyphica have the conventional start codons ATN (six ATT, five ATG and one ATC), with the exception of nad1 (TTG). Except for three genes (cox1, nad4 and nad5) end with the incomplete stop codon T--, all other PCGs terminated with the stop codon TAA or TAG. The whole mitogenome exhibited heavy AT nucleotide bias (80.0%). Phylogenetic analysis positioned M. hieroglyphica in a well-supported clade within the subfamily Galerucinae with Monolepta occifluvis, Monolepta sp. and Paleosepharia posticata. These results provided an important basis for further studies on mitochondrial genome and phylogenetics of Galerucinae.

Prathapanius fortis, a new genus and new species of Galerucinae from Ecuador (Coleoptera, Chrysomelidae)

Prathapanius Viswajyothi & Clark, gen. nov., is described and illustrated. The genus is placed in the chrysomelid section Diabroticites Chapuis (subtribe Diabroticina Chapuis, tribe Luperini Chapuis, subfamily Galerucinae Latreille). It is monobasic, containing Prathapanius fortis Viswajyothi & Clark, sp. nov., from Ecuador. The new genus is briefly compared with Acalymma Barber, Isotes Weise, and Zischkaita Bechyné.

The mitochondrial genomes of ladybird beetles and implications for evolution and phylogeny

Ladybirds formed the most familiar beetle group, namely the family Coccinellidae, whose internal relationships remain unclear. In particular, the subfamily relationships could not be well resolved in previous studies based on the conventional nuclear and/or mitochondrial gene fragments. In this study, we used next-generation sequencing to obtain new mitochondrial genomes (mitogenomes) from 13 species representing four coccinellid subfamilies (i.e., Coccinellinae, Epilachninae, Coccidulinae and Chilocorinae). Together with 24 existing mitogenome sequences of Cucujoidea, we conducted phylogenetic analyses to investigate the deep phylogenetic relationships in Coccinellidae, under maximum likelihood and Bayesian inference criteria. The analyses from nucleotide datasets resulted in a largely identical tree topology, where Epilachninae and Coccinellinae were monophyletic groups. The Scymninae and Coccidulinae were recovered as non-monophyletic. Amino acids differed from nucleotides in that the Epilachninae was retrieved as paraphyletic, with respect to Epilachna admirabilis. Ancestral state reconstruction suggested that the plant eating ladybird beetles arose within an aphidophagous/coccidophagous clade. In addition, three independent shifts toward coccidophagy and one shift toward mycophagy occurred in Coccinellidae.

Comparative Mitogenomic Analysis of the Eurydema Genus in the Context of Representative Pentatomidae (Hemiptera: Heteroptera) Taxa

The family Pentatomidae, the largest within the superfamily Pentatomoidae, comprises about 5,000 species; many of which are economically important pests. Although the phylogeny of Pentatomidae species has been studied using various molecular markers, their phylogenetic relationships remain controversial. Recently, mitochondrial genomes (mitogenomes) have been extensively employed to examine the phylogenetics and evolution of different insects, and in this study, we sequenced complete/near-complete mitochondrial genomes from five shield bug species of Eurydema to gain a better understanding of phylogenetic relationships in the Pentatomidae. The five mitogenomes ranged in length from 15,500 to 16,752 bp and comprised 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region. We compared mitogenomic characteristics of the Pentatomidae and constructed phylogenetic trees using Bayesian inference and maximum likelihood methods. Our results showed that gene arrangements, base composition, start/stop codons, gene overlaps, and RNA structures were conserved within the Pentatomidae and that congeneric species shared more characteristics. Saturation and heterogeneity analyses revealed that our PCGs and PCGRNA datasets were valid for phylogenetic analysis. Phylogenetic analyses showed consistent topologies based on BI and ML methods. These analyses strongly supported that Eurydema species belong to the tribe Strachiini, and formed a sister group with Pentatomini. The relationships among Eurydema species were shown to be consistent with their morphological features. (Strachiini + Pentatomini) was found to be a stable sibling of the clade comprising Cappaeini, Graphosomini, and Carpocorini. Furthermore, our results indicated that Graphosoma rubrolineatum (Heteroptera: Pentatomidae) belongs to the Pentatominae and not the Podopinae.

Identification and analysis of the complete mitochondrial genome of Phyllotreta striolata (Coleoptera, Chrysomelidae)

In this study, the complete mitochondrial genome of Phyllotreta striolata (Coleoptera, Chrysomeloidea, Chrysomelidae) was first determined. The complete genome is 15,689 bp in length. It contains 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a control region (A + T-rich region). The gene organization, nucleotide composition, and codon usage are similar to other Chrysomelidae mitogenomes. The overall nucleotide composition was 39.90% A, 35.94% T, 15.27% C, and 8.89% G, respectively. Phylogenetic analysis both highly supported that P. striolata showed a close relationship with P. undulata. The measure of complete mitogenome sequence of P. striolata will provide fundamental data for the phylogenetic and biogeographic studies of the Chrysomeloidea and Coleoptera.

Mitochondrial Architecture Rearrangements Produce Asymmetrical Nonadaptive Mutational Pressures That Subvert the Phylogenetic Reconstruction in Isopoda

The phylogeny of Isopoda, a speciose order of crustaceans, remains unresolved, with different data sets (morphological, nuclear, mitochondrial) often producing starkly incongruent phylogenetic hypotheses. We hypothesized that extreme diversity in their life histories might be causing compositional heterogeneity/heterotachy in their mitochondrial genomes, and compromising the phylogenetic reconstruction. We tested the effects of different data sets (mitochondrial, nuclear, nucleotides, amino acids, concatenated genes, individual genes, gene orders), phylogenetic algorithms (assuming data homogeneity, heterogeneity, and heterotachy), and partitioning; and found that almost all of them produced unique topologies. As we also found that mitogenomes of Asellota and two Cymothoida families (Cymothoidae and Corallanidae) possess inversed base (GC) skew patterns in comparison to other isopods, we concluded that inverted skews cause long-branch attraction phylogenetic artifacts between these taxa. These asymmetrical skews are most likely driven by multiple independent inversions of origin of replication (i.e., nonadaptive mutational pressures). Although the PhyloBayes CAT-GTR algorithm managed to attenuate some of these artifacts (and outperform partitioning), mitochondrial data have limited applicability for reconstructing the phylogeny of Isopoda. Regardless of this, our analyses allowed us to propose solutions to some unresolved phylogenetic debates, and support Asellota are the most likely candidate for the basal isopod branch. As our findings show that architectural rearrangements might produce major compositional biases even on relatively short evolutionary timescales, the implications are that proving the suitability of data via composition skew analyses should be a prerequisite for every study that aims to use mitochondrial data for phylogenetic reconstruction, even among closely related taxa.

Host plant associations in Western Palaearctic Longitarsus flea beetles (Chrysomelidae, Galerucinae, Alticini): a preliminary phylogenetic assessment

Longitarsus Latreille (Chrysomelidae, Galerucinae, Alticini) is a very large genus of phytophagous insects, with more than 700 species distributed in all zoogeographical regions. Patterns of host use have been a central topic in phytophagous insect research. In this study a first assessment is provided to test the hypothesis that host-plant association is phylogenetically conserved in Western Palaearctic Longitarsus species. Maximum Likelihood and Bayesian Inference methods were used to infer a phylogeny based on DNA sequence data from two mitochondrial genes from 52 Longitarsus species from the Western Palaearctic. In agreement with the host phylogenetic conservatism hypothesis, a strict association between most of the recovered clades and specific plant families was found, except for species associated with Boraginaceae. Low phylogenetic resolution at deep nodes limited the evaluation of whether closely related Longitarsus clades are associated with the same plant family or to closely related plant families.

Adaptation of flea beetles to Brassicaceae: host plant associations and geographic distribution of Psylliodes Latreille and Phyllotreta Chevrolat (Coleoptera, Chrysomelidae)

The cosmopolitan flea beetle genera Phyllotreta and Psylliodes (Galerucinae, Alticini) are mainly associated with host plants in the family Brassicaceae and include economically important pests of crucifer crops. In this review, the host plant associations and geographical distributions of known species in these genera are summarised from the literature, and their proposed phylogenetic relationships to other Alticini analysed from published molecular phylogenetic studies of Galerucinae. Almost all Phyllotreta species are specialised on Brassicaceae and related plant families in the order Brassicales, whereas Psylliodes species are associated with host plants in approximately 24 different plant families, and 50% are specialised to feed on Brassicaceae. The current knowledge on how Phyllotreta and Psylliodes are adapted to the characteristic chemical defence in Brassicaceae is reviewed. Based on our findings we postulate that Phyllotreta and Psylliodes colonised Brassicaceae independently from each other.

The complete mitochondrial genomes of two sibling species of camellia weevils (Coleoptera: Curculionidae) and patterns of Curculionini speciation

Complete mitochondrial genomes contain large and diverse datasets for species delineation. To better understand the divergence of the two morphologically indistinguishable weevil species in Curculionini, we first sequenced and compared their complete mitochondrial genomes. The complete mitochondrial genomes of Curculio chinensis and Curculio sp. were 19,713 bp with an A + T content of 76.61% and 19,216 bp with an A + T content of 76.85%, respectively. All 37 of the typical mitochondrial genes were determined in both species. The 13 protein sequences of the two species shared high homology (about 90%) except for ATP8 (73.08%). The differences in secondary structure of ATP8 were the number of possible proteins and nucleic acid binding sites. There were 22 and 15 mismatched base-pairs in the tRNA secondary structures from C. chinensis and Curculio sp., respectively. Maximum Likelihood and Bayesian analyses indicated that Curculio sp. is a novel species closely related to C. chinensis. The divergence time estimation suggests that Cryptorhynchinae and Curculionini lines diverged in the Cenozoic Period, while C. chinensis and Curculio sp. diverged at 6.7079 (95% CI 5-13) Mya. This study demonstrates the utility of using complete mitochondrial gene sets for phylogenetic analysis and enhances our understanding of the genetic basis for the evolution of the Curculionini.

One Pathway Is Not Enough: The Cabbage Stem Flea Beetle Psylliodes chrysocephala Uses Multiple Strategies to Overcome the Glucosinolate-Myrosinase Defense in Its Host Plants

The cabbage stem flea beetle (Psylliodes chrysocephala) is a key pest of oilseed rape in Europe, and is specialized to feed on Brassicaceae plants armed with the glucosinolate-myrosinase defense system. Upon tissue damage, the β-thioglucosidase enzyme myrosinase hydrolyzes glucosinolates (GLS) to form toxic isothiocyanates (ITCs) which deter non-adapted herbivores. Here, we show that P. chrysocephala selectively sequester GLS from their host plants and store these throughout their life cycle. In addition, P. chrysocephala metabolize GLS to desulfo-GLS, which implies the evolution of GLS sulfatase activity in this specialist. To assess whether P. chrysocephala can largely prevent GLS hydrolysis in ingested plant tissue by sequestration and desulfation, we analyzed the metabolic fate of 4-methylsulfinylbutyl (4MSOB) GLS in adults. Surprisingly, intact and desulfo-GLS together accounted for the metabolic fate of only 26% of the total ingested GLS in P. chrysocephala, indicating that most ingested GLS are nevertheless activated by the plant myrosinase. The presence of 4MSOB-ITC and the corresponding nitrile in feces extracts confirmed the activation of ingested GLS, but the detected amounts of unmetabolized ITCs were low. P. chrysocephala partially detoxifies ITCs by conjugation with glutathione via the conserved mercapturic acid pathway. In addition to known products of the mercapturic acid pathway, we identified two previously unknown cyclic metabolites derived from the cysteine-conjugate of 4MSOB-ITC. In summary, the cabbage stem flea beetle avoids ITC formation by specialized strategies, but also relies on and extends the conserved mercapturic acid pathway to prevent toxicity of formed ITCs.

The Mitochondrial Genomes of Phytophagous Scarab Beetles and Systematic Implications

In this study, we newly sequenced five mitogenomes of representatives of phytophagous scarab beetles (Coleoptera: Scarabaeidae) by using next-generation sequencing technology. Two species have complete (or nearly complete) mitogenome sequences, namely Popillia mutans Newman (Coleoptera: Scarabaeidae) and Holotrichia oblita Faldermann (Coleoptera: Scarabaeidae). The remaining three species have the partial mitogenomes, and the missing genes are mainly located adjacent to the control region. The complete (or nearly complete) mitogenomes have the same genome structure as most of the existing Scarabaeidae mitogenomes. We conducted phylogenetic analyses together with 24 published mitogenomes of Scarabaeoidea. The results supported a basal split of coprophagous and phytophagous Scarabaeidae. The subfamily Sericinae was recovered as sister to all other phytophagous scarab beetles. All analyses supported a non-monophyletic Melolonthinae, which included two different non-sister clades. The Cetoniinae was recovered as sister to a clade including Rutelinae and Dynastinae. Although the Rutelinae was rendered paraphyletic by Dynastinae in the Bayesian trees inferred under the site-heterogeneous CAT-GTR or CAT-MTART model, discordant patterns were given in some of ML trees estimated using the homogeneous GTR model.

How many genera and species of Galerucinaes. str. do we know? Updated statistics (Coleoptera, Chrysomelidae)

Galerucinae s. str. is a rich group of leaf beetles. A new, up-to date checklist of Galerucinae genera in the world is provided, including the number of valid species of each genus. Genera and species were counted in literature published before the end of 2016. In summary, 7145 species (7132 recent, 13 fossils) and 192 subspecies from 543 genera (542 recent, 1 fossil) were quantified in Galerucinae s. str. In comparison with the previous catalogue of worldwide Galerucinae (Wilcox 1971-1973), an additional 91 valid genera, 1341 valid species (1337 recent, 4 fossils) and 38 subspecies have been published; 43 genera were synonymized, four genera were transferred into Alticini, two subgenera were elevated to genus rank, and one genus was downgraded to subgenus rank. The updated list of references to taxonomic publications on Galerucinae s. str. from the period 1971-2016 is provided.