Positive selection and sequence rearrangements generate extensive polymorphism in the gamete recognition protein bindin

Bindin is a gamete recognition protein of sea urchins that mediates species-specific attachment of sperm to an egg-surface receptor during fertilization. Sequences of bindin from closely related urchins show fixed species-specific differences. Within species, highly polymorphic bindin alleles result from point substitution, insertion/deletion, and recombination. Since speciation, positive selection favoring allelic variants has generated diversity in bindin polypeptides. Intraspecific bindin variation can be tolerated by the egg receptor, which suggests functional parallels between this system and other flexible recognition systems, including immune recognition. These results show that polymorphism in mate recognition loci required for rapid evolution of sexual isolation can arise within natural populations.      

Bindin from a sea star

SUMMARY The genetic basis for the evolution of development includes genes that encode proteins expressed on the surfaces of sperm and eggs. Previous studies of the sperm acrosomal protein bindin have helped to characterize the adaptive evolution of gamete compatibility and speciation in sea urchins. The absence of evidence for bindin expression in taxa other than the Echinoidea has limited such studies to sea urchins, and led to the suggestion that bindin might be a sea urchin-specific molecule. Here we characterize the gene that encodes bindin in a broadcast-spawning asterinid sea star ( Patiria miniata ). We describe the sequence and domain structure of a full-length bindin cDNA and its single intron. In comparison with sea urchins, P. miniata bindin is larger but the two molecules share several general features of their domain structure and some sequence features of two domains. Our results extend the known evolutionary history of bindin from the Mesozoic (among the crown group sea urchins) into the early Paleozoic (and the common ancestor of eleutherozoans), and present new opportunities for understanding the role of bindin molecular evolution in sexual selection, life history evolution, and speciation among sea stars.     

Low rates of bindin codon evolution in lecithotrophic Heliocidaris sea urchins

Life-history variables including egg size affect the evolutionary response to sexual selection in broadcast-spawning sea urchins and other marine animals. Such responses include high or low rates of codon evolution at gamete recognition loci that encode sperm- and egg-surface peptides. Strong positive selection on such loci affects intraspecific mating success and interspecific reproductive divergence (and may play a role in speciation). Here, we analyze adaptive codon evolution in the sperm acrosomal protein bindin from a brooding sea urchin (Heliocidaris bajulus, with large eggs and nonfeeding or lecithotrophic larval development) and compare our results to previously published data for two closely related congeners. Purifying selection and low relative rates of bindin nonsynonymous substitution in H. bajulus were significantly different from selectively neutral bindin evolution in H. erythrogramma despite similar large egg size in those two species, but were similar to the background rate of nonsynonymous bindin substitution for other closely related sea urchins (including H. tuberculata, all with small egg size and feeding planktonic larval development). Bindin evolution is not driven by egg size variation among Heliocidaris species, but may be more consistent with an alternative mechanism based on the effects of high or low spatial density of conspecific mates.     

The evolution of sea urchin sperm bindin

Sea urchins have been model organisms for the study of fertilization for more than a century. Fertilization in sea urchins happens externally, which facilitates the study of sperm-egg attachment and fusion, and means that all of the molecules involved in gamete recognition and fusion are associated with the gametes. Sea urchin sperm bindin was the first "gamete recognition protein" to be isolated and characterized (Vacquier and Moy 1977), and bindin has since been studied by developmental biologists interested in fertilization, by biochemists interested in membrane fusion and by evolutionary biologists interested in reproductive isolation and speciation. Research on bindin was last reviewed thirteen years ago by Vacquier et al. (1995) in an article titled "What have we learned about sea urchin sperm bindin?" in which the authors reviewed the identification, isolation and early molecular examinations of bindin. Research since then has focused on bindin's potential role in fusing egg and sperm membranes, comparisons of bindin between distantly related species, studies within genera linking bindin evolution to reproductive isolation, and studies within species looking at fertilization effects of individual bindin alleles. In addition, the egg receptor for bindin has been cloned and sequenced. I review this recent research here.     

Reproductive character displacement and the genetics of gamete recognition in tropical sea urchins

Reproductive character displacement occurs when sympatric and allopatric populations of a species differ in traits crucial to reproduction, and it is commonly thought of as a signal of selection acting to limit hybridization. Most documented cases of reproductive character displacement involve characters that are poorly understood at the genetic level, and rejecting alternative hypotheses for biogeographic shifts in reproductive traits is often very difficult. In sea urchins, the gamete recognition protein bindin evolves under positive selection when species are broadly sympatric, suggesting character displacement may be operating in this system. We sampled sympatric and allopatric populations of two species in the sea urchin genus Echinometra for variation in bindin and for the mitochondrial cytochrome oxidase I to examine patterns of population differentiation and molecular evolution at a reproductive gene. We found a major shift in bindin alleles between central Pacific (allopatric) and western Pacific (sympatric) populations of E. oblonga. Allopatric populations of E. oblonga are polyphyletic with E. sp. C at bindin, whereas sympatric populations of the two species are reciprocally monophyletic. There is a strong signal of positive selection (P(N)/P(S) = 4.5) in the variable region of the first exon of bindin, which is associated with alleles found in sympatric populations of E. oblonga. These results indicate that there is a strong pattern of reproductive character displacement between E. oblonga and E. sp. C and that the divergence is driven by selection. There is much higher population structure in sympatric populations at the bindin locus than at the neutral mitochondrial locus, but this difference is not seen in allopatric populations. These data suggest a pattern of speciation driven by selection for local gamete coevolution as a result of interactions between sympatric species. Although this pattern is highly suggestive of speciation by reinforcement, further research into hybrid fitness and egg-sperm interactions is required to address this potential mechanism for character displacement.     

SIMULTANEOUS POSITIVE AND NEGATIVE FREQUENCY‐DEPENDENT SELECTION ON SPERM BINDIN,A GAMETE RECOGNITION PROTEIN IN THE SEA URCHIN STRONGYLOCENTROTUS PURPURATUS

Gamete‐recognition proteins often, but not always, evolve rapidly. We explored how variation in sperm bindin influences reproductive success of the sea urchin Strongylocentrotus purpuratus during group spawning in the sea. Despite large variation in male and female abundance and neighbor distances, males with common genotypes had higher reproductive success than males with rare genotypes. However, males with a relatively uncommon proline‐for‐serine substitution were the most successful. Females also showed a fitness consequence of sperm‐bindin genotype, suggesting linkage disequilibrium between the sperm‐bindin locus and the egg receptor locus. Females with common genotypes had higher reproductive success than rare genotypes, but females with relatively uncommon insertions were most successful. Overall, these results suggest that rare male proteins are selected against, as supported by molecular evidence of purifying selection and probably caused by poor matches to the female receptor protein. Within the pool of moderately common to common alleles, however, individuals with less‐common functional variants were favored and probably maintained by negative frequency‐dependent selection. These results support the hypothesis that sperm availability and sexual conflict influence the evolution of gamete recognition systems in broadcast spawners and highlight the benefits of combining fitness measures with molecular signatures for estimation of patterns of selection.     

Genetic structure of populations of the red sea urchin, Strongylocentrotus franciscanus

Population subdivision was evaluated in the red sea urchin, Strongylocentrotus franciscanus, using DNA sequence data from 134 adult individuals collected in 1995 and 1996. On average 22 individuals were sequenced from six geographic locations between Alaska and Baja California (N=134), nearly the full extent of the species range. DNA sequence data was obtained from direct sequencing of a 273 base pair region of the bindin gene, which encodes a sperm fertilization protein. Results indicate that bindin is sufficiently polymorphic to serve as a genetic marker. We identified 14 unique alleles present in the entire range sampled with a maximum of eight alleles at a specific site. Mean pairwise comparison of the 14 unique alleles indicates moderate sequence diversity (p-distance=1.06). Although there is conflicting evidence to suggest that Alaska populations may deviate from the Hardy-Weinberg expectations, analysis of bindin genotype frequencies indicate that it is not possible to reject the null hypothesis of random mating throughout the species range. The results of a chi-square test with pooling conform to Hardy-Weinberg expectations for all populations (P>0.05) except for the Alaska population (P=0.037). Inbreeding coefficients are consistent with this result and suggest that for the bindin locus, there is high gene flow. These results are compared with previously published results of genetic substructuring in sea urchins to examine relationships among population structure, dispersal potential and biogeography.     

The evolution of gametic compatibility and compatibility groups in the sea urchin Mesocentrotus franciscanus: An avenue for speciation in the sea

The generation of reproductive incompatibility between groups requires a rare genotype with low compatibility to increase in frequency. We tested the hypothesis that sexual conflict driven by the risk of polyspermy can generate compatibility groups in gamete recognition proteins (GRPs) in the sea urchin Mesocentrotus franciscanus. We examined variation in the sperm (bindin) and egg (EBR1) GRPs, how this variation influences fertilization success and how allele frequencies shift in these GRPs over time. The EBR1 gene is a large, 4595 amino acid protein made up of 27 thrombospondin type 1 domain (TSP) and 20 C1s/C1r, uEGF and bone morphogenic protein subdomain (CUB) repeats. Two TSP and two CUB repeats each demonstrate two common non‐synonymous haplotypes (alleles). Sperm bindin and one of these EBR1 repeats (TSP8) shift allele frequencies from one common to two common types over an approximate 200 year interval associated with the removal of predatory sea otters and rising sea urchin abundances; the egg receptor shifts first, followed by the sperm ligand. Laboratory crosses indicate that the historically common sperm and egg gamete recognition proteins have high compatibility as do the new common proteins, with mismatches having lower compatibility. This process of creating compatibility groups sets the stage for reproductive isolation and speciation.     

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     

Ancient Origin of the Null Allele se 428 of the Human ABO-Secretor Locus (FUT2)

In human populations, a null allele having several nucleotide differences from the wild-type allele is segregating at the FUT2 locus (the ABO-Secretor locus) encoding α(1,2)fucosyltransferase. To estimate the age of the most recent common ancestor (MRCA) of these two alleles, we sequenced FUT2 homologues from chimpanzee, gorilla, orangutan, and green monkey. Since we did not detect acceleration or any heterogeneity in the substitution rate at this locus among these species, the age of the MRCA was estimated to be around 3 MYA, assuming the divergence time of human and chimpanzee to be 5 MYA. We developed a simple test to examine whether or not the old age of the MRCA of the FUT2 is consistent with that expected for two divergent neutral alleles sampled from a random mating population. An application of the test to the data at FUT2 indicated that the age of the MRCA is too old to be explained by the simple neutral assumptions, although our test depends on accurate estimation of the divergence time of human and chimpanzee in units of twice the human population size. Various possibilities including balancing selection are discussed to explain this old age of the MRCA. Received: 9 May 1999 / Accepted: 20 September 1999     

Positive Darwinian selection in gamete recognition proteins of Strongylocentrotus sea urchins

Gamete recognition proteins commonly experience positive Darwinian selection and evolve more rapidly than nonreproductive proteins, but the selective forces responsible for their adaptive diversification remain unclear. We examined the patterns of positive selection in the cognate interacting pair of proteins formed by sperm bindin and its egg receptor (EBR1) and in two regions of the sea urchin sperm receptor for egg jelly suREJ3 gene (exons 22 and 26) among four species of Strongylocentrotus sea urchins (S. purpuratus, S. droebachiensis, S. pallidus and S. franciscanus). The signatures of selection differed at each reproductive protein. A strong signal of positive selection was detected at bindin in all lineages even though the species compared had highly variable gamete traits and experience different intensities and forms of sexual selection and sexual conflict in nature. Weaker selection was observed at EBR1 but the small region studied precluded a clear understanding of the extent of sexual conflict between bindin and the EBR1 protein. At the suREJ3 locus, diversifying selection was observed in exon 22 but not exon 26, suggesting that these regions experience different selective pressures and evolutionary constraints. Positive selection was also detected within S. pallidus at suREJ‐22 because of the presence of 12 amino acid replacement mutations segregating at frequencies >0.10. Our results suggest that sexual conflict may be the predominant evolutionary mechanism driving the rapid diversification of reproductive proteins between, and polymorphism within, strongylocentrotid sea urchins.     

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

Abstract 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.     

Special evolutionary properties of genes encoding a protein with a simple amino acid repeat

We have examined the evolution of a gene, SM50, encoding a component of the spicule matrix, which plays an integral role in the formation of the echinoderm skeleton. This gene was originally characterized in Strongylocentrotus purpuratus and encodes an imperfect tandem repeat of six or seven amino acids. We have analyzed the sequence of this repeat in a number of sea urchin species and have determined that the repeat regions have undergone concerted evolution. There are differences in the repeat region between species, but the overall repeat structure is conserved, suggesting the repeat forms a structural domain important in biomineralization. The inherent conserved amino acid repeat structure promotes concerted evolution due to the high probability of misreplication and unequal crossing-over in the repeated segment of the gene. While there are constraints on the amino acids allowed in the repeat region, there are also variations, so that the sequences observed illustrate the balance between amino acid substitutions and concerted evolution. We have evidence that substitutions can alter the mechanisms of unequal crossing-over, altering the way concerted evolution occurs. The way in which concerted evolution occurred appears to be determined by the degree of sequence similarity between the repeats in a given gene, which influences how unequal crossing over may occur. We have mapped the differences in repeat regions on existing phylogenetic trees and indicate where concerted evolution has taken place. We also confirm an earlier report that Hemicentrotus pulcherrimus fits into the Strongylocentrotus genus and examine the evolution of the H. pulcherrimus SM50 repeat relative to other members of this genus. Received: 31 October 2000 / Accepted: 20 March 2001     

Comparative Genomics of Mitochondrial DNA in Drosophila simulans

The current study compares the nucleotide variation among 22 complete mitochondrial genomes of the three distinct Drosophila simulans haplotypes with intron 1 of the alcohol dehydrogenase-related locus. This is the first study to investigate the sequence variation of multiple complete mitochondrial genomes within distinct mitochondrial haplotypes of a single species. Patterns of variation suggest distinct forces are influencing the evolution of mitochondrial DNA (mtDNA) and autosomal DNA in D. simulans. First, there is little variation within each mtDNA haplotype but strong differentiation among them. In contrast, there is no support for differentiation of the mitochondrial haplotypes at the autosomal locus. Second, there is a significant deficiency of mitochondrial variation in each haplotype relative to the autosomal locus. Third, the ratio of nonsynonymous to synonymous substitutions is not equal in all branches of the well-resolved phylogeny. There is an excess of nonsynonymous substitutions relative to synonymous substitutions within each D. simulans haplotype. This result is similar to that previously observed within the mtDNA of distinct species. A single evolutionary force may be causally linked to the observed patterns of mtDNA variation—a rickettsia-like microorganism, Wolbachia pipientis, which is known to directly influence mitochondrial evolution but have a less direct influence on autosomal loci. Received: 16 September 1999 / Accepted: 14 March 2000     

The molecular evolution of sperm bindin in six species of sea urchins (Echinoida: Strongylocentrotidae)

The acrosomal protein bindin attaches sperm to eggs during sea urchin fertilization. Complementary to ongoing functional biochemical studies, I take a comparative approach to explore the molecular evolution of bindin in a group of closely related free-spawning echinoid species. Two alleles of the mature bindin gene were sequenced for each of six species in the sea urchin family Strongylocentrotidae. The nucleotide sequences diverged by at least 1% per Myr at both silent and replacement sites. Two short sections flanking the conserved block show an excess of nonsynonymous substitutions. Each is homologous to a region that had been identified as a target of selection in other sea urchin comparisons. A large proportion of the bindin-coding sequence consists of a highly variable repeat region. Bindin sequences, even including the large intron, could not resolve the branching order among five of the species.      

Evidence for multiple functional copies of the male sex-determining locus, sry, in African murine rodents

Southern hybridization data suggest that the male sex-determining locus, Sry, is often duplicated in rodents. Here we explore DNA sequence evolution of orthologous and paralogous copies of Sry isolated from six species of African murines. PCR amplification followed by direct sequencing revealed from two to four copies of Sry per species. All copies include a long open reading frame, with a stop codon that coincides closely with the stop codon of the house mouse, Mus musculus, a species known to have a single copy of Sry. A phylogenetic analysis suggests that there are at least seven paralogous copies of Sry in this group of rodents. Putative orthologues are identical; sequence divergence among putative paralogues ranges from 1 to 8% (excluding the CAG repeat), with much lower levels of divergence in the high-mobility group (HMG-box) region than in the C-terminal region. A high proportion of nucleotide substitutions in both regions result in amino-acid replacement. The long open reading frame, conserved HMG-box, and pattern of evolution of the putative paralogues suggest that they are functional. Received: 4 October 1996 / Accepted: 17 January 1997     

Molecular Evolution of Cytochrome c Oxidase Subunit IV: Evidence for Positive Selection in Simian Primates

Cytochrome c oxidase (COX) is a multi-subunit enzyme complex that catalyzes the final step of electron transfer through the respiratory chain on the mitochondrial inner membrane. Up to 13 subunits encoded by both the mitochondrial (subunits I, II, and III) and nuclear genomes occur in eukaryotic organisms ranging from yeast to human. Previously, we observed a high number of amino acid replacements in the human COX IV subunit compared to mouse, rat, and cow orthologues. Here we examined COX IV evolution in the two groups of anthropoid primates, the catarrhines (hominoids, cercopithecoids) and platyrrhines (ceboids), as well as one prosimian primate (lorisiform), by sequencing PCR-amplified portions of functional COX4 genes from genomic DNAs. Phylogenetic analysis of the COX4 sequence data revealed that accelerated nonsynonymous substitution rates were evident in the early evolution of both catarrhines and, to a lesser extent, platyrrhines. These accelerated rates were followed later by decelerated rates, suggesting that positive selection for adaptive amino acid replacement became purifying selection, preserving replacements that had occurred. The evidence for positive selection was especially pronounced along the catarrhine lineage to hominoids in which the nonsynonymous rate was first faster than the synonymous rate, then later much slower. The rates of three types of ``neutral DNA' nucleotide substitutions (synonymous substitutions, pseudogene nucleotide substitutions, and intron nucleotide substitutions) are similar and are consistent with previous observations of a slower rate of such substitutions in the nuclear genomes of hominoids than in the nuclear genomes of other primate and mammalian lineages. Received: 22 May 1996 / Accepted: 24 November 1996     

Selection and demographic history shape the molecular evolution of the gamete compatibility protein bindin in Pisaster sea stars

Reproductive compatibility proteins have been shown to evolve rapidly under positive selection leading to reproductive isolation, despite the potential homogenizing effects of gene flow. This process has been implicated in both primary divergence among conspecific populations and reinforcement during secondary contact; however, these two selective regimes can be difficult to discriminate from each other. Here, we describe the gene that encodes the gamete compatibility protein bindin for three sea star species in the genus Pisaster. First, we compare the full‐length bindin‐coding sequence among all three species and analyze the evolutionary relationships between the repetitive domains of the variable second bindin exon. The comparison suggests that concerted evolution of repetitive domains has an effect on bindin divergence among species and bindin variation within species. Second, we characterize population variation in the second bindin exon of two species: We show that positive selection acts on bindin variation in Pisaster ochraceus but not in Pisaster brevispinus, which is consistent with higher polyspermy risk in P. ochraceus. Third, we show that there is no significant genetic differentiation among populations and no apparent effect of sympatry with congeners that would suggest selection based on reinforcement. Fourth, we combine bindin and cytochrome c oxidase 1 data in isolation‐with‐migration models to estimate gene flow parameter values and explore the historical demographic context of our positive selection results. Our findings suggest that positive selection on bindin divergence among P. ochraceus alleles can be accounted for in part by relatively recent northward population expansions that may be coupled with the potential homogenizing effects of concerted evolution.     

Duplication and Diversification of the Apolipoprotein CI (APOCI) Genomic Segment in Association with Retroelements

We have previously shown that several multicopy gene families within the major histocompatibility complex (MHC) arose from a process of segmental duplication. It has also been observed that retroelements play a role in generating diversity within these duplicated segments. The objective of this study was to compare the genomic organization of a gene duplication within another multicopy gene family outside the MHC. Using new continuous genomic sequence encompassing the APOE-CII gene cluster, we show that APOCI and its pseudogene, APOCI′, are contained within large duplicated segments which include sequences from the hepatic control region (HCR). Flanking Alu sequences are observed at both ends of the duplicated unit, suggesting a possible role in the integration of these segments. As observed previously within the MHC, the major differences between the segments are the insertion of sequences (approximately 200–1000 bp in length), consisting predominantly of Alu sequences. Ancestral retroelements also contribute to the generation of sequence diversity between the segments, especially within the 3′ poly(A) tract of Alu sequences. The exonic and regulatory sequences of the APOCI and HCR loci show limited sequence diversity, with exon 3 being an exception. Finally, the typing of pre- and postduplication Alus from both segments indicates an estimated time of duplication of approximately 37 million years ago (mya), some time prior to the separation of Old and New World monkeys. Received: 17 July 1999 / Accepted: 6 November 1999