Structural alignment of 18S and 28S rDNA sequences provides insights into phylogeny of Phytophaga (Coleoptera: Curculionoidea and Chrysomeloidea)

The sperm ultrastructure of Spermophagus kuesteri Schilsky, 1905 (Chrysomelidae, Bruchinae) and considerations on the relationships of Cucujiformia superfamilies

The sperm ultrastructure of the bean-weevil Spermophagus kuesteri (Bruchinae) was studied to verify the congruence of the new position of the subfamily within Chrysomelidae. The results indicated a positive answer to the question supporting a close relationship between Chrysomelidae and Curculionidae, a finding confirmed also by molecular data. Moreover, the sperm morphology of Divales cinctus, a member of Melyridae (Cleroidea) allowed to confirm the different sperm organization between members of this superfamily and Phytophaga (Chrysomeloidea + Curculionoidea). While studying the spermiogenesis of S. kuesteri, some sperm cysts showed aberrant cells provided with two flagella in the same plasma membrane. These aberrant sperm could be the result, during early spermiogenesis, of irregular processes involving the canal rings between spermatids.

Three chemosensory proteins enriched in antennae and tarsi of Rhaphuma horsfieldi differentially contribute to the binding of insecticides

Antennae and legs (primarily the tarsal segments) of insects are the foremost sensory organs that contact a diverse range of toxic chemicals including insecticides. Binding proteins expressed in the two tissues are potential molecular candidates serving as the binding and sequestering of insecticides, like chemosensory proteins (CSPs). Insect CSPs endowed with multiple roles have been suggested to participate in insecticide resistance, focusing mainly on moths, aphids and mosquitos. Yet, the molecular underpinnings underlying the interactions of cerambycid CSPs and insecticides remain unexplored. Here, we present binding properties of three antenna- and tarsus-enriched RhorCSPs (RhorCSP1, CSP2 and CSP3) in Rhaphuma horsfieldi to eight insecticide classes totaling 15 chemicals. From the transcriptome of this beetle, totally 16 CSP-coding genes were found, with seven full-length sequences. In phylogeny, these RhorCSPs were distributed dispersedly in different clades. Expression profiles revealed the abundant expression of RhorCSP1, CSP2 and CSP3 in antennae and tarsi, thus as representatives for studying the protein-insecticide interactions. Binding assays showed that the three RhorCSPs were tuned differentially to insecticides but exhibited the highest affinities with hexaflumuron, chlorpyrifos and rotenone (dissociation constants <13 μM). In particular, RhorCSP3 could interact strongly with 10 of tested insecticides, of which four residues (Tyr25, Phe42, Val65 and Phe68) contributed significantly to the binding of six, four, three and four ligands, respectively. Of these, the binding of four mutated RhorCSP3s to a botanical insecticide rotenone was significantly weakened compared to the wildtype protein. Furthermore, we also evidenced that RhorCSP3 was a broadly-tuned carrier protein in response to a wide variety of plant odorants outside insecticides. Altogether, our findings shed light on different binding mechanisms and odorant-tuning profiles of three RhorCSPs in R. horsfieldi and identify key residues of the RhorCSP3-insecticide interactions.

Taxonomy, Biology, Symbionts, Omics, and Management of Rhynchophorus Palm Weevils (Coleoptera: Curculionidae: Dryophthorinae)

Palm weevils, Rhynchophorus spp., are destructive pests of native, ornamental, and agricultural palm species. Of the 10 recognized species, two of the most injurious species, Rhynchophorus ferrugineus and Rhynchophorus palmarum, both of which have spread beyond their native range, are the best studied. Due to its greater global spread and damage to edible date industries in the Middle East, R. ferrugineus has received more research interest. Integrated pest management programs utilize traps baited with aggregation pheromone, removal of infested palms, and insecticides. However, weevil control is costly, development of resistance to insecticides is problematic, and program efficacy can be impaired because early detection of infestations is difficult. The genome of R. ferrugineus has been sequenced, and omics research is providing insight into pheromone communication and changes in volatile and metabolism profiles of weevil-infested palms. We outline how such developments could lead to new control strategies and early detection tools.

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

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

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.

Genome sequencing provides insights into the evolution of gene families encoding plant cell wall-degrading enzymes in longhorned beetles

With more than 36,000 species, the longhorned beetles (family Cerambycidae) are a mega-diverse lineage of mostly xylophagous insects, all of which are represented by the sole sequenced genome of the Asian longhorned beetle (Anoplophora glabripennis; Lamiinae). Their successful radiation has been linked to their ability to degrade plant cell wall components using a range of so-called plant cell wall-degrading enzymes (PCWDEs). Our previous analysis of larval gut transcriptomes demonstrated that cerambycid beetles horizontally acquired genes encoding PCWDEs from various microbial donors; these genes evolved through multiple duplication events to form gene families. To gain further insights into the evolution of these gene families during the Cerambycidae radiation, we assembled draft genomes for four beetle species belonging to three subfamilies using long-read nanopore sequencing. All the PCWDE-encoding genes we annotated from the corresponding larval gut transcriptomes were present in these draft genomes. We confirmed that the newly discovered horizontally acquired glycoside hydrolase family 7 (GH7), subfamily 26 of GH43 (GH43_26), and GH53 (all of which are absent from the A. glabripennis genome) were indeed encoded by these beetles' genome. Most of the PCWDE-encoding genes of bacterial origin gained introns after their transfer into the beetle genome. Altogether, we show that draft genome assemblies generated from nanopore long-reads offer meaningful information to the study of the evolution of gene families in insects. We anticipate that our data will support studies aiming to better understand the biology of the Cerambycidae and other beetles in general.

Phylogenomics of weevils revisited: data curation and modelling compositional heterogeneity

Weevils represent one of the most prolific radiations of beetles and the most diverse group of herbivores on land. The phylogeny of weevils (Curculionoidea) has received extensive attention, and a largely satisfactory framework for their interfamilial relationships has been established. However, a recent phylogenomic study of Curculionoidea based on anchored hybrid enrichment (AHE) data yielded an abnormal placement for the family Belidae (strongly supported as sister to Nemonychidae + Anthribidae). Here we reanalyse the genome-scale AHE data for Curculionoidea using various models of molecular evolution and data filtering methods to mitigate anticipated systematic errors and reduce compositional heterogeneity. When analysed with the infinite mixture model CAT-GTR or using appropriately filtered datasets, Belidae are always recovered as sister to the clade (Attelabidae, (Caridae, (Brentidae, Curculionidae))), which is congruent with studies based on morphology and other sources of molecular data. Although the relationships of the 'higher Curculionidae' remain challenging to resolve, we provide a consistent and robust backbone phylogeny of weevils. Our extensive analyses emphasize the significance of data curation and modelling across-site compositional heterogeneity in phylogenomic studies.

Metabolic novelty originating from horizontal gene transfer is essential for leaf beetle survival

Horizontal gene transfer (HGT) provides an evolutionary shortcut for recipient organisms to gain novel functions. Although reports of HGT in higher eukaryotes are rapidly accumulating, in most cases the evolutionary trajectory, metabolic integration, and ecological relevance of acquired genes remain unclear. Plant cell wall degradation by HGT-derived enzymes is widespread in herbivorous insect lineages. Pectin is an abundant polysaccharide in the walls of growing parts of plants. We investigated the significance of horizontally acquired pectin-digesting polygalacturonases (PGs) of the leaf beetle Phaedon cochleariae. Using a CRISPR/Cas9-guided gene knockout approach, we generated a triple knockout and a quadruple PG-null mutant in order to investigate the enzymatic, biological, and ecological effects. We found that pectin-digestion 1) is exclusively linked to the horizontally acquired PGs from fungi, 2) became fixed in the host genome by gene duplication leading to functional redundancy, 3) compensates for nutrient-poor diet by making the nutritious cell contents more accessible, and 4) facilitates the beetles development and survival. Our analysis highlights the selective advantage PGs provide to herbivorous insects and demonstrate the impact of HGT on the evolutionary success of leaf-feeding beetles, major contributors to species diversity.

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.

Phylogenetic analysis of the family Megalopodidae (Coleoptera: Chrysomeloidea): better taxon-sampling facilitates detection of new relationships and new taxa

The systematics of Megalopodidae is not adequately known, in spite of it being a relatively small group of phytophagous beetles. The first phylogenetic analysis of Megalopodidae with a comprehensive generic representation (25 genera of 30 described, 10 subgenera and 77 species) is undertaken. A parsimony analysis under equal and implied weights was carried out based on 147 adult and larval morphological characters. Subfamilies Palophaginae and Zeugophorinae were recovered as monophyletic, by contrast with Megalopodinae, which proved to be paraphyletic. Atelederinae are proposed as a new subfamily. Also, three tribes and three subtribes within Megalopodinae are proposed: Leucasteini trib.n., Sphondyliini trib.n. and Megalopodini, the latter including Macrolophina subtrib.n., Temnaspidina subtrib.n. and Megalopodina. The genera Macrolopha, Kuilua, Poecilomorpha, Temnaspis, Antonaria, Agathomerus, Megalopus and Bothromegalopus were recovered as non-monophyletic. New delimitations of the polyphyletic genera Poecilomorpha and Macrolopha are proposed, Clythraxeloma is resurrected, and the subgenera of Agathomerus are suppressed. The following new combinations are proposed: Kuilua apicata (Fairmaire), K. nyassae (Jacoby), Poecilomorpha cribricollis (Pic), P. minuta (Pic), Clythraxeloma assamensis (Jacoby), C. bipartita (Lacordaeri), C. discolineata (Pic), C. downesii (Baly), C. gerstaeckeri (Westwood), C. laosensis (Pic), C. maculata (Pic), C. mouhoti (Baly), C. nigrocyanea (Motschulsky), C. pretiosa (Reineck), Temnaspis tricoloripes (Pic) and Barticaria faciatus (Dalman). Clythraxeloma cyanipennis Kraatz is a restored combination. Distribution patterns of Megalopodidae largely conform to the breakup of Gondwanaland, with its main clades having particular distributions: Andean-Australian (Palophaginae), Ethiopian (Leucasteini, Sphondyliini, and Macrolophina), Neotropical (Ateledrinae and Megalopodina) and Ethiopian-Oriental-Palaearctic (Temnaspidina the result of a secondary expansion. Zeugophorinae present a worldwide distribution, except for the Neotropical and Andean regions, which may be the result of geodispersal. The findings of the present study also shed light on groups with taxonomic issues, where phylogenetic analyses are strongly needed.

Comparative Mitogenomic Analysis of Two Longhorn Beetles (Coleoptera: Cerambycidae: Lamiinae) with Preliminary Investigation into Phylogenetic Relationships of Tribes of Lamiinae

The subfamily Lamiinae is the most taxonomically diverse subfamily of Cerambycidae, but relationships between tribes of Lamiinae are still unresolved. In order to study the characteristics of the mitogenomes of Lamiinae and the tribal-level phylogenetic relationships, we sequenced the mitogenomes of two species representing two tribes, Agapanthia amurensis (Agapanthiini) and Moechotypa diphysis (Ceroplesini), with a total length of 15,512 bp and 15,493 bp, respectively. The gene arrangements of these two new mitogenomes were consistent with the inferred ancestral insect mitogenomes. Each species contained 37 typical mitochondrial genes and a control region (A + T-rich region), including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs). All PCGs initiated with the standard start codon ATN, and terminated with the complete stop codons of TAA and TAG, or incomplete stop codon T. All tRNAs could be folded into a clover-leaf secondary structure except for trnS1, in which the dihydrouridine (DHU) arm was reduced. Moreover, we studied the phylogenetic relationships between some tribes of Lamiinae based in mitochondrial PCGs in nucleotides; our results show that the relationships were as follows: (Onciderini + ((Apomecynini + Acanthocinini) + ((Ceroplesini + Agapanthiini) + ((Mesosini + Pteropliini) + ((Dorcaschematini + (Saperdini 1 + (Phytoeciini + Saperdini 2))) + (Batocerini + Lamiini)))))).

First phylogenetic analysis of Dryophthorinae (Coleoptera, Curculionidae) based on structural alignment of ribosomal DNA reveals Cenozoic diversification

Dryophthorinae is an economically important, ecologically distinct, and ubiquitous monophyletic group of pantropical weevils with more than 1,200 species in 153 genera. This study provides the first comprehensive phylogeny of the group with the aim to provide insights into the process and timing of diversification of phytophagous insects, inform classification and facilitate predictions. The taxon sampling is the most extensive to date and includes representatives of all five dryophthorine tribes and all but one subtribe. The phylogeny is based on secondary structural alignment of 18S and 28S rRNA totaling 3,764 nucleotides analyzed under Bayesian and maximum likelihood inference. We used a fossil-calibrated relaxed clock model with two approaches, node-dating and fossilized birth-death models, to estimate divergence times for the subfamily. All tribes except the species-rich Rhynchophorini were found to be monophyletic, but higher support is required to ascertain the paraphyly of Rhynchophorini with more confidence. Nephius is closely related to Dryophthorini and Stromboscerini, and there is strong evidence for paraphyly of Sphenophorina. We find a large gap between the divergence of Dryophthorinae from their sister group Platypodinae in the Jurassic-Cretaceous boundary and the diversification of extant species in the Cenozoic, highlighting the role of coevolution with angiosperms in this group.

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))).

Analyzing the Substrate Specificity of a Class of Long-Horned-Beetle-Derived Xylanases by Using Synthetic Arabinoxylan Oligo- and Polysaccharides

Xylophagous long-horned beetles thrive in challenging environments. To access nutrients, they secrete plant-cell-wall-degrading enzymes in their gut fluid; among them are cellulases of the subfamily 2 of glycoside hydrolase family 5 (GH5_2). Recently, we discovered that several beetle-derived GH5_2s use xylan as a substrate instead of cellulose, which is unusual for this family of enzymes. Here, we analyze the substrate specificity of a GH5_2 xylanase from the beetle Apriona japonica (AJAGH5_2-1) using commercially available substrates and synthetic arabinoxylan oligo- and polysaccharides. We demonstrate that AJAGH5_2-1 processes arabinoxylan polysaccharides in a manner distinct from classical xylanase families such as GH10 and GH11. AJAGH5_2-1 is active on long oligosaccharides and cleaves at the non-reducing end of a substituted xylose residue (position +1) only if: 1) three xylose residues are present upstream and downstream of the cleavage site, and 2) xylose residues at positions -1, -2, +2 and +3 are not substituted.

Plants use identical inhibitors to protect their cell wall pectin against microbes and insects

As fundamentally different as phytopathogenic microbes and herbivorous insects are, they enjoy plant-based diets. Hence, they encounter similar challenges to acquire nutrients. Both microbes and beetles possess polygalacturonases (PGs) that hydrolyze the plant cell wall polysaccharide pectin. Countering these threats, plant proteins inhibit PGs of microbes, thereby lowering their infection rate. Whether PG-inhibiting proteins (PGIPs) play a role in defense against herbivorous beetles is unknown. To investigate the significance of PGIPs in insect-plant interactions, feeding assays with the leaf beetle Phaedon cochleariae on Arabidopsis thaliana pgip mutants were performed. Fitness was increased when larvae were fed on mutant plants compared to wild-type plants. Moreover, PG activity was higher, although PG genes were downregulated in larvae fed on PGIP-deficient plants, strongly suggesting that PGIPs impair PG activity. As low PG activity resulted in delayed larval growth, our data provide the first in vivo correlative evidence that PGIPs act as defense against insects.

Molecular phylogenetic assessment of the tribal classification of Lamiinae (Coleoptera: Cerambycidae)

Lamiinae is the most diverse subfamily of longhorned beetles, with about 20,000 described species classified into 80 tribes. Most of the tribes of Lamiinae were proposed during the 19th century and the suprageneric classification of the subfamily has never been assessed under phylogenetic criteria. In this study, we present the first tribal-level phylogeny of Lamiinae, inferred from 130 terminals (representing 46 tribes, prioritizing generic type species of the tribes) and fragments of two mitochondrial and three nuclear markers (cox1, rrnL, Wg, CPS and LSU; 5,024 aligned positions in total). Analyses were performed under Maximum Likelihood and Bayesian methods based on two datasets: a dataset including all taxa available for the study, and a reduced dataset with 111 terminals where taxa only contributing with mitochondrial markers were excluded from the matrix. The monophyly of Lamiinae was corroborated in three of the four analyses and 11 of the 35 tribes with more than one species represented in the analyses were consistently recovered as monophyletic. However, 15 tribes were not retrieved as monophyletic, requiring a revision of their boundaries: Acanthocinini, Acanthoderini, Agapanthiini, Apomecynini, Desmiphorini, Dorcaschematini, Enicodini, Hemilophini, Monochamini, Onciderini, Parmenini, Phytoeciini, Pogonocherini, Pteropliini and Saperdini. Based on these results, when strong support values for paraphyly were recovered, we argue a number of tribe synonymies, including Moneilemini as synonym of Acanthocinini; Onocephalini of Onciderini; Dorcadionini, Gnomini, Monochamini and Rhodopinini of Lamiini; and Obereini and Phytoeciini of Saperdini. Other taxonomic changes proposed in this study based on the criterion of monophyly and supported by morphological characters include the transfer of Tricondyloides and Stenellipsis to Enicodini, and of Dylobolus stat. rest., which is removed as subgenus of Mecas and restituted as genus, to Hemilophini. Furthermore, our analyses suggest that Ostedes and Neohoplonotus should be removed from Acanthocinini and Parmenini, respectively, and Colobotheini should be redefined to encompass several genera currently placed in Acanthocinini.

The mitochondrial genome of Apion squamigerum (Coleoptera, Curculionoidea, Brentidae) and the phylogenetic implications

In this article, we present the nearly complete mitochondrial genome (mitogenome) of the weevil beetle Apion squamigerum (Curculionoidea, Brentidae), assembled using data from Illumina next generation sequencing (NGS). This mitogenome was found to be very large, with the total length of 18,562 bp. Two trnM genes were identified. A large non-coding intergenic spacer spanning 1,949 bp occurred between trnI and trnM2. Combined with 111 existing weevil mitogenomes, we conducted phylogenetic reconstructions based on various datasets under maximum likelihood and Bayesian inference. Firstly, phylogenetic analyses robustly supported a sister group of A. squamigerum and Rhopalapion longirostre, namely, that two species of Apioninae (Brentidae) formed a clade. Within the entire Curculionoidea, the Nemonychidae diverged firstly, following the families Anthribidae and Attelabidae. In addition, a large clade comprising the sister families Brentidae and Curculionidae was strongly supported in all trees.

Pectin Digestion in Herbivorous Beetles: Impact of Pseudoenzymes Exceeds That of Their Active Counterparts

Many protein families harbor pseudoenzymes that have lost the catalytic function of their enzymatically active counterparts. Assigning alternative function and importance to these proteins is challenging. Because the evolution toward pseudoenzymes is driven by gene duplication, they often accumulate in multigene families. Plant cell wall-degrading enzymes (PCWDEs) are prominent examples of expanded gene families. The pectolytic glycoside hydrolase family 28 (GH28) allows herbivorous insects to break down the PCW polysaccharide pectin. GH28 in the Phytophaga clade of beetles contains many active enzymes but also many inactive counterparts. Using functional characterization, gene silencing, global transcriptome analyses, and recordings of life history traits, we found that not only catalytically active but also inactive GH28 proteins are part of the same pectin-digesting pathway. The robustness and plasticity of this pathway and thus its importance for the beetle is supported by extremely high steady-state expression levels and counter-regulatory mechanisms. Unexpectedly, the impact of pseudoenzymes on the pectin-digesting pathway in Phytophaga beetles exceeds even the influence of their active counterparts, such as a lowered efficiency of food-to-energy conversion and a prolongation of the developmental period.

Functional diversification of horizontally acquired glycoside hydrolase family 45 (GH45) proteins in Phytophaga beetles

BACKGROUND

Cellulose, a major polysaccharide of the plant cell wall, consists of β-1,4-linked glucose moieties forming a molecular network recalcitrant to enzymatic breakdown. Although cellulose is potentially a rich source of energy, the ability to degrade it is rare in animals and was believed to be present only in cellulolytic microbes. Recently, it has become clear that some animals encode endogenous cellulases belonging to several glycoside hydrolase families (GHs), including GH45. GH45s are distributed patchily among the Metazoa and, in insects, are encoded only by the genomes of Phytophaga beetles. This study aims to understand both the enzymatic functions and the evolutionary history of GH45s in these beetles.

RESULTS

To this end, we biochemically assessed the enzymatic activities of 37 GH45s derived from five species of Phytophaga beetles and discovered that beetle-derived GH45s degrade three different substrates: amorphous cellulose, xyloglucan and glucomannan. Our phylogenetic and gene structure analyses indicate that at least one gene encoding a putative cellulolytic GH45 was present in the last common ancestor of the Phytophaga, and that GH45 xyloglucanases evolved several times independently in these beetles. The most closely related clade to Phytophaga GH45s was composed of fungal sequences, suggesting this GH family was acquired by horizontal gene transfer from fungi. Besides the insects, other arthropod GH45s do not share a common origin and appear to have emerged at least three times independently.

CONCLUSION

The rise of functional innovation from gene duplication events has been a fundamental process in the evolution of GH45s in Phytophaga beetles. Both, enzymatic activity and ancestral origin suggest that GH45s were likely an essential prerequisite for the adaptation allowing Phytophaga beetles to feed on plants.

Phylogenomic Data Yield New and Robust Insights into the Phylogeny and Evolution of Weevils

The phylogeny and evolution of weevils (the beetle superfamily Curculionoidea) has been extensively studied, but many relationships, especially in the large family Curculionidae (true weevils; > 50,000 species), remain uncertain. We used phylogenomic methods to obtain DNA sequences from 522 protein-coding genes for representatives of all families of weevils and all subfamilies of Curculionidae. Most of our phylogenomic results had strong statistical support, and the inferred relationships were generally congruent with those reported in previous studies, but with some interesting exceptions. Notably, the backbone relationships of the weevil phylogeny were consistently strongly supported, and the former Nemonychidae (pine flower snout beetles) were polyphyletic, with the subfamily Cimberidinae (here elevated to Cimberididae) placed as sister group of all other weevils. The clade comprising the sister families Brentidae (straight-snouted weevils) and Curculionidae was maximally supported and the composition of both families was firmly established. The contributions of substitution modeling, codon usage and/or mutational bias to differences between trees reconstructed from amino acid and nucleotide sequences were explored. A reconstructed timetree for weevils is consistent with a Mesozoic radiation of gymnosperm-associated taxa to form most extant families and diversification of Curculionidae alongside flowering plants-first monocots, then other groups-beginning in the Cretaceous.

Cellulose degradation in Gastrophysa viridula (Coleoptera: Chrysomelidae): functional characterization of two CAZymes belonging to glycoside hydrolase family 45 reveals a novel enzymatic activity

Cellulose is a major component of the primary and secondary cell walls in plants. Cellulose is considered to be the most abundant biopolymer on Earth and represents a large potential source of metabolic energy. Yet, cellulose degradation is rare and mostly restricted to cellulolytic microorganisms. Recently, various metazoans, including leaf beetles, have been found to encode their own cellulases, giving them the ability to degrade cellulose independently of cellulolytic symbionts. Here, we analyzed the cellulosic capacity of the leaf beetle Gastrophysa viridula, which typically feeds on Rumex plants. We identified three putative cellulases member of two glycoside hydrolase (GH) families, namely GH45 and GH9. Using heterologous expression and functional assays, we demonstrated that both GH45 proteins are active enzymes, in contrast to the GH9 protein. One GH45 protein acted on amorphous cellulose as an endo-β-1,4-glucanase, whereas the other evolved to become an endo-β-1,4-xyloglucanase. We successfully knocked down the expression of both GH45 genes using RNAi, but no changes in weight gain or mortality were observed compared to control insects. Our data indicated that the breakdown of these polysaccharides in G. viridula may facilitate access to plant cell content, which is rich in nitrogen and simple sugars.

A Combined Molecular and Morphological Approach to Explore the Higher Phylogeny of Entimine Weevils (Coleoptera: Curculionidae), with Special Reference to South American Taxa

The Entiminae are broad-nosed weevils constituting the most diverse subfamily of Curculionidae, with over 50 tribes. We performed Bayesian and Maximum Parsimony combined phylogenetic analyses with the main objective of testing higher-level relationships and the naturalness of the major Neotropical and Southern South American (Patagonia and Andes) tribes, including some members from other regions. We compiled a data matrix of 67 terminal units with 63 Entiminae species, as well as four outgroup taxa from Cyclominae, by 3522 molecular (from nuclear 18S rDNA and 28S rDNA, and mitochondrial 16S rDNA and COI gene sequences) and 70 morphological characters. The resulting trees recover a clade Entiminae with a monophyletic Cylydrorhinini and Premnotrypes branching off early. The tree resulting from parsimony analysis shows a clade of Leptopiini from the Australian region and another clade including taxa mainly distributed in the Palaearctic and Neotropical regions, but in the Bayesian tree the South American and Australian Leptopiini are grouped together. The mainly Palaearctic Entiminae (e.g., Brachyderini, Laparocerini, Otiorhynchini, Peritelini, Polydrusini, Phyllobiini and Sciaphylini) form a subclade separated from Southern Hemisphere taxa. Among the latter, the well-supported Naupactini are the sister group of the South American Tanymecini, excluding Platyaspistes, herein transferred to Leptopiini (new placement). Another well-justified clade is Eustylini–Geonemini, which also includes the enigmatic Galapagonotus, and the genus Artipus, thus corroborating its recent exclusion from Naupactini.

Morphological and Molecular Perspectives on the Phylogeny, Evolution, and Classification of Weevils (Coleoptera: Curculionoidea): Proceedings from the 2016 International Weevil Meeting

The 2016 International Weevil Meeting was held immediately after the International Congress of Entomology (ICE). It built on the topics and content of the 2016 ICE weevil symposium Phylogeny and Evolution of Weevils (Coleoptera: Curculionoidea): A Symposium in Honor of Dr. Guillermo "Willy” Kuschel. Beyond catalyzing research and collaboration, the meeting was intended to serve as a forum for identifying priorities and goals for those who study weevils. The meeting consisted of 46 invited and contributed lectures, discussion sessions and introductory remarks presented by 23 speakers along with eight contributed research posters. These were organized into three convened sessions, each lasting one day: (1) weevil morphology; (2) weevil fossils, biogeography and host/habitat associations; and (3) molecular phylogenetics and classification of weevils. Some of the topics covered included the 1K Weevils Project, major morphological character systems of adult and larval weevils, weevil morphological terminology, prospects for future morphological character discovery, phylogenetic analysis of morphological character data, the current status of weevil molecular phylogenetics and evolution, resources available for phylogenetic and comparative genomic studies of weevils, the weevil fossil record, weevil biogeography and evolution, weevil host plants, evolutionary development of the weevil rostrum, resources available for weevil identification and the current status of and challenges in weevil classification.

Topological support and data quality can only be assessed through multiple tests in reviewing Blattodea phylogeny

Assessing support for molecular phylogenies is difficult because the data is heterogeneous in quality and overwhelming in quantity. Traditionally, node support values (bootstrap frequency, Bayesian posterior probability) are used to assess confidence in tree topologies. Other analyses to assess the quality of phylogenetic data (e.g. Lento plots, saturation plots, trait consistency) and the resulting phylogenetic trees (e.g. internode certainty, parameter permutation tests, topological tests) exist but are rarely applied. Here we argue that a single qualitative analysis is insufficient to assess support of a phylogenetic hypothesis and relate data quality to tree quality. We use six molecular markers to infer the phylogeny of Blattodea and apply various tests to assess relationship support, locus quality, and the relationship between the two. We use internode-certainty calculations in conjunction with bootstrap scores, alignment permutations, and an approximately unbiased (AU) test to assess if the molecular data unambiguously support the phylogenetic relationships found. Our results show higher support for the position of Lamproblattidae, high support for the termite phylogeny, and low support for the position of Anaplectidae, Corydioidea and phylogeny of Blaberoidea. We use Lento plots in conjunction with mutation-saturation plots, calculations of locus homoplasy to assess locus quality, identify long branch attraction, and decide if the tree's relationships are the result of data biases. We conclude that multiple tests and metrics need to be taken into account to assess tree support and data robustness.

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

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

Comparative morphology of the spermatheca in Megalopodidae (Coleoptera, Chrysomeloidea)

The spermatheca is an organ that stores and maintains viability of sperm until fertilization. It has an important role in copulation and oviposition, and it is highly informative in species delimitation. Here, we present a comparative study of the spermathecal morphology in the coleopteran family Megalopodidae. The spermathecae of 34 species, representing 13 genera and all three subfamilies, were studied. Illustrations are newly provided for all species, except in 14 cases in which illustrations were reproduced from previously published literature. Our results show that each subfamily of Megalopodidae can be effectively differentiated based on the particular spermathecal anatomy. In addition, the spermathecal anatomy presents a range of variation within each subfamily, useful for diagnosing species and, in some cases, identifying groups of genera. For instance, the “American group” is thus recognized in this study.

Evolution and functional characterization of CAZymes belonging to subfamily 10 of glycoside hydrolase family 5 (GH5_10) in two species of phytophagous beetles

Hemicelluloses, such as xyloglucan, xylan and mannans, consist of a heterogeneous array of plant-derived polysaccharides that form the plant cell wall. These polysaccharides differ from each other in their structure and physiochemical properties, but they share a β-(1,4)-linked sugar backbone. Hemicelluloses can be hydrolyzed by plant-cell-wall-degrading enzymes (PCWDEs), which are widely distributed in phytopathogenic microbes. Recently, it has become apparent that phytophagous beetles also produce their own PCWDEs. Our previous work identified genes encoding putative mannanases belonging to the subfamily 10 of glycoside hydrolase (GH) family 5 (GH5_10) in the genomes of the leaf beetle, Gastrophysa viridula (Chrysomelidae, Chrysomelinae; one gene), and of the bean beetle, Callosobruchus maculatus (Chrysomelidae, Bruchinae; four genes). In contrast to proteins from other GH5 subfamilies, GH5_10 proteins are patchily distributed within the tree of life and have so far hardly been investigated. We addressed the following questions: Are beetle-derived GH5_10s active PCWDEs? How did they evolve? What is their physiological function? Using heterologous protein expression and enzymatic assays, we show that the G. viridula GH5_10 protein is an endo-β-1,4-mannanase. We also demonstrate that only one out of four C. maculatus GH5_10 proteins is an endo-β-1,4-mannanase, which has additional activity on carboxymethyl cellulose. Unexpectedly, another C. maculatus GH5_10 protein has evolved to use xylan instead of mannans as a substrate. RNAi experiments in G. viridula indicate (i) that the sole GH5_10 protein is responsible for breaking down mannans in the gut and (ii) that this breakdown may rather be accessory and may facilitate access to plant cell content, which is rich in nitrogen and simple sugars. Phylogenetic analyses indicate that coleopteran-derived GH5_10 proteins cluster together with Chelicerata-derived ones. Interestingly, other insect-derived GH5_10 proteins cluster elsewhere, suggesting insects have several independent evolutionary origins.

Horizontal Gene Transfer Contributes to the Evolution of Arthropod Herbivory

Within animals, evolutionary transition toward herbivory is severely limited by the hostile characteristics of plants. Arthropods have nonetheless counteracted many nutritional and defensive barriers imposed by plants and are currently considered as the most successful animal herbivores in terrestrial ecosystems. We gather a body of evidence showing that genomes of various plant feeding insects and mites possess genes whose presence can only be explained by horizontal gene transfer (HGT). HGT is the asexual transmission of genetic information between reproductively isolated species. Although HGT is known to have great adaptive significance in prokaryotes, its impact on eukaryotic evolution remains obscure. Here, we show that laterally transferred genes into arthropods underpin many adaptations to phytophagy, including efficient assimilation and detoxification of plant produced metabolites. Horizontally acquired genes and the traits they encode often functionally diversify within arthropod recipients, enabling the colonization of more host plant species and organs. We demonstrate that HGT can drive metazoan evolution by uncovering its prominent role in the adaptations of arthropods to exploit plants.

Family-Level Sampling of Mitochondrial Genomes in Coleoptera: Compositional Heterogeneity and Phylogenetics

Mitochondrial genomes are readily sequenced with recent technology and thus evolutionary lineages can be densely sampled. This permits better phylogenetic estimates and assessment of potential biases resulting from heterogeneity in nucleotide composition and rate of change. We gathered 245 mitochondrial sequences for the Coleoptera representing all 4 suborders, 15 superfamilies of Polyphaga, and altogether 97 families, including 159 newly sequenced full or partial mitogenomes. Compositional heterogeneity greatly affected 3rd codon positions, and to a lesser extent the 1st and 2nd positions, even after RY coding. Heterogeneity also affected the encoded protein sequence, in particular in the nad2, nad4, nad5, and nad6 genes. Credible tree topologies were obtained with the nhPhyML ("nonhomogeneous") algorithm implementing a model for branch-specific equilibrium frequencies. Likelihood searches using RAxML were improved by data partitioning by gene and codon position. Finally, the PhyloBayes software, which allows different substitution processes for amino acid replacement at various sites, produced a tree that best matched known higher level taxa and defined basal relationships in Coleoptera. After rooting with Neuropterida outgroups, suborder relationships were resolved as (Polyphaga (Myxophaga (Archostemata + Adephaga))). The infraorder relationships in Polyphaga were (Scirtiformia (Elateriformia ((Staphyliniformia + Scarabaeiformia) (Bostrichiformia (Cucujiformia))))). Polyphagan superfamilies were recovered as monophyla except Staphylinoidea (paraphyletic for Scarabaeiformia) and Cucujoidea, which can no longer be considered a valid taxon. The study shows that, although compositional heterogeneity is not universal, it cannot be eliminated for some mitochondrial genes, but dense taxon sampling and the use of appropriate Bayesian analyses can still produce robust phylogenetic trees.

The sperm ultrastructure of Stictoleptura cordigera (Füssli, 1775) (Insecta, Coleoptera, Cerambycidae)

The spermatozoa of the longhorn beetles Stictoleptura cordigera were ultrastructurally described in this paper. They have an apical bilayered acrosome, an elongated nucleus, a centriole with star-shape links, two asymmetric mitochondrial derivatives partially crystallized and a 9+9+2 flagellar axoneme with accessory tubules provided with 16 protofilaments in their wall. A centriole adjunct is present and gives rise to two thick laminae as accessory bodies, also asymmetrical, to which two relatively small puff-like structures of different size are connected. These features were previously found in the sperm of the cerambycid Morimus asper. The strict similarity of the cerambycid sperm characters with those of curculionoids indicates a clear phylogenetic relationship between Chrysomeloidea and Curculionoidea.

Horizontal gene transfer and functional diversification of plant cell wall degrading polygalacturonases: Key events in the evolution of herbivory in beetles

Plant cell walls are the largest reservoir of organic carbon on earth. To breach and utilize this carbohydrate-rich protective barrier, microbes secrete plant cell wall degrading enzymes (PCWDEs) targeting pectin, cellulose and hemicelluloses. There is a growing body of evidence that genomes of some herbivorous insects also encode PCWDEs, raising questions about their evolutionary origins and functions. Among herbivorous beetles, pectin-degrading polygalacturonases (PGs) are found in the diverse superfamilies Chrysomeloidea (leaf beetles, long-horn beetles) and Curculionoidea (weevils). Here our aim was to test whether these arose from a common ancestor of beetles or via horizontal gene transfer (HGT), and whether PGs kept their ancestral function in degrading pectin or evolved novel functions. Transcriptome data derived from 10 beetle species were screened for PG-encoding sequences and used for phylogenetic comparisons with their bacterial, fungal and plant counterparts. These analyses revealed a large family of PG-encoding genes of Chrysomeloidea and Curculionoidea sharing a common ancestor, most similar to PG genes of ascomycete fungi. In addition, 50 PGs from beetle digestive systems were heterologously expressed and functionally characterized, showing a set of lineage-specific consecutively pectin-degrading enzymes, as well as conserved but enzymatically inactive PG proteins. The evidence indicates that a PG gene was horizontally transferred ∼200 million years ago from an ascomycete fungus to a common ancestor of Chrysomeloidea and Curculionoidea. This has been followed by independent duplications in these two lineages, as well as independent replacement in two sublineages of Chrysomeloidea by two other subsequent HGTs. This origin, leading to subsequent functional diversification of the PG gene family within its new hosts, was a key event promoting the evolution of herbivory in these beetles.

Identification and characterization of plant cell wall degrading enzymes from three glycoside hydrolase families in the cerambycid beetle Apriona japonica

Xylophagous insects have evolved to thrive in a highly challenging environment. For example, wood-boring beetles from the family Cerambycidae feed exclusively on woody tissues, and to efficiently access the nutrients present in this sub-optimal environment, they have to cope with the lignocellulose barrier. Whereas microbes of the insect's gut flora were hypothesized to be responsible for the degradation of lignin, the beetle itself depends heavily on the secretion of a range of enzymes, known as plant cell wall degrading enzymes (PCWDEs), to efficiently digest both hemicellulose and cellulose networks. Here we sequenced the larval gut transcriptome of the Mulberry longhorn beetle, Apriona japonica (Cerambycidae, Lamiinae), in order to investigate the arsenal of putative PCWDEs secreted by this species. We combined our transcriptome with all available sequencing data derived from other cerambycid beetles in order to analyze and get insight into the evolutionary history of the corresponding gene families. Finally, we heterologously expressed and functionally characterized the A. japonica PCWDEs we identified from the transcriptome. Together with a range of endo-β-1,4-glucanases, we describe here for the first time the presence in a species of Cerambycidae of (i) a xylanase member of the subfamily 2 of glycoside hydrolase family 5 (GH5 subfamily 2), as well as (ii) an exopolygalacturonase from family GH28. Our analyses greatly contribute to a better understanding of the digestion physiology of this important group of insects, many of which are major pests of forestry worldwide.

Early Cretaceous angiosperms and beetle evolution

The Coleoptera (beetles) constitute almost one-fourth of all known life-forms on earth. They are also among the most important pollinators of flowering plants, especially basal angiosperms. Beetle fossils are abundant, almost spanning the entire Early Cretaceous, and thus provide important clues to explore the co-evolutionary processes between beetles and angiosperms. We review the fossil record of some Early Cretaceous polyphagan beetles including Tenebrionoidea, Scarabaeoidea, Curculionoidea, and Chrysomeloidea. Both the fossil record and molecular analyses reveal that these four groups had already diversified during or before the Early Cretaceous, clearly before the initial rise of angiosperms to widespread floristic dominance. These four beetle groups are important pollinators of basal angiosperms today, suggesting that their ecological association with angiosperms probably formed as early as in the Early Cretaceous. With the description of additional well-preserved fossils and improvements in phylogenetic analyses, our knowledge of Mesozoic beetle-angiosperm mutualisms will greatly increase during the near future.

Do diet and taxonomy influence insect gut bacterial communities?

Many insects contain diverse gut microbial communities. While several studies have focused on a single or small group of species, comparative studies of phylogenetically diverse hosts can illuminate general patterns of host-microbiota associations. In this study, we tested the hypotheses that (i) host diet and (ii) host taxonomy structure intestinal bacterial community composition among insects. We used published 16S rRNA gene sequence data for 58 insect species in addition to four beetle species sampled from the Sevilleta National Wildlife Refuge to test these hypotheses. Overall, gut bacterial species richness in these insects was low. Decaying wood xylophagous insects harboured the richest bacterial gut flora (102.8 species level operational taxonomic units (OTUs)/sample ± 71.7, 11.8 ± 5.9 phylogenetic diversity (PD)/sample), while bees and wasps harboured the least rich bacterial communities (11.0 species level OTUs/sample ± 5.4, 2.6 ± 0.8 PD/sample). We found evidence to support our hypotheses that host diet and taxonomy structure insect gut bacterial communities (P < 0.001 for both). However, while host taxonomy was important in hymenopteran and termite gut community structure, diet was an important community structuring factor particularly for insect hosts that ingest lignocellulose-derived substances. Our analysis provides a baseline comparison of insect gut bacterial communities from which to test further hypotheses concerning proximate and ultimate causes of these associations.

Molecular barcoding for central-eastern European Crioceris leaf-beetles (Coleoptera: Chrysomelidae)

Among Crioceris leaf-beetles, the two most widespread species (Crioceris asparagi and C. duodecimpunctata) are serious invasive plant pests, while another two (C. quatuordecimpunctata and C. quinquepunctata) are rare species restricted to steppe-like habitats in Eurasia. The aim of the research was to check the genetic distinctiveness of these four species and develop barcodes for their molecular identification using the mitochondrial Cytochrome Oxidase I (COI) gene and two nuclear markers: Elongation Factor 1-α (EF1-α) and Internal Transcribed Spacer 1 (ITS1). The identification of each species was possible and reliable with the use of COI and ITS1 markers. EF1-α was omitted in analyses due to its high level of heterozygosity (presence of multiple PCR products). C. duodecimpunctata and C. quatuordecimpunctata were shown to be sister taxa, but the similar genetic distances between all of the species indicate that these species originated almost simultaneously from a common ancestor. Identification of two separate clades in populations of C. quatuordecimpunctata suggested that the clades are isolated and can be considered as separate conservation units.