Skip to main content
Download PDF

Transcript expression plasticity as a response to alternative larval host plants in the speciation process of corn and rice strains of Spodoptera frugiperda

BMC Genomics volume 18, Article number: 792 (2017) Cite this article

Abstract

Background

Our main purpose was to evaluate the expression of plastic and evolved genes involved in ecological speciation in the noctuid moth Spodoptera frugiperda, the fall armyworm (FAW); and to demonstrate how host plants might influence lineage differentiation in this polyphagous insect. FAW is an important pest of several crops worldwide, and it is differentiated into host plant-related strains, corn (CS) and rice strains (RS). RNA-Seq and transcriptome characterization were applied to evaluate unbiased genetic expression differences in larvae from the two strains, fed on primary (corn) and alternative (rice) host plants. We consider that genes that are differently regulated by the same FAW strain, as a response to different hosts, are “plastic”. Otherwise, differences in gene expression between the two strains fed on the same host are considered constitutive differences.

Results

Individual performance parameters (larval and pupal weight) varied among conditions (strains vs. hosts). A total of 3657 contigs was related to plastic response, and 2395 contigs were differentially regulated in the two strains feeding on preferential and alternative hosts (constitutive contigs). Three molecular functions were present in all comparisons, both down- and up-regulated: oxidoreductase activity, metal-ion binding, and hydrolase activity.

Conclusions

Metabolization of foreign chemicals is among the key functions involved in the phenotypic variation of FAW strains. From an agricultural perspective, high plasticity in families of detoxifying genes indicates the capacity for a rapid response to control compounds such as insecticides.

Background

The intimate interaction between herbivorous insects and their host plants has an enormous influence on the dynamics and evolution of their populations. First, the ability to feed on plants has driven diversification and speciation processes in insects [1, 2]. Populations of oligophagous or polyphagous insects associated with host plants belonging to the same plant family or to different families, respectively [3], can become isolated into distinct lineages in response to spatial or temporal preferential use of hosts [4,5,6]. Should assortative mating of insects on the preferential host plant occur, pre-zygotic isolation may arise [7]. In short, the initial divergence due to host-plant preference can lead to speciation within insect lineages.

The study of ecological speciation – that is, adaptive divergence between populations due to ecological factors [7, 8] – encompasses the investigation of several mechanisms through which divergent selection can cause reproductive isolation [7, 9, 10]. The process of lineage differentiation is continuous and can lead to increasing reproductive isolation [5], designated as a “stage” of speciation by some authors [10].

How does the initial diversification arise in divergent populations that use alternative host plants? Recent studies have emphasized the role of phenotypic plasticity in diversification and speciation events, including numerous species of phytophagous insects that use different host plants [11,12,13,14]. If the initial preference for a new host plant is a plastic response, in the sense of the ability of a single-genotype organism to react to a novel environment and to produce different phenotypes [11, 13], the evolutionary mechanism behind the use of alternative host plants in the distributional range of an insect, leading to an increase in reproductive isolation, can follow the developmental-plasticity hypothesis proposed by West-Eberhard [12]. According to this model, adaptive selection of plastic phenotypes occurs in two steps: 1) a novel environmental factor affects plastic phenotypes, leading to novel variants; 2) novel variants affected by environmental recurrences of the initial stimulus are selected to produce evolutionary genetic change that can finally become fixed by genetic assimilation [13]. In this way, phenotypic plasticity can facilitate evolutionary change and speciation by enabling subsequent adaptation [13]. Accordingly, phenotypic plasticity allows diverse environments to perform a dual role in evolution, since it generates phenotypic variation that can potentially respond to selection. In addition, plasticity in one trait can influence the variation and selection in linked or correlated traits [11].

Another important matter related to divergent selection and the speciation process in the interactions between phytophagous insects and their host plants is whether adaptation is the product of many loci with small effects or of a few loci with a large effect [15]. A prevalent idea is that adaptation to a new environment, or to a host plant in the present scenario, would involve a few key loci with large effects on the organism genome [16, 17]. How then can we identify the key genes involved in the speciation process as a response to host-plant preference?

Our main purpose here was to identify plastic and evolved genes that are involved in ecological speciation toward primary and alternative host plants, attempting to answer how the host plants influence the differentiation of lineages with increasing reproductive isolation in polyphagous insect populations. The fall armyworm (FAW), the noctuid moth Spodoptera frugiperda (J. E. Smith), is an appropriate model to answer these questions. S. frugiperda is a polyphagous species and one of the most important pest in several crops worldwide, including corn (Zea mays L.), rice (Oryza sativa L.), and, more recently, cotton (Gossypium hirsutum L.). The most impressive feature of S. frugiperda, however, is its differentiation into host-plant-related strains, one that feeds preferentially on corn, cotton, and sorghum (corn strain, CS), and another that is found usually on rice and several pasture grasses (rice strain, RS) [18]. This distinction is indeed an ongoing process, as suggested by evidence of hybridization between the strains [19]. Several terms have been used in the literature to describe these strains, sometimes as synonyms in the same manuscript, such as biotypes [20], host strains [18, 19, 21,22,23,24,25,26], host races [19], host assemblages [26], ecological races [21], genetically differentiated forms [26], host forms [21, 27], genetic groups [19], and even sibling species [18, 21]. Numerous pieces of evidence confirm the differentiation of FAW into strains, including habitat, temporal and behavioral isolation, unidirectional mating (revised as pre- and post-zygotic barriers in [22] and in [19]). Although morphologically identical, the two strains possess at least two main reproductive isolation mechanisms that are responsible for their maintenance: differences in the composition of female pheromones [23], and in the period of reproductive activities [25]. Although there is no strong evidence of complete reproductive isolation between the strains [28], interbreeding between them results in loss of fertility [29].

The effect of host plants on the differentiation of these FAW strains has long been recognized [19], but how much the preferential use of a host plant contributes to reproductive isolation between the strains is still undetermined. Likewise, the genetic mechanism of adaptation and the preferential use of host plants in the speciation process of FAW are unknown. In agreement with the hypothesis that a few genes with key functions should be responsible for the adaptation to one host plant as opposed to another [16], few AFLP loci specific to each strain have been identified [20, 26].

Currently, the most commonly applied next-generation sequencing technique is transcriptome characterization, although few studies have used this approach to study the interactions between phytophagous insects and their host plants [30, 31]. Unbiased studies of the transcriptional-profile resulting from RNA-sequencing (RNA-Seq) are a powerful tool for studying the adaptation and speciation processes in such interactions, since they allow differentially regulated genes to be functionally characterized under diverse conditions. They are also the initial stage in understanding genetic, ecological and evolutionary mechanisms involved in the speciation process as a result of adaptive divergent selection in response to host-plant use [15, 31]. In an agricultural context, the mechanisms of insect adaptation to different alternative hosts offer a useful model for understanding how evolution may proceed in genetically engineered crops [32]. We applied RNA-Seq and transcriptome characterization to evaluate unbiased genetic-expression differences in larvae of the two strains of S. frugiperda fed on primary (corn) or alternative (rice) host plants. We consider corn as the primary FAW host plant because currently the species is more easily found on this host, while it is rarely found on rice in Brazil; however, the preference between these two host plants is considered ambiguous for the FAW [33].

According to the West-Eberhard model [12], the genes that underlie a plastic response to new environments, exhibiting differences in expression, can be the same genes that are differently regulated in recently diverged lineages [34]. Thus, we consider that genes that are differently regulated by the same FAW strain in response to exposure to different host plants, to be “plastic” in the sense that they are a reaction of the organisms to two distinct environments. Otherwise, differences in gene expression in the two strains fed on the same host plant are considered constitutive differences, and should be involved in the process of speciation between FAW strains.

We can consider that host-plant recognition and use-induced phenotypic changes involve multiple regulatory genes and processes through different hierarchies, as in other environment-induced phenotypic changes [13]. For this reason, an unbiased approach to estimate the gene expression is not only desirable, but also essential to better understand the phenotypic and genetic responses of phytophagous insects to preferential and alternative host plants, and how reproductive isolation evolves in these circumstances. Therefore, we first evaluated differences in the larval performance of the FAW strains fed on primary and alternative host plants under laboratory conditions. We then assembled and functionally characterized the unbiased transcriptome profile of S. frugiperda, with emphasis on differentially regulated genes in larvae reared in different conditions. Finally, we investigated the hypothesis that differentially regulated genes of the transcriptional plastic response to host plants are highly represented as constitutive differentially regulated genes in each FAW strain reared on different hosts.

Methods

Insect rearing

Populations of Spodoptera frugiperda for either CS or RS were kept on a white bean-based artificial diet [35] under laboratory conditions. The rice strain (RS) colony was originated from a cornfield collection of 170 larvae at Santa Helena de Goiás, Goiás, Brazil, in 2011-winter season. This population was separated in single pair mating in laboratory and the adults were genotyped to determine the strain. A single couple was genotyped as RS and the colony was established. The colony was maintained on artificial diet in mass mating cages for 2 years. In 2013, we collected insect samples from this population and established single-pair matings that were strain genotyped for the present study. The corn strain (CS) was obtained from a cornfield collection of 317 larvae at Campo Mourão, Paraná, Brazil in 2013-winter season. The population was kept in artificial diet and mass mating cages under laboratory conditions for ~ 5 generations. Then, we established single pair matings and genotyped both male and female to check the strain.

For each strain, one couple was placed in a cylindrical plastic cage (23 cm height × 10 cm diameter) for mating, and immediately after the female oviposited, the adults were removed from the cage to extract DNA for strain genotyping, as explained below. Experiments with both CS and RS were conducted with sibling larvae from the same family line.

One hundred twenty-eight to 144 neonate larvae of CS and RS (Table 1) were transferred to individual plates with a white bean-based artificial diet or with fresh leaves of one of the two host plants: i) corn (Z. mays), cultivar 2B688 Dow AgroSciences, hereafter referred to as the primary host plant; and ii) rice (O. sativa), cultivar Puitá INTA-CL, hereafter the alternative host plant. The six combinations were: 1) CS on diet; 2) RS on diet; 3) CS on corn leaves; 4) CS on rice leaves; 5) RS on corn leaves; and 6) RS on rice leaves. The leaves were replaced with fresh leaves every two days throughout the larval development period. The larvae were kept in climate-controlled chambers at 27 ± 1 °C, relative humidity of 60 ± 10% and a 14-h photophase. Part of the larvae were reared until pupation, and some 5th-instar larvae were stored in RNAlater® (Life Technologies, Carlsbad, CA, USA) for RNA extraction, in the proportion of 1:5 of mass:RNAlater (Table 1).

Table 1 Number of replicates in each larval stage during rearing on primary and alternative host plants
Full size table

Strain genotyping

Total genomic DNA was extracted from each couple (both male and female) presumed to be of either CS or RS. DNA was extracted using the DNeasy Blood & Tissue Kit (Qiagen, Dusseldorf, Germany) from the thoracic tissue of each individual. The mitochondrial gene cytochrome c oxidase I (COI) (ca. 569 bp) was amplified using the primers JM76 and JM77, under the same PCR program described in [36]. Amplification reactions were conducted with 25 μL total volume, using 1 μL of total DNA, 2 mM of 25 M MgCl2, 40 μM of dNTPs, 0.2 μM of each primer, 1 U of GoTaq DNA Polymerase (Promega, Madison, WI, USA), 10% of 10X Taq buffer and 10% volume of 5% dimethyl sulfoxide (DMSO). PCR products were run in 1% agarose gel in TAE 1X buffer (40 mM Tris, 20 mM acetic acid, 1 mM EDTA (pH 8.0)) to confirm amplification. After amplification, 1.0 μL of FastDigest MspI (Thermo Scientific, Waltham, MA, USA) was added to 10 μL of each reaction, incubated at 37 °C for 10 min, and the complete volume was loaded in 2% agarose gel in TAE buffer.

A repeated DNA sequence known as FR [37] was also amplified for each male and female to confirm their strains, using primers FR-c and FR-2, and conditions described elsewhere [36], using the same PCR procedure as above. PCR products were run in 2% agarose gel in TAE buffer to observe the band patterns related to each strain.

Larval performance

Larval performance in each rearing condition was evaluated by measuring the weight (in mg) of 10-day-old larvae and 24-h-old pupae, and the time (in days) to complete larval development (from first instar to pupa). Statistical differences were evaluated by log-transforming weight values and comparing mean differences among treatments, using the Tukey test [38].

RNA extraction, library preparation, and sequencing

Twelve larvae from each feeding condition in the early 5th instar were stored in RNAlater. Later, larvae were removed from the storage reagent, their gut contents were washed with 0.9% NaCl physiological solution to remove residual food, and the whole bodies were immediately immersed in liquid nitrogen and ground together. To increase the power of the post-sequencing statistical analyses and efficiently use sequencing resources [39, 40], three independent biological replicates for each condition, with 12 larvae each, were conducted using this procedure, totaling 36 individuals from each treatment. Individual replicates were stored at −80 °C until RNA extraction.

RNA was extracted using TRIzol™ (Life Technologies) combined with Direct-zol™ RNA MiniPrep (Zymo Research, Irvine, CA, USA). Each sample was eluted in 40 μL of Ultrapure water (Qiagen). RNA quality was evaluated using a Nanodrop UV spectrophotometer (Techno Scientific, Wilmington, DE, USA), and quantified in a Qubit® 2.0 Fluorometer (Life Technologies). RNA validation was completed in an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) in the Central Laboratory of High Performance Technologies (LaCTAD), University of Campinas, Brazil, where libraries were prepared using the Denaturing and Diluting Libraries for the HiSeq® and GAIIx (using 15 pM of sample per lane), TruSeq® RNA Sample Preparation v2, and cBot (Illumina, San Diego, CA, USA). Each of the three biological replicates per condition (library) was run in a different lane in an Illumina HiSeq 2500 System (Illumina), with 12 libraries per lane (Table 2).

Table 2 De novo assembly descriptive metrics
Full size table

Transcriptome assembly

Illumina reads from all replicates were processed using Illumina pipeline v. 1.8 or later. Prior to assembly, reads from all replicates were merged into a single dataset and the resulting FASTQ files (Illumina) were filtered to remove low-quality bases (score < 30) and adapters, using the SeqyClean pipeline (https://github.com/ibest/seqyclean). The remaining data were normalized by a single-pass digital normalization using the “normalize_by_kmer_coverage” procedure of the Trinity assembler suite [41]. A de novo hybrid assembly using ESTs data from the SpodoBase project (http://bioweb.ensam.inra.fr/spodobase/) and in-house Illumina sequencing reads was performed with the MIRA assembler v. 4.9.3 [42] with two passes and the following settings: i) spoiler detection on (−AS:sd = yes), ii) 70% relative percentage of exact word matches (−SK:pr = 70), iii) maximum megahub ratio = 1 (−SK:mmhr = 1), iv) stepping increment = 2 (−SK:kss = 2), v) quality clipping on (−CL:qc = yes), vi) minimum base quality = 5 (−CL:qcmq = 5), vii) length of window for quality clipping = 5 (−CL:qcwl = 5), and viii) elimination of sequences that form a contig with <3 reads (−AS:mrpc = 3). The final transcriptome was constructed by filtering the CAF (Common Assembly Format) file generated by MIRA for contigs larger than 200 bases, which shows more than five times coverage. Chimeric assembled contigs were identified and removed using a suite of customized Perl script (blast CHECK).

Annotation

Gene ontology (GO) annotation and sequence descriptions were indicated by a multi-step process, using a customized set of Perl scripts and local databases constructed with publicly available data (this database is available upon request to the authors). All assembled contigs were searched by similarity against NCBI REFSEQ [43] (updated on September 3, 2015) and MEROPS v. 9.12 [44], a specific database for peptidases, using an e-value cut-off of 10e−5 and HSP similarity threshold of 80%. Patterns of RNA families were indicated by HMMScan [45] using the RFAM database v. 12 [46]. Patterns of protein families from the PFAM database [47] were proposed by HMMScan using a set of translated peptides from candidate coding regions within the assembled transcriptome sequences indicated by transdecoder (https://transdecoder.github.io/).

The sequence description was achieved by integrating all database searches. Blast best hit results, from both databases previously described, were designated by a restrictive e-value and HSP similarity cut-off (1e−10 and 90% respectively) sorted by the latter; RNA and protein families from HMMScan were filtered by the Expectation Value (1e−10) in the “full sequence” column from the resulting analyses. Sequences that did not recover any information from all the databases used were tagged as “UNKNOWN DESCRIPTION”. The sequence description can be presented by all or any information from each of the databases. The gene ontology controlled vocabulary terms were assigned to all sequences by a custom Perl script that searches the GO local database (updated on September 5, 2015) using all sequence-similarity results from previous described searches, following the thresholds described above and removing all obsolete terms and redundant ontologies for the same sequence.

Differential regulation in S. frugiperda strains feeding on primary or alternative host plants

Differential expression analyses were carried out with scripts from the RSEM pipeline [48] to prepare the transcriptome reference index and to calculate the relative expression within each library. Differential gene expression analysis was conducted using edgeR [49] with TMM normalization [50], as suggested in the Trinity pipeline procedure [41], and the p-values were corrected for multiple testing by the false discovery rate [50]. The RPKM values for the most differentially regulated genes (corrected p-value <0.001 and log2 fold change >2 or < −2) were submitted to a MySQL database for faster consulting and comparisons among conditions. Differential gene regulation was compared in four tests: two comparisons were run to test the response of each strain to the same primary or alternative host plant (CS on corn vs. RS on corn, and CS on rice vs. RS on rice) (constitutive response), and two comparisons tested the response of the same strain to primary and alternative host plants (CS on corn vs. CS on rice, and RS on corn vs. RS on rice) (plastic response).

To investigate the hypothesis that differentially regulated genes of a transcriptional plastic response to host plants are also highly represented as differentially regulated evolved genes, we surveyed contigs that were simultaneously differentially regulated in both response types. First we constructed two lists of differentially regulated genes, the first composed of genes present in the two plastic comparisons, and the second composed of genes present in the two constitutive comparisons. Then we examined whether overregulated genes in the first list were also present in the second one. These contigs were gathered by a set of SQL scripts, consulting the previously generated database, and gene ontology annotation was summarized using REVIGO [51]. The 10 most enriched GO terms were identified by the weight01 method implemented in topGO [52].

All scripts used in this study are available upon request to the authors.

Results

Strain genotyping

The identity of each parental couple was confirmed as CS and RS by both the PCR-RFLP analysis of the mitochondrial COI gene and the FR repeated DNA sequence (data not shown).

Larval performance

Larval weight at 10 days was similar in both CS and RS when they were reared on corn leaves (Fig. 1). CS larvae were significantly heavier than RS larvae when they were fed on rice and on the artificial diet. However, the pupal weight was not significantly different for CS and RS fed on rice and on the artificial diet, but RS pupae were significantly heavier than CS pupae when they were fed on corn leaves. Although the larval and pupal weights did not differ significantly between the two strains, RS larvae required at least two more days to pupate than CS larvae when both were reared on rice leaves, indicating phenotypic differences between CS and RS as a response to the same host plant (Table 1). Larval development time ranged from 12 to 25 days, depending on feeding conditions (Table 1). Both CS and RS showed the longest larval development times when reared on the alternative host plant (rice leaves), and the shortest times on the primary host plant (corn leaves).

Fig. 1
figure 1

Larval performance described as 10-day-old larval weight (above) and 24-h-old pupal weight (below) in CS and RS reared on corn and rice leaves and on an artificial diet. * indicates significant difference (P < 0.05) between strains reared in the same conditions

Full size image

RNA sequencing

As expected, the number of reads varied among replicates. However, the quality of sequencing did not vary widely, with more than 80% of the reads with a quality index above Q30.

Transcriptome assembly and annotation

After assembly, 71,425 contigs, with a mean size of 937 nt, were annotated by similarity (Table 2), and 26,389 contigs (ca. 37%) resulted in an unknown annotation. The size of non-annotated contigs varied from 199 to 7950 nt (mean size 695 nt). Most of the annotation information (98.6%) came from the phylum Arthropoda, class Insecta. The order Lepidoptera represented the source of annotation for 98.9% of the contigs within Insecta, but other orders were represented as well in lower numbers, including Hymenoptera, Diptera, Hemiptera, Coleoptera, and Phythiraptera. Within the order Lepidoptera, the pyralid Amyelois transitella was the most represented species (33%), followed by Bombyx mori (31.8%) and Papilio xuthus (16.7%) (Fig. 2).

Fig. 2
figure 2

Contribution of species of Lepidoptera for annotation. The number of annotated contigs is noted on the bars

Full size image

Annotation by Gene Ontology resulted in 1038 different processes and 9558 unigenes, and 853 different molecular functions and 13,601 unigenes (the complete databank is available upon request). Metabolic process GOs were associated with 12% of annotated contigs, followed by oxidation-reduction process, with 8% (Fig. 3A). Sequences annotated as metabolic process included mainly glutathione transferases and UDP-glucosyltransferases and UDP-glucuronosyltransferases. Oxidation-reduction process sequences included mainly cytochrome P450, and many glucose dehydrogenase (PFAM: GMC oxidoreductase) and NADH dehydrogenase sequences.

Fig. 3
figure 3

Gene Ontology (GO) assignments for FAW transcriptome, showing the 20 most-represented GO categories, presented as biological processes (a) and molecular functions (b). The number of unigenes within each GO category is shown after the comma

Full size image

The most represented categories of GOs in molecular functions were metal-ion binding, nucleotide binding, oxidoreductase activity, and hydrolase activity (Fig. 3B). Unigenes annotated as metal-ion binding and oxidoreductase activity included several alcohol dehydrogenases, several kinds of cytochromes, including cytochrome b and cytochrome c, and many entries of cytochrome P450. Nucleotide binding included many unigenes annotated as serine/threonine protein kinase, as well as T-complex protein, ras-related protein, and multidrug resistance-associated protein, among many others. Unigenes annotated as hydrolase activity function included many peptidases (such as trypsin, esterase, serine protease, and proteasome), juvenile hormones, and serine/threonine-protein phosphatase, among others.

Differential regulation on primary and alternative host plants

Comparisons of gene expression between pairs of feeding conditions revealed 225 contigs that are differentially regulated (at least two-fold) in all conditions (Fig. 4). The comparisons between the two host strains feeding on primary or alternative host plants resulted in 3657 contigs that are differentially regulated, and we considered them constitutive responses to the host plant. The main results of the two comparisons designed to investigate evolved or constitutive response of the two strains to the same host plant are described below.

Fig. 4
figure 4

Number of at least two-fold differentially regulated genes in each comparison between pairs of feeding conditions

Full size image

1) RS-corn vs. CS-corn: A total of 2395 contigs were differentially expressed when comparing the transcriptome of the two host strains fed on corn (Fig. 4). The number of contigs differently expressed was highest in this pairwise comparison contrasted with all others (Fig. 4). Gene Ontology annotation for molecular function was assigned for 768 differentially regulated contigs; 482 contigs were up-regulated in CS and 286 were down-regulated in CS in relation to RS. The top five most differentially regulated contigs, both up- and down-regulated, include molecular functions such as metal-ion binding (GO:0046872), oxidoreductase activity (GO:0016491), nucleotide binding (GO:0000166), ATP binding (GO:0005524), and hydrolase activity (GO:0016787) (Fig. 5A). The most common unigenes annotated as metal ion-binding and oxidoreductase activity functions are cytochrome P450 (4d1-like, 4 g15, 9120, CYP4L6 and CYP6AE9), and many cytochrome b and c unigenes. Unigenes with oxidoreductase activity also included aldo-keto reductase proteins, among others. The aldo-keto reductase (AKR) superfamily comprises several enzymes that catalyze redox transformations involved in biosynthesis, intermediary metabolism, and detoxification. Substrates of the family include glucose, steroids, glycosylation end products, lipid peroxidation products, and environmental pollutants [53].

Fig. 5
figure 5

Those Gene Ontologies (GOs) are the 10 most-represented Gene Ontologies categories of up- and down-regulated contigs in each pairwise comparison presented as molecular function. a and b: constitutive response; c and d: plastic response. Up- and down-expression is a result of the comparison of the second condition in relation to the first condition. X-axis represents a different expression, and is presented on different scales for better visualization. Illustration: Dadi (www.ilustradoradadi.com)

Full size image

The most common unigenes with nucleotide and ATP binding functions included different subunits of the T-complex protein 1 and serine/threonine-protein kinases. Unigenes with hydrolase activity included mainly serine proteases (MEROPS: S01). Highly differentially regulated contigs also included other functions such as peptidase activity, monooxygenase activity, serine-type endopeptidase activity, and even an odorant-binding function, described as pheromone-binding protein.

2) RS-rice vs. CS-rice: A total of 1748 contigs were differentially regulated when both strains were fed on fresh rice leaves (Fig. 4). Six hundred fourteen (614) contigs were assigned to specific molecular functions, from which 290 were up-regulated in CS and 324 were down-regulated in CS in comparison to RS reared on rice. Again, the three most differentially regulated contigs, both up- and down-regulated, included molecular functions such as oxidoreductase activity, metal-ion binding, and hydrolase activity (Fig. 5B). The most common unigenes annotated as oxidoreductase activity and metal ion-binding functions are cytochrome P450 (e.g., 302A1, 4 g15, 9a20, CYP4L6, CYP6AE9, 6AB4, CYP332A1, CYP9A21), and many cytochrome b, acyl-CoA, aldo-keto reductase, sorbitol dehydrogenase, NADH dehydrogenase unigenes; many unigenes were annotated as zinc-finger protein under metal ion-binding function. Unigenes annotated as catalytic and ATP-binding functions were also up-regulated in RS, and include unigenes such as kynureninase, p270, serine hydroxylmethyltransferase, transketolase and C-1-tetrahydrofolate synthase (catalytic function), 26S protease regulatory, C-1-tetrahydrofolate synthase and T-complex protein 1 (ATP binding function). Down-regulated unigenes in CS (up-regulated in RS when reared on rice) include nucleotide binding and transferase activity functions; under nucleotide binding category the most common unigenes were annotated as 26S protease regulatory, ABC transporter, serine/threonine-protein kinase, T-complex protein 1, and tubulin chain. Transferase activity unigenes included mainly glutathione transferases, followed by UDP-glucosyltransferases and serine/threonine-protein kinase.

A total of 1730 contigs were differentially regulated as a response of the same strain to primary or alternative host plants, and were considered as plastic responses; 1404 contigs were differentially regulated in CS fed on corn and rice, and 1312 in RS fed on the two host plants (Fig. 4). The main results of the two comparisons designed to investigate the plastic response of each strain to the two host plants are described next:

3) CS-rice vs. CS-corn: A total of 381 contigs were appointed to particular molecular functions. Of these, 124 were up-regulated in CS fed on rice, and 257 were down-regulated in the same condition, in comparison to CS reared on corn. The three most differentially regulated functions, presented as up- and down-regulated, include oxidoreductase activity, metal-ion binding, and hydrolase activity (Fig. 5C). The most common unigenes annotated as oxidoreductase activity include several cytochrome P450 (e.g., 4 g15, 6AB4, 302A1, CYP6AE9, CYP4L6) and apoptosis-inducing factor. Metal ion-binding function also comprised cytochromes P450, and several zinc-finger proteins. Unigenes annotated as hydrolase activity comprised serine proteases. Up-regulated contigs also include functions such as transferase activity (mainly glutathione transferase) and structural constituent of cuticle (mainly cuticular protein RR-1 and RR-2 motifs). Down-regulated contigs also comprised nucleotide and ATP-binding functions (mostly multidrug resistance proteins).

4) RS-rice vs. RS-corn: A total of 420 contigs were assigned to specific molecular functions. Of these, 183 were up-regulated in RS reared on rice, and 237 were down-regulated in the same condition, in comparison to RS fed on corn. The top five most differentially up-regulated contigs include molecular functions such as metal-ion binding, structural constituent of cuticle, structural constituent of ribosome, hydrolase activity, and oxidoreductase activity (Fig. 5D). The five most differentially down-regulated contigs also include oxidoreductase activity, metal-ion binding, and hydrolase activity, and also comprised catalytic activity and nucleotide binding. The most common unigenes annotated as metal-ion binding and oxidoreductase activity include several cytochrome P450 (e.g., 9a20, 4 g15, 6AE8, CYP6AE9, CYP332A1, CYP4M5), zinc-finger proteins, sorbitol dehydrogenase, and malate dehydrogenase. Unigenes annotated as hydrolase activity comprise beta-N-acetylglucosaminidases, integument esterase, and alpha-amylases. Unigenes under the structural constituent of cuticle function include mainly cuticular proteins RR-1 and RR-2 motifs, and unigenes under the structural constituent of ribosome function comprise mainly the 60S ribosomal protein. Finally, down-regulated unigenes under catalytic-activity functions include kynurenine/alpha-aminoadipate aminotransferase and malate dehydrogenase, and those under nucleotide-binding function contain mainly a tubulin beta chain.

Functional-ontology groups enrichment of genes that simultaneously show plastic and constitutive responses as differentially regulated contigs indicated that ca. 21% (1541 in 7219) of the plastic-response genes were also present as evolved-response contigs. These contigs include several different molecular functions (Fig. 6), and the 10 most-enriched genes in both types of comparisons are involved in oxidoreductase activity, constitution of cuticle, transport of lipids and transcription, structural integrity of a cytoskeletal structure, amino-acid attachment and synthesis of peptides, pyruvate kinase activity, and polymerase activity (Table 3). GO terms identified as significant terms have considerably large numbers of annotated genes, which gives additional confidence in the significance of our results.

Fig. 6
figure 6

Gene Ontology (GO) annotation summarized using REVIGO of simultaneously plastic and constitutive differentially regulated genes as a response to primary and alternative host plants

Full size image
Table 3 Ten most-enriched GO terms identified by the weight01 method implemented in topGO. Annotaded: number of genes annotated under the GO term; Significant: GOs identified as significant terms; Expected: number of interesting genes mapped to the GO term if the interesting genes were randomly distributed over all GO terms; Classic Fisher: test statistics based on gene counts [52]
Full size table

Discussion

The study of the speciation process in Spodoptera frugiperda strains that exploit different hosts is valuable for many reasons, one of them worthy of special mention: pest-insect problems are also evolutionary problems, in the sense that it is evolution that creates genetic modification in pest populations [32]. In doing so, pest-insects can offer simpler models of study for understanding evolution and speciation in herbivorous insects.

The differences in development time of FAW host strains reared on corn and rice found here were also found in other studies, in which RS larvae required longer to pupate than CS larvae on corn [54]. Groot and collaborators [19], however, summarized several performance experiments with S. frugiperda reared on corn and rice, which gave variable results. In general, RS is recorded as outperforming CS on rice (under the same parameters that we measured here), while CS tended to perform better than RS when reared on corn. In our experiment, both strains showed poorer performance when fed on rice, mainly in relation to the time needed for pupation. In field conditions, taking longer to pupate can mean that the larvae are more susceptible to predation and parasitism [55], and the FAW would have an ecological advantage when feeding on corn in relation to rice. In general, performance on a new or alternative host can be lower than the performance on the native host plant [34]. Our results suggest that corn is currently the preferred host plant of S. frugiperda, although both rice [27] and dicots [33] have been suggested as likely ancestral host plants of the FAW.

Although we found several plastic-response genes that were also constitutive-response genes, we were not able to show that the first kind of response is overrepresented in relation to the second, as found by other studies [34]. This result could be a consequence of the high number of differentially regulated genes that we found in each rearing comparison, due to the unbiased approach that we used. Despite this, a substantial proportion of contigs are simultaneously overregulated within the same host strain feeding on different host plants and in the two host strains feeding on the same host, and we can consider that the molecular functions that they have in this insect-plant system are important as both kinds of responses. The presence of the same genes as a plastic and evolved response, as we defined here, agrees with the model proposed by West-Eberhard [12] in which the use of a new environment would arise through a plastic response, and would be eventually selected and fixed in the population by repetition of the initial stimulus and by genetic assimilation.

In summary, three molecular functions were present in all four comparisons between pairs of feeding conditions, both down- and up-regulated: oxidoreductase activity, metal-ion binding, and hydrolase activity, as also found for other insect-plant systems [30]. Oxidoreductase activity and metal-ion binding include predominantly unigenes of several kinds of cytochrome P450, involved mainly in xenobiotics metabolization [56] and zinc-finger proteins. Unigenes under hydrolase activity comprise serine proteases and trypsins, involved in protein digestion in insects, as a response to plant protease inhibitors [31, 57, 58], and sometimes used as an anti-herbivore bypass mechanism. Hydrolase activity also includes fumarylacetoacetate hydrolases, proteasome subunits, and venom serine carboxypeptidase, which has the ability to release amino acids from the C-terminus of a peptide chain, can be employed in the determination of amino-acid sequences [59], and is used as a digestive enzyme by some insects [60].

Generalist herbivores have a range of detoxifying enzymes that enable them to feed on a diversity of available host plants [31]. These enzymes include many of the unigenes described as up- and down-regulated in our comparisons, such as cytochrome P450 monooxygenases (P450s), glutathione transferases, and UDP-glycosyltransferases.

P450 enzymes were the most commonly up- and down-regulated unigenes in all of our comparisons between conditions, both as plastic and as evolved responses by the FAW strains to their primary and alternative host plants, as found in other studies [54, 61,62,63]. In a very comprehensive review of P450s, Feyereisen [64] stated that most transcriptomics studies have shown that one or more P450s genes were differentially regulated, but in addition to being a direct relationship, this may be a cascade effect. Cytochrome P450, or CYP genes, comprises one of the largest families of genes. These enzymes have a monooxygenase function, catalyzing the transfer of one atom of molecular oxygen to a substrate and reducing the other atom to water, and also show several other catalytic activities [64]. An important role of P450s in insects is related to the detoxification of xenobiotics, and P450s are involved in many events of resistance to insecticides. P450 gene expression is regulated by chemicals that allow insects to respond to new conditions directly by building a detoxification defense, and indirectly by adapting their basal metabolism, including the hormone balance, rate of development and reproduction [64]. New hosts that contain new, and possibly toxic, compounds can plastically inhibit or enhance P450 gene expression, as we found here. Among several other functions performed by P450s [65], an important function related to reproductive isolation is the role of P450 enzymes in the biosynthesis of many insect pheromones and allomones, and evidence is accumulating for pheromone catabolism by P450 enzymes [64].

Serine proteases, also highly differentiated in our comparisons, are mostly related to digestion in insects, since they produce abundant proteases for the digestion of dietary proteins [58]. Many lepidopteran larvae use serine proteases for protein digestion [31], and other studies have found that the expression of these proteases depends on the host plant [30]. These digestive enzymes allow insects to overcome plant protein-inhibitor defenses, either by overexpressing existing proteases that are not a target of the inhibitors, or by expressing new ones [58, 66]. Serine proteases, however, are associated with several non-digestive functions in insects, and often function in cascade pathways; this can occur in insect embryonic development and immune responses [67]. Serine proteases can also have a physiological function in defense against infection, when they are present in the hemolymph [67].

Plant secondary metabolites can trigger these metabolic and/or digestive responses in insects. Corn, the primary host plant tested here, contains benzoxazinoids (BXDs) or hydroxamic acids, one of the main secondary compounds in many grasses. BXDs confer resistance on herbivorous insects and pathogens because of their antifeeding, insecticidal, antimicrobial and allelopathic activities [68,69,70]. These compounds are present in several cereal crops such as corn, wheat, and rye, but are absent in rice, cultivated barley and oats. DIMBOA (2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one) is the major BXD in the aerial parts of corn [71]. Although its mode of action is not fully elucidated, this hydroxamic acid inhibits digestive proteases in the midgut of larvae of the lepidopteran Ostrinia nubilalis (Pyralidae) [72] and Sesamia nonagrioides (Noctuidae) [73], suggesting that it acts as a digestive toxin [74]. When 5th instar larvae of O. nubilalis were fed on leaves of corn, several detoxification enzymes in the midgut, such as cytochrome b5, NADPH-cytochrome c reductase, NADPH oxidase and O-demethylase, increased activity [75]; while in in vitro tests with S. nonagrioides, esterases and glutathione transferases were strongly inhibited by DIMBOA [73]. Additionally, larvae of Ostrinia furnacalis fed on cabbage dipped in DIMBOA showed increases in the activity of cytochrome P450 monooxygenase and glutathione transferase [76]. Interestingly, although DIMBOA has a feeding-deterrent and/or toxic action for several lepidopteran pest species [68, 73], it acts as a feeding stimulant for S. frugiperda, and enhances FAW larval growth at low concentrations [69, 77]. All things considered, transcriptional differences in the responses of CS and RS to corn and rice leaves can be related to the presence of DIMBOA in corn, and/or to its absence in rice. Additional work with feeding of FAW with purified BXDs from corn seedlings added to an artificial diet could help to clarify this hypothesis.

Conclusions

As expected, differentially regulated genes in all feeding comparisons involved multiple regulatory genes and processes [13]. Using a much more informative molecular approach to evaluate the unbiased transcriptome profile of S. frugiperda, we were able to show that phenotypic plasticity and subsequent selection in response to alternative host plants is the product of actions of many loci, with diverse molecular functions through different hierarchies, and small individual effects. Metabolism, however, is suggested as the most important function, and the variable regulation of this molecular function indicates that metabolization of foreign chemicals is among the key players in the phenotypic variation in FAW strains. From an agricultural perspective, high plasticity in detoxifying-gene families indicates the possibility of a rapid response to control substances such as insecticides.

Abbreviations

AKR:

Aldo-keto Reductase

BXDs:

Benzoxazinoids

CAF:

Common Assembly Format

COI:

Cytochrome c Oxidase Subunit I

CS:

Corn Strain

CYP:

Cytochrome P450

DIMBOA:

2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one

EST:

Expressed Sequence Tag

FAW:

Fall Armyworm

GO:

Gene Ontology

HSP:

High-scoring Segment Pair

LaCTAD:

Central Laboratory of Higher Performance Technologies

NaCl:

Sodium Chloride

PFAM:

Protein families

RFAM:

RNA families

RNA-Seq:

RNA-sequencing

RPKM:

Reads Per Kilobase Million

RS:

Rice Strain

SQL:

Structured Query Language

TAE:

Tris-acetate-EDTA buffer

TMM:

Trimmed Mean of M-values

References

  1. Mitter C, Farrell B, Wiegmann B. The phylogenetic study of adaptive zones: has phytophagy promoted insect diversification? Am Nat. 1988;132:107–28.

    Article  Google Scholar 

  2. Mopper S. Adaptive genetic structure in phytophagous insect populations. Trends Ecol Evol. 1996;11:235–8.

    Article  CAS  PubMed  Google Scholar 

  3. Schoonhoven LM, van Loon JJA, Dicke M. Insect-plant biology. New York: Oxford University Press; 2008.

    Google Scholar 

  4. Machado V, Wunder M, Baldissera VD, Oliveira JV, Fiuza LM, Nagoshi RN. Molecular characterization of host strains of Spodoptera frugiperda (Lepidoptera: Noctuidae) in southern Brazil. Ann Entomol Soc Am. 2008;101:619–26.

    Article  CAS  Google Scholar 

  5. Dres M, Mallet J. Host races in plant-feeding insects and their importance in sympatric speciation. P Roy Soc B-Biol Sci. 2002;357:471–92.

    Google Scholar 

  6. Via S. Sympatric speciation in animals: the ugly duckling grows up. Trends Ecol Evol. 2001;16:381–90.

    Article  CAS  PubMed  Google Scholar 

  7. Rundle HD, Nosil P. Ecological speciation. Ecol Lett. 2005;8:336–52.

    Article  Google Scholar 

  8. Matsubayashi KW, Ohshima I, Nosil P. Ecological speciation in phytophagous insects. Entomol Exp Appl. 2010;134:1–27.

    Article  Google Scholar 

  9. Funk DJ, Nosil P. Comparative analyses of ecological speciation. In: Tilmon KJ, editor. Specialization, speciation, and radiation. Berkeley and Los Angeles: University of California Press; 2008. p. 117–35.

    Google Scholar 

  10. Nosil P, Harmon LJ, Seehausen O. Ecological explanations for (incomplete) speciation. Trends Ecol Evol. 2009;24:145–56.

    Article  PubMed  Google Scholar 

  11. Pfennig DW, Wund MA, Snell-Rood EC, Cruickshank T, Schlichting CD, Moczek AP. Phenotypic plasticity's impacts on diversification and speciation. Trends Ecol Evol. 2010;25:459–67.

    Article  PubMed  Google Scholar 

  12. West-Eberhard MJ. Developmental plasticity and evolution. New York: Oxford University Press; 2003.

    Google Scholar 

  13. Whitman DW, Agrawal AA. What is phenotypic plasticity and why is it important? In: Phenotypic Plasticity of Insects: Mechanisms and Consequences. Edited by Whitman DW, Ananthakrishnan TN: CRC press; 2009. p. 1–63.

    Chapter  Google Scholar 

  14. Moczek AP, Sultan S, Foster S, Ledon-Rettig C, Dworkin I, Nijhout HF, Abouheif E, Pfennig DW. The role of developmental plasticity in evolutionary innovation. P Roy Soc B-Biol Sci. 2011;278:2705–13.

    Article  Google Scholar 

  15. Stapley J, Reger J, Feulner PGD, Smadja C, Galindo J, Ekblom R, Bennison C, Ball AD, Beckerman AP, Slate J. Adaptation genomics: the next generation. Trends Ecol Evol. 2010;25:705–12.

    Article  PubMed  Google Scholar 

  16. Beaumont MA, Balding DJ. Identifying adaptive genetic divergence among populations from genome scans. Mol Ecol. 2004;13:969–80.

    Article  CAS  PubMed  Google Scholar 

  17. Orr HA. The genetic theory of adaptation: a brief history. Nat Rev Genet. 2005;6:119–27.

    Article  CAS  PubMed  Google Scholar 

  18. Pashley DP. Host-associated genetic differentiation in fall armyworm (Lepidoptera: Noctuidae): a sibling species complex? Ann Entomol Soc Am. 1986;79:898–904.

    Article  Google Scholar 

  19. Groot AT, Marr M, Heckel DG, Schofl G. The roles and interactions of reproductive isolation mechanisms in fall armyworm (Lepidoptera: Noctuidae) host strains. Ecol Entomol. 2010;35:105–18.

    Article  Google Scholar 

  20. Busato GR, Grutzmacher AD, de Oliveira AC, Vieira EA, Zimmer PD, Kopp MM, Bandeira JD, Magalhães TR. Analysis of the molecular structure and diversity of Spodoptera frugiperda (JE smith) (Lepidoptera: Noctuidae) populations associated to the corn and rice crops in Rio Grande do Sul state, Brazil. Neotrop Entomol. 2004;33:709–16.

    Article  CAS  Google Scholar 

  21. Dumas P, Barbut J, Le Ru B, Silvain JF, Clamens AL, d'Alencon E, Kergoat GJ. Phylogenetic molecular species delimitations unravel potential new species in the pest genus Spodoptera Guenee, 1852 (Lepidoptera, Noctuidae). PLoS One. 2015;10:e0122407. doi:10.1371/journal.pone.0122407

  22. Dumas P, Legeai F, Lemaitre C, Scaon E, Orsucci M, Labadie K, Gimenez S, Clamens AL, Henri H, Vavre F, et al. Spodoptera frugiperda (Lepidoptera: Noctuidae) host-plant variants: two host strains or two distinct species? Genetica. 2015;143:305–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Groot AT, Marr M, Schofl G, Lorenz S, Svatos A, Heckel DG. Host strain specific sex pheromone variation in Spodoptera frugiperda. Front Zool. 2008;5 doi:10.1186/1742–9994–1185-1120.

  24. Juarez ML, Murua MG, Garcia MG, Ontivero M, Vera MT, Vilardi JC, Groot AT, Castagnaro AP, Gastaminza G, Willink E. Host association of Spodoptera frugiperda (Lepidoptera: Noctuidae) corn and rice strains in Argentina, Brazil, and Paraguay. J Econ Entomol. 2012;105:573–82.

    Article  PubMed  Google Scholar 

  25. Pashley DP, Hammond AM, Hardy TN. Reproductive isolating mechanisms in fall armyworm host strains (Lepidoptera, Noctuidae). Ann Entomol Soc Am. 1992;85:400–5.

    Article  Google Scholar 

  26. Prowell DP, McMichael M, Silvain JF. Multilocus genetic analysis of host use, introgression, and speciation in host strains of fall armyworm (Lepidoptera: Noctuidae). Ann Entomol Soc Am. 2004;97:1034–44.

    Article  CAS  Google Scholar 

  27. Juarez ML, Schofl G, Vera MT, Vilardi JC, Murua MG, Willink E, Hanniger S, Heckel DG, Groot AT. Population structure of Spodoptera frugiperda maize and rice host forms in South America: are they host strains? Entomol Exp Appl. 2014;152:182–99.

    Article  CAS  Google Scholar 

  28. Pashley DP. Causes of host-associated variation in insect herbivores: an example from fall armyworm. In: Kim KC, York MPBAN, editors. Evolution of insect pests: patterns of variation. Inc.: John Wiley and Sons; 1993. p. 351–9.

    Google Scholar 

  29. Busato GR, Loeck AE, Garcia MS, Bernardi O, Zart M, Nunes AM, Zazycki LCF. Compatibilidade reprodutiva entre os biótipos "milho" e "arroz" de Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae). Rev Bras Agro. 2008;14:273–8.

    Google Scholar 

  30. Eyres I, Jaquiery J, Sugio A, Duvaux L, Gharbi K, Zhou JJ, Legeai F, Nelson M, Simon JC, Smadja CM, et al. Differential gene expression according to race and host plant in the pea aphid. Mol Ecol. 2016;25:4197–215.

    Article  CAS  PubMed  Google Scholar 

  31. Vogel H, Musser RO. Celorio-Mancera MdlP. Transcriptome responses in herbivorous insects towards host plant and toxin feeding. Annual Plant Reviews. 2014;47:197–234.

    Article  CAS  Google Scholar 

  32. Via S. Ecological genetics and host adaptation in herbivorous insects - the experimental study of evolution in natural and agricultural systems. Annu Rev Entomol. 1990;35:421–46.

    Article  CAS  PubMed  Google Scholar 

  33. Kergoat GJ, Prowell DP, Le Ru BP, Mitchell A, Dumas P, Clamens AL, Condamine FL, Silvain JF. Disentangling dispersal, vicariance and adaptive radiation patterns: a case study using armyworms in the pest genus Spodoptera (Lepidoptera: Noctuidae). Mol Phylogenet Evol. 2012;65:855–70.

    Article  PubMed  Google Scholar 

  34. Ragland GJ, Almskaar K, Vertacnik KL, Gough HM, Feder JL, Hahn DA, Schwarz D. Differences in performance and transcriptome-wide gene expression associated with Rhagoletis (Diptera: Tephritidae) larvae feeding in alternate host fruit environments. Mol Ecol. 2015;24:2759–76.

    Article  CAS  PubMed  Google Scholar 

  35. Kasten Jr. P, Precetti AACM, Parra JRP. Dados biológicos comparativos de Spodoptera frugiperda (J.E. Smith, 1797) em duas dietas artificiais e substrato natural. Rev Agri 1978;53:68–78.

  36. Nagoshi RN, Meagher R. Fall armyworm FR sequences map to sex chromosomes and their distribution in the wild indicate limitations in interstrain mating. Insect Mol Biol. 2003;12:453–8.

    Article  CAS  PubMed  Google Scholar 

  37. Y-J L, Kochert GD, Isenhourt DJ, Adang MJ. Molecular characterization of a strain-specific repeated DNA sequence in the fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae). Insect Mol Biol. 1994;3:123–30.

    Article  Google Scholar 

  38. Zar JH. Biostatistical analysis. New Jersey: Pearson; 2010.

    Google Scholar 

  39. Hansen KD, Wu Z, Irizarry RA, Leek JT. Sequencing technology does not eliminate biological variability. Nat Biotechnol. 2011;29:572–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Conesa A, Madrigal P, Tarazona S, Gomez-Cabrero D, Cervera A, McPherson A, Szczesniak MW, Gaffney DJ, Elo LL, Zhang X, et al. A survey of best practices for RNA-seq data analysis. Genome Biol. 2016;17:13.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, Couger MB, Eccles D, Li B, Lieber M, et al. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.

    Article  CAS  PubMed  Google Scholar 

  42. Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Muller WEG, Wetter T, Suhai S. Using the miraEST assembler for reliable and automated mRNA transcript assembly and SNP detection in sequenced ESTs. Genome Res. 2004;14:1147–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Pruitt KD, Tatusova T, Maglott DRNCBI. Reference sequences (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007;35:D61–5.

    Article  CAS  PubMed  Google Scholar 

  44. Rawlings ND, Waller M, Barrett AJ, Bateman AMEROPS. The database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res. 2014;42:D503–9.

    Article  CAS  PubMed  Google Scholar 

  45. Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011;7:e1002195.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Nawrocki EP, Burge SW, Bateman A, Daub J, Eberhardt RY, Eddy SR, Floden EW, Gardner PP, Jones TA, Tate J, et al. Rfam 12.0: updates to the RNA families database. Nucleic Acids Res. 2015;43:D130–7.

    Article  CAS  PubMed  Google Scholar 

  47. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, Heger A, Hetherington K, Holm L, Mistry J, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42:D222–30.

    Article  CAS  PubMed  Google Scholar 

  48. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12: doi:10.1186/1471–2105–1112-1323.

  49. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.

    Article  CAS  PubMed  Google Scholar 

  50. Kadota K, Nishiyama T, Shimizu K. A normalization strategy for comparing tag count data. Algorithm Mol Biol. 2012;7:5. doi:10.1186/1748-7188-1187-1185.

  51. Supek F, Bosnjak M, Skunca N, Smuc TREVIGO. Summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6:e21800.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Alexa A, Rahnenfuhrer J, Lengauer T. Improved scoring of functional groups from gene expression data by decorrelating GO graph structure. Bioinformatics. 2006;22:1600–7.

    Article  CAS  PubMed  Google Scholar 

  53. Barski OA, Tipparaju SM, Bhatnagar A. The aldo-keto reductase superfamily and its role in drug metabolism and detoxification. Drug Metab Rev. 2008;40:553–624.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Roy A, Walker WB, Vogel H, Chattington S, Larsson MC, Anderson P, Heckel DG, Schlyter F. Diet dependent metabolic responses in three generalist insect herbivores Spodoptera spp. Insect Biochem Molec. 2016;71:91–105.

    Article  CAS  Google Scholar 

  55. Lima S, Dill LM. Behavioral decisions made under the risk of predation: a review and prospectus. Can J Zool. 1990;68:619–40.

    Article  Google Scholar 

  56. Heidel-Fischer HM, Vogel H. Molecular mechanisms of insect adaptation to plant secondary compounds. Curr Opin Insect Sci. 2015;8:8–14.

    Article  Google Scholar 

  57. Terra WR, Ferreira C, Jordão BP, Dillon RJ. Digestive Enzymes. In: Biology of the insect midgut. Edited by Lehane MJ, Billingsley PF. London: Chapman & Hall; 1996.

  58. Terra WR, Ferreira C. Biochemistry and molecular biology of digestion. In: Insect Molecular Biology and Biochemistry. Edited by Gilbert LI. San Diego: academic press; 2012. p. 365–418.

    Chapter  Google Scholar 

  59. Breddam K. Serine carboxypeptidases. A review. Carlberg research. Communication. 1986;51:83–128.

    CAS  Google Scholar 

  60. Waniek PJ, Araújo CAC, Momoli MM, Azambuja P, Hansen AM, Genta FA. Serine carboxypeptidases of Triatoma brasiliensis (Hemiptera, Reduviidae): sequence characterization, expression pattern and activity localization. J Insect Physiol. 2014;63:9–20.

    Article  CAS  PubMed  Google Scholar 

  61. Celorio-Mancera MD, Heckel DG, Vogel H. Transcriptional analysis of physiological pathways in a generalist herbivore: responses to different host plants and plant structures by the cotton bollworm, Helicoverpa armigera. Entomol Exp Appl. 2012;144:123–33.

    Article  Google Scholar 

  62. Koenig C, Bretschneider A, Heckel DG, Grosse-Wilde E, Hansson BS, Vogel H. The plastic response of Manduca sexta to host and non-host plants. Insect Biochem Molec. 2015;63:72–85.

    Article  CAS  Google Scholar 

  63. Perera OP, Shelby KS, Popham HJR, Gould F, Adang MJ, Jurat-Fuentes JL. Generation of a transcriptome in a model lepidopteran pest, Heliothis virescens, using multiple sequencing strategies for profiling midgut gene expression. PLoS One. 2015;10:e0128563. doi:10.1371/journal.pone.0128563

  64. Feyereisen R, Insect CYP. genes and P450 enzymes. In: Insect Molecular Biology and Biochemistry. Edited by Gilbert LI. San Diego: academic press; 2012. p. 236–316.

    Chapter  Google Scholar 

  65. Feyereisen R. Insect cytochrome P450. In: Science CMI, editor. Gilbert LE, Iatrou K, gill SS, vol. 4. New York: Elsevier Pergamon; 2005. p. 1–77.

    Google Scholar 

  66. Souza TP, Dias RO, Castelhano EC, Brandão MM, Moura DS, Silva-Filho MC. Comparative analysis of expression profiling of the trypsin and chymotrypsin genes from Lepidoptera species with different levels of sensitivity to soybean peptidase inhibitors. Comp Biochem Physiol B: Biochem Mol Biol. 2016;196–197:67–73.

    Article  Google Scholar 

  67. Kanost MR, Clem RJ. Insect proteases. In: Insect Molecular Biology and Biochemistry. Edited by Gilbert LI. San Diego: academic press; 2012. p. 347–64.

    Google Scholar 

  68. Klun JA, Tipton CL, Brindley TA. 2,4-Dihydroxy-7-Methoxy-1,4-Benzoxazin-3-1 (Dimboa) an active agent in resistance of maize to European corn borer. J Econ Entomol. 1967;60:1529–33.

    Article  CAS  Google Scholar 

  69. Rostas M. The effects of 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one on two species of Spodoptera and the growth of Setosphaeria turcica in vitro. J Pest Sci. 2007;80:35–41.

    Article  Google Scholar 

  70. Wouters FC, Gershenzon J, Vassão DG. Benzoxazinoids: reactivity and modes of action of a versatile class of plant chemical defenses. J Braz Chem Soc. 2016;27:1379–97.

    CAS  Google Scholar 

  71. Niemeyer HM. Hydroxamic acids (4-Hydroxy-1,4-Benzoxazin-3-ones), defense chemicals in the Gramineae. Phytochemistry. 1988;27:3349–58.

    Article  CAS  Google Scholar 

  72. Houseman JG, Campos F, Thie NMR, Philogene BJR, Atkinson J, Morand P, Arnason JT. Effect of the maize-derived compounds DIMBOA and MBOA on growth and digestive processes of European corn borer (Lepidoptera, Pyralidae). J Econ Entomol. 1992;85:669–74.

    Article  CAS  Google Scholar 

  73. Ortego F, Ruiz M, Castanera P. Effect of DIMBOA on growth and digestive physiology of Sesamia nonagrioides (Lepidoptera : Noctuidae) larvae. J Insect Physiol. 1998;44:95–101.

    Article  CAS  PubMed  Google Scholar 

  74. Niemeyer HM. Hydroxamic acids derived from 2-hydroxy-2H-1,4-benzoxazin-3(4H)-one: key defense chemicals of cereals. J Agri Food Chem. 2009;57:1677–96.

    Article  CAS  Google Scholar 

  75. Feng R, Houseman JG, Downe AER, Atkinson J, Arnason JT. Effects of 2,4-Dihydroxy-7-methoxy-l,4-benzoxazin-3-one (DIMBOA) and 6-Methoxybenzoxazolinone (MBOA) on the detoxification processes in the larval midgut of the European corn borer. Pestic Biochem Phys. 1992;44:147–54.

    Article  CAS  Google Scholar 

  76. Yan F, Xu C, Li S, Lin C, Li J. Effects of DIMBOA on several enzymatic systems in Asian corn borer, Ostrinia furnacalis (Guenee). J Chem Ecol. 1995;21:2047–56.

    Article  CAS  PubMed  Google Scholar 

  77. Glauser G, Marti G, Villard N, Doyen GA, Wolfender JL, Turlings TCJ, Erb M. Induction and detoxification of maize 1,4-benzoxazin-3-ones by insect herbivores. Plant J. 2011;68:901–11.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank the editors and the reviewers, in special Dr. Astrid Groot, for their valuable comments, which helped us to improve the manuscript.

Funding

This work was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (Fapesp, process 2012/16266–7). Analyzes were performed in the Thunder Computer Cluster funded by Fapesp (process 2011/00417–3). AF and CO are recipients of Brazilian CNPq fellowships.

Availability of data and materials

RNA sequencing data are available at Bioproject (PRJNA408280) and Biosample (SAMN07678153). All scripts generated during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

  1. Laboratório de Melhoramento de Plantas, Centro de Energia Nuclear na Agricultura, Universidade de São Paulo, Campus “Luiz de Queiroz”, Av. Centenário, 303, Piracicaba, SP, 13400-970, Brazil

    Karina Lucas Silva-Brandão & Antonio Figueira

  2. Present address: Centro de Biologia Molecular e Engenharia Genética, Universidade Estadual de Campinas, Av. Cândido Rondon, 400, Campinas, SP, 13083-875, Brazil

    Karina Lucas Silva-Brandão

  3. Programa de Pós-graduação em Entomologia, Departamento de Entomologia e Acarologia, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Av. Pádua Dias, 11, Piracicaba, SP, 13418-900, Brazil

    Renato Jun Horikoshi & Daniel Bernardi

  4. Departamento de Entomologia e Acarologia, Escola Superior de Agricultura “Luiz de Queiroz”, Universidade de São Paulo, Av. Pádua Dias, 11, Piracicaba, SP, 13418-900, Brazil

    Celso Omoto

  5. Centro de Biologia Molecular e Engenharia Genética, Universidade Estadual de Campinas, Av. Cândido Rondon, 400, Campinas, SP, 13083-875, Brazil

    Marcelo Mendes Brandão

Authors
  1. Karina Lucas Silva-Brandão
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Renato Jun Horikoshi
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Daniel Bernardi
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Celso Omoto
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Antonio Figueira
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Marcelo Mendes Brandão
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

KLSB idealized the study, completed lab work, interpreted the data, prepared figures and wrote the manuscript; MMB analyzed and annotated the transcriptome, and was a major contributor in idealizing the study and in writing the manuscript; RJH and DB reared all insects and completed larval performance studies; CO and AF helped idealizing the study, in discussing the results and in writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Karina Lucas Silva-Brandão.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Silva-Brandão, K.L., Horikoshi, R.J., Bernardi, D. et al. Transcript expression plasticity as a response to alternative larval host plants in the speciation process of corn and rice strains of Spodoptera frugiperda . BMC Genomics 18, 792 (2017). https://doi.org/10.1186/s12864-017-4170-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-017-4170-z

Keywords

  • Digestive enzyme
  • Ecological speciation
  • Fall armyworm
  • P450
  • Transcriptome

BMC Genomics

ISSN: 1471-2164

Contact us