REPRODUCTIVE CHARACTER DISPLACEMENT AND THE GENETICS OF GAMETE RECOGNITION IN TROPICAL SEA URCHINS

Assortative mating drives linkage disequilibrium between sperm and egg recognition protein loci in the sea urchin Strongylocentrotus purpuratus

Sperm and eggs have interacting proteins on their surfaces that influence their compatibility during fertilization. These proteins are often polymorphic within species, producing variation in gamete affinities. We first demonstrate the fitness consequences of various sperm bindin protein (Bindin) variants in the sea urchin Strongylocentrotus purpuratus, and assortative mating between males and females based on their sperm Bindin genotype. This empirical finding of assortative mating based on sperm Bindin genotype could arise by linkage disequilibrium (LD) between interacting sperm and egg recognition loci. We then examine sequence variation in eight exons of the sea urchin egg receptor for sperm Bindin (EBR1). We find little evidence of LD among the eight exons of EBR1, yet strong evidence for LD between sperm Bindin and EBR1 overall, and varying degrees of LD between sperm Bindin among the eight exons. We reject the alternate hypotheses of LD driven by shared evolutionary histories, population structure, or close physical linkage between these interacting loci on the genome. The most parsimonious explanation for this pattern of LD is that it represents selection driven by assortative mating based on interactions among these sperm and egg loci. These findings indicate the importance of ongoing sexual selection in the maintenance of protein polymorphisms and LD, and more generally highlight how LD can be used as an indication of current mate choice, as opposed to historic selection.

Extreme positive selection on a new highly-expressed larval glycoprotein (LGP) gene in Galaxias fishes (Osmeriformes: Galaxiidae)

We describe the intron-exon structure and DNA/protein sequences of a new larval glycoprotein (LGP) gene from nine species of galaxiid fish. The gene has a distant similarity to Danio THP (Tamm-Horsfall urinary glycoprotein; uromodulin) and cichlid SPP120 (seminal plasma glycoprotein) due to conserved features of its zona pellucida (ZP) domain, including eight highly conserved cysteines and a consensus furin cleavage site. Using a combination of 454 sequencing of cDNA and exon-primed intron-spanning sequencing of genomic DNA, we obtained full sequences of the coding region (996 bp) and its intervening sequences (1,459 bp). LGP shows an exceptionally strong signal of positive selection over the entire coding region, as evidenced by d(N)/d(S) values >1. Across nine species of Galaxias, 87/332 (26%) amino acid residues are variable, compared with 9/386 (2%) for mitochondrial cytochrome b (cytb) in the same group of species. Across 36 interspecific pairwise comparisons, genetic distances are in all cases larger for coding region than for introns, by a factor of 2.4-fold on average. Reading frame, gene structure, splice sites, and many ZP motifs are conserved across all species. Together with the fact that the gene is expressed in all species, these results argue clearly against the possibility of a pseudogene. We show by 454 sequencing and quantitative polymerase chain reaction that the transcript is abundant (ca. 0.5%) in newly hatched larvae and appears to be almost absent from a range of adult tissues. We postulate that the strong Darwinian evolution exhibited by this protein may reflect some type of immunoprotection at this vulnerable larval stage.

Extraordinary intraspecific diversity in oyster sperm bindin

In free-spawning invertebrates sperm-egg incompatibility is a barrier to mating between species, and divergence of gamete recognition proteins (GRPs) can result in reproductive isolation. Of interest are processes that create reproductive protein diversity within species, because intraspecific variants are potentially involved in mate choice and early speciation. Sperm acrosomes of the Pacific oyster Crassostrea gigas contain the protein bindin that bonds sperm to egg during fertilization. Oyster bindin is a single-copy gene encoding a diversity of protein variants. Oyster bindins have a conserved N-terminal region followed by one to five tandem fucose-binding lectin (F-lectin) domains. These repeats have diversified by positive selection at eight sites clustered on the F-lectin's fucose binding face. Additional bindin variants result from recombination in an intron in each F-lectin repeat. Males also express alternatively spliced bindin cDNAs with one to five repeats, but typically translate only one or two isoforms into protein. Thus, positive selection, alternative splicing, and recombination can create thousands of bindin variants within C. gigas. Models of sexual conflict predict high male diversity when females are diverse and sexual conflict is strong. The amount of intraspecific polymorphism in male GRPs may be a consequence of the relative efficiency of local (molecular recognition) and global (electrical, cortical, and physical) polyspermy blocks that operate during fertilization.

Adaptive gamete-recognition divergence in a hybridizing Mytilus population

Gamete-recognition proteins often evolve rapidly, but it is not known if their divergence occurs within species and corresponds with the evolution of reproductive isolation, or if divergence typically accumulates between already isolated lineages. We examined the evolution of a candidate gamete-recognition protein in several sympatric and allopatric populations of Mytilus blue mussels, species that hybridize in nature. Within a single species, Mytilus galloprovincialis, we found adaptive divergence of Lysin-M7, a sperm acrosomal protein that dissolves the egg vitelline envelope during fertilization. Mytilus galloprovincialis Lysin-M7 alleles group into two distinct clades (termed G and G(D)), and individual alleles in these clades are separated from each other by at least three and up to eleven amino-acid substitutions. Maximum-likelihood estimates of selective pressure (dN/dS =omega) implicate selection in the divergence between M. galloprovincialis Lysin-M7 clades, and within the G(D) clade. Exact tests of population differentiation indicate that the relative frequency of G and G(D) Lysin-M7 alleles differs significantly among M. galloprovincialis populations. Compared with allopatric Mediterranean samples, Lysin-M7 alleles in the G(D) clade are found at elevated frequency in samples from the East Atlantic and California, areas of secondary contact and hybridization between Mytilus species, and Australia, an area of unknown species composition. Adaptive divergence between the alleles most common in allopatry and those found at elevated frequency in samples from sympatry suggests that selection pressures acting in hybridizing populations, likely following Pleistocene secondary contact with M. edulis in the East Atlantic, drove the divergence of Lysin-M7 in M. galloprovincialis.

The evolution of conspecific gamete precedence and its effect on reinforcement

Conspecific gamete precedence, the usage of conspecific sperm by a female that mates with both a conspecific and a heterospecific male, has been found in many taxa. We construct a population genetic model to examine the evolution of conspecific gamete precedence and its coevolution with premating isolation in the process of reinforcement. Our findings suggest that conspecific gamete precedence can evolve via a process very similar to reinforcement. We explore the nature of the selection against hybridization necessary to drive this evolution. Moreover, our results confirm the prediction of Marshall et al. (Trends Ecol. Evol. 2002;17:558-563) that conspecific gamete precedence will inhibit the evolution of reinforcement between two species. We further find that reinforcement will inhibit the evolution of conspecific gamete precedence. Both reinforcement and conspecific gamete precedence increase reproductive isolation and contribute to the process of speciation. We discuss factors that may affect which of these phenomena are likely to become predominant between incipient species.

Evolutionary animation: how do molecular phylogenies compare to Mayr's reconstruction of speciation patterns in the sea?

Ernst Mayr used the geography of closely related species in various stages of increasing divergence to "animate" the process of geographic, or allopatric, speciation. This approach was applied to a wide set of taxa, and a seminal paper by Mayr used it to explore speciation patterns in tropical sea urchins. Since then, taxonomic information in several of these genera has been augmented by detailed molecular phylogenies. We compare Mayr's animation with the phylogenies of eight sea urchin genera placed by Mayr into four speciation groups. True to Mayr's predictions, early-stage genera have on average lower species divergence and more polytypic species than genera in later stages. For six of these genera, we also have information about the evolution of the gamete recognition protein bindin, which is critical to reproductive isolation. These comparisons show that later-stage genera with many sympatric species tend to be those with rapid bindin evolution. By contrast, early-stage genera with few sympatric species are not necessarily earlier in the divergence process; they happen to be those with slow rates of bindin evolution. These results show that the rate of speciation in sea urchins does not only depend on the steady accumulation of genome divergence over time, but also on the rate of evolution of gamete recognition proteins. The animation method used by Mayr is generally supported by molecular phylogenies. However, the existence of multiple rates in the acquisition of reproductive isolation complicates placement of different genera in an evolutionary series.

Speciation on the coasts of the new world: phylogeography and the evolution of bindin in the sea urchin genus Lytechinus

Beginning with E. Mayr's study in 1954, tropical sea urchins have played an important role in studies of speciation in the sea, but what are the processes of cladogenesis and divergence that give rise to new species in this group? We attempt to answer this question in the genus Lytechinus. Unlike the majority of other tropical sea urchin genera, which have circumtropical distributions, Lytechinus is mostly confined to the tropics and subtropics of the New World. We sequenced a region of mitochondrial cytochrome oxidase I and the entire molecule of nuclear bindin (a sperm gamete recognition protein) of nearly all species in the genus, and we assayed isozymes of three partially sympatric closely related species and subspecies. We found that in both mitochondrial DNA (mtDNA) and in bindin the genus Lytechinus is paraphyletic, encompassing Sphaerechinus granularis as the sister species of L. euerces. The rest of the species are arranged in an Atlantic clade composed of L. williamsi and L. variegatus, and a Pacific clade containing L. anamesus, L. pictus, L. semituberculatus, and L. panamensis. Divergence between these clades suggests that they were separated no later than the closure of the Isthmus of Panama, and possibly before this time. Our data confirm that L. anamesus and L. pictus from California are a single species, and provide no evidence of differentiation between L. variegatus variegatus from the Caribbean and L. variegatus atlanticus from Bermuda. Lytechinus variegatus variegatus mtDNA is distinct from that of L. variegatus carolinus from the North American seaboard and the Gulf of Mexico, whereas their bindins are very similar. However, there is clear evidence of introgression of mtDNA between the two subspecies and they share alleles in all sampled isozyme loci. Lytechinus williamsi from the Caribbean shares mtDNA haplotypes with L. variegatus variegatus, and they also share isozymes in all assayed loci. Their bindin, however, is distinct and coalesces within each morphospecies. A private clade of mtDNA in L. williamsi may be indicative of former differentiation in the process of being swamped by introgression, or of recent speciation. Recent sudden expansions in effective population size may explain the predominance of a few mitochondrial haplotypes common to the two species. Despite the high divergence of bindin (relative to differentiation of mtDNA) between L. variegatus and L. williamsi, comparison of amino acid replacement to silent substitutions by various methods uncovered no evidence for positive selection on the bindin of any clade of Lytechinus. With the possible exception of L. williamsi and L. variegatus, our results are consistent with a history of allopatric speciation in Lytechinus. The molecular results from Lytechinus, along with those of other similar studies of sea urchins, suggest that the general speciation patterns deduced in the middle of last century by Mayr from morphology and geography have held up, but also have uncovered peculiarities in the evolution of each genus.