A high‐quality chromosome‐level genome assembly of a generalist herbivore, Trichoplusia ni

Downregulation of APN1 and ABCC2 mutation in Bt Cry1Ac-resistant Trichoplusia ni are genetically independent

The resistance to the insecticidal protein Cry1Ac from the bacterium Bacillus thuringiensis (Bt) in the cabbage looper, Trichoplusia ni, has previously been identified to be associated with a frameshift mutation in the ABC transporter ABCC2 gene and with altered expression of the aminopeptidase N (APN) genes APN1 and APN6, shown as missing of the 110-kDa APN1 (phenotype APN1¯) in larval midgut brush border membrane vesicles (BBMV). In this study, genetic linkage analysis identified that the APN1¯ phenotype and the ABCC2 mutation in Cry1Ac-resistant T. ni segregated independently, although they were always associated under Cry1Ac selection. The ABCC2 mutation and APN1¯ phenotype were separated into two T. ni strains respectively. Bioassays of the T. ni strains with Cry1Ac determined that the T. ni with the APN1¯ phenotype showed a low level resistance to Cry1Ac (3.5-fold), and the associated resistance is incompletely dominant in the background of the ABCC2 mutation. Whereas the ABCC2 mutation-associated resistance to Cry1Ac is at a moderate level, and the resistance is incompletely recessive in the genetic background of downregulated APN1. Analysis of Cry1Ac binding to larval midgut BBMV indicated that the midgut in larvae with the APN1¯ phenotype had reduced binding affinity for Cry1Ac, but the number of binding sites remained unchanged, and the midgut in larvae with the ABCC2 mutation had both reduced binding affinity and reduced number of binding sites for Cry1Ac. The reduced Cry1Ac binding to BBMV from larvae with the ABCC2 mutation or APN1¯ phenotype correlated with the lower levels of resistance.IMPORTANCEThe soil bacterium Bacillus thuringiensis (Bt) is an important insect pathogen used as a bioinsecticide for pest control. Bt genes coding for insecticidal proteins are the primary transgenes engineered into transgenic crops (Bt crops) to confer insect resistance. However, the evolution of resistance to Bt proteins in insect populations in response to exposure to Bt threatens the sustainable application of Bt biotechnology. Cry1Ac is a major insecticidal toxin utilized for insect control. Genetic mechanisms of insect resistance to Cry1Ac are complex and require to be better understood. The resistance to Cry1Ac in Trichoplusia ni is associated with a mutation in the ABCC2 gene and also associated with the APN expression phenotype APN1¯. This study identified the genetic independence of the APN1¯ phenotype from the ABCC2 mutation and isolated and analyzed the ABCC2 mutation-associated and APN1¯ phenotype-associated resistance traits in T. ni to provide new insights into the genetic mechanisms of Cry1Ac resistance in insects.

Effect of sublethal concentrations of the bioinsecticide spinosyn treatment of Trichoplusia ni eggs on the caterpillar and its parasitoid, Trichogramma brassicae

BACKGROUND

Integrated Pest Management (IPM) seeks to combine multiple management strategies for optimal pest control. One method that is successfully employed in IPM is the use of beneficial organisms. However, in severe circumstances when pest insects exceed threshold limits, insecticides may still need to be implemented. Thus, understanding the effects of insecticides on biocontrol agents, such as parasitoid wasps, is paramount to ensure sustainable agroecosystems. Sublethal effects of the bioinsecticide spinosyn, a mixture of the bacterial Saccharopolyspora spinosa (Mertz and Yao) fermentation products spinosyn A and D, on eggs of Trichoplusia ni (Hübner), a cruciferous crop pest, and its egg parasitoid Trichogramma brassicae (Bezdenko) was investigated.

RESULTS

The LC50 for spinosyn A and D (dissolved in ethanol) on T. ni eggs is 54 ng mL-1. Transcriptomics on caterpillars (1st and 3rd instars) that hatched from eggs treated with sublethal concentrations of spinosyn identified the upregulation of several genes encoding proteins that may be involved in insecticide resistance including detoxification enzymes, such as cytochrome P450s, glutathione S-transferases and esterases. Sublethal T. ni egg treatments did not affect parasitoid emergence, however, there was a marked increase in the size of T. brassicae hind tibia and wings that emerged from spinosyn-treated eggs.

CONCLUSIONS

For the caterpillar, treatment of eggs with sublethal concentrations of spinosyn may induce insecticide resistance mechanisms. For the parasitoids, their increased size when reared in spinosyn-treated eggs suggests that the emerged wasps may have higher performance. © 2024 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Gene gain and loss from the Asian corn borer W chromosome

BACKGROUND

Sex-limited chromosomes Y and W share some characteristics, including the degeneration of protein-coding genes, enrichment of repetitive elements, and heterochromatin. However, although many studies have suggested that Y chromosomes retain genes related to male function, far less is known about W chromosomes and whether they retain genes related to female-specific function.

RESULTS

Here, we built a chromosome-level genome assembly of the Asian corn borer, Ostrinia furnacalis Guenée (Lepidoptera: Crambidae, Pyraloidea), an economically important pest in corn, from a female, including both the Z and W chromosome. Despite deep conservation of the Z chromosome across Lepidoptera, our chromosome-level W assembly reveals little conservation with available W chromosome sequence in related species or with the Z chromosome, consistent with a non-canonical origin of the W chromosome. The W chromosome has accumulated significant repetitive elements and experienced rapid gene gain from the remainder of the genome, with most genes exhibiting pseudogenization after duplication to the W. The genes that retain significant expression are largely enriched for functions in DNA recombination, the nucleosome, chromatin, and DNA binding, likely related to meiotic and mitotic processes within the female gonad.

CONCLUSIONS

Overall, our chromosome-level genome assembly supports the non-canonical origin of the W chromosome in O. furnacalis, which experienced rapid gene gain and loss, with the retention of genes related to female-specific function.

A chromosome-level genome assembly of Sesamia inferens

The pink stem borer, Sesamia inferens (Walker), is a significant polyphagous pest historically restricted to regions south of N34° latitude. However, with changes in global climate and farming practices, the distribution of this moth has progressively exceeded its traditional limit of 34° N and encompassed most regions in North China. The genetic adaptations of S. inferens remain incompletely understood due to the lack of high-quality genome resources. Here, we sequenced the genome of S. inferens using PacBio and Hi-C technology, yielding a genome assembly of 865.04 Mb with contig N50 of 1.23 Mb. BUSCO analysis demonstrated this genome assembly has a high-level completeness of 96.1% gene coverage. In total, 459.72 Mb repeat sequences (53.14% of the assembled genome) and 20858 protein-coding genes were identified. We used the Hi-C technique to anchor 1135 contigs to 31 chromosomes, yielding a chromosome-level genome assembly with a scaffold N50 of 29.99 Mb. In conclusion, our high-quality genome assembly provided valuable resource that exploring the genetic characteristics of local adaptation and developing an efficient control strategy.

Epigenetic weapons in plant-herbivore interactions: Sulforaphane disrupts histone deacetylases, gene expression, and larval development in Spodoptera exigua while the specialist feeder Trichoplusia ni is largely resistant to these effects

Cruciferous plants produce sulforaphane (SFN), an inhibitor of nuclear histone deacetylases (HDACs). In humans and other mammals, the consumption of SFN alters enzyme activities, DNA-histone binding, and gene expression within minutes. However, the ability of SFN to act as an HDAC inhibitor in nature, disrupting the epigenetic machinery of insects feeding on these plants, has not been explored. Here, we demonstrate that SFN consumed in the diet inhibits the activity of HDAC enzymes and slows the development of the generalist grazer Spodoptera exigua, in a dose-dependent fashion. After consuming SFN for seven days, the activities of HDAC enzymes in S. exigua were reduced by 50%. Similarly, larval mass was reduced by 50% and pupation was delayed by 2-5 days, with no additional mortality. Similar results were obtained when SFN was applied topically to eggs. RNA-seq analyses confirm that SFN altered the expression of thousands of genes in S. exigua. Genes associated with energy conversion pathways were significantly downregulated while those encoding for ribosomal proteins were dramatically upregulated in response to the consumption of SFN. In contrast, the co-evolved specialist feeder Trichoplusia ni was not negatively impacted by SFN, whether it was consumed in their diet at natural concentrations or applied topically to eggs. The activities of HDAC enzymes were not inhibited and development was not disrupted. In fact, SFN exposure sometimes accelerated T. ni development. RNA-seq analyses revealed that the consumption of SFN alters gene expression in T. ni in similar ways, but to a lesser degree, compared to S. exigua. This apparent resistance of T. ni can be overwhelmed by unnaturally high levels of SFN or by exposure to more powerful pharmaceutical HDAC inhibitors. These results demonstrate that dietary SFN interferes with the epigenetic machinery of insects, supporting the hypothesis that plant-derived HDAC inhibitors serve as "epigenetic weapons" against herbivores.

Chromosomal-level Genome Assembly of the Coffee Bee Hawk Moth Reveals the Evolution of Chromosomes and the Molecular Basis of Distinct Phenotypes

Cephonodes hylas, the coffee bee hawk moth is a hawk moth species with unique characteristics, such as larvae feeding on gardenia, overcoming the toxicity of its iridoid glycosides, diurnal adults, and transparent wings. Although C. hylas is a fascinating model for molecular biological research, genome sequence analysis-based genetic approaches to elucidate these peculiarities have not yet been undertaken. We successfully achieved de novo genome assembly at the chromosome level of C. hylas comparable to the Lepidoptera model organism, silkworm. Additionally, 16,854 protein-coding genes were annotated, and the constructed genome sequence and annotated genes were of the highest quality BUSCO completion compared to closely related species. Comparative genome analysis revealed the process of chromosomal evolution from the Bombycoidea ancestral (n = 31) genome and changes in turnover at the chromosome level associated with chromosomal fusion events, such as the rate of repetitive sequence insertion. These analyses were only possible because the genome was constructed at the chromosome level. Additionally, increased the nonsynonymous/synonymous rate (dN/dS) ratios were observed in multiple photoreceptor-related genes that were strongly associated with the acquisition of diurnal activity. Furthermore, tandemly duplicated expanded genes containing many digestive and other enzymes and larval midgut-specific expression were also confirmed. These genes may be involved in the metabolism of genipin, a toxin found in gardenias. Using the genome sequence of C. hylas determined at the chromosome level, we have successfully identified new insights into the chromosomal evolution of Bombycoidea, as well as the relationship between the genome sequence and its characteristic traits.

Impact of dual-baculovirus infection on the Sf9 insect cell transcriptome during rAAV production using single-cell RNA-seq

The insect cell-baculovirus expression vector system (IC-BEVS) has shown to be a powerful platform to produce complex biopharmaceutical products, such as recombinant proteins and virus-like particles. More recently, IC-BEVS has also been used as an alternative to produce recombinant adeno-associated virus (rAAV). However, little is known about the variability of insect cell populations and the potential effect of heterogeneity (e.g., stochastic infection process and differences in infection kinetics) on product titer and/or quality. In this study, transcriptomics analysis of Sf9 insect cells during the production of rAAV of serotype 2 (rAAV2) using a low multiplicity of infection, dual-baculovirus system was performed via single-cell RNA-sequencing (scRNA-seq). Before infection, the principal source of variability in Sf9 insect cells was associated with the cell cycle. Over the course of infection, an increase in transcriptional heterogeneity was detected, which was linked to the expression of baculovirus genes as well as to differences in rAAV transgenes (rep, cap and gfp) expression. Noteworthy, at 24 h post-infection, only 29.4% of cells enclosed all three necessary rAAV transgenes to produce packed rAAV2 particles, indicating limitations of the dual-baculovirus system. In addition, the trajectory analysis herein performed highlighted that biological processes such as protein folding, metabolic processes, translation, and stress response have been significantly altered upon infection. Overall, this work reports the first application of scRNA-seq to the IC-BEVS and highlights significant variations in individual cells within the population, providing insight into the rational cell and process engineering toward improved rAAV2 production in IC-BEVS.

Ultra-deep sequencing of 45S rDNA to discern intragenomic diversity in three Chrysodeixis species for molecular identification

Species identification is necessary to prevent introductions of exotic plant pests through global trade. Many of these pests are understudied and lack publicly available DNA sequence data on which rapid molecular identification methods can be based. One such lineage is the genus Chrysodeixis, which includes three species of potential concern for United States trade initiatives: C. includens, C. chalcites, and C. eriosoma. Here we describe a method to generate robust 45S rDNA profiles using long read sequencing in order to clarify evolutionary relationships and develop a real-time PCR identification technique. Such an identification tool will be useful in rapidly differentiating between Chrysodeixis species of quarantine concern where traditional morphological identification methods are insufficient. Molecular methods such as this greatly reduce the time spent identifying each specimen, allow for detection of eDNA, vastly increase throughput, and increase the probability of detection. The methods presented here will be generally adaptable to any understudied lepidopteran taxa that necessitates a molecular diagnostic assay and, with adjustment or testing of the primers, could be applied to any group for which development of rDNA profiles in a benchtop setting would prove useful.

Increased Phenoloxidase Activity Constitutes the Main Defense Strategy of Trichoplusia ni Larvae against Fungal Entomopathogenic Infections

The cabbage looper Trichoplusia ni is an important agricultural pest worldwide and is frequently used as a model organism for assessing entomopathogenic fungi virulence, though few studies have measured the host response repertoire to fungal biocontrol agents. Here, we quantified the immune response of T. ni larvae following exposure to two entomopathogenic fungal species: Beauveria bassiana and Cordyceps javanica. Results from our study demonstrate that T. ni larvae exposed to fungal entomopathogens had higher total phenoloxidase activity compared to controls, indicating that the melanization cascade is one of the main immune components driving defense against fungal infection and contrasting observations from other insect-fungi interaction studies. We also observed differences in host response depending on the species of entomopathogenic fungi, with significantly higher induction observed during infections with B. bassiana than with C. javanica. Larvae exposed to B. bassiana had an increased expression of genes involved in prophenoloxidase response and the Imd, JNK, and Jak/STAT immune signaling pathways. Our results indicate a notable absence of Toll pathway-related responses, further contrasting results to other insect-fungi pathosystems. Important differences were also observed in the induction of antimicrobial effectors, with B. bassiana infections eliciting three antimicrobial effectors (lysozyme, gloverin, and cecropin), while C. javanica only induced cecropin expression. These results provide insight into the host response strategies employed by T. ni for protection against entomopathogenic fungi and increase our understanding of insect-fungal entomopathogen interactions, aiding in the design of more effective microbial control strategies for this important agricultural pest.

Chloride intracellular channel gene knockdown induces insect cell lines death and level increases of intracellular calcium ions

Chloride intracellular channel (CLIC) is a member of the chloride channel protein family for which growing evidence supports a pivotal role in fundamental cellular events. However, the physiological function of CLIC in insects is still rarely uncovered. The ovary-derived High Five (Hi-5) cell line isolated from the cabbage looper (Trichoplusia ni) is widely used in laboratories. Here, we studied both characteristics and functions of CLIC in Hi-5 cells (TnCLIC). We identified the TnCLIC gene in Hi-5 cells and annotated highly conserved CLIC proteins in most insect species. After RNA interference of TnCLIC, the phenomenon of significantly increased cell death suggests that the TnCLIC protein is essential for the survival of Hi-5 cells. The same lethal effect was also observed in Spodoptera frugiperda 9 and Drosophila melanogaster Schneider 2 cells after CLIC knockdown. Furthermore, we found that this kind of cell death was accompanied by increases in intracellular calcium ions after TnCLIC knockdown with the transcriptomic analyses and the detection of calcium levels. Our results provide insights into insect CLIC as a key factor for cell survival and lay the foundation for the cell death mechanism.

Reference Genome Sequences of the Oriental Armyworm, Mythimna separata (Lepidoptera: Noctuidae)

Lepidopteran insects are an important group of animals, including those used as biochemical and physiological model species in the insect and silk industries as well as others that are major agricultural pests. Therefore, the genome sequences of several lepidopteran insects have been reported. The oriental armyworm, Mythimna separata, is an agricultural pest commonly used to study insect immune reactions and interactions with parasitoid wasps as hosts. To improve our understanding of these research topics, reference genome sequences were constructed in the present study. Using long-read and short-read sequence data, de novo assembly and polishing were performed and haplotigs were purged. Subsequently, gene predictions and functional annotations were performed. To search for orthologs of the Toll and Immune Deficiency (IMD) pathways and for C-type lectins, annotation data analysis, BLASTp, and Hummer scans were performed. The M. separata genome is 682 Mbp; its contig N50 was 2.7 Mbp, with 21,970 genes and 24,452 coding sites predicted. All orthologs of the core components of the Toll and IMD pathways and 105 C-type lectins were identified. These results suggest that the genome data were of sufficient quality for use as reference genome data and could contribute to promoting M. separata and lepidopteran research at the molecular and genome levels.

A database of crop pest cell lines

We have developed an online database describing the known cell lines from Coleoptera, Diptera, Hemiptera, Hymenoptera, and Lepidoptera that were developed from agricultural pests. Cell line information has been primarily obtained from previous compilations of insect cell lines. We conducted in-depth Internet literature searches and drew on Internet sources such as the Cellosaurus database (https://web.expasy.org/cellosaurus/), and inventories from cell line depositories. Here, we report on a new database of insect cell lines, which covers 719 cell lines from 86 species. We have not included cell lines developed from Drosophila because they are already known from published databases, such as https://dgrc.bio.indiana.edu/cells/Catalog. We provide the designation, tissue and species of origin, cell line developer, unique characteristics, its use in various applications, publications, and patents, and, when known, insect virus susceptibility. This information has been assembled and organized into a searchable database available at the link https://entomology.ca.uky.edu/aginsectcellsdatabase which will be updated on an ongoing basis.

Novel genetic basis of resistance to Bt toxin Cry1Ac in Helicoverpa zea

Crops genetically engineered to produce insecticidal proteins from the bacterium Bacillus thuringiensis have advanced pest management, but their benefits are diminished when pests evolve resistance. Elucidating the genetic basis of pest resistance to Bacillus thuringiensis toxins can improve resistance monitoring, resistance management, and the design of new insecticides. Here, we investigated the genetic basis of resistance to Bacillus thuringiensis toxin Cry1Ac in the lepidopteran Helicoverpa zea, one of the most damaging crop pests in the United States. To facilitate this research, we built the first chromosome-level genome assembly for this species, which has 31 chromosomes containing 375 Mb and 15,482 predicted proteins. Using a genome-wide association study, fine-scale mapping, and RNA-seq, we identified a 250-kb quantitative trait locus on chromosome 13 that was strongly associated with resistance in a strain of Helicoverpa zea that had been selected for resistance in the field and lab. The mutation in this quantitative trait locus contributed to but was not sufficient for resistance, which implies alleles in more than one gene contributed to resistance. This quantitative trait locus contains no genes with a previously reported role in resistance or susceptibility to Bacillus thuringiensis toxins. However, in resistant insects, this quantitative trait locus has a premature stop codon in a kinesin gene, which is a primary candidate as a mutation contributing to resistance. We found no changes in gene sequence or expression consistently associated with resistance for 11 genes previously implicated in lepidopteran resistance to Cry1Ac. Thus, the results reveal a novel and polygenic basis of resistance.

The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture

Simple Summary

In this study, we produced a genomic resource for the light brown apple moth, Epiphyas postvittana, to understand the biological basis of adaptation to a high number of hosts (polyphagy) and the invasive nature of this and other lepidopteran pests. The light brown apple moth is an invasive pest of horticultural plants, with over 500 recorded plant hosts. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, causing significant economic losses for fruit producers. Comparative genomic analyses with other lepidopteran genomes indicate that a high proportion of the genome is made up of repetitive sequences, with the majority of the known elements being DNA transposable elements and retrotransposons. Twenty gene families show significant expansions, including some likely to have a role in its pest status. Finally, population genomics, investigated by a RAD-tag approach, indicated likely patterns of invasion and admixture, with Californian moths most probably being derived from Australia.

Abstract

The light brown apple moth, Epiphyas postvittana is an invasive, polyphagous pest of horticultural systems around the world. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, where it has been found on over 500 plants, including many horticultural crops. We have produced a genomic resource, to understand the biological basis of the polyphagous and invasive nature of this and other lepidopteran pests. The assembled genome sequence encompassed 598 Mb and has an N50 of 301.17 kb, with a BUSCO completion rate of 97.9%. Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb). Of the 31,389 predicted genes, 28,714 (91.5%) were assigned to 11,438 orthogroups across the Lepidoptera, of which 945 were specific to E. postvittana. Twenty gene families showed significant expansions in E. postvittana, including some likely to have a role in its pest status, such as cytochrome p450s, glutathione-S-transferases and UDP-glucuronosyltransferases. Finally, using a RAD-tag approach, we investigated the population genomics of this pest, looking at its likely patterns of invasion.

A glimpse on metazoan ZNFX1 helicases, ancient players of antiviral innate immunity

The human zinc finger NFX1-type containing 1 (ZNFX1) is an interferon-stimulated protein associated to the outer mitochondrial membrane, able to bind dsRNAs and interact with MAVS proteins, promoting type I IFN response in the early stage of viral infection. An N-terminal Armadillo (ARM)-type fold and a large helicase core (P-loop) and zinc fingers confer RNA-binding and ATPase activities to ZNFX1. We studied the phylogenetic distribution of metazoan ZNFX1s, ZNFX1 gene expression trends and genomic and protein signatures during viral infection of invertebrates. Based on 221 ZNFX1 sequences, we obtained a polyphyletic tree with a taxonomy-consistent branching at the phylum-level only. In metazoan genomes, ZNFX1 genes were found either in single copy, with up to some tens of exons in vertebrates, or in multiple copies, with one or a few exons and one of them sometimes encompassing most of the coding sequence, in invertebrates like sponges, sea urchins and mollusks. Structural analyses of selected ZNFX1 proteins showed high conservation of the helicase region (P-loop), an overall conserved region and domain architecture, an ARM-fold mostly traceable, and the presence of intrinsically disordered regions of varying length and position. The remarkable over-expression of ZNFX1 in bivalve and gastropod mollusks infected with dsDNA viruses underscores the antiviral role of ZNFX1, whereas nothing similar was found in virus-infected nematodes and corals. Whether the functional diversification reported in the C. elegans ZNFX1 occurs in other metazoan proteins remains to be established.

Bt Cry1Ac resistance in Trichoplusia ni is conferred by multi-gene mutations

The three-domain Cry toxin Cry1Ac from Bacillus thuringiensis （Bt） is an important insecticidal toxin in Bt sprays and has been used in transgenic Bt-crops to confer insect resistance. The cabbage looper, Trichoplusia ni, has developed resistance to Bt sprays in commercial greenhouses, and the resistance to Cry1Ac has been previously identified to be associated with altered expression of the APN1 and APN6 genes and be genetically linked to a locus on chromosome 15. In this study, the Cry1Ac resistance locus in T. ni was further finely mapped, and the specific Cry1Ac resistance-conferring mutation in the resistance locus was identified to be a 4 bp frameshift insertion in the ABCC2 gene by whole genome resequencing, midgut transcriptome analysis, candidate gene cDNA sequencing and mutation site genomic DNA sequencing. By CRISPR/Cas9 mutagenesis, a series of ABCC2 and ABCC3 mutant T. ni strains were generated, and the role of ABCC2 in the toxicity of Cry1Ac in T. ni was confirmed. The results from this study also showed that knockout of ABCC2 in T. ni conferred resistance to Cry1Ac at a level lower than that in the greenhouse-derived resistant T. ni strain and that the Cry1Ac resistance-associated alteration of APN1 and APN6 expression was independent of ABCC2 gene mutations, indicating that the altered expression of APN1 and APN6 was controlled by another gene mutation in Cry1Ac resistant T. ni. Furthermore, T. ni larval bioassays showed that the level of Cry1Ac resistance in F₁ families from reciprocal crosses of the Cry1Ac resistant strain with an ABCC2 knockout CRISPR strain was significantly higher than that in ABCC2 knockout strain, indicating the presence of additional Cry1Ac resistance-conferring mutation(s) in the Cry1Ac resistant strain. Therefore, the resistance to Cry1Ac in T. ni is conferred by a mutation in ABCC2 and an additional mutation (or mutations) which leads to altered expression of APN1 and APN6. The additional Cry1Ac resistance mutation or mutations remain to be identified.

Assessing the Impact of a Viral Infection on the Expression of Transposable Elements in the Cabbage Looper Moth (Trichoplusia ni)

Most studies of stress-induced transposable element (TE) expression have so far focused on abiotic sources of stress. Here, we analyzed the impact of an infection by the AcMNPV baculovirus on TE expression in a cell line (Tnms42) and midgut tissues of the cabbage looper moth (Trichoplusia ni). We find that a large fraction of TE families (576/636 in Tnms42 cells and 503/612 in midgut) is lowly expressed or not expressed at all [≤ 4 transcripts per million (TPM)] in the uninfected condition (median TPM of 0.37 in Tnms42 and 0.46 in midgut cells). In the infected condition, a total of 62 and 187 TE families were differentially expressed (DE) in midgut and Tnms42 cells, respectively, with more up- (46) than downregulated (16) TE families in the former and as many up- (91) as downregulated (96) TE families in the latter. Expression log2 fold changes of DE TE families varied from -4.95 to 9.11 in Tnms42 cells and from -4.28 to 7.66 in midgut. Large variations in expression profiles of DE TEs were observed depending on the type of cells and on time after infection. Overall, the impact of AcMNPV on TE expression in T. ni is moderate but potentially sufficient to affect TE activity and genome architecture. Interestingly, one host-derived TE integrated into AcMNPV genomes is highly expressed in infected Tnms42 cells. This result shows that virus-borne TEs can be expressed, further suggesting that they may be able to transpose and that viruses may act as vectors of horizontal transfer of TEs in insects.

Genome and transcriptome analysis of the beet armyworm Spodoptera exigua reveals targets for pest control

The genus Spodoptera (Lepidoptera: Noctuidae) includes some of the most infamous insect pests of cultivated plants including Spodoptera frugiperda, Spodoptera litura, and Spodoptera exigua. To effectively develop targeted pest control strategies for diverse Spodoptera species, genomic resources are highly desired. To this aim, we provide the genome assembly and developmental transcriptome comprising all major life stages of S. exigua, the beet armyworm. Spodoptera exigua is a polyphagous herbivore that can feed on > 130 host plants, including several economically important crops. The 419 Mb beet armyworm genome was sequenced from a female S. exigua pupa. Using a hybrid genome sequencing approach (Nanopore long-read data and Illumina short read), a high-quality genome assembly was achieved (N50 = 1.1 Mb). An official gene set (18,477 transcripts) was generated by automatic annotation and by using transcriptomic RNA-seq datasets of 18 S. exigua samples as supporting evidence. In-depth analyses of developmental stage-specific expression combined with gene tree analyses of identified homologous genes across Lepidoptera genomes revealed four potential genes of interest (three of them Spodoptera-specific) upregulated during first- and third-instar larval stages for targeted pest-outbreak management. The beet armyworm genome sequence and developmental transcriptome covering all major developmental stages provide critical insights into the biology of this devastating polyphagous insect pest species worldwide. In addition, comparative genomic analyses across Lepidoptera significantly advance our knowledge to further control other invasive Spodoptera species and reveals potential lineage-specific target genes for pest control strategies.

Extended in vivo transcriptomes of two ascoviruses with different tissue tropisms reveal alternative mechanisms for enhancing virus reproduction in hemolymph

Ascoviruses are large dsDNA viruses characterized by the extraordinary changes they induce in cellular pathogenesis and architecture whereby after nuclear lysis and extensive hypertrophy, each cell is cleaved into numerous vesicles for virion reproduction. However, the level of viral replication and transcription in vesicles compared to other host tissues remains uncertain. Therefore, we applied RNA-Sequencing to compare the temporal transcriptome of Spodoptera frugiperda ascovirus (SfAV) and Trichoplusia ni ascovirus (TnAV) at 7, 14, and 21 days post-infection (dpi). We found most transcription occurred in viral vesicles, not in initial tissues infected, a remarkably novel reproduction mechanism compared to all other viruses and most other intracellular pathogens. Specifically, the highest level of viral gene expression occurred in hemolymph, for TnAV at 7 dpi, and SfAV at 14 dpi. Moreover, we found that host immune genes were partially down-regulated in hemolymph, where most viral replication occurred in highly dense accumulations of vesicles.

Monitoring insect transposable elements in large double-stranded DNA viruses reveals host-to-virus and virus-to-virus transposition

The mechanisms by which transposable elements (TEs) can be horizontally transferred between animals are unknown, but viruses are possible candidate vectors. Here, we surveyed the presence of host-derived TEs in viral genomes in 35 deep sequencing data sets produced from 11 host–virus systems, encompassing nine arthropod host species (five lepidopterans, two dipterans, and two crustaceans) and six different double-stranded (ds) DNA viruses (four baculoviruses and two iridoviruses). We found evidence of viral-borne TEs in 14 data sets, with frequencies of viral genomes carrying a TE ranging from 0.01% to 26.33% for baculoviruses and from 0.45% to 7.36% for iridoviruses. The analysis of viral populations separated by a single replication cycle revealed that viral-borne TEs originating from an initial host species can be retrieved after viral replication in another host species, sometimes at higher frequencies. Furthermore, we detected a strong increase in the number of integrations in a viral population for a TE absent from the hosts’ genomes, indicating that this TE has undergone intense transposition within the viral population. Finally, we provide evidence that many TEs found integrated in viral genomes (15/41) have been horizontally transferred in insects. Altogether, our results indicate that multiple large dsDNA viruses have the capacity to shuttle TEs in insects and they underline the potential of viruses to act as vectors of horizontal transfer of TEs. Furthermore, the finding that TEs can transpose between viral genomes of a viral species sets viruses as possible new niches in which TEs can persist and evolve.

Proteomic Identification of Bombyx mori Organelles Using the Engineered Ascorbate Peroxidase APEX and Development of Silkworm Organelle Proteome Database (SilkOrganPDB)

Silkworm Bombyx mori is an economically important insect and a lepidopteran model. Organelle proteome is vital to understanding gene functions; however, it remains to be identified in silkworm. Here, using the engineered ascorbate peroxidase APEX, we constructed transgenic B. mori embryo cells (BmE) expressing APEX-NLS, COX4-APEX, APEX-Rev, and APEX-KDEL in nucleus, mitochondrial matrix (MM), cytosol, and endoplasmic reticulum (ER), and isolated the biotin-labeled proteins using streptavidin-affinity purification, respectively. The isolated proteins were determined using LC-MS/MS and annotated by searching B. mori genomes downloaded from GenBank, SilkBase, SilkDB 2.0, and SilkDB 3.0, resulting in 842, 495, 311, and 445 organelle proteins identified, respectively. We mapped the 296 MM proteins annotated in the GenBank data to mitochondrial protein databases of the fly, human, and mouse, and found that 140 (47%) proteins are homologous to 80 fly proteins, and 65 (22%) proteins match to 31 and 29 human and mouse proteins, respectively. Protein orthology was predicted in multiple insects using OrthoMCL, producing 460 families containing 839 proteins we identified. Out of 460 families, 363 were highly conserved and found in all insects, leaving only three proteins without orthology in other insects, indicating that the identified proteins are highly conserved and probably play important roles in insects. A gene ontology enrichment analysis by clusterProfiler revealed that the nucleus proteins significantly enriched in cellular component terms of nucleus and nucleolus, the MM proteins markedly enriched in molecular function terms of nucleotide binding, and the cytosol proteins mainly enriched in biological process terms of small molecule metabolism. To facilitate the usage and analysis of our data, we developed an open-access database, Silkworm Organelle Proteome Database (SilkOrganPDB), which provides multiple modules for searching, browsing, downloading, and analyzing these proteins, including BLAST, HMMER, Organelle Proteins, Protein Locations, Sequences, Gene Ontology, Homologs, and Phylogeny. In summary, our work revealed the protein composition of silkworm BmE organelles and provided a database resource helpful for understanding the functions and evolution of these proteins.

Draft nuclear genome and complete mitogenome of the Mediterranean corn borer, Sesamia nonagrioides, a major pest of maize

The Mediterranean corn borer (Sesamia nonagrioides, Noctuidae, Lepidoptera) is a major pest of maize in Europe and Africa. Here, we report an assembly of the nuclear and mitochondrial genome of a pool of inbred males and females third-instar larvae, based on short- and long-read sequencing. The complete mitochondrial genome is 15,330 bp and contains all expected 13 and 24 protein-coding and RNA genes, respectively. The nuclear assembly is 1021 Mb, composed of 2553 scaffolds and it has an N50 of 1105 kb. It is more than twice larger than that of all Noctuidae species sequenced to date, mainly due to a higher repeat content. A total of 17,230 protein-coding genes were predicted, including 15,776 with InterPro domains. We provide detailed annotation of genes involved in sex determination (doublesex, insulin-like growth factor 2 mRNA-binding protein, and P-element somatic inhibitor) and of alpha-amylase genes possibly involved in interaction with parasitoid wasps. We found no evidence of recent horizontal transfer of bracovirus genes from parasitoid wasps. These genome assemblies provide a solid molecular basis to study insect genome evolution and to further develop biocontrol strategies against S. nonagrioides.

A Tale of Two Transcriptomic Responses in Agricultural Pests via Host Defenses and Viral Replication

Autographa californica Multiple Nucleopolyhedrovirus (AcMNPV) is a baculovirus that causes systemic infections in many arthropod pests. The specific molecular processes underlying the biocidal activity of AcMNPV on its insect hosts are largely unknown. We describe the transcriptional responses in two major pests, Spodoptera frugiperda (fall armyworm) and Trichoplusia ni (cabbage looper), to determine the host-pathogen responses during systemic infection, concurrently with the viral response to the host. We assembled species-specific transcriptomes of the hemolymph to identify host transcriptional responses during systemic infection and assessed the viral transcript abundance in infected hemolymph from both species. We found transcriptional suppression of chitin metabolism and tracheal development in infected hosts. Synergistic transcriptional support was observed to suggest suppression of immune responses and induction of oxidative stress indicating disease progression in the host. The entire AcMNPV core genome was expressed in the infected host hemolymph with a proportional high abundance detected for viral transcripts associated with replication, structure, and movement. Interestingly, several of the host genes that were targeted by AcMNPV as revealed by our study are also targets of chemical insecticides currently used commercially to control arthropod pests. Our results reveal an extensive overlap between biological processes represented by transcriptional responses in both hosts, as well as convergence on highly abundant viral genes expressed in the two hosts, providing an overview of the host-pathogen transcriptomic landscape during systemic infection.

In vitro screening of 65 mycotoxins for insecticidal potential

The economic losses and threats to human and animal health caused by insects and the pathogens transmitted by them require effective and environmentally-friendly methods of controlling them. One such group of natural biocontrol agents which may be used as biopesticides is that of the entomopathogenic fungi and their toxic secondary metabolites (mycotoxins). The present in vitro work examined the insecticidal potential of 65 commercially-available mycotoxins against the insect Sf-9 cell line. Mammalian Caco-2 and THP-1 cell lines served as reference controls to select insecticidal mycotoxins harmless to mammalian cells. All tested mycotoxins significantly reduced the in vitro proliferation of the Sf-9 cells and evoked morphological changes. Ten of the mycotoxins found to strongly inhibit Sf-9 proliferation also had moderate or no effect on Caco-2 cells. The THP-1 cells were highly resistant to the tested mycotoxins: doses 10³ times higher were needed to affect viability and morphology (1 μg/ml for THP-1 versus 1 ng/ml for Sf-9 and Caco-2). Nine mycotoxins significantly decreased Sf-9 cell proliferation with minor effects on mammalian cells: cyclosporins B and D, cytochalasin E, gliotoxin, HC toxin, paxilline, penitrem A, stachybotrylactam and verruculogen. These may be good candidates for future biopesticide formulations.

InSexBase: an annotated genomic resource of sex chromosomes and sex-biased genes in insects

Sex determination and the regulation of sexual dimorphism are among the most fascinating topics in modern biology. As the most species-rich group of sexually reproducing organisms on Earth, insects have multiple sex determination systems. Though sex chromosomes and sex-biased genes are well-studied in dozens of insects, their gene sequences are scattered in various databases. Moreover, a shortage of annotation hinders the deep mining of these data. Here, we collected the chromosome-level sex chromosome data of 49 insect species, including 34 X chromosomes, 15 Z chromosomes, 5 W chromosomes and 2 Y chromosomes. We also obtained Y-linked contigs of four insects species—Anopheles gambiae, Drosophila innubila, Drosophila yakuba and Tribolium castaneum. The unannotated chromosome-level sex chromosomes were annotated using a standard pipeline, yielding a total of 123 030 protein-coding genes, 2 159 427 repeat sequences, 894 miRNAs, 1574 rRNAs, 5105 tRNAs, 395 snoRNAs (small nucleolar RNA), 54 snRNAs (small nuclear RNA) and 5959 other ncRNAs (non-coding RNA). In addition, 36 781 sex-biased genes were identified by analyzing 62 RNA-seq (RNA sequencing) datasets. Together with 5707 sex-biased genes from the Drosophila genus collected from the Sex-Associated Gene Database, we obtained a total of 42 488 sex-biased genes from 13 insect species. All these data were deposited into InSexBase, a new user-friendly database of insect sex chromosomes and sex-biased genes.

Database URL:http://www.insect-genome.com/Sexdb/.

Chromosome-level genome of the peach fruit moth Carposina sasakii (Lepidoptera: Carposinidae) provides a resource for evolutionary studies on moths

The peach fruit moth (PFM), Carposina sasakii Matsumura, is a major phytophagous orchard pest widely distributed across Northeast Asia. Here, we report the chromosome-level genome for the PFM, representing the first genome for the family Carposinidae, from the lepidopteran superfamily Copromorphoidea. The genome was assembled into 404.83 Mb sequences using PacBio long-read and Illumina short-read sequences, including 275 contigs, with a contig N50 length of 2.62 Mb. All contigs were assembled into 31 linkage groups assisted by the Hi-C technique, including 30 autosomes and a Z chromosome. BUSCO analysis showed that 98.3% of genes were complete and 0.4% of genes were fragmented, while 1.3% of genes were missing in the assembled genome. In total, 21,697 protein-coding genes were predicted, of which 84.80% were functionally annotated. Because of the importance of diapause triggered by photoperiod in PFM, five circadian genes in the PFM as well as in the other related species were annotated, and potential genes related to diapause and photoperiodic reaction were also identified from transcriptome sequencing. In addition, manual annotation of detoxification gene families was undertaken and showed a higher number of glutathione S-transferase (GST) gene in PFM than in most other lepidopterans, in contrast to a lower number of uridine diphosphate (UDP)-glycosyltransferase (UGT) gene, carboxyl/cholinesterases (CCE) gene and cytochrome P450 monooxygenase (P450) gene, suggesting different detoxication pathways in this moth. The high-quality genome provides a resource for comparative evolutionary studies of this moth and its relatives within the context of radiations across Lepidoptera.

Structures of diverse poxin cGAMP nucleases reveal a widespread role for cGAS-STING evasion in host-pathogen conflict

DNA viruses in the family Poxviridae encode poxin enzymes that degrade the immune second messenger 2'3'-cGAMP to inhibit cGAS-STING immunity in mammalian cells. The closest homologs of poxin exist in the genomes of insect viruses suggesting a key mechanism of cGAS-STING evasion may have evolved outside of mammalian biology. Here we use a biochemical and structural approach to discover a broad family of 369 poxins encoded in diverse viral and animal genomes and define a prominent role for 2'3'-cGAMP cleavage in metazoan host-pathogen conflict. Structures of insect poxins reveal unexpected homology to flavivirus proteases and enable identification of functional self-cleaving poxins in RNA-virus polyproteins. Our data suggest widespread 2'3'-cGAMP signaling in insect antiviral immunity and explain how a family of cGAS-STING evasion enzymes evolved from viral proteases through gain of secondary nuclease activity. Poxin acquisition by poxviruses demonstrates the importance of environmental connections in shaping evolution of mammalian pathogens.

Joyce A, Pereira de Souza A, Gabriela Murúa M, J. Clough S, Imaculada Zucchi M. A first draft genome of the Sugarcane borer, Diatraea saccharalis

Background: The sugarcane borer (Diatraea saccharalis), a widely distributed moth throughout the Americas, is a pest that affects economically important crops such as sugarcane, sorghum, wheat, maize and rice. Given its significant impact on yield reduction, whole-genome information of the species is needed. Here, we report the first draft assembly of the D. saccharalis genome. Methods: The genomic sequences were obtained using the Illumina HiSeq 2500 whole-genome sequencing of a single adult male specimen. We assembled the short-reads using the SPAdes software and predicted protein-coding genes using MAKER. Genome assembly completeness was assessed through BUSCO and the repetitive content by RepeatMasker. Results: The 453 Mb assembled sequences contain 1,445 BUSCO gene orthologs and 1,161 predicted gene models identified based on homology evidence to the domestic silk moth, Bombyx mori. The repeat content composes 41.18% of the genomic sequences which is in the range of other lepidopteran species. Conclusions: Functional annotation reveals that predicted gene models are involved in important cellular mechanisms such as metabolic pathways and protein synthesis. Thus, the data generated in this study expands our knowledge on the genomic characteristics of this devastating pest and provides essential resources for future genetic studies of the species.

Functional redundancy of structural proteins of the peritrophic membrane in Trichoplusia ni

The peritrophic membrane (or peritrophic matrix) (PM) in insects is formed by binding of PM proteins with multiple chitin binding domains (CBDs) to chitin fibrils. Multi-CBD chitin binding proteins (CBPs) and the insect intestinal mucin (IIM) are major PM structural proteins. To understand the biochemical and physiological role of IIM in structural formation and physiological function of the PM, Trichoplusia ni mutant strains lacking IIM were generated by CRISPR/Cas9 mutagenesis. The mutant T. ni larvae were confirmed to lack IIM, but PM formation was observed as in wild type larvae and lacking IIM in the PM did not result in changes of protease activities in the larval midgut. Larval growth and development of the mutant strains were similar to the wild type strain on artificial diet and cabbage leaves, but had a decreased survival in the 5th instar. The larvae of the mutant strains with the PM formed without IIM did not have a change of susceptibility to the infection of the baculovirus AcMNPV and the Bacillus thuringiensis (Bt) formulation Dipel, to the toxicity of the Bt toxins Cry1Ac and Cry2Ab and the chemical insecticide sodium aluminofluoride. Treatment of the mutant T. ni larvae with Calcofluor reduced the larval susceptibility to the toxicity of Bt Cry1Ac, as similarly observed in the wild type larvae. Overall, in the mutant T. ni larvae, the PM was formed without IIM and the lacking of IIM in the PM did not drastically impact the performance of larvae on diet or cabbage leaves under the laboratory conditions.

A haplotype-resolved, de novo genome assembly for the wood tiger moth (Arctia plantaginis) through trio binning

BACKGROUND

Diploid genome assembly is typically impeded by heterozygosity because it introduces errors when haplotypes are collapsed into a consensus sequence. Trio binning offers an innovative solution that exploits heterozygosity for assembly. Short, parental reads are used to assign parental origin to long reads from their F1 offspring before assembly, enabling complete haplotype resolution. Trio binning could therefore provide an effective strategy for assembling highly heterozygous genomes, which are traditionally problematic, such as insect genomes. This includes the wood tiger moth (Arctia plantaginis), which is an evolutionary study system for warning colour polymorphism.

FINDINGS

We produced a high-quality, haplotype-resolved assembly for Arctia plantaginis through trio binning. We sequenced a same-species family (F1 heterozygosity ∼1.9%) and used parental Illumina reads to bin 99.98% of offspring Pacific Biosciences reads by parental origin, before assembling each haplotype separately and scaffolding with 10X linked reads. Both assemblies are contiguous (mean scaffold N50: 8.2 Mb) and complete (mean BUSCO completeness: 97.3%), with annotations and 31 chromosomes identified through karyotyping. We used the assembly to analyse genome-wide population structure and relationships between 40 wild resequenced individuals from 5 populations across Europe, revealing the Georgian population as the most genetically differentiated with the lowest genetic diversity.

CONCLUSIONS

We present the first invertebrate genome to be assembled via trio binning. This assembly is one of the highest quality genomes available for Lepidoptera, supporting trio binning as a potent strategy for assembling heterozygous genomes. Using our assembly, we provide genomic insights into the geographic population structure of A. plantaginis.

Characterization of the First W-Specific Protein-Coding Gene for Sex Identification in Helicoverpa armigera

Helicoverpa armigera is a globally-important crop pest with a WZ (female)/ZZ (male) sex chromosome system. The absence of discernible sexual dimorphism in its egg and larval stages makes it impossible to address any sex-related theoretical and applied questions before pupation unless a W-specific sequence marker is available for sex diagnosis. To this end, we used one pair of morphologically pre-sexed pupae to PCR-screen 17 non-transposon transcripts selected from 4855 W-linked candidate reads identified by mapping a publicly available egg transcriptome of both sexes to the male genome of this species and detected the read SRR1015458.67499 only in the female pupa. Subsequent PCR screenings of this read and the previously reported female-specific RAPD (random amplified polymorphic DNA) marker AF18 with ten more pairs of pre-sexed pupae and different annealing positions and/or temperatures as well as its co-occurrence with the female-specific transcript splicing isoforms of doublesex gene of H. armigera (Hadsx) and amplification and sequencing of their 5′ unknown flanking sequences in three additional pairs of pre-sexed pupae verified that SRR1015458.67499 is a single copy protein-coding gene unique to W chromosome (named GUW1) while AF18 is a multicopy MITE transposon located on various chromosomes. Test application of GUW1 as a marker to sex 30 neonates of H. armigera yielded a female/male ratio of 1.14: 1.00. Both GUW1 and Hadsx splicing isoforms assays revealed that the H. armigera embryo cell line QB-Ha-E-1 is a male cell line. Taken together, GUW1 is not only a reliable DNA marker for sexing all stages of H. armigera and its cell lines, but also represents the first W-specific protein-coding gene in lepidopterans.

The genetic adaptations of fall armyworm Spodoptera frugiperda facilitated its rapid global dispersal and invasion

The fall armyworm (Spodoptera frugiperda) is a lepidopteran insect pest that causes huge economic losses. This notorious insect pest has rapidly spread over the world in the past few years. However, the mechanisms of rapid dispersal are not well understood. Here, we report a chromosome-level assembled genome of the fall armyworm, named the ZJ-version, using PacBio and Hi-C technology. The sequenced individual was a female collected from the Zhejiang province of China and had high heterozygosity. The assembled genome size of ZJ-version was 486 Mb, containing 361 contigs with an N50 of 1.13 Mb. Hi-C scaffolding further assembled the genome into 31 chromosomes and a portion of W chromosome, representing 97.4% of all contigs and resulted in a chromosome-level genome with scaffold N50 of 16.3 Mb. The sex chromosomes were identified by genome resequencing of a single male pupa and a single female pupa. About 28% of the genome was annotated as repeat sequences, and 22,623 protein-coding genes were identified. Comparative genomics revealed the expansion of the detoxification-associated gene families, chemoreception-associated gene families, nutrition metabolism and transport system gene families in the fall armyworm. Transcriptomic and phylogenetic analyses focused on these gene families revealed the potential roles of the genes in polyphagia and invasion of fall armyworm. The high-quality of the fall armyworm genome provides an important genomic resource for further explorations of the mechanisms of polyphagia and insecticide resistance, as well as for pest management of fall armyworm.

Stimulation of Insect Herbivory by Elevated Temperature Outweighs Protection by the Jasmonate Pathway

Rising global temperatures are associated with increases in the geographic range, population size, and feeding voracity of insect herbivores. Although it is well established that the plant hormone jasmonate (JA) promotes durable resistance to many ectothermic herbivores, little is known about how JA-mediated defense is influenced by rising temperatures. Here, we used the Arabidopsis-Trichoplusia ni (cabbage looper) interaction to investigate the relative contribution of JA and elevated temperature to host resistance. Video monitoring of T. ni larval behavior showed that elevated temperature greatly enhanced defoliation by increasing the bite rate and total time spent feeding, whereas loss of resistance in a JA-deficient mutant did not strongly affect these behaviors. The acceleration of insect feeding at elevated temperature was not attributed to decreases in wound-induced JA biosynthesis, expression of JA-responsive genes, or the accumulation of defensive glucosinolates prior to insect challenge. Quantitative proteomic analysis of insect frass, however, provided evidence for a temperature-dependent increase in the production of T. ni digestive enzymes. Our results demonstrate that temperature-driven stimulation of T. ni feeding outweighs the protective effects of JA-mediated resistance in Arabidopsis, thus highlighting a potential threat to plant resilience in a warming world.

Chromosomal-level reference genome of Chinese peacock butterfly (Papilio bianor) based on third-generation DNA sequencing and Hi-C analysis

BACKGROUND

Papilio bianor Cramer, 1777 (commonly known as the Chinese peacock butterfly) (Insecta, Lepidoptera, Papilionidae) is a widely distributed swallowtail butterfly with a wide number of geographic populations ranging from the southeast of Russia to China, Japan, India, Vietnam, Myanmar, and Thailand. Its wing color consists of both pigmentary colored scales (black, reddish) and structural colored scales (iridescent blue or green dust). A high-quality reference genome of P. bianor is an important foundation for investigating iridescent color evolution, phylogeography, and the evolution of swallowtail butterflies.

FINDINGS

We obtained a chromosome-level de novo genome assembly of the highly heterozygous P. bianor using long Pacific Biosciences sequencing reads and high-throughput chromosome conformation capture technology. The final assembly is 421.52 Mb on 30 chromosomes (29 autosomes and 1 Z sex chromosome) with 13.12 Mb scaffold N50. In total, 15,375 protein-coding genes and 233.09 Mb of repetitive sequences were identified. Phylogenetic analyses indicated that P. bianor separated from a common ancestor of swallowtails ∼23.69-36.04 million years ago. Demographic history suggested that the population expansion of this species from the last interglacial period to the last glacial maximum possibly resulted from its decreased natural enemies and its adaptation to climate change during the glacial period.

CONCLUSIONS

We present a high-quality chromosome-level reference genome of P. bianor using long-read single-molecule sequencing and Hi-C-based chromatin interaction maps. Our results lay the foundation for exploring the genetic basis of special biological features of P. bianor and also provide a useful data source for comparative genomics and phylogenomics among butterflies and moths.

A chromosome-level genome assembly of Cydia pomonella provides insights into chemical ecology and insecticide resistance

The codling moth Cydia pomonella, a major invasive pest of pome fruit, has spread around the globe in the last half century. We generated a chromosome-level scaffold assembly including the Z chromosome and a portion of the W chromosome. This assembly reveals the duplication of an olfactory receptor gene (OR3), which we demonstrate enhances the ability of C. pomonella to exploit kairomones and pheromones in locating both host plants and mates. Genome-wide association studies contrasting insecticide-resistant and susceptible strains identify hundreds of single nucleotide polymorphisms (SNPs) potentially associated with insecticide resistance, including three SNPs found in the promoter of CYP6B2. RNAi knockdown of CYP6B2 increases C. pomonella sensitivity to two insecticides, deltamethrin and azinphos methyl. The high-quality genome assembly of C. pomonella informs the genetic basis of its invasiveness, suggesting the codling moth has distinctive capabilities and adaptive potential that may explain its worldwide expansion.

The codling moth, Cydia pomonella, is one of the major pests of pome fruit (apples and pears) and walnuts. Here, the authors sequence and analyze its genome, providing insights on olfactory and detoxification processes that may underlie its worldwide expansion.

Mutation of ABC transporter ABCA2 confers resistance to Bt toxin Cry2Ab in Trichoplusia ni

Insecticidal proteins from Bacillus thuringiensis (Bt) are the primary recombinant proteins expressed in transgenic crops (Bt-crops) to confer insect resistance. Development of resistance to Bt toxins in insect populations threatens the sustainable application of Bt-crops in agriculture. The Bt toxin Cry2Ab is a major insecticidal protein used in current Bt-crops, and resistance to Cry2Ab has been selected in several insects, including the cabbage looper, Trichoplusia ni. In this study, the Cry2Ab resistance gene in T. ni was mapped to Chromosome 17 by genetic linkage analyses using a whole genome resequencing approach, and was then finely mapped using RNA-seq-based bulked segregant analysis (BSA) and amplicon sequencing (AmpSeq)-based fine linkage mapping to a locus containing two genes, ABCA1 and ABCA2. Mutations in ABCA1 and ABCA2 in Cry2Ab resistant T. ni were identified by both genomic DNA and cDNA sequencing. Analysis of the expression of ABCA1 and ABCA2 in T. ni larvae indicated that ABCA2 is abundantly expressed in the larval midgut, but ABCA1 is not a midgut-expressed gene. The mutation in ABCA2 in Cry2Ab resistant T. ni was identified to be an insertion of a transposon Tntransib in ABCA2. For confirmation of ABCA2 as the Cry2Ab-resistance gene, T. ni mutants with frameshift mutations in ABCA1 and ABCA2 were generated by CRISPR/Cas9 mutagenesis. Bioassays of the T. ni mutants with Cry2Ab verified that the mutations of ABCA1 did not change larval susceptibility to Cry2Ab, but the ABCA2 mutants were highly resistant to Cry2Ab. Genetic complementation test of the ABCA2 allele in Cry2Ab resistant T. ni with an ABCA2 mutant generated by CRISPR/Cas9 confirmed that the ABCA2 mutation in the Cry2Ab resistant strain confers the resistance. The results from this study confirmed that ABCA2 is essential for the toxicity of Cry2Ab in T. ni and mutation of ABCA2 confers the resistance to Cry2Ab in the resistant T. ni strain derived from a Bt resistant greenhouse population.

Genetic engineering approaches to improve posttranslational modification of biopharmaceuticals in different production platforms

The number of approved biopharmaceuticals, where product quality attributes remain of major importance, is increasing steadily. Within the available variety of expression hosts, the production of biopharmaceuticals faces diverse limitations with respect to posttranslational modifications (PTM), while different biopharmaceuticals demand different forms and specifications of PTMs for proper functionality. With the growing toolbox of genetic engineering technologies, it is now possible to address general as well as host- or biopharmaceutical-specific product quality obstacles. In this review, we present diverse expression systems derived from mammalians, bacteria, yeast, plants, and insects as well as available genetic engineering tools. We focus on genes for knockout/knockdown and overexpression for meaningful approaches to improve biopharmaceutical PTMs and discuss their applicability as well as future trends in the field.

Transcriptional Responses of the Trichoplusia ni Midgut to Oral Infection by the Baculovirus Autographa californica Multiple Nucleopolyhedrovirus

Baculoviruses are large double-stranded DNA viruses that are virulent pathogens of certain insect species. In a natural host, Trichoplusia ni, infection by the model baculovirus Autographa californica multiple nucleopolyhedrovirus (AcMNPV) begins when the occluded form of the virus disassembles in the midgut and virions infect midgut epithelial cells to establish the primary phase of the infection. To better understand the primary phase of the AcMNPV infection cycle, newly molted 5th-instar T. ni larvae were orally infected with AcMNPV occlusion bodies and the transcriptional responses of the T. ni midgut were analyzed at various times from 0 to 72 h postinfection, using transcriptome sequencing analysis and a T. ni reference genome. The numbers of differentially expressed host genes increased as the infection progressed, and we identified a total of 3,372 differentially expressed T. ni transcripts in the AcMNPV-infected midgut. Genes encoding orthologs of HMG176, atlastin, and CPH43 were among the most dramatically upregulated in response to AcMNPV infection. A number of cytochrome P450 genes were downregulated in response to infection. We also identified the effects of AcMNPV infection on a large variety of genes associated with innate immunity. This analysis provides an abundance of new and detailed information on host responses to baculovirus infection during the primary phase of the infection in the midgut and will be important for understanding how baculoviruses establish productive infections in the organism.IMPORTANCE Baculoviruses are virulent pathogens of a number of important insect pest species. In the host Trichoplusia ni, infection begins in the midgut when infectious virions of the occlusion-derived virus (ODV) phenotype enter and subsequently replicate in cells of the midgut epithelium. A second virion phenotype (budded virus [BV]) is produced there, and BV mediates systemic infection of the animal. Most prior detailed studies of baculovirus infections have focused on BV infections of cultured cells. In this study, we examined the transcriptional responses of the T. ni midgut to infection by ODV of the baculovirus AcMNPV and identified a variety of host genes that respond dramatically to viral infection. Understanding the transcriptional responses of the host midgut to viral infection is critically important for understanding the biphasic infection in the animal as a whole.