Genetic differentiation among geographic populations of Gonatocerus ashmeadi, the predominant egg parasitoid of the glassy-winged sharpshooter, Homalodisca coagulata

Jesse H. de León and Walker A. Jones

Introduction
Materials and Methods
Results and Discussion
Acknowledgements
References

The aim of genetically comparing different populations of the same species of natural enemies is to identify the strain that is most adapted to the environment where it will be released. In the present study, Inter-Simple Sequence Repeat-Polymerase Chain Reaction (ISSR–PCR) was utilized to estimate the population genetic structure of Gonatocerus ashmeadi (Girault) (Hymenoptera: Mymaridae), the predominant egg parasitoid of Homalodisca coagulata (Say) (Homoptera: Cicadellidae), the glassy-winged sharpshooter. Six populations from throughout the U.S. and a population from Argentina identified as near G. ashmeadi were analyzed. Four populations (California; San Antonio, Texas; Weslaco, Texas [WTX-2]; and Florida) were field collected and two (Louisiana and Weslaco, Texas [WTX-1]) were reared. Three ISSR–PCR reactions were pooled to generate 41 polymorphic markers among the six U.S. populations. Nei's expected heterozygosity values (h), including the reared population from Louisiana, were high (9.01–14.3%) for all populations, except for a reared population from WTX-1 (2.9%). The total genetic diversity value (Ht) for the field populations was high (23%). Interestingly, the Florida population that was collected from one egg mass (siblings) generated the greatest number of polymorphic markers (20) and was observed with the highest gene diversity value (14.3%). All populations, except WTX-2 generated population-specific markers. Comparison of genetic differentiation estimates, which evaluate the degree of genetic subdivision, demonstrated good agreement between G_ST and values, 0.38 and 0.50, respectively for field populations, and 0.44 and 0.50, respectively for all populations. Genetic divergence (D) indicated that the WTX-1 population was the most differentiated. Average D results from the Argentina population support the taxonomic data that it is a different species. The present results estimate the population genetic structure of G. ashmeadi, demonstrating genetic divergence and restricted gene flow (Nm = 0.83) among populations. These results are of interest to the Pierce's disease/glassy-winged sharpshooter biological control program because the key to successful biological control may not be in another species, but instead in different geographic races or biotypes.

Abstract
Materials and Methods
Results and Discussion
Acknowledgements
References

Gonatocerus ashmeadi (Girault) (Hymenoptera: Mymaridae) is a primary egg parasitoid of Homalodisca coagulata (Say) (Homoptera: Cicadellidae), the glassy-winged sharpshooter (Huber 1998). A survey conducted by Triapitsyn et al. (1998) in California, Florida, and Louisiana concluded that this mymarid wasp was the predominant or the most common natural enemy of H. coagulata. A biological control program is currently in progress in California against H. coagulata because this xylem feeding leafhopper is a serious economic pest that vectors a strain of Xylella fastidiosa, a bacterium that causes Pierce's disease in grapevines (Vitis vinifera and V. labrusca) (Hopkins and Mollenbauer 1973). H. coagulata are native to the southern United States, ranging from Florida to Texas and northern Mexico (Young 1958; Turner and Pollard 1959; Nielsen 1968; Brlansky et al. 1983). Within the last 10 years, H. coagulata have established in southern California where they pose a serious threat to the wine and table grape industry (Sorensen and Gill 1996).

Studies of allele or marker frequencies in naturally occurring parasitoid populations are important, not only for identifying genetic variation of potential benefit in the selection and screening of biological control organisms, but also for the detection of genetic markers indicative of specific biological traits or geographic origins. In addition, the recognition of intraspecific variation can be as crucial for the success of biological control programs as is sound species determination. Populations of parasitoids from distinct geographical regions may differ in relevant biological characteristics of importance to biological control (Powell and Walton 1989; Narang et al. 1993b; Unruh and Woolley 1999). An aim of genetically comparing different populations of the same species of natural enemies is to identify the strain that is most adapted to the environment where it will be released (Messenger and van den Bosch 1971); in other words, the key to successful biological control may not be in another species, but instead in different geographic races or biotypes of one species (Diehl and Bush 1984). Laboratory populations of parasitoids are maintained for three main purposes; to provide material for research into their basic biology and behavior, to meet quarantine requirements, and to mass-rear insects for release in biological control programs (Unruh et al. 1983; Powell and Walton 1989; Hopper et al. 1993). Reliable methods are needed for distinguishing various exotic strains of these biological control agents from those indigenous to the U.S., including parasitoids from different states within the U.S. Release of unidentified and uncharacterized strains can make it difficult to document their establishment and dispersal. Therefore, genetic typing of strains prior to their release in the field is highly desirable (Narang et al. 1993b). Molecular methods are also needed to monitor the quality of laboratory cultures, to detect development of major deviations in genetic composition from that existing in field populations, and for the detection of cross-contamination between cultures of sibling species (Powell and Walton 1989; Menken and Ulenberg 1987; Narang et al. 1993b; Hopper et al. 1993; Unruh and Woolley 1999).

A sensitive approach for obtaining polymorphic DNA markers is based on the use of simple sequence repeats (SSR) or microsatellites. Microsatellites are ubiquitous in eukaryotic genomes and can be found in both protein-coding and noncoding regions (Tóth et al. 2000). Microsatellites are widely used as genetic markers because they are co-dominant, multiallelic, easily scored, and highly polymorphic, although drawbacks for their use are the time and cost required to characterize them (Karp and Edwards 1997). However, a DNA fingerprinting procedure, Inter-Simple Sequence Repeat-Polymerase Chain Reaction (ISSR–PCR) (Zietkiewicz et al. 1994), including adding a primer pair to the reaction (pp-ISSR–PCR) (Prevost and Wilkinson 1999; Cekic et al. 2001), permits detection of DNA variation in microsatellites without the need to isolate and sequence specific DNA fragments. The approach for this technique involves amplification with oligonucleotide primers corresponding directly to random SSR sites. This involves the use of 5'-anchored or compound ISSR primers where the anchor serves to fix the annealing of the primer to a single position of the target site, thus resulting in a low level of slippage during amplification. Many classes of microsatellite repeat motifs have been identified, though the class most abundant in eukaryotic genomes is the CA-repeat. The presence of these repeat motifs in high copy number and their dispersion throughout the genome of all eukaryotes tested has been demonstrated by earlier studies (Hamada et al. 1982; Tautz and Renz 1984; Tóth et al. 2000). Therefore, because of their high density, a few oligonucleotides complementary to these CA-repeat motifs can be used to target a significant portion of the genome and reveal highly polymorphic banding patterns (Zietkiewicz et al. 1994).

To gain an understanding of the population genetics of G. ashmeadi, molecular genetic markers were developed by the ISSR–PCR DNA fingerprinting procedure. Recently, we developed DNA markers for H. coagulata for the purpose of estimating genetic variation in natural populations (de León and Jones 2004). We screened and compared four DNA fingerprinting procedures and determined that procedures incorporating SSR or microsatellites were the most sensitive with H. coagulata template. To this end, the specific objectives of the present study were to 1) estimate genetic variation or gene diversity within and among populations, 2) estimate the population genetic structure, 3) determine whether ISSR–PCR was sensitive enough to identify diagnostic markers in geographic populations, and 4) confirm the species identification of a population of egg parasitoids from Argentina identified as near G. ashmeadi. We demonstrated the ability of the ISSR–PCR procedure with compound and 5'-anchored ISSR primers, including combining them, to generate polymorphic markers and estimate geographic variation in G. ashmeadi populations.

Abstract
Introduction
Results and Discussion
Acknowledgements
References

Insect collection

G. ashmeadi were collected from H. coagulata egg masses from various regions of the United States by the individuals listed on Table 1. A population of wasps from Argentina classified as near G. ashmeadi (M02012; 30 males) was also included. S. Triapitsyn from the University of California Riverside identified this egg parasitoid. This egg parasitoid was imported to the USDA, APHIS, Arthropod Quaratine Facility, Edinburg, Texas. M02012 was the unique code associated with this imported colony. Specimens were then shipped to Weslaco in 95% non-denatured ethanol. G. ashmeadi from California (CA; ~10 egg masses; 17 females and 13 males), Weslaco, Texas (WTX-2; ~4 egg masses; 22 females and 9 males), San Antonio, Texas (SATX; ~10 egg masses; 23 females and 7 males), and Quincy, Florida (QFL; 1 egg mass; 10 females and 3 males) were field collected. Populations from Weslaco, Texas (WTX-1; F1/F2 generation; 23 females and 7 males) and Louisiana (LA; F5 generation; 11 females and 19 males) were reared. The Weslaco populations (WTX-1 and -2) were collected a year apart, with WTX-1 collected in 2002 and WTX-2 collected in 2003.

Genomic DNA isolation

Genomic DNA was extracted according to standard methods (Sambrook and Russell 2001). Individual wasps were homogenized on ice in 1.5 ml microfuge tubes in 60 µl of lysis buffer [10 mM Tris-HCl (pH 7.5), 1 mM EDTA (pH 8.0), 1% IGEPAL CA-630)] with two 20 s bursts with 10 min intervals on ice (Pellet Pestle Motor, Bel-Art Products, www.bel-art.com). To avoid cross-contamination between samples, a sterile plastic pestle was used per insect. The final DNA pellet was resuspended in 61 µl of TE [Tris-HCl (pH 7.5), 1 mM EDTA (pH 7.5)]. To confirm for the present of genomic DNA, amplification reactions were performed with 1 µl of stock DNA and 28S primers at a Tm of 65 °C (forward CCCTGTTGAGCTTGACTCTAGTCTGGC and reverse AAGAGCCGACATCGAAGGATC) with 1.5 mM MgCl₂ with the assay conditions described below.

ISSR–PCR DNA fingerprinting

ISSR–PCR amplification reactions and modification of ISSR primers are described in de León and Jones (2004). Three reactions for this study included: (a) ISSR–PCR p-6, (b) pp-ISSR–PCR #4 and (c) #9 (Table 1). The reactions were performed with 5'-anchored or compound ISSR primers that included either a single primer or a primer pair: (a) HVH(TG)₇T, (b) KKVRVRV(TG)₆ and CCAG(GT)₇, and (c) CCAG(GT)₇ and T(GT)₇(AT)₂. The use of these primers in their original forms are reported in several sources (Wu et al. 1994; Zietkiewicz et al. 1994; Fisher et al. 1996; Witsenboer et al. 1997). Notations: K = G/T, V = G/C/A, R = G/A, H = A/T/C, and V = G/C/A. The reactions were performed in a final volume of 20 µl with the following components: 1X PCR buffer (50 mM KCl, 20 mM Tris-HCl [pH 8.4], 2.0 mM MgCl₂, and 0.01% gelatin), 0.25 mM deoxynucleotide triphosphates, 0.25 µM ISSR primer(s), 1.0 µl of stock genomic DNA and 0.05 U/µl Taq DNA Polymerase (New England Biolabs, Beverly, Massachusetts). The cycling parameters were as follows: 1 cycle at 94 °C for 2 min followed by 45 cycles at 94 °C for 1 min, Tm (Table 1) for 1 min, and 72 °C for 2 min. Reactions were optimized for amount of genomic DNA, annealing temperature, MgCl2 concentration, and cycle number. Negative control reactions were performed in the absence of genomic DNA. Amplification products were loaded onto 2% agarose gels and submitted to electrophoresis in 1X TBE buffer (90 mM Tris-borate, 2 mM EDTA) in the presence of 0.2 µg/ml ethidium bromide. Gels were photographed with the Chemi Doc System and markers/bands were scored via the Quantity One Software™ (Bio-Rad Laboratories, www.bio-rad.com).

Data analysis

Data were analyzed with the pooled products of the three amplification reactions. Bands observed in each lane were compared with all the other lanes of the same gel and only reproducible bands were scored as present (1) or absent (0). Each ISSR–PCR reaction generated unique banding patterns. Only bands unique to each reaction were scored and in cases where a similar sized band was generated by more than one reaction, a band from only one reaction was chosen. ISSR–PCR markers are treated as dominant markers (reviewed in Wolfe and Liston 1998). Fragment sizes were estimated based on the 1.0 Kb Plus DNA Ladder (Invitrogen Life Technologies, www.invitrogen.com) according to the algorithm provided in the Quantity One software. Genetic variation was analyzed with the programs POPGENE 3.2 (Yen et al. 1996) and Tools for Population Genetic Analyses (TFPGA, ver. 1.3, Miller 1997) as previously performed (de León et al. 2004). Marker frequencies were estimated based on Lynch and Milligan's (1994) Taylor expansion estimate. Different approaches to estimate population differentiation- exact tests (Raymond and Rousset 1995), G_ST (Nei 1987), (Weir 1990, 1996), and dendrograms based on genetic distance (Nei 1987; Reynolds et al. 1983) by the UPGMA method (Sneath and Sokal 1973) were applied and compared.

Abstract
Introduction
Materials and Methods
Acknowledgements
References

ISSR–PCR marker heterozygosity

A total of 41 polymorphic markers were generated in the six populations of G. ashmeadi (163 individuals) from the U.S. with the ISSR–PCR reactions (Table 2), 34 or 83% of these markers tested neutral based on the Ewens-Watterson analysis (Manly 1985). Reactions a, b, and c generated several polymorphic markers each. Results of G₂-contingency tests indicated significant heterogeneity of marker frequency across all U.S. populations for 31 of 41 markers and for 25 of 34 markers for the field populations. Also shown on Table 2 are the number of population-specific or rare markers identified in each geographic population, the ISSR–PCR reaction that generated the markers, the estimated size of the markers (bp), and the estimated frequency of each marker. All populations, except WTX-2, were associated with population-specific markers. Reared populations, WTX-1 (F1/F2) and LA (F5) produced two and five population-specific markers, respectively. Rare markers can be an indication of population subdivision. It was surprising that WTX-2 was not found associated with any population-specific markers since it was field collected from about ten different egg masses. Absence of rare or population-specific markers can indicate a genetic bottleneck in recent history of a population. This occurs when a population has undergone a drastic contraction, and its numbers rebuilt from a small number of individuals, a common occurrence in insect colonies (Munstermann 1994); however none of the populations were colonized in our study. Only seven (17%) of the 41 markers were shared among all G. ashmeadi populations. None of the same fixed markers were identified across all G. ashmeadi populations, though each population was associated with fixed or monomorphic markers. The following markers were present at high frequency: #10, 25, and 36 and 37 for CA, WTX-1, and LA, respectively. Together, markers # 36, 37, and 48 can be considered diagnostic for the LA population. Seven fixed markers (#41–47) were identified in ARG that can be considered diagnostic.

Genetic diversity

Within populations, gene diversity values (h) were observed ranging from 2.9 to 14.3% with WTX-1 (F1/F2) having the lowest and QFL having the highest value (Table 3). In general, the two Weslaco populations (WTX-1 and -2) were found to have the lowest h values. No significant differences in h were seen between the two Weslaco populations (t = 1.49, df = 58, P > 0.05), but significant differences (P < 0.05) were observed between WTX-1 and the rest of the U.S. populations. Interestingly, no significant differences in h were observed between the reared LA (F5) and the rest of the field populations. The fact that QFL was associated with an h value of 14.3% was surprising since this population was from a single egg mass. Overall, the field populations and all the U.S. G. ashmeadi populations together had an h value of 23 and 20.8%, respectively. The number of polymorphic markers ranged from 12 to 20 with WTX-1 and -2 having the lowest and QFL the highest. Percentage of polymorphic markers (%P) ranged from 29.3 to 58.8%, but overall, 100% of the ISSR–PCR markers were polymorphic, including the field populations analyzed separately. The two Weslaco populations were associated with the lowest %P and QFL with the highest. It is interesting to note that even though both LA and WTX-1 were reared, WTX-1 is presented with a significantly (P < 0.05) lower h value. On the other hand, the Weslaco populations contain the same number of monomorphic markers (7 vs. 6), polymorphic markers (12 vs. 13), total number of markers (19 each), and %P (29.3 vs. 31.7%), but a three-fold difference in h (2.9 vs. 9.0%) (non-significant) is seen. If a decrease in total number of markers occurs, a genuine loss of genetic material is indicated (Munstermann 1994), but this situation is not seen between the two Weslaco populations. These results may indicate a real genetic difference between the two populations, including the possibility of sympatric strains. The highest polymorphic ratio (1.54) was seen with QFL, this value was about three- and four-fold higher than CA, SATX, LA and WTX-1 and -2, respectively. This ratio also demonstrates the high diversity in QFL; this is an interesting result since all individuals are assumed to be siblings (from the same egg mass).

ISSR–PCR differentiation among U.S. G. ashmeadi populations

Tables 4 and 5 present the results from the different approaches used to apportion variation into within- and among-populations levels. Simultaneous exact tests for population differentiation (Table 4) indicate that highly significant differences in marker frequencies exist among the six populations; significant differences are also seen among the field populations analyzed separately. In addition, pairwise exact tests showed highly significant differences in marker frequencies between populations, including between WTX-1 and -2. These statistically significant tests suggest that discrete subpopulations exist. Total genetic diversity (Ht) values, 23% for field and 20.9% for all populations (Table 5) are in agreement with the h values on Table 3. The average genetic diversity within populations (Hs) value for the field populations is 14.4%. Table 5 also shows a comparison of other genetic differentiation estimates, G_ST and , which evaluate the degree of genetic subdivision among populations. Good agreement was seen between G_ST and values, respectively for field and for all populations. The G_ST values, for field and all populations indicate that about 38 and 44% of the variance is distributed among populations, and 62 and 56% is distributed within populations, respectively. The values show that about 50% of the variance is seen among populations in both field and all populations. The indirect estimate of gene flow, Nm base on G_ST, demonstrated low values for both field and all U.S. populations. These values indicate restricted gene flow among the populations. Overall, genetic differentiation measurements (G2 tests, exact tests, G_ST, , and Nm) indicate genetic divergence among G. ashmeadi populations from the U.S.

Genetic relatedness among G. ashmeadi populations from the U.S.

Since the levels of genetic divergence among populations can also be summarized by calculating pairwise estimates for genetic distance, we used the procedures of Nei (1978) and Reynolds et al. (1983) (Table 6). Average genetic divergence (D) among both field [Nei = 0.1702 (0.1021–0.2230); Reynolds = 0.6208 (0.4069–0.8138)] and all populations [Nei = 0.1304 (0.0715–0.2024); Reynolds = 0.6512 (0.3705–0.8890)] was high. We compared the level of genetic divergence between the field populations and the WTX-1 and LA reared populations and found mean D values of 0.1806 (Nei) and 0.8589 (Reynolds) and 0.1065 (Nei) and 0.5371 (Reynolds), respectively. These results indicate that WTX-1 is more diverged than LA. A comparison of Nei's genetic distance within the Texas populations, WTX-2 vs. WTX-1 (0.1391) and WTX-2 vs. SATX (0.1286), showed that divergence is slightly higher between the Weslaco populations. Sympatric species tend to have higher levels of genetic differentiation; more work is needed to confirm this possibility. The divergence between ARG (near G. ashmeadi) and all U.S. G. ashmeadi populations was very high, 0.3633 (Nei) and 1.6093 (Reynolds), respectively. These results support the taxonomic data that ARG (near G. ashmeadi) is another species. Dendrograms based on Nei's genetic distance are shown on Fig. 1 with all populations including ARG (Fig. 1A) and the field populations analyzed separately (Fig. 1B). At least two main clusters are identified on the dendrogram with ARG clustered as an outlier (Fig. 1A). Within a second cluster or all G. ashmeadi from the U.S., WTX-1 appears to be the most differentiated (Fig. 1A). The CA population appears to form a second subcluster and the two southeastern populations, LA and QFL form a single cluster showing their genetic similarity. The WTX-1 and -2 populations are distributed in different clusters. Also within Texas (Fig. 1B), WTX-2 and SATX show divergence as they appear on a separate cluster. It is interesting to note that this same pattern of differentiation is seen with H. coagulata within Texas (de León et al. 2004).

In general, genetic variation of Hymenoptera as measured by isozyme electrophoresis is approximately one-third of that found in diploid insects (Berkelhamer 1983; Graur 1985; Unruh et al. 1986; Powell and Walton 1989; Owen 1993; Unruh and Woolley 1999), though there have been limited reports of high genetic diversity (Sheppard and Heydon 1986; Boato and Battisti 1996). Berkelhamer (1983) demonstrated that Hymenoptera were associated with low gene diversity values 0.037 (3.7%) vs. 51 diploid insects 0.135 (13.5%). On the other hand, Sheppard and Heydon (1986) demonstrated mean expected gene diversity values of 0.121 to 0.144 (12.1 to 14.4%) in sawflies and Boato and Battisti (1996) demonstrated a mean expected gene diversity value of 0.197 (19.7%) in Cephalcia populations. Though a direct comparison of mean expected gene diversity values can not be made between enzyme electrophoreses studies and PCR methods incorporating SSR, in the present study with G. ashmeadi we observed individual mean expected gene diversity values ranging from 2.9 to 14.3%, but a total genetic diversity value (Ht) of 23.1% for the field populations. This Ht value falls within the range described by Boato and Battisti (1996) for Cephalcia populations. PCR-based techniques, such as, RAPDs have been described to have generated many polymorphisms in four parasitic hymenoptera (Narang et al. 1993a) and in 40 Trichogramma isofemale lines that revealed several RAPD polymporphic markers (Vanlerberghe-Masutti 1994); but these studies did not report expected mean gene diversity values and therefore, direct comparisons are difficult to make. In a social wasp, Polistes annularis, the SSR technique identified polymorphic microsatellite loci that produced a mean observed heterozygous value of 0.62 (Hughes and Queller 1993); the SSR and ISSR–PCR methods are different techniques though are similar in that they both target SSR.

The ISSR primers that were utilized in this study contained CA-repeat motifs in their sequence (see Materials and Methods). Recently, in a genetic study of H. coagulata, we demonstrated high numbers of polymorphic markers with various ISSR primers containing CA-repeat motifs in their sequences (de León and Jones 2004). Based on a survey of SSR or microsatellites in different genomes by Tóth et al. (2000) and earlier studies (Hamada et al. 1982; Tautz and Renz 1984), the most frequent repeat motif in Arthropoda (mainly Drosophila) in all regions of the genome, which include intergenic regions, introns, and exons is the dinucleotide repeat AC. This dinucleotide repeat motif is found in greater frequency in intergenic regions and introns than in exons. The number of polymorphic markers produced in the present study therefore confirmed the utility of using ISSR primers containing these repeat motifs. The ISSR–PCR method has been utilized successfully for genetic characterization of the silkworm moth, Bombyx mori (Reddy et al. 1999), for genetic variability studies of Biomphalaria straminea snail complexes (Caldeira et al. 2001), and in many plant studies (reviewed by Wolfe and Liston 1998). pp-ISSR–PCR has been used successfully to fingerprint potato cultivars (Prevost and Wilkinson 1999) and for genetic linkage analysis in the seasonal flowering locus in Fragaria (Cekic et al. 2001).

To our knowledge, the present report represents the first study of the population genetics of the Hymenopteran G. ashmeadi performed with the ISSR–PCR DNA fingerprinting method, including adding a second ISSR primer to the reaction (pp-ISSR–PCR). The ISSR–PCR procedure is a sensitive method and requires no prior knowledge of DNA sequence information (Zietkiewicz et al. 1994; reviewed in Karp and Edwards 1997). Amplification of G. ashmeadi template with the ISSR–PCR procedure (Table 2) produced a high quantity of polymorphic markers. Within the six U.S. populations of G. ashmeadi, 41 polymorphic markers were generated with the three pooled reactions. The high number of polymorphic markers produced is a good indication that a significant portion of the genome was targeted.

In summary, the major observations of this study were that 1) among G. ashmeadi populations, based on genetic differentiation measurements (exact test, G_ST, ), pronounced genetic structure was identified; 2) the mean expected gene diversity value for LA (F5) did not differ from field populations, whereas WTX-1 was observed with a significantly lower mean expected gene diversity value as compared to field populations (except WTX-2); 3) QFL generated the most polymorphic markers (20) with only 13 individuals, even though they were all siblings or from one egg mass. This is an interesting result, since it may be assumed that siblings are not associated with high variability or have isofemale line characteristics. These results indicate that G. ashmeadi parasitoid siblings somehow manage to maintain their genetic diversity. Further studies are required to confirm this observation in this species and other Gonatocerus species. Variation within 10 male individuals (Anaphes sp. nov.) was demonstrated with RAPD markers by Landry et al. (1993), but they were not from the same egg mass; 4) based on genetic distance or average divergence, WTX-1 appeared to be the most differentiated population. Within Texas, field populations WTX-2 and SATX appeared on separate clusters, indicating that these populations are differentiated even though they are within the same state; and 5) The ARG population (near G. ashmeadi, M02012) is confirmed to be a different species. More research is required to confirm these results, sequencing of standard genes (e. g., mitochrondrial cytochrome oxidase [COI and II]) and ITS fragments are in progress.

Abstract
Introduction
Materials and Methods
Results and Discussion
References

We thank Marissa González and Lisa A. Ledezma for their excellent technical assistance. We are thankful to our following colleagues for samples of G. ashmeadi, David J. W. Morgan from Agricultural Operations, University of California, Riverside, California, to Russell F. Mizell III from the University of Florida, Quincy, Florida, to William Warfield from USDA, ARS Weslaco, Texas, to Isabelle Lauziere from Mission Plant Protection Center, Moore Air Field, Edinburg, Texas, and to Eduardo Virla and Guillermo Logarzo from USDA, ARS, SABCL. We thank Kyung S. Kim from USDA, ARS, CICGRU, Iowa State University, Ames, Iowa and Thomas R. Unruh from the USDA, ARS in Wapato, Washington for reading a previous version of this manuscript. Lastly, we thank the anonymous reviewers for improving the manuscript.

Disclaimer

Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. Department of Agriculture.

Abstract
Introduction
Materials and Methods
Results and Discussion
Acknowledgements

Berkelhamer RC. 1983. Intraspecificgenetic variation and haplodiploidy, eusociality, and polygyny in the Hymenoptera. Evolution 37: 540–545.

Boato A, Battisti A. 1996. High genetic variability despite haplodipoidy in primitive sawflies of the genus Celpalcia (Hymenoptera: Pamphiliidae). Experientia 52: 516–521.

Brlansky RH, Timmer LW, French WJ, McCoy RE. 1983. Colonization of the sharpshooter vectors, Oncometopia nigricans and Homalodisca coagulata, by xylem-limited bacteria carriers of xylem-limited bacterial disease. Phytopathology 73: 530–535.

Caldeira RL, Vidigal THDA, Simpson AJG, Carvalho OS. 2001. Genetic variability in Brazilian populations of Biomphalaria straminea complex detected by Simple Sequence Repeat Anchored Polymerase Chain Reaction Amplification. Memorias do Instituto Oswaldo Cruz 96: 535–544.

Cekic C, Battey NH, Wilkinson MJ. 2001. The potential of ISSR–PCR primer-pair combinations for genetic linkage analysis using the seasonal flowering locus in Fragaria as a model. Theoretical and Applied Genetics 103: 540–546.

de León JH, Jones WA. 2004. Detection of DNA Polymorphisms in Homalodisca coagulata (Homoptera: Cicadellidae) by Polymerase Chain Reaction-based DNA Fingerprinting Methods. Annals of the Entomological Society of America 97: 574–585.

de León JH, Jones WA, Morgan DJW. 2004. Population genetic structure of Homalodisca coagulata (Homoptera: Cicadellidae), the vector of the bacterium Xylella fastidiosa causing Pierce's disease in grapevines. Annals of the Entomological Society of America 97: 809–818.

Diehl, SR, Bush GL. 1984. An evolutionary and applied perspective of insect biotypes. Annual Review of Entomology 29: 471–504.

Fisher PJ, Gardner RC, Richardson TE. 1996. Single locus microsatellites isolated using 5' anchored PCR. Nucleic Acids Research 24: 4369–4371.

Graur D. 1985. Gene diversity in Hymenoptera. Evolution 39: 190–199.

Hamada H, Petrino MG, Kakunaga T. 1982. A novel repeated element with Z-DNA forming potential is widely found in evolutionary diverse eukaryotic genomes. Proceedings of the National Academy of Science USA 79: 6465–6469.

Hopkins DL, Mollenhauer HH. 1973. Rickettsia-like bacterium associated with Pierce's disease of grapes. Science 179: 298–300.

Hopper KR, Roush RT, Powell W. 1993. Management of genetics of biological control introductions. Annual Reviews in Entomology 38: 27–51.

Huber JT. 1998. The species groups of Gonatocerus Nees in North America with a revision of the sulphuripes an ater groups (Hymenoptera: Mymaridae). Memoirs of the Entomological Society of Canada 141: 1–109.

Hughes C, Queller DC. 1993. Detection of highly polymorphic microsatellite loci in a species with little allozyme polymorphism. Molecular Ecology 2: 131–137.

Karp A, Edwards J. 1997. DNA markers: A global overview. In: Caetano-Anolles G, Gresshoff PM, editors. DNA markers-protocols, applications, and overviews, 1-13. New York: Wiley-Liss, Inc.

Landry BS, Dextraze L, Biovin G. 1993. Random amplified polymorphic DNA markers for DNA fingerprinting and genetic variability assessment of minute parasitic wasp species (Hymenoptera: Mymaridae and Trichogrammatidae) used in biological control programs of phytophagous insects. Genome 36: 580–587.

Lynch M, Milligan BG. 1994. Analysis of population genetic structure with RAPD markers. Molecular Ecology 3: 91–99.

Manly BFJ. 1985. The statistics of natural selection on animal populations / Bryan F.J. Manly. London; New York: Chapman and Hall, 1987, c1985, p. 272–282.

Menken SBJ, Ulenberg SA. 1987. Biochemical characters in agricultural entomology. Agricultural Zoology Reviews. 2: 305–353.

Messenger PS, van den Bosch R. 1971. The adaptability of introduced biological control agents. In: Huffaker, CB editor. Biological Control, 68–92. New York: Plenum Press.

Miller MP. 1997. Tools for population genetic analysis (TFPGA) 1.3: A windows program for the analysis of allozyme and molecular population genetic data. Computer software distributed by author.

Munstermann LE. 1994. Unexpected Genetic Consequences of Colonization and Inbreed: Allozyme Tracking in Culicidae (Diptera). Annals of the Entomological Society of America 87: 157–164.

Narang SK, Leopold RA, Krueger CM, DeVault JD. 1993a. Dichomotomus RAPD-PCR Key for Identification of Four Species of Parasitic Hymenoptera. In: Narang SK, Barlett AC, Faust RM, eds. Applications of Genetics to Arthropods of Biological Control Significance, 53–67. CRC Press Inc, Boca Raton, Florida.

Narang SK, Tabachnick WJ, Faust RM. 1993b. Complexities of population genetic structure and implications for biological control programs. In: Narang SK, Barlett AC, Faust RM, editors. Applications of Genetics to Arthropods of Biological Control Significance, 19–52. CRC Press Inc, Boca Raton, Florida.

Nei M. 1973. Analysis of gene diversity in subdivided populations. Proceeding of the National Academy of Science U.S.A. 70: 3321–3323.

Nei M. 1977. F-statistics and analysis of gene diversity in subdivided populations. Annals of Human Genetics 41: 225–233.

Nei M. 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 89: 583–590.

Nei M. 1987. Molecular evolutionary genetics / Masatoshi Nei. New York : Columbia.

Owen RE. 1993. Genetics of Parasitic Hymenoptera. In: Narang SK, Barlett AC, Faust RM, editors. Applications of Genetics to Arthropods of Biological Control Significance, 69–89. Boca Raton, Florida: CRC Press Inc.

Nielson MW. 1968. The leafhopper vectors of phytopathogenic viruses (Homoptera, Cicadellidae) taxonomy, biology, and virus transmission. USDA Technical Bulletin. 1382: 81–84.

Powell W, Walton MP. 1989. The use of electrophoresis in the study of hymenopteran parasitoids of agricultural pest. In: Loxdale HD, den Hollander J, editors. Electrophoretic Studies on Agricultural Pests, 443–465. Oxford: Clarendon.

Prevost A, Wilkinson MJ. 1999. A new system of comparing PCR primers applied to ISSR fingerprinting of potato cultivars. Theoretical and Applied Genetics 98: 107–112.

Raymond ML, Rousset F. 1995. An exact test for population differentiation. Evolution 49: 1280–1283.

Reddy KD, Nagaraju J, Abraham EG. 1999. Genetic characterization of the silkworm Bombyx mori by simple sequence repeat (SSR)-anchored PCR. Heredity 83: 681–687.

Reynolds J, Weir BS, Cockerham CC. 1983. Estimation of the coancestry coeffient: basis for a short-term genetic distance. Genetics 105: 767–779.

Sambrook J, Russell DW. 2001. Molecular Cloning: A Laboratory Manual, 3rd ed. Cold Spring Harbor, New York: Cold Spring Harbor Laboratory Press.

Sheppard WS, Heydon SL. 1986. High levels of genetic variability in three male-haploid species (Hymenoptera: Argidae, Tenthredinidae). Evolution 40: 1350–1353.

Slatkin M, Barton NH. 1989. A comparison of three indirect methods for estimating average levels of gene flow. Evolution 43: 1349–1368.

Sneath PHA, Sokal RR. 1973. Numerical Taxonomy. San Francisco: W. H. Freeman. 573 p.

Sorensen JT, Gill RJ. 1996. A range extension of Homalodisca coagulata (Say) (Hemiptera: Clypeorrhyncha: Cicadellidae) to southern California. Pan-Pacific Entomology 72: 160–161.

Tautz D, Renz M. 1984. Simple sequences are ubiquitous repetitive components of eukaryotic genomes. Nucleic Acids Research 12: 4127–4138.

Tóth G, Gáspári Z, Jurka J. 2000. Microsatellites in different eukaryotic genomes: Survey and Analysis. Genome Research 10: 967–981.

Triapitsyn SV, Mizell RF, Bosart JL, Carlton CE. 1998. Egg parasitoids of Homalodisca coagulata (Homoptera: Cicadellidae). Florida Entomologist 81: 241–243.

Turner WF, Pollard HN. 1959. Life histories and behavior of five insect vectors of phony peach disease. USDA Technical Bulletin. 1188: 28 pp.

Unruh TR, White W, Gonzalez D, Gordh G, Luck RF. 1983. Heterozygosity and effective size in laboratory populations of Aphidius ervi (Hymn.: Aphidiidae). Entomophaga 28: 245–58.

Unruh TR, White W, Gonzalez D, Luck RF. 1986. Electrophoretic studies of parasitic hymenoptera and implications for biological control. In: Patterson RS, Riuz DA, editors. Biological Control of Muscoid Flies. Misc Publ, 61:150–63. College Park, Massachusetts: Entomological Society of America.

Unruh TR, Woolley JB. 1999. Molecular Methods in Classical Biological Control In: Van Driesche RG, Bellows TS, Jr., editors. Biological Control, 57–85. New York: Chapman and Hall.

Vanlerberghe-Masutti F. 1994. Detection of genetic variability in Trichogramma populations using molecular markers. Norwegian Journal of Agricultural Sciences (Suppl.) 16: 171–176.

Weir BS. 1990. Genetic Data Analysis: Methods for discrete population genetic data. Sunderland, Massachusetts: Sinauer Associates, Inc. 377 pp.

Weir BS. 1996. Genetic Data Analysis II: Methods for discrete population genetic data. Sunderland, Massachusetts: Sinauer Associates, Inc. 445 pp.

Weir BS, Cockerham CC. 1984. Estimating F-statistics for the analysis of population structure. Evolution 38:1358–1370.

Witsenboer H, Vogel J, Michelmore RW. 1997. Identification, genetic localisation and allelic diversity of amplified microsatellite polymorphic loci in lettuce and wild relatives (Lactuca ssp.). Genome 40: 923–936.

Wolfe AD, Liston A. 1998. Contributions of PCR-based methods to plant systematic and evolutionary biology. In Soltis DE, Soltis PS, Doyle JJ, editors. Molecular Systematics of Plants II: DNA Sequencing, 43–86. New York: Kluwer, New York.

Wu K-S, Jones R, Dannenberg L, Scolnik PA. 1994. Detection of microsatellites without cloning. Nucleic Acids Research 22: 3257–3258.

Yen FC, Yang RC, Mao J, Ye Z, Boyle TJB. 1996. POPGENE, the microsoft windows-based user-friendly software for population genetic analysis of co-dominant and dominant markers and quantitative traits. Dept. Renew. Res., Univ. Alberta, Edmonton, Canada.

Young, DA. 1958. A synopsis of the species of Homalodisca in the United States (Homoptera: Cicadellidae). Bulletin of the Brooklyn Entomological Society. 53: 7–13.

Zietkiewicz E, Rafalski A, Labuda D. 1994. Genomic fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification. Genomics 20:176–183.