Duplication of the Rdl GABA receptor subunit gene in an insecticide-resistant aphid, Myzus persicae

Resistance to cyclodiene insecticides is associated with replacements of a single amino acid (alanine 302) in a γ-aminobutyric acid (GABA) receptor subunit encoded by the single-copy gene Resistance to dieldrin (Rdl). Alanine 302 is predicted to reside within the second membrane-spanning region of the Rdl receptor, a region that is thought to line the integral chloride ion channel pore. In all cyclodiene-resistant insects studied to date, this same alanine residue is replaced either by a serine, or, in some resistant strains of Drosophila simulans, a glycine residue. Therefore, individuals can carry only two different Rdl alleles. In contrast, here we report the presence of up to four different Rdl-like alleles in individual clones of the green peach aphid, Myzus persicae. In addition to the wild-type copy of Rdl gene (encoding A302 or allele A), M. persicae carries three other alleles with the following amino acid replacements: A302 → Glycine (allele G), A302 → Serine^TCG (allele S) and A302 → Serine^AGT (allele S′). Evidence from direct nucleotide sequencing and Single Stranded Conformational Polymorphism (SSCP) analysis shows that at least three of these different Rdl alleles (i.e. A, G and S) are commonly present in individual aphids or aphid clones. Southern analysis using allele-specific probes and analysis of sequences downstream of the exon containing the resistance-associated mutation confirm the presence of two independent Rdl-like loci in M. persicae. One locus carries the susceptible alanine (A) and/or resistant glycine (G) allele while the other carries the two serine alleles (S or S′). Whereas resistance levels are correlated with the glycine replacement, the S allele was present in all aphid clones, regardless of their resistance status. These results suggest that target site insensitivity is associated with replacements at the first (A/G) but not the second (S/S′) locus. Phylogenetic analysis of nucleotide sequences indicates that both putative aphid Rdl loci are monophyletic with respect to other insect Rdl genes and may have arisen through a recent gene duplication event. The implications of this duplication with respect to insecticide resistance and insect GABA receptor subunit diversity are discussed. Received: 10 March 1998 / Accepted: 21 July 1998     

Duplication of the Rdl GABA receptor subunit gene in an insecticide-resistant aphid, Myzus persicae

Resistance to cyclodiene insecticides is associated with replacements of a single amino acid (alanine 302) in a γ-aminobutyric acid (GABA) receptor subunit encoded by the single-copy gene Resistance to dieldrin (Rdl). Alanine 302 is predicted to reside within the second membrane-spanning region of the Rdl receptor, a region that is thought to line the integral chloride ion channel pore. In all cyclodiene-resistant insects studied to date, this same alanine residue is replaced either by a serine, or, in some resistant strains of Drosophila simulans, a glycine residue. Therefore, individuals can carry only two different Rdl alleles. In contrast, here we report the presence of up to four different Rdl-like alleles in individual clones of the green peach aphid, Myzus persicae. In addition to the wild-type copy of Rdl gene (encoding A302 or allele A), M. persicae carries three other alleles with the following amino acid replacements: A302?→?Glycine (allele G), A302?→ Serine^TCG (allele S) and A302?→?Serine^AGT (allele S′). Evidence from direct nucleotide sequencing and Single Stranded Conformational Polymorphism (SSCP) analysis shows that at least three of these different Rdl alleles (i.e. A, G and S) are commonly present in individual aphids or aphid clones. Southern analysis using allele-specific probes and analysis of sequences downstream of the exon containing the resistance-associated mutation confirm the presence of two independent Rdl-like loci in M. persicae. One locus carries the susceptible alanine (A) and/or resistant glycine (G) allele while the other carries the two serine alleles (S or S′). Whereas resistance levels are correlated with the glycine replacement, the S allele was present in all aphid clones, regardless of their resistance status. These results suggest that target site insensitivity is associated with replacements at the first (A/G) but not the second (S/S′) locus. Phylogenetic analysis of nucleotide sequences indicates that both putative aphid Rdl loci are monophyletic with respect to other insect Rdl genes and may have arisen through a recent gene duplication event. The implications of this duplication with respect to insecticide resistance and insect GABA receptor subunit diversity are discussed.     

The role of Rdl in resistance to phenylpyrazoles in Drosophila melanogaster

Extensive use of older generation insecticides may result in pre-existing cross-resistance to new chemical classes acting at the same target site. Phenylpyrazole insecticides block inhibitory neurotransmission in insects via their action on ligand-gated chloride channels (LGCCs). Phenylpyrazoles are broad-spectrum insecticides widely used in agriculture and domestic pest control. So far, all identified cases of target site resistance to phenylpyrazoles are based on mutations in the Rdl (Resistance to dieldrin) LGCC subunit, the major target site for cyclodiene insecticides. We examined the role that mutations in Rdl have on phenylpyrazole resistance in Drosophila melanogaster, exploring naturally occurring variation, and generating predicted resistance mutations by mutagenesis. Natural variation at the Rdl locus in inbred strains of D. melanogaster included gene duplication, and a line containing two Rdl mutations found in a highly resistant line of Drosophila simulans. These mutations had a moderate impact on survival following exposure to two phenylpyrazoles, fipronil and pyriprole. Homology modelling suggested that the Rdl chloride channel pore contains key residues for binding fipronil and pyriprole. Mutagenesis of these sites and assessment of resistance in vivo in transgenic lines showed that amino acid identity at the Ala³⁰¹ site influenced resistance levels, with glycine showing greater survival than serine replacement. We confirm that point mutations at the Rdl 301 site provide moderate resistance to phenylpyrazoles in D. melanogaster. We also emphasize the beneficial aspects of testing predicted mutations in a whole organism to validate a candidate gene approach.     

Sequence of PRAT Satellite DNA ``Frozen' in Some Coleopteran Species

The intriguing diversity of highly abundant satellite repeats found even among closely related species can result from processes leading to dramatic changes in copy number of a particular sequence in the genome and not from rapid accumulation of mutations. To test this hypothesis, we investigated the distribution of the PRAT satellite DNA family, a highly abundant major satellite in the coleopteran species Palorus ratzeburgii, in eight species belonging to the related genera (Tribolium, Tenebrio, Latheticus), the subfamily (Pimeliinae), and the family (Chrysomelidae). Dot blot analysis and PCR assay followed by Southern hybridization revealed that the PRAT satellite, in the form of low-copy number repeats, was present in all tested species. The PRAT satellite detected in the species Pimelia elevata has been sequenced, and compared with previously cloned PRAT monomers from Palorus ratzeburgii and Palorus subdepressus. Although the two Palorus species diverged at least 7 Myr ago, and the subfamily Pimeliinae separated from the genus Palorus 50–60 Myr ago, all PRAT clones exhibit high mutual homology, with average variability relative to the common consensus sequence of 1.3%. The presence of ancestral mutations found in PRAT clones from all three species as well as the absence of species diagnostic mutations illustrate extremely slow sequence evolution. This unexpectedly high conservation of PRAT satellite DNA sequence might be induced by a small bias of turnover mechanisms favoring the ancestral sequence in the process of molecular drive.     

Evolution of Microsatellite Alleles in Four Species of Mice (Genus Apodemus)

Microsatellite length variation was investigated at a highly variable microsatellite locus in four species of Apodemus. Information obtained from microsatellite allele sequences was contrasted with allele sizes, which included 18 electromorphs. Additional analysis of a 400-bp unique sequence in the flanking region identified 26 different haplotype sequences or ``true' alleles in the sample. Three molecular mechanisms, namely, (1) addition/deletion of repeats, (2) substitutions and indels in the flanking region, and (3) mutations interrupting the repeat, contributed to the generation of allelic variation. Size homoplasy can be inferred for alleles within populations, from different populations of the same species, and from different species. We propose that microsatellite flanking sequences may be informative markers for investigating mutation processes in microsatellite repeats as well as phylogenetic relationships among alleles, populations, and species. Received: 3 November 1999 / Accepted: 2 May 2000     

A system to efficiently maintain embryonic lethal mutations in the flour beetle Tribolium castaneum

 Due to its small size, short life cycle, and easy maintenance, the flour beetle Tribolium castaneum is well suited for the genetic analysis of development. One drawback of Tribolium as a genetic system is, however, the difficulty of keeping embryonic lethal lines. Presently, only few lethal mutations can be kept as balanced stocks. Therefore, heterozygous carriers must be identified anew in every generation in order to maintain a recessive embryonic mutation. To alleviate this problem we have devised a block system that allows the simultaneous processing of many mutant lines or test crosses for visual inspection of larval cuticle phenotypes. Using this technique, one person can maintain about 100 embryonic lethal stocks, which makes feasible the thorough genetic analysis of embryogenesis in this species. Received: 4 November 1998 / Accepted: 7 February 1999     

Allelic variation in the squirrel monkey x-linked color vision gene: biogeographical and behavioral correlates

Most Neotropical primate species possess a polymorphic X-linked and a monomorphic autosomal color vision gene. Consequently, populations are composed of both dichromatics and trichromatics. Most theories on the maintenance of this genetic system revolve around possible advantages for foraging ecology. To examine the issue from a different angle, we compared the numbers and relative frequencies of alleles at the X-linked locus among three species of Saimiri representing a wide range of geographical and behavioral variation in the genus. Exons 3, 4, and 5 of the X-linked opsin gene were sequenced for a large number of X chromosomes for all three species. Several synonymous mutations were detected in exons 4 and 5 for the originally reported alleles but only a single nonsynonymous change was detected. Two alleles were found that appeared to be the result of recombination events. The low occurrence of recombinant alleles and absence of mutations in the amino acids critical for spectral tuning indicates that stabilizing selection acts to maintain the combinations of critical sites specific to each allele. Allele frequencies were approximately the same for all Saimiri species, with a slight but significant difference between S. boliviensis and S. oerstedii. No apparent correlation exists between allele frequencies and behavioral or biogeographical differences between species, casting doubt on the speculation that the spectral sensitivities of the alleles have been maintained because they are specifically well-tuned to Saimiri visual ecology. Rather, the spectral tuning peaks might have been maintained because they are as widely spaced as possible within the limited range of middlewave to longwave spectra useful to all primates. This arrangement creates a balance between maximizing the distance between spectral tuning peaks (allowing the color opponency of the visual system to distinguish between peaks) and maximizing the number of alleles within a limited range (yielding the greatest possible frequency of heterozygotes).     

The Conserved Motif in Hydrophilic Loop 2/3 and Loop 8/9 of the Lactose Permease of Escherichia coli. Analysis of Suppressor Mutations

The major facilitator superfamily (MFS) of transport proteins, which includes the lactose permease of Escherichia coli, contains a conserved motif G-X-X-X-D/E-R/K-X-G-R/K-R/K in the loops that connect transmembrane segments 2 and 3, and transmembrane segments 8 and 9. In three previous studies (Jessen-Marshall, A.E., & Brooker, R.J. 1996. J. Biol. Chem. 271:1400–1404; Jessen-Marshall, A.E., Parker, N., & Brooker, R.J. 1997. J. Bacteriol. 179:2616–2622; and Pazdernik, N., Cain, S.M., & Brooker, R.J. 1997. J. Biol. Chem. 272:26110–26116), suppressor mutations at twenty different sites were identified which restore function to mutant permeases that have deleterious mutations in the conserved loop 2/3 or loop 8/9 motif. In the current study, several of these second-site suppressor mutations have been separated from the original mutation in the conserved motif. The loop 2/3 suppressors were then coupled to a loop 8/9 mutation (P280L) and the loop 8/9 suppressors were coupled to a loop 2/3 mutation (i.e., G64S) to determine if the suppressors could restore function only to a loop 2/3 mutation, a loop 8/9 mutation, or both. The single parent mutations changing the first position in loop 2/3 (i.e., G64S) and loop 8/9 (i.e., P280L) had less than 4% lactose transport activity. Interestingly, most of the suppressors were very inhibitory when separated from the parent mutation. Two suppressors, A50T and G370V, restored substantial transport activity when individually coupled to the mutation in loop 2/3 and also when coupled to the corresponding mutation in loop 8/9. In other words, these suppressors could alleviate a defect imposed by mutations in either half of the permease. From a kinetic analysis, these suppressors were shown to exert their effects by increasing the V _max values for lactose transport compared with the single G64S and P280L strains. These results are discussed within the context of our model in which the two halves of the lactose permease interact at a rotationally symmetrical interface, and that lactose transport is mediated by conformational changes at the interface. Received: 18 November 1999/Revised: 11 April 2000     

Target site mediated insecticide resistance: what questions remain?

Recent advances in the molecular genetics of insecticide resistance have identified the point mutations associated with target site insensitivity in the genes encoding the three major insecticide targets: the Rdl GABA receptor, the para voltage gated sodium channel and insect acetylcholinesterase. However central questions relating to the origin, selection and fitness of these mutations in natural populations remain. This review examines the extent to which we understand how a specific subset of potential resistance associated mutations are selected, how often they may arise and/or recombine, and whether we can explain any potential fitness disadvantages based on our knowledge of the molecular mechanisms of resistance involved.     

Forty years of erratic insecticide resistance evolution in the mosquito Culex pipiens

One view of adaptation is that it proceeds by the slow and steady accumulation of beneficial mutations with small effects. It is difficult to test this model, since in most cases the genetic basis of adaptation can only be studied a posteriori with traits that have evolved for a long period of time through an unknown sequence of steps. In this paper, we show how ace-1, a gene involved in resistance to organophosphorous insecticide in the mosquito Culex pipiens, has evolved during 40 years of an insecticide control program. Initially, a major resistance allele with strong deleterious side effects spread through the population. Later, a duplication combining a susceptible and a resistance ace-1 allele began to spread but did not replace the original resistance allele, as it is sublethal when homozygous. Last, a second duplication, (also sublethal when homozygous) began to spread because heterozygotes for the two duplications do not exhibit deleterious pleiotropic effects. Double overdominance now maintains these four alleles across treated and nontreated areas. Thus, ace-1 evolution does not proceed via the steady accumulation of beneficial mutations. Instead, resistance evolution has been an erratic combination of mutation, positive selection, and the rearrangement of existing variation leading to complex genetic architecture.     

Enhanced Evolvability in Immunoglobulin V Genes Under Somatic Hypermutation

Darwinian theory requires that mutations be produced in a nonanticipatory manner; it is nonetheless consistent to suggest that mutations that have repeatedly led to nonviable phenotypes would be introduced less frequently than others—if under appropriate genetic control. Immunoglobulins produced during infection acquire point mutations that are subsequently selected for improved binding to the eliciting antigen. We and others have speculated that an enhancement of mutability in the complementarity-determining regions (CDR; where mutations have a greater chance of being advantageous) and/or decrement of mutability in the framework regions (FR; where mutations are more likely to be lethal) may be accomplished by differential codon usage in concert with the known sequence specificity of the hypermutation mechanism. We have examined 115 nonproductively rearranged human Ig sequences. The mutation patterns in these unexpressed genes are unselected and therefore directly reflect inherent mutation biases. Using a χ² test, we have shown that the number of mutations in the CDRs is significantly higher than the number of mutations found in the FRs, providing direct evidence for the hypothesis that mutations are preferentially targeted into the CDRs.     

Mutation in ace1 associated with an insecticide resistant population of Plutella xylostella

Insensitive acetylcholinesterase (AChE) was determined to be involved in an EPN-resistant (ER) strain and a contaminated susceptible (CS) strain of diamondback moth (DBM, Plutella xylostella L.), as estimated by AChE inhibition assay using DDVP as a inhibitor in a nondenaturing electrophoresis gel. The ER strain exhibited very high AChE insensitivity, high resistance ratio, and two point mutations (G324A, A298S) in ace1-type AChE gene (Pxace1). The CS strain showed low AChE insensitivity, low resistance ratio, and it has only one point mutation (G324A). These findings suggest that the A298S mutation, along with reported G324A mutation (Baek et al, 2005), can be important in the development of organophosphate resistance. These results also suggest that the A298S mutation could be a good candidate for a molecular diagnosis marker for resistance monitoring. Three molecular diagnosis methods (Quantitative Sequencing; QS, PCR amplification of specific alleles; PASA and restriction fragment length polymorphism; RFLP) were developed which successfully detected specific resistance associated point mutations. Seven local population DBMs were surveyed and showed high insecticide resistance levels and a A298S mutation in Pxace1. These methods can be used to monitor the resistance allele in field population of DBMs and resistance management strategy.     

The Parthenogenetic Rock Lizard Lacerta unisexualis: An Example of Limited Genetic Polymorphism

Protein electrophoresis of Lacerta unisexualis from three populations found that 21 of 36 allozyme loci were homozygous, while 14 expressed fixed heterozygotes and one locus was variable. Three clones were detected at the locus Cat-A. Two individuals represent two rare clones while all others form a common clone. Our favored explanation is the mutation of a preexisting common clone rather than multiple origins. Received: 27 January 1997 / Accepted: 29 April 1997     

Expected Relationship Between the Silent Substitution Rate and the GC Content: Implications for the Evolution of Isochores

The relationship between the silent substitution rate (K _s) and the GC content along the genome is a focal point of the debate about the origin of the isochore structure in vertebrates. Recent estimation of the silent substitution rate showed a positive correlation between K _s and GC content, in contradiction with the predictions of both the regional mutation bias model and the selection or biased gene conversion model. The aim of this paper is to help resolve this contradiction between theoretical studies and data. We analyzed the relationship between K _s and GC content under (1) uniform mutation bias, (2) a regional mutation bias, and (3) mutation bias and selection. We report that an increase in K _s with GC content is expected under mutation bias because of either nonequilibrium of the isochore structure or an increasing mutation rate from AT toward GC nucleotides in GC-richer isochores. We show by simulations that CpG deamination tends to increase the mutation rate with GC content in a regional mutation bias model. We also demonstrate that the relationship between K _s and GC under the selectionist or biased gene conversion model is positive under weak selection if the mutation selection equilibrium GC frequency is less than 0.5. Received: 28 March 2001 / Accepted: 16 May 2001     

Distinct roles of two RDL GABA receptors in fipronil action in the diamondback moth (Plutella xylostella)

The phenylpyrazole insecticide fipronil blocks resistance to dieldrin (RDL) γ-aminobutyric acid (GABA) receptors in insects, thereby impairing inhibitory neurotransmission. Some insect species, such as the diamondback moth (Plutella xylostella), possess more than one Rdl gene. The involvement of multiple Rdls in fipronil toxicity and resistance remains largely unknown. In this study, we investigated the roles of two Rdl genes, PxRdl1 and PxRdl2, in P. xylostella fipronil action. In Xenopus oocytes, PxRDL2 receptors were 40 times less sensitive to fipronil than PxRDL1. PxRDL2 receptors were also less sensitive to GABA compared with PxRDL1. Knockout of the fipronil-sensitive PxRdl1 reduced the fipronil potency 10-fold, whereas knockout of the fipronil-resistant PxRdl2 enhanced the fipronil potency 4.4-fold. Furthermore, in two fipronil-resistant diamondback moth field populations, PxRdl2 expression was elevated 3.7- and 4.1-fold compared with a susceptible strain, whereas PxRdl1 expression was comparable among the resistant and susceptible strains. Collectively, our results indicate antagonistic effects of PxRDL1 and PxRDL2 on fipronil action in vivo and suggest that enhanced expression of fipronil-resistant PxRdl2 is potentially a new mechanism of fipronil resistance in insects.     

Evolutionary Studies on Uricases of Fungal Endosymbionts of Aphids and Planthoppers

Aphids belonging to the three genera Tuberaphis, Glyphinaphis, and Cerataphis contain extracellular fungal symbionts that resemble endocellular yeast-like symbionts of planthoppers. Whereas the symbiont of planthoppers has a uricase (urate oxidase; EC 1.7.3.3) and recycles uric acid that the host stores, no uric acid was found in Tuberaphis styraci, and its fungal symbiont did not exhibit the uricase activity. However, the fungal symbionts of these aphids, including that of T. styraci, were shown to have putative uricase genes, or pseudogenes, for the uricase. Sequence analysis of these genes revealed that deleterious mutations occurred independently on each lineage of Glyphinaphis and Tuberaphis, while no such mutation was found in the lineage of Cerataphis. These genes were almost identical to those cloned from the symbionts of planthoppers, though the host aphids and planthoppers are phylogenetically distant. To estimate the phylogenetic relationship in detail between the fungal symbionts of aphids and those of planthoppers, a gene tree was constructed based on the sequences of the uricase genes including their flanking regions. As a result, the symbionts of planthoppers and Tuberaphis aphids formed a sister group against those of Glyphinaphis and Cerataphis aphids with high bootstrap confidence levels, which strongly suggests that symbionts have been horizontally transferred from the aphids' lineage to the planthoppers'. Received: 29 March 2000 / Accepted: 31 May 2000     

Selection of dieldrin-resistant strains of Lucilia cuprina (Diptera: Calliphoridae) after ethyl methanesulfonate mutagenesis of a susceptible strain.

After ethyl methanesulfonate (EMS) mutagenesis of a susceptible strain (SWT), selective screening of Lucilia cuprina (Wiedemann) resulted in four strains that were resistant to the insecticide dieldrin. Concentrations used for selection were greater than LC99 of susceptible phenotypes. No resistant variants were screened from the standard laboratory strain (SWT) not treated with EMS. The resistance phenotypes of the four resistant strains were similar to each other and to that of a field-selected resistant strain. The genetic basis of resistance is monogenic in all strains and the data are consistent with the same locus, Rdl, determining resistance status in each strain. The Rdl locus maps to chromosome V, approximately 3.5 map units distal to the Sut locus. Dieldrin resistance may be caused by less effective blocking of insect neuronal GABA receptors by the chemical in resistant strains. The data indicate that the evolution of resistance to an insecticide in the field may be constrained by a limited number of genetical and biochemical options if a monogenic response is selected for and that the spontaneous mutation rate to the Rdl allele is less than 1 in 10(6) in the laboratory.     

Divergent Human Y-Chromosome Microsatellite Evolution Rates

In this work, we analyze several characteristics influencing the low variability of the microsatellite DYS19 in the major founder Amerindian Y chromosome lineage containing the point mutation DYS199-T. Variation of DYS19 was compared with that of five other Y-linked tetranucleotide repeat loci (DYS389A, DYS389B, DYS390, DYS391, and DYS393) in the DYS199-T lineage. All the other microsatellites showed significantly higher levels of variability than DYS19 as measured by gene diversity and repeat number variance. Moreover, we had previously shown that DYS19 had high diversity in Brazilians and in several other populations worldwide. Thus, the slow DYS19 evolution in the DYS199-T lineage seems to be both locus and allele specific. To understand the slow DYS19 evolutionary rate, the microsatellite loci were compared according to their mapping on the Y chromosome and also on the basis of structural aspects such as the base composition of the repeat motif and flanking regions and the degree of perfection and size (repeat number) of the variable blocks. The only observed difference that might be related to the low DYS19 variability is its small average number of repeats, a value expected to be closer to the founder DYS19 allele in the DYS199-T lineage. These data were also compared to other derived Y lineages. The Tat-C lineage displayed a lower DYS19 variability correlated to a small average repeat number, while in the DYS234-G lineage, a high DYS19 variability was found associated to a larger average repeat number. This approach reveals that evolution of Y microsatellites in lineages defined by slowly evolving markers, such as point mutations, can be greatly influenced by the size (number of repeats of the variable block) of the founder allele in each microsatellite locus. Thus lineage-dating methods using microsatellite variation should be practiced with great care. Received: 7 November 1998 / Accepted: 9 April 1999     

The <Emphasis Type="Italic">Tribolium spineless</Emphasis> ortholog specifies both larval and adult antennal identity

The morphology of insect antennae varies widely among species, but our understanding of antennal development comes almost solely from studies of one species—the fruit fly, Drosophila melanogaster. Moreover, this knowledge applies mostly to adult structures, since Drosophila lacks external larval appendages. In contrast to Drosophila, the red flour beetle, Tribolium castaneum, has both larval and adult antennae, which are very different from one another in morphology. Thus, Tribolium provides an ideal system to compare modes of antennal development both within and between species. Here, we report that the Tribolium ortholog of spineless (Tc-ss) is required in both the larval and adult antennae. Knockdown of Tc-ss by RNAi during either larval or imaginal development causes transformation of the distal portion of the antennae to legs. Thus, the function of ss is conserved between Drosophila and Tribolium with respect to adult antennal specification and also between Tribolium larval and adult antennal development. The similarity of the Tc-ss RNAi phenotype to that of a classically described Tribolium mutation, antennapedia (ap) (of no relationship to the Drosophila Hox gene of the same name), led us to characterize the original ap mutation and two newly identified ap alleles. Our mapping and phenotypic data suggest that Tc-ss is the best candidate for the ap locus. These results represent a first step in characterizing larval and adult antennal patterning in Tribolium, which should provide important insights into the evolution of insect antennal development.     

Evaluation of Sequence Variation and Selection in the Bindin Locus of the Red Sea Urchin, Strongylocentrotus franciscanus

Recent evidence suggests that gamete recognition proteins may be subjected to directed evolutionary pressure that enhances sequence variability. We evaluated whether diversity enhancing selection is operating on a marine invertebrate fertilization protein by examining the intraspecific DNA sequence variation of a 273-base pair region located at the 5′ end of the sperm bindin locus in 134 adult red sea urchins (Strongylocentrotus franciscanus). Bindin is a sperm recognition protein that mediates species-specific gamete interactions in sea urchins. The region of the bindin locus examined was found to be polymorphic with 14 alleles. Mean pairwise comparison of the 14 alleles indicates moderate sequence diversity (p-distance = 1.06). No evidence of diversity enhancing selection was found. It was not possible to reject the null hypothesis that the sequence variation observed in S. franciscanus bindin is a result of neutral evolution. Statistical evaluation of expected proportions of replacement and silent nucleotide substitutions, observed versus expected proportions of radical replacement substitutions, and conformance to the McDonald and Kreitman test of neutral evolution all indicate that random mutation followed by genetic drift created the polymorphisms observed in bindin. Observed frequencies were also highly similar to results expected for a neutrally evolving locus, suggesting that the polymorphism observed in the 5′ region of S. franciscanus bindin is a result of neutral evolution. Received: 19 June 1998 / Accepted: 2 August 2000