Spodoptera frugiperda: Ecology, Evolution, and Management Options of an Invasive Species.

The fall armyworm (FAW), Spodoptera frugiperda (Lepidoptera, Noctuidae), is a well-known agricultural pest in its native range, North and South America, and has become a major invasive pest around the globe in the past decade. In this review, we provide an overview to update what is known about S. frugiperda in its native geographic ranges. This is followed by discussion of studies from the invaded areas to gain insights into S. frugiperda's ecology, specifically its reproductive biology, host plant use, status of insecticide resistance alleles, and biocontrol methods in native and invasive regions. We show that reference to host strains is uninformative in the invasive populations because multidirectional introduction events likely underpinned its recent rapid spread. Given that recent genomic analyses show that FAW is much more diverse than was previously assumed, and natural selection forces likely differ geographically, region-specific approaches will be needed to control this global pest.

Virulence of Beauveria sp. and Metarhizium sp. fungi towards fall armyworm (Spodoptera frugiperda).

The development of effective pest management strategies for Spodoptera frugiperda is a high priority for crop protection across its invasive ranges. Here, we examined six Beauveria and five Metarhizium fungal isolates against this pest. Two Beauveria isolates (B-0571, B-1311) induced high mortality toward 3rd and 6th instar caterpillars and adults. For B-0571 mortality was 82.81 ± 5.75%, 61.46 ± 6.83%, and 93.75 ± 3.61%, and 73.72 ± 2.51%, 71.88 ± 5.41%, and 97.92 ± 2.08% for B-1311, with deaths in caterpillars largely occurring under 24 h (3rd instar control 0.74 ± 0.33%, B-0571 73.96 ± 7.85% and B-1311 62.08 ± 3.67%; 6th instar control 0%, B-0571 66.67% ± 11.02% and B-1311 62.5% ± 9.55%). Infection from both Beauveria isolates fully prevented reproduction in surviving S. frugiperda females. In contrast, all five Metarhizium isolates tested and the remaining four Beauveria isolates exhibited lower virulence. The discovery of two highly virulent Beauveria fungal isolates to S. frugiperda opens avenues to develop novel biological control tools against this highly invasive pest.

Adaptive Introgression across Semipermeable Species Boundaries between Local Helicoverpa zea and Invasive Helicoverpa armigera Moths.

Hybridization between invasive and native species has raised global concern, given the dramatic increase in species range shifts and pest outbreaks due to anthropogenic dispersal. Nevertheless, secondary contact between sister lineages of local and invasive species provides a natural laboratory to understand the factors that determine introgression and the maintenance or loss of species barriers. Here, we characterize the early evolutionary outcomes following secondary contact between invasive Helicoverpa armigera and native H. zea in Brazil. We carried out whole-genome resequencing of Helicoverpa moths from Brazil in two temporal samples: during the outbreak of H. armigera in 2013 and 2017. There is evidence for a burst of hybridization and widespread introgression from local H. zea into invasive H. armigera coinciding with H. armigera expansion in 2013. However, in H. armigera, the admixture proportion and the length of introgressed blocks were significantly reduced between 2013 and 2017, suggesting selection against admixture. In contrast to the genome-wide pattern, there was striking evidence for adaptive introgression of a single region from the invasive H. armigera into local H. zea, including an insecticide resistance allele that increased in frequency over time. In summary, despite extensive gene flow after secondary contact, the species boundaries are largely maintained except for the single introgressed region containing the insecticide-resistant locus. We document the worst-case scenario for an invasive species, in which there are now two pest species instead of one, and the native species has acquired resistance to pyrethroid insecticides through introgression.

Origin, behaviour, and genetics of reproductive workers in an invasive ant.

BACKGROUND: Worker reproduction has an important influence on the social cohesion and efficiency of social insect colonies, but its role in the success of invasive ants has been neglected. We used observations of 233 captive colonies, laboratory experiments, and genetic analyses to investigate the conditions for worker reproduction in the invasive Anoplolepis gracilipes (yellow crazy ant) and its potential cost on interspecific defence. We determined the prevalence of worker production of males and whether it is triggered by queen absence; whether physogastric workers with enlarged abdomens are more likely to be reproductive, how normal workers and physogastric workers compare in their contributions to foraging and defence; and whether worker-produced males and males that could have been queen- or worker-produced differ in their size and heterozygosity. RESULTS: Sixty-six of our 233 captive colonies produced males, and in 25 of these, some males could only have been produced by workers. Colonies with more workers were more likely to produce males, especially for queenless colonies. The average number of days between the first appearance of eggs and adult males in our colonies was 54.1 ± 10.2 (mean ± SD, n = 20). In our laboratory experiment, queen removal triggered an increase in the proportion of physogastric workers. Physogastric workers were more likely to have yolky oocytes (37-54.9%) than normal workers (2-25.6%), which is an indicator of fertile or trophic egg production. Physogastric workers were less aggressive during interspecific aggression tests and foraged less than normal workers. The head width and wing length of worker-produced males were on average 4.0 and 4.3% greater respectively than those of males of undetermined source. Our microsatellite DNA analyses indicate that 5.5% of worker-produced males and 14.3% of males of undetermined source were heterozygous, which suggests the presence of diploid males and/or genetic mosaics in A. gracilipes. CONCLUSIONS: Our experimental work provides crucial information on worker reproduction in A. gracilipes and its potential cost to colony defence. The ability of A. gracilipes workers to produce males in the absence of queens may also contribute to its success as an invasive species if intranidal mating can take place between virgin queens and worker-produced males.

Hybridization and gene flow in the mega-pest lineage of moth, Helicoverpa.

Within the mega-pest lineage of heliothine moths are a number of polyphagous, highly mobile species for which the exchange of adaptive traits through hybridization would affect their properties as pests. The recent invasion of South America by one of the most significant agricultural pests, Helicoverpa armigera, raises concerns for the formation of novel combinations of adaptive genes following hybridization with the closely related Helicoverpa zea To investigate the propensity for hybridization within the genus Helicoverpa, we carried out whole-genome resequencing of samples from six species, focusing in particular upon H. armigera population structure and its relationship with H. zea We show that both H. armigera subspecies have greater genetic diversity and effective population sizes than do the other species. We find no signals for gene flow among the six species, other than between H. armigera and H. zea, with nine Brazilian individuals proving to be hybrids of those two species. Eight had largely H. armigera genomes with some introgressed DNA from H. zea scattered throughout. The ninth resembled an F1 hybrid but with stretches of homozygosity for each parental species that reflect previous hybridization. Regions homozygous for H. armigera-derived DNA in this individual included one containing a gustatory receptor and esterase genes previously associated with host range, while another encoded a cytochrome P450 that confers insecticide resistance. Our data point toward the emergence of novel hybrid ecotypes and highlight the importance of monitoring H. armigera genotypes as they spread through the Americas.

Movement Ecology of Pest Helicoverpa: Implications for Ongoing Spread.

The recent introduction and spread of Helicoverpa armigera throughout South America highlight the invasiveness and adaptability of moths in the Helicoverpa genus. Long-range movement in three key members, H. armigera, H. zea, and H. punctigera, occurs by migration and international trade. These movements facilitate high population admixture and genetic diversity, with important economic, biosecurity, and control implications in today's agricultural landscape. This is particularly true for the spread of resistance alleles to transgenic crops expressing Bacillus thuringiensis (Bt) toxins that are planted over vast areas to suppress Helicoverpa spp. The ability to track long-distance movement through radar technology, population genetic markers, and/or long-distance dispersal modeling has advanced in recent years, yet we still know relatively little about the population trajectories or migratory routes in Helicoverpa spp. Here, we consider how experimental and theoretical approaches can be integrated to fill key knowledge gaps and assist management practices.

Insect Resistance to Bacillus thuringiensis Toxin Cry2Ab Is Conferred by Mutations in an ABC Transporter Subfamily A Protein.

The use of conventional chemical insecticides and bacterial toxins to control lepidopteran pests of global agriculture has imposed significant selection pressure leading to the rapid evolution of insecticide resistance. Transgenic crops (e.g., cotton) expressing the Bt Cry toxins are now used world wide to control these pests, including the highly polyphagous and invasive cotton bollworm Helicoverpa armigera. Since 2004, the Cry2Ab toxin has become widely used for controlling H. armigera, often used in combination with Cry1Ac to delay resistance evolution. Isolation of H. armigera and H. punctigera individuals heterozygous for Cry2Ab resistance in 2002 and 2004, respectively, allowed aspects of Cry2Ab resistance (level, fitness costs, genetic dominance, complementation tests) to be characterised in both species. However, the gene identity and genetic changes conferring this resistance were unknown, as was the detailed Cry2Ab mode of action. No cross-resistance to Cry1Ac was observed in mutant lines. Biphasic linkage analysis of a Cry2Ab-resistant H. armigera family followed by exon-primed intron-crossing (EPIC) marker mapping and candidate gene sequencing identified three independent resistance-associated INDEL mutations in an ATP-Binding Cassette (ABC) transporter gene we named HaABCA2. A deletion mutation was also identified in the H. punctigera homolog from the resistant line. All mutations truncate the ABCA2 protein. Isolation of further Cry2Ab resistance alleles in the same gene from field H. armigera populations indicates unequal resistance allele frequencies and the potential for Bt resistance evolution. Identification of the gene involved in resistance as an ABC transporter of the A subfamily adds to the body of evidence on the crucial role this gene family plays in the mode of action of the Bt Cry toxins. The structural differences between the ABCA2, and that of the C subfamily required for Cry1Ac toxicity, indicate differences in the detailed mode-of-action of the two Bt Cry toxins.

A genomic approach to identify and monitor a novel pyrethroid resistance mutation in the redlegged earth mite, Halotydeus destructor.

Resistance mechanisms are typically uncovered by identifying sequence variation in known candidate genes, however this strategy can be problematic for species with no reference data in known relatives. Here we take a genomic approach to identify resistance to pyrethroids in the redlegged earth mite, Halotydeus destructor, a member of the Penthalidae family of mites that are virtually uncharacterized genetically. Based on shallow genome sequencing followed by a genome assembly, we first identified contigs of the H. destructor parasodium channel gene. By linking variation in this gene to known resistant phenotypes, we located a single nucleotide polymorphism in resistant mites. This polymorphism results in a leucine (L) to phenylalanine (F) amino acid substitution in the II6 region (predicted) of the gene (L1024F). This novel mutation has not previously been linked to pyrethroid resistance, although other polymorphisms have been identified in the two-spotted spider mite, Tetranychus urticae at the same locus (L1024V). The sequencing approach was successful in generating a candidate polymorphism that was then validated using laboratory bioassays and field surveys. A high throughput Illumina-based sequencing diagnostic was developed to rapidly assess resistance allele frequencies in pools of mites sourced from hundreds of populations across Australia. Resistance was confirmed to be widespread in the southern wheatbelt region of Western Australia. Two different resistance mutations were identified in field populations, both resulting in the same amino acid substitution. The frequency and distribution of resistance amplicon haplotypes suggests at least two, and probably more independent origins of resistance.

Complex multiple introductions drive fall armyworm invasions into Asia and Australia.

The fall armyworm (FAW) Spodoptera frugiperda is thought to have undergone a rapid 'west-to-east' spread since 2016 when it was first identified in western Africa. Between 2018 and 2020, it was recorded from South Asia (SA), Southeast Asia (SEA), East Asia (EA), and Pacific/Australia (PA). Population genomic analyses enabled the understanding of pathways, population sources, and gene flow in this notorious agricultural pest species. Using neutral single nucleotide polymorphic (SNP) DNA markers, we detected genome introgression that suggested most populations in this study were overwhelmingly C- and R-strain hybrids (n = 252/262). SNP and mitochondrial DNA markers identified multiple introductions that were most parsimoniously explained by anthropogenic-assisted spread, i.e., associated with international trade of live/fresh plants and plant products, and involved 'bridgehead populations' in countries to enable successful pest establishment in neighbouring countries. Distinct population genomic signatures between Myanmar and China do not support the 'African origin spread' nor the 'Myanmar source population to China' hypotheses. Significant genetic differentiation between populations from different Australian states supported multiple pathways involving distinct SEA populations. Our study identified Asia as a biosecurity hotspot and a FAW genetic melting pot, and demonstrated the use of genome analysis to disentangle preventable human-assisted pest introductions from unpreventable natural pest spread.

Screening for insecticide resistance in Australian field populations of Bemisia tabaci (Hemiptera: Aleyrodidae) using bioassays and DNA sequencing.

BACKGROUND: Species within the Bemisia tabaci cryptic species complex can cause significant crop damage. We used high-throughput amplicon sequencing to identify the species composition and resistance allele genotypes in field populations from cotton fields in Australia. For selected populations, the resistance phenotype was determined in bioassays and compared with sequencing data. RESULTS: A metabarcoding approach was used to analyse the species composition in 144 field populations collected between 2013 and 2021. Two mixed AUS I and MEAM1 populations were detected, whereas the remaining 142 populations consisted of MEAM1 only. High-throughput sequencing of organophosphate and pyrethroid resistance gene amplicons showed that the organophosphate resistance allele F331W was fixed (> 99%) in all MEAM1 populations, whereas the pyrethroid resistance allele L925I in the voltage-gated sodium channel gene was detected at varying frequencies [1.0%-7.0% (43 populations); 27.7% and 42.1% (two populations); 95%-97.5% (three populations)]. Neither organophosphate nor pyrethroid resistance alleles were detected in the AUS I populations. Pyrethroid bioassays of 85 MEAM1 field-derived populations detected no resistance in 51 populations, whereas 32 populations showed low frequency resistance, and 2 populations were highly resistant. CONCLUSIONS: We demonstrate that high-throughput sequencing and bioassays are complementary approaches. The detection of target site mutations and the phenotypic provides a comprehensive analysis of the low-level resistance to pyrethroids that is present in Australian cotton farms. By contrast, a limited survey of whitefly populations from horticulture found evidence of high-level resistance against pyrethroids. Furthermore, we found that the F331W allele (linked to organophosphate resistance) is ubiquitous in Australian MEAM1. © 2022 Commonwealth of Australia. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

A diagnostic LAMP assay for rapid identification of an invasive plant pest, fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae).

Fall armyworm (FAW), Spodoptera frugiperda (Lepidoptera: Noctuidae), is a highly polyphagous invasive plant pest that has expanded its global geographic distribution, including recently into much of Australia. Rapid diagnostic tests are required for identification of FAW to assist subsequent management and control. We developed a new loop-mediated isothermal amplification (LAMP) assay based on the mitochondrial cytochrome c oxidase subunit I (COI) gene for accurate and timely diagnosis of FAW in the field. The specificity of the new assay was tested against a broad panel of twenty non-target noctuids, including eight other Spodoptera species. Only S. frugiperda samples produced amplification within 20 min, with an anneal derivative temperature of 78.3 ± 0.3 °C. A gBlock dsDNA fragment was developed and trialled as a synthetic positive control, with a different anneal derivative of 81 °C. The new FAW LAMP assay was able to detect FAW DNA down to 2.4 pg, similar to an existing laboratory-based real-time PCR assay. We also trialled the new FAW assay with a colorimetric master mix and found it could successfully amplify positive FAW samples in half the time compared to an existing FAW colorimetric LAMP assay. Given the high sensitivity and rapid amplification time, we recommend the use of this newly developed FAW LAMP assay in a portable real-time fluorometer for in-field diagnosis of FAW.

A high-throughput amplicon sequencing approach for population-wide species diversity and composition survey.

Management of agricultural pests requires an understanding of pest species diversity, their interactions with beneficial insects and spatial-temporal patterns of pest abundance. Invasive and agriculturally important insect pests can build up very high populations, especially in cropping landscapes. Traditionally, sampling effort for species identification involves small sample sizes and is labour intensive. Here, we describe a multiprimer high throughput sequencing (HTS) metabarcoding method and associated analytical workflow for a rapid, intensive, high-volume survey of pest species compositions. We demonstrate our method using the taxonomically challenging Bemisia pest cryptic species complex as examples. The whiteflies Bemisia including the"tabaci" species are agriculturally important capable of vectoring diverse plant viruses that cause diseases and crop losses. Our multiprimer metabarcoding HTS amplicon approach simultaneously process high volumes of whitefly individuals, with efficiency to detect rare (i.e., 1%) test-species, while our improved whitefly primers for metabarcoding also detected beneficial hymenopteran parasitoid species from whitefly nymphs. Field-testing our redesigned Bemisia metabarcoding primer sets across the Tanzania, Uganda and Malawi cassava cultivation landscapes, we identified the sub-Saharan Africa 1 Bemisia putative species as the dominant pest species, with other cryptic Bemisia species being detected at various abundances. We also provide evidence that Bemisia species compositions can be affected by host crops and sampling techniques that target either nymphs or adults. Our multiprimer HTS metabarcoding method incorporated two overlapping amplicons of 472 bp and 518 bp that spanned the entire 657 bp 3' barcoding region for Bemisia, and is particularly suitable to molecular diagnostic surveys of this highly cryptic insect pest species complex that also typically exhibited high population densities in heavy crop infestation episodes. Our approach can be adopted to understand species biodiversity across landscapes, with broad implications for improving transboundary biosecurity preparedness, thus contributing to molecular ecological knowledge and the development of control strategies for high-density, cryptic, pest-species complexes.

Population genomics provides insights into lineage divergence and local adaptation within the cotton bollworm.

The cotton bollworm Helicoverpa armigera is a cosmopolitan pest and its diverse habitats plausibly contribute to the formation of diverse lineages. Despite the significant threat it poses to economic crops worldwide, its evolutionary history and genetic basis of local adaptation are poorly understood. In this study, we de novo assembled a high-quality chromosome-level reference genome of H. a. armigera (contig N50 = 7.34 Mb), with 99.13% of the HaSCD2 assembly assigned to 31 chromosomes (Z-chromosome + 30 autosomes). We constructed an ultradense variation map across 14 cotton bollworm populations and identified a novel lineage in northwestern China. Historical inference showed that effective population size changes coincided with global temperature fluctuation. We identified nine differentiated genes in the three H. armigera lineages (H. a. armigera, H. a. conferta and the new northwestern Chinese lineage), of which per and clk genes are involved in circadian rhythm. Selective sweep analyses identified a series of Gene Ontology categories related to climate adaptation, feeding behaviour and insecticide tolerance. Our findings reveal fundamental knowledge of the local adaptation of different cotton bollworm lineages and will guide the formulation of cotton bollworm management measures at different scales.

Population Genetic Structure of the Bean Leaf Beetle Ootheca mutabilis (Coleoptera: Chrysomelidae) in Uganda.

Bean leaf beetle (BLB) (Ootheca mutabilis) has emerged as an important bean pest in Uganda, leading to devastating crop losses. There is limited information on the population genetic structure of BLB despite its importance. In this study, novel microsatellite DNA markers and the partial mitochondrial cytochrome oxidase subunit I (mtCOI) gene sequences were used to analyze the spatial population genetic structure, genetic differentiation and haplotype diversity of 86 O. mutabilis samples from 16 (districts) populations. We identified 19,356 simple sequence repeats (SSRs) (mono, di-, tri-, tetra-, penta-, and hexa-nucleotides) of which 81 di, tri and tetra-nucleotides were selected for primer synthesis. Five highly polymorphic SSR markers (4-21 alleles, heterozygosity 0.59-0.84, polymorphic information content (PIC) 50.13-83.14%) were used for this study. Analyses of the 16 O. mutabilis populations with these five novel SSRs found nearly all the genetic variation occurring within populations and there was no evidence of genetic differentiation detected for both types of markers. Also, there was no evidence of isolation by distance between geographical and genetic distances for SSR data and mtCOI data except in one agro-ecological zone for mtCOI data. Bayesian clustering identified a signature of admixture that suggests genetic contributions from two hypothetical ancestral genetic lineages for both types of markers, and the minimum-spanning haplotype network showed low differentiation in minor haplotypes from the most common haplotype with the most common haplotype occurring in all the 16 districts. A lack of genetic differentiation indicates unrestricted migrations between populations. This information will contribute to the design of BLB control strategies.

Global population genomic signature of Spodoptera frugiperda (fall armyworm) supports complex introduction events across the Old World.

Native to the Americas, the invasive Spodoptera frugiperda (fall armyworm; FAW) was reported in West Africa in 2016, followed by its chronological detection across the Old World and the hypothesis of an eastward Asia expansion. We explored population genomic signatures of American and Old World FAW and identified 12 maternal mitochondrial DNA genome lineages across the invasive range. 870 high-quality nuclear single nucleotide polymorphic DNA markers identified five distinct New World population clusters, broadly reflecting FAW native geographical ranges and the absence of host-plant preferences. We identified unique admixed Old World populations, and admixed and non-admixed Asian FAW individuals, all of which suggested multiple introductions underpinning the pest's global spread. Directional gene flow from the East into eastern Africa was also detected, in contrast to the west-to-east spread hypothesis. Our study demonstrated the potential of population genomic approaches via international partnership to address global emerging pest threats and biosecurity challenges.

The evolutionary process of invasion in the fall armyworm (Spodoptera frugiperda).

The fall armyworm (FAW; Spodoptera frugiperda) is one of the major agricultural pest insects. FAW is native to the Americas, and its invasion was first reported in West Africa in 2016. Then it quickly spread through Africa, Asia, and Oceania, becoming one of the main threats to corn production. We analyzed whole genome sequences of 177 FAW individuals from 12 locations on four continents to infer evolutionary processes of invasion. Principal component analysis from the TPI gene and whole genome sequences shows that invasive FAW populations originated from the corn strain. Ancestry coefficient and phylogenetic analyses from the nuclear genome indicate that invasive populations are derived from a single ancestry, distinct from native populations, while the mitochondrial phylogenetic tree supports the hypothesis of multiple introductions. Adaptive evolution specific to invasive populations was observed in detoxification, chemosensory, and digestion genes. We concluded that extant invasive FAW populations originated from the corn strain with potential contributions of adaptive evolution.

Within-Season Changes in Land-Use Impact Pest Abundance in Smallholder African Cassava Production Systems.

Cassava (Manihot esculenta Crantz), an important commercial and food security crop in East and Central Africa, continues to be adversely affected by the whitefly Bemisia tabaci. In Uganda, changes in smallholder farming landscapes due to crop rotations can impact pest populations but how these changes affect pest outbreak risk is unknown. We investigated how seasonal changes in land-use have affected B. tabaci population dynamics and its parasitoids. We used a large-scale field experiment to standardize the focal field in terms of cassava age and cultivar, then measured how Bemisia populations responded to surrounding land-use change. Bemisia tabaci Sub-Saharan Africa 1 (SSA1) was identified using molecular diagnostics as the most prevalent species and the same species was also found on surrounding soybean, groundnut, and sesame crops. We found that an increase in the area of cassava in the 3-7-month age range in the landscape resulted in an increase in the abundance of the B. tabaci SSA1 on cassava. There was a negative relationship between the extent of non-crop vegetation in the landscape and parasitism of nymphs suggesting that these parasitoids do not rely on resources in the non-crop patches. The highest abundance of B. tabaci SSA1 nymphs in cassava fields occurred at times when landscapes had large areas of weeds, low to moderate areas of maize, and low areas of banana. Our results can guide the development of land-use strategies that smallholder farmers can employ to manage these pests.

Geographic Monitoring of Insecticide Resistance Mutations in Native and Invasive Populations of the Fall Armyworm.

Field evolved resistance to insecticides is one of the main challenges in pest control. The fall armyworm (FAW) is a lepidopteran pest species causing severe crop losses, especially corn. While native to the Americas, the presence of FAW was confirmed in West Africa in 2016. Since then, the FAW has been detected in over 70 countries covering sub-Saharan Africa, the Middle East, North Africa, South Asia, Southeast Asia, and Oceania. In this study, we tested whether this invasion was accompanied by the spread of resistance mutations from native to invasive areas. We observed that mutations causing Bt resistance at ABCC2 genes were observed only in native populations where the mutations were initially reported. Invasive populations were found to have higher gene numbers of cytochrome P450 genes than native populations and a higher proportion of multiple resistance mutations at acetylcholinesterase genes, supporting strong selective pressure for resistance against synthetic insecticides. This result explains the susceptibility to Bt insecticides and resistance to various synthetic insecticides in Chinese populations. These results highlight the necessity of regular and standardized monitoring of insecticide resistance in invasive populations using both genomic approaches and bioassay experiments.

Whole-genome sequencing to detect mutations associated with resistance to insecticides and Bt proteins in Spodoptera frugiperda.

The fall armyworm (FAW), Spodoptera frugiperda, is a major pest native to the Americas that has recently invaded the Old World. Point mutations in the target-site proteins acetylcholinesterase-1 (ace-1), voltage-gated sodium channel (VGSC) and ryanodine receptor (RyR) have been identified in S. frugiperda as major resistance mechanisms to organophosphate, pyrethroid and diamide insecticides respectively. Mutations in the adenosine triphosphate-binding cassette transporter C2 gene (ABCC2) have also been identified to confer resistance to Cry1F protein. In this study, we applied a whole-genome sequencing (WGS) approach to identify point mutations in the target-site genes in 150 FAW individuals collected from China, Malawi, Uganda and Brazil. This approach revealed three amino acid substitutions (A201S, G227A and F290V) of S. frugiperda ace-1, which are known to be associated with organophosphate resistance. The Brazilian population had all three ace-1 point mutations and the 227A allele (mean frequency = 0.54) was the most common. Populations from China, Malawi and Uganda harbored two of the three ace-1 point mutations (A201S and F290V) with the 290V allele (0.47-0.58) as the dominant allele. Point mutations in VGSC (T929I, L932F and L1014F) and RyR (I4790M and G4946E) were not detected in any of the 150 individuals. A novel 12-bp insertion mutation in exon 15 of the ABCC2 gene was identified in some of the Brazilian individuals but absent in the invasive populations. Our results not only demonstrate robustness of the WGS-based genomic approach for detection of resistance mutations, but also provide insights for improvement of resistance management tactics in S. frugiperda.

Generation of microsatellite repeat families by RTE retrotransposons in lepidopteran genomes.

BACKGROUND: Developing lepidopteran microsatellite DNA markers can be problematical, as markers often exhibit multiple banding patterns and high frequencies of non-amplifying "null" alleles. Previous studies identified sequences flanking simple sequence repeat (SSR) units that are shared among many lepidopteran species and can be grouped into microsatellite-associated DNA families. These families are thought to be associated with unequal crossing-over during DNA recombination or with transposable elements (TEs). RESULTS: We identified full-length lepidopteran non-LTR retrotransposable elements of the RTE clade in Heliconius melpomene and Bombyx mori. These retroelements possess a single open reading frame encoding the Exonuclease/Endonuclease/Phosphatase and the Reverse Transcriptase/nLTR domains, a 5' UTR (untranslated region), and an extremely short 3' UTR that regularly consists of SSR units. Phylogenetic analysis supported previous suggestions of horizontal transfer among unrelated groups of organisms, but the diversity of lepidopteran RTE elements appears due to ancient divergence of ancestral elements rather than introgression by horizontal transfer. Similarity searches of lepidopteran genomic sequences in GenBank identified partial RTE elements, usually consisting of the 3' terminal region, in 29 species. Furthermore, we identified the C-terminal end of the Reverse Transcriptase/nLTR domain and the associated 3' UTR in over 190 microsatellite markers from 22 lepidopteran species, accounting for 10% of the lepidopteran microsatellites in GenBank. Occasional retrotransposition of autonomous elements, frequent retrotransposition of 3' partial elements, and DNA replication slippage during retrotransposition offers a mechanistic explanation for the association of SSRs with RTE elements in lepidopteran genomes. CONCLUSIONS: Non-LTR retrotransposable elements of the RTE clade therefore join a diverse group of TEs as progenitors of SSR units in various organisms. When microsatellites are isolated using standard SSR enrichment protocols and primers designed at complementary repeated regions, amplification from multiple genomic sites can cause scoring difficulties that compromise their utility as markers. Screening against RTE elements in the isolation procedure provides one strategy for minimizing this problem.