Draft genome of the mountain pine beetle, Dendroctonus ponderosae Hopkins, a major forest pest

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.

Chromosomal-level genome assembly of Hylurgus ligniperda: insights into host adaptation and environmental tolerance

BACKGROUND

Hylurgus ligniperda (Coleoptera: Curculionidae) is a worldwide forest quarantine pest. It is widely distributed, has many host tree species, and possesses strong adaptability. To explore its environmental adaptability and the related molecular mechanisms, we conducted chromosome-level genome sequencing and analyzed the transcriptome under different environmental factors, identifying key expressed genes.

RESULTS

We employed PacBio, Illumina, and Hi-C sequencing techniques to assemble a 520 Mb chromosomal-level genome of H. ligniperda, obtaining an N50 of 39.97 Mb across 138 scaffolds. A total of 10,765 protein-coding genes were annotated after repeat masking. Fourteen chromosomes were identified, among which Hyli14 was determined to be the sex chromosome. Survival statistics were tested over various growth periods under high temperature and low humidity conditions. The maximum survival period of adults reached 292 days at 25 °C, 65% relative humidity. In comparison, the maximum survival period was 14 days under 35 °C, 65% relative humidity, and 106 days under 25°C, 40% relative humidity. This indicated that environmental stress conditions significantly reduced adults' survival period. We further conducted transcriptome analysis to screen for potentially influential differentially expressed genes, such as CYP450 and Histone. Subsequently, we performed gene family analysis to gain insights into their functions and interactions, such as CYP450 and Histone. CYP450 genes affected the detoxification metabolism of enzymes in the Cytochrome P450 pathway to adapt to different environments. Histone genes are involved in insect hormone biosynthesis and longevity-regulating pathways in H. ligniperda to adapt to environmental stress.

CONCLUSIONS

The genome at the chromosome level of H. ligniperda was assembled for the first time. The mortality of H. ligniperda increased significantly at 35 ℃, 65% RH, and 25 ℃, 40% RH. CYP450 and Histone genes played an important role in response to environmental stress. This genome offers a substantial genetic resource for investigating the molecular mechanisms behind beetle invasion and spread.

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.

Exploring the contribution of the salivary gland and midgut to digestion in the swede midge (Contarinia nasturtii) through a genomics-guided approach

The larvae of Contarinia nasturtii (Kieffer) (Diptera: Cecidomyiidae), the swede midge, targets the meristem of brassica crops where they induce the formation of galls and disrupt seed and vegetable production. Previously, we examined the salivary gland transcriptome of newly-hatched first instar larvae as they penetrated the host and initiated gall formation. Here we examine the salivary gland and midgut transcriptome of third instar larvae and provide evidence for cooperative nutrient acquisition beginning with secretion of enzymes and feeding facilitators followed by gastrointestinal digestion. Sucrose, presumably obtained from the phloem, appeared to be a major nutrient source as several α-glucosidases (sucrases, maltases) and β-fructofuranosidases (invertases) were identified. Genes encoding β-fructofuranosidases/invertases were among the most highly expressed in both tissues and represented two distinct gene families that may have originated via horizontal gene transfer from bacteria. The importance of the phloem as a nutrient source is underscored by the expression of genes encoding regucalcin and ARMET (arginine-rich mutated in early stages of tumor) which interfere with calcium signalling and prevent sieve tube occlusion. Lipids, proteins, and starch appear to serve as a secondary nutrient sources. Genes encoding enzymes involved in the detoxification of glucosinolates (myrosinases, arylsulfatases, and glutathione-S-transferases) were expressed indicative of Brassicaceae host specialization. The midgut expressed simple peritrophins and mucins typical of those found in Type II peritrophic matrices, the first such description for a gall midge.

Genome assembly of the southern pine beetle (Dendroctonus frontalis Zimmerman) reveals the origins of gene content reduction in Dendroctonus

Dendroctonus frontalis, also known as southern pine beetle (SPB), represents the most damaging forest pest in the southeastern United States. Strategies to predict, monitor and suppress SPB outbreaks have had limited success. Genomic data are critical to inform on pest biology and to identify molecular targets to develop improved management approaches. Here, we produced a chromosome-level genome assembly of SPB using long-read sequencing data. Synteny analyses confirmed the conservation of the core coleopteran Stevens elements and validated the bona fide SPB X chromosome. Transcriptomic data were used to obtain 39,588 transcripts corresponding to 13,354 putative protein-coding loci. Comparative analyses of gene content across 14 beetle and 3 other insects revealed several losses of conserved genes in the Dendroctonus clade and gene gains in SPB and Dendroctonus that were enriched for loci encoding membrane proteins and extracellular matrix proteins. While lineage-specific gene losses contributed to the gene content reduction observed in Dendroctonus, we also showed that widespread misannotation of transposable elements represents a major cause of the apparent gene expansion in several non-Dendroctonus species. Our findings uncovered distinctive features of the SPB gene complement and disentangled the role of biological and annotation-related factors contributing to gene content variation across beetles.

Comparative transcriptomics reveals the conservation and divergence of reproductive genes across three sympatric Tomicus bark beetles

Three tree-killing bark beetles belonging to the genus Tomicus, Tomicus yunnanensis, Tomicus brevipilosus and Tomicus minor (Coleoptera; Curculionidae, Scolytinae), are serious wood-borers with larvae feeding on the phloem tissues of Pinus yunnanensis. The three Tomicus beetles, in some cases, coexist in a same habitat, providing a best system for exploring the conservation and divergence of reproductive genes. Here, we applied comparative transcriptomics and molecular biology approaches to characterize reproductive-related genes in three sympatric Tomicus species. Illumina sequencing of female and male reproductive systems and residual bodies generated a large number of clean reads, representing 185,920,232 sequences in T. yunnanensis, 169,153,404 in T. brevipilosus and 178,493,176 in T. minor that were assembled into 32,802, 56,912 and 33,670 unigenes, respectively. The majority of the genes had detectable expression in reproductive tissues (FPKM >1), particularly those genes in T. brevipilosus accounting for 76.61 % of the total genes. From the transcriptomes, totally 838 genes encoding 463 detoxification enzymes, 339 chemosensory membrane proteins and 36 ionotropic glutamate receptors (iGluRs) were identified, including 622 reproductive tissue-expressed genes. Of these, members of carboxylesterases (COEs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs) and iGluRs were highly conserved in gene numbers and sequence identities across three Tomicus species. Further, expression profiling analyses revealed a number of genes expressed in reproductive tissues and the diverse expression characteristics in these beetles. The results provide evidence for the conservation and differences of reproductive genes among three sympatric closely related beetles, helping understand their different reproductive strategies and the maximization of the reproductive success.

Chromosome-level genome assembly of the Japanese sawyer beetle Monochamus alternatus

The Japanese sawyer beetle Monochamus alternatus (Coleoptera: Cerambycidae) is a pest in pine forests and acts as a vector for the pine wood nematode Bursaphelenchus xylophilus, which causes the pine wilt disease. We assembled a high-quality genome of M. alternatus at the chromosomal level using Illumina, Nanopore, and Hi-C sequencing technologies. The assembled genome is 767.12 Mb, with a scaffold N50 of 82.0 Mb. All contigs were assembled into ten pseudo-chromosomes. The genome contains 63.95% repeat sequences. We identify 16, 284 protein-coding genes in the genome, of which 11,244 were functionally annotated. The high-quality genome of M. alternatus provides an invaluable resource for the biological, ecological, and genetic study of this beetle and opens new avenues for understanding the transmission of pine wood nematode by insect vectors.

The genome of the invasive and broadly polyphagous Diaprepes root weevil, Diaprepes abbreviatus (Coleoptera), reveals an arsenal of putative polysaccharide-degrading enzymes

The Diaprepes root weevil (DRW), Diaprepes abbreviatus, is a broadly polyphagous invasive pest of agriculture in the southern United States and the Caribbean. Its genome was sequenced, assembled, and annotated to study genomic correlates of specialized plant-feeding and invasiveness and to facilitate the development of new methods for DRW control. The 1.69 Gb D. abbreviatus genome assembly was distributed across 653 contigs, with an N50 of 7.8 Mb and the largest contig of 62 Mb. Most of the genome was comprised of repetitive sequences, with 66.17% in transposable elements, 5.75% in macrosatellites, and 2.06% in microsatellites. Most expected orthologous genes were present and fully assembled, with 99.5% of BUSCO genes present and 1.5% duplicated. One hundred and nine contigs (27.19 Mb) were identified as putative fragments of the X and Y sex chromosomes, and homology assessment with other beetle X chromosomes indicated a possible sex chromosome turnover event. Genome annotation identified 18,412 genes, including 43 putative horizontally transferred (HT) loci. Notably, 258 genes were identified from gene families known to encode plant cell wall degrading enzymes and invertases, including carbohydrate esterases, polysaccharide lyases, and glycoside hydrolases (GH). GH genes were unusually numerous, with 239 putative genes representing 19 GH families. Interestingly, several other beetle species with large numbers of GH genes are (like D. abbreviatus) successful invasive pests of agriculture or forestry.

Genome and transcriptome of Ips nitidus provide insights into high-altitude hypoxia adaptation and symbiosis

Ips nitidus is a well-known conifer pest that has contributed significantly to spruce forest disturbance in the Qinghai-Tibet Plateau and seriously threatens the ecological balance of these areas. We report a chromosome-level genome of I. nitidus determined by PacBio and Hi-C technology. Phylogenetic inference showed that it diverged from the common ancestor of I. typographus ∼2.27 mya. Gene family expansion in I. nitidus was characterized by DNA damage repair and energy metabolism, which may facilitate adaptation to high-altitude hypoxia. Interestingly, differential gene expression analysis revealed upregulated genes associated with high-altitude hypoxia adaptation and downregulated genes associated with detoxification after feeding and tunneling in fungal symbiont Ophiostoma bicolor-colonized substrates. Our findings provide evidence of the potential adaptability of I. nitidus to conifer host, high-altitude hypoxia and insight into how fungal symbiont assist in this process. This study enhances our understanding of insect adaptation, symbiosis, and pest management.

Evolution of five environmentally responsive gene families in a pine-feeding sawfly, Neodiprion lecontei (Hymenoptera: Diprionidae)

A central goal in evolutionary biology is to determine the predictability of adaptive genetic changes. Despite many documented cases of convergent evolution at individual loci, little is known about the repeatability of gene family expansions and contractions. To address this void, we examined gene family evolution in the redheaded pine sawfly Neodiprion lecontei, a noneusocial hymenopteran and exemplar of a pine-specialized lineage evolved from angiosperm-feeding ancestors. After assembling and annotating a draft genome, we manually annotated multiple gene families with chemosensory, detoxification, or immunity functions before characterizing their genomic distributions and molecular evolution. We find evidence of recent expansions of bitter gustatory receptor, clan 3 cytochrome P450, olfactory receptor, and antimicrobial peptide subfamilies, with strong evidence of positive selection among paralogs in a clade of gustatory receptors possibly involved in the detection of bitter compounds. In contrast, these gene families had little evidence of recent contraction via pseudogenization. Overall, our results are consistent with the hypothesis that in response to novel selection pressures, gene families that mediate ecological interactions may expand and contract predictably. Testing this hypothesis will require the comparative analysis of high-quality annotation data from phylogenetically and ecologically diverse insect species and functionally diverse gene families. To this end, increasing sampling in under-sampled hymenopteran lineages and environmentally responsive gene families and standardizing manual annotation methods should be prioritized.

Rearranging and completing the puzzle: Phylogenomic analysis of bark beetles Dendroctonus reveals new hypotheses about genus diversification

Studies carried out on bark beetles within Dendroctonus have been extensive and revealed diverse information in different areas of their natural history, taxonomy, evolution, and interactions, among others. Despite these efforts, phylogenetic hypotheses have remained obscured mainly due to limited information analyzed (taxonomic, gene sampling, or both) in studies focused on obtaining evolutionary hypotheses for this genus. With the aim of filling these gaps in the evolutionary history for Dendroctonus, we analyzed ∼1800 loci mapped to a reference genome obtained for 20 of the 21 species recognized to date, minimizing the impact of missing information and improving the assumption of orthology in a phylogenomic framework. We obtained congruent phylogenetic topologies from two phylogenomic inference strategies: loci concatenation (ML framework) and a multispecies coalescent model (MSC) through the analysis of site pattern frequencies (SNPs). Dendroctonus is composed of two major clades (A and B), each containing five and four subclades, respectively. According to our divergence dating analysis, the MRCA for Dendroctonus dates back to the early Eocene, while the MRCA for each major clade diverged in the mid-Eocene. Interestingly, most of the speciation events of extant species occurred during the Miocene, which could be correlated with the diversification of pine trees (Pinus). The MRCA for Dendroctonus inhabited large regions of North America, with all ancestors and descendants of clade A having diversified within this region. The Mexican Transition Zone is important in the diversification processes for the majority of clade A species. For clade B, we identified two important colonization events to the Old World from America: the first in the early Oligocene from the Arctic to Asia (via Beringia), and the second during the Miocene from the Arctic-Western-Alleghany region to Europe and Siberia (also via Beringia). Our genomic analyses also supported the existence of hidden structured lineages within the frontalis complex, and also that D. beckeri represent a lineage independent from D. valens, as previously suggested. The information presented here updates the knowledge concerning the diversification of a genus with remarkable ecological and economic importance.

Omics in the Red Palm Weevil Rhynchophorus ferrugineus (Olivier) (Coleoptera: Curculionidae): A Bridge to the Pest

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is the most devastating pest of palm trees worldwide. Mitigation of the economic and biodiversity impact it causes is an international priority that could be greatly aided by a better understanding of its biology and genetics. Despite its relevance, the biology of the RPW remains poorly understood, and research on management strategies often focuses on outdated empirical methods that produce sub-optimal results. With the development of omics approaches in genetic research, new avenues for pest control are becoming increasingly feasible. For example, genetic engineering approaches become available once a species's target genes are well characterized in terms of their sequence, but also population variability, epistatic interactions, and more. In the last few years alone, there have been major advances in omics studies of the RPW. Multiple draft genomes are currently available, along with short and long-read transcriptomes, and metagenomes, which have facilitated the identification of genes of interest to the RPW scientific community. This review describes omics approaches previously applied to RPW research, highlights findings that could be impactful for pest management, and emphasizes future opportunities and challenges in this area of research.

Insight into weevil biology from a reference quality genome of the boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae)

The boll weevil, Anthonomus grandis grandis Boheman, is one of the most historically impactful insects due to its near destruction of the US cotton industry in the early 20th century. Contemporary efforts to manage this insect primarily use pheromone baited traps for detection and organophosphate insecticides for control, but this strategy is not sustainable due to financial and environmental costs. We present a high-quality boll weevil genome assembly, consisting of 306 scaffolds with approximately 24,000 annotated genes, as a first step in the identification of gene targets for novel pest control. Gene content and transposable element distribution are similar to those found in other Curculionidae genomes; however, this is the most contiguous and only assembly reported to date for a member in the species-rich genus Anthonomus. Transcriptome profiles across larval, pupal, and adult life stages led to identification of several genes and gene families that could present targets for novel control strategies.

Odorant-Binding and Chemosensory Proteins in Anthonomus eugenii (Coleoptera: Curculionidae) and Their Tissue Expression

The pepper weevil Anthonomus eugenii is one of the most damaging pests to the pepper crop. To offer alternative management strategies to insecticides, several studies have identified the semiochemicals that are involved in the pepper weevil's aggregation and mating behavior; however, there is no information on its perireceptor molecular mechanism, to date. In this study, bioinformatics tools were used to functionally annotate and characterize the A. eugenii head transcriptome and their probable coding proteins. We identified twenty-two transcripts belonging to families related to chemosensory processes, seventeen corresponding to odorant-binding proteins (OBP), and six to chemosensory proteins (CSP). All results matched with closely related Coleoptera: Curculionidae homologous proteins. Likewise, twelve OBP and three CSP transcripts were experimentally characterized by RT-PCR in different female and male tissues. The results by sex and tissue display the different expression patterns of the AeugOBPs and AeugCSPs; some are present in both sexes and all tissues, while others show expressions with higher specificity, which suggests diverse physiological functions in addition to chemo-detection. This study provides information to support the understanding of odor perception in the pepper weevil.

Role of Insect and Mammal Glutathione Transferases in Chemoperception

Glutathione transferases (GSTs) are ubiquitous key enzymes with different activities as transferases or isomerases. As key detoxifying enzymes, GSTs are expressed in the chemosensory organs. They fulfill an essential protective role because the chemosensory organs are located in the main entry paths of exogenous compounds within the body. In addition to this protective function, they modulate the perception process by metabolizing exogenous molecules, including tastants and odorants. Chemosensory detection involves the interaction of chemosensory molecules with receptors. GST contributes to signal termination by metabolizing these molecules. By reducing the concentration of chemosensory molecules before receptor binding, GST modulates receptor activation and, therefore, the perception of these molecules. The balance of chemoperception by GSTs has been shown in insects as well as in mammals, although their chemosensory systems are not evolutionarily connected. This review will provide knowledge supporting the involvement of GSTs in chemoperception, describing their localization in these systems as well as their enzymatic capacity toward odorants, sapid molecules, and pheromones in insects and mammals. Their different roles in chemosensory organs will be discussed in light of the evolutionary advantage of the coupling of the detoxification system and chemosensory system through GSTs.

Evolutionary ecology of the bark beetles Ips typographus and Pityogenes chalcographus

Ips typographus (L.) and Pityogenes chalcographus (L.) (Coleoptera: Curculionidae) are two common bark beetle species on Norway spruce in Eurasia. Multiple biotic and abiotic factors affect the life cycles of these two beetles, shaping their ecology and evolution. In this article, we provide a comprehensive and comparative summary of selected life-history traits. We highlight similarities and differences in biotic factors, like host range, interspecific competition, host colonization, reproductive behaviour and fungal symbioses. Moreover, we focus on the species' responses to abiotic factors and compare their temperature-dependent development and flight behaviour, cold adaptations and diapause strategies. Differences in biotic and abiotic traits might be the result of recent, species-specific evolutionary histories, particularly during the Pleistocene, with differences in glacial survival and postglacial recolonization. Finally, we discuss future research directions to understand ecological and evolutionary pathways of the two bark beetle species, for both basic research and applied forest management.

Genome-wide screening and expression of glutathione S-transferase genes reveal that GSTe4 contributes to sensitivity against β-cypermethrin in Zeugodacus cucurbitae

Glutathione S-transferases (GSTs) are an essential multifunctional protein family with common detoxifying enzymes. In this study, 34 GST genes were identified from the melon fly, Zeugodacus cucurbitae, one of the most destructive pests worldwide. These GSTs include 32 cytosolic genes and two microsomal genes. Furthermore, these cytosolic GSTs were classified into six classes: 11 delta, 13 epsilon, three theta, one sigma, two zeta, and two omega. Most of these showed dynamic expression during the developmental stage, some of which showed stage-specific expression. The expression in various adult tissues showed that most of them were expressed in anti-stress-related tissues. The transcriptional response of the delta and epsilon families was determined when Z. cucurbitae was exposed to three insecticides, abamectin, dinotefuran, and β-cypermethrin. Seven genes were significantly up-regulated by abamectin exposure. Moreover, five and four genes were significantly up-regulated with dinotefuran and β-cypermethrin exposure, respectively, demonstrating their involvement in the detoxification of these such toxic substances in Z. cucurbitae. One example of these genes, ZcGSTe4, was randomly selected to explore its function in response to β-cypermethrin exposure. Over-expressed ZcGSTe4 in E. coli showed significant tolerance to β-cypermethrin, and RNAi-mediated suppression of ZcGSTe4 also increased the sensitivity of melon fly to this agent. This study provides a foundation for further studies on the mechanism of detoxification metabolism in the melon fly.

Tolerance to dietary linalool primarily involves co-expression of cytochrome P450s and cuticular proteins in Pagiophloeus tsushimanus (Coleoptera: Curculionidae) larvae using SMRT sequencing and RNA-seq

BACKGROUND

Pagiophloeus tsushimanus (Coleoptera: Curculionidae), an emerging forest pest exclusively infesting camphor trees, has recently caused severe ecological and economic damage in localized areas in China. Its population outbreak depends largely on the capacity to overcome the pressure of terpenoid-derived metabolites (e.g. linalool) from camphor trees. At present, the molecular basis of physiological adaptation of P. tsushimanus to dietary linalool is poorly understood, and there is no available reference genome or transcriptome.

RESULTS

Herein, we constructed the transcriptome profiling of P. tsushimanus larvae reared on linalool-infused diets using RNA sequencing and single-molecule real-time sequencing. A total of 20,325 high-quality full-length transcripts were identified as a reference transcriptome, of which 14,492 protein-coding transcripts including 130 transcription factors (TFs), and 5561 long non-coding RNAs (lncRNAs) were detected. Also, 30 alternative splicing events and 8049 simple sequence repeats were captured. Gene ontology enrichment of differential expressed transcripts revealed that overall up-regulation of both cytochrome P450s (CYP450s) and cuticular proteins (CPs), was the primary response characteristic against dietary linalool. Other physiological effects possibly caused by linalool exposure, such as increase in Reactive Oxygen Species (ROS) and hormetic stimulation, were compensated by a handful of induced genes encoding antioxidases, heat shock proteins (HSPs), juvenile hormone (JH) epoxide hydrolases, and digestive enzymes. Additionally, based on co-expression networks analysis, a diverse array of hub lncRNAs and TFs co-expressed with CYP450s and CPs were screened as the potential gene regulators. Temporal expression of candidate transcripts determined by quantitative real-time PCR also indicated a cooperative relationship between the inductions of CYP450s and CPs upon exposure to linalool.

CONCLUSIONS

Our present study provides an important transcriptome resource of P. tsushimanus, and lays a valuable foundation for understanding how this specialist pest copes with chemical challenges in its specific host environments.

Experimental Life History Evolution Results in Sex-specific Evolution of Gene Expression in Seed Beetles

The patterns of reproductive timing and senescence vary within and across species owing to differences in reproductive strategies, but our understanding of the molecular underpinnings of such variation is incomplete. This is perhaps particularly true for sex differences. We investigated the evolution of sex-specific gene expression associated with life history divergence in replicated populations of the seed beetle Acanthoscelides obtectus, experimentally evolving under (E)arly or (L)ate life reproduction for >200 generations which has resulted in strongly divergent life histories. We detected 1,646 genes that were differentially expressed in E and L lines, consistent with a highly polygenic basis of life history evolution. Only 30% of differentially expressed genes were similarly affected in males and females. The evolution of long life was associated with significantly reduced sex differences in expression, especially in non-reproductive tissues. The expression differences were overall more pronounced in females, in accordance with their greater phenotypic divergence in lifespan. Functional enrichment analysis revealed differences between E and L beetles in gene categories previously implicated in aging, such as mitochondrial function and defense response. The results show that divergent life history evolution can be associated with profound changes in gene expression that alter the transcriptome in a sex-specific way, highlighting the importance of understanding the mechanisms of aging in each sex.

Genome-wide and expression-profiling analyses of the cytochrome P450 genes in Tenebrionidea

Cytochrome P450 monooxygenases (CYPs) are present in almost all areas of the tree of life. As one of the largest and most diverse superfamilies of multifunctional enzymes, they play important roles in the metabolism of xenobiotics and biosynthesis of endogenous compounds, shaping the success of insects. In this study, the CYPome (an omics term for all the CYP genes in a genome) diversification was examined in the four Tenebrionidea species through genome-wide analysis. A total of 483 CYP genes were identified, of which 103, 157, 122, and 101 were respectively deciphered from the genomes of Tebebrio molitor, Asbolus verucosus, Hycleus cichorii and Hycleus phaleratus. These CYPs were classified into four major clans (mitochondrial, CYP2, CYP3, and CYP4), and clans CYP3 and CYP4 are most diverse. Phylogenetic analysis showed that most CYPs of these Tenebrionidea beetles from each clan had a very close 1:1 orthology to each other, suggesting that they originate closely and have evolutionally conserved function. Expression analysis at different developmental stages and in various tissues showed the life stage-, gut-, salivary gland-, fat body-, Malpighian tubule-, antennae-, ovary- and testis-specific expression patterns of T. molitor CYP genes, implying their various potential roles in development, detoxification, immune response, digestion, olfaction, and reproduction. Our studies provide a platform to understand the evolution of Tenebrionidea CYP gene superfamily, and a basis for further functional investigation of the T. molitor CYPs involved in various biological processes.

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

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

Identification and characterization of candidate detoxification genes in Pharsalia antennata Gahan (Coleoptera: Cerambycidae)

The wood-boring beetles, including the majority of Cerambycidae, have developed the ability to metabolize a variety of toxic compounds derived from host plants and the surrounding environment. However, detoxification mechanisms underlying the evolutionary adaptation of a cerambycid beetle Pharsalia antennata to hosts and habitats are largely unexplored. Here, we characterized three key gene families in relation to detoxification (cytochrome P450 monooxygenases: P450s, carboxylesterases: COEs and glutathione-S-transferases: GSTs), by combinations of transcriptomics, gene identification, phylogenetics and expression profiles. Illumina sequencing generated 668,701,566 filtered reads in 12 tissues of P. antennata, summing to 100.28 gigabases data. From the transcriptome, 215 genes encoding 106 P450s, 77 COEs and 32 GSTs were identified, of which 107 relatives were differentially expressed genes. Of the identified 215 genes, a number of relatives showed the orthology to those in Anoplophora glabripennis, revealing 1:1 relationships in 94 phylogenetic clades. In the trees, P. antennata detoxification genes mainly clustered into one or two subfamilies, including 64 P450s in the CYP3 clan, 33 COEs in clade A, and 20 GSTs in Delta and Epsilon subclasses. Combining transcriptomic data and PCR approaches, the numbers of detoxification genes expressed in abdomens, antennae and legs were 188, 148 and 141, respectively. Notably, some genes exhibited significantly sex-biased levels in antennae or legs of both sexes. The findings provide valuable reference resources for further exploring xenobiotics metabolism and odorant detection in P. antennata.

BmSuc1 Affects Silk Properties by Acting on Sericin1 in Bombyx mori

BmSuc1, a novel animal-type β-fructofuranosidase (β-FFase, EC 3.2.1.26) encoding gene, was cloned and identified for the first time in the silkworm, Bombyx mori. BmSuc1 was specifically and highly expressed in the midgut and silk gland of Bombyx mori. Until now, the function of BmSuc1 in the silk gland was unclear. In this study, it was found that the expression changes of BmSuc1 in the fifth instar silk gland were consistent with the growth rate of the silk gland. Next, with the aid of the CRISPR/Cas9 system, the BmSuc1 locus was genetically mutated, and homozygous mutant silkworm strains with truncated β-FFase (BmSUC1) proteins were established. BmSuc1 mutant larvae exhibited stunted growth and decreased body weight. Interestingly, the molecular weight of part of Sericin1 (Ser1) in the silk gland of the mutant silkworms was reduced. The knockout of BmSuc1 reduced the sericin content in the silkworm cocoon shell, and the mechanical properties of the mutant line silk fibers were also negatively affected. These results reveal that BmSUC1 is involved in the synthesis of Ser1 protein in silk glands and helps to maintain the homeostasis of silk protein content in silk fibers and the mechanical properties of silk fibers, laying a foundation for the study of BmSUC1 regulation of silk protein synthesis in silk glands.

Chromosome-level genome assembly and population genomic analyses provide insights into adaptive evolution of the red turpentine beetle, Dendroctonus valens

BACKGROUND

Biological invasions are responsible for substantial environmental and economic losses. The red turpentine beetle (RTB), Dendroctonus valens LeConte, is an important invasive bark beetle from North America that has caused substantial tree mortality in China. The lack of a high-quality reference genome seriously limits deciphering the extent to which genetic adaptions resulted in a secondary pest becoming so destructive in its invaded area.

RESULTS

Here, we present a 322.41 Mb chromosome-scale reference genome of RTB, of which 98% of assembled sequences are anchored onto fourteen linkage groups including the X chromosome with a N50 size of 24.36 Mb, which is significantly greater than other Coleoptera species. Repetitive sequences make up 45.22% of the genome, which is higher than four other Coleoptera species, i.e., Mountain pine beetle Dendroctonus ponderosae, red flour beetle Tribolium castaneum, blister beetle Hycleus cichorii, and Colorado potato beetle Leptinotarsa decemlineata. We identify rapidly expanded gene families and positively selected genes in RTB, which may be responsible for its rapid environmental adaptation. Population genetic structure of RTB was revealed by genome resequencing of geographic populations in native and invaded regions, suggesting substantial divergence of the North American population and illustrates the possible invasion and spread route in China. Selective sweep analysis highlighted the enhanced ability of Chinese populations in environmental adaptation.

CONCLUSIONS

Overall, our high-quality reference genome represents an important resource for genomics study of invasive bark beetles, which will facilitate the functional study and decipher mechanism underlying invasion success of RTB by integrating the Pinus tabuliformis genome.

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 the forest insect pest Pissodes strobi reveals genome expansion and evidence of a Wolbachia endosymbiont

The highly diverse insect family of true weevils, Curculionidae, includes many agricultural and forest pests. Pissodes strobi, commonly known as the spruce weevil or white pine weevil, is a major pest of spruce and pine forests in North America. Pissodes strobi larvae feed on the apical shoots of young trees, causing stunted growth and can destroy regenerating spruce or pine forests. Here, we describe the nuclear and mitochondrial Pissodes strobi genomes and their annotations, as well as the genome of an apparent Wolbachia endosymbiont. We report a substantial expansion of the weevil nuclear genome, relative to other Curculionidae species, possibly driven by an abundance of class II DNA transposons. The endosymbiont observed belongs to a group (supergroup A) of Wolbachia species that generally form parasitic relationships with their arthropod host.

Glutathione S-Transferase Genes are Involved in Lambda-Cyhalothrin Resistance in Cydia pomonella via Sequestration

Pest management is mostly accomplished by the use of insecticides. However, the overuse of insecticides has led to the development of resistance. Glutathione S-transferases (GSTs) are vital detoxification enzymes involved in insecticide resistance in insects. In this study, we report the involvement of GSTs in insecticide resistance to lambda-cyhalothrin in Cydia pomonella, a globally quarantined fruit pest. A total of 25 GST, including 22 cytosolic genes and 3 microsomal genes, are identified from the genome database of C. pomonella. These cytosolic genes are further classified into six classes, including four in delta, eight in epsilon, three in omega, three in sigma, one in theta, and one in zeta class, as well as two unclassified genes. The real-time quantitative polymerase chain reaction (RT-qPCR) shows that the majority of these genes are mainly expressed throughout the larval stage and in the midgut of the fourth-instar larvae. Exposure to an LD₁₀ dose of lambda-cyhalothrin resulted in the upregulation of 17 GST genes. Moreover, mRNA levels of most GST genes, with the exception of CpGSTe6, CpGSTd2, CpGSTd4, and CpGSTz1, are considerably higher in a lambda-cyhalothrin-resistant population (ZW_R) than those of susceptible strains. Recombinant CpGSTd1, CpGSTd3, CpGSTe3, and CpGSTs2 can bind and metabolize lambda-cyhalothrin, with the highest metabolic rate observed for CpGSTd3 but no metabolite(s) was detected, supporting the role of GSTs in sequestration of lambda-cyhalothrin. Molecular dynamics simulation analysis indicates that key residues of hydrophobic pocket-derived lipophilic energy S(lipo) interactions with a hydrophobic pharmacophore of lambda-cyhalothrin are crucial for metabolism by CpGSTd3 and further lead to resistance. Our study is the first to experimentally confirm the involvement of GSTs in lambda-cyhalothrin resistance via sequestration and provides new insights into resistance management in C. pomonella.

Genomic insight into the scale specialization of the biological control agent Novius pumilus (Weise, 1892)

BACKGROUND

Members of the genus Novius Mulsant, 1846 (= Rodolia Mulsant, 1850) (Coleoptera, Coccinellidae), play important roles in the biological control of cotton cushion scale pests, especially those belonging to Icerya. Since the best-known species, the vedalia beetle Novius cardinalis (Mulsant, 1850) was introduced into California from Australia, more than a century of successful use in classical biological control, some species of Novius have begun to exhibit some field adaptations to novel but related prey species. Despite their economic importance, relatively little is known about the underlying genetic adaptations associated with their feeding habits. Knowledge of the genome sequence of Novius is a major step towards further understanding its biology and potential applications in pest control.

RESULTS

We report the first high-quality genome sequence for Novius pumilus (Weise, 1892), a representative specialist of Novius. Computational Analysis of gene Family Evolution (CAFE) analysis showed that several orthogroups encoding chemosensors, digestive, and immunity-related enzymes were significantly expanded (P < 0.05) in N. pumilus compared to the published genomes of other four ladybirds. Furthermore, some of these orthogroups were under significant positive selection pressure (P < 0.05). Notably, transcriptome profiling demonstrated that many genes among the significantly expanded and positively selected orthogroups, as well as genes related to detoxification were differentially expressed, when N. pumilus feeding on the nature prey Icerya compared with the no feeding set. We speculate that these genes are vital in the Icerya adaptation of Novius species.

CONCLUSIONS

We report the first Novius genome thus far. In addition, we provide comprehensive transcriptomic resources for N. pumilus. The results from this study may be helpful for understanding the association of the evolution of genes related to chemosensing, digestion, detoxification and immunity with the prey adaptation of insect predators. This will provide a reference for future research and utilization of Novius in biological control programs. Moreover, understanding the possible molecular mechanisms of prey adaptation also inform mass rearing of N. pumilus and other Novius, which may benefit pest control.

Genome resequencing reveals rapid, repeated evolution in the Colorado potato beetle

Insecticide resistance and rapid pest evolution threatens food security and the development of sustainable agricultural practices, yet the evolutionary mechanisms that allow pests to rapidly adapt to control tactics remains unclear. Here we examine how a global super-pest, the Colorado potato beetle (CPB), Leptinotarsa decemlineata, rapidly evolves resistance to insecticides. Using whole genome resequencing and transcriptomic data focused on its ancestral and pest range in North America, we assess evidence for three, non-mutually exclusive models of rapid evolution: pervasive selection on novel mutations, rapid regulatory evolution, and repeated selection on standing genetic variation. Population genomic analysis demonstrates that CPB is geographically structured, even among recently established pest populations. Pest populations exhibit similar levels of nucleotide diversity, relative to non-pest populations, and show evidence of recent expansion. Genome scans provide clear signatures of repeated adaptation across CPB populations, with especially strong evidence of selection on insecticide resistance genes in different populations. Analyses of gene expression show that constitutive upregulation of candidate insecticide resistance genes drives distinctive population patterns. CPB evolves insecticide resistance repeatedly across agricultural regions, leveraging similar genetic pathways but different genes, demonstrating a polygenic trait architecture for insecticide resistance that can evolve from standing genetic variation. Despite expectations, we do not find support for strong selection on novel mutations, or rapid evolution from selection on regulatory genes. These results suggest that integrated pest management practices must mitigate the evolution of polygenic resistance phenotypes among local pest populations, in order to maintain the efficacy and sustainability of novel control techniques.

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.

Phylogenomic and mitogenomic data can accelerate inventorying of tropical beetles during the current biodiversity crisis

Conservation efforts must be evidence-based, so rapid and economically feasible methods should be used to quantify diversity and distribution patterns. We have attempted to overcome current impediments to the gathering of biodiversity data by using integrative phylogenomic and three mtDNA fragment analyses. As a model, we sequenced the Metriorrhynchini beetle fauna, sampled from ~700 localities in three continents. The species-rich dataset included ~6,500 terminals, ~1,850 putative species delimited at 5% uncorrected pairwise threshold, possibly ~1,000 of them unknown to science. Neither type of data could alone answer our questions on biodiversity and phylogeny. The phylogenomic backbone enabled the integrative delimitation of robustly defined natural genus-group units that will inform future research. Using constrained mtDNA analysis, we identified the spatial structure of species diversity, very high species-level endemism, and a biodiversity hotspot in New Guinea. We suggest that focused field research and subsequent laboratory and bioinformatic workflow steps would substantially accelerate the inventorying of any hyperdiverse tropical group with several thousand species. The outcome would be a scaffold for the incorporation of further data from environmental sequencing and ecological studies. The database of sequences could set a benchmark for the spatiotemporal evaluation of biodiversity, would support evidence-based conservation planning, and would provide a robust framework for systematic, biogeographic, and evolutionary studies.

Identification of genes and gene expression associated with dispersal capacity in the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae)

Dispersal flights by the mountain pine beetle have allowed range expansion and major damage to pine stands in western Canada. We asked what the genetic and transcriptional basis of mountain pine beetle dispersal capacity is. Using flight mills, RNA-seq and a targeted association study, we compared strong-flying, weak-flying, and non-flying female beetles from the recently colonized northern end of their range. Nearly 3,000 genes were differentially expressed between strong and weak flying beetles, while weak fliers and nonfliers did not significantly differ. The differentially expressed genes were mainly associated with lipid metabolism, muscle maintenance, oxidative stress response, detoxification, endocrine function, and flight behavior. Three variant loci, two in the coding region of genes, were significantly associated with flight capacity but these genes had no known functional link to flight. Several differentially expressed gene systems may be important for sustained flight, while other systems are downregulated during dispersal and likely to conserve energy before host colonization. The candidate genes and SNPs identified here will inform further studies and management of mountain pine beetle, as well as contribute to understanding the mechanisms of insect dispersal flights.

The transposable element-rich genome of the cereal pest Sitophilus oryzae

BACKGROUND

The rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

RESULTS

We sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino acids required for insect development and cuticle biosynthesis.

CONCLUSIONS

Here we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.

Draft genome of the blister beetle, Epicauta chinensis

Existence of cantharidin (CTD) in blister beetles is a significant ecological adaptive mechanism that defends against predators and regulates courtship and mating behaviors. To better understand CTD biosynthetic information as well as its biology and pharmacology, we assembled a genome of 151.88 Mb for Epicauta chinensis using PacBio sequencing technology. Gene annotation yielded 249,238 repeats, 527 non-coding RNAs and 12,520 protein-coding genes. Compared to other 11 insects, expansions of gene families in E. chinensis for most core gene families likely associated with environmental adaptation, such as chemoreception, immunity, and detoxification. We further annotated P450s and immune-related genes, a total of 117 putative P450s comprising 7 CYP2, 67 CYP3, 36 CYP4, and 7 mitochondrial P450s and 281 immune-related genes were identified. Comparative analysis of the insect immune repertoires indicated presence of immune genes detected only from Coleopteran insects such as MD2-like. This suggested a lineage-specific gene evolution for Coleopteran insects. Based on the gene family evolution analysis, we identified two probable candidate genes including CYP4TT1 and phytanoyl-CoA dioxygenase for CTD biosynthesis. The high-quality reference genome of E. chinensis provides the genetic basis for further investigation of CTD biosynthesis and in-depth studies of the development and evolution of blister beetles.

Mutually Exclusive Expression of Closely Related Odorant-Binding Proteins 9A and 9B in the Antenna of the Red Flour Beetle Tribolium castaneum

Olfaction is crucial for insects to find food sources, mates, and oviposition sites. One of the initial steps in olfaction is facilitated by odorant-binding proteins (OBPs) that translocate hydrophobic odorants through the aqueous olfactory sensilla lymph to the odorant receptor complexes embedded in the dendritic membrane of olfactory sensory neurons. The Tribolium castaneum (Coleoptera, Tenebrionidae) OBPs encoded by the gene pair TcasOBP9A and TcasOBP9B represent the closest homologs to the well-studied Drosophila melanogaster OBP Lush (DmelOBP76a), which mediates pheromone reception. By an electroantennographic analysis, we can show that these two OBPs are not pheromone-specific but rather enhance the detection of a broad spectrum of organic volatiles. Both OBPs are expressed in the antenna but in a mutually exclusive pattern, despite their homology and gene pair character by chromosomal location. A phylogenetic analysis indicates that this gene pair arose at the base of the Cucujiformia, which dates the gene duplication event to about 200 Mio years ago. Therefore, this gene pair is not the result of a recent gene duplication event and the high sequence conservation in spite of their expression in different sensilla is potentially the result of a common function as co-OBPs.

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.

Whole-Genome Sequence of the Trypoxylus dichotomus Japanese rhinoceros beetle

The draft whole-genome sequence of the Japanese rhinoceros beetle, Trypoxylus dichotomus was obtained using long-read PacBio sequence technology. The final assembled genome consisted of 739 Mbp in 2,347 contigs, with 24.5× mean coverage and a G+C content of 35.99%.

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.

Chromosome-Level Genome Assembly of Cyrtotrachelus buqueti and Mining of Its Specific Genes

Background: The most severe insect damage to bamboo shoots is the bamboo-snout beetle (Cyrtotrachelus buqueti). Bamboo is a perennial plant that has significant economic value. C. buqueti also plays a vital role in the degradation of bamboo lignocellulose and causing damage. The genome sequencing and functional gene annotation of C. buqueti are of great significance to reveal the molecular mechanism of its efficient degradation of bamboo fiber and the development of the bamboo industry.

Results: The size of C. buqueti genome was close to 600.92 Mb by building a one paired-end (PE) library and k-mer analysis. Then, we developed nine 20-kb SMRTbell libraries for genome sequencing and got a total of 51.12 Gb of the original PacBio sequel reads. Furthermore, after filtering with a coverage depth of 85.06×, clean reads with 48.71 Gb were obtained. The final size of C. buqueti genome is 633.85 Mb after being assembled and measured, and the contig N50 of C. buqueti genome is 27.93 Mb. The value of contig N50 shows that the assembly quality of C. buqueti genome exceeds that of most published insect genomes. The size of the gene sequence located on chromosomes reaches 630.86 Mb, accounting for 99.53% of the genome sequence. A 1,063 conserved genes were collected at this assembled genome, comprising 99.72% of the overall genes with 1,066 using the Benchmark Uniform Single-Copy Orthology (BUSCO). Moreover, 63.78% of the C. buqueti genome is repetitive, and 57.15% is redundant with long-term elements. A 12,569 protein-coding genes distributed on 12 chromosomes were acquired after function annotation, of which 96.18% were functional genes. The comparative genomic analysis results revealed that C. buqueti was similar to D. ponderosae. Moreover, the comparative analysis of specific genes in C. buqueti genome showed that it had 244 unique lignocellulose degradation genes and 240 genes related to energy production and conversion. At the same time, 73 P450 genes and 30 GST genes were identified, respectively, in the C. buqueti genome.

Conclusion: The high-quality C. buqueti genome has been obtained in the present study. The assembly level of this insect’s genome is higher than that of other most reported insects’ genomes. The phylogenetic analysis of P450 and GST gene family showed that C. buqueti had a vital detoxification function to plant chemical components.

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.

A highly-contiguous genome assembly of the Eurasian spruce bark beetle, Ips typographus, provides insight into a major forest pest

Conifer-feeding bark beetles are important herbivores and decomposers in forest ecosystems. These species complete their life cycle in nutritionally poor substrates and some can kill enormous numbers of trees during population outbreaks. The Eurasian spruce bark beetle (Ips typographus) can destroy >100 million m3 of spruce in a single year. We report a 236.8 Mb I. typographus genome assembly using PacBio long-read sequencing. The final phased assembly has a contig N50 of 6.65 Mb in 272 contigs and is predicted to contain 23,923 protein-coding genes. We reveal expanded gene families associated with plant cell wall degradation, including pectinases, aspartyl proteases, and glycosyl hydrolases. This genome sequence from the genus Ips provides timely resources to address questions about the evolutionary biology of the true weevils (Curculionidae), one of the most species-rich animal families. In forests of today, increasingly stressed by global warming, this draft genome may assist in developing pest control strategies to mitigate outbreaks.

Comparative genome and transcriptome analyses reveal innate differences in response to host plants by two color forms of the two-spotted spider mite Tetranychus urticae

Background

The two-spotted spider mite, Tetranychus urticae, is a major agricultural pest with a cosmopolitan distribution, and its polyphagous habits provide a model for investigating herbivore-plant interactions. There are two body color forms of T. urticae with a different host preference. Comparative genomics and transcriptomics are used here to investigate differences in responses of the forms to host plants at the molecular level. Biological responses of the two forms sourced from multiple populations are also presented.

Results

We carried out principal component analysis of transcription changes in three red and three green T. urticae populations feeding on their original host (common bean), and three hosts to which they were transferred: cotton, cucumber and eggplant. There were differences among the forms in gene expression regardless of their host plant. In addition, different changes in gene expression were evident in the two forms when responding to the same host transfer. We further compared biological performance among populations of the two forms after feeding on each of the four hosts. Fecundity of 2-day-old adult females showed a consistent difference between the forms after feeding on bean. We produced a 90.1-Mb genome of the red form of T. urticae with scaffold N50 of 12.78 Mb. Transcriptional profiles of genes associated with saliva, digestion and detoxification showed form-dependent responses to the same host and these genes also showed host-specific expression effects.

Conclusions

Our research revealed that forms of T. urticae differ in host-determined transcription responses and that there is form-dependent plasticity in the transcriptomic responses. These differences may facilitate the extreme polyphagy shown by spider mites, although fitness differences on hosts are also influenced by population differences unrelated to color form.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07894-7.

Genome Assembly of the Cold-Tolerant Leaf Beetle Gonioctena quinquepunctata, an Important Resource for Studying Its Evolution and Reproductive Barriers between Species

Coleoptera is the most species-rich insect order, yet is currently underrepresented in genomic databases. An assembly was generated for ca. 1.7 Gb genome of the leaf beetle Gonioctena quinquepunctata by first assembling long-sequence reads (Oxford Nanopore; ± 27-fold coverage) and subsequently polishing the resulting assembly with short sequence reads (Illumina; ± 85-fold coverage). The unusually large size (most Coleoptera species are associated with a reported size below 1 Gb) was at least partially attributed to the presence of a large fraction of repeated elements (73.8%). The final assembly was characterized by an N50 length of 432 kb and a BUSCO score of 95.5%. The heterozygosity rate was ± 0.6%. Automated genome annotation informed by RNA-Seq resulted in 40,568 predicted proteins, which is much larger than the typical range 17,000–23,000 predicted for other Coleoptera. However, no evidence of a genome duplication was detected. This new reference genome will contribute to our understanding of genetic variation in the Coleoptera. Among others, it will also allow exploring reproductive barriers between species, investigating introgression in the nuclear genome, and identifying genes involved in resistance to extreme climate conditions.

Comparative analysis of mite genomes reveals positive selection for diet adaptation

Diet is a powerful evolutionary force for species adaptation and diversification. Acari is one of the most abundant clades of Arachnida, exhibiting diverse dietary types, while the underlying genetic adaptive mechanisms are not fully understood. Based on comparative analyses of 15 Acari genomes, we found genetic bases for three specialized diets. Herbivores experienced stronger selection pressure than other groups; the olfactory genes and gene families involving metabolizing toxins showed strong adaptive signals. Genes and gene families related to anticoagulation, detoxification, and haemoglobin digestion were found to be under strong selection pressure or significantly expanded in the blood-feeding species. Lipid metabolism genes have a faster evolutionary rate and been subjected to greater selection pressures in fat-feeding species; one positively selected site in the fatty-acid amide hydrolases 2 gene was identified. Our research provides a new perspective for the evolution of Acari and offers potential target loci for novel pesticide development.

Comparative analyses of transcriptional responses of Dectes texanus LeConte (Coleoptera: Cerambycidae) larvae fed on three different host plants and artificial diet

Dectes texanus is an important coleopteran pest of soybeans and cultivated sunflowers in the Midwestern United States that causes yield losses by girdling stems of their host plants. Although sunflower and giant ragweed are primary hosts of D. texanus, they began colonizing soybeans approximately 50 years ago and no reliable management method has been established to prevent or reduce losses by this pest. To identify genes putatively involved when feeding soybean, we compared gene expression of D. texanus third-instar larvae fed soybean to those fed sunflower, giant ragweed, or artificial diet. Dectes texanus larvae differentially expressed 514 unigenes when fed on soybean compared to those fed the other diet treatments. Enrichment analyses of gene ontology terms from up-regulated unigenes in soybean-fed larvae compared to those fed both primary hosts highlighted unigenes involved in oxidoreductase and polygalacturonase activities. Cytochrome P450s, carboxylesterases, major facilitator superfamily transporters, lipocalins, apolipoproteins, glycoside hydrolases 1 and 28, and lytic monooxygenases were among the most commonly up-regulated unigenes in soybean-fed larvae compared to those fed their primary hosts. These results suggest that D. texanus larvae differentially expressed unigenes involved in biotransformation of allelochemicals, digestion of plant cell walls and transport of small solutes and lipids when feeding in soybean.

De novo transcriptome analysis and identification of reproduction control genes from the red palm weevil Rhynchophorus ferrugineus

Recent attacks by the red palm weevil, Rhynchophorus ferrugineus (Olivier), have become a severe problem for palm species. In present work, fat body transcriptome of adult female red palm weevil was analyzed, focusing on the identification of reproduction control genes. Transcriptome study was completed by means of next-generation sequencing (NGS) using Illumina Hiseq 2000 sequencing system. A total of 105,938,182 raw reads, 102,645,544 clean reads, and 9,238,098,960 clean nucleotides with a guanine-cytosine content of 40.31%, were produced. The processed transcriptome data resulted in 43,789 unique transcripts (with mean lengths of 1,172 bp). It was found that 20% of total unique transcripts shared up to 80%-100% sequence identity with homologous species, mainly the mountain pine beetle Dendroctonus ponderosae (59.9%) and red flour beetle Tribolium castaneum (26.9%). Nearly 25 annotated genes were predicted to be involved in red palm weevil reproduction, including five vitellogenin (Vg) transcripts. Among the five Vg gene transcripts, one was highly expressed compared with the other four (FPKM values of 1.963, 1.471, 1.028, and 1.017, respectively), and the five Vg gene transcripts were designated as RfVg, RfVg-equivalent1, RfVg-equivalent2, RfVg-equivalent3, and RfVg-equivalent4, respectively. The high expression level of RfVg verified by RT-polymerase chain reaction analysis suggested that RfVg is the primary functional Vg gene in red palm weevil. A high similarity of RfVg with other Coleopterans was also reflected in a phylogenetic tree, where RfVg was placed within the clade of the order Coleoptera. Awareness of the major genes that play critical roles in reproduction and proliferation of red palm weevil is valuable to understand their reproduction mechanism at a molecular level. In addition, for future molecular studies, the NGS dataset obtained will be useful and will promote the exploration of biotech-based control strategies against red palm weevil, a primary pest of palm trees.

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 (~147 Mb) 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.

Genomic Signature of Shifts in Selection in a Subalpine Ant and Its Physiological Adaptations

Understanding how organisms adapt to extreme environments is fundamental and can provide insightful case studies for both evolutionary biology and climate-change biology. Here, we take advantage of the vast diversity of lifestyles in ants to identify genomic signatures of adaptation to extreme habitats such as high altitude. We hypothesized two parallel patterns would occur in a genome adapting to an extreme habitat: 1) strong positive selection on genes related to adaptation and 2) a relaxation of previous purifying selection. We tested this hypothesis by sequencing the high-elevation specialist Tetramorium alpestre and four other phylogenetically related species. In support of our hypothesis, we recorded a strong shift of selective forces in T. alpestre, in particular a stronger magnitude of diversifying and relaxed selection when compared with all other ants. We further disentangled candidate molecular adaptations in both gene expression and protein-coding sequence that were identified by our genome-wide analyses. In particular, we demonstrate that T. alpestre has 1) a higher level of expression for stv and other heat-shock proteins in chill-shock tests and 2) enzymatic enhancement of Hex-T1, a rate-limiting regulatory enzyme that controls the entry of glucose into the glycolytic pathway. Together, our analyses highlight the adaptive molecular changes that support colonization of high-altitude environments.

Jekyll or Hyde? The genome (and more) of Nesidiocoris tenuis, a zoophytophagous predatory bug that is both a biological control agent and a pest

Nesidiocoris tenuis (Reuter) is an efficient predatory biological control agent used throughout the Mediterranean Basin in tomato crops but regarded as a pest in northern European countries. From the family Miridae, it is an economically important insect yet very little is known in terms of genetic information and no genomic or transcriptomic studies have been published. Here, we use a linked-read sequencing strategy on a single female N. tenuis. From this, we assembled the 355 Mbp genome and delivered an ab initio, homology-based and evidence-based annotation. Along the way, the bacterial "contamination" was removed from the assembly. In addition, bacterial lateral gene transfer (LGT) candidates were detected in the N. tenuis genome. The complete gene set is composed of 24 688 genes; the associated proteins were compared to other hemipterans (Cimex lectularis, Halyomorpha halys and Acyrthosiphon pisum). We visualized the genome using various cytogenetic techniques, such as karyotyping, CGH and GISH, indicating a karyotype of 2n = 32. Additional analyses include the localization of 18S rDNA and unique satellite probes as well as pooled sequencing to assess nucleotide diversity and neutrality of the commercial population. This is one of the first mirid genomes to be released and the first of a mirid biological control agent.

Horizontal Gene Transfer and Gene Duplication of β-Fructofuranosidase Confer Lepidopteran Insects Metabolic Benefits

Horizontal gene transfer (HGT) is a potentially critical source of material for ecological adaptation and the evolution of novel genetic traits. However, reports on posttransfer duplication in organism genomes are lacking, and the evolutionary advantages conferred on the recipient are generally poorly understood. Sucrase plays an important role in insect physiological growth and development. Here, we performed a comprehensive analysis of the evolution of insect β-fructofuranosidase transferred from bacteria via HGT. We found that posttransfer duplications of β-fructofuranosidase were widespread in Lepidoptera and sporadic occurrences of β-fructofuranosidase were found in Coleoptera and Hymenoptera. β-fructofuranosidase genes often undergo modifications, such as gene duplication, differential gene loss, and changes in mutation rates. Lepidopteran β-fructofuranosidase gene (SUC) clusters showed marked divergence in gene expression patterns and enzymatic properties in Bombyx mori (moth) and Papilio xuthus (butterfly). We generated SUC1 mutations in B. mori using CRISPR/Cas9 to thoroughly examine the physiological function of SUC. BmSUC1 mutant larvae were viable but displayed delayed growth and reduced sucrase activities that included susceptibility to the sugar mimic alkaloid found in high concentrations in mulberry. BmSUC1 served as a critical sucrase and supported metabolic homeostasis in the larval midgut and silk gland, suggesting that gene transfer of β-fructofuranosidase enhanced the digestive and metabolic adaptation of lepidopteran insects. These findings highlight not only the universal function of β-fructofuranosidase with a link to the maintenance of carbohydrate metabolism but also an underexplored function in the silk gland. This study expands our knowledge of posttransfer duplication and subsequent functional diversification in the adaptive evolution and lineage-specific adaptation of organisms.