De novo genome assembly of the tobacco hornworm moth (Manduca sexta)

A long-read draft assembly of the Chinese mantis (Mantodea: Mantidae: Tenodera sinensis) genome reveals patterns of ion channel gain and loss across Arthropoda

Praying mantids (Mantodea: Mantidae) are iconic insects that have captivated biologists for decades, especially the species with cannibalistic copulatory behavior. This behavior has been cited as evidence that insects lack nociceptive capacities and cannot feel pain; however, this behaviorally driven hypothesis has never been rigorously tested at the genetic or functional level. To enable future studies of nociceptive capabilities in mantids, we sequenced and assembled a draft genome of the Chinese praying mantis (Tenodera sinensis) and identified multiple classes of nociceptive ion channels by comparison to orthologous gene families in Arthropoda. Our assembly-produced using PacBio HiFi reads-is fragmented (total size = 3.03 Gb; N50 = 1.8 Mb; 4,966 contigs), but is highly complete with respect to gene content (BUSCO complete = 98.7% [odb10_insecta]). The size of our assembly is substantially larger than that of most other insects, but is consistent with the size of other mantid genomes. We found that most families of nociceptive ion channels are present in the T. sinensis genome; that they are most closely related to those found in the damp-wood termite (Zootermopsis nevadensis); and that some families have expanded in T. sinensis while others have contracted relative to nearby lineages. Our findings suggest that mantids are likely to possess nociceptive capabilities and provide a foundation for future experimentation regarding ion channel functions and their consequences for insect behavior.

Starvation-Induced Changes to the Midgut Proteome and Neuropeptides in Manduca sexta

Starvation is a complex physiological state that induces changes in protein expression to ensure survival. The insect midgut is sensitive to changes in dietary content as it is at the forefront of communicating information about incoming nutrients to the body via hormones. Therefore, a DIA proteomics approach was used to examine starvation physiology and, specifically, the role of midgut neuropeptide hormones in a representative lepidopteran, Manduca sexta. Proteomes were generated from midguts of M. sexta fourth-instar caterpillars, starved for 24 h and 48 h, and compared to fed controls. A total of 3047 proteins were identified, and 854 of these were significantly different in abundance. KEGG analysis revealed that metabolism pathways were less abundant in starved caterpillars, but oxidative phosphorylation proteins were more abundant. In addition, six neuropeptides or related signaling cascade proteins were detected. Particularly, neuropeptide F1 (NPF1) was significantly higher in abundance in starved larvae. A change in juvenile hormone-degrading enzymes was also detected during starvation. Overall, our results provide an exploration of the midgut response to starvation in M. sexta and validate DIA proteomics as a useful tool for quantifying insect midgut neuropeptide hormones.

An enteric ultrastructural surface atlas of the model insect Manducasexta

The tobacco hornworm is a laboratory model that is particularly suitable for analyzing gut inflammation, but a physiological reference standard is currently unavailable. Here, we present a surface atlas of the healthy hornworm gut generated by scanning electron microscopy and nano-computed tomography. This comprehensive overview of the gut surface reveals morphological differences between the anterior, middle, and posterior midgut, allowing the screening of aberrant gut phenotypes while accommodating normal physiological variations. We estimated a total resorptive midgut surface of 0.42 m2 for L5d6 larvae, revealing its remarkable size. Our data will support allometric scaling and dose conversion from Manduca sexta to mammals in preclinical research, embracing the 3R principles. We also observed non-uniform gut colonization by enterococci, characterized by dense biofilms in the pyloric cone and downstream of the pylorus associated with pore and spine structures in the hindgut intima, indicating a putative immunosurveillance function in the lepidopteran hindgut.

High-quality haploid genomes corroborate 29 chromosomes and highly conserved synteny of genes in Hyles hawkmoths (Lepidoptera: Sphingidae)

BACKGROUND

Morphological and traditional genetic studies of the young Pliocene genus Hyles have led to the understanding that despite its importance for taxonomy, phenotypic similarity of wing patterns does not correlate with phylogenetic relationship. To gain insights into various aspects of speciation in the Spurge Hawkmoth (Hyles euphorbiae), we assembled a chromosome-level genome and investigated some of its characteristics.

RESULTS

The genome of a male H. euphorbiae was sequenced using PacBio and Hi-C data, yielding a 504 Mb assembly (scaffold N50 of 18.2 Mb) with 99.9% of data represented by the 29 largest scaffolds forming the haploid chromosome set. Consistent with this, FISH analysis of the karyotype revealed n = 29 chromosomes and a WZ/ZZ (female/male) sex chromosome system. Estimates of chromosome length based on the karyotype image provided an additional quality metric of assembled chromosome size. Rescaffolding the published male H. vespertilio genome resulted in a high-quality assembly (651 Mb, scaffold N50 of 22 Mb) with 98% of sequence data in the 29 chromosomes. The larger genome size of H. vespertilio (average 1C DNA value of 562 Mb) was accompanied by a proportional increase in repeats from 45% in H. euphorbiae (measured as 472 Mb) to almost 55% in H. vespertilio. Several wing pattern genes were found on the same chromosomes in the two species, with varying amounts and positions of repetitive elements and inversions possibly corrupting their function.

CONCLUSIONS

Our two-fold comparative genomics approach revealed high gene synteny of the Hyles genomes to other Sphingidae and high correspondence to intact Merian elements, the ancestral linkage groups of Lepidoptera, with the exception of three simple fusion events. We propose a standardized approach for genome taxonomy using nucleotide homology via scaffold chaining as the primary tool combined with Oxford plots based on Merian elements to infer and visualize directionality of chromosomal rearrangements. The identification of wing pattern genes promises future understanding of the evolution of forewing patterns in the genus Hyles, although further sequencing data from more individuals are needed. The genomic data obtained provide additional reliable references for further comparative studies in hawkmoths (Sphingidae).

A salivary GMC oxidoreductase of Manduca sexta re-arranges the green leaf volatile profile of its host plant

Green leaf volatiles (GLVs) are short-chain oxylipins that are emitted from plants in response to stress. Previous studies have shown that oral secretions (OS) of the tobacco hornworm Manduca sexta, introduced into plant wounds during feeding, catalyze the re-arrangement of GLVs from Z-3- to E-2-isomers. This change in the volatile signal however is bittersweet for the insect as it can be used by their natural enemies, as a prey location cue. Here we show that (3Z):(2E)-hexenal isomerase (Hi-1) in M. sexta's OS catalyzes the conversion of the GLV Z-3-hexenal to E-2-hexenal. Hi-1 mutants that were raised on a GLV-free diet showed developmental disorders, indicating that Hi-1 also metabolizes other substrates important for the insect's development. Phylogenetic analysis placed Hi-1 within the GMCβ-subfamily and showed that Hi-1 homologs from other lepidopterans could catalyze similar reactions. Our results indicate that Hi-1 not only modulates the plant's GLV-bouquet but also functions in insect development.

A quantitative micro-tomographic gut atlas of the lepidopteran model insect Manduca sexta

The tobacco hornworm is used extensively as a model system for ecotoxicology, immunology and gut physiology. Here, we established a micro-computed tomography approach based on the oral application of the clinical contrast agent iodixanol, allowing for a high-resolution quantitative analysis of the Manduca sexta gut. This technique permitted the identification of previously unknown and understudied structures, such as the crop or gastric ceca, and revealed the underlying complexity of the hindgut folding pattern, which is involved in fecal pellet formation. The acquired data enabled the volume rendering of all gut parts, the reliable calculation of their volumes, and the virtual endoscopy of the entire alimentary tract. It can provide information for accurate orientation in histology uses, enable quantitative anatomical phenotyping in three dimensions, and allow the calculation of locally effective midgut concentrations of applied chemicals. This atlas will provide critical insights into the evolution of the alimentary tract in lepidopterans.

A high-quality, long-read genome assembly of the whitelined sphinx moth (Lepidoptera: Sphingidae: Hyles lineata) shows highly conserved melanin synthesis pathway genes

The sphinx moth genus Hyles comprises 29 described species inhabiting all continents except Antarctica. The genus diverged relatively recently (40-25 MYA), arising in the Americas and rapidly establishing a cosmopolitan distribution. The whitelined sphinx moth, Hyles lineata, represents the oldest extant lineage of this group and is one of the most widespread and abundant sphinx moths in North America. Hyles lineata exhibits the large body size and adept flight control characteristic of the sphinx moth family (Sphingidae), but it is unique in displaying extreme larval color variation and broad host plant use. These traits, in combination with its broad distribution and high relative abundance within its range, have made H. lineata a model organism for studying phenotypic plasticity, plant-herbivore interactions, physiological ecology, and flight control. Despite being one of the most well-studied sphinx moths, little data exist on genetic variation or regulation of gene expression. Here, we report a high-quality genome showing high contiguity (N50 of 14.2 Mb) and completeness (98.2% of Lepidoptera BUSCO genes), an important first characterization to facilitate such studies. We also annotate the core melanin synthesis pathway genes and confirm that they have high sequence conservation with other moths and are most similar to those of another, well-characterized sphinx moth, the tobacco hornworm (Manduca sexta).

Genome-wide identification, classification, and expression profiling of serine esterases and other esterase-related proteins in the tobacco hornworm, Manduca sexta

Serine esterases (SEs) are hydrolases that catalyze the conversion of carboxylic esters into acids and alcohols. Lipases and carboxylesterases constitute two major groups of SEs. Although over a hundred of insect genomes are known, systematic identification and classification of SEs are rarely performed, likely due to large size and complex composition of the gene family in each species. Considering their key roles in lipid metabolism and other physiological processes, we have categorized 144 M. sexta SEs and SE homologs (SEHs), 114 of which contain a motif of GXSXG. Multiple sequence alignment and phylogenetic tree analysis have revealed 39 neutral lipases (NLs), 3 neutral lipase homologs (NLHs), 11 acidic lipases (ALs), 3 acidic lipase homologs (ALHs), a lipase-3, a triglyceride lipase, a monoglyceride lipase, a hormone-sensitive lipase, and a GDSL lipase. Eighty-three carboxylesterase genes encode 29 α-esterases (AEs), 12 AEHs (e.g., SEH4-1-3), 20 feruloyl esterases (FEs), 2 FEHs, 2 β-esterases (BEs), 2 integument esterases (IEs), 1 IEH, 4 juvenile hormone esterases, 2 acetylcholinesterases, gliotactin, 6 neuroligins, neurotactin, and an uncharacteristic esterase homolog. In addition to these GXSXG proteins, we have identified 26 phospholipases and 13 thioesterases. Expression profiling of these genes in specific tissues and stages has provided insights into their functions including digestion, detoxification, hormone processing, neurotransmission, reproduction, and developmental regulation. In summary, we have established a framework of information on SEs and related proteins in M. sexta to stimulate their research in the model species and comparative investigations in agricultural pests or disease vectors.

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.

Parallel evolution of opsin visual pigments in hawkmoths by tuning of spectral sensitivities during transition from a nocturnal to a diurnal ecology

Light environments differ dramatically between day and night. The transition between diurnal and nocturnal visual ecology has happened repeatedly throughout evolution in many species. However, the molecular mechanism underlying the evolution of vision in recent diurnal-nocturnal transition is poorly understood. Here, we focus on hawkmoths (Lepidoptera: Sphingidae) to address this question by investigating five nocturnal and five diurnal species. We performed RNA-sequencing analysis and identified opsin genes corresponding to the ultraviolet (UV), short-wavelength (SW) and long-wavelength (LW)-absorbing visual pigments. We found no significant differences in the expression patterns of opsin genes between the nocturnal and diurnal species. We then constructed the phylogenetic trees of hawkmoth species and opsins. The diurnal lineages had emerged at least three times from the nocturnal ancestors. The evolutionary rates of amino acid substitutions in the three opsins differed between the nocturnal and diurnal species. We found an excess number of parallel amino acid substitutions in the opsins in three independent diurnal lineages. The numbers were significantly more than those inferred from neutral evolution, suggesting that positive selection acted on these parallel substitutions. Moreover, we predicted the visual pigment absorption spectra based on electrophysiologically determined spectral sensitivity in two nocturnal and two diurnal species belonging to different clades. In the diurnal species, the LW pigments shift 10 nm towards shorter wavelengths, and the SW pigments shift 10 nm in the opposite direction. Taken together, our results suggest that parallel evolution of opsins may have enhanced the colour discrimination properties of diurnal hawkmoths in ambient light.

Genome sequence and silkomics of the spindle ermine moth, Yponomeuta cagnagella, representing the early diverging lineage of the ditrysian Lepidoptera

Many lepidopteran species produce silk, cocoons, feeding tubes, or nests for protection from predators and parasites for caterpillars and pupae. Yet, the number of lepidopteran species whose silk composition has been studied in detail is very small, because the genes encoding the major structural silk proteins tend to be large and repetitive, making their assembly and sequence analysis difficult. Here we have analyzed the silk of Yponomeuta cagnagella, which represents one of the early diverging lineages of the ditrysian Lepidoptera thus improving the coverage of the order. To obtain a comprehensive list of the Y. cagnagella silk genes, we sequenced and assembled a draft genome using Oxford Nanopore and Illumina technologies. We used a silk-gland transcriptome and a silk proteome to identify major silk components and verified the tissue specificity of expression of individual genes. A detailed annotation of the major genes and their putative products, including their complete sequences and exon-intron structures is provided. The morphology of silk glands and fibers are also shown. This study fills an important gap in our growing understanding of the structure, evolution, and function of silk genes and provides genomic resources for future studies of the chemical ecology of Yponomeuta species.

Insect Models in Nutrition Research

Insects are the most diverse organisms on earth, accounting for ~80% of all animals. They are valuable as model organisms, particularly in the context of genetics, development, behavior, neurobiology and evolutionary biology. Compared to other laboratory animals, insects are advantageous because they are inexpensive to house and breed in large numbers, making them suitable for high-throughput testing. They also have a short life cycle, facilitating the analysis of generational effects, and they fulfil the 3R principle (replacement, reduction and refinement). Many insect genomes have now been sequenced, highlighting their genetic and physiological similarities with humans. These factors also make insects favorable as whole-animal high-throughput models in nutritional research. In this review, we discuss the impact of insect models in nutritional science, focusing on studies investigating the role of nutrition in metabolic diseases and aging/longevity. We also consider food toxicology and the use of insects to study the gut microbiome. The benefits of insects as models to study the relationship between nutrition and biological markers of fitness and longevity can be exploited to improve human health.

Sex- and tissue-specific expression of chemosensory receptor genes in a hawkmoth

For the nocturnal hawkmoth Manduca sexta, olfactory and gustatory cues are essential for finding partners, food, and oviposition sites. Three chemosensory receptor families, odorant receptors (ORs), ionotropic receptors (IRs), and gustatory receptors (GRs) are involved in the detection of these stimuli. While many chemosensory receptor genes have been identified, knowledge of their expression profile in potentially chemoreceptive organs is incomplete. Here, we studied the expression patterns of chemosensory receptors in different tissues including the antennae, labial palps, proboscis, legs, wings and ovipositor. We compared the receptors’ expression in female and male moths both before and after mating by using the NanoString platform. This tool allowed us to measure expression levels of chemosensory receptor genes in a single reaction using probes designed against 71 OR, 29 IR and 49 GR transcripts. In all tissues investigated, we detected expression of genes from all three receptor families. The highest number of receptors was detected in the antennae (92), followed by the ovipositor (59), while the least number was detected in the hindlegs (21). The highest number of OR genes were expressed in the antennae (63), of which 24 were specific to this main olfactory organ. The highest number of IRs were also expressed in the antennae (16), followed by the ovipositor (15). Likewise, antennae and ovipositor expressed the highest number of GRs (13 and 14). Expression of the OR co-receptor MsexORCo, presumably a prerequisite for OR function, was found in the antennae, labial palps, forelegs and ovipositor. IR co-receptors MsexIR25a and MsexIR76b were expressed across all tested tissues, while expression of the IR co-receptor MsexIR8a was restricted to antennae and ovipositor. Comparing the levels of all 149 transcripts across the nine tested tissues allowed us to identify sex-biased gene expression in the antennae and the legs, two appendages that are also morphologically different between the sexes. However, none of the chemosensory receptors was differentially expressed based on the moths’ mating state. The observed gene expression patterns form a strong base for the functional characterization of chemosensory receptors and the understanding of olfaction and gustation at the molecular level in M. sexta.

Population differentiation and structural variation in the Manduca sexta genome across the United States

Many species that are extensively studied in the laboratory are less well characterized in their natural habitat, and laboratory strains represent only a small fraction of the variation in a species’ genome. Here we investigate genomic variation in 3 natural North American populations of an agricultural pest and a model insect for many scientific disciplines, the tobacco hornworm (Manduca sexta). We show that hornworms from Arizona, Kansas, and North Carolina are genetically distinct, with Arizona being particularly differentiated from the other 2 populations using Illumina whole-genome resequencing. Peaks of differentiation exist across the genome, but here, we focus in on the most striking regions. In particular, we identify 2 likely segregating inversions found in the Arizona population. One inversion on the Z chromosome may enhance adaptive evolution of the sex chromosome. The larger, 8 Mb inversion on chromosome 12 contains a pseudogene which may be involved in the exploitation of a novel hostplant in Arizona, but functional genetic assays will be required to support this hypothesis. Nevertheless, our results reveal undiscovered natural variation and provide useful genomic data for both pest management and evolutionary genetics of this insect species.

Dynamic monitoring of vital functions and tissue re-organization in Saturnia pavonia (Lepidoptera, Saturniidae) during final metamorphosis by non-invasive MRI

Magnetic resonance imaging (MRI) is the key whole-body imaging technology for observing processes within a living object providing excellent resolution and contrast between soft tissues. In the present work, we exploited the non-destructive properties of MRI to track longitudinally the dynamic changes that take place in developing pupae of the Emperor Moth (Saturnia pavonia) during the last days before eclosion. While in diapause pupae, body fluid was almost homogeneously distributed over the internal compartments, as soon as wings, legs, flight muscles and the head region were fully developed, a significant redistribution of water levels occurred between thoracic and abdominal regions. During the last two days before eclosion, the developing moths transferred substantial amounts of liquid into the gut and the labial gland, and in case of females, into developing eggs. Concomitantly, the volume of the air sacs increased drastically and their expansion/compression became clearly visible in time-resolved MR images. Furthermore, besides ventilation of the tracheal system, air sacs are likely to serve as volume reservoir for liquid transfer during development of the moths inside their pupal case. In parallel, we were able to monitor noninvasively lipid consumption, cardiac activity and haemolymph circulation during final metamorphosis.

Comparative Analysis of PacBio and Oxford Nanopore Sequencing Technologies for Transcriptomic Landscape Identification of Penaeus monodon

With the advantages that long-read sequencing platforms such as Pacific Biosciences (Menlo Park, CA, USA) (PacBio) and Oxford Nanopore Technologies (Oxford, UK) (ONT) can offer, various research fields such as genomics and transcriptomics can exploit their benefits. Selecting an appropriate sequencing platform is undoubtedly crucial for the success of the research outcome, thus there is a need to compare these long-read sequencing platforms and evaluate them for specific research questions. This study aims to compare the performance of PacBio and ONT platforms for transcriptomic analysis by utilizing transcriptome data from three different tissues (hepatopancreas, intestine, and gonads) of the juvenile black tiger shrimp, Penaeus monodon. We compared three important features: (i) main characteristics of the sequencing libraries and their alignment with the reference genome, (ii) transcript assembly features and isoform identification, and (iii) correlation of the quantification of gene expression levels for both platforms. Our analyses suggest that read-length bias and differences in sequencing throughput are highly influential factors when using long reads in transcriptome studies. These comparisons can provide a guideline when designing a transcriptome study utilizing these two long-read sequencing technologies.

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.