Determinants of Insecticide Resistance Evolution: Comparative Analysis Among Heliothines.

It is increasingly clear that pest species vary widely in their propensities to develop insecticide resistance. This review uses a comparative approach to analyze the key pest management practices and ecological and biochemical or genetic characteristics of the target that contribute to this variation. We focus on six heliothine species, three of which, Helicoverpa armigera, Heliothis virescens, and Helicoverpa zea, have developed resistances to many pesticide classes. The three others, Helicoverpa punctigera, Helicoverpa assulta, and Helicoverpa gelotopoeon, also significant pests, have developed resistance to very few pesticide classes. We find that host range and movement between alternate hosts are key ecological traits that influence effective selection intensities for resistance. Operational issues are also critical; area-wide, cross-pesticide management practices that account for these ecological factors are key to reducing selection intensity. Without such management, treatment using broad-spectrum chemicals serves to multiply the effects of host plant preference, preadaptive detoxification ability, and high genetic diversity to create a pesticide treadmill for the three high-propensity species.Without rigorous ongoing management, such a treadmill could still develop for newer, more selective chemistries and insecticidal transgenic crops.

A Generalist Herbivore Copes with Specialized Plant Defence: the Effects of Induction and Feeding by Helicoverpa armigera (Lepidoptera: Noctuidae) Larvae on Intact Arabidopsis thaliana (Brassicales) Plants.

Plants of the Brassicaceae are defended from feeding by generalist insects by constitutively-expressed and herbivory-induced glucosinolates (GS). We induced Arabidopsis plants 1, 16 and 24 h prior to allowing neonate larvae of the generalist Helicoverpa armigera to feed on whole plants for 72 h. These plants were subsequently retested with another group of neonates for a further 72 h. We used wild-type A. thaliana Col-0, and mutant lines lacking indolic GS, aliphatic GS or all GS. We hypothesized that larvae would not grow well on defended plants (WT) compared to those lacking GS, and would not grow well if plants had been primed or fed on for longer, due to the expected induced GS. There was survivorship on all lines suggesting H. armigera is a suitable generalist for these experiments. Larvae performed less well on wild-type and no indolic lines than on no aliphatic and no GS lines. Larvae distributed feeding damage extensively in all lines, more so on wild type and no-indolic lines. Contrary to expectations, larvae grew better on plants that had been induced for 1 to 16 h than on un-induced plants suggesting they moved to and selected less toxic plant parts within a heterogeneously defended plant. Performance declined on all lines if plants had been induced for 24 h, or had been fed upon for a further 72 h. However, contrary to expectation, individual and total GS did not increase after these two treatments. This suggests that Arabidopsis plants induce additional (not GS) defenses after longer induction periods.

Studying the organization of genes encoding plant cell wall degrading enzymes in Chrysomela tremula provides insights into a leaf beetle genome.

The ability of herbivorous beetles from the superfamilies Chrysomeloidea and Curculionoidea to degrade plant cell wall polysaccharides has only recently begun to be appreciated. The presence of plant cell wall degrading enzymes (PCWDEs) in the beetle's digestive tract makes this degradation possible. Sequences encoding these beetle-derived PCWDEs were originally identified from transcriptomes and strikingly resemble those of saprophytic and phytopathogenic microorganisms, raising questions about their origin; e.g. are they insect- or microorganism-derived? To demonstrate unambiguously that the genes encoding PCWDEs found in beetle transcriptomes are indeed of insect origin, we generated a bacterial artificial chromosome library from the genome of the leaf beetle Chrysomela tremula, containing 18 432 clones with an average size of 143 kb. After hybridizing this library with probes derived from 12 C. tremula PCWDE-encoding genes and sequencing the positive clones, we demonstrated that the latter genes are encoded by the insect's genome and are surrounded by genes possessing orthologues in the genome of Tribolium castaneum as well as in three other beetle genomes. Our analyses showed that although the level of overall synteny between C. tremula and T. castaneum seems high, the degree of microsynteny between both species is relatively low, in contrast to the more closely related Colorado potato beetle.

One quantitative trait locus for intra- and interspecific variation in a sex pheromone.

Even though premating isolation is hypothesized to be a major driving force in speciation, its genetic basis is poorly known. In the noctuid moth Heliothis subflexa, one group of sex pheromone components, the acetates, emitted by the female, plays a crucial isolating role in preventing interspecific matings to males of the closely related Heliothis virescens, in which females do not produce acetates and males are repelled by them. We previously found intraspecific variation in acetates in H. subflexa: females in eastern North America contain significantly more acetates than females in Western Mexico. Here we describe the persistence of this intraspecific variation in laboratory-reared strains and the identification of one major quantitative trait locus (QTL), explaining 40% of the variance in acetate amounts. We homologized this intraspecific QTL to our previously identified interspecific QTL using restriction-associated DNA (RAD) tags. We found that a major intraspecific QTL overlaps with one of the two major interspecific QTL. To identify candidate genes underlying the acetate variation, we investigated a number of gene families with known or suspected acetyl- or acyltransferase activity. The most likely candidate genes did not map to our QTL, so that we currently hypothesize that a transcription factor underlies this QTL. Finding a single, large QTL that impacts variation in pheromone blends between and within species is, to our knowledge, the first such example for traits that have been demonstrated to affect premating isolation.

Coevolutionary fine-tuning: evidence for genetic tracking between a specialist wasp parasitoid and its aphid host in a dual metapopulation interaction.

In the interaction between two ecologically-associated species, the population structure of one species may affect the population structure of the other. Here, we examine the population structures of the aphid Metopeurum fuscoviride, a specialist on tansy Tanacetum vulgare, and its specialist primary hymenopterous parasitoid Lysiphlebus hirticornis, both of which are characterized by multivoltine life histories and a classic metapopulation structure. Samples of the aphid host and the parasitoid were collected from eight sites in and around Jena, Germany, where both insect species co-occur, and then were genotyped using suites of polymorphic microsatellite markers. The host aphid was greatly differentiated in terms of its spatial population genetic patterning, while the parasitoid was, in comparison, only moderately differentiated. There was a positive Mantel test correlation between pairwise shared allele distance (DAS) of the host and parasitoid, i.e. if host subpopulation samples were more similar between two particular sites, so were the parasitoid subpopulation samples. We argue that while the differences in the levels of genetic differentiation are due to the differences in the biology of the species, the correlations between host and parasitoid are indicative of dependence of the parasitoid population structure on that of its aphid host. The parasitoid is genetically tracking behind the aphid host, as can be expected in a classic metapopulation structure where host persistence depends on a delay between host and parasitoid colonization of the patch. The results may also have relevance to the Red Queen hypothesis, whereupon in the 'arms race' between parasitoid and its host, the latter 'attempts' to evolve away from the former.

Genetic differentiation across North America in the generalist moth Heliothis virescens and the specialist H. subflexa.

The two moth species Heliothis virescens (Hv) and H. subflexa (Hs) are closely related, but have vastly different feeding habits. Hv is a generalist and an important pest in many crops in the USA, while Hs is a specialist feeding only on plants in the genus Physalis. In this study, we conducted a comparative population genetic analysis to assess whether and how generalist and specialist life styles are reflected in differences in population structures. In Hv 98% of the total variation occurred within populations. The overall differentiation (F(ST) ) between regions was 0.006 and even lower between years (0.0039) and hosts (0.0028). Analyses of population structure suggest that all individuals form one genetically homogeneous population, except for at most 12 individuals (6%) that diverged from this cluster. Population homogeneity likely results from the high mobility of Hv and its generalist feeding behaviour. Hs exhibited substantially more population structure. Even though 96% of the total variation was attributable to within-population variability, F(ST) -values between Hs populations were 10 times higher than between Hv populations. Hs populations showed significant isolation by distance. Analyses of Hs population structure suggest at least two subpopulations and thus some degree of metapopulation structure. We speculate that the patchy distribution of Physalis- the exclusive food source of Hs - contributes to differences in population structure between these closely related species. The finding that the specialist shows more population differentiation than the generalist corroborates the notion that host specialization is not an evolutionary dead end but a dynamic trait.

Complex inheritance of larval adaptation in Plutella xylostella to a novel host plant.

Studying the genetics of host shifts and range expansions in phytophagous insects contributes to our understanding of the evolution of host plant adaptation. We investigated the recent host range expansion to pea, in the pea-adapted strain (P-strain) of the crucifer-specialist diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae). Larval survivorship on the novel host plant pea and a typical crucifer host (kale) was measured in reciprocal F(1), F(2) and backcrosses between the P-strain and a strain reared only on crucifers (C-strain). Reciprocal F(1) hybrids differed: offspring from P-strain mothers survived better on pea, indicating a maternal effect. However, no evidence for sex-linkage was found. Backcrosses to the P-strain produced higher survivorship on pea than C-strain backcrosses, suggesting recessive inheritance. In a linkage analysis with amplified fragment length polymorphism markers using P-strain backcrosses, two, four and five linkage groups contributing to survival on pea were identified in three different families respectively, indicating oligogenic inheritance. Thus, the newly evolved ability to survive on pea has a complex genetic basis, and the P-strain is still genetically heterogeneous and not yet fixed for all the alleles enabling it to survive on pea. Survivorship on kale was variable, but not related to survivorship on pea. This pattern may characterize the genetic inheritance of early host plant adaptation in oligophagous insect species.

Oviposition of diamondback moth in the presence and absence of a novel host plant.

The diamondback moth (DBM, Plutella xylostella L. (Lepidoptera: Plutellidae)) consumes a wide variety of brassicaceous host plants and is a common pest of crucifer crops worldwide. A highly unusual infestation of a sugar pea crop was recorded in Kenya in 1999, which persisted for two consecutive years. A strain (DBM-P) from this population was established in the laboratory and is the only one of several strains tested that can complete larval development on sugar peas. The oviposition acceptance and preference of the DBM-P strain was assessed in the presence of cabbage plants, sugar pea plants or both, in comparison to another strain (DBM-Cj) that was collected from cabbage and is unable to grow on pea plants. As expected, DBM-Cj females preferred to oviposit on cabbage plants. Surprisingly, DBM-P females also laid most eggs on cabbage and very few on peas. However, they laid significantly more eggs on the cabbage plant when pea plants were present. Our findings suggest that DBM-P manifested the initial stages of an evolutionary host range expansion, which is incomplete due to lack of oviposition fidelity on pea plants.

Widespread pyrethroid resistance in Australian diamondback moth, Plutella xylostella (L.), is related to multiple mutations in the para sodium channel gene.

Populations of Plutella xylostella, extending over 3800 km in southern Australia, show no genetic structure as assessed by microsatellite markers; yet outbreaks of pyrethroid resistance occur sporadically in cropping areas. Since mutations in the para voltage-gated sodium channel gene have been implicated in pyrethroid resistance, we looked for DNA sequence variation at this target among Australian moths. We found two resistance mutations previously reported for this species (L1014F and T929I), as well as a novel substitution (F1020S). Of the eight possible haplotypes formed by combinations of these three biallelic polymorphisms, only four were found in Australian populations: the wild-type allele (w), the kdr mutation allele (kdr) with only L1014F, the super-kdr-like combination of L1014F and T929I (skdrl), and the crashdown allele with only F1020S (cdr). Comparison of genotype frequencies among survivors of permethrin assays with those from untreated controls identified three resistant genotypes: skdrl homozygotes, cdr homozygotes and the corresponding heterozygote, cdr/skrdl - the heterozygote being at least as resistant as either homozygote. Spatial heterogeneity of allele frequencies was conspicuous, both across the continent and among local collections, consistent with reported spatial heterogeneity of pyrethroid resistance. Further, high resistance samples were sometimes associated with high frequency of cdr, sometimes high frequency of skdrl, or sometimes with a high combined cdr+skdrl frequency. The skdrl and cdr alleles explain a high proportion of the Australia-wide resistance variation. These data add to evidence that nerve insensitivity by mutations in the para-sodium channel gene is a common pyrethroid resistance mechanism in P. xylostella.

A haploid diamondback moth (Plutella xylostella L.) genome assembly resolves 31 chromosomes and identifies a diamide resistance mutation.

The diamondback moth, Plutella xylostella (L.), is a highly mobile brassica crop pest with worldwide distribution and can rapidly evolve resistance to insecticides, including group 28 diamides. Reference genomes assembled using Illumina sequencing technology have provided valuable resources to advance our knowledge regarding the biology, origin and movement of diamondback moth, and more recently with its sister species, Plutella australiana. Here we apply a trio binning approach to sequence and annotate a chromosome level reference genome of P. xylostella using PacBio Sequel and Dovetail Hi-C sequencing technology and identify a point mutation that causes resistance to commercial diamides. A P. xylostella population collected from brassica crops in the Lockyer Valley, Australia (LV-R), was reselected for chlorantraniliprole resistance then a single male was crossed to a P. australiana female and a hybrid pupa sequenced. A chromosome level 328 Mb P. xylostella genome was assembled with 98.1% assigned to 30 autosomes and the Z chromosome. The genome was highly complete with 98.4% of BUSCO Insecta genes identified and RNAseq informed protein prediction annotated 19,002 coding genes. The LV-R strain survived recommended field application doses of chlorantraniliprole, flubendiamide and cyclaniliprole. Some hybrids also survived these doses, indicating significant departure from recessivity, which has not been previously documented for diamides. Diamide chemicals modulate insect Ryanodine Receptors (RyR), disrupting calcium homeostasis, and we identified an amino acid substitution (I4790K) recently reported to cause diamide resistance in a strain from Japan. This chromosome level assembly provides a new resource for insect comparative genomics and highlights the emergence of diamide resistance in Australia. Resistance management plans need to account for the fact that resistance is not completely recessive.

Pyrosequencing the Manduca sexta larval midgut transcriptome: messages for digestion, detoxification and defence.

The tobacco hornworm Manduca sexta is an important model for insect physiology but genomic and transcriptomic data are currently lacking. Following a recent pyrosequencing study generating immune related expressed sequence tags (ESTs), here we use this new technology to define the M. sexta larval midgut transcriptome. We generated over 387,000 midgut ESTs, using a combination of Sanger and 454 sequencing, and classified predicted proteins into those involved in digestion, detoxification and immunity. In many cases the depth of 454 pyrosequencing coverage allowed us to define the entire cDNA sequence of a particular gene. Many new M. sexta genes are described including up to 36 new cytochrome P450s, some of which have been implicated in the metabolism of host plant-derived nicotine. New lepidopteran gene families such as the beta-fructofuranosidases, previously thought to be restricted to Bombyx mori, are also described. An unexpectedly high number of ESTs were involved in immunity, for example 39 contigs encoding serpins, and the increasingly appreciated role of the midgut in insect immunity is discussed. Similar studies of other tissues will allow for a tissue by tissue description of the M. sexta transcriptome and will form an essential complimentary step on the road to genome sequencing and annotation.

Phenotypic plasticity in sexual communication signal of a noctuid moth.

Variability within sex pheromone signalling systems is generally believed to be low because of strong stabilizing selection; yet the noctuid moth Heliothis subflexa (Hs) shows significant intraspecific variation. One possible explanation is that females may alter their sex pheromone blend depending on prevailing olfactory cues in the habitat, which we termed the 'experience hypothesis'. This could be adaptive if Hs females experiencing the pheromone of another species, Heliothis virescens (Hv), responded to reduce the frequency of heterospecific matings. We exposed Hs females to no pheromone, Hs pheromone or Hv pheromone in the first 3 days of their adult lives. Hs females in the latter treatment produced significantly more of the acetate Z11-16:OAc, which inhibits the attraction of Hv males. To our knowledge, this is the first study showing adaptive phenotypic plasticity in a moth sex pheromone and suggests that behavioural differentiation may precede genetic divergence in the sexual signals of moths.

Spatial population dynamics of a specialist aphid parasitoid, Lysiphlebus hirticornis Mackauer (Hymenoptera: Braconidae: Aphidiinae): evidence for philopatry and restricted dispersal.

Within insect communities, the population ecology of organisms representing higher trophic levels, for example, hymenopterous parasitoids, may be influenced by the structure of their insect hosts. Using microsatellite markers and ecological data, we investigated the population structure of the specialist braconid wasp parasitoid, Lysiphlebus hirticornis Mackauer attacking Metopeurum fuscoviride, a specialist aphid feeding on tansy, Tanacetum vulgare. Previous studies revealed that M. fuscoviride has a classic metapopulation structure with high subpopulation turnover. In this study, up to 100% of ramets within a host plant genet colonized by aphids were colonized by the parasitoid, yet plants with aphids but no parasitoids were also observed. Genetic differentiation measured by F(ST), actual differentiation (D) and relative differentiation (G(ST)) indicated highly structured parasitoid population demes, with restricted gene flow among and between parasitoid subpopulations at the various sites. Interestingly, both field data and population assignment analysis showed that the parasitoid is highly philopatric. Thus, despite the frequent local extinctions of the aphid host, the parasitoid continuously exploits its aphid host and contributes to the demise of local aphid subpopulations, rather than spreading its genes over many aphid populations. F(ST) values for the haplodiploid parasitoid were similar to those found in an independent study of the diploid aphid host, M. fuscoviride, hence supporting the view that an insect herbivore's population structure directly influences the ecology and genetics of the higher trophic level, in this case the wasp parasitoid.

Time-shifted reproductive behaviours among fall armyworm (Noctuidae: Spodoptera frugiperda) host strains: evidence for differing modes of inheritance.

The noctuid moth Spodoptera frugiperda consists of two strains associated with different larval host plants (most notably corn and rice). These strains exhibit differential temporal patterns of female calling and copulation during scotophase, with the corn strain more active earlier in the night. We investigated strain-specific constraints in reproductive timing, mating interactions between the two strains, and the mode of inheritance of timing of female calling, male calling, copulation and oviposition. We observed an allochronic shift of all reproductive behaviours by approximately 3 h and a parallel shift of nonreproductive locomotor activity, suggesting involvement of the circadian clock. The corn strain was more variable in the timing of calling and copulation than the rice strain. Rice strain females were more restricted in the timing of copulation than rice strain males, while such differences between the sexes were not apparent in the corn strain. There were significant interactions between the strains affecting onset times of copulation and male calling. The four investigated reproductive traits differed in their modes of inheritance: timing of female and male calling exhibited strong maternal effects, timing of copulation was controlled by a combination of maternal effects and corn strain dominant autosomal factors, and timing of oviposition was inherited in a corn strain dominant fashion. We conclude that the allochronic separation of reproduction between fall armyworm strains is asymmetric, less pronounced than previously thought, and under complex genetic control.

Identification of a gene associated with Bt resistance in Heliothis virescens.

Transgenic crops producing insecticidal toxins from Bacillus thuringiensis (Bt) are widely used for pest control. Bt-resistant insect strains have been studied, but the molecular basis of resistance has remained elusive. Here, we show that disruption of a cadherin-superfamily gene by retrotransposon-mediated insertion was linked to high levels of resistance to the Bt toxin Cry1Ac in the cotton pest Heliothis virescens. Monitoring the early phases of Bt resistance evolution in the field has been viewed as crucial but extremely difficult, especially when resistance is recessive. Our findings enable efficient DNA-based screening for resistant heterozygotes by directly detecting the recessive allele.

Identification of candidate genes for fenvalerate resistance in Helicoverpa armigera using cDNA-AFLP.

A cDNA-amplified fragment length polymorphisms approach was undertaken to screen for candidate genes associated with fenvalerate resistance in the AN02 strain of Helicoverpa armigera. Larvae and adults of this strain manifest approximately 50-fold resistance, which is suppressible by piperonyl butoxide and controlled by the semidominant gene RFen1 previously mapped to AFLP Linkage Group 13. Two cytochrome P450s (CYP337B1 and CYP4S1), one carboxylesterase-like protein and one glutathione transferase were found to be constitutively upregulated in resistant insects. Mapping of these potential detoxification genes showed that one of them, the novel P450 CYP337B1, was tightly linked to the resistance locus. This suggests that the RFen1(R) allele has a cis-acting effect on CYP337B1 expression, and possible trans-acting effects on expression of other genes.

Data set for diet specific differential gene expression analysis in three Spodoptera moths.

Examination of closely related species pairs is suggested for evolutionary comparisons of different degrees of polyphagy, which we did here with three taxa of lepidopteran herbivores, Spodoptera spp (S. littoralis, S. frugiperda maize (C) and rice (R) strains) for a RNAseq analysis of the midguts from the 3rd instar insect larvae for differential metabolic responses after feeding on pinto bean based artificial diet vs maize leaves. Paired-end (2×100 bp) Illumina HiSeq2500 sequencing resulted in a total of 24, 23, 24, and 21 million reads for the SF-C-Maize, SF-C-Pinto, SF-R-Maize, SF-R Pinto, and a total of 35 and 36 million reads for the SL-Maize and SL-Pinto samples, respectively. After quality control measures, a total of 62.2 million reads from SL and 71.7 million reads from SF were used for transcriptome assembly (TA). The resulting final de novo reference TA (backbone) for the SF taxa contained 37,985 contigs with a N50 contig size of 1030 bp and a maximum contig length of 17,093 bp, while for SL, 28,329 contigs were generated with a N50 contig size of 1980 bp and a maximum contig length of 18,267 bp. The data presented herein contains supporting information related to our research article Roy et al. (2016) http://dx.doi.org/10.1016/j.ibmb.2016.02.006[1].

Diet dependent metabolic responses in three generalist insect herbivores Spodoptera spp.

Adaption to dietary changes is critical in the evolution of host plant ranges in polyphagous insects. We compared three taxa of lepidopteran herbivores from the predominantly generalist genus Spodoptera showing different degrees of polyphagy: Spodoptera littoralis, with a broad host range including both mono- and dicotyledonous plants, and two Spodoptera frugiperda strains [Corn (i.e. maize) (C) and Rice (R)] adapted primarily to different grass species. When feeding on maize we show a lower performance in the broad generalist taxon compared to the grass adapted taxa. Among these taxa, the maize adapted S. frugiperda C-strain generally performed better than the R-strain on maize leaves. On artificial pinto diet, all taxa performed well. Our RNA-Seq analysis of midgut transcriptomes from 3rd instar larvae feeding on maize showed broader transcriptional readjustments in the generalist S. littoralis compared to grass adapted S. frugiperda strains. Substantial alteration in the expression levels of midgut physiological function related transcripts, such as digestive and detoxifying enzymes, transporters, immunity, and peritrophic membrane associated transcripts, existed in all taxa. We found high background expression of UDP-glucosyl transferases, which are known to neutralize maize leaf toxins, in the maize adapted S. frugiperda C-strain, contributing to its fitness on maize compared to the R-strain. Our findings provide evidence for divergent diet specific response of digestive physiology within these Spodoptera taxa. Unexpectedly, the C- and R-strains of S. frugiperda fed on the same diet showed large differences in expression patterns between these two closely related taxa.

Interactions between pesticide genes: model and experiment.

In response to years of intense selection pressure by organophosphate insecticides, several different insecticide resistance mechanisms have evolved in natural populations of the mosquito Culex pipiens. We examined interactions between two of the most important mechanisms using a four-compartment model of insecticide pharmacokinetics. The joint effect of different mechanisms of resistance can be expressed in terms of epistasis at the physiological level in this model. The type of epistasis predicted by the model depends on the particular physiological mechanisms of resistance involved. Resistance due to a reduced penetration of the insecticide combines multiplicatively with other resistance factors, but resistance due to detoxicative processes and to insensitivity of the target site combines additively. How the pattern of epistasis at the physiological level is translated into fitness epistasis in natural populations of this mosquito depends on the intensity and pattern of insecticide selection in the field.

A direct assessment of the role of genetic drift in determining allele frequency variation in populations of Euphydryas editha.

Estimates of allele frequencies at six polymorphic loci were collected over eight generations in two populations of Euphydryas editha. We have estimated, in addition, the effective population size for each generation for both populations with results from mark-recapture and other field data. The variation in allele frequencies generated by random genetic drift was then studied using computer simulations and our direct estimates of effective population size. Substantial differences between observed values and computer-generated expected values assuming drift alone were found for three loci (Got, Hk, Pgi) in one population. These observations are consistent with natural selection in a variable environment.