Limited information exists on molecular genetic variation and distribution of the corn and rice strains of the fall armyworm, Spodoptera frugiperda (J.E. Smith). This study was conducted to investigate the genetic structure of S. frugiperda across a part of its range in the United States. A 608-base-pair portion of the mitochondrial cytochrome oxidase I and II genes was sequenced from 71 individuals resulting in three corn and four rice strain haplotypes. Genetic divergence between the two strains ranged from 0.66 to 0.99%. A 562-base-pair region of the nuclear ITS-1 gene was also amplified and sequenced from 17 individuals representing both corn and rice strains. No variation was detected in any of the samples for the ITS-1 region. Analysis of molecular variance was conducted on the resulting mtDNA haplotypes from the Arkansas and Florida populations and as a hierarchical analysis between populations in the two states. Results indicate a significant overall ΦST for all populations with the hierarchical analysis revealing that this significant ΦST is due to structuring of the populations between states. The observed genetic structure is possibly due to the distribution of fall armyworm strains.
The fall armyworm, Spodoptera frugiperda (J.E. Smith), is a major pest on corn, sorghum, and bermudagrass in the southeastern United States (Knipling 1980; Pashley 1986; Sparks 1979). The preferred host plants of the fall armyworm came under new scrutiny in 1986 when Pashley proposed that the fall armyworm consists of two morphologically undistinguishable strains, a corn strain that prefers corn, cotton, and sorghum, and a rice strain that prefers rice and bermudagrass (Pashley 1986, 1988a). The range of S. frugiperda is known to cover most of the Western hemisphere, and the range of each strain, however, has been examined from Louisiana down through Central America and in the Caribbean to Brazil (Pashley et al. 1985; Pashley 1986, 1988b).
Despite the possible benefits that population genetic analysis of the fall armyworm may provide towards understanding dispersal, monitoring the spread of insecticide resistance, and the implementation of area-wide control programs, relatively little research in this area has been conducted. A survey of 22 allozyme loci by Pashley et al. (1985) indicated significant heterogeneity between populations at five of 11 polymorphic loci, due in large part to the distinctness of a single Puerto Rican population collected from rice. The phylogenetic relationships between the two strains were further examined with three of these polymorphic allozymes (Hbdh, PepF, and Est3; Pashley 1988b). The majority of the genetic studies have focused on differentiating the rice and corn strains with polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP), strain specific PCR, RFLP, amplified fragment length polymorphism (AFLP) and allozyme markers (Meagher & Gallo-Meagher 2003; Levy et al. 2002; Nagoshi & Meagher 2003; McMichael & Prowell 1999; Pashley et al. 1985; Lu et al. 1992; Adamczyk 1993; Lu & Adang 1996; Pashley 1989). A genetic variation study by Lu et al. (1992) involving RFLP of a random genomic library from six populations (five of which were lab colonies) from Louisiana, Mississippi, and Georgia revealed high levels of genetic variation within and among populations. However, no population genetic analysis was conducted in that study, which focused on finding diagnostic markers for the corn and rice strains.
Mitochondrial-DNA (mtDNA) analysis is generally assumed to be more powerful than allozyme analysis for revealing population structure, and has been used for numerous population genetic studies (Avise 1994). The cytochrome oxidase I (COI) and cytochrome oxidase II (COII) regions of the mtDNA genome have proved useful for measuring genetic variation in numerous insect taxa (Szalanski & Owens 2003; Austin et al. 2002; Taylor et al. 1997; Brower & Jeansonne 2004). Comparison of mtDNA variation with a nuclear genetic variation can provide insight into current versus historical gene flow in a species. For example, high levels of mtDNA variation combined with a lack of nuclear DNA variation may indicate unidirectional mating between strains.
We investigated the extent of genetic variation within and between races of fall armyworm using DNA sequences of a portion of the mitochondrial COI and COII genes, and the nuclear rRNA first internal transcribed spacer (ITS-1) region.
Materials and Methods
Larval fall armyworm samples were collected from sorghum and cotton in Raymond, MS and Colfax, LA, respectively (Table 1). Additional larval samples were obtained from southern Florida and Altheimer, Arkansas, and larval and pupal samples from lab colonies maintained at the University of Mississippi and the University of Florida also were obtained. Larval species identification was confirmed with morphological keys of Peterson (1962), and samples were designated as corn or rice strain based on the host from which they were collected (Table 1). Fall armyworm adults were collected with pheromone traps through summer and fall of 2001 to 2003 from three locations in Arkansas: Tillar, Foreman, and Fayetteville (Table 1). The traps at Tillar were located on the border of experimental research plots of different field crops (cotton, corn, soybean, and sorghum). The adjacent landscape was predominantly cotton with limited acreages of soybean, rice, and corn. A large commercial field of coastal bermudagrass was located within ¼ mile of the traps. The location at Foreman was on a grain farm and the predominant crops were corn, soybean, peanuts, and sorghum. Some limited areas of commercial pasture were near the sample areas. The location at Fayetteville was on an agricultural research farm located in an urban/suburban area. Diverse crops and grasslands were located nearby. Adult fall armyworm identification was confirmed by comparing DNA sequences to larval fall armyworm and other noctuid DNA sequences (unpublished data).
DNA was extracted from individual moths, larvae, and pupae with the Puregene DNA isolation kit D-5000A (Gentra, Minneapolis, MN). Voucher specimens are maintained at the Arthropod Museum, Department of Entomology, University of Arkansas, Fayetteville, AR. DNA vouchers, preserved on filter paper according to Owens & Szalanski (2005), are maintained at the Insect Genetics Laboratory, Department of Entomology, University of Arkansas, Fayetteville, AR.
PCR reactions were conducted with 1 μl of the extracted DNA with New England Biolabs (Ipswich, MA) Taq DNA polymerase with thermopol buffer. Approximately 608 bp of a mtDNA region containing the COI, tRNA leucine, and COII genes was amplified with the primers C1-J-2797 (5'-CCTCGACGTTATTCAGATTACC-3') (Simon et al. 1994) and C2-N-3400 (5'-TCAATATCATTGATGACCAAT-3') (Taylor et al. 1997). The mtDNA marker was amplified with a thermal cycler profile consisting of 35 cycles of 94°C for 45 s, 46°C for 45 s and 72°C for 45 s according to Szalanski et al. (2000). A 562-bp section of the nuclear 3' portion of 18S rDNA, all of ITS-1, and the 5' portion of 5.8S were amplified with the primers rDNA2 (5'-TTGATTACGTCCCTGCCCTTT-3') (Vrain et al. 1992) and rDNA 1.58S (5'-GCCACCTAGTGAGCCGAGCA-3') (Cherry et al. 1997) with a thermal cycler profile consisting of 40 cycles of 94°C for 45 s, 53°C for 1 min and 72°C for 1 min as described by Szalanski & Owens (2003). Amplified DNA from individual moths was purified and concentrated with minicolumns according to the manufacturer's instructions (Wizard PCRpreps, Promega). Samples were sent to The University of Arkansas Medical School DNA Sequencing Facility (Little Rock, AR) for direct sequencing in both directions.
Consensus sequences were derived from both of DNA sequences from an individual with Bioedit 5.09 (Hall 1999) to verify nucleotide polymorphisms, and sequences were aligned by CLUSTAL W (Thompson et al. 1994) for both mtDNA and nDNA sequences. Mitochondrial DNA haplotypes were aligned by MacClade v4 (Sinauer Associates, Sunderland, MA). GenBank accession numbers were AY714298 to AY714304 for the different fall armyworm haplotypes. Genealogical relationships among mtDNA haplotypes were constructed with TCS (Clement et al. 2000) and the method described by Templeton et al. (1992). The distance matrix option of PAUP* 4.0b10 (Swofford 2001) was used to calculate genetic distances according to the Kimura 2-parameter model of sequence evolution (Kimura 1980). Tests for population differentiation were conducted by AMOVA as implemented in Arlequin v. 2.0 (Schneider et al. 2000). An analog of FST, ΦST, was calculated from the haplotypes frequencies and Tajima and Nei (1984) genetic distances (Excoffier et al. 1992). Initially, AMOVA was used to test mtDNA genetic differentiation among all Arkansas and Florida populations sampled (ΦST). Subsequently, a hierarchical AMOVA was conducted in which populations were grouped into states to determine differentiation between states (ΦCT) and among populations within states (ΦST). Pairwise comparisons, calculated independently for all Arkansas and Florida population pairs of ΦST also were calculated. Permutations of the data set were used to determine statistical significance of the pairwise comparisons (P < 0.05).
Mitochondrial DNA sequencing of 71 fall armyworm samples revealed an amplicon size of 608 bp. Nucleotide positions 1 to 222 were COI, 223 to 289 tRNA-leu, and 290 to 608 COII. The average base frequencies were A = 0.36, C = 0.13, G = 0.09, and T = 0.42. Corn haplotype C1 was the most common haplotype for the corn strain and occurred in all of the sampled locations where the corn strain occurred (Table 1). The other two corn strain haplotypes were found only in Arkansas. Rice strain haplotype R2 was the most common haplotype and occurred in every location where the rice strain was found. Rice strains R1 and R4 were found only in Florida, while strain R3 was found in both Arkansas and Florida.
Nine nucleotide sites were variable among the observed three corn and four rice strain haplotypes (Table 2). Three variable nucleotide sites were located in the COI gene and the remainder were located in the COII gene. Tajima-Nei distances (Tajima & Nei 1984) among the fall armyworm haplotypes ranged from 0.164 to 0.329% for the corn strain, 0.164 to 0.329% for the rice strain, and 0.658 to 0.987% between strains. Fig. 1 shows the 95% parsimony network for the seven haplotypes (Posada & Crandall 2001). Missing haplotypes probably represent sampling gaps.
DNA sequencing of the nuclear ITS-1 region from 17 FAW samples (Table 3) revealed an amplicon size of 562 base pairs. No sequence variation was detected in any of the 17 individuals and the base frequencies were A = 0.23, C = 0.24, G = 0.26, and T = 0.27.
AMOVA detected a significant overall ΦST (0.493, P < 0.001) when comparing mtDNA genetic variation among populations (Table 4). The amount of variation was almost equal within versus among populations (within 50.74%, among 49.26%). Hierarchical AMOVA conducted between Arkansas and Florida populations detected a significant ΦCT (0.387, P < 0.005) between the two states (Table 4). The comparison among groups accounted for 38.70% of the observed variation.
This genetic investigation of the fall armyworm mtDNA revealed significant levels of genetic differentiation among populations both within and between the two fall armyworm strains. This research also represents the first attempt to determine the geographical distribution of fall armyworm haplotypes from mtDNA sequence data as well as determining the extent of genetic variation within each strain. A haplotype or allele is defined by one unique form of the gene and differs from any other gene by at least one nucleotide. Haplotype diversity or gene diversity quantifies the number of haplotypes in relation to their relative frequency to each other, and haplotype diversity is described as the probability that two sequences randomly selected from a population are different (Nei 1987).
Four haplotypes were observed for the rice strain and three haplotypes were found for the corn strain, although it is likely more haplotypes may be discovered for each strain. Observed genetic variation between strains was approximately 0.66%. Estimated time of divergence between corn and rice strain is approximately 287,000 years based on a molecular clock rate of 2.3% divergence per million years (Brower 1994). Populations of nearly all species, social or otherwise, exhibit at least some degree of genetic differentiation among geographic locales (Ehrlich & Raven 1969). This observation becomes more difficult to accurately discern when dealing with a migratory species such as the fall armyworm; however, more studies such as this one could help determine the migratory paths of the insect.
One of the purposes of the research presented herein was to estimate the baseline genetic variation which occurs both within and between fall armyworm strains. As with other animal populations, additional genetic structure normally is to be expected over increasing spatial scales, where populations can show additional differentiation due to spatial habitat structure, isolation by distance, or other factors (Avise 1994). There may be temporal differences in the occurrence of the rice and corn. Temporal data, obtained by sampling the same area throughout a season and over a period of years, also may provide insight into the specific migratory patterns of the fall armyworm.
Comparing mtDNA sequences with nuclear markers can provide evidence of inter-strain mating within a species. The lack of variation in the nuclear rDNA ITS-1 region combined with previously conducted laboratory-based mating studies (Pashley & Martin 1987; Whitford et al. 1988; Nagoshi & Meagher 2003) suggests that inter-strain mating does occur in the field. However, the lack of genetic variation in the rDNA ITS-1 region must be approached with caution, because this marker has no power to detect gene flow and this invariant region may be ancestral to strain subdivision. Prowell et al. (2004) also reported a lack of variation in the ITS-1 region of the fall armyworm, but it was cited as unpublished data.
Based on this study, there appears to be sufficient genetic variation both within and between populations to substantiate a more comprehensive population genetics study on this species, and we would recommend also that temporal data be taken into consideration.
We thank Ralph Bagwell, Don Parker, and Clint Allen for providing samples, and John W. Jones for technical assistance. Research was supported by a University of Arkansas, Arkansas Agricultural Experiment Station research initiation grant.
- J. Adamczyk 1993. The UBS DNA Isolation Kit is a reliable method for extracting total nucleic acid from insect tissue. USBiochemical Comments 20:25–26. Google Scholar
- J. W. Austin, A. L. Szalanski, P. Uva, A. Bagneres, and A. Kence . 2002. A comparative genetic analysis of the subterranean termite genus Reticulitermes (Isoptera: Rhinotermitidae). Ann. Entomol. Soc. Amer 95:753–760. Google Scholar
- J. C. Avise 1994. Molecular Markers, Natural History and Evolution. Chapman & Hall, New York, NY. 511 pp. Google Scholar
- A. V Z. Brower 1994. Rapid morphological radiation and convergence among races of the butterfly Heliconuius erato inferred from patterns of mitochondrial DNA evolution. Proc. Natl. Acad. Sci. USA 91:6491–6495. Google Scholar
- A. V Z. Brower and M. M. Jeansonne . 2004. Geographical populations and “subspecies” of new world monarch butterflies (Nymphalidae) share a recent origin and are not phylogenetically distinct. Ann. Entomol. Soc. Amer 97:519–527. Google Scholar
- T. Cherry, A. L. Szalanski, T. C. Todd, and T. O. Powers . 1997. The internal transcribed spacer region of Belonolaimus (Nemata: Belonolaimidae). J. Nematol 29:21–29. Google Scholar
- M. Clement, D. Posada, and K. A. Crandall . 2000. TCS: a computer program to estimate 6 genealogies. Mol. Ecol 9:1657–1659. Google Scholar
- P. R. Ehrlich and P. H. Raven . 1969. Differentiation of populations. Science 165:1228–1232. Google Scholar
- L. Excoffier, P. E. Smouse, and J. M. Quattro . 1992. Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics 131:479–491. Google Scholar
- T. A. Hall 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl. Acids. Symp. Ser 41:95–98. Google Scholar
- M. Kimura 1980. A simple method for estimating evolutionary rate of base substitutions through comparative study of nucleotide sequences. J. Molec. Evol 16:111–120. Google Scholar
- E. F. Knipling 1980. Regional management of the fall armyworm--a realistic approach? Florida Entomol 63:468–480. Google Scholar
- H. C. Levy, A. Garcia-Maruniak, and J. E. Maruniak . 2002. Strain identification of Spodoptera frugiperda (Lepidoptera: Noctuidae) insects and cell line: PCR-RFLP of cytochrome oxidase subunit I gene. Florida Entomol 85:186–190. Google Scholar
- Y. Lu and M. J. Adang . 1996. Distinguishing fall armyworm (Lepidoptera: Noctuidae) strains using a diagnostic mitochondrial DNA marker. Florida Entomol 79 49–55. Google Scholar
- Y. Lu, M. J. Adang, D. J. Eisenhour, and G. D. Kochert . 1992. RFLP analysis of genetic variation in North American populations of the fall armyworm moth Spodoptera frugiperda (Lepidoptera: Noctuidae). Mol. Ecology 1:199–208. Google Scholar
- M. McMichael and D. P. Prowell . 1999. Differences in amplified fragment-length polymorphisms in fall armyworm (Lepidoptera: Noctuidae) host strains. Ann. Entomol. Soc. Amer 92:175–181. Google Scholar
- R. L. Meagher Jr. and M. Gallo-Meagher . 2003. Identifying host strains of fall armyworm (Lepidoptera: Noctuidae) in Florida using mitochondrial markers. Florida Entomol 86:450–455. Google Scholar
- R. N. Nagoshi and R. L. Meagher . 2003. FR tandem-repeat sequence in fall armyworm (Lepidoptera: Noctuidae) host strains. Ann. Entomol. Soc Amer 96:329–335. Google Scholar
- M. Nei 1987. Molecular evolutionary genetics. Columbia University Press, New York. Google Scholar
- C. Owens and A. L. Szalanski . 2005. Filter paper for preservation, storage and distribution of insect and pathogen DNA samples. J. Med. Entomol. (in press). Google Scholar
- D. P. Pashley 1986. Host associated genetic differentiation in fall armyworm (Lepidoptera: Noctuidae): a sibling species complex? Ann. Entomol. Soc. Amer 79:898–904. Google Scholar
- D. P. Pashley 1988a. Quantitative genetics, development and physiological adaptation in sympatric host strains of fall armyworm. Evolution 42:93–102. Google Scholar
- D. P. Pashley 1988b. Current status of fall armyworm host strains. Florida Entomol 71:227–234. Google Scholar
- D. P. Pashley 1989. Host-associated differentiation in armyworms (Lepidoptera: Noctuidae): An allozymic and mitochondrial DNA perspective. pp. 103-114 In H. D. Loxdale and J. den Hollander [eds.], Systematics Association Special Volume No. 39. Claredon Press, Oxford. Google Scholar
- D. P. Pashley and J. A. Martin . 1987. Reproductive incompatibility between host strains of the fall armyworm (Lepidoptera: Noctuidae). Ann. Entomol. Soc. Amer 80:731–733. Google Scholar
- D. P. Pashley, S. J. Johnson, and A. N. Sparks . 1985. Genetic population structure of migratory moths: the fall armyworm (Lepidoptera: Noctuidae). Ann. Entomol. Soc. Amer 78:756–762. Google Scholar
- A. Peterson 1962. Larvae of Insects. Edward Bros., Inc., Ann Arbor, MI. Google Scholar
- D. Posada and K. A. Crandall . 2001. Performance of methods for detecting recombination from DNA sequences: computer simulations. Proc. Natl. Acad. Sci. USA 98:13757–13762. Google Scholar
- D. P. Prowell, M. McMichael, and J. F. Silvain . 2004. Multilocus genetic analysis of host use, introgression, and speciation in host strains of fall armyworm (Lepidoptera: Noctuidae). Ann. Entomol. Soc. Amer 97:1034–1044. Google Scholar
- S. Schneider, D. Roessli, and L. Excoffier . 2000. Arlequin: A software for population genetic data. Genetics and Biometry Laboratory, University of Geneva, Switzerland. Google Scholar
- C. Simon, F. Frati, A. Beckenbach, B. Crespi, H. Liu, and P. Flook . 1994. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann. Entomol. Soc. Amer 87:651–701. Google Scholar
- A. N. Sparks 1979. A review of the biology of the fall armyworm. Florida Entomol 62:82–87. Google Scholar
- D. L. Swofford 2001. PAUP*: Phylogenetic analysis using parsimony (*and other methods), ver. 4.0b10. Sinauer, Sunderland, MA. Google Scholar
- A. L. Szalanski and C. B. Owens . 2003. Genetic variation of the southern corn rootworm, Diabrotica undecimpunctata howardi (Coleoptera: Chrysomelidae). Florida Entomol 86:329–333. Google Scholar
- A. L. Szalanski, D. S. Sikes, R. Bischof, and M. Fritz . 2000. Population genetics and phylogenetics of the endangered American burying beetle, Nicrophorus americanus (Coleoptera: Silphidae). Ann. Entomol. Soc. Amer 93:589–594. Google Scholar
- F. Tajima and M. Nei . 1984. Estimation of evolutionary distance between nucleotide sequences. Mol. Biol. Evol 1:269–285. Google Scholar
- D. B. Taylor, R. D. Peterson II, A. L. Szalanski, and J. J. Petersen . 1997. Mitochondrial DNA variation among Muscidifurax spp. (Hymenoptera: Pteromalidae), pupal parasitoids of filth flies. Ann. Entomol. Soc. Am 90:814–824. Google Scholar
- A. R. Templeton, K. A. Crandall, and C. F. Sing . 1992. A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. Cladogram estimation. Genetics 132:619–633. Google Scholar
- J. D. Thompson, D. G. Higgins, and T. J. Gibson . 1994. CLUSTAL W: improving the sensitivity of progressive multiples sequence alignments through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 22:4673–4680. Google Scholar
- T. C. Vrain, D. C. Wakarchuk, A. C. Levesque, and R. I. Hamilton . 1992. Intraspecific rDNA restriction fragment length polymorphism in the Xiphinema americanum group. Fundam. Appl. Nematol 15:563–573. Google Scholar
- F. Whitford, S. S. Quisenberry, J. J. Riley, and J. W. Lee . 1988. Oviposition preference, mating compatibility and development of two fall armyworm strains. Florida Entomol 71:234–243. Google Scholar
Sampling locations and frequency of fall armyworm corn “C” and rice “R” haplotypes.
Genetic variation at nine nucleotide sites among fall armyworm haplotypes.
Sampling locations of fall armyworm with nuclear rDNA ITS-1 marker.
Estimates of genetic differentiation calculated among Arkansas and Florida populations as one hierarchical group (One Group) and with an additional hierarchical group between states (Two Groups). In both, the amount of variation occurring among categories and the estimate of genetic differentiation is provided. An asterisk denotes statistical significance at α = 0.05.