Transcriptome and full-length cDNA resources for the mountain pine beetle, Dendroctonus ponderosae Hopkins, a major insect pest of pine forests

Efficient Biosynthesis of (+)-α-Pinene and de Novo Synthesis of (+)-cis-Verbenol in Escherichia coli

Bark beetles, major pests that bore into forest stems, cause significant economic damage to forests globally. (+)-α-Pinene is the precursor to (+)-cis-verbenol, a crucial component of the aggregation pheromones produced by bark beetles. This paper describes the de novo synthesis of (+)-cis-verbenol in Escherichia coli. Initially, the truncation position of (+)-α-pinene synthase (PtPS30 from Pinus taeda) and monoterpene precursor (geranyl diphosphate/neryl diphosphate) synthases were evaluated. Neryl diphosphate synthase from Solanum lycopersicum (SlNPPS1) and truncated (+)-α-pinene synthase (PtPS30-39) were selected as promising candidates. Subsequently, the titer of (+)-α-pinene was significantly increased 8.9-fold by using the fusion tag CM29, which enhanced the solubility of PtPS30-39. In addition, by optimizing expression elements (ribosomal binding sites, linkers, and up elements) and overexpressing CM29*PtPS30-39, a yield of 134.12 mg/L (+)-α-pinene was achieved. Finally, the first de novo synthesis of enantiopure (+)-cis-verbenol was achieved by introducing a cytochrome P450 mutant from Pseudomonas putida (P450camF89W,Y98F,L246A), resulting in a yield of 11.13 mg/L. This study lays the groundwork for developing verbenol-based trapping technology for controlling bark beetles.

Identification and Functional Analysis of an Epsilon Class Glutathione S-Transferase Gene Associated with α-Pinene Adaptation in Monochamus alternatus

Alpha-pinene is one of the main defensive components in conifers. Monochamus alternatus (Coleoptera: Cerambycidae), a wood borer feeding on Pinaceae plants, relies on its detoxifying enzymes to resist the defensive terpenoids. Here, we assayed the peroxide level and GST activity of M. alternatus larvae treated with different concentrations of α-pinene. Meanwhile, a gst gene (MaGSTe3) was isolated and analyzed. We determined its expression level and verified its function. The results showed that α-pinene treatment led to membrane lipid peroxidation and thus increased the GST activity. Expression of MaGSTe3 was significantly upregulated in guts following exposure to α-pinene, which has a similar pattern with the malonaldehyde level. In vitro expression and disk diffusion assay showed that the MaGSTe3 protein had high antioxidant capacity. However, RNAi treatment of MaGSTe3 did not reduce the hydrogen peroxide and malonaldehyde levels, while GST activity was significantly reduced. These results suggested MaGSTe3 takes part in α-pinene adaptation, but it does not play a great role in the resistance of M. alternatus larvae to α-pinene.

Identification of metabolizing enzyme genes associated with xenobiotics and odorants in the predatory stink bug Arma custos based on transcriptome analysis

The predatory stink bug, Arma custos, is a highly effective beneficial predator of crop pests. The lack of gene information related to xenobiotic detoxification and odorant degrading enzymes in the predator stink bugs to date has limited our ability for more in-depth studies of biological control. Hence, we conducted de novo assembly of the A. custos transcriptome from guts, antennae, and other tiussue samples of 5th instar larvae using Illumina sequencing technology. A total of 91, 50 and 23 genes of cytochrome P450 monooxygenases (CYPs), carboxyl/choline esterases (CCEs) and glutathione S-transferases (GSTs) genes were identified, respectively. Gene expansions of CYP3 and CYP4 clans and the hormone and pheromone processing CCE class were found in A. custos. Analysis of tissue-specific expression patterns showed that 37 CYPs, 14 CCEs and 8 GSTs were enriched in guts, and 6 CYPs, 5 CCEs and 2 GSTs were up-regulated in antennae, suggesting their potential roles on xenobiotics detoxification and ordorant degradation. Gene information data presented here could be useful for a deeper understanding of the ecology, physiology and behavior of this beneficial species and could be helpful to improve their bio-control efficiency.

Escaping drought: The pectin methylesterase inhibitor gene Slpmei27 can significantly change drought resistance in tomato

Drought stress will lead to a decrease in tomato yield and poor flavour, yield and quality, resulting in economic losses in agricultural production. Mining the key genes regulating tomato drought resistance is of great significance to improve the drought resistance of tomato plants. The cell wall can directly participate in the plant drought stress response as one of the main components of the cell wall, and the regulation of pectin content in plant drought resistance is still unclear. Here, the candidate gene Solyc08g006690 (Slpmei27) was obtained by fine mapping based on genome sequencing technology (BSA-seq) of late-maturing stress-resistant tomato mutants found in the field. Slpmei27 is expressed in the cell wall. The transient silencing of Slpmei27 by VIGS significantly improved the drought resistance of tomato. Meanwhile, Slpmei27 silencing could significantly change the cell wall structure of plants, change the stomatal pass rate, reduce the water loss rate of plants, improve the scavenging ability of reactive oxygen species, change the redox balance in plants, and thus improve the drought resistance of tomato. The promoter region of this gene contains a large number of hormone-response and stress-response binding sites. The promoter region of the Slpmei27 gene in the mutant could lower the expression of downstream genes. Through this study, the mechanism by which Slpmei27 improves tomato drought resistance was revealed, and the relationship between pectin methyl ester metabolism and plant drought resistance was established, providing a theoretical basis for the production of high-quality tomato materials with high drought resistance.

A full-length transcriptome and gene expression analysis of three detoxification gene families in a predatory stink bug, Picromerus lewisi

The predatory stink bug P. Lewisi shows potential for Integrated Pest Management programs for controlling Lepidoptera pest insects in crops and forests. The importance of this insect for biological control has stimulated several studies into its biology and ecology. However, P. lewisi has little genetic information available. In the present study, PacBio single-molecule real-time (SMRT) sequencing and Illumina RNA-seq sequencing technologies were used to reveal the full-length transcriptome profiling and tissue-specific expression patterns of P. lewisi. A total of 12,997 high-quality transcripts with an average length of 2,292 bp were obtained from different stages of P. lewisi using SMRT sequencing. Among these, 12,101 were successfully annotated in seven public databases. A total of 67 genes of cytochrome P450 monooxygenases, 43 carboxylesterase genes, and 18 glutathione S-transferase genes were identified, most of which were obtained with full-length ORFs. Then, tissue-specific expression patterns of 5th instar nymphs were analyzed using Illumina sequencing. Several candidate genes related to detoxification of insecticides and other xenobiotics as well as the degradation of odors, were identified in the guts and antennae of P. lewisi. The current study offered in-depth knowledge to understand the biology and ecology of this beneficial predator and related species.

Mountain Pine Beetle Epidemic: An Interplay of Terpenoids in Host Defense and Insect Pheromones

The mountain pine beetle epidemic has highlighted the complex interactions of bark beetles with conifer host defenses. In these interactions, oleoresin terpenoids and volatiles, produced and released by the host tree, can be both harmful and beneficial to the beetle's success in colonizing a tree and completing its life cycle. The insect spends almost its entire life, from egg to adult, within the bark and phloem of a pine host, exposed to large quantities of complex mixtures of oleoresin terpenoids. Conifer oleoresin comprises mostly monoterpenes and diterpene resin acids as well as many different sesquiterpenes. It functions as a major chemical and physical defense system. However, the insect has evolved host colonization behavior and enzymes for terpenoid metabolism and detoxification that allow it to overcome some of the terpenoid defenses and, importantly, to co-opt pine monoterpenes as cues for host search and as a precursor for its own pheromone system. The insect-associated microbiome also plays a role in the metabolism of conifer terpenoids.

The Genome of Rhyzopertha dominica (Fab.) (Coleoptera: Bostrichidae): Adaptation for Success

The lesser grain borer, Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae), is a major global pest of cereal grains. Infestations are difficult to control as larvae feed inside grain kernels, and many populations are resistant to both contact insecticides and fumigants. We sequenced the genome of R. dominica to identify genes responsible for important biological functions and develop more targeted and efficacious management strategies. The genome was assembled from long read sequencing and long-range scaffolding technologies. The genome assembly is 479.1 Mb, close to the predicted genome size of 480.4 Mb by flow cytometry. This assembly is among the most contiguous beetle assemblies published to date, with 139 scaffolds, an N₅₀ of 53.6 Mb, and L₅₀ of 4, indicating chromosome-scale scaffolds. Predicted genes from biologically relevant groups were manually annotated using transcriptome data from adults and different larval tissues to guide annotation. The expansion of carbohydrase and serine peptidase genes suggest that they combine to enable efficient digestion of cereal proteins. A reduction in the copy number of several detoxification gene families relative to other coleopterans may reflect the low selective pressure on these genes in an insect that spends most of its life feeding internally. Chemoreceptor genes contain elevated numbers of pseudogenes for odorant receptors that also may be related to the recent ontogenetic shift of R. dominica to a diet consisting primarily of stored grains. Analysis of repetitive sequences will further define the evolution of bostrichid beetles compared to other species. The data overall contribute significantly to coleopteran genetic research.

Gut transcriptome of two bark beetle species stimulated with the same kairomones reveals molecular differences in detoxification pathways

Dendroctonus bark beetles are the most destructive agents in coniferous forests. These beetles come into contact with the toxic compounds of their host's chemical defenses throughout their life cycle, some of which are also used by the insects as kairomones to select their host trees during the colonization process. However, little is known about the molecular mechanisms by which the insects counteract the toxicity of these compounds. Here, two sibling species of bark beetles, D. valens and D. rhizophagus, were stimulated with vapors of a blend of their main kairomones (α-pinene, β-pinene and 3-carene), in order to compare the transcriptional response of their gut. A total of 48 180 unigenes were identified in D. valens and 43 704 in D. rhizophagus, in response to kairomones blend. The analysis of differential gene expression showed a transcriptional response in D. valens (739 unigenes, 0.58–10.36 Log2FC) related to digestive process and in D. rhizophagus (322 unigenes 0.87–13.08 Log2FC) related to xenobiotics metabolism. The expression profiles of detoxification genes mainly evidenced the up-regulation of COEs and GSTs in D. valens, and the up-regulation of P450s in D. rhizophagus. Results suggest that terpenes metabolism comes accompanied by an integral hormetic response, result of compensatory mechanisms, including the activation of other metabolic pathways, to ensure the supply of energy and the survival of organisms which is specific for each species, according to its life history and ecological strategy.

Chromosome-level genome assembly reveals genomic architecture of northern range expansion in the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae)

Genome sequencing methods and assembly tools have improved dramatically since the 2013 publication of draft genome assemblies for the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae). We conducted proximity ligation library sequencing and scaffolding to improve contiguity, and then used linkage mapping and recent bioinformatic tools for correction and further improvement. The new assemblies have dramatically improved contiguity and gaps compared to the originals: N50 values increased 26- to 36-fold, and the number of gaps were reduced by half. Ninety per cent of the content of the assemblies is now contained in 12 and 11 scaffolds for the female and male assemblies, respectively. Based on linkage mapping information, the 12 largest scaffolds in both assemblies represent all 11 autosomal chromosomes and the neo-X chromosome. These assemblies now have nearly chromosome-sized scaffolds and will be instrumental for studying genomic architecture, chromosome evolution, population genomics, functional genomics, and adaptation in this and other pest insects. We also identified regions in two chromosomes, including the ancestral-X portion of the neo-X chromosome, with elevated differentiation between northern and southern Canadian populations.

Expression Levels of Detoxification Enzyme Genes from Dendroctonus armandi (Coleoptera: Curculionidae) Fed on a Solid Diet Containing Pine Phloem and Terpenoids

Simple Summary

The bark beetle is the most well-known pest in coniferous trees worldwide. These insects only leave the host pine bark when they disperse to locate a new host. Determining how Dendroctonus armandi overcome the trees’ terpene-based defense systems has been the key problem in the study of bark beetles. Therefore, the aim of this study was to discover the molecular mechanism of insect detoxification enzymes’ ability to confer resistance to terpenes. For this purpose, the genes of cytochrome P450s, glutathione S-transferases, and carboxylesterases were studied in beetles given diets containing terpenes. The results suggest that beetles express different genes in response to terpenoids, and the responses of multiple detoxifying enzymes indicate these insects’ adaption to their chemical environment.

Abstract

Bark beetles overcome the toxic terpenoids produced by pine trees by both detoxifying and converting them into a pheromone system. Detoxification enzymes such as cytochrome P450s, glutathione S-transferases, and carboxylesterases are involved in the ability of Dendroctonus armandi to adapt to its chemical environment. Ten genes from these three major classes of detoxification enzymes were selected to study how these enzymes help D. armandi to respond to the host defenses. The expression profile of these detoxification enzyme genes was observed in adult beetles after feeding on different types of diet. Significant differences were observed between two types of seminatural diet containing the phloem of pines, and a purely artificial diet containing five monoterpenes ((−)-α-pinene, (−)-β-pinene, (+)-3-carene, (±)-limonene, and turpentine oil) also caused differential transcript levels in the detoxification enzyme genes. The results suggest that monoterpenes enter the beetles through different routes (i.e., respiratory and digestive systems) and cause the expression of different genes in response, which might be involved in pheromone metabolism. In addition, the xenobiotic metabolism in bark beetles should be considered as a system comprising multiple detoxifying enzymes.

RNA Interference-Based Forest Protection Products (FPPs) Against Wood-Boring Coleopterans: Hope or Hype?

Forest insects are emerging in large extension in response to ongoing climatic changes, penetrating geographic barriers, utilizing novel hosts, and influencing many hectares of conifer forests worldwide. Current management strategies have been unable to keep pace with forest insect population outbreaks, and therefore novel and aggressive management strategies are urgently required to manage forest insects. RNA interference (RNAi), a Noble Prize-winning discovery, is an emerging approach that can be used for forest protection. The RNAi pathway is triggered by dsRNA molecules, which, in turn, silences genes and disrupts protein function, ultimately causing the death of the targeted insect. RNAi is very effective against pest insects; however, its proficiency varies significantly among insect species, tissues, and genes. The coleopteran forest insects are susceptible to RNAi and can be the initial target, but we lack practical means of delivery, particularly in systems with long-lived, endophagous insects such as the Emerald ash borer, Asian longhorn beetles, and bark beetles. The widespread use of RNAi in forest pest management has major challenges, including its efficiency, target gene selection, dsRNA design, lack of reliable dsRNA delivery methods, non-target and off-target effects, and potential resistance development in wood-boring pest populations. This review focuses on recent innovations in RNAi delivery that can be deployed against forest pests, such as cationic liposome-assisted (lipids), nanoparticle-enabled (polymers or peptides), symbiont-mediated (fungi, bacteria, and viruses), and plant-mediated deliveries (trunk injection, root absorption). Our findings guide future risk analysis of dsRNA-based forest protection products (FPPs) and risk assessment frameworks incorporating sequence complementarity-based analysis for off-target predictions. This review also points out barriers to further developing RNAi for forest pest management and suggests future directions of research that will build the future use of RNAi against wood-boring coleopterans.

Transcriptomic Characterization of Odorant Binding Proteins in Cacia cretifera thibetana and Their Association with Different Host Emitted Volatiles

Simple Summary

The odorant binding proteins (OBPs) interact with host chemical compounds to elicit olfactory responses. Transcriptome analysis of six different tissues of male and female Cacia cretifera thibetana was performed to unravel the interaction of OBPs with host compounds. In both sexes, differentially expressed genes were associated with the KEGG pathways such as cutin, suberine and wax biosynthesis, glycerophospholipid metabolism, choline metabolism in cancer, and the chemokine signaling pathway. The expression of 11 out of 31 OBPs were confirmed by quantitative RT-PCR and seven were found to be specifically expressed in antennae. CcreOBP6 and CcreOBP10 showed strong affinity for terpineol and trans-2-hexenal exhibiting their potential role as an attractant or repellent to control C. cretifera thibetana.

Abstract

This study characterized the transcriptome of Cacia cretifera thibetana and explored odorant binding proteins (OBPs) and their interaction with host-specific compounds. A total of 36 samples from six different organs including antennae, head, thorax, abdomen, wings, and legs (12 groups with 3 replicates per group) from both male and female insects were collected for RNA extraction. Transcriptomic analysis revealed a total of 89,897 transcripts as unigenes, with an average length of 1036 bp. Between male and female groups, 31,095 transcripts were identified as differentially expressed genes (DEGs). The KEGG pathway analysis revealed 26 DEGs associated with cutin, suberine, and wax biosynthesis and 70, 48, and 62 were linked to glycerophospholipid metabolism, choline metabolism in cancer, and chemokine signaling pathways, respectively. A total of 31 OBP genes were identified. Among them, the relative expression of 11 OBP genes (OBP6, 10, 12, 14, 17, 20, 22, 26, 28, 30, and 31) was confirmed by quantitative RT-PCR in different tissues. Seven OBP genes including CcreOBP6 and CcreOBP10 revealed antennae-specific expression. Further, we selected two OBPs (CcreOBP6 and CcreOBP10) for functional analysis to evaluate their binding affinity with 20 host odorant compounds. The CcreOBP6 and CcreOBP10 exhibited strong binding affinities with terpineol and trans-2-hexenal revealing their potential as an attractant or repellent for controlling C. cretifera thibetana.

Isolation, Expression Profiling, and Regulation via Host Allelochemicals of 16 Glutathione S-Transferases in the Chinese White Pine Beetle, Dendroctonus armandi

The Chinese white pine beetle (Dendroctonus armandi) is undoubtedly one of the most important pests causing ecological damage in the Qinling Mountains. When bark beetles invade conifers, they must overcome host tree defenses, including primary resistance and induced resistance responses. Moreover, this induced resistance occurs following herbivory by bark beetles. Bark beetles have a corresponding defense mechanism for degrading toxic compounds, and glutathione S-transferases (GSTs) can catalyze the binding of endogenous substances that reduce glutathione (GSH) to various harmful electrophilic substrates, increasing their solubility and facilitating their excretion from cells. In this experiment, we successfully obtained sixteen full-length sequences of D. armandi, which belonged to four GST categories (delta, epsilon, sigma, and theta). The transcript levels of sixteen GSTs in D. armandi were compared at four developmental stages (larvae, pupae, teneral adults, and adults), three different tissues (antennae, gut, and reproductive organs), and under various levels of terpenoid stimuli during feeding on phloem tissue to evaluate the various relevant modes of action. This study aids in the understanding of the interaction between monoterpenes and beetles, and beetles’ detoxification through GSTs.

Expression Profiling of Plant Cell Wall-Degrading Enzyme Genes in Eucryptorrhynchus scrobiculatus Midgut

In China, the wood-boring weevil Eucryptorrhynchus scrobiculatus damages and eventually kills the tree of heaven Ailanthus altissima. To feed and digest the cell wall of A. altissima, E. scrobiculatus requires plant cell wall-degrading enzymes (PCWDEs). In the present study, we used next-generation sequencing to analyze the midgut transcriptome of E. scrobiculatus. Using three midgut transcriptomes, we assembled 21,491 unigenes from 167,714,100 clean reads. We identified 25 putative PCWDEs, including 11 cellulases and 14 pectinases. We constructed phylogenetic trees with a maximum likelihood algorithm to elucidate the relationships between sequences of the PCWDE protein families and speculate the functions of the PCWDE genes in E. scrobiculatus. The expression patterns of 17 enzymes in the midgut transcriptome were analyzed in various tissues by quantitative real-time PCR (RT-qPCR). The relative expression levels of 12 genes in the midgut and two genes in the proboscis were significantly higher than those in the other tissues. The proboscis and midgut are the digestive organs of insects, and the high expression level indirectly indicates that these genes are related to digestion. The present study has enabled us to understand the types and numbers of the PCWDEs of E. scrobiculatus and will be helpful for research regarding other weevils’ PCWDEs in the future.

Survival physiology and sex ratio of the Chinese white pine beetle Dendroctonus armandi (Coleoptera: Scolytinae) during host colonization and overwintering

The Chinese white pine beetle Dendroctonus armandi (Coleoptera: Scolytinae) typically displays bivoltinism at altitudes below 1700 m in the Qinling Mountains, China. The periods of host colonization and larval overwintering are two important phases in the life cycle of bark beetles, as it is during these periods that they have to contend with host plant defences and periods of intense cold, respectively. Although during different seasons, the females and males of Chinese white pine beetles show varying tolerances to host plant terpenoids, the sex ratio and survival physiology condition of the two beetle generations are unknown. We investigated the sex ratio of individuals, and also examined the body mass, energy stores, and detoxication enzymes of males and females in each of the two generations in order to determine the overall population stability of each generation. We identified a female-biased sex ratio among adults in both generations. Furthermore, patterns of body mass, energy stores, and detoxication enzymes were found to differ between the two sexes and two seasons. Compared with the males, the females have a larger body mass and higher amounts of stored lipids, which are assumed to be adaptations designed to overcome host resistance and facilitate subsequent oviposition.

Genomic dissection of an extended phenotype: Oak galling by a cynipid gall wasp

Galls are plant tissues whose development is induced by another organism for the inducer's benefit. 30,000 arthropod species induce galls, and in most cases the inducing effectors and target plant systems are unknown. Cynipid gall wasps are a speciose monophyletic radiation that induce structurally complex galls on oaks and other plants. We used a model system comprising the gall wasp Biorhiza pallida and the oak Quercus robur to characterise inducer and host plant gene expression at defined stages through the development of galled and ungalled plant tissues, and tested alternative hypotheses for the origin and type of galling effectors and plant metabolic pathways involved. Oak gene expression patterns diverged markedly during development of galled and normal buds. Young galls showed elevated expression of oak genes similar to legume root nodule Nod factor-induced early nodulin (ENOD) genes and developmental parallels with oak buds. In contrast, mature galls showed substantially different patterns of gene expression to mature leaves. While most oak transcripts could be functionally annotated, many gall wasp transcripts of interest were novel. We found no evidence in the gall wasp for involvement of third-party symbionts in gall induction, for effector delivery using virus-like-particles, or for gallwasp expression of genes coding for plant hormones. Many differentially and highly expressed genes in young larvae encoded secretory peptides, which we hypothesise are effector proteins exported to plant tissues. Specifically, we propose that host arabinogalactan proteins and gall wasp chitinases interact in young galls to generate a somatic embryogenesis-like process in oak tissues surrounding the gall wasp larvae. Gall wasp larvae also expressed genes encoding multiple plant cell wall degrading enzymes (PCWDEs). These have functional orthologues in other gall inducing cynipids but not in figitid parasitoid sister groups, suggesting that they may be evolutionary innovations associated with cynipid gall induction.

Plant galls are induced by organisms that manipulate host plant development to produce novel structures. The organisms involved range from mutualistic (such as nitrogen fixing bacteria) to parasitic. In the case of parasites, the gall benefits only the gall-inducing partner. A wide range of organisms can induce galls, but the processes involved are understood only for some bacterial and fungal galls. Cynipid gall wasps induce diverse and structurally complex galls, particularly on oaks (Quercus). We used transcriptome and genome sequencing for one gall wasp and its host oak to identify genes active in gall development. On the plant side, when compared to normally developing bud tissues, young gall tissues showed elevated expression of loci similar to those found in nitrogen-fixing root nodules of leguminous plants. On the wasp side, we found no evidence for involvement of viruses or microorganisms carried by the insects in gall induction or delivery of inducing stimuli. We found that gall wasps express many genes whose products may be secreted to the host, including enzymes that degrade plant cell walls. Genome comparisons between galling and non-galling relatives showed cell wall-degrading enzymes are restricted to gall inducers, and hence potentially key components of a gall inducing lifestyle.

Transcriptome and gene expression analysis of three developmental stages of the coffee berry borer, Hypothenemus hampei

Coffee production is a global industry valued at approximately 173 billion US dollars. One of the main challenges facing coffee production is the management of the coffee berry borer (CBB), Hypothenemus hampei, which is considered the primary arthropod pest of coffee worldwide. Current control strategies are inefficient for CBB management. Although biotechnological alternatives, including RNA interference (RNAi), have been proposed in recent years to control insect pests, characterizing the genetics of the target pest is essential for the successful application of these emerging technologies. In this study, we employed RNA-seq to obtain the transcriptome of three developmental stages of the CBB (larva, female and male) to increase our understanding of the CBB life cycle in relation to molecular features. The CBB transcriptome was sequenced using Illumina Hiseq and assembled de novo. Differential gene expression analysis was performed across the developmental stages. The final assembly produced 29,434 unigenes, of which 4,664 transcripts were differentially expressed. Genes linked to crucial physiological functions, such as digestion and detoxification, were determined to be tightly regulated between the reproductive and nonreproductive stages of CBB. The data obtained in this study help to elucidate the critical roles that several genes play as regulatory elements in CBB development.

Genomic content of chemosensory genes correlates with host range in wood-boring beetles (Dendroctonus ponderosae, Agrilus planipennis, and Anoplophora glabripennis)

Background

Olfaction and gustation underlie behaviors that are crucial for insect fitness, such as host and mate selection. The detection of semiochemicals is mediated via proteins from large and rapidly evolving chemosensory gene families; however, the links between a species’ ecology and the diversification of these genes remain poorly understood. Hence, we annotated the chemosensory genes from genomes of select wood-boring coleopterans, and compared the gene repertoires from stenophagous species with those from polyphagous species.

Results

We annotated 86 odorant receptors (ORs), 60 gustatory receptors (GRs), 57 ionotropic receptors (IRs), 4 sensory neuron membrane proteins (SNMPs), 36 odorant binding proteins (OBPs), and 11 chemosensory proteins (CSPs) in the mountain pine beetle (Dendroctonus ponderosae), and 47 ORs, 30 GRs, 31 IRs, 4 SNMPs, 12 OBPs, and 14 CSPs in the emerald ash borer (Agrilus planipennis). Four SNMPs and 17 CSPs were annotated in the polyphagous wood-borer Anoplophora glabripennis. The gene repertoires in the stenophagous D. ponderosae and A. planipennis are reduced compared with those in the polyphagous A. glabripennis and T. castaneum, which is largely manifested through small gene lineage expansions and entire lineage losses. Alternative splicing of GR genes was limited in D. ponderosae and apparently absent in A. planipennis, which also seems to have lost one carbon dioxide receptor (GR1). A. planipennis has two SNMPs, which are related to SNMP3 in T. castaneum. D. ponderosae has two alternatively spliced OBP genes, a novel OBP “tetramer”, and as many as eleven IR75 members. Simple orthology was generally rare in beetles; however, we found one clade with orthologues of putative bitter-taste GRs (named the “GR215 clade”), and conservation of IR60a from Drosophila melanogaster.

Conclusions

Our genome annotations represent important quantitative and qualitative improvements of the original datasets derived from transcriptomes of D. ponderosae and A. planipennis, facilitating evolutionary analysis of chemosensory genes in the Coleoptera where only a few genomes were previously annotated. Our analysis suggests a correlation between chemosensory gene content and host specificity in beetles. Future studies should include additional species to consolidate this correlation, and functionally characterize identified proteins as an important step towards improved control of these pests.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6054-x) contains supplementary material, which is available to authorized users.

Functions of mountain pine beetle cytochromes P450 CYP6DJ1, CYP6BW1 and CYP6BW3 in the oxidation of pine monoterpenes and diterpene resin acids

The mountain pine beetle (MPB; Dendroctonus ponderosae) is a forest insect pest that attacks several different pine (Pinus) species in its native range of distribution in western North America. MPB are exposed for most of their life cycle to the chemical defenses of their hosts. These defenses are dominated by oleoresin secretions containing mostly various monoterpenes and diterpene resin acids (DRAs). Cytochrome P450 enzymes (P450s) of the MPB are thought to be involved in the metabolism of at least some of these defense compounds. Here we describe the cloning and characterization of three MPB P450s, CYP6DJ1, CYP6BW1 and CYP6BW3, and their functions in the oxidation of various monoterpenes and diterpene resin acids. CYP6DJ1 oxidizes the monoterpenes (+)-(4R)-limonene, (-)-(4S)-limonene and terpinolene and produces (4R,8R)-limonene-8,9-epoxide, (4R,8S)-limonene-8,9-epoxide, (4S,8S)-limonene-8,9-epoxide, (4S,8R)-limonene-8,9-epoxide, perilla alcohol and several unidentified oxidized compounds. These products of CYP6DJ1 were also identified in extracts of MPB treated with the same monoterpenes. CYP6BW1 and CYP6BW3 both oxidize the DRAs abietic acid, dehydroabietic acid, neoabietic acid, levopimaric acid, palustric acid, and isopimaric acid, producing hydroxylated and epoxidized DRAs. CYP6DJ1, CYP6BW1 and CYP6BW3 appear to contribute to the metabolism of oleoresin terpenes as part of the MPB's ability to cope with host defenses.

Time-Course of CYP450 Genes Expression From Dendroctonus rhizophagus (Curculionidae: Scolytinae) During Early Hours of Drilling Bark and Settling Into the Host Tree

Dendroctonus bark beetles (Scolytinae) are one of the most important disturbance agents of coniferous forests in North and Central America. These beetles spend their lives almost entirely under the tree bark, and their survival and reproductive success depend on their ability to overcome the toxic effect of the trees' oleoresin. The cytochromes P450 (CYPs) are associated with the detoxification process of xenobiotics, as well as other physiological processes. Different cytochromes (families 4, 6, and 9) in the Dendroctonus species have been expressed under several experimental conditions; nevertheless, the expression time-course of these genes is unknown. To explore the induction speed of CYPs, we evaluated the relative expression of the CYP6BW5, CYP6DG1, CYP6DJ2, CYP9Z18, and CYP9Z20 genes at the early hours of drilling and settling into a tree (1, 2, 4, 6, 8, 12, 18 h) both in females and males, solitary or paired, of the bark beetle Dendroctonus rhizophagus Thomas and Bright. Our findings show that the five genes were rapidly overexpressed in the early hours (1 to 6 h) in both sexes and in solitary and paired conditions, suggesting their participation in the detoxification process. Additionally, the CYPs expression shows up- and down-regulation patterns through these short times, suggesting their probable participation in other physiological processes as the biosynthesis of hormones, pheromones or compounds related to reproduction.

The cytochrome P450 CYP6DE1 catalyzes the conversion of α-pinene into the mountain pine beetle aggregation pheromone trans-verbenol

The recent outbreak of the mountain pine beetle (Dendroctonus ponderosae; MPB) has affected over 20 M hectares of pine forests in western North America. During the colonization of host trees, female MPB release the aggregation pheromone (-)-trans-verbenol. (-)-trans-Verbenol is thought to be produced from the pine defense compound (-)-α-pinene by cytochrome P450 (P450) dependent hydroxylation. MPB may also use P450s for the detoxification of other monoterpenes of the pine defense system. Here we describe the functional characterization of MPB CYP6DE1. CYP6DE1, but not the closely related CYP6DE2, used the bicyclic monoterpenes (-)-α-pinene, (+)-α-pinene, (-)-β-pinene, (+)-β-pinene and (+)-3-carene as substrates. CYP6DE1 was not active with other monoterpenes or diterpene resin acids that were tested as substrates. trans-Verbenol is the major product of CYP6DE1 activity with (-)-α-pinene or (+)-α-pinene as substrates. When tested with blends of different ratios of (-)-α-pinene and (+)-α-pinene, CYP6DE1 produced trans-verbenol with an enantiomeric profile that was similar to that produced by female MPB exposed to the α-pinene enantiomers.

Functional characterization of CYP4G11-a highly conserved enzyme in the western honey bee Apis mellifera

Determining the functionality of CYP4G11, the only CYP4G in the genome of the western honey bee Apis mellifera, can provide insight into its reduced CYP4 inventory. Toward this objective, CYP4G11 transcripts were quantified, and CYP4G11 was expressed as a fusion protein with housefly CPR in Sf9 cells. Transcript levels varied with age, task, and tissue type in a manner consistent with the need for cuticular hydrocarbon production to prevent desiccation or with comb wax production. Young larvae, with minimal need for desiccation protection, expressed CYP4G11 at very low levels. Higher levels were observed in nurses, and even higher levels in wax producers and foragers, the latter of which risk desiccation upon leaving the hive. Recombinant CYP4G11 readily converted octadecanal to n-heptadecane in a time-dependent manner, demonstrating its functions as an oxidative decarbonylase. CYP4G11 expression levels are high in antennae; heterologously expressed CYP4G11 converted tetradecanal to n-tridecane, demonstrating that it metabolizes shorter-chain aldehydes. Together, these findings confirm the involvement of CYP4G11 in cuticular hydrocarbon production and suggest a possible role in clearing pheromonal and phytochemical compounds from antennae. This possible dual functionality of CYP4G11, i.e., cuticular hydrocarbon and comb wax production and antennal odorant clearance, may explain how honey bees function with a reduced CYP4G inventory.

Monoterpenyl esters in juvenile mountain pine beetle and sex-specific release of the aggregation pheromone trans-verbenol

trans-Verbenol is a critical aggregation pheromone of the female mountain pine beetle. It is derived from a compound of the host defense, the monoterpene α-pinene. We found that beetles accumulate monoterpenyl esters during early life stages. These may serve as a previously unknown reservoir for the sex-specific release of aggregation pheromone when female beetles synchronize the mass colonization of a new host tree. The improved knowledge of the mountain pine beetle pheromone system can inform prediction of bark beetle outbreaks and invasion of new habitats.

A recent outbreak of mountain pine beetle (MPB) has spread over more than 25 million hectares of pine forests in western North America, affecting pine species of sensitive boreal and mountain ecosystems. During initial host colonization, female MPB produce and release the aggregation pheromone trans-verbenol to coordinate a mass attack of individual trees. trans-Verbenol is formed by hydroxylation of α-pinene, a monoterpene of the pine oleoresin defense. It is thought that adult females produce and immediately release trans-verbenol when encountering α-pinene on a new host tree. Here, we show that both sexes of MPB accumulate the monoterpenyl esters verbenyl oleate and verbenyl palmitate during their development in the brood tree. Verbenyl oleate and verbenyl palmitate were retained in adult female MPB until the time of emergence from brood trees, but were depleted in males. Adult females released trans-verbenol in response to treatment with juvenile hormone III (JHIII). While both sexes produced verbenyl esters when exposed to α-pinene, only females responded to JHIII with release of trans-verbenol. Accumulation of verbenyl esters at earlier life stages may allow adult females to release the aggregation pheromone trans-verbenol upon landing on a new host tree, independent of access to α-pinene. Formation of verbenyl esters may be part of a general detoxification system to overcome host monoterpene defenses in both sexes, from which a specialized and female-specific system of pheromone biosynthesis and release may have evolved.

Expansion of cytochrome P450 and cathepsin genes in the generalist herbivore brown marmorated stink bug

BACKGROUND

The brown marmorated stink bug (Halyomorpha halys) is an invasive pest in North America which causes severe economic losses on tree fruits, ornamentals, vegetables, and field crops. The H. halys is an extreme generalist and this feeding behaviour may have been a major contributor behind its establishment and successful adaptation in invasive habitats of North America. To develop an understanding into the mechanism of H. halys' generalist herbivory, here we specifically focused on genes putatively facilitating its adaptation on diverse host plants.

RESULTS

We generated over 142 million reads via sequencing eight RNA-Seq libraries, each representing an individual H. halys adult. The de novo assembly contained 79,855 high quality transcripts, totalling 39,600,178 bases. Following a comprehensive transcriptome analysis, H. halys had an expanded suite of cytochrome P450 and cathepsin-L genes compared to other insects. Detailed characterization of P450 genes from the CYP6 family, known for herbivore adaptation on host plants, strongly hinted towards H. halys-specific expansions involving gene duplications. In subsequent RT-PCR experiments, both P450 and cathepsin genes exhibited tissue-specific or distinct expression patterns which supported their principal roles of detoxification and/or digestion in a particular tissue.

CONCLUSIONS

Our analysis into P450 and cathepsin genes in H. halys offers new insights into potential mechanisms for understanding generalist herbivory and adaptation success in invasive habitats. Additionally, the large-scale transcriptomic resource developed here provides highly useful data for gene discovery; functional, population and comparative genomics as well as efforts to assemble and annotate the H. halys genome.

Toxicity of Pine Monoterpenes to Mountain Pine Beetle

The mountain pine beetle (Dendroctonus ponderosae; MPB) is an eruptive bark beetle species affecting pine forests of western North America. MPB are exposed to volatile monoterpenes, which are important host defense chemicals. We assessed the toxicity of the ten most abundant monoterpenes of lodgepole pine (Pinus contorta), a major host in the current MPB epidemic, against adult MPB from two locations in British Columbia, Canada. Monoterpenes were tested as individual volatiles and included (-)-β-phellandrene, (+)-3-carene, myrcene, terpinolene, and both enantiomers of α-pinene, β-pinene and limonene. Dose-mortality experiments identified (-)-limonene as the most toxic (LC50: 32 μL/L), and (-)-α-pinene (LC50: 290 μL/L) and terpinolene (LC50: >500 μL/L) as the least toxic. MPB body weight had a significant positive effect on the ability to survive most monoterpene volatiles, while sex did not have a significant effect with most monoterpenes. This study helps to quantitatively define the effects of individual monoterpenes towards MPB mortality, which is critical when assessing the variable monoterpene chemical defense profiles of its host species.

Seasonal shifts in accumulation of glycerol biosynthetic gene transcripts in mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae), larvae

Winter mortality is a major factor regulating population size of the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae). Glycerol is the major cryoprotectant in this freeze intolerant insect. We report findings from a gene expression study on an overwintering mountain pine beetle population over the course of 35 weeks. mRNA transcript levels suggest glycerol production in the mountain pine beetle occurs through glycogenolytic, gluconeogenic and potentially glyceroneogenic pathways, but not from metabolism of lipids. A two-week lag period between fall glycogen phosphorylase transcript and phosphoenolpyruvate carboxykinase transcript up-regulation suggests that gluconeogenesis serves as a secondary glycerol-production process, subsequent to exhaustion of the primary glycogenolytic source. These results provide a first look at the details of seasonal gene expression related to the production of glycerol in the mountain pine beetle.

Comparative transcriptomics of mountain pine beetle pheromone-biosynthetic tissues and functional analysis of CYP6DE3

Background

The mountain pine beetle (MPB, Dendroctonus ponderosae Hopkins) is a highly destructive pest of pine forests in western North America. During flight to a new host tree and initiation of feeding, mountain pine beetles release aggregation pheromones. The biosynthetic pathways of these pheromones are sex-specific and localized in the midgut and fat body, but the enzymes involved have not all been identified or characterized.

Results

We used a comparative RNA-Seq analysis between fed and unfed male and female MPB midguts and fat bodies to identify candidate genes involved in pheromone biosynthesis. The 13,407 potentially unique transcripts showed clear separation based on feeding state and gender. Gene co-expression network construction and examination using petal identified gene groups that were tightly connected. This, as well as other co-expression and gene ontology analyses, identified all four known pheromone biosynthetic genes, confirmed the tentative identification of four others from a previous study, and suggested nine novel candidates. One cytochrome P450 monooxygenase, CYP6DE3, identified as a possible exo-brevicomin-biosynthetic enzyme in this study, was functionally characterized and likely is involved in resin detoxification rather than pheromone biosynthesis.

Conclusions

Our analysis supported previously characterized pheromone-biosynthetic genes involved in exo-brevicomin and frontalin biosynthesis and identified a number of candidate cytochrome P450 monooxygenases and a putative cyclase for further studies. Functional analyses of CYP6DE3 suggest its role in resin detoxification and underscore the limitation of using high-throughput data to tentatively identify candidate genes. Further functional analyses of candidate genes found in this study should lead to the full characterization of MPB pheromone biosynthetic pathways and the identification of molecular targets for possible pest management strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3696-4) contains supplementary material, which is available to authorized users.

Transcriptomic alterations in Sitophilus zeamais in response to allyl isothiocyanate fumigation

To study the fumigation mechanisms of Allyl isothiocyanate (AITC) a promising biorational alternative to present fumigants (phosphine and methyl bromide), and provide theoretical basis for its further development in the control of stored grain pests, this research presents a transcriptome analysis of Sitophilus zeamais fumigated with AITC at the concentration of LC₅₀ (5.69μg/mL) and control over 8h. 21,869,022 and 23,873,110 clean reads in insects fumigated with AITC and control were gained, respectively. The results of RNA-seq were confirmed by qRT-PCR determination of the expression levels of NADH dehydrogenase subunit 6 and Vacuolar ATP synthase subunit B in the insects fumigated with AITC at different concentrations. After enrichment analysis of differentially expressed genes, 117 over-expressed and 271 down-regulated transcripts were gained. Following GO enrichment, these transcripts were classified into 38 GO subgroups (at level 2), and the majority enriched GO terms were "Binding" "Cell process" and "metabolic". KEGG enrichment analysis showed that the majority enriched pathway were "Folding, sorting and degradation", "Transport and catabolism", "Energy metabolism", and "Carbohydrate metabolism". Connected with previous researches on mechanisms of isothiocyanates, cytoskeleton collapse and mitochondria dysfunction are proposed to be significant lethal mechanisms of AITC.

Deep, Staged Transcriptomic Resources for the Novel Coleopteran Models Atrachya menetriesi and Callosobruchus maculatus

Despite recent efforts to sample broadly across metazoan and insect diversity, current sequence resources in the Coleoptera do not adequately describe the diversity of the clade. Here we present deep, staged transcriptomic data for two coleopteran species, Atrachya menetriesi (Faldermann 1835) and Callosobruchus maculatus (Fabricius 1775). Our sampling covered key stages in ovary and early embryonic development in each species. We utilized this data to build combined assemblies for each species which were then analysed in detail. The combined A. menetriesi assembly consists of 228,096 contigs with an N50 of 1,598 bp, while the combined C. maculatus assembly consists of 128,837 contigs with an N50 of 2,263 bp. For these assemblies, 34.6% and 32.4% of contigs were identified using Blast2GO, and 97% and 98.3% of the BUSCO set of metazoan orthologs were present, respectively. We also carried out manual annotation of developmental signalling pathways and found that nearly all expected genes were present in each transcriptome. Our analyses show that both transcriptomes are of high quality. Lastly, we performed read mapping utilising our timed, stage specific RNA samples to identify differentially expressed contigs. The resources presented here will provide a firm basis for a variety of experimentation, both in developmental biology and in comparative genomic studies.

petal: Co-expression network modelling in R

Background

Networks provide effective models to study complex biological systems, such as gene and protein interaction networks. With the advent of new sequencing technologies, many life scientists are grasping for user-friendly methods and tools to examine biological components at the whole-systems level. Gene co-expression network analysis approaches are frequently used to successfully associate genes with biological processes and demonstrate great potential to gain further insights into the functionality of genes, thus becoming a standard approach in Systems Biology. Here the objective is to construct biologically meaningful and statistically strong co-expression networks, the identification of research dependent subnetworks, and the presentation of self-contained results.

Results

We introduce petal, a novel approach to generate gene co-expression network models based on experimental gene expression measures. petal focuses on statistical, mathematical, and biological characteristics of both, input data and output network models. Often over-looked issues of current co-expression analysis tools include the assumption of data normality, which is seldom the case for hight-throughput expression data obtained from RNA-seq technologies. petal does not assume data normality, making it a statistically appropriate method for RNA-seq data. Also, network models are rarely tested for their known typical architecture: scale-free and small-world. petal explicitly constructs networks based on both these characteristics, thereby generating biologically meaningful models. Furthermore, many network analysis tools require a number of user-defined input variables, these often require tuning and/or an understanding of the underlying algorithm; petal requires no user input other than experimental data. This allows for reproducible results, and simplifies the use of petal. Lastly, this approach is specifically designed for very large high-throughput datasets; this way, petal’s network models represent as much of the entire system as possible to provide a whole-system approach.

Conclusion

petal is a novel tool for generating co-expression network models of whole-genomics experiments. It is implemented in R and available as a library. Its application to several whole-genome experiments has generated novel meaningful results and has lead the way to new testing hypothesizes for further biological investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-016-0298-8) contains supplementary material, which is available to authorized users.

Gene expression analysis of overwintering mountain pine beetle larvae suggests multiple systems involved in overwintering stress, cold hardiness, and preparation for spring development

Cold-induced mortality has historically been a key aspect of mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae), population control, but little is known about the molecular basis for cold tolerance in this insect. We used RNA-seq analysis to monitor gene expression patterns of mountain pine beetle larvae at four time points during their overwintering period—early-autumn, late-autumn, early-spring, and late-spring. Changing transcript profiles over the winter indicates a multipronged physiological response from larvae that is broadly characterized by gene transcripts involved in insect immune responses and detoxification during the autumn. In the spring, although transcripts associated with developmental process are present, there was no particular biological process dominating the transcriptome.

Transcriptome and Gene Expression Analysis of Cylas formicarius (Coleoptera: Brentidae) During Different Development Stages

The sweet potato weevil, Cylas formicarius (F.) (Coleoptera: Brentidae), is an important pest of sweet potato worldwide. However, there is limited knowledge on the molecular mechanisms underlying growth and differentiation of C. formicarius The transcriptomes of the eggs, second instar larvae, third instar larvae (L3), pupae, females, and males of C. formicarius were sequenced using Illumina sequencing technology for obtaining global insights into developing transcriptome characteristics and elucidating the relative functional genes. A total of 54,255,544 high-quality reads were produced, trimmed, and de novo assembled into 115,281 contigs. 61,686 unigenes were obtained, with an average length of 1,009 nt. Among these unigenes, 17,348 were annotated into 59 Gene Ontology (GO) terms and 12,660 were assigned to 25 Cluster of Orthologous Groups classes, whereas 24,796 unigenes were mapped to 258 pathways. Differentially expressed unigenes between various developmental stages of C. formicarius were detected. Higher numbers of differentially expressed genes (DEGs) were recorded in the eggs versus L3 and eggs versus male samples (2,141 and 2,058 unigenes, respectively) than the others. Genes preferentially expressed in each stage were also identified. GO and pathway-based enrichment analysis were used to further investigate the functions of the DEGs. In addition, the expression profiles of ten DEGs were validated by quantitative real-time PCR. The transcriptome profiles presented in this study and these DEGs detected by comparative analysis of different developed stages of C. formicarius will facilitate the understanding of the molecular mechanism of various living process and will contribute to further genome-wide research.

Structure and Properties of a Non-processive, Salt-requiring, and Acidophilic Pectin Methylesterase from Aspergillus niger Provide Insights into the Key Determinants of Processivity Control

Many pectin methylesterases (PMEs) are expressed in plants to modify plant cell-wall pectins for various physiological roles. These pectins are also attacked by PMEs from phytopathogens and phytophagous insects. The de-methylesterification by PMEs of the O6-methyl ester groups of the homogalacturonan component of pectin, exposing galacturonic acids, can occur processively or non-processively, respectively, describing sequential versus single de-methylesterification events occurring before enzyme-substrate dissociation. The high resolution x-ray structures of a PME from Aspergillus niger in deglycosylated and Asn-linked N-acetylglucosamine-stub forms reveal a 10⅔-turn parallel β-helix (similar to but with less extensive loops than bacterial, plant, and insect PMEs). Capillary electrophoresis shows that this PME is non-processive, halophilic, and acidophilic. Molecular dynamics simulations and electrostatic potential calculations reveal very different behavior and properties compared with processive PMEs. Specifically, uncorrelated rotations are observed about the glycosidic bonds of a partially de-methyl-esterified decasaccharide model substrate, in sharp contrast to the correlated rotations of processive PMEs, and the substrate-binding groove is negatively not positively charged.

Cytochrome P450s from the Chinese white pine beetle, Dendroctonus armandi (Curculionidae: Scolytinae): Expression profiles of different stages and responses to host allelochemicals

Bark beetles oxidize the defensive allelochemicals from their host trees to both detoxify and convert these materials into components of their pheromone system. The ability of this insect to adapt to its chemical environment might be explained by the action of major detoxification enzymes such as cytochrome P450s (or CYPs). Sixty-four sequences coding for P450s were identified, and most of the transcripts were found to be expressed in the larvae, pupae and adults of Dendroctonus armandi. To gain information on how these genes help D. armandi overcome the host defense, differential transcript levels of the CYP genes were observed between sexes and within the sexes. Significant differences were observed among developmental stages, in feeding on the phloem of Pinus armandi and in exposure to stimuli ((±)-α-pinene, (S)-(-)-α-pinene, (S)-(-)-β-pinene, (+)-3-carene, (±)-limonene and turpentine oil) for 8 h. We investigated the effect of sex and generations on the survivorship of individual D. armandi that were exposed to host volatiles at concentrations comparable to constitutive and induced levels of defense using fumigant exposure to understand the ability of the beetles to tolerate host defensive chemicals. The differential transcript accumulation patterns of CYP genes of these bark beetle provided insight into the ecological interactions of D. armandi with its host pine.

Endogenous cellulolytic enzyme systems in the longhorn beetle Mesosa myops (Insecta: Coleoptera) studied by transcriptomic analysis

The Cerambycidae (longhorn beetle) is a large family of Coleoptera with xylophagous feeding habits. Cellulose digestion plays an important role in these wood-feeding insects. In this study, transcriptomic technology was used to obtain one glycoside hydrolase family 45 (GH45) cellulase and seven GH5 cellulases from Mesosa myops, a typical longhorn beetle. Analyses of expression dynamics and evolutionary relationships provided a complete description of the cellulolytic system. The expression dynamics related to individual development indicated that endogenous GH45 and GH5 cellulases dominate cellulose digestion in M. myops. Evolutionary analyses suggested that GH45 cellulase gene is a general gene in the Coleoptera Suborder Polyphaga. Evolutionary analyses also indicated that the GH5 cellulase group in Lamiinae longhorn beetles is closely associated with wood feeding. This study demonstrated that there is a complex endogenous cellulolytic system in M. myops that is dominated by cellulases belonging to two glycoside hydrolase families.

Metagenome Skimming of Insect Specimen Pools: Potential for Comparative Genomics

Metagenomic analyses are challenging in metazoans, but high-copy number and repeat regions can be assembled from low-coverage sequencing by "genome skimming," which is applied here as a new way of characterizing metagenomes obtained in an ecological or taxonomic context. Illumina shotgun sequencing on two pools of Coleoptera (beetles) of approximately 200 species each were assembled into tens of thousands of scaffolds. Repeated low-coverage sequencing recovered similar scaffold sets consistently, although approximately 70% of scaffolds could not be identified against existing genome databases. Identifiable scaffolds included mitochondrial DNA, conserved sequences with hits to expressed sequence tag and protein databases, and known repeat elements of high and low complexity, including numerous copies of rRNA and histone genes. Assemblies of histones captured a diversity of gene order and primary sequence in Coleoptera. Scaffolds with similarity to multiple sites in available coleopteran genome sequences for Dendroctonus and Tribolium revealed high specificity of scaffolds to either of these genomes, in particular for high-copy number repeats. Numerous "clusters" of scaffolds mapped to the same genomic site revealed intra- and/or intergenomic variation within a metagenome pool. In addition to effect of taxonomic composition of the metagenomes, the number of mapped scaffolds also revealed structural differences between the two reference genomes, although the significance of this striking finding remains unclear. Finally, apparently exogenous sequences were recovered, including potential food plants, fungal pathogens, and bacterial symbionts. The "metagenome skimming" approach is useful for capturing the genomic diversity of poorly studied, species-rich lineages and opens new prospects in environmental genomics.

Pectinases from Sphenophorus levis Vaurie, 1978 (Coleoptera: Curculionidae): putative accessory digestive enzymes

The cell wall in plants offers protection against invading organisms and is mainly composed of the polysaccharides pectin, cellulose, and hemicellulose, which can be degraded by plant cell wall degrading enzymes (PCWDEs). Such enzymes are often synthesized by free living microorganisms or endosymbionts that live in the gut of some animals, including certain phytophagous insects. Thus, the ability of an insect to degrade the cell wall was once thought to be related to endosymbiont enzyme activity. However, recent studies have revealed that some phytophagous insects are able to synthesize their own PCWDEs by endogenous genes, although questions regarding the origin of these genes remain unclear. This study describes two pectinases from the sugarcane weevil, Sphenophorus levis Vaurie, 1978 (Sl-pectinases), which is considered one of the most serious agricultural pests in Brazil. Two cDNA sequences identified in a cDNA library of the insect larvae coding for a pectin methylesterase (PME) and an endo-polygalacturonase (endo-PG)-denominated Sl-PME and Sl-endoPG, respectively-were isolated and characterized. The quantitative real-time reverse transcriptase polymerase chain reaction expression profile for both Sl-pectinases showed mRNA production mainly in the insect feeding stages and exclusively in midgut tissue of the larvae. This analysis, together Western blotting data, suggests that Sl-pectinases have a digestive role. Phylogenetic analyses indicate that Sl-PME and Sl-endoPG sequences are closely related to bacteria and fungi, respectively. Moreover, the partial genomic sequences of the pectinases were amplified from insect fat body DNA, which was certified to be free of endosymbiotic DNA. The analysis of genomic sequences revealed the existence of two small introns with 53 and 166 bp in Sl-endoPG, which is similar to the common pattern in fungal introns. In contrast, no intron was identified in the Sl-PME genomic sequence, as generally observed in bacteria. These data support the theory of horizontal gene transfer proposed for the origin of insect pectinases, reinforcing the acquisition of PME genes from bacteria and endo-PG genes from fungi.

Transcriptome analysis and systemic RNAi response in the African sweetpotato weevil (Cylas puncticollis, Coleoptera, Brentidae)

The African sweetpotato weevil (SPW) Cylas puncticollis Boheman is one of the most important constraints of sweetpotato production in Sub-Saharan Africa and yet is largely an uncharacterized insect pest. Here, we report on the transcriptome analysis of SPW generated using an Illumina platform. More than 213 million sequencing reads were obtained and assembled into 89,599 contigs. This assembly was followed by a gene ontology annotation. Subsequently, a transcriptome search showed that the necessary RNAi components relevant to the three major RNAi pathways, were found to be expressed in SPW. To address the functionality of the RNAi mechanism in this species, dsRNA was injected into second instar larvae targeting laccase2, a gene which encodes an enzyme involved in the sclerotization of insect exoskeleton. The body of treated insects showed inhibition of sclerotization, leading eventually to death. Quantitative Real Time PCR (qPCR) confirmed this phenotype to be the result of gene silencing. Together, our results provide valuable sequence data on this important insect pest and demonstrate that a functional RNAi pathway with a strong and systemic effect is present in SPW and can further be explored as a new strategy for controlling this important pest.

Proteomics indicators of the rapidly shifting physiology from whole mountain pine beetle, Dendroctonus ponderosae (Coleoptera: Curculionidae), adults during early host colonization

We developed proteome profiles for host colonizing mountain pine beetle adults, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae). Adult insects were fed in pairs on fresh host lodgepole pine, Pinus contorta Dougl. ex Loud, phloem tissue. The proteomes of fed individuals were monitored using iTRAQ and compared to those of starved beetles, revealing 757 and 739 expressed proteins in females and males, respectively, for which quantitative information was obtained. Overall functional category distributions were similar for males and females, with the majority of proteins falling under carbohydrate metabolism (glycolysis, gluconeogenesis, citric acid cycle), structure (cuticle, muscle, cytoskeleton), and protein and amino acid metabolism. Females had 23 proteins with levels that changed significantly with feeding (p<0.05, FDR<0.20), including chaperones and enzymes required for vitellogenesis. In males, levels of 29 proteins changed significantly with feeding (p<0.05, FDR<0.20), including chaperones as well as motor proteins. Only two proteins, both chaperones, exhibited a significant change in both females and males with feeding. Proteins with differential accumulation patterns in females exhibited higher fold changes with feeding than did those in males. This difference may be due to major and rapid physiological changes occurring in females upon finding a host tree during the physiological shift from dispersal to reproduction. The significant accumulation of chaperone proteins, a cytochrome P450, and a glutathione S-transferase, indicate secondary metabolite-induced stress physiology related to chemical detoxification during early host colonization. The females' activation of vitellogenin only after encountering a host indicates deliberate partitioning of resources and a balancing of the needs of dispersal and reproduction.

Two CYP4 genes of the Chinese white pine beetle, Dendroctonus armandi (Curculionidae: Scolytinae), and their transcript levels under different development stages and treatments

Bark beetles oxidize the defensive monoterpenes of their host trees both to detoxify them and to convert them into components of their pheromone system. This oxidation is catalysed by cytochrome P450 (CYP) enzymes and occurs in different stages of the insect. We identified two new CYP4 genes in the Chinese white pine beetle (Dendroctonus armandi), and carried out bioinformatic analysis one the full-length nucleic acid sequences and deduced amino acid sequences. Differential expression of the CYP4 genes was observed between sexes, and within these significant differences amongst development stages, fed on phloem of Pinus armandi and exposed to stimuli((±)- α-pinene, (R)-(+)- α-pinene, (S)-(-)-α-pinene, (S)-(-)-β-pinene and (+)-3-carene) at 8 and 24 h, and their interactions were found upon exposure to host monoterpenes. Increased expression of CYP4 genes suggested that they play a role in the detoxification of monoterpenes released by the host trees. The differential transcript accumulation patterns of these bark beetle CYP4 genes provides insight into the ecological interactions of D. armandi with its host pine.

TransPS: A Transcriptome Post Scaffolding Method for Assembling High Quality Contigs

MOTIVATION

As the development of the high throughput sequencing technologies, transcriptome can be sequenced with a low price and high efficiency. Sequence assembly approaches have been renewed to meet the new requirements from new sequencing technologies. Assembly strategies are important for biologists who need to assemble the transcriptome generated in their experiments. However, some modern de novo assembly strategies generate a large section of redundant contigs due to sequence variations, which greatly affect downstream analysis and experiments. This work proposed TransPS, a post transcriptome scaffolding method to generate high quality transcriptomes.

RESULTS

TransPS shows promising results on the test transcriptome data sets where the redundancy is greatly reduced by at least 50%, while the coverage is improved considerately.

AVAILABILITY

The web server and source code are available at https://bioinformatics.cs.vt.edu/zhanglab/transps/.

How the mountain pine beetle (Dendroctonus ponderosae) breached the Canadian Rocky Mountains

The mountain pine beetle (MPB; Dendroctonus ponderosae Hopkins), a major pine forest pest native to western North America, has extended its range north and eastward during an ongoing outbreak. Determining how the MPB has expanded its range to breach putative barriers, whether physical (nonforested prairie and high elevation of the Rocky Mountains) or climatic (extreme continental climate where temperatures can be below −40 °C), may contribute to our general understanding of range changes as well as management of the current epidemic. Here, we use a panel of 1,536 single nucleotide polymorphisms (SNPs) to assess population genetic structure, connectivity, and signals of selection within this MPB range expansion. Biallelic SNPs in MPB from southwestern Canada revealed higher genetic differentiation and lower genetic connectivity than in the northern part of its range. A total of 208 unique SNPs were identified using different outlier detection tests, of which 32 returned annotations for products with putative functions in cholesterol synthesis, actin filament contraction, and membrane transport. We suggest that MPB has been able to spread beyond its previous range by adjusting its cellular and metabolic functions, with genome scale differentiation enabling populations to better withstand cooler climates and facilitate longer dispersal distances. Our study is the first to assess landscape-wide selective adaptation in an insect. We have shown that interrogation of genomic resources can identify shifts in genetic diversity and putative adaptive signals in this forest pest species.

Molecular evolution of glycoside hydrolase genes in the Western corn rootworm (Diabrotica virgifera virgifera)

Cellulose is an important nutritional resource for a number of insect herbivores. Digestion of cellulose and other polysaccharides in plant-based diets requires several types of enzymes including a number of glycoside hydrolase (GH) families. In a previous study, we showed that a single GH45 gene is present in the midgut tissue of the western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae). However, the presence of multiple enzymes was also suggested by the lack of a significant biological response when the expression of the gene was silenced by RNA interference. In order to clarify the repertoire of cellulose-degrading enzymes and related GH family proteins in D. v. virgifera, we performed next-generation sequencing and assembled transcriptomes from the tissue of three different developmental stages (eggs, neonates, and third instar larvae). Results of this study revealed the presence of seventy-eight genes that potentially encode GH enzymes belonging to eight families (GH45, GH48, GH28, GH16, GH31, GH27, GH5, and GH1). The numbers of GH45 and GH28 genes identified in D. v. virgifera are among the largest in insects where these genes have been identified. Three GH family genes (GH45, GH48, and GH28) are found almost exclusively in two coleopteran superfamilies (Chrysomeloidea and Curculionoidea) among insects, indicating the possibility of their acquisitions by horizontal gene transfer rather than simple vertical transmission from ancestral lineages of insects. Acquisition of GH genes by horizontal gene transfers and subsequent lineage-specific GH gene expansion appear to have played important roles for phytophagous beetles in specializing on particular groups of host plants and in the case of D. v. virgifera, its close association with maize.

Transcriptome immune analysis of the invasive beetle Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae) parasitized by Tetrastichus brontispae Ferrière (Hymenoptera: Eulophidae)

The beetle Octodonta nipae (Maulik) (Coleoptera: Chrysomelidae) is a serious invasive insect pest of palm plants in southern China, and the endoparasitoid Tetrastichus brontispae Ferrière (Hymenoptera: Eulophidae) is a natural enemy of this pest that exhibits great ability in the biocontrol of O. nipae. For successful parasitism, endoparasitoids often introduce or secrete various virulence factors to suppress host immunity. To investigate the effects of parasitization by T. brontispae on the O. nipae immune system, the transcriptome of O. nipae pupae was analyzed with a focus on immune-related genes through Illumina sequencing. De novo assembly generated 49,919 unigenes with a mean length of 598 bp. Of these genes, 27,490 unigenes (55.1% of all unigenes) exhibited clear homology to known genes in the NCBI nr database. Parasitization had significant effects on the transcriptome profile of O. nipae pupae, and most of these differentially expressed genes were down-regulated. Importantly, the expression profiles of immune-related genes were significantly regulated after parasitization. Taken together, these transcriptome sequencing efforts shed valuable light on the host (O. nipae) manipulation mechanisms induced by T. brontispae, which will pave the way for the development of novel immune defense-based management strategies of O. nipae, and provide a springboard for further molecular analyses, particularly of O. nipae invasion.

Processive pectin methylesterases: the role of electrostatic potential, breathing motions and bond cleavage in the rectification of Brownian motions

Pectin methylesterases (PMEs) hydrolyze the methylester groups that are found on the homogalacturonan (HG) chains of pectic polysaccharides in the plant cell wall. Plant and bacterial PMEs are especially interesting as the resulting de-methylesterified (carboxylated) sugar residues are found to be arranged contiguously, indicating a so-called processive nature of these enzymes. Here we report the results of continuum electrostatics calculations performed along the molecular dynamics trajectory of a PME-HG-decasaccharide complex. In particular it was observed that, when the methylester groups of the decasaccharide were arranged in order to mimic the just-formed carboxylate product of de-methylesterification, a net unidirectional sliding of the model decasaccharide was subsequently observed along the enzyme’s binding groove. The changes that occurred in the electrostatic binding energy and protein dynamics during this translocation provide insights into the mechanism by which the enzyme rectifies Brownian motions to achieve processivity. The free energy that drives these molecular motors is thus demonstrated to be incorporated endogenously in the methylesterified groups of the HG chains and is not supplied exogenously.

RNA-sequencing analysis reveals abundant developmental stage-specific and immunity-related genes in the pollen beetle Meligethes aeneus

The pollen beetle (Meligethes aeneus) is a major pest of oilseed rape (Brassica napus) and other cruciferous crops in Europe. Pesticide-resistant pollen beetle populations are emerging, increasing the economic impact of this species. We isolated total RNA from the larval and adult stages, the latter either naïve or immunized by injection with bacteria and yeast. High-throughput RNA sequencing (RNA-Seq) was carried out to establish a comprehensive transcriptome catalogue and to screen for developmental stage-specific and immunity-related transcripts. We assembled the transcriptome de novo by combining sequence tags from all developmental stages and treatments. Gene expression data based on normalized read counts revealed several functional gene categories that were differentially expressed between larvae and adults, particularly genes associated with digestion and detoxification that were induced in larvae, and genes associated with reproduction and environmental signalling that were induced in adults. We also identified many genes associated with microbe recognition, immunity-related signalling and defence effectors, such as antimicrobial peptides (AMPs) and lysozymes. Digital gene expression analysis revealed significant differences in the profile of AMPs expressed in larvae, naïve adults and immune-challenged adults, providing insight into the steady-state differences between developmental stages and the complex transcriptional remodelling that occurs following the induction of immunity. Our data provide insight into the adaptive mechanisms used by phytophagous insects and could lead to the development of more effective control strategies for insect pests.

Expression of the mevalonate pathway enzymes in the Lutzomyia longipalpis (Diptera: Psychodidae) sex pheromone gland demonstrated by an integrated proteomic approach

In Latin America, Lutzomyia longipalpis is the main vector of the protozoan parasite Leishmania infantum, which is the causal agent of American Visceral Leishmaniasis. This insect uses male-produced pheromones for mate recognition. Elucidation of pheromone biogenesis or its regulation may enable molecular strategies for mating disruption and, consequently, the vector's population management. Motivated by our recent results of the transcriptomic characterization of the L. longipalpis pheromone gland, we performed a proteomic analysis of this tissue combining SDS-PAGE, and mass spectrometry followed by an integrative data analysis. Considering that annotated genome sequences of this sand fly are not available, we designed an alternative workflow searching MS/MS data against two customized databases using three search engines: Mascot, OMSSA and ProLuCID. A total of 542 proteins were confidently characterized, 445 of them using a Uniref100-insect protein database, and 97 using a transcript translated database. In addition, use of PEAKS for de novo peptide sequencing of MS/MS data confirmed ~90% identifications made with the combination of the three search engines. Our results include the identification of six of the seven enzymes of the mevalonate-pathway, plus the enzymes involved in sesquiterpenoid biosynthesis, all of which are proposed to be involved in pheromone production in L. longipalpis.

BIOLOGICAL SIGNIFICANCE

L. longipalpis is the main vector of the protozoan parasite L. infantum, which is the causal agent of American Visceral Leishmaniasis. One of the control measures of such disease is focused on vector population control. As this insect uses male-produced pheromones for mate recognition, the elucidation of pheromone biogenesis or its regulating process may enable molecular strategies for mating disruption and, consequently, this vector's population management. On this regard, in this manuscript we report expression evidence, at the protein level, of several molecules potentially involved in the pheromone production of L. longipalpis. Our results include the identification of the mevalonate-pathway enzymes, plus the enzymes involved in sesquiterpenoid biosynthesis, all of which are proposed to be involved in pheromone production in L. longipalpis. In addition, considering that the annotated genome sequences of this sand fly are not yet available, we designed an alternative workflow searching MS/MS data against proteomic and transcript translated customized databases, using three search engines: Mascot, OMSSA, and ProLuCID. In addition, a de novo peptide sequencing software (PEAKS) was used to further analyze the MS/MS data. This approach made it possible to identify and annotate 542 proteins for the pheromone gland of L. longipalpis. Importantly, all annotated protein sequences and raw data are available for the research community in protein repositories that provide free access to the data.

De novo transcriptome assembly from fat body and flight muscles transcripts to identify morph-specific gene expression profiles in Gryllus firmus

Wing polymorphism is a powerful model for examining many aspects of adaptation. The wing dimorphic cricket species, Gryllus firmus, consists of a long-winged morph with functional flight muscles that is capable of flight, and two flightless morphs. One (obligately) flightless morph emerges as an adult with vestigial wings and vestigial flight muscles. The other (plastic) flightless morph emerges with fully-developed wings but later in adulthood histolyzes its flight muscles. Importantly both flightless morphs have substantially increased reproductive output relative to the flight-capable morph. Much is known about the physiological and biochemical differences between the morphs with respect to adaptations for flight versus reproduction. In contrast, little is known about the molecular genetic basis of these morph-specific adaptations. To address this issue, we assembled a de novo transcriptome of G. firmus using 141.5 million Illumina reads generated from flight muscles and fat body, two organs that play key roles in flight and reproduction. We used the resulting 34,411 transcripts as a reference transcriptome for differential gene expression analyses. A comparison of gene expression profiles from functional flight muscles in the flight-capable morph versus histolyzed flight muscles in the plastic flight incapable morph identified a suite of genes involved in respiration that were highly expressed in pink (functional) flight muscles and genes involved in proteolysis highly expressed in the white (histolyzed) flight muscles. A comparison of fat body transcripts from the obligately flightless versus the flight-capable morphs revealed differential expression of genes involved in triglyceride biosynthesis, lipid transport, immune function and reproduction. These data provide a valuable resource for future molecular genetics research in this and related species and provide insight on the role of gene expression in morph-specific adaptations for flight versus reproduction.