Mechanisms of macroevolution: polyphagous plasticity in butterfly larvae revealed by RNA-Seq

Evolved and Plastic Gene Expression in Adaptation of a Specialist Fly to a Novel Niche

How gene expression evolves to enable divergent ecological adaptation and how changes in gene expression relate to genomic architecture are pressing questions for understanding the mechanisms enabling adaptation and ecological speciation. Furthermore, how plasticity in gene expression can both contribute to and be affected by the process of ecological adaptation is crucial to understanding gene expression evolution, colonisation of novel niches and response to rapid environmental change. Here, we investigate the role of constitutive and plastic gene expression differences between host races, or host-specific ecotypes, of the peacock fly Tephritis conura, a thistle bud specialist. By cross-fostering larvae to new buds of their natal host plant or the alternative, novel host plant, we uncover extensive constitutive differences in gene expression between the host races, especially genes associated with processing of host plant chemicals. However, evidence for expression plasticity was minimal and limited to the ancestral host race. Genes with host race-specific expression are found more often than expected within a large inversion in the T. conura genome, adding to evidence that inversions are important for enabling diversification in the face of gene flow and underscores that altered gene expression may be key to understanding the evolutionary consequences of inversions.

Chromosome-level genome assembly of cotton thrips Thrips tabaci (Thysanoptera: Thripidae)

Cotton thrip, Thrips tabaci is a major polyphagous pest widely distributed on a variety of crops around the world, causing huge economic losses to agricultural production. Due to its biological and genomic characteristics, this pest can reproduce quickly and develop resistance to various pesticides in a very short time. However, the lack of high-quality reference genomes has hindered deeper gene function exploration and slows down the development of new management strategies. Here, we assembled a high-quality genome of T. tabaci at the chromosome level for the first time by using Illumina, PacBio long reads, and Hi-C technologies. The 329.59 Mb genome was obtained from 320 contigs, with a contig N50 of 1.53 Mb, and 94.21% of the assembly was anchored to 18 chromosomes. In total, 17,816 protein-coding genes were annotated, and 96.78% of BUSCO genes were fully represented. In conclusion, this high-quality genome provides a valuable genetic basis for our understanding of the biology of T. tabaci and contributes to the development of management strategies for cotton thrip.

Plasticity for the win: Flexible transcriptional response to host plant switches in the comma butterfly (Polygonia c-album)

Generalist plant-feeding insects are characterised by a broad host repertoire that can comprise several families or even different orders of plants. The genetic and physiological mechanisms underlying the use of such a wide host range are still not fully understood. Earlier studies indicate that the consumption of different host plants is associated with host-specific gene expression profiles. It remained, however, unclear if and how larvae can alter these profiles in the case of a changing host environment. Using the polyphagous comma butterfly (Polygonia c-album) we show that larvae can adjust their transcriptional profiles in response to a new host plant. The switch to some of the host plants, however, resulted in a larger transcriptional response and, thus, seems to be more challenging. At a physiological level, no correspondence for these patterns could be found in larval performance. This suggests that a high transcriptional but also phenotypic flexibility are essential for the use of a broad and diverse host range. We furthermore propose that host switch tests in the laboratory followed by transcriptomic investigations can be a valuable tool to examine not only plasticity in host use but also subtle and/or transient trade-offs in the evolution of host plant repertoires.

Transcriptional signature of host shift in the seed beetle Zabrotes subfasciatus

In phytophagous insects, adaptation to a new host is a dynamic process, in which early and later steps may be underpinned by different features of the insect genome. Here, we tested the hypothesis that early steps of this process are underpinned by a shift in gene expression patterns. We set up a short-term artificial selection experiment (10 generations) for the use of an alternative host (Cicer arietinum) on populations of the bean beetle Zabrotes subfasciatus. Using Illumina sequencing on young adult females, we show the selected populations differ in the expression of genes associated to stimuli, signalling, and developmental processes. Particularly, the "C. arietinum" population shows upregulation of histone methylation genes, which may constitute a strategy for fine-tuning the insect global gene expression network. Using qPCR on body regions, we demonstrated that the "Phaseolus vulgaris" population upregulates the genes polygalacturonase and egalitarian and that the expression of an odorant receptor transcript variant changes over generations. Moreover, in this population we detected the existence of vitellogenin (Vg) variants in both males and females, possibly harbouring canonical reproductive function in females and extracellular unknown functions in males. This study provides the basis for future genomic investigations seeking to shed light on the nature of the proximate mechanisms involved in promoting differential gene expression associated to insect development and adaptation to new hosts.

Gene expression plasticity facilitates different host feeding in Ips sexdentatus (Coleoptera: Curculionidae: Scolytinae)

Host shift is ecologically advantageous and a crucial driver for herbivore insect speciation. Insects on the non-native host obtain enemy-free space and confront reduced competition, but they must adapt to survive. Such signatures of adaptations can often be detected at the gene expression level. It is astonishing how bark beetles cope with distinct chemical environments while feeding on various conifers. Hence, we aim to disentangle the six-toothed bark beetle (Ips sexdentatus) response against two different conifer defences upon host shift (Scots pine to Norway spruce). We conducted bioassay and metabolomic analysis followed by RNA-seq experiments to comprehend the beetle's ability to surpass two different terpene-based conifer defence systems. Beetle growth rate and fecundity were increased when reared exclusively on spruce logs (alternative host) compared to pine logs (native host). Comparative gene expression analysis identified differentially expressed genes (DEGs) related to digestion, detoxification, transporter activity, growth, signalling, and stress response in the spruce-feeding beetle gut. Transporter genes were highly abundant during spruce feeding, suggesting they could play a role in pumping a wide variety of endogenous and xenobiotic compounds or allelochemicals out. Trehalose transporter (TRET) is also up-regulated in the spruce-fed beetle gut to maintain homeostasis and stress tolerance. RT-qPCR and enzymatic assays further corroborated some of our findings. Taken together, the transcriptional plasticity of key physiological genes plays a crucial role after the host shift and provides vital clues for the adaptive potential of bark beetles on different conifer hosts.

Chromosome level genome assembly of oriental armyworm Mythimna separata

The oriental armyworm, Mythimna separata, is an extremely destructive polyphagous pest with a broad host range that seriously threatens the safety of agricultural production. Here, a high-quality chromosome-level genome was assembled using Illumina, PacBio HiFi long sequencing, and Hi-C scaffolding technologies. The genome size was 706.30 Mb with a contig N50 of 22.08 Mb, and 99.2% of the assembled sequences were anchored to 31 chromosomes. In addition, 20,375 protein-coding genes and 258.68 Mb transposable elements were identified. The chromosome-level genome assembly of M. separata provides a significant genetic resource for future studies of this insect and contributes to the development of management strategies.

Midgut serine proteinases participate in dietary adaptations of the castor (Eri) silkworm Samia ricini Anderson transferred from Ricinus communis to an ancestral host, Ailanthus excelsa Roxb

Dietary change influenced the life-history traits, nutritional utilization, and midgut serine proteinases in the larvae of the domesticated polyphagous S. ricini, transferred from R. communis (common name: castor; family Euphorbiaceae; the host plant implicated in its domestication) to A. excelsa (common name: Indian tree of heaven; family Simaroubaceae; an ancestral host of wild Samia species). Significantly higher values for fecundity and body weight were observed in larvae feeding on R. communis (Scr diet), and they took less time to reach pupation than insects feeding on A. excelsa (Scai diet). Nevertheless, the nutritional index for efficiency of conversion of digested matter (ECD) was similar for larvae feeding on the two plant species, suggesting the physiological adaptation of S. ricini (especially older instars) to an A. excelsa diet. In vitro protease assays and gelatinolytic zymograms using diagnostic substrates and protease inhibitors revealed significantly elevated levels (p ≤ 0.05) of digestive trypsins, which may be associated with the metabolic costs influencing slow growth in larvae feeding on A. excelsa. RT-PCR with semidegenerate serine proteinase gene-specific primers, and cloning and sequencing of 3' cDNA ends identified a large gene family comprising at least two groups of putative chymotrypsins (i.e., Sr I and Sr II) resembling invertebrate brachyurins/collagenases with wide substrate specificities, and five groups of putative trypsins (i.e., Sr III, Sr IV, Sr V, Sr VII, and Sr VIII). Quantitative RT-PCR indicated that transcripts belonging to the Sr I, Sr III, Sr IV, and Sr V groups, especially the Sr IV group (resembling achelase I from Lonomia achelous), were expressed differentially in the midguts of fourth instars reared on the two plant species. Sequence similarity indicated shared lineages with lepidopteran orthologs associated with expression in the gut, protein digestion, and phytophagy. The results obtained are discussed in the context of larval serine proteinases in dietary adaptations, domestication, and exploration of new host plant species for commercial rearing of S. ricini.

Metabolic consequences of various fruit-based diets in a generalist insect species

Most phytophagous insect species exhibit a limited diet breadth and specialize on a few or a single host plant. In contrast, some species display a remarkably large diet breadth, with host plants spanning several families and many species. It is unclear, however, whether this phylogenetic generalism is supported by a generic metabolic use of common host chemical compounds ('metabolic generalism') or alternatively by distinct uses of diet-specific compounds ('multi-host metabolic specialism')? Here, we simultaneously investigated the metabolomes of fruit diets and of individuals of a generalist phytophagous species, Drosophila suzukii, that developed on them. The direct comparison of metabolomes of diets and consumers enabled us to disentangle the metabolic fate of common and rarer dietary compounds. We showed that the consumption of biochemically dissimilar diets resulted in a canalized, generic response from generalist individuals, consistent with the metabolic generalism hypothesis. We also showed that many diet-specific metabolites, such as those related to the particular color, odor, or taste of diets, were not metabolized, and rather accumulated in consumer individuals, even when probably detrimental to fitness. As a result, while individuals were mostly similar across diets, the detection of their particular diet was straightforward. Our study thus supports the view that dietary generalism may emerge from a passive, opportunistic use of various resources, contrary to more widespread views of an active role of adaptation in this process. Such a passive stance towards dietary chemicals, probably costly in the short term, might favor the later evolution of new diet specializations.

A draft Diabrotica virgifera virgifera genome: insights into control and host plant adaption by a major maize pest insect

BACKGROUND

Adaptations by arthropod pests to host plant defenses of crops determine their impacts on agricultural production. The larval host range of western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is restricted to maize and a few grasses. Resistance of D. v. virgifera to crop rotation practices and multiple insecticides contributes to its status as the most damaging pest of cultivated maize in North America and Europe. The extent to which adaptations by this pest contributes to host plant specialization remains unknown.

RESULTS

A 2.42 Gb draft D. v. virgifera genome, Dvir_v2.0, was assembled from short shotgun reads and scaffolded using long-insert mate-pair, transcriptome and linked read data. K-mer analysis predicted a repeat content of ≥ 61.5%. Ortholog assignments for Dvir_2.0 RefSeq models predict a greater number of species-specific gene duplications, including expansions in ATP binding cassette transporter and chemosensory gene families, than in other Coleoptera. A majority of annotated D. v. virgifera cytochrome P450s belong to CYP4, 6, and 9 clades. A total of 5,404 transcripts were differentially-expressed between D. v. virgifera larvae fed maize roots compared to alternative host (Miscanthus), a marginal host (Panicum virgatum), a poor host (Sorghum bicolor) and starvation treatments; Among differentially-expressed transcripts, 1,908 were shared across treatments and the least number were between Miscanthus compared to maize. Differentially-expressed transcripts were enriched for putative spliceosome, proteosome, and intracellular transport functions. General stress pathway functions were unique and enriched among up-regulated transcripts in marginal host, poor host, and starvation responses compared to responses on primary (maize) and alternate hosts.

CONCLUSIONS

Manual annotation of D. v. virgifera Dvir_2.0 RefSeq models predicted expansion of paralogs with gene families putatively involved in insecticide resistance and chemosensory perception. Our study also suggests that adaptations of D. v. virgifera larvae to feeding on an alternate host plant invoke fewer transcriptional changes compared to marginal or poor hosts. The shared up-regulation of stress response pathways between marginal host and poor host, and starvation treatments may reflect nutrient deprivation. This study provides insight into transcriptomic responses of larval feeding on different host plants and resources for genomic research on this economically significant pest of maize.

Adaptation of Helicoverpa armigera to Soybean Peptidase Inhibitors Is Associated with the Transgenerational Upregulation of Serine Peptidases

Molecular phenotypes induced by environmental stimuli can be transmitted to offspring through epigenetic inheritance. Using transcriptome profiling, we show that the adaptation of Helicoverpa armigera larvae to soybean peptidase inhibitors (SPIs) is associated with large-scale gene expression changes including the upregulation of genes encoding serine peptidases in the digestive system. Furthermore, approximately 60% of the gene expression changes induced by SPIs persisted in the next generation of larvae fed on SPI-free diets including genes encoding regulatory, oxidoreductase, and protease functions. To investigate the role of epigenetic mechanisms in regulating SPI adaptation, the methylome of the digestive system of first-generation larvae (fed on a diet with and without SPIs) and of the progeny of larvae exposed to SPIs were characterized. A comparative analysis between RNA-seq and Methyl-seq data did not show a direct relationship between differentially methylated and differentially expressed genes, while trypsin and chymotrypsin genes were unmethylated in all treatments. Rather, DNA methylation potential epialleles were associated with transcriptional and translational controls; these may play a regulatory role in the adaptation of H. armigera to SPIs. Altogether, our findings provided insight into the mechanisms of insect adaptation to plant antiherbivore defense proteins and illustrated how large-scale transcriptional reprograming of insect genes can be transmitted across generations.

Sublethal toxicity, transgenerational effects, and transcriptome expression of the neonicotinoid pesticide cycloxaprid on demographic fitness of Coccinella septempunctata

Evaluating side effects of new neonicotinoids in terms of sublethal doses and transcriptome expression is a crucial but challenging part of integrated pest management (IPM) approaches. To this end, a study of lethal and sublethal effects on Coccinella septempunctata larvae was conducted, and an age-stage, two-sex life table procedure was performed to investigate life-table parameters. Cycloxaprid (CYC) was shown to have adverse effects on survival, development, total longevity, reproductive capacity, and predation ability in C. septempunctata. In addition, demographic growth parameters of the F1 generation such as net reproductive rate, and the intrinsic and finite rates of increase were significantly decreased under sublethal dosage LR₃₀ (1.91 g ai/hm²). These results demonstrated that the population growth of C. septempunctata was impacted by a sublethal dosage of CYC. For transcriptome expression, 544 up- and 338 down-regulated significantly differentially expressed genes (DEGs), were observed between LR₃₀ treatment and control groups. Moreover, pathways related to metabolism of retinol, carcinogenesis, biosynthesis of steroid hormone, P450 metabolism, and metabolism of xenobiotics were identified in KEGG pathway analysis. Ten DEGs were chosen and confirmed with quantitative real-time PCR analysis. Based on these findings, CYC should be considered as a component of IPM strategies in the field.

De Novo Transcription Responses Describe Host-Related Differentiation of Paracoccus marginatus (Hemiptera: Pseudococcidae)

Paracoccus marginatus (Hemiptera: Pseudococcidae) is an invasive pest with a diverse host range, strong diffusion, and high fecundity. It has been observed that P. marginatus feeding on Carica papaya have a higher survival rate, fecundity, and longer lifespan than P. marginatus feeding on Solanum tuberosum, indicating their successful adaptation to C. papaya; however, the mechanisms underlying host plant adaptation remain unclear. Therefore, RNA-seq was performed to study the transcriptional responses of P. marginatus feeding on C. papaya and S. tuberosum plants. A total of 408 genes with significant differential expression were defined; most of them were downregulated in S. tuberosum, including those of digestive enzymes, detoxifying enzymes, ribosomes, and reproductive-related genes, which may result from the adaptation of the host to nutritional needs and changes in toxic chemical levels. Enrichment analysis of the Kyoto Encyclopedia of Genes and Genomes showed that lysosome and longevity regulating pathways related to digestion, detoxification, and longevity were enriched. We suggest that C. papaya is a more suitable host than S. tuberosum, and downregulated target genes may have important effects on the adaptation of P. marginatus to host transfer.

Family dinner: Transcriptional plasticity of five Noctuidae (Lepidoptera) feeding on three host plant species

Polyphagous insects often show specialization in feeding on different host plants in terms of survival and growth and, therefore, can be considered minor or major pests of particular hosts. Whether polyphagous insects employ a common transcriptional response to cope with defenses from diverse host plants is under-studied. We focused on patterns of transcriptional plasticity in polyphagous moths (Noctuidae), of which many species are notorious pests, in relation to herbivore performance on different host plants. We compared the transcriptional plasticity of five polyphagous moth species feeding and developing on three different host plant species. Using a comparative phylogenetic framework, we evaluated if successful herbivory, as measured by larval performance, is determined by a shared or lineage-specific transcriptional response. The upregulated transcriptional activity, or gene expression pattern, of larvae feeding on the different host plants and artificial control diet was highly plastic and moth species-specific. Specialization, defined as high herbivore success for specific host plants, was not generally linked to a lower number of induced genes. Moths that were more distantly related and showing high herbivore success for certain host plants showed shared expression of multiple homologous genes, indicating convergence. We further observed specific transcriptional responses within phylogenetic lineages. These expression patterns for specific host plant species are likely caused by shared evolutionary histories, for example, symplesiomorphic patterns, and could therefore not be associated with herbivore success alone. Multiple gene families, with roles in plant digestion and detoxification, were widely expressed in response to host plant feeding but again showed highly moth species-specific. Consequently, high herbivore success for specific host plants is also driven by species-specific transcriptional plasticity. Thus, potential pest moths display a complex and species-specific transcriptional plasticity.

Evolution of resource generalism via generalized stress response confers increased reproductive thermal tolerance in a pest beetle

Generalism should be favoured evolutionarily when there is no genetic constraint or loss of fitness across alternative environments. However, evolution of generalism can require substantial evolutionary change, which can confer a general stress response to other aspects of the environment. We created generalist lineages from an ancestral, resource-specialized laboratory population of seed beetles (Callosobruchus maculatus) by rearing lines over 60 generations on a mixture of both ancestral and novel host species to test for costs associated with the evolution of generalism involving evolutionary changes in gene expression and correlated phenotypic responses during a shift to generalism. Evolved lines had higher fitness on the novel resource, with no loss of fitness on the ancestral resource, indicating that they overcame initial fitness trade-offs. This involved upregulation of major stress response (heat shock protein) genes and genes coding for metabolic enzymes, suggesting an underpinning metabolic and physiological cost. Resource generalist populations also evolved greater thermal tolerance breadth, highlighting that the evolution of resource generalism might pre-adapt species to respond favourably to other environmental stressors, following selection for generalized stress response gene upregulation. The rapid gain of novel hosts during a pest invasion might also confer greater thermal resilience to ongoing climate change.

Evolved transcriptional responses and their trade-offs after long-term adaptation of Bemisia tabaci to a marginally-suitable host

Although generalist insect herbivores can migrate and rapidly adapt to a broad range of host plants, they can face significant difficulties when accidentally migrating to novel and marginally-suitable hosts. What happens, both in performance and gene expression regulation, if these marginally-suitable hosts must be used for multiple generations before migration to a suitable host can take place, largely remains unknown. In this study, we established multigenerational colonies of the whitefly Bemisia tabaci, a generalist phloem-feeding species, adapted to a marginally-suitable host (habanero pepper) or an optimal host (cotton). We used reciprocal host tests to estimate the differences in performance of the populations on both hosts under optimal (30 oC) and mild-stressful (24 oC) temperature conditions, and documented the associated transcriptomic changes. The habanero pepper-adapted population greatly improved its performance on habanero pepper but did not reach its performance level on cotton, the original host. It also showed reduced performance on cotton, relative to the non-adapted population, and an antagonistic effect of the lower-temperature stressor. The transcriptomic data revealed that most of the expression changes, associated with long-term adaptation to habanero pepper, can be categorized as "evolved" with no initial plastic response. Three molecular functions dominated: enhanced formation of cuticle structural constituents, enhanced activity of oxidation-reduction processes involved in neutralization of phytotoxins and reduced production of proteins from the cathepsin B family. Taken together, these findings indicate that generalist insects can adapt to novel host plants by modifying the expression of a relatively small set of specific molecular functions.

Transcriptional responses are oriented towards different components of the rearing environment in two Drosophila sibling species

BACKGROUND

The chance to compare patterns of differential gene expression in related ecologically distinct species can be particularly fruitful to investigate the genetics of adaptation and phenotypic plasticity. In this regard, a powerful technique such as RNA-Seq applied to ecologically amenable taxa allows to address issues that are not possible in classic model species. Here, we study gene expression profiles and larval performance of the cactophilic siblings Drosophila buzzatii and D. koepferae reared in media that approximate natural conditions and evaluate both chemical and nutritional components of the diet. These closely related species are complementary in terms of host-plant use since the primary host of one is the secondary of the other. D. koepferae is mainly a columnar cactus dweller while D. buzzatii prefers Opuntia hosts.

RESULTS

Our comparative study shows that D. buzzatii and D. koepferae have different transcriptional strategies to face the challenges posed by their natural resources. The former has greater transcriptional plasticity, and its response is mainly modulated by alkaloids of its secondary host, while the latter has a more canalized genetic response, and its transcriptional plasticity is associated with the cactus species.

CONCLUSIONS

Our study unveils a complex pleiotropic genetic landscape in both species, with functional links that relate detox responses and redox mechanisms with developmental and neurobiological processes. These results contribute to deepen our understanding of the role of host plant shifts and natural stress driving ecological specialization.

Role of Genes in Regulating Host Plants Expansion in Tephritid Fruit Flies (Diptera) and Potential for RNAi-Based Control

Host plant expansion is an important survival strategy for tephritids as they expand their range. Successful host expansion requires tephritids to adapt to the chemical and nonchemical properties of a novel host fruit, such as fruit color, phenology, and phytochemicals. These plant properties trigger a series of processes in tephritids, with each process having its own genetic basis, which means that various genes are involved in regulating host plant expansion by tephritids. This review summarizes current knowledge on the categories and roles of genes involved in host plant expansion in several important tephritid species, including genes related to chemoreception (olfactory and gustation), vision, digestion, detoxification, development, ribosomal and energy metabolism. Chemoreception- and detoxification- and digestion-related genes are stimulated by volatile chemicals and secondary chemicals of different hosts, respectively, which are involved in the regulation of nervous signal transduction that triggers behavioral, physical, and chemical responses to the novel host fruit. Vision-, nerve-, and development-related genes and metabolism-associated genes are activated in response to nonchemical stimuli from different hosts, such as color and phenology, to regulate a comprehensive adaptation of the extending host for tephritids. The chemical and nonchemical signals of hosts activate ribosomal and energy-related genes that result in the basic regulation of many processes of host expansion, including detoxification and development. These genes do not regulate novel host use individually, but multiple genes regulate multilevel adaptation to novel host fruits via multiple mechanisms. These genes may also be potential target genes for RNAi-based control of tephritid pests.

Transcriptomic modulation in response to an intoxication with deltamethrin in a population of Triatoma infestans with low resistance to pyrethroids

Background

Triatoma infestans is the main vector of Chagas disease in the Southern Cone. The resistance to pyrethroid insecticides developed by populations of this species impairs the effectiveness of vector control campaigns in wide regions of Argentina. The study of the global transcriptomic response to pyrethroid insecticides is important to deepen the knowledge about detoxification in triatomines.

Methodology and findings

We used RNA-Seq to explore the early transcriptomic response after intoxication with deltamethrin in a population of T. infestans which presents low resistance to pyrethroids. We were able to assemble a complete transcriptome of this vector and found evidence of differentially expressed genes belonging to diverse families such as chemosensory and odorant-binding proteins, ABC transporters and heat-shock proteins. Moreover, genes related to transcription and translation, energetic metabolism and cuticle rearrangements were also modulated. Finally, we characterized the repertoire of previously uncharacterized detoxification-related gene families in T. infestans and Rhodnius prolixus.

Conclusions and significance

Our work contributes to the understanding of the detoxification response in vectors of Chagas disease. Given the absence of an annotated genome from T. infestans, the analysis presented here constitutes a resource for molecular and physiological studies in this species. The results increase the knowledge on detoxification processes in vectors of Chagas disease, and provide relevant information to explore undescribed potential insecticide resistance mechanisms in populations of these insects.

Climate-driven variation in biotic interactions provides a narrow and variable window of opportunity for an insect herbivore at its ecological margin

Climate-driven geographic range shifts have been associated with transitions between dietary specialism and generalism at range margins. The mechanisms underpinning these often transient niche breadth modifications are poorly known, but utilization of novel resources likely depends on phenological synchrony between the consumer and resource. We use a climate-driven range and host shift by the butterfly Aricia agestis to test how climate-driven changes in host phenology and condition affect phenological synchrony, and consider implications for host use. Our data suggest that the perennial plant that was the primary host before range expansion is a more reliable resource than the annual Geraniaceae upon which the butterfly has become specialized in newly colonized parts of its range. In particular, climate-driven phenological variation in the novel host Geranium dissectum generates a narrow and variable 'window of opportunity' for larval productivity in summer. Therefore, although climatic change may allow species to shift hosts and colonise novel environments, specialization on phenologically limited hosts may not persist at ecological margins as climate change continues. We highlight the potential role for phenological (a)synchrony in determining lability of consumer-resource associations at range margins and the importance of considering causes of synchrony in biotic interactions when predicting range shifts. This article is part of the theme issue 'Species' ranges in the face of changing environments (Part II)'.

Unraveling the role of male reproductive tract and haemolymph in cantharidin-exuding Lydus trimaculatus and Mylabris variabilis (Coleoptera: Meloidae): a comparative transcriptomics approach

Background

Meloidae (blister beetles) are known to synthetize cantharidin (CA), a toxic and defensive terpene mainly stored in male accessory glands (MAG) and emitted outward through reflex-bleeding. Recent progresses in understanding CA biosynthesis and production organ(s) in Meloidae have been made, but the way in which self-protection is achieved from the hazardous accumulation and release of CA in blister beetles has been experimentally neglected. To provide hints on this pending question, a comparative de novo assembly transcriptomic approach was performed by targeting two tissues where CA is largely accumulated and regularly circulates in Meloidae: the male reproductive tract (MRT) and the haemolymph. Differential gene expression profiles in these tissues were examined in two blister beetle species, Lydus trimaculatus (Fabricius, 1775) (tribe Lyttini) and Mylabris variabilis (Pallas, 1781) (tribe Mylabrini). Upregulated transcripts were compared between the two species to identify conserved genes possibly involved in CA detoxification and transport.

Results

Based on our results, we hypothesize that, to avoid auto-intoxication, ABC, MFS or other solute transporters might sequester purported glycosylated CA precursors into MAG, and lipocalins could bind CA and mitigate its reactivity when released into the haemolymph during the autohaemorrhaging response. We also found an over-representation in haemolymph of protein-domains related to coagulation and integument repairing mechanisms that likely reflects the need to limit fluid loss during reflex-bleeding.

Conclusions

The de novo assembled transcriptomes of L. trimaculatus and M. variabilis here provided represent valuable genetic resources to further explore the mechanisms employed to cope with toxicity of CA in blister beetle tissues. These, if revealed, might help conceiving safe and effective drug-delivery approaches to enhance the use of CA in medicine.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08118-8.

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.

Behavioral and Physiological Plasticity Provides Insights into Molecular Based Adaptation Mechanism to Strain Shift in Spodoptera frugiperda

How herbivorous insects adapt to host plants is a key question in ecological and evolutionary biology. The fall armyworm, (FAW) Spodoptera frugiperda (J.E. Smith), although polyphagous and a major pest on various crops, has been reported to have a rice and corn (maize) feeding strain in its native range in the Americas. The species is highly invasive and has recently established in China. We compared behavioral changes in larvae and adults of a corn population (Corn) when selected on rice (Rice) and the molecular basis of these adaptational changes in midgut and antennae based on a comparative transcriptome analysis. Larvae of S. frugiperda reared on rice plants continuously for 20 generations exhibited strong feeding preference for with higher larval performance and pupal weight on rice than on maize plants. Similarly, females from the rice selected population laid significantly more eggs on rice as compared to females from maize population. The most highly expressed DEGs were shown in the midgut of Rice vs. Corn. A total of 6430 DEGs were identified between the populations mostly in genes related to digestion and detoxification. These results suggest that potential adaptations for feeding on rice crops, may contribute to the current rapid spread of fall armyworm on rice crops in China and potentially elsewhere. Consistently, highly expressed DEGs were also shown in antennae; a total of 5125 differentially expressed genes (DEGs) s were identified related to the expansions of major chemosensory genes family in Rice compared to the Corn feeding population. These results not only provide valuable insight into the molecular mechanisms in host plants adaptation of S. frugiperda but may provide new gene targets for the management of this pest.

Host-specific gene expression as a tool for introduction success in Naupactus parthenogenetic weevils

Food resource access can mediate establishment success in invasive species, and generalist herbivorous insects are thought to rely on mechanisms of transcriptional plasticity to respond to dietary variation. While asexually reproducing invasives typically have low genetic variation, the twofold reproductive capacity of asexual organisms is a marked advantage for colonization. We studied host-related transcriptional acclimation in parthenogenetic, invasive, and polyphagous weevils: Naupactus cervinus and N. leucoloma. We analyzed patterns of gene expression in three gene categories that can mediate weevil-host plant interactions through identification of suitable host plants, short-term acclimation to host plant defenses, and long-term adaptation to host plant defenses and their pathogens. This approach employed comparative transcriptomic methods to investigate differentially expressed host detection, detoxification, immune defense genes, and pathway-level gene set enrichment. Our results show that weevil gene expression responses can be host plant-specific, and that elements of that response can be maintained in the offspring. Some host plant groups, such as legumes, appear to be more taxing as they elicit a complex gene expression response which is both strong in intensity and specific in identity. However, the weevil response to taxing host plants shares many differentially expressed genes with other stressful situations, such as host plant cultivation conditions and transition to novel host, suggesting that there is an evolutionarily favorable shared gene expression regime for responding to different types of stressful situations. Modulating gene expression in the absence of other avenues for phenotypic adaptation may be an important mechanism of successful colonization for these introduced insects.

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Transcriptomic analysis and molecular docking reveal genes involved in the response of Aedes aegypti larvae to an essential oil extracted from Eucalyptus

BACKGROUND

Aedes aegypti (L.) is an urban mosquito, vector of several arboviruses that cause severe diseases in hundreds of million people each year. The resistance to synthetic insecticides developed by Ae. aegypti populations worldwide has contributed to failures in vector control campaigns, increasing the impact of arbovirus diseases. In this context, plant-derived essential oils with larvicidal activity could be an attractive alternative for vector control. However, the mode of action and the detoxificant response of mosquitoes to plant derived compounds have not been established, impairing the optimization of their use.

METHODS AND FINDINGS

Here we compare gene expression in Ae. aegypti larvae after 14 hrs of exposure to Eucalyptus camaldulensis essential oil with a control group exposed to vehicle (acetone) for the same lapse, by using RNA-Seq. We found differentially expressed genes encoding for cuticle proteins, fatty-acid synthesis, membrane transporters and detoxificant related gene families (i.e. heat shock proteins, cytochromes P450, glutathione transferases, UDP-glycosyltransferases and ABC transporters). Finally, our RNA-Seq and molecular docking results provide evidence pointing to a central involvement of chemosensory proteins in the detoxificant response in mosquitoes.

CONCLUSIONS AND SIGNIFICANCE

Our work contributes to the understanding of the physiological response of Ae. aegypti larvae to an intoxication with a natural toxic distilled from Eucalyptus leafs. The results suggest an involvement of most of the gene families associated to detoxification of xenobiotics in insects. Noteworthy, this work provides important information regarding the implication of chemosensory proteins in the detoxification of a natural larvicide. Understanding the mode of detoxification of Eucalyptus distilled compounds could contribute to their implementation as a tool in mosquito control.

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

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

Chromosome Level Assembly of the Comma Butterfly (Polygonia c-album)

The comma butterfly (Polygonia c-album, Nymphalidae, Lepidoptera) is a model insect species, most notably in the study of phenotypic plasticity and plant-insect coevolutionary interactions. In order to facilitate the integration of genomic tools with a diverse body of ecological and evolutionary research, we assembled the genome of a Swedish comma using 10X sequencing, scaffolding with matepair data, genome polishing, and assignment to linkage groups using a high-density linkage map. The resulting genome is 373 Mb in size, with a scaffold N50 of 11.7 Mb and contig N50 of 11,2Mb. The genome contained 90.1% of single-copy Lepidopteran orthologs in a BUSCO analysis of 5,286 genes. A total of 21,004 gene-models were annotated on the genome using RNA-Seq data from larval and adult tissue in combination with proteins from the Arthropoda database, resulting in a high-quality annotation for which functional annotations were generated. We further documented the quality of the chromosomal assembly via synteny assessment with Melitaea cinxia. The resulting annotated, chromosome-level genome will provide an important resource for investigating coevolutionary dynamics and comparative analyses in Lepidoptera.

Analyses of microRNAs and transcriptomes in the midgut of Spodoptera litura feeding on Brassica juncea

Spodoptera litura is a destructive agricultural pest in tropical and subtropical areas. Understanding the molecular mechanisms of S. litura adaptation to its preferred host plants may help identify target genes useful for pest control. We used high-throughput sequencing to characterize the expression patterns of messenger RNAs (mRNAs) and microRNAs (miRNAs) in the midgut of S. litura fed on Brassica juncea for 6 h and 48 h. A total of 108 known and 134 novel miRNAs were identified, 29 miRNAs and 237 mRNAs were differentially expressed at 6 h of B. juncea feeding, 26 miRNAs and 433 mRNAs were differentially expressed at 48 h. For the mRNAs, the up-regulated genes were mostly enriched in detoxification enzymes (cytochrome P450, esterase, glutathione S-transferase, uridine diphosphate-glucuronosyl transferase), while the down-regulated genes were mostly enriched in proteinases and immune-related genes. Furthermore, most detoxification enzymes begin to up-regulate at 6 h, while most digestion and immune-related genes begin to up- or down-regulate at 48 h. Eighteen and 37 differently expressed transcription factors were identified at 6 h and 48 h, which may regulate the functional genes. We acquired 136 and 41 miRNA versus mRNA pairs at 6 h and 48 h, respectively. Some down-regulated and up-regulated miRNAs were predicted to target detoxification enzymes and proteinases, respectively. Real-time quantitative polymerase chain reaction of nine randomly selected miRNAs and 28 genes confirmed the results of RNA-seq. This analyses of miRNA and mRNA transcriptomes provides useful information about the molecular mechanisms of S. litura response to B. juncea.

Comparative transcriptome analysis of the newly discovered insect vector of the pine wood nematode in China, revealing putative genes related to host plant adaptation

BACKGROUND

In many insect species, the larvae/nymphs are unable to disperse far from the oviposition site selected by adults. The Sakhalin pine sawyer Monochamus saltuarius (Gebler) is the newly discovered insect vector of the pine wood nematode (Bursaphelenchus xylophilus) in China. Adult M. saltuarius prefers to oviposit on the host plant Pinus koraiensis, rather than P. tabuliformis. However, the genetic basis of adaptation of the larvae of M. saltuarius with weaken dispersal ability to host environments selected by the adult is not well understood.

RESULTS

In this study, the free amino and fatty acid composition and content of the host plants of M. saltuarius larvae, i.e., P. koraiensis and P. tabuliformis were investigated. Compared with P. koraiensis, P. tabuliformis had a substantially higher content of various free amino acids, while the opposite trend was detected for fatty acid content. The transcriptional profiles of larval populations feeding on P. koraiensis and P. tabuliformis were compared using PacBio Sequel II sequencing combined with Illumina sequencing. The results showed that genes relating to digestion, fatty acid synthesis, detoxification, oxidation-reduction, and stress response, as well as nutrients and energy sensing ability, were differentially expressed, possibly reflecting adaptive changes of M. saltuarius in response to different host diets. Additionally, genes coding for cuticle structure were differentially expressed, indicating that cuticle may be a potential target for plant defense. Differential regulation of genes related to the antibacterial and immune response were also observed, suggesting that larvae of M. saltuarius may have evolved adaptations to cope with bacterial challenges in their host environments.

CONCLUSIONS

The present study provides comprehensive transcriptome resource of M. saltuarius relating to host plant adaptation. Results from this study help to illustrate the fundamental relationship between transcriptional plasticity and adaptation mechanisms of insect herbivores to host plants.

The molecular mechanisms that determine different degrees of polyphagy in the Bemisia tabaci species complex

The whitefly Bemisia tabaci is a closely related group of >35 cryptic species that feed on the phloem sap of a broad range of host plants. Species in the complex differ in their host-range breadth, but the mechanisms involved remain poorly understood. We investigated, therefore, how six different B. tabaci species cope with the environmental unpredictability presented by a set of four common and novel host plants. Behavioral studies indicated large differences in performances on the four hosts and putative specialization of one of the species to cassava plants. Transcriptomic analyses revealed two main insights. First, a large set of genes involved in metabolism (>85%) showed differences in expression between the six species, and each species could be characterized by its own unique expression pattern of metabolic genes. However, within species, these genes were constitutively expressed, with a low level of environmental responsiveness (i.e., to host change). Second, within each species, sets of genes mainly associated with the super-pathways "environmental information processing" and "organismal systems" responded to the host switching events. These included genes encoding for proteins involved in sugar homeostasis, signal transduction, membrane transport, and immune, endocrine, sensory and digestive responses. Our findings suggested that the six B. tabaci species can be divided into four performance/transcriptomic "Types" and that polyphagy can be achieved in multiple ways. However, polyphagy level is determined by the specific identity of the metabolic genes/pathways that are enriched and overexpressed in each species (the species' individual metabolic "tool kit").

Consequences of 'no-choice, fixed time' reciprocal host plant switches on nutrition and gut serine protease gene expression in Pieris brassicae L. (Lepidoptera: Pieridae)

Rapid adaptive responses were evident from reciprocal host-plant switches on performance, digestive physiology and relative gene expression of gut serine proteases in larvae of crucifer pest P. brassicae transferred from cauliflower (CF, Brassica oleracea var. botrytis, family Brassicaceae) to an alternate host, garden nasturtium, (GN, Tropaeolum majus L., family Tropaeolaceae) and vice-versa under laboratory conditions. Estimation of nutritional indices indicated that larvae of all instars tested consumed the least food and gained less weight on CF-GN diet (significant at p≤0.05) as compared to larvae feeding on CF-CF, GN-GN and GN-CF diets suggesting that the switch to GN was nutritionally less favorable for larval growth. Nevertheless, these larvae, especially fourth instars, were adroit in utilizing and digesting GN as a new host plant type. In vitro protease assays conducted to understand associated physiological responses within twelve hours indicated that levels and properties of gut proteases were significantly influenced by type of natal host-plant consumed, change in diet as well as larval age. Activities of gut trypsins and chymotrypsins in larvae feeding on CF-GN and GN-CF diets were distinct, and represented shifts toward profiles observed in larvae feeding continuously on GN-GN and CF-CF diets respectively. Results with diagnostic protease inhibitors like TLCK, STI and SBBI in these assays and gelatinolytic zymograms indicated complex and contrasting trends in gut serine protease activities in different instars from CF-GN diet versus GN-CF diet, likely due to ingestion of plant protease inhibitors present in the new diet. Cloning and sequencing of serine protease gene fragments expressed in gut tissues of fourth instar P. brassicae revealed diverse transcripts encoding putative trypsins and chymotrypsins belonging to at least ten lineages. Sequences of members of each lineage closely resembled lepidopteran serine protease orthologs including uncharacterized transcripts from Pieris rapae. Differential regulation of serine protease genes (Pbr1-Pbr5) was observed in larval guts of P. brassicae from CF-CF and GN-GN diets while expression of transcripts encoding two putative trypsins (Pbr3 and Pbr5) were significantly different in larvae from CF-GN and GN-CF diets. These results suggested that some gut serine proteases that were differentially expressed in larvae feeding on different species of host plants were also involved in rapid adaptations to dietary switches. A gene encoding nitrile-specifier protein (nsp) likely involved in detoxification of toxic products from interactions of ingested host plant glucosinolates with myrosinases was expressed to similar levels in these larvae. Taken together, these snapshots reflected contrasts in physiological and developmental plasticity of P. brassicae larvae to nutritional challenges from wide dietary switches in the short term and the prominent role of gut serine proteases in rapid dietary adaptations. This study may be useful in designing novel management strategies targeting candidate gut serine proteases of P. brassicae using RNA interference, gene editing or crops with transgenes encoding protease inhibitors from taxonomically-distant host plants.

Host plant diet affects growth and induces altered gene expression and microbiome composition in the wood white (Leptidea sinapis) butterfly

In a time with decreasing biodiversity, especially among insects, a detailed understanding about specific resource utilization strategies is crucial. The physiological and behavioural responses to host switches in phytophagous insects are poorly understood. Earlier studies indicate that a host plant switch might be associated with distinctive molecular and physiological responses in different lineages. Expanding the assessment of such associations across Lepidoptera will reveal if there are general patterns in adaptive responses, or if each switch event is more of a unique character. We investigated host plant preference, fitness consequences, effects on expression profiles and gut microbiome composition in two common wood white (Leptidea sinapis) populations with different host plant preferences from the extremes of the species distribution area (Sweden and Catalonia). Our results show that female Catalonian wood whites lack preference for either host plant (Lotus corniculatus or L. dorycnium), while Swedish females laid significantly more eggs on L. corniculatus. Individuals from both populations reared on L. dorycnium had longer developmental times and smaller body size as adults. This indicates that both environmental and genetic factors determine the choice to use a specific host plant. Gene expression analysis revealed a more pronounced response to host plant in the Catalonian compared to the Swedish population. In addition, host plant treatment resulted in a significant shift in microbiome community structure in the Catalonian population. Together, this suggests that population specific plasticity associated with local conditions underlies host plant utilisation in wood whites.

Genome-Wide Identification and Expression Analysis of Udp-Glucuronosyltransferases in the Whitefly Bemisia Tabaci (Gennadius) (HemipterA: Aleyrodidae)

Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is an important agricultural pest worldwide. Uridine diphosphate (UDP)-glucuronosyltransferases (UGTs) are one of the largest and most ubiquitous groups of proteins. Because of their role in detoxification, insect UGTs are attracting increasing attention. In this study, we identified and analyzed UGT genes in B. tabaci MEAM1 to investigate their potential roles in host adaptation and reproductive capacity. Based on phylogenetic and structural analyses, we identified 76 UGT genes in the B. tabaci MEAM1 genome. RNA-seq and real-time quantitative PCR (RT-qPCR) revealed differential expression patterns of these genes at different developmental stages and in association with four host plants (cabbage, cucumber, cotton and tomato). RNA interference results of selected UGTs showed that, when UGT352A1, UGT352B1, and UGT354A1 were respectively silenced by feeding on dsRNA, the fecundity of B. tabaci MEAM1 was reduced, suggesting that the expressions of these three UGT genes in this species may be associated with host-related fecundity. Together, our results provide detailed UGTs data in B.tabaci and help guide future studies on the mechanisms of host adaptation by B.tabaci.

Genome streamlining in a minute herbivore that manipulates its host plant

The tomato russet mite, Aculops lycopersici, is among the smallest animals on earth. It is a worldwide pest on tomato and can potently suppress the host's natural resistance. We sequenced its genome, the first of an eriophyoid, and explored whether there are genomic features associated with the mite's minute size and lifestyle. At only 32.5 Mb, the genome is the smallest yet reported for any arthropod and, reminiscent of microbial eukaryotes, exceptionally streamlined. It has few transposable elements, tiny intergenic regions, and is remarkably intron-poor, as more than 80% of coding genes are intronless. Furthermore, in accordance with ecological specialization theory, this defense-suppressing herbivore has extremely reduced environmental response gene families such as those involved in chemoreception and detoxification. Other losses associate with this species' highly derived body plan. Our findings accelerate the understanding of evolutionary forces underpinning metazoan life at the limits of small physical and genome size.

Physiological adaptations to sugar-mimic alkaloids: Insights from Bombyx mori for long-term adaption and short-term response

Insects evolved adaptive plasticity to minimize the effects of the chemical defenses of their host plants. Nevertheless, the expressional response and adaptation of phytophagous specialists for long-term adaption and short-term response to host phytochemicals remains largely unexplored. The mulberry (Morus alba)-silkworm (Bombyx mori) interaction is an old and well-known model of plant-insect interaction. In this study, we examined the long-term adaption and short-term response of the mulberry-specialist silkworm to two sugar-mimic alkaloids in mulberry: the commonly encountered 1-deoxynojirimycin (1-DNJ) and occasionally encountered 1,4-dideoxy-1,4-imino-D-arabinitol (D-AB1), respectively. Global transcriptional patterns revealed that the physiological responses induced by the selective expression of genes involved in manifold cellular processes, including detoxification networks, canonical digestion processes, target enzymes, and other fundamental physiological processes, were crucial for regulating metabolic homeostasis. Comparative network analysis of the effects of exposure to D-AB1 and 1-DNJ supported the contention that B. mori produced similar and specific trajectories of changed gene expression in response to different sugar-mimic alkaloids. D-AB1 elicited a substantial proportion of downregulated genes relating to carbohydrate metabolism, catabolic process, lipid metabolism, and glycan biosynthesis and metabolism. This study dramatically expands our knowledge of the physiological adaptations to dietary sugar-mimic alkaloid intake and uncovered both metabolic evolutionarily responses and unique adaptive mechanisms previously unknown in insects.

Mechanisms of Novel Host Use by Bactrocera tau (Tephritid: Diptera) Revealed by RNA Transcriptomes

Use of novel plant hosts can facilitate the establishment and range expansion of herbivorous invasive species. However, the inherent mechanisms of novel host use are still unclear in many herbivorous species. Here, we examine mechanisms of novel host use in the invasive tephritid fruit fly Bactrocera tau (Walker)(Diptera: Tephritidae) by documenting changes in the RNA transcriptomes associated with a novel host. RNA transcripts of B. tau were obtained with high-throughput sequencing from samples continuously reared on two traditional Cucurbitaceae hosts and a novel host (banana). We found transcriptome variation was strongly associated with feeding on banana. Moreover, B. tau feeding on banana contained more differentially expressed genes (DEGs) and more annotated categories of DEGs in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database with 1,595 DEGs and 21 major annotated pathways. The annotated categories of DEGs in individuals reared on banana differed with from those individuals feeding on other hosts and were enriched in oxidative phosphorylation, citrate cycle pathway, and four other carbohydrate pathways. For B. tau feeding on banana, the predominant numbers of upregulated genes in the mitochondrial NADH (56 on average) and a relatively higher numbers of upregulated genes (13 on average) were found in oxidative phosphorylation and the TCA pathway, respectively. Changes in RNA transcriptomes associated with novel host use, especially for genes related to energy and carbohydrate metabolism, help to explain how B. tau can be successful in use of novel hosts and may be useful in developing novel strategies for control of tephritid flies.

Comparative analysis of diet-associated responses in two rice planthopper species

Background

Host adaptation is the primary determinant of insect diversification. However, knowledge of different host ranges in closely related species remains scarce. The brown planthopper (Nilaparvata lugens, BPH) and the small brown planthopper (Laodelphax striatellus, SBPH) are the most destructive insect pests within the family Delphacidae. These two species differ in their host range (SBPH can well colonize rice and wheat plants, whereas BPH survives on only rice plants), but the underlying mechanism of this difference remains unknown. High-throughput sequencing provides a powerful approach for analyzing the association between changes in gene expression and the physiological responses of insects. Therefore, gut transcriptomes were performed to elucidate the genes associated with host adaptation in planthoppers. The comparative analysis of planthopper responses to different diets will improve our knowledge of host adaptation regarding herbivorous insects.

Results

In the present study, we analyzed the change in gene expression of SBPHs that were transferred from rice plants to wheat plants over the short term (rSBPH vs tSBPH) or were colonized on wheat plants over the long term (rSBPH vs wSBPH). The results showed that the majority of differentially expressed genes in SBPH showed similar changes in expression for short-term transfer and long-term colonization. Based on a comparative analysis of BPH and SBPH after transfer, the genes associated with sugar transporters and heat-shock proteins showed similar variation. However, most of the genes were differentially regulated between the two species. The detoxification-related genes were upregulated in SBPH after transfer from the rice plants to the wheat plants, but these genes were downregulated in BPH under the same conditions. In contrast, ribosomal-related genes were downregulated in SBPH after transfer, but these genes were upregulated in BPH under the same conditions.

Conclusions

The results of this study provide evidence that host plants played a dominant role in shaping gene expression and that the low fitness of BPH on wheat plants might be determined within 24 h after transfer. This study deepens our understanding of different host ranges for the two planthopper species, which may provide a potential strategy for pest management.

Host-specialized transcriptome of plant-associated organisms

Living organisms respond to their immediate environment by modulating their genetic programme to perform adapted functions. Eukaryotic organisms that associate with plants (fungi, oomycetes, insects, …) alter their transcriptome in a host-specific manner. Recent comparative transcriptomic studies revealed that host-specialized transcriptomes consist of a limited set of genes. Such a set typically encodes proteins that modulate host structures and functions (predicted effectors and other secreted proteins), control nutrient assimilation (proteases, transporters), and maintain cellular homeostasis (oxidoreductases, detoxification enzymes). We conclude by discussing open mechanistic and evolutionary questions and integrated approaches to move beyond descriptive studies.

Legacy of a Butterfly's Parental Microbiome in Offspring Performance

An insect's phenotype can be influenced by the experiences of the parental generation. However, the effects of the parental symbiotic microbiome and host plant use on the offspring are unclear. We addressed this gap of knowledge by studying Pieris brassicae, a multivoltine butterfly species feeding on different brassicaceous plants across generations. We investigated how disturbance of the parental bacterial community by antibiotic treatment affects F1 larval traits. We tested the effects depending on whether F1 larvae are feeding on the same plant species as their parents or on a different one. The parental treatment alone had no impact on the biomass of F1 larvae feeding on the parental plant species. However, the parental treatment had a detrimental effect on F1 larval biomass when F1 larvae had a different host plant than their parents. This effect was linked to higher larval prophenoloxidase activity and greater downregulation of the major allergen gene (MA), a glucosinolate detoxification gene of P. brassicae Bacterial abundance in untreated adult parents was high, while it was very low in F1 larvae from either parental type, and thus unlikely to directly influence larval traits. Our results suggest that transgenerational effects of the parental microbiome on the offspring's phenotype become evident when the offspring is exposed to a transgenerational host plant shift.IMPORTANCE Resident bacterial communities are almost absent in larvae of butterflies and thus are unlikely to affect their host. In contrast, adult butterflies contain conspicuous amounts of bacteria. While the host plant and immune state of adult parental butterflies are known to affect offspring traits, it has been unclear whether also the parental microbiome imposes direct effects on the offspring. Here, we show that disturbance of the bacterial community in parental butterflies by an antibiotic treatment has a detrimental effect on those offspring larvae feeding on a different host plant than their parents. Hence, the study indicates that disturbance of an insect's parental microbiome by an antibiotic treatment shapes how the offspring individuals can adjust themselves to a novel host plant.

Host use diversification during range shifts shapes global variation in Lepidopteran dietary breadth

Niche breadths tend to be greater at higher latitudes. This pattern is frequently assumed to emerge from the cumulative effects of multiple, independent local adaptation events along latitudinal environmental gradients, although evidence that generalization is more beneficial at higher-latitude locations remains equivocal. Here I propose an alternative hypothesis: that latitudinal variation in niche breadths emerges as a non-adaptive consequence of range shift dynamics. Based on analysis of a global dataset comprising more than 6,934 globally distributed dietary records from 4,410 Lepidopteran species, this hypothesis receives robust support. Population-level dietary niche breadths are better explained by the relative position of the population within its geographic range and the species' poleward range extent than by the latitude of diet observation. Broader diets are observed closer to poleward range limits and in species that have attained higher latitudes. Moreover, latitudinal variation in diet breadth is more prominent within and among species undergoing rapid, contemporary range shifts than for species with more stable ranges. Together these results suggest that latitudinal patterns in niche breadth represent a transient and emergent property of recent geographic range dynamics and need not require underlying gradients in selective agents or fitness trade-offs. The results have wide-ranging implications for global ecology and for anticipating changes in host use during ongoing distributional shifts of pests and disease vectors.

Host plant adaptation in the polyphagous whitefly, Trialeurodes vaporariorum, is associated with transcriptional plasticity and altered sensitivity to insecticides

BACKGROUND

The glasshouse whitefly, Trialeurodes vaporariorum, is a damaging crop pest and an invasive generalist capable of feeding on a broad range of host plants. As such this species has evolved mechanisms to circumvent the wide spectrum of anti-herbivore allelochemicals produced by its host range. T. vaporariorum has also demonstrated a remarkable ability to evolve resistance to many of the synthetic insecticides used for control.

RESULTS

To gain insight into the molecular mechanisms that underpin the polyphagy of T. vaporariorum and its resistance to natural and synthetic xenobiotics, we sequenced and assembled a reference genome for this species. Curation of genes putatively involved in the detoxification of natural and synthetic xenobiotics revealed a marked reduction in specific gene families between this species and another generalist whitefly, Bemisia tabaci. Transcriptome profiling of T. vaporariorum upon transfer to a range of different host plants revealed profound differences in the transcriptional response to more or less challenging hosts. Large scale changes in gene expression (> 20% of genes) were observed during adaptation to challenging hosts with a range of genes involved in gene regulation, signalling, and detoxification differentially expressed. Remarkably, these changes in gene expression were associated with significant shifts in the tolerance of host-adapted T. vaporariorum lines to natural and synthetic insecticides.

CONCLUSIONS

Our findings provide further insights into the ability of polyphagous insects to extensively reprogram gene expression during host adaptation and illustrate the potential implications of this on their sensitivity to synthetic insecticides.

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

Characterizing potential repelling volatiles for "push-pull" strategy against stem borer: a case study in Chilo auricilius

Background

Massive techniques have been evaluated for developing different pest control methods to minimize fertilizer and pesticide inputs. As “push-pull” strategy utilizes generally non-toxic chemicals to manipulate behaviors of insects, such strategy is considered to be environmentally friendly. “Push-pull” strategy has been extraordinarily effective in controlling stem borers, and the identification of new “pushing” or “pull” components against stem borers could be significantly helpful.

Results

In this study, the results of field trapping assay and behavioral assay showed the larvae of C.auricilius, one kind of stem borers, could be deterred by rice plant under tilling stage, its main host crop. The profiles of volatiles were compared between rice plants under two different developmental stages, and α-pinene was identified as a key differential component. The repelling activity of α-pinene against C.auricilius was confirmed by Y-tube olfactometer. For illuminating the olfactory recognition mechanism, transcriptome analysis was carried out, and 13 chemosensory proteins (CSPs) were identified in larvae and 19 CSPs were identified in adult of C.auriciliu, which was reported for the first time in this insect. Among these identified CSPs, 4 CSPs were significantly regulated by α-pinene treatment, and CSP8 showed good binding affinity with α-pinene in vitro.

Conclusions

Overall, C.auricilius could be repelled by rice plant at tilling stage, and our results highlighted α-pinene as a key component in inducing repelling activity at this specific stage and confirmed the roles of some candidate chemosensory elements in this chemo-sensing process. The results in this study could provide valuable information for chemosensory mechanism of C.auricilius and for identification of “push” agent against rice stem borers.

Divergent Switchgrass Cultivars Modify Cereal Aphid Transcriptomes

Schizaphis graminum Rondani (Hemiptera: Aphididae) and Sipha flava Forbes (Hemiptera: Aphididae) are two common pests of bioenergy grasses. Despite the fact that they are both considered generalists, they differ in their ability to colonize Panicum virgatum cultivars. For example, S. flava colonizes both P. virgatum cv. Summer and P. virgatum cv. Kanlow whereas S. graminum can only colonize Summer. To study the molecular responses of these aphids to these two switchgrass cultivars, we generated de novo transcriptome assemblies and compared the expression profiles of aphids feeding on both cultivars to profiles associated with feeding on a highly susceptible sorghum host and a starvation treatment. Transcriptome assemblies yielded 8,428 and 8,866 high-quality unigenes for S. graminum and S. flava, respectively. Overall, S. graminum responded strongly to all three treatments after 12 h with an upregulation of unigenes coding for detoxification enzymes while major transcriptional changes were not observed in S. flava until 24 h. Additionally, while the two aphids responded to the switchgrass feeding treatment by downregulating unigenes linked to growth and development, their responses to Summer and Kanlow diverged significantly. Schizaphis graminum upregulated more unigenes coding for stress-responsive enzymes in the Summer treatment compared to S. flava; however, many of these unigenes were actually downregulated in the Kanlow treatment. In contrast, S. flava appeared capable of overcoming host defenses by upregulating a larger number of unigenes coding for detoxification enzymes in the Kanlow treatment. Overall, these findings are consistent with previous studies on the interactions of these two cereal aphids to divergent switchgrass hosts.

Differential regulation of host plant adaptive genes in Pieris butterflies exposed to a range of glucosinolate profiles in their host plants

Specialist herbivores have often evolved highly sophisticated mechanisms to counteract defenses mediated by major plant secondary-metabolites. Plant species of the herbivore host range often display high chemical diversity and it is not well understood how specialist herbivores respond to this chemical diversity. Pieris larvae overcome toxic products from glucosinolate hydrolysis, the major chemical defense of their Brassicaceae hosts, by expressing nitrile-specifier proteins (NSP) in their gut. Furthermore, Pieris butterflies possess so-called major allergen (MA) proteins, which are multi-domain variants of a single domain major allergen (SDMA) protein expressed in the guts of Lepidopteran larvae. Here we show that Pieris larvae fine-tune NSP and MA gene expression depending on the glucosinolate profiles of their Brassicaceae hosts. Although the role of MA is not yet fully understood, the expression levels of NSP and MA in larvae that fed on plants whose glucosinolate composition varied was dramatically changed, whereas levels of SDMA expression remained unchanged. In addition, we found a similar regulation pattern among these genes in larvae feeding on Arabidopsis mutants with different glucosinolate profiles. Our results demonstrate that Pieris larvae appear to use different host plant adaptive genes to overcome a wide range of glucosinolate profiles in their host plants.

Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution

Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies—swallowtails of the family Papilionidae—that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6–142.7 Ma) and Papilionidae around 71.4 Ma (49.8–103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.

De novo comparative transcriptome analysis of a rare cicada, with identification of candidate genes related to adaptation to a novel host plant and drier habitats

Background

Although the importance of host plant chemistry in plant–insect interactions is widely recognized, our understanding about the genetic basis underlying the relationship between changes in midgut proteins and adaptation of plant-feeding insects to novel host plants and habitats is very limited. To address this knowledge gap, the transcriptional profiles of midguts among three populations of the cicada Subpsaltria yangi Chen were compared. Among which, the Hancheng (HC) and Fengxiang (FX) populations occurring in the Loess Plateau feed on Ziziphus jujuba Mill. var. spinosa (Bunge) Hu ex H. F. Chow, while the population occurring in a much drier habitat in the Helan (HL) Mountains is locally specialized on a chemically divergent plant, Ephedra lepidosperma C. Y. Cheng.

Results

Based on comparative analysis, 1826 (HL vs HC) differentially expressed genes (DEGs) and 723 DEGs (HL vs FX) were identified between the populations utilizing different host plants, including 20, 36, 2, 5 and 2 genes related to digestion, detoxification, oxidation-reduction, stress response and water-deprivation response, respectively, and 35 genes presumably associated with osmoregulation. However, only 183 DEGs were identified between the HC and FX populations, including two genes related to detoxification, two genes related to stress response, and one gene presumably associated with osmoregulation. These results suggest that the weakest expression differences were between the populations utilizing the same host plant and occurring in the closest habitats, which may help explain the metabolic mechanism of adaptation in S. yangi populations to novel host plants and new niches.

Conclusions

The observed differences in gene expression among S. yangi populations are consistent with the hypothesis that the host plant shift and habitat adaptation in the HL population was facilitated by differential regulation of genes related to digestion, detoxification, oxidation-reduction, stress response, water-deprivation response and osmoregulation. The results may inform future studies on the molecular mechanisms underlying the relationship between changes in midgut proteins and adaptation of herbivorous insects to novel host plants and new niches.

Electronic supplementary material

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

Transcriptomic Plasticity in the Arthropod Generalist Tetranychus urticae Upon Long-Term Acclimation to Different Host Plants

The two-spotted spider mite Tetranychus urticae is an important pest with an exceptionally broad host plant range. This generalist rapidly acclimatizes and adapts to a new host, hereby overcoming nutritional challenges and a novel pallet of constitutive and induced plant defenses. Although recent studies reveal that a broad transcriptomic response upon host plant transfer is associated with a generalist life style in arthropod herbivores, it remains uncertain to what extent these transcriptional changes are general stress responses or host-specific. In the present study, we analyzed and compared the transcriptomic changes that occur in a single T. urticae population upon long-term transfer from Phaseolus vulgaris to a similar, but chemically defended, host (cyanogenic Phaseolus lunatus) and to multiple economically important crops (Glycine max, Gossypium hirsutum, Solanum lycopersicum and Zea mays). These long-term host plant transfers were associated with distinct transcriptomic responses with only a limited overlap in both specificity and directionality, suggestive of a fine-tuned transcriptional plasticity. Nonetheless, analysis at the gene family level uncovered overlapping functional processes, recruiting genes from both well-known and newly discovered detoxification families. Of note, our analyses highlighted a possible detoxification role for Tetranychus-specific short-chain dehydrogenases and single PLAT domain proteins, and manual genome annotation showed that both families are expanded in T. urticae Our results shed new light on the molecular mechanisms underlying the remarkable adaptive potential for host plant use of generalist arthropods and set the stage for functional validation of important players in T. urticae detoxification of plant secondary metabolites.

Characterization and expression profiling of microRNAs in response to plant feeding in two host-plant strains of the lepidopteran pest Spodoptera frugiperda

Background

A change in the environment may impair development or survival of living organisms leading them to adapt to the change. The resulting adaptation trait may reverse, or become fixed in the population leading to evolution of species. Deciphering the molecular basis of adaptive traits can thus give evolutionary clues. In phytophagous insects, a change in host-plant range can lead to emergence of new species. Among them, Spodoptera frugiperda is a major agricultural lepidopteran pest consisting of two host-plant strains having diverged 3 MA, based on mitochondrial markers. In this paper, we address the role of microRNAs, important gene expression regulators, in response to host-plant change and in adaptive evolution.

Results

Using small RNA sequencing, we characterized miRNA repertoires of the corn (C) and rice (R) strains of S. frugiperda, expressed during larval development on two different host-plants, corn and rice, in the frame of reciprocal transplant experiments. We provide evidence for 76 and 68 known miRNAs in C and R strains and 139 and 171 novel miRNAs. Based on read counts analysis, 34 of the microRNAs were differentially expressed in the C strain larvae fed on rice as compared to the C strain larvae fed on corn. Twenty one were differentially expressed on rice compared to corn in R strain. Nine were differentially expressed in the R strain compared to C strain when reared on corn. A similar ratio of microRNAs was differentially expressed between strains on rice. We could validate experimentally by QPCR, variation in expression of the most differentially expressed candidates. We used bioinformatics methods to determine the target mRNAs of known microRNAs. Comparison with the mRNA expression profile during similar reciprocal transplant experiment revealed potential mRNA targets of these host-plant regulated miRNAs.

Conclusions

In the current study, we performed the first systematic analysis of miRNAs in Lepidopteran pests feeding on host-plants. We identified a set of the differentially expressed miRNAs that respond to the plant diet, or differ constitutively between the two host plant strains. Among the latter, the ones that are also deregulated in response to host-plant are molecular candidates underlying a complex adaptive trait.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5119-6) contains supplementary material, which is available to authorized users.