Biocontrol Science and Technology, August 2005; 15(5): 481/495

A single-step multiplex PCR assay for the detection of European Peristenus spp., parasitoids of Lygus spp.

TARA D. GARIEPY13, ULRICH KUHLMANN2, TIM HAYE2, CEDRIC GILLOTT3, & MARTIN ERLANDSON1

1Agriculture and Agri-Food Canada, Saskatoon Research Centre, Saskatoon, SK, Canada, 2CABI Bioscience Switzerland Centre, Dele´mont, Switzerland, and 3Department of Biology, University of Saskatchewan, Saskatoon, SK, Canada

(Received 8 August 2004; returned 14 October 2004; accepted 13 November 2004)

Abstract Lygus Hahn (Hemiptera: Miridae) are serious pests of agricultural and greenhouse crops throughout North America. In Europe, bivoltine Peristenus Fo¨rster (: ) species have a significant impact on Lygus populations. Release and establishment of European P. digoneutis Loan in Lygus lineolaris (Palisot de Beauvois) populations in northeastern USA has renewed interest in the intended liberation of European parasitoids for Lygus control in Canada. Accurate identification of natural enemies is the cornerstone of biological control but conventional methods for identifying Peristenus species and estimating parasitism rates rely on tedious and time-consuming dissection and rearing methods. The present study describes species-specific PCR primers for three species of Peristenus, and the use of a multiplex PCR assay to detect P.digoneutis and P.stygicus Loan eggs and larvae from Lygus rugulipennis Poppius nymphs. Results indicate that the primers amplify uniquely sized, species-specific PCR products for the three species and are capable of detecting single eggs in parasitized nymphs within 3 days post-parasitism. Using a multiplex PCR assay, the primers maintain specificity and sensitivity, and allow detection of each of the three species in a single reaction. Although molecular diagnostics have previously been used in the identification of parasitoids and estimation of parasitism rates, this is the first time a single-step multiplex PCR protocol has been described.

Keywords: Peristenus digoneutis, P. stygicus, P. pallipes, Lygus, biological control, molecular diagnostics, multiplex PCR

Introduction Species of Lygus Hahn (Hemiptera: Miridae) are serious pests of crops grown in agricultural and greenhouse environments throughout North America. Feeding damage by Lygus nymphs and adults can cause flower bud and fruit abortion, side- shoot abscission, leaf tissue perforation and deformation or death of meristem tissue (Strong 1970; Broadbent et al. 2002). In North America, Lygus have historically been recorded as destructive pests (Craig & Loan 1984; Young 1986; Coulson 1987), and recent literature demonstrates their continued pest status on a variety of crops,

Correspondence: M. Erlandson, Agriculture and Agri-Food Canada, Saskatoon Research Centre, Saskatoon, SK, Canada S7N 0X2. Tel: 1 306 956 7276. Fax: 1 306 956 7247. E-mail: erlandsonm@ agr.gc.ca

ISSN 0958-3157 print/ISSN 1360-0478 online # 2005 Taylor & Francis Group Ltd DOI: 10.1080/09583150500086771 482 T. D. Gariepy et al. including wheat (Wise et al. 2000), oilseed flax (Wise & Lamb 2000), cotton (Layton 2000), canola (Boyd & Lentz 1999; Braun et al. 2001; Carcamo et al. 2002), alfalfa (Braun et al. 2001; May et al. 2003;) and strawberry (Rhainds & English-Loeb 2003). In North America, native Lygus parasitoids in the genus Peristenus Fo¨rster (Hymenoptera: Braconidae) are primarily univoltine, attacking only the first Lygus generation. As Lygus species are generally multivoltine, subsequent generations are not attacked by native parasitoids (Kuhlmann et al. 1998; Coutinot & Hoelmer 1999). Several bivoltine and multivoltine Peristenus species have a significant impact on populations of Lygus rugulipennis Poppius in Europe (Loan & Bilewicz-Pawinska 1973; Bilewicz-Pawinska 1976; Kuhlmann et al. 1998). Thus, the introduction of bivoltine European Peristenus species as classical biological control agents into Canada may fill an unoccupied niche and suppress native Lygus populations throughout the season. Attempts have been made since the 1960s to introduce and establish populations of P. digoneutis Loan and P. stygicus Loan in the USA (Coulson 1987). Peristenus digoneutis was released in the 1980s in New Jersey and has become established in Lygus lineolaris (Palisot de Beauvois) populations in alfalfa and strawberry fields in northeastern USA (Day 1999; Day et al. 1990, 2000; Tilmon & Hoffmann 2003). The success of this program has renewed interest in the intended liberation of European parasitoids for Lygus control in Canada. Accurate identification of potential natural enemies is the cornerstone of biological control, and methods that can definitively identify potential biological control agents are essential, especially when morphological variation among species is slight. Members of the genus Peristenus can be difficult to distinguish because species differences are small (Loan & Bilewicz-Pawinska 1973; Bilewicz-Pawinska & Pankanin 1974). As such, PCR-based methods for distinguishing between species would be advantageous for identifying adult Peristenus species and for detecting immature stages of Peristenus in parasitized Lygus nymphs. Recently, molecular techniques have been investigated for their utility in identifica- tion of some species of Peristenus (Tilmon et al. 2000; Erlandson et al. 2003; Ashfaq et al. 2004; Zhu et al. 2004). These molecular marker systems include PCR-based methods that potentially preclude the need for tedious and time-consuming dissection and rearing procedures. This would be advantageous as dissection provides no information as to which parasitoid species is present (Carignan et al. 1995), and rearing methods often require several months because of diapause prior to adult emergence and identification (Loan & Bilewicz-Pawinska 1973). Tilmon et al. (2000) developed molecular markers based on the cytochrome oxidase 1 (CO1) gene for P. pallipes Curtis, P. digoneutis, and P. conradi Marsh in order to assess parasitism rates in L. lineolaris in the field. Similarly, Erlandson et al. (2003) developed molecular markers based on sequences from the internal transcribed spacer (ITS) regions for P. digoneutis, P. stygicus, and P. pallipes. However, the molecular markers developed in both these studies amplify fragment sizes that are identical for all Peristenus species investigated, and therefore require the additional step of restriction enzyme digestion or additional PCRs in order to distinguish between species. Recently, Zhu et al. (2004) screened a collection of single, short oligonucleotide primers and selected several which produced PCR product patterns that distinguish between P. pallipes, Single-step multiplex PCR for detection of European Peristenus spp. 483

P. pseudopallipes and P. howardi without the need for additional digestion with restriction enzymes. Multiplex PCR involves the use of multiple species-specific primer pairs in a single reaction to amplify unique fragment sizes for different DNA targets, allowing samples to be screened for several organisms simultaneously (Bej et al. 1991; Henson & French 1993). Although this type of system has yet to be implemented in host- parasitoid studies to identify multiple parasitoid species, multiplex PCR is routinely used in the diagnosis of medical and veterinary diseases (Rossiter & Caskey 1993; Zarlena & Higgins 2001; Favia et al. 2003) as well as in the detection of plant pathogens (Hamelin et al. 1996; Cullen et al. 2000; Taylor et al. 2001; Bertolini et al. 2003; Gariepy et al. 2003; Lopez et al. 2003). The purpose of this study was to design PCR primers for P. digoneutis, P. stygicus, and P. pallipes that produce unique fragment sizes, and that could be applied in multiplex to improve their diagnostic utility. Two of these species, P. digoneutis and P. stygicus, have been released as biological control agents for Lygus spp. in the USA (Coulson 1987; Day 1999), and are currently being considered for release in Canada. Peristenus pallipes, however, is an holarctic species occuring in Europe and North America (Loan 1965, 1974; Loan & Bilewicz-Pawinska 1973; Hormchan 1977) where it is known to parasitize several mirid species (Craig 1963; Loan 1965; Bilewicz- Pawinka 1977, 1982). A multiplex PCR protocol for these species would be particularly useful in North America and Europe for detecting immature stages of Peristenus in parasitized Lygus nymphs in order to determine Peristenus-host associations, which are particularly important for non-target risk assessment studies.

Materials and methods Acquisition of voucher specimens Cocoons of the Lygus parasitoids P. digoneutis and P. stygicus were reared from L. rugulipennis nymphs collected from geographically dispersed red clover (Trifolium pratense L.) and camomile (Matricaria recutita L.) fields in Schleswig-Holstein, northern Germany. Peristenus pallipes cocoons were obtained from Liocoris tripustulatus Fabricius (Hemiptera: Miridae) and Closterotomus norvegicus Gmelin (Hemiptera: Miridae) collected from stinging nettles (Urtica dioica L.) in the same region. Emerged adults were put in 95% ethanol and identified based on morphological characteristics by H. Goulet (Systematic Entomology Unit, Agriculture and Agri-Food Canada, Ottawa, Canada). Twenty identified specimens for each of the Peristenus species investigated were randomly selected from collections made in different habitats and locations in northern Germany. These specimens were used as voucher samples for molecular analysis, and were kept in 95% ethanol.

DNA extraction Peristenus adults preserved in 95% ethanol were rinsed in sterile, distilled water, and placed individually in 1.5-mL centrifuge tubes. Individual Peristenus adults were homogenized using disposable pestles (Kimble Kontes, Vineland, NJ, USA) in 250 mL of Lifton buffer [0.2 M sucrose, 50 mM ethylene-diaminetetraacetic acid, 100 mM

Tris /HCl (pH 7.5), and 0.5% sodium dodecyl sulfate], and incubated for 2 h at room temperature. Proteinase K (1.25 mLof20mg/mL) was added, and samples were 484 T. D. Gariepy et al. incubated for 4 h at 658C. Samples were put on ice for 30 min following the addition of 35 mL of 5 M potassium acetate, and centrifuged for 20 min at 20 000/g.

The supernatant was extracted once with 200 mL of phenol/chloroform/isoamyl alcohol (25:24:1), and once with 200 mL chloroform/isoamyl alcohol (24:1). The aqueous phase was precipitated with two volumes of cold 95% ethanol by overnight incubation at /208C. Precipitated DNA was pelleted by centrifugation (30 min at 20 000/g), washed with cold 70% ethanol, and air-dried for 15 min. DNA samples were resuspended in 50 mLTris/EDTA (pH 7) with 0.5 mL of 10 mg/mL RNase A.

DNA sequencing and primer design The oligonucleotide primers inDNA44 (5?-TCCTCCGCTTATTGATATGC-3?) and inDNA45 (5?-GGAAGTAAAAGTCGTAACAAGG-3?) (Nucleic Acid / Protein Service, University of British Columbia, Vancouver, BC, Canada) were used in PCR to amplify conserved rRNA gene sequences and ITS regions of nuclear DNA from voucher specimens. PCR conditions were as follows: 25-mL reactions contained

2.5 mLof10/ reaction buffer, 2.5 mM MgCl2, 0.25 mM dNTPs, 0.5 mM of each primer, 1 U of Ta q DNA polymerase (Invitrogen, Burlington, ON, Canada), and approximately 20 ng of DNA template. PCR reactions were carried out in a Strategene Robocycler thermocycler (LaJolla, CA, USA) with an initial denaturation at 948C for 120 s, followed by 30 cycles at 948C for 60 s, 558C for 60 s, and 728C for 90 s. A final extension period at 728C for 5 min followed the 30 cycles. PCR products were electrophoresed on 1% agarose gels containing 0.5 mg/mL ethidium bromide at

75 V for 45/60 min and products were visualized by UV transillumination. PCR † products were directly cloned into pGem / T Easy vectors (Promega, Madison, WI, USA) using the manufacturer’s protocol. Cloned PCR products were sequenced at the DNA Services Laboratory, Plant Biotechnology Institute (NRCC, Saskatoon, SK, Canada) using universal primers to sequence from both ends of the cloned insert. DNA sequence data were assembled and aligned using SeqMan (DNAStar Inc., Madison, WI, USA) and Vector NTI (InfoMax, Fredrick, MD, USA). Areas of sequence variation determined by visual inspection of the sequence data were used to design primers specific for P. digoneutis, P. stygicus and P. pallipes.

Specificity of Peristenus species-specific primers The specificity of the primer pairs designed for P. digoneutis, P. stygicus and P. pallipes was evaluated by amplifying 20 voucher specimens of each species using each of the primer pairs. Adult Lygus DNA was used as a negative control to ensure that the primer sets did not amplify Lygus sequences. Amplification of DNA was performed in a Stratagene Robocycler thermocycler using the same reagents as above, except the primer pairs used were Per R1 (a presumed universal Peristenus primer) in combination with either dig F1096 for P. digoneutis, sty F1230 for P. stygicus, or pal F517 for P. pallipes. Primer concentrations were as described above. Amplification conditions were: initial denaturation at 948C for 120 s, followed by 35 cycles of 948C for 45 s, 548C for 45 s, and 728C for 60 s. A final extension period of 5 min at 728C followed. Single-step multiplex PCR for detection of European Peristenus spp. 485

Multiplex PCR assay and detection sensitivity In order to enhance the diagnostic utility of Peristenus spp. primers, a multiplex PCR protocol was tested. For this approach, the reagents mentioned above were used with the following modifications: 0.2 mM dig F1096, 0.2 mM sty F1230, 0.2 mM pal F517 and 0.4 mM Per R1 were combined in the reaction mixture. Amplification conditions were as described above for specificity testing of the Peristenus species- specific primers. DNA from adult P. digoneutis, P. stygicus, and P. pallipes voucher specimens was serially diluted 10-fold, from 20 to 0.02 ng, and used in a multiplex PCR assay to determine the sensitivity of the primers.

Detection of immature stages of Peristenus in parasitized Lygus nymphs To determine the ability of the molecular markers to detect immature stages of the biological control agents, P. digoneutis and P. stygicus inside their host, L. rugulipennis nymphs were parasitized in the laboratory. Between 20 and 25 experienced, mated P. digoneutis and P. stygicus females (5/10 days old) were provided with second or third instar Lygus nymphs in a 5-cm diameter Petri dish. Each wasp was transferred into a Petri dish containing five nymphs, and attacks were observed. When an attack occurred, the nymph was immediately removed to prevent superparasitism. In total, 100/125 nymphs attacked by P. digoneutis and

100/125 nymphs attacked by P. stygicus were obtained. Attacked nymphs were either placed immediately in 95% ethanol (0 h) or reared for fixed periods at approximately 229/28C, 16L:8D photoperiod, and 50% RH in rearing cylinders with sprouting potatoes, lettuce, and green beans. Rearing periods were chosen to cover the documented range of developmental stages of P. digoneutis and P. stygicus within Lygus spp. nymphs (Drea et al. 1973; Carignan et al. 1995). Rearing periods included: 3 days (egg), 6 days (first instar), 9 days (second instar), and 12 days (third instar). After the specified rearing period, attacked nymphs were placed in vials containing 95% ethanol. DNA was extracted from the attacked nymphs (20 nymphs per rearing period for each Peristenus species), and amplified using the multiplex PCR assay.

Application of the multiplex PCR assay to field-collected Lygus nymphs Lygus rugulipennis nymphs were collected from sweep-net samples from three different sites [two camomile, Matricaria recutita L. (GPS coordinates: N 54824.757? E 010809.919?;N54825.484? E 010801.349?) and one red clover, Trifolium pratense L. (GPS coordinates: N 54805.191? E 010804.353?)] in Schleswig-Holstein, Germany. From each site, 300 late-instar Lygus nymphs were collected: 100 nymphs from each site were reared, 100 from each site were dissected, and 100 from each site were used for molecular analysis using the multiplex PCR assay. DNA extraction and amplification using the multiplex assay was as described above, and PCR products from individual nymphs were electrophoresed at 65 V for 60 min on 1% agarose gels containing 0.5 mL/mL ethidium bromide. Parasitism and/or species composition based on rearing, dissection and molecular analysis were compared. 486 T. D. Gariepy et al.

Results DNA sequencing and primer design DNA sequences were obtained for cloned PCR products generated using inDNA44 and inDNA45 from at least four individuals (one to three clones per individual) for each of the Peristenus species investigated. Sequences were obtained for the ITS1 regions of P. digoneutis (GenBank accession no. AY605234) and P. stygicus (GenBank accession no. AY605233). ITS1 and ITS2 sequences were obtained for P. pallipes (GenBank accession no. AY608602). Sequences were aligned to analyse sequence variation among species. Several areas of variation were identified, allowing the design of species-specific primers for each species (Table I). Primer sets consisted of a species-specific reverse primer (dig F1096, sty F1230, or pal F517) used in combination with a genus-specific forward primer (Per R1) common to all three species investigated. The combination of Per R1 with dig F1096 for P. digoneutis, with sty F1230 for P. stygicus, and with pal F517 for P. pallipes produced PCR products of 515, 330, and 1060 bp, respectively (Figure 1).

Specificity of Peristenus species-specific primers Each primer set selectively amplified the DNA of the species for which it was designed, and yielded the expected product sizes when using the PCR conditions described above. Species-specific amplification was consistent for all voucher speci- mens tested (20 specimens per species), and a representative for each species is shown (Figure 2). In addition to the 20 voucher specimens of P.pallipes from European mirid populations, the P. pallipes primers (Per R1 and palF517) also amplified P. pallipes (25 individuals) obtained from three widely dispersed geographic regions in western Canada. All samples yielded the expected 1060-bp PCR product for P. pallipes (Erlandson, Ashfaq and Gariepy, data not shown), indicating that the primers designed for European P. pallipes were also specific for representative North American populations. None of the primer sets amplified Lygus DNA.

Multiplex PCR assay and detection sensitivity Combining the three species-specific reverse primers (dig F1096, sty F1230, pal F517) and the genus-specific forward primer (Per R1) in a multiplex PCR protocol distinguished each of the three species, and amplified the appropriate species-specific fragment; 515 bp for P. digoneutis, 330 bp for P. stygicus, and 1060 bp for P. pallipes (Figure 3). These results were consistent with those obtained using the species- specific primer sets individually (Figure 2). When used in multiplex, the lower limit of

Table I. Peristenus species-specific primers, their respective sequences, and fragment sizes when used in combination with a genus-specific reverse primer (Per R1).

Species Primer Primer sequence Fragment size

Peristenus spp. Per R1 5?-ACAAGGTTTCCGTAGGTG-3?/ P. digoneutis dig F1096 5?-GAACATAAAAACTTCTTCTACGC-3? 515 P. stygicus sty F1230 5?-CAGGTAGAGATACATGGCTGT-3? 330 P. pallipes pal F517 5?-TAAACTTTGGCCAGATAAATG-3? 1060 Single-step multiplex PCR for detection of European Peristenus spp. 487

Figure 1. Diagrammatic representation of the ITS region of P. digoneutis, P. stygicus and P. pallipes, showing location of species-specific primers. detection appears to be 0.2 ng for P. digoneutis primers, and 0.02 ng for P. stygicus and P. pallipes primers (Figure 4).

Detection of immature stages of Peristenus in parasitized Lygus nymphs Using the multiplex PCR protocol, primers successfully detected Peristenus eggs and larvae in L. rugulipennis nymphs 3, 6, 9, and 12 days post-parasitism. For these periods, detection of both parasitoid species was consistent among Lygus nymphs (20 attacked nymphs per rearing period for P. digoneutis and P. stygicus), with 90/100% of observed attacks resulting in confirmed parasitism; representa- tives for each rearing period are presented (Figure 5). Detection of Peristenus within nymphs placed in 95% ethanol immediately after an attack (0 h) was inconsistent. Faint bands of the expected size were occasionally present, but required 4/5 times more DNA template than was required for samples from the other post-parasitism periods. The primer sets did not amplify Lygus DNA.

Figure 2. Diagnostic PCR using species-specific primers for P.digoneutis (dig F1096 and Per R1), P.stygicus (sty F1230 and Per R1), and P.pallipes (pal F517 and Per R1). D, S, and P designate P.digoneutis, P.stygicus, and P. pallipes, respectively. M, 100-bp marker (Invitrogen); NEG, negative control (no DNA). 488 T. D. Gariepy et al.

Figure 3. Diagnostic multiplex PCR of P.digoneutis, P.stygicus and P. pallipes and Lygus adult DNA (NEG, negative control) using a combination of dig F1096, sty F1230, pal F517, and Per R1 in each reaction. M, 100-bp marker (Invitrogen).

Application of the multiplex PCR assay to field-collected Lygus nymphs Results from rearing, dissection, and molecular analysis of Lygus nymphs collected from three sites are summarized in Table II. For each site, the percent parasitism based on molecular analysis, rearing, and dissection were not significantly different 2 2 2 (Chi-square test, site 1: x /0.36, P B/0.05; site 2: x /1.5, P B/0.05; site 3: x /1.02,

P B/0.05). Multiplex PCR results indicated the occurrence of multiparasitism in some Lygus nymphs, as demonstrated by the presence of P. digoneutis and P. stygicus specific PCR products in single nymphs (Figure 6). When the incidence of multiple parasitism detected by PCR analysis was taken into account, the relative proportions of P. digoneutis and P. stygicus in the overall species composition were similar when determined by rearing and multiplex PCR.

Discussion The PCR primers developed for P.digoneutis, P.stygicus,andP.pallipes are specific and sensitive when used individually or in multiplex with DNA extracted from Peristenus adults or from parasitized Lygus nymphs. When used individually or in multiplex, the primers yield unique fragment sizes for each species, allowing the identification of

Figure 4. Multiplex PCR using 20, 2, 0.2, and 0.02 ng of DNA template from P. digoneutis, P. stygicus, and P. pallipes adults. M, 100-bp marker (Invitrogen). Single-step multiplex PCR for detection of European Peristenus spp. 489

Figure 5. Multiplex PCR of Lygus nymphs parasitized by P.digoneutis and P.stygicus and killed 0 h and 3, 6, 9, and 12 days after attack. NEG 1 and NEG 2 are negative controls containing no DNA and Lygus DNA, respectively. M, 100-bp marker (Invitrogen).

Peristenus eggs (3 days post-parasitism), larvae, and adults based on the size of the PCR fragment generated. The inability of molecular markers to detect newly laid parasitoid eggs in their host is not uncommon. Primers developed by Amornsak et al. (1998) for detection of Trichogramma australicum Girault in Helicoverpa armigera Hubner and H. punctigera Wallengren were unable to detect the parasitoid immediately after oviposition; PCR products were only visible 12 h after parasitism. Within the first 24 h after oviposition, Peristenus eggs swell and further embryonic development occurs (Carignan et al.

1995). Presumably within 12/24 h post-parasitism, our multiplex PCR system would be capable of detecting eggs, but this remains to be tested. The fact that electrophoresis of PCR product from the 3-day post-parasitism stage produced very strong bands indicates that this is certainly not the lower limit of detection. Amplification of serial dilutions of adult Peristenus DNA suggests that the lower limit of detection is between 0.2 and 0.02 ng for P.digoneutis and approximately 0.02 ng for P. stygicus. However, this represents the ability of the primers to detect pure Peristenus DNA, and may not reflect the minimum quantity detectable in samples containing parasitoid and host DNA. In order to determine the lower-limit of parasitoid detection within Lygus nymphs, it would be prudent to expand these experiments to include 12, 24, and 48 h post-parasitism periods. In addition to being used previously to distinguish between Peristenus species (Tilmon et al. 2000; Erlandson et al. 2003; Zhu et al. 2004), molecular diagnostic methods have been used to separate morphologically similar parasitoid species using PCR or PCR followed by restriction enzyme digestion. For example, these techniques have been used by Silva et al. (1999) and Stouthamer et al. (1999) to distinguish between species of Trichogramma Westwood (Hymenoptera: Trichogrammatidae), by Zhu and Greenstone (1999) and Zhu et al. (2000) to detect different species and strains of Aphelinus Dalman (Hymenoptera: Aphelinidae) in Diuraphis noxia Mordvilko (Homoptera: Aphididae), and by Ratcliffe et al. (2002) to detect eight pteromalid (Hymenoptera: Pteromalidae) species parasitizing Musca domestica L. and Stomoxys calcitrans L. However, none of these studies use techniques that allow detection and identification of more than one species in a single reaction. 490 .D aip tal. et Gariepy D. T.

Table II. Total percent parasitism and parasitoid species composition based on molecular analysis, rearing, and dissection of field-collected Lygus rugulipennis nymphs from northern Germany. From each site, 100 nymphs were tested per method (i.e. 100 for molecular analysis, 100 for rearing, and 100 for dissection).

Site 1 Site 2 Site 3

Molecular Rearing Dissections Molecular Rearing Dissections Molecular Rearing Dissection

Total% parasitism: 68 68 62 70 70 58 60 67 56 Species composition:

P. digoneutis (%) 53 66 / 83 100 / 50 58 / P. stygicus (%) 29 34 / 60/ 47 42 / Multiparasitism (%) 18 / / 11 / / 3 / /

/, information not available using this technique. Single-step multiplex PCR for detection of European Peristenus spp. 491

Figure 6. Application of the multiplex PCR assay to field collected Lygus rugulipennis nymphs. PCR products specific to P. digoneutis (515 bp) and P. stygicus (330 bp) are indicated by D and S, respectively. Lanes with no bands indicate nymphs that are not parasitized (NP), and lanes showing both 515- and 330-bp products simultaneously indicate multiparasitism (DS). M, 100-bp marker (Invitrogen).

Lopez et al. (2003) list multiplex PCR methods among the innovative tools currently used in pathogen diagnosis based on the sensitivity of the technique and the ability to simultaneously detect different DNA targets in a single reaction. In parasitology, molecular diagnostic tools incorporating multiplex PCR have been used to identify different mosquito species and strains that vector malaria (Kengne et al. 2001; Phuc et al. 2003), to detect viruses and plant pathogens in insect vectors (Wilson & Chase 1993; Rodrigues et al. 2003), to identify microsporidian infection in ants (Valles et al. 2002) and bacterial infection in insect cadavers (Bourque et al. 1993). However, the implementation of multiplex assays to detect and identify insect parasitoid species within their hosts has not previously been reported. The application of the multiplex PCR protocol to field-collected Lygus samples demonstrates the applicability and potential of this method for screening field- collected samples for the presence of Peristenus spp. Percent parasitism by rearing, dissection and molecular analysis were not significantly different; however, molecular methods provide additional information pertaining to multiple parasitism. Such information is not evident by rearing or dissection, but could become important to assess inter-specific competition between indigenous and exotic biological control agents in target and non-target hosts. Furthermore, the use of a multiplex PCR assay for Peristenus species improves upon the conventional and molecular techniques currently used to identify Peristenus spp. and quantify parasitism levels. It allows rapid and accurate identification of multiple Peristenus species and eliminates mortality issues encountered during the rearing process. Furthermore, unlike dissections, the multiplex assay allows the determination of species composition, and unlike rearing, it does so without the time-delay due to the lengthy diapause period of the parasitoid prior to adult emergence. This is a novel approach that could be valuable in insect biological control programs, particularly for pre-release and post-release studies on potential biological control agents. Pre-release studies, including the determination of ecological and physiological host ranges, are strongly recommended by van Lenteren et al. (2003) prior to importation and release of exotic parasitoids for inundative biological control programs. Furthermore, recommendations by Louda et al. (2003) suggest that pre-release field research generates better ecological data for biological control agents in order to reduce the risk of non-target effects. Insect molecular systematics has supplemented and improved the value of morphological and ecological data (Caterino et al. 2000), 492 T. D. Gariepy et al. and in fact, molecular data may complement and enhance ecological and physiological studies on parasitoid host range. The single-step multiplex PCR assay described here would be useful in conducting pre-release studies in the area of origin; nymphs can be screened to determine

Peristenus/ host associations as well as the distribution of Peristenus species. In addition, following introduction, molecular tools could facilitate post-release mon- itoring studies. In fact, recent post-introduction studies in New York State using molecular markers for the COI gene have been used to map the distribution of P. digoneutis in Lygus populations in alfalfa and strawberry fields (Tilmon et al. 2000; Tilmon & Hoffmann 2003). However, as previously mentioned, the markers used in these studies produce identical fragment sizes for each species and require additional processing to obtain species-specific identifications. The use of the single-step multiplex PCR protocol described here could facilitate similar post-introduction studies in the future. Furthermore, it would allow samples to be screened for the introduced Peristenus species and indigenous North American P. pallipes simulta- neously to assess establishment, distribution, potential displacement of native P. pallipes, and realized host range in the area of introduction.

Acknowledgements This work was supported by the Biocontrol Network (Canada), CABI Bioscience, and Agriculture and Agri-Food Canada. We thank Dr H. Goulet (AAFC, ECORC, Ottawa, ON, Canada) for voucher specimen identification, and Dr P. Mason (AAFC, ECORC, Ottawa, ON, Canada) for useful discussions and advice.

References Amornsak W, Gordh G, Graham G. 1998. Detecting parasitized eggs with polymerase chain reaction and DNA sequence of Trichogramma australicum Girault (Hymenoptera: Trichogrammatidae). Australian Journal of Entomology 37:174/179. Ashfaq M, Braun L, Hegedus D, Erlandson M. 2004. Estimating parasitism levels in Lygus spp. (Hemiptera: Miridae) field populations using standard and molecular techniques. Biocontrol Science

and Technology 14:731/735. Bej AK, McCarty SC, Atlas RM. 1991. Detection of coliform bacteria and Escherichia coli by multiplex polymerase chain reaction, comparison with defined substrate and plating methods for water quality

monitoring. Applied Environmental Microbiology 57:2429/2432. Bertolini E, Olmos A, Lopez MM, Cambra M. 2003. Multiplex nested reverse transcription-polymerase chain reaction in a single tube for sensitive and simultaneous detection of four RNA viruses and

Pseudomonas savastanoi pv. savastanoi in olive trees. Phytopathology 93:286/292. Bilewicz-Pawinska T. 1976. Natural enemies of some plant bugs (Miridae) in Poland (English Summary). Roczniki Nauk Rolniczych 6:125/135. Bilewicz-Pawinska T. 1977. Parasitism of Adelphocoris lineolatus Goeze and Lygus rugulipennis Popp.

(Heteroptera) by braconids and their occurence on alfalfa. Ekologia Polska 25:539/550. Bilewicz-Pawinska T. 1982. Plant bugs (Heteroptera: Miridae) and their parasitoids (Hymenoptera: Braconidae) on cereal crops. Polish Ecological Studies 8:113/191. Bilewicz-Pawinska T, Pankanin M. 1974. Larvae of Peristenus Foerster (Hym., Braconidae) parasitic on Lygus rugulipennis Popp. (Heteroptera: Miridae) (English Summary). Polskie Pismo Entomologiczne 44:759/764. Bourque SN, Valero JR, Mercier J, Lavoie MC, Levesque RC. 1993. Multiplex polymerase chain reaction for detection and differentiation of the microbial insecticide Bacillus thuringiensis. Applied Environmental

Microbiology 59:523/527. Single-step multiplex PCR for detection of European Peristenus spp. 493

Boyd ML, Lentz GL. 1999. Seasonal occurence and abundance of the tarnished plant bug (Hemiptera: Miridae) and thrips (Thysanoptera: Thripidae) on rapeseed in West Tennessee. Journal of Agricultural

and Urban Entomology 16:171/178. Braun L, Erlandson M, Baldwin D, Soroka J, Mason P, Foottit R, Hegedus D. 2001. Seasonal occurrence, species composition, and parasitism of Lygus spp. in alfalfa, canola, and mustard. Canadian Entomologist

133:565/577. Broadbent AB, Mason PG, Lachance S, Whistlecraft JW, Soroka JJ, Kuhlmann U. 2002. Lygus spp. plant bugs (Hemiptera: Miridae). In: Mason PG, Huber JT, editors. Biological control programmes in Canada

1981/2000. Wallingford: CAB International Publishing. pp 152/159. Carcamo H, Otani J, Herle C, Dolinski M, Dosdall L, Mason P, Butts R, Kaminski L, Olfert O. 2002. Variation of Lygus species assemblages in canola agroecosystems in relation to ecoregion and crop stage.

Canadian Entomologist 134:97/111. Carignan S, Boivin G, Stewart RK. 1995. Developmental biology and morphology of Peristenus digoneutis

Loan (Hymenoptera: Braconidae: ). Biological Control 5:553/560. Caterino MS, Cho S, Sperling FAH. 2000. The current state of insect molecular systematics: a thriving

tower of Babel. Annual Review of Entomology 45:1/54. Coulson JR. 1987. Studies on the biological control of plant bugs (Heteroptera: Miridae): An introduction

and history, 1961/1983. In: Hedlund RC, Graham HM, editors. Economic importance and biological control of Lygus and Adelphocoris in North America. Maryland: United States Department of

Agriculture, Agricultural Research Service. pp 1/12. Coutinot D, Hoelmer K. 1999. Parasitoids of Lygus spp. in Europe and their potential for biological control

of Lygus spp. North America, in ANPP / Fifth International Conference on Pests in Agriculture,

Montpelier. pp 641/648. Craig CH. 1963. The alfalfa plant bug, Adelphocoris lineolatus (Goeze) in Northern Saskatchewan. Canadian

Entomologist 95:6/13. Craig CH, Loan CC. 1984. Lygus spp. plant bugs (Heteroptera: Miridae). In: Kelleher JS, Hulme MA,

editors. Biological control programmes against and weeds in Canada, 1969/1980. Norwich: Page

Bros. Ltd. pp 45/47. Cullen DW, Lees AK, Toth IK, Bell KS, Duncan JM. 2000. Detection and quantification of fungal and

bacterial potato pathogens in plants and soil. Bulletin OEPP 30:485/488. Day WH. 1999. Host preferences of introduced and native parasites (Hymenoptera: Braconidae) of phytophagous plant bugs (Hemiptera: Miridae) in alfalfa-grass fields in the northeastern USA.

BioControl 44:249/261. Day WH, Hedlund RC, Saunders LB, Coutinot D. 1990. Establishment of Peristenus digoneutis (Hymenoptera: Braconidae), a parasite of the tarnished plant bug (Hemiptera: Miridae), in the United

States. Environmental Entomology 19:1528/1533. Day WH, Tilmon KJ, Romig RF, Eaton AT, Murray KD. 2000. Recent range expansions of Peristenus digoneutis (Hymenoptera: Braconidae), a parasite of the tarnished plant bug (Hemiptera: Miridae), and high summer temperatures as a factor limiting its geographic distribution in North America. Journal of

the New York Entomological Society 108:326/331. Drea JJ, Dureseau L, Rivet E. 1973. Biology of Peristenus stygicus from Turkey, a potential natural enemy of

Lygus bugs in North America. Environmental Entomology 2:278/280. Erlandson M, Braun L, Baldwin D, Soroka J, Ashfaq M, Hegedus D. 2003. Molecular markers for Peristenus spp. (Hymenoptera: Braconidae) parasitoids associated with Lygus spp. (Hemiptera: Miridae). Canadian

Entomologist 35:71/83. Favia G, Cancrini G, Ricci I, Bazzocchi C, Magi M, Pietrobelli M, Genchi C, Bandi C. 2000. 5S ribosomal spacer sequences of some filarial parasites: comparative analysis and diagnostic applications. Molecular

and Cellular Probes 14:285/290. Gariepy TD, Levesque CA, de Jong SN, Rahe JE. 2003. Species-specific identification of the Neofabraea pathogen complex associated with pome fruits using PCR and multiplex DNA amplification. Mycological

Research 107:528/536. Hamelin RC, Berube P, Gignac M, Bourassa M. 1996. Identification of root rot fungi in nursery seedlings

by nested multiplex PCR. Applied Environmental Microbiology 62:4026/4031. Henson JM, French R. 1993. The polymerase chain reaction and plant disease diagnosis. Annual Review of

Phytopathology 31:81/109. Hormchan P. 1977. Biology of three exotic species, and role of native species of the genus Peristenus / Parasites of tarnished plant bug, Lygus lineolaris, in Mississippi. PhD thesis, Mississippi State University. 494 T. D. Gariepy et al.

Kengne P, Trung HD, Baimai V, Coosemans M, Manguin S. 2001. A multiplex PCR-based method derived from random amplified polymorphic DNA (RAPD) markers for the identification of species of the

Anopheles minimus group in Southeast Asia. Insect Molecular Biology 10:427/435. Kuhlmann U, Mason PG, Greathead D. 1998. Assessment of potential risks for introducing European Peristenus species as biological control agents of Lygus species in North America: A cooperative approach.

Biocontrol News and Information 19:83N/90N. Layton MB. 2000. Biology and damage of the tarnished plant bug, Lygus lineolaris, in cotton. Southwestern

Entomologist Supplement 23:7/20. Loan CC. 1965. Life cycle and development of Leiophron pallipes Curtis (Hymenoptera: Braconidae, Euphorinae) in five mirid hosts in the Belleville District. Proceedings of the Entomological Society of

Ontario 95:115/121. Loan CC. 1974. The European species of Leiophron Nees and Peristenus Foerster (Hymenoptera:

Braconidae, Euphorinae). Transactions of the Royal Entomological Society of London 126:207/238. Loan CC, Bilewicz-Pawinska T. 1973. Systematics and biology of four Polish species of Peristenus Foerster

(Hymenoptera: Braconidae: Euphorinae). Environmental Entomology 2:271/278. Lopez MM, Bertolini E, Olmos A, Caruso P, Gorris MT, Llop P, Panyalver R, Cambra M. 2003. Innovative

tools for detection of plant pathogenic viruses and bacteria. International Microbiology 6:233/243. Louda SM, Pemberton RW, Johnson MT, Follett PA. 2003. Non-target effects / the Achilles’ heel of biological control? Retrospective analysis to reduce risk associated with biocontrol introductions. Annual

Review of Entomology 48:365/396. May WE, Soroka JJ, Loeppky HA, Murrell DC. 2003. The effects of trichlorfon and deltamethrin on alfalfa plant bug and Lygus bug (Heteroptera: Miridae) populations in alfalfa grown in Canada. Crop Protection

22:883/889. Phuc HK, Ball AJ, Son L, Hanh NV, Tu ND, Lien NG, Verardi A, Townson H. 2003. Multiplex PCR assay for malaria vector Anopheles minimus and four related species in the Myzomyia series from Southeast Asia.

Medical and Veterinary Entomology 17:423/428. Ratcliffe ST, Robertson HM, Jones CJ, Bollero GA, Weinzierl RA. 2002. Assessment of parasitism of house fly and stable fly (Diptera: Muscidae) pupae by pteromalid (Hymenoptera: Pteromalidae) parasitoids

using a polymerase chain reaction assay. Journal of Medical Entomology 39:52/60. Rhainds M, English-Loeb G. 2003. Variation in abundance and feeding impact of tarnished plant bug (Hemiptera: Miridae) for different cultivars of strawberry: role of flowering phenology and yield

attributes. Journal of Economic Entomology 96:433/440. Rodrigues JLM, Silva-Stenico ME, Gomes JE, Lopes JRS, Tsai SM. 2003. Detection and diversity assessment of Xylella fastidiosa in field-collected plant and insect samples using 16S rRNA and gyrB

sequences. Applied Environmental Microbiology 69:4249/55. Rossiter BJF, Caskey CT. 1993. The polymerase chain reaction. In: Rapley R, Walker MR, editors. Molecular diagnostics. Oxford: Blackwell Scientific. p 76. Silva IMM, Honda J, van Kan F, Hu J, Neto L, Pintureau B, Stouthamer R. 1999. Molecular differentiation

of five Trichogramma species occuring in Portugal. Biological Control 16:177/184. Stouthamer R, Hu J, van Kan FJPM, Platner GR, Pinto JD. 1999. The utility of internally transcribed spacer 2 DNA sequences of the nuclear ribosomal gene for distinguishing sibling species of

Trichogramma. BioControl 43:421/440. Strong FE. 1970. Physiological injury caused by Lygus hesperus. Journal of Economic Entomology 63:808/ 814. Taylor E, Bates J, Kenyon D, Maccaferri M, Thomas J. 2001. Modern molecular methods for

characterisation and diagnosis of seed-borne fungal pathogens. Journal of Plant Pathology 83:75/81. Tilmon KJ, Hoffmann MP. 2003. Biological control of Lygus lineolaris by Peristenus spp. in strawberry.

Biological Control 26:287/292. Tilmon KJ, Danforth BN, Day WH, Hoffmann MP. 2000. Determining parasitoid species composition in a

host population: A molecular approach. Annals of the Entomological Society of America 93:640/647. Valles SM, Oi DH, Perera OP, Williams DF. 2002. Detection of Thelohanis solenopsae (Microsporidia: Thelohaniidae) in Solenopsis invicta (Hymenoptera: Formicidae) by multiplex PCR. Journal of

Invertebrate Pathology 81:196/201. van Lenteren JC, Babendreier D, Bigler F, Burgio G, Hokkanen HMT, Kuske S, Loomans AJM, Menzler- Hokkanen I, van Rijn PCJ, Thomas MB, Tommasini MG, Zeng Q-Q. 2003. Environmental risk

assessment of exotic natural enemies used in inundative biological control. BioControl 48:3/38. Single-step multiplex PCR for detection of European Peristenus spp. 495

Wilson WC, Chase CCL. 1993. Nested and multiplex polymerase chain reactions for the identification of bluetongue virus infection in the biting midge, Culicoides variipennis. Journal of Virological Methods 45:39/47. Wise IL, Lamb RJ. 2000. Seasonal occurrence of plant bugs (Hemiptera: Miridae) on oilseed flax (Linaceae)

and their effect on yield. Canadian Entomologist 132:369/371. Wise IL, Tucker JR, Lamb RJ. 2000. Damage to wheat seeds caused by a plant bug, Lygus lineolaris L. Canadian Journal of Plant Science 80:459/461. Young OP. 1986. Host plants of the tarnished plant bug, Lygus lineolaris (Heteroptera: Miridae). Annals of

the Entomological Society of America 79:747/762. Zarlenga DS, Higgins J. 2001. PCR as a diagnostic and quantitative technique in veterinary parasitology.

Veterinary Parasitology 101:215/230. Zhu Y, Greenstone MH. 1999. Polymerase chain reaction techniques for distinguishing three species and two strains of Aphelinus (Hymenoptera: Aphelinidae) from Diuraphis noxia and Schizaphis graminum (Homoptera: Aphididae). Annals of the Entomological Society of America 92:71/79. Zhu YC, Burd JD, Elliott NC, Greenstone MH. 2000. Specific ribosomal DNA marker for early polymerase chain reaction detection of Aphelinus hordei (Hymenoptera: Aphelinidae) and Aphidius colemani (Hymenoptera: Aphidiidae) from Diuraphis noxia (Homoptera: Aphididae). Annals of the Entomological Society of America 93:486/491. Zhu YC, Riddick EW, Williams L, Schotzko DJ, Logarzo GA, Jackson CG. 2004. Potential of detection and identification of nymphal parasitoids (Hymenoptera: Braconidae) of Lygus bugs (Heteroptera: Miridae) by using polymerase chain reaction and ITS2 sequence analysis techniques. Annals of the Entomological Society of America 97:743/752.