Beetle Genome Research Articles

Functional and evolutionary insights into chemosensation and specialized herbivory from the genome of the red milkweed beetle, Tetraopes tetrophthalmus (Cerambycidae: Lamiinae).

Tetraopes are aposematic longhorn beetles (Cerambycidae) that feed primarily on toxic plants in the genus Asclepias (milkweeds). Studies of Tetraopes and their host plants have revealed compelling evidence for insect-plant coevolution and cospeciation. We sequenced, assembled and annotated the genome of the common red milkweed beetle, Tetraopes tetrophthalmus, and explored gene content and evolution, focusing on annotated genes putatively involved in chemosensation, allelochemical detoxification, and phytophagy. Comparisons were made to the Asian longhorned beetle (Anoplophora glabripennis) genome. The genome assembly comprised 779Mb distributed across 1057 contigs, with an N50 of 2.21Mb and 13,089 putative genes, including 97.3% of expected single-copy orthologs. Manual curation identified 122 putative odorant receptors (OR) and 162 gustatory receptors (GR), the former number similar to A. glabripennis but the latter only 69% of the A. glabripennis suite. We also documented a greater percentage of pseudogenic GRs and ORs compared to A. glabripennis, suggesting an ongoing reduction in chemosensory function, perhaps related to host specialization. We found lower diversity within certain well-studied gene families predicted to encode putative plant cell wall degrading enzymes in the T. tetrophthalmus genome, perhaps also due to host specialization. Exploring genes relevant to stress and allelochemical detoxification revealed evidence of an abundance of ABC-family genes in the T. tetrophthalmus genome, which may be related to sequestering toxic cardiac glycosides. Our studies further illuminate the genomic basis and evolution of chemosensation in longhorn beetles and provide a new vantage point from which to explore the ecology and evolution of specialized plant-feeding in Tetraopes and other phytophagous beetles.

Comparative analyses of the banded alder borer (Rosalia funebris) and Asian longhorned beetle (Anoplophora glabripennis) genomes reveal significant differences in genome architecture and gene content among these and other Cerambycidae.

Rosalia funebris (RFUNE; Cerambycidae), the banded alder borer, is a longhorn beetle whose larvae feed on the wood of various economically and ecologically significant trees in western North America. Adults are short-lived and not known to consume plant material substantially. We sequenced, assembled, and annotated the RFUNE genome using HiFi and RNASeq data. We documented genome architecture and gene content, focusing on genes putatively involved in plant feeding (phytophagy). Comparisons were made to the well-studied genome of the Asian longhorned beetle (AGLAB; Anoplophora glabripennis) and other Cerambycidae. The 814 Mb RFUNE genome assembly was distributed across 42 contigs, with an N50 of 30.18 Mb. Repetitive sequences comprised 60.27% of the genome, and 99.0% of expected single-copy orthologous genes were fully assembled. We identified 12,657 genes, fewer than in the four other species studied, and 46.4% fewer than for Aromia moschata (same subfamily as RFUNE). Of the 7,258 orthogroups shared between RFUNE and AGLAB, 1,461 had more copies in AGLAB and 1,023 had more copies in RFUNE. We identified 240 genes in RFUNE that putatively arose via horizontal transfer events. The RFUNE genome encoded substantially fewer putative plant cell wall degrading enzymes than AGLAB, which may relate to the longer-lived plant-feeding adults of the latter species. The RFUNE genome provides new insights into cerambycid genome architecture and gene content and provides a new vantage point from which to study the evolution and genomic basis of phytophagy in beetles.

The genomes of Hercules beetles reveal putative adaptive loci and distinct demographic histories in pristine North American forests.

Beetles, despite their remarkable biodiversity and a long history of research, remain lacking in reference genomes annotated with structural variations in loci of adaptive significance. We sequenced and assembled high-quality chromosome-level genomes of four Hercules beetles which exhibit divergence in male horn size and shape and body colouration. The four Hercules beetle genomes were assembled to 11 pseudo-chromosomes, where the three genomes assembled using Nanopore data (Dynastes grantii, D. hyllus and D. tityus) were mapped to the genome assembled using PacBio + Hi-C data (D. maya). We demonstrated a striking similarity in genome structure among the four species. This conservative genome structure may be attributed to our use of the D. maya assembly as the reference; however, it is worth noting that such a conservative genome structure is a recurring phenomenon among scarab beetles. We further identified homologues of nine and three candidate-gene families that may be associated with the evolution of horn structure and body colouration respectively. Structural variations in Scr and Ebony2 were detected and discussed for their putative impacts on generating morphological diversity in beetles. We also reconstructed the demographic histories of the four Hercules beetles using heterozygosity information from the diploid genomes. We found that the demographic histories of the beetles closely recapitulated historical changes in suitable forest habitats driven by climate shifts.

First evidence of a horizontally-acquired GH-7 cellobiohydrolase from a longhorned beetle genome.

Xylophagous larvae of longhorned beetles (Coleoptera; Cerambycidae) efficiently break down polysaccharides of the plant cell wall, which make the bulk of their food, using a range of carbohydrate-active enzymes (CAZymes). In this study, we investigated the function and evolutionary history of the first identified example of insect-encoded members of glycoside hydrolase family 7 (GH7) derived from the Lamiinae Exocentrus adspersus. The genome of this beetle contained two genes encoding GH7 proteins located in tandem and flanked by transposable elements. Phylogenetic analysis revealed that the GH7 sequences of E. adspersus were closely related to those of Ascomycete fungi, suggesting that they were acquired through horizontal gene transfer (HGT) from fungi. However, they were more distantly related to those encoded by genomes of Crustacea and of protist symbionts of termites and cockroaches, supporting that the same enzyme family was recruited several times independently in Metazoa during the course of their evolution. The recombinant E. adspersus GH7 was found to primarily break down cellulose polysaccharides into cellobiose, indicating that it is a cellobiohydrolase, and could also use smaller cellulose oligomers as substrates. Additionally, the cellobiohydrolase activity was boosted by the presence of calcium chloride. Our findings suggest that the combination of GH7 cellobiohydrolases with other previously characterized endo-β-1,4-glucanases and β-glucosidases allows longhorned beetles like E. adspersus to efficiently break down cellulose into monomeric glucose.

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.

Phylogenetic Analysis of Small Hive Beetles From Native to Introduced Populations.

The small hive beetle (SHB), a social parasite of beehives, is native to sub-Saharan Africa and has spread to America, Europe, and Australia. Recently, these beetles invaded China, causing widespread colony collapses in the honeybee, Apis cerana. In this study, single nucleotide polymorphisms (SNPs) were identified in the beetle genome from its native range (Africa), a region that was invaded by SHBs nearly 30 years ago (America), and more recent invasions (Asia). The beetles in the United States formed the earliest branch and show signs of two decades of gene flow and local adaptation to differentiate this population from the native ones. The beetles in China were deep branched and showed the highest fixation index when compared to the US populations. The number of SNPs in overexpressed genes was significantly higher than the transcriptome. Gene-expression profiles presented here distinguish the characters between adult and larvae SHBs.

The Wnt gene family in Tenebrio molitor and other coleopterans.

The Wnt gene family is involved in a wide range of developmental processes. Despite its significance, the evolution and function of Wnt genes remain largely unclear. Here, an exhaustive survey of Wnt genes was conducted in Tenebrio molitor and 17 other beetle genomes. A total of 146 Wnt genes were identified, creating a comprehensive coleopteran Wnt gene catalog. Comparative genomics indicates that dynamic evolutionary patterns of Wnt gene loss and duplication occurred in Coleoptera, leading to the diverse Wnt gene repertoire in various beetles. A striking loss of particular Wnt gene subfamilies occurs in Coleoptera. Remarkably, Wnt gene duplication was discovered for the first time in insects. Further analysis of Wnt gene expression in T. molitor indicates that each Wnt gene, including the duplicated ones, has a unique spatial or temporal expression pattern. The current study provides valuable insight into the evolution and functional validation of Wnt genes in Coleoptera.

Firefly genomes illuminate the evolution of beetle bioluminescent systems.

Fireflies are one of the best-known bioluminescent organisms, and the reaction mechanism and ecological utility of bioluminescence have been well-studied. Genome assemblies of six species of bioluminescent beetles have recently been published. These studies have focused on the evolution of novelties; luciferase, and the biosynthesis of luciferin and defensive chemicals. For example, clustering of the luciferase gene with acyl-CoA synthetase genes on a chromosome in luminous beetle genomes suggests the involvement of tandem gene duplications and neofunctionalization during the evolution of beetle bioluminescence. Several candidate genes for critical roles in beetle bioluminescence have been identified, but their functional analyses are still ongoing. The establishment of a long-term mass-rearing system and strain will be the key for the post-genome research on bioluminescent beetles. Lastly, the application of contemporary chromosome-scale genome assembly techniques to luminous beetles will help resolve outstanding evolutionary questions, such as how many times bioluminescence evolved in this clade.

A high-quality carabid genome assembly provides insights into beetle genome evolution and cold adaptation.

The hyperdiverse order Coleoptera comprises a staggering ~25% of known species on Earth. Despite recent breakthroughs in next generation sequencing, there remains a limited representation of beetle diversity in assembled genomes. Most notably, the ground beetle family Carabidae, comprising more than 40,000 described species, has not been studied in a comparative genomics framework using whole genome data. Here we generate a high-quality genome assembly for Nebria riversi, to examine sources of novelty in the genome evolution of beetles, as well as genetic changes associated with specialization to high-elevation alpine habitats. In particular, this genome resource provides a foundation for expanding comparative molecular research into mechanisms of insect cold adaptation. Comparison to other beetles shows a strong signature of genome compaction, with N.riversi possessing a relatively small genome (~147Mb) compared to other beetles, with associated reductions in repeat element content and intron length. Small genome size is not, however, associated with fewer protein-coding genes, and an analysis of gene family diversity shows significant expansions of genes associated with cellular membranes and membrane transport, as well as protein phosphorylation and muscle filament structure. Finally, our genomic analyses show that these high-elevation beetles have endosymbiotic Spiroplasma, with several metabolic pathways (e.g., propanoate biosynthesis) that might complement N.riversi, although its role as a beneficial symbiont or as a reproductive parasite remains equivocal.

Draft genome of a biparental beetle species, Lethrus apterus

BackgroundThe lack of an understanding about the genomic architecture underpinning parental behaviour in subsocial insects displaying simple parental behaviours prevents the development of a full understanding about the evolutionary origin of sociality. Lethrus apterus is one of the few insect species that has biparental care. Division of labour can be observed between parents during the reproductive period in order to provide food and protection for their offspring.ResultsHere, we report the draft genome of L. apterus, the first genome in the family Geotrupidae. The final assembly consisted of 286.93 Mbp in 66,933 scaffolds. Completeness analysis found the assembly contained 93.5% of the Endopterygota core BUSCO gene set. Ab initio gene prediction resulted in 25,385 coding genes, whereas homology-based analyses predicted 22,551 protein coding genes. After merging, 20,734 were found during functional annotation. Compared to other publicly available beetle genomes, 23,528 genes among the predicted genes were assigned to orthogroups of which 1664 were in species-specific groups. Additionally, reproduction related genes were found among the predicted genes based on which a reduction in the number of odorant- and pheromone-binding proteins was detected.ConclusionsThese genes can be used in further comparative and functional genomic researches which can advance our understanding of the genetic basis and hence the evolution of parental behaviour.

Distinct Regulation of the Expression of Satellite DNAs in the Beetle Tribolium castaneum.

In the flour beetle, Tribolium castaneum (peri)centromeric heterochromatin is mainly composed of a major satellite DNA TCAST1 interspersed with minor satellites. With the exception of heterochromatin, clustered satellite repeats are found dispersed within euchromatin. In order to uncover a possible satellite DNA function within the beetle genome, we analysed the expression of the major TCAST1 and a minor TCAST2 satellite during the development and upon heat stress. The results reveal that TCAST1 transcription was strongly induced at specific embryonic stages and upon heat stress, while TCAST2 transcription is stable during both processes. TCAST1 transcripts are processed preferentially into piRNAs during embryogenesis and into siRNAs during later development, contrary to TCAST2 transcripts, which are processed exclusively into piRNAs. In addition, increased TCAST1 expression upon heat stress is accompanied by the enrichment of the silent histone mark H3K9me3 on the major satellite, while the H3K9me3 level at TCAST2 remains unchanged. The transcription of the two satellites is proposed to be affected by the chromatin state: heterochromatin and euchromatin, which are assumed to be the prevalent sources of TCAST1 and TCAST2 transcripts, respectively. In addition, distinct regulation of the expression might be related to diverse roles that major and minor satellite RNAs play during the development and stress response.

Genomic content of chemosensory receptors in two sister blister beetles facilitates characterization of chemosensory evolution

BackgroundMore than 2500 species belong to the Meloidae family (Coleoptera: Tenebrionoidea), members of which produce the potent defensive blistering agent cantharidin and are commonly known as blister beetles or Spanishflies. Cantharidin has recently been used for cancer therapy. Hycleus cichorii and Hycleus phaleratus have been used in traditional Chinese medicine for more than 2000 years due to their ability to biosynthesize cantharidin. To understand the role of the chemosensory system in beetle evolution, we comparatively analysed the chemosensory receptor families of both blister beetle species and compared them with those of other beetles.ResultsWe identified 89 odorant receptors (ORs), 86 gustatory receptors (GRs), and 45 ionotropic receptors (IRs) in H. phaleratus and 149 ORs, 102 GRs and 50 IRs in H. cichorii. Nine groups of beetle ORs were recovered, and a similar pattern of ORs in Coleoptera emerged. Two evident expanded clades in Hycleus (Groups 5A and 3) were reconstructed in the phylogenetic tree. Four of eight genes with evidence of positive selection were clustered in the expanded clades of Group 5A. Three, eight and three orthologous pairs of CO2, sugar and fructose receptors, respectively, were identified in both blister beetles. Two evident expanded clades of putative bitter GRs in Hycleus were also found, and the GR in one clade had notably low divergence. Interestingly, IR41a was specifically expanded in blister beetles compared to other insects identified to date, and IR75 was also clearly expanded in both blister beetles based on our phylogenetic tree analysis. Moreover, evidence of positive selection was detected for eight ORs, three GRs and two IRs, half of which were from five duplicate clades.ConclusionsWe first annotated the chemosensory receptor families in a pair of sister beetle genomes (Meloidae: Hycleus), which facilitated evolutionary analysis of chemosensory receptors between sibling species in the Coleoptera group. Our analysis suggests that changes in chemosensory receptors have a possible role in chemical-based species evolution in blister beetles. Future studies should include more species to verify this correlation, which will help us understand the evolution of blister beetles.

Comparison of microsatellite distribution patterns in twenty-nine beetle genomes.

Simple sequence repeats (SSRs) represent an important source of genetic variation that provides a basis for adaptation to different environments in organisms. In this study, we examined the distribution patterns of SSRs in twenty-nine beetle genomes and carried out Gene Ontology (GO) analysis of CDSs embedded with perfect SSRs (P-SSRs). The results demonstrated that imperfect SSRs (I-SSRs) represented the most abundant SSR category in beetle genomes and in different genomic regions (CDS, exon, and intron regions). The numbers of P-SSRs, I-SSRs, compound SSRs, and variable number tandem repeats were positively correlated with beetle genome size, whereas neither the frequency nor the density of the SSRs was correlated with genome size. Moreover, our results demonstrated that common genomic features of P-SSRs within the same suborder or family of Coleoptera were rare. Mono-, di-, tri-, or tetranucleotide SSRs were the most abundant P-SSR categories in beetle genomes. The preferred predominant repeat motif among the mononucleotide P-SSRs was (A)n, but the most frequent repeat motifs for other length classes varied differentially among these genomes. Furthermore, the P-SSR type with the highest GC content differed in the beetle genomes and in different genomic regions. CV (coefficient of variability) analysis demonstrated that the repeat copy numbers of P-SSRs presented relatively higher variation in introns than in CDSs and exons. The GO terms of CDSs containing P-SSRs for molecular functions were mainly enriched in "binding" and "transcription". Our findings will be useful for studying the functional roles of microsatellite heterogeneity in beetle adaptation.

In silico identification and assessment of insecticide target sites in the genome of the small hive beetle, Aethina tumida

BackgroundThe small hive beetle, Aethina tumida, is a rapidly emerging global pest of honey bee colonies. Small hive beetle infestation can be extremely destructive, which may cause honey bees to abscond and render colony infrastructure unusable. Due to the impacts small hive beetles have on honey bees, a wide variety of physical, cultural, and chemical control measures have been implemented to manage small hive beetle infestations. The use of insecticides to control small hive beetle populations is an emerging management tactic. Currently, very little genomic information exists on insecticide target sites in the small hive beetle. Therefore, the objective of this study is to utilize focused in silico comparative genomics approaches to identify and assess the potential insecticide sensitivity of the major insecticide target sites in the small hive beetle genome.ResultsNo previously described resistance mutations were identified in any orthologs of insecticide target sites. Alternative exon use and A-to-I RNA editing were absent in AtumSC1. The ryanodine receptor in small hive beetle (Atum_Ryr) was highly conserved and no previously described resistance mutations were identified. A total of 12 nAChR subunits were identified with similar alternative exon use in other insects. Alternative exon use and critical structural features of the GABA-gated chloride channel subunits (Atum_RDL, Atum_GRD, and Atum_LCCH3) were conserved. Five splice variants were found for the glutamate-gated chloride channel subunit. Exon 3c of Atum_GluCl may be a beetle-specific alternative exon. The co-occurrence of exons 9a and 9b in the pH-sensitive chloride channel (Atum_pHCl) is a unique combination that introduces sites of post-translational modification. The repertoire and alternative exon use for histamine-gated chloride channels (Atum-HisCl), octopamine (Atum_OctR) and tyramine receptors (Atum_TAR) were conserved.ConclusionsThe recently published small hive beetle genome likely serves as a reference for insecticide-susceptible versions of insecticide target sites. These comparative in silico studies are the first step in discovering targets that can be exploited for small hive beetle-specific control as well as tracking changes in the frequency of resistance alleles as part of a resistance monitoring program. Comparative toxicity alongside honey bees is required to verify these in silico predictions.

Draft Genome of a Blister Beetle Mylabris aulica.

Mylabris aulica is a widely distributed blister beetle of the Meloidae family. It has the ability to synthesize a potent defensive secretion that includes cantharidin, a toxic compound used to treat many major illnesses. However, owing to the lack of genetic studies on cantharidin biosynthesis in M. aulica, the commercial use of this species is less extensive than that of other blister beetle species in China. This study reports a draft assembly and possible genes and pathways related to cantharidin biosynthesis for the M. aulica blister beetle using nanopore sequencing data. The draft genome assembly size was 288.5 Mb with a 467.8 Kb N50, and a repeat content of 50.62%. An integrated gene finding pipeline performed for assembly obtained 16,500 protein coding genes. Benchmarking universal single-copy orthologs assessment showed that this gene set included 94.4% complete Insecta universal single-copy orthologs. Over 99% of these genes were assigned functional annotations in the gene ontology, Kyoto Encyclopedia of Genes and Genomes, or Genbank non-redundant databases. Comparative genomic analysis showed that the completeness and continuity of our assembly was better than those of Hycleus cichorii and Hycleus phaleratus blister beetle genomes. The analysis of homologous orthologous genes and inference from evolutionary history imply that the Mylabris and Hycleus genera are genetically close, have a similar genetic background, and have differentiated within one million years. This M. aulica genome assembly provides a valuable resource for future blister beetle studies and will contribute to cantharidin biosynthesis.

Draft genomes of two blister beetles Hycleus cichorii and Hycleus phaleratus.

BackgroundCommonly known as blister beetles or Spanish fly, there are more than 1500 species in the Meloidae family (Hexapoda: Coleoptera: Tenebrionoidea) that produce the potent defensive blistering agent cantharidin. Cantharidin and its derivatives have been used to treat cancers such as liver, stomach, lung, and esophageal cancers. Hycleus cichorii and Hycleus phaleratus are the most commercially important blister beetles in China due to their ability to biosynthesize this potent vesicant. However, there is a lack of genome reference, which has hindered development of studies on the biosynthesis of cantharidin and a better understanding of its biology and pharmacology.ResultsWe report 2 draft genomes and quantified gene sets for the blister beetles H. cichorii and H. phaleratus, 2 complex genomes with >72% repeats and approximately 1% heterozygosity, using Illumina sequencing data. An integrated assembly pipeline was performed for assembly, and most of the coding regions were obtained. Benchmarking universal single-copy orthologs (BUSCO) assessment showed that our assembly obtained more than 98% of the Endopterygota universal single-copy orthologs. Comparison analysis showed that the completeness of coding genes in our assembly was comparable to other beetle genomes such as Dendroctonus ponderosae and Agrilus planipennis. Gene annotation yielded 13 813 and 13 725 protein-coding genes in H. cichorii and H. phaleratus, of which approximately 89% were functionally annotated. BUSCO assessment showed that approximately 86% and 84% of the Endopterygota universal single-copy orthologs were annotated completely in these 2 gene sets, whose completeness is comparable to that of D. ponderosae and A. planipennis.ConclusionsAssembly of both blister beetle genomes provides a valuable resource for future biosynthesis of cantharidin and comparative genomic studies of blister beetles and other beetles.

Genome of the Asian longhorned beetle (Anoplophora glabripennis), a globally significant invasive species, reveals key functional and evolutionary innovations at the beetle-plant interface.

BackgroundRelatively little is known about the genomic basis and evolution of wood-feeding in beetles. We undertook genome sequencing and annotation, gene expression assays, studies of plant cell wall degrading enzymes, and other functional and comparative studies of the Asian longhorned beetle, Anoplophora glabripennis, a globally significant invasive species capable of inflicting severe feeding damage on many important tree species. Complementary studies of genes encoding enzymes involved in digestion of woody plant tissues or detoxification of plant allelochemicals were undertaken with the genomes of 14 additional insects, including the newly sequenced emerald ash borer and bull-headed dung beetle.ResultsThe Asian longhorned beetle genome encodes a uniquely diverse arsenal of enzymes that can degrade the main polysaccharide networks in plant cell walls, detoxify plant allelochemicals, and otherwise facilitate feeding on woody plants. It has the metabolic plasticity needed to feed on diverse plant species, contributing to its highly invasive nature. Large expansions of chemosensory genes involved in the reception of pheromones and plant kairomones are consistent with the complexity of chemical cues it uses to find host plants and mates.ConclusionsAmplification and functional divergence of genes associated with specialized feeding on plants, including genes originally obtained via horizontal gene transfer from fungi and bacteria, contributed to the addition, expansion, and enhancement of the metabolic repertoire of the Asian longhorned beetle, certain other phytophagous beetles, and to a lesser degree, other phytophagous insects. Our results thus begin to establish a genomic basis for the evolutionary success of beetles on plants.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-016-1088-8) contains supplementary material, which is available to authorized users.

Unexpected functional diversity in the fatty acid desaturases of the flour beetle Tribolium castaneum and identification of key residues determining activity

Desaturases catalyse modifications to fatty acids which are essential to homeostasis and for pheromone and defensive chemical production. All desaturases of the flour beetle Tribolium castaneum were investigated via query of the sequenced genome which yielded 15 putative acyl-Coenzyme A genes. Eleven desaturase mRNA were obtained in full length and functionally expressed in yeast. Phylogenetic analysis separated the desaturases into 4 distinct clades; one clade contained conserved beetle Δ9 desaturases, second clade was Tribolium-specific having diverse activities including Δ5, Δ9 and Δ12 desaturation and the other 2 clades had mixed insect representatives. Three members of this clade contained unusual inserted sequences of ∼20 residues in the C-terminal region and were related to desaturases that all contained similar inserts. Deletion of the entirety of the insert in the flour beetle Δ12 desaturase abolished its activity but this was partially restored by the reintroduction of two histidine residues, suggesting the histidine(s) are required for activity but the full length insert is not. Five new desaturase activities were discovered: Δ9 desaturation of C12:0–C16:0 substrates; two unprecedented Δ5 enzymes acting on C18:0 and C16:0; Δ9 activity exclusively on C16:0 and a further stearate Δ9 desaturase. qPCR analysis ruled out a role in sex pheromone synthesis for the Δ5 and Δ9/C16:0 desaturases. The flour beetle genome has underpinned an examination of all transcribed desaturases in the organism and revealed a diversity of novel and unusual activities, an improved understanding of the evolutionary relationships among insect desaturases and sequence determinants of activity.

Studying the organization of genes encoding plant cell wall degrading enzymes inChrysomela tremulaprovides insights into a leaf beetle genome

The ability of herbivorous beetles from the superfamilies Chrysomeloidea and Curculionoidea to degrade plant cell wall polysaccharides has only recently begun to be appreciated. The presence of plant cell wall degrading enzymes (PCWDEs) in the beetle's digestive tract makes this degradation possible. Sequences encoding these beetle-derived PCWDEs were originally identified from transcriptomes and strikingly resemble those of saprophytic and phytopathogenic microorganisms, raising questions about their origin; e.g. are they insect- or microorganism-derived? To demonstrate unambiguously that the genes encoding PCWDEs found in beetle transcriptomes are indeed of insect origin, we generated a bacterial artificial chromosome library from the genome of the leaf beetle Chrysomela tremula, containing 18 432 clones with an average size of 143 kb. After hybridizing this library with probes derived from 12 C. tremula PCWDE-encoding genes and sequencing the positive clones, we demonstrated that the latter genes are encoded by the insect's genome and are surrounded by genes possessing orthologues in the genome of Tribolium castaneum as well as in three other beetle genomes. Our analyses showed that although the level of overall synteny between C. tremula and T. castaneum seems high, the degree of microsynteny between both species is relatively low, in contrast to the more closely related Colorado potato beetle.

A remarkably stable TipE gene cluster: evolution of insect Para sodium channel auxiliary subunits.

BackgroundFirst identified in fruit flies with temperature-sensitive paralysis phenotypes, the Drosophila melanogaster TipE locus encodes four voltage-gated sodium (NaV) channel auxiliary subunits. This cluster of TipE-like genes on chromosome 3L, and a fifth family member on chromosome 3R, are important for the optional expression and functionality of the Para NaV channel but appear quite distinct from auxiliary subunits in vertebrates. Here, we exploited available arthropod genomic resources to trace the origin of TipE-like genes by mapping their evolutionary histories and examining their genomic architectures.ResultsWe identified a remarkably conserved synteny block of TipE-like orthologues with well-maintained local gene arrangements from 21 insect species. Homologues in the water flea, Daphnia pulex, suggest an ancestral pancrustacean repertoire of four TipE-like genes; a subsequent gene duplication may have generated functional redundancy allowing gene losses in the silk moth and mosquitoes. Intronic nesting of the insect TipE gene cluster probably occurred following the divergence from crustaceans, but in the flour beetle and silk moth genomes the clusters apparently escaped from nesting. Across Pancrustacea, TipE gene family members have experienced intronic nesting, escape from nesting, retrotransposition, translocation, and gene loss events while generally maintaining their local gene neighbourhoods. D. melanogaster TipE-like genes exhibit coordinated spatial and temporal regulation of expression distinct from their host gene but well-correlated with their regulatory target, the Para NaV channel, suggesting that functional constraints may preserve the TipE gene cluster. We identified homology between TipE-like NaV channel regulators and vertebrate Slo-beta auxiliary subunits of big-conductance calcium-activated potassium (BKCa) channels, which suggests that ion channel regulatory partners have evolved distinct lineage-specific characteristics.ConclusionsTipE-like genes form a remarkably conserved genomic cluster across all examined insect genomes. This study reveals likely structural and functional constraints on the genomic evolution of insect TipE gene family members maintained in synteny over hundreds of millions of years of evolution. The likely common origin of these NaV channel regulators with BKCa auxiliary subunits highlights the evolutionary plasticity of ion channel regulatory mechanisms.