Natural hybridization in the sea urchin genus Pseudoboletia between species without apparent barriers to gamete recognition

Marine species with high dispersal potential often have huge ranges and minimal population structure. Combined with the paucity of geographic barriers in the oceans, this pattern raises the question as to how speciation occurs in the sea. Over the past 20 years, evidence has accumulated that marine speciation is often linked to the evolution of gamete recognition proteins. Rapid evolution of gamete recognition proteins in gastropods, bivalves, and sea urchins is correlated with gamete incompatibility and contributes to the maintenance of species boundaries between sympatric congeners. Here, we present a counterexample to this general pattern. The sea urchins Pseudoboletia indiana and P. maculata have broad ranges that overlap in the Indian and Pacific oceans. Cytochrome oxidase I sequences indicated that these species are distinct, and their 7.3% divergence suggests that they diverged at least 2 mya. Despite this, we suspected hybridization between them based on the presence of morphologically intermediate individuals in sympatric populations at Sydney, Australia. We assessed the opportunity for hybridization between the two species and found that (1) individuals of the two species occur within a meter of each other in nature, (2) they have overlapping annual reproductive cycles, and (3) their gametes cross-fertilize readily in the laboratory and in the field. We genotyped individuals with intermediate morphology and confirmed that many were hybrids. Hybrids were fertile, and some female hybrids had egg sizes intermediate between the two parental species. Consistent with their high level of gamete compatibility, there is minimal divergence between P. indiana and P. maculata in the gamete recognition protein bindin, with a single fixed amino acid difference between the two species. Pseudoboletia thus provides a well-characterized exception to the idea that broadcast spawning marine species living in sympatry develop and maintain species boundaries through the divergence of gamete recognition proteins and the associated evolution of gamete incompatibility.     

Sea urchin bindin divergence predicts gamete compatibility

Studies on the evolution of reproductive proteins have shown that they tend to evolve more rapidly than other proteins, frequently under positive selection. Progress on understanding the implications of these patterns is possible for marine invertebrates, where molecular evolution can be linked to gamete compatibility. In this study, we surveyed data from the literature from five genera of sea urchins for which there was information on gamete compatibility, divergence of the sperm-egg recognition protein bindin, and mitochondrial divergence. We draw three conclusions: (1) bindin divergence at nonsynonymous sites predicts gamete compatibility, whereas (2) bindin divergence at synonymous sites and mitochondrial DNA divergence do not, and (3) as few as 10 amino acid changes in bindin can lead to complete gamete incompatibility between species. Using mitochondrial divergence as a proxy for time, we find that complete gamete incompatibility can evolve in approximately one and a half million years, whereas sister species can maintain complete gamete compatibility for as long as five million years.     

Multiple gene genealogies reveal asymmetrical hybridization and introgression among strongylocentrotid sea urchins

The evolution of incompatibilities between eggs and sperm is thought to play important roles in establishing and maintaining reproductive isolation among species of broadcast-spawning marine invertebrates. However, the effectiveness of gametic isolation in initiating the speciation process and/or in limiting the introgression of genes among species at later stages of divergence remains largely unknown. In the present study, we collected DNA sequence data from five loci in four species of Strongylocentrotus sea urchins ( S. droebachiensis , S. pallidus , S. purpuratus , and S. franciscanus ) to test whether the susceptibility of S. droebachiensis eggs to fertilization by heterospecific sperm results in gene flow between species. Despite the potential for introgression, a small but statistically significant signal of introgression was observed only between the youngest pair of sister taxa ( S. pallidus and S. droebachiensis ) that was strongly asymmetrical (from the former into the latter). No significant gene flow was observed for either S. purpuratus or S. franciscanus despite the ability of their sperm to readily fertilize the eggs of S. droebachiensis . Our results demonstrate that asymmetrical gamete compatibilities in strongylocentrotids can give rise to asymmetrical patterns of introgression but suggest that gamete traits alone cannot be responsible for maintaining species integrities. The genetic boundaries between strongylocentrotid urchin species in the northeast Pacific appear to be related to postzygotic isolating mechanisms that scale with divergence times and not intrinsic gametic incompatibilities per se .     

HYBRIDIZATION OF SYMPATRIC PATIRIELLA SPECIES (ECHINODERMATA: ASTEROIDEA) IN NEW SOUTH WALES

Three species of the asteroid genus Patiriella occur sympatrically in New South Wales and the possibility for hybridization among them was examined through a series of cross-fertilization experiments. Patiriella calcar and P. gunnii are morphologically distinct as adults but indistinguishable as larvae. Patiriella exigua is morphologically distinct in both its adult and larval morphologies. The gametes of P. calcar and P. gunnii were reciprocally compatible: laboratory crosses between these species produced viable hybrid juveniles. In crosses between female P. calcar and male P. gunnii, most of the juveniles metamorphosed with an arm number intermediate between that of the parents, whereas crosses between female P. gunnii and male P. calcar produced juveniles with an arm number more similar to the maternal phenotype. Heterospecific crosses with P. exigua resulted in low fertilization rates, and viable hybrids were not produced. This species appears capable of self-fertilization. Because hybrids between P. calcar and P. gunnii were viable, neither gametic incompatibility nor hybrid inviability appears to ensure reproductive isolation between these species. Ecological or habitat segregation and temporal separation in breeding may isolate these species in the field. The results demonstrate that if gamete surface recognition molecules are involved in fertilization of P. calcar and P. gunnii, then they are not strongly species specific, at least at the sperm concentrations used in this study. Reproductive isolation between these species has evolved despite their gametic compatibility. In contrast, P. exigua is isolated from its congeners because of gametic incompatibility and several features characteristic of its reproduction and development. The implications of these findings for reproductive isolation and speciation of Patiriella and for the evolution of reproductive isolation in free-spawning marine organisms are discussed.     

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat‐related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model‐based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.     

Cross-incompatibility between two sympatric polyploid Solanum species

Summary Preliminary results from a large number of reciprocal crosses between the closely related sympatric species S. gourlayi Hawkes (2n=4x=48) and S. oplocense Hawkes (2n=6x=72) indicated that they are difficult to hybridize. Pollen-pistil incompatibility barriers were detected via fluorescent microscopy. The cross incompatibility reaction occurred at three sites in 6x×4x crosses; on the stigma, in the first one-third of the style, and in the first two-thirds of the style. In the reciprocal 4x×6x crosses the incompatibility reaction invariably occurred in the ovary. Backcrosses of interspecific pentaploid hybrids (that were occasionally formed) to both parental populations were fully compatible, partially compatible, and fully incompatible with three sites of cross-incompatibility reaction similar to those observed in 6x×4x crosses, respectively. Both polyploid species were found to be selfcompatible, whereas their F1 hybrids were found to be self-incompatible. An hypothesis based on interactions of dominant cross-incompatibility (CI) genes in pistils and dominant specific complementary genes in pollen grains is postulated to explain these observations. The cross-incompatibility system that appears to be operating in nature between 4x S. gourlayi and 6x S. oplocense provides a way for gene exchange between sympatric populations without threatening the identity of either species.     

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.     

Isolation‐driven divergence: speciation in a widespread North American songbird (Aves: Certhiidae)

Lineage, or true ‘species’, trees may differ from gene trees because of stochastic processes in molecular evolution leading to gene‐tree heterogeneity. Problems with inferring species trees because of excessive incomplete lineage sorting may be exacerbated in lineages with rapid diversification or recent divergences necessitating the use of multiple loci and individuals. Many recent multilocus studies that investigate divergence times identify lineage splitting to be more recent than single‐locus studies, forcing the revision of biogeographic scenarios driving divergence. Here, we use 21 nuclear loci from regional populations to re‐evaluate hypotheses identified in an mtDNA phylogeographic study of the Brown Creeper (Certhia americana), as well as identify processes driving divergence. Nuclear phylogeographic analyses identified hierarchical genetic structure, supporting a basal split at approximately 32°N latitude, splitting northern and southern populations, with mixed patterns of genealogical concordance and discordance between data sets within the major lineages. Coalescent‐based analyses identify isolation, with little to no gene flow, as the primary driver of divergence between lineages. Recent isolation appears to have caused genetic bottlenecks in populations in the Sierra Madre Oriental and coastal mountain ranges of California, which may be targets for conservation concerns.     

Geographic variation and positive selection on M7 lysin, an acrosomal sperm protein in mussels (Mytilus spp.)

Successful fertilization in free-spawning marine organisms depends on the interactions between genes expressed on the surfaces of eggs and sperm. Positive selection frequently characterizes the molecular evolution of such genes, raising the possibility that some common deterministic process drives the evolution of gamete recognition genes and may even be important for understanding the evolution of prezygotic isolation and speciation in the marine realm. One hypothesis is that gamete recognition genes are subject to selection for prezygotic isolation, namely, reinforcement. In a previous study, positive selection on the gene coding for the acrosomal sperm protein M7 lysin was demonstrated among allopatric populations of mussels in the Mytilus edulis species group (M. edulis, Mytilus galloprovincialis, and Mytilus trossulus). Here, we expand sampling to include M7 lysin haplotypes from populations where mussel species are sympatric and hybridize to determine whether there is a pattern of reproductive character displacement (RCD), which would be consistent with reinforcement driving selection on this gene. We do not detect a strong pattern of RCD; neither are there unique haplotypes in sympatry nor is there consistently greater population structure in comparisons involving sympatric populations. One distinct group of haplotypes, however, is strongly affected by natural selection, and this group of haplotypes is found within M. galloprovincialis populations throughout the Northern Hemisphere concurrent with haplotypes common to M. galloprovincialis and M. edulis. We suggest that balancing selection, perhaps resulting from sexual conflicts between sperm and eggs, maintains old allelic diversity within M. galloprovincialis.     

Asymmetric conspecific sperm precedence in relation to spawning times in the Montastraea annularis species complex (Cnidaria: Scleractinia)

In broadcast spawners, prezygotic reproductive isolation depends on differences in the spatial and temporal patterns of gamete release and gametic incompatibility. Typically, gametic incompatibility is measured in no‐choice crosses, but conspecific sperm precedence (CSP) can prevent hybridization in gametes that are compatible in the absence of sperm competition. Broadcast spawning corals in the Montastraea annularis species complex spawn annually on the same few evenings. Montastraea franksi spawns an average of 110 min before M. annularis, with a minimum gap of approximately 40 min. Gametes are compatible in no‐choice heterospecific assays, but it is unknown whether eggs exhibit choice when in competition. Hybridization depends on either M. franksi eggs remaining unfertilized and in proximity to M. annularis when the latter species spawns or M. franksi sperm remaining in sufficient viable concentrations when M. annularis spawns. We found that the eggs of the early spawning M. franksi demonstrate strong CSP, whereas CSP appears to be lacking for M. annularis eggs. This study provides evidence of diverging gamete affinities between these recently separated species and suggests for the first time that selection may favour CSP in earlier spawning species when conspecific sperm is diluted and aged and is otherwise at a numeric and viability disadvantage with heterospecific sperm.     

Genetic studies of self incompatibility in the garden chrysanthemum,Chrysanthemum morifolium ramat

Summary Self incompatibility was investigated in the hexaploid garden chrysanthemum, a member of Compositae. Nine sibling clones selected from a highly compatible cross were all self incompatible. 14.8% of the crosses between these sibs in diallel were compatible, but one sib, 67-111-42, accounted for 10 of the 12 compatible crosses. 67-111-42 was also more compatible than the remaining 8 sibs in crosses to other closely related plants. Crosses of the 9 sibs to 12 unrelated tester clones indicated that none were male or female sterile. Inbreeding via pseudocompatibility was successful in increasing homozygosity at the S loci. The percentage of compatible crosses obtained in 3 sib diallels of I ₂ clones from crosses of 67-111-42I ₁ plants approached that of the original 9 × 9 diallel, but no one individual accounted for most of the compatible crosses. It was possible to separate the 9 sibs into 9 incompatibility patterns from the pollinations made in this study. The evidence suggests that the self-incompatibility reaction in the garden chrysanthemum is sporophytic and involves more than 1 locus.Scientific Journal Series Paper Number 7882 of the Minnesota Agricultural Experiment Station.Lyndon W. Drewlow was a National Science Foundation Trainee.     

MORE THAN BINDIN DIVERGENCE: REPRODUCTIVE ISOLATION BETWEEN SYMPATRIC SUBSPECIES OF A SEA URCHIN BY ASYNCHRONOUS SPAWNING

The evolution of reproductive barriers is crucial to the process of speciation. In the Echinoidea, studies have focused on divergence in the gamete recognition protein, bindin, as the primary isolating mechanism among species. As such, the capacity of alternate mechanisms to be effective reproductive barriers and the phylogenetic context in which they arise is unclear. Here, we examine the evolutionary histories and factors limiting gene exchange between two subspecies of Heliocidaris erythrogramma that occur sympatrically in Western Australia. We found low, but significant differentiation between the subspecies in two mitochondrial genes. Further, coalescent analyses suggest that they diverged in isolation on the east and west coasts of Australia, with a subsequent range expansion of H. e. erythrogramma into Western Australia. Differentiation in bindin was minimal, indicating gamete incompatibility is an unlikely reproductive barrier. We did, however, detect strong asynchrony in spawning seasons; H. e. erythrogramma spawned over summer whereas H. e. armigera spawned in autumn. Taken together, we provide compelling evidence for a recent divergence of these subspecies and their reproductive isolation without gamete incompatibility. Western Australian H. erythrogramma may therefore present an intriguing case of incipient speciation, which depends on long‐term persistence of the factors underlying this spawning asynchrony.     

Adaptive evolution of sperm bindin tracks egg incompatibility in neotropical sea urchins of the genus Echinometra

Bindin is a gamete recognition protein known to control species-specificsperm-egg adhesion and membrane fusion in sea urchins. Previousanalyses have shown that diversifying selection on bindin aminoacid sequence is found when gametically incompatible speciesare compared, but not when species are compatible. The presentstudy analyzes bindin polymorphism and divergence in the threeclosely related species of Echinometra in Central America: E.lucunter and E. viridis from the Caribbean, and E. vanbruntifrom the eastern Pacific. The eggs of E. lucunter have evolveda strong block to fertilization by sperm of its neotropicalcongeners, whereas those of the other two species have not.As in the Indo-West Pacific (IWP) Echinometra, the neotropicalspecies show high intraspecific bindin polymorphism in the samegene regions as in the IWP species. Maximum likelihood analysisshows that many of the polymorphic codon sites are under mildpositive selection. Of the fixed amino acid replacements, mosthave accumulated along the bindin lineage of E. lucunter. Weanalyzed the data with maximum likelihood models of variationin positive selection across lineages and codon sites, and withmodels that consider sites and lineages simultaneously. Ourresults show that positive selection is concentrated along theE. lucunter bindin lineage, and that codon sites with aminoacid replacements fixed in this species show by far the highestsignal of positive selection. Lineage-specific positive selectionparalleling egg incompatibility provides support that adaptiveevolution of sperm proteins acts to maintain recognition ofbindin by changing egg receptors. Because both egg incompatibilityand bindin divergence are greater between allopatric speciesthan between sympatric species, the hypothesis of selectionagainst hybridization (reinforcement) cannot explain why adaptiveevolution has been confined to a single lineage in the AmericanEchinometra. Instead, processes acting to varying degrees withinspecies (e.g., sperm competition, sexual selection, and sexualconflict) are more promising explanations for lineage-specificpositive selection on bindin.     

Marine speciation on a small planet

The scale of population structure in many marine species is on the order of thousands to tens of thousands of kilometers. How does speciation take place in oceans that are only about this same size? Recent results suggest an important role for transient isolation, gamete ecology and molecular evolution at gamete recognition loci. These factors have long been appreciated by plant biologists, and are likely to be a fruitful area of research for marine biologists as well.     

Divergence in coloration and ecological speciation in the Anolis marmoratus species complex

Adaptive divergence in coloration is expected to produce reproductive isolation in species that use colourful signals in mate choice and species recognition. Indeed, many adaptive radiations are characterized by differentiation in colourful signals, suggesting that divergent selection acting on coloration may be an important component of speciation. Populations in the Anolis marmoratus species complex from the Caribbean island of Guadeloupe display striking divergence in the colour and pattern of adult males that occurs over small geographic distances, suggesting strong divergent selection. Here we test the hypothesis that divergence in coloration results in reduced gene flow among populations. We quantify variation in adult male coloration across a habitat gradient between mesic and xeric habitats, use a multilocus coalescent approach to infer historical demographic parameters of divergence, and examine gene flow and population structure using microsatellite variation. We find that colour variation evolved without geographic isolation and in the face of gene flow, consistent with strong divergent selection and that both ecological and sexual selection are implicated. However, we find no significant differentiation at microsatellite loci across populations, suggesting little reproductive isolation and high levels of contemporary gene exchange. Strong divergent selection on loci affecting coloration probably maintains clinal phenotypic variation despite high gene flow at neutral loci, supporting the notion of a porous genome in which adaptive portions of the genome remain fixed whereas neutral portions are homogenized by gene flow and recombination. We discuss the impact of these findings for studies of colour evolution and ecological speciation.     

Introgression and selection shaped the evolutionary history of sympatric sister‐species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister‐species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site‐associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister‐species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174–0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.     

Evolutionary consequences of introgression at M7 lysin, a gamete recognition locus, following secondary contact between blue mussel species

Hybridizing populations of blue mussels, Mytilus edulis and Mytilus trossulus, in Cobscook Bay (eastern Maine) have been used by our laboratory to study the evolution of gamete incompatibility and molecular evolution of the vitelline coat lysin proteins expressed in sperm. The M7 lysin locus has been the most studied of the three lysins, but evidence for positive selection necessary to help confirm its role in gamete recognition in western Atlantic hybrid zones is contradictory. We developed an alternative test, based on rates of introgression at M7 lysin. Contrary to expectations, introgression at this locus is much higher (instead of much lower) than is introgression at neutral markers. In this article, we present simulations, constructed using synthetic populations of combinations of admixed genotypes, representing various hybrid offspring categories. Simulations produced variation in introgression across loci, but did not generate the massive introgression at M7 lysin observed in natural populations in Cobscook Bay. We consider these results in the context of selection operating on gamete recognition loci, both within and between species, during the third stage of allopatric speciation in Mytilus.     

Genetic divergence with ongoing gene flow is maintained by the use of different hosts in phytophagous ladybird beetles genus Henosepilachna

Adaptation to different environments can promote population divergence via natural selection even in the presence of gene flow – a phenomenon that typically occurs during ecological speciation. To elucidate how natural selection promotes and maintains population divergence during speciation, we investigated the population genetic structure, degree of gene flow and heterogeneous genomic divergence in three closely related Japanese phytophagous ladybird beetles: Henosepilachna pustulosa, H. niponica and H. yasutomii. These species act as a generalist, a wild thistle (Cirsium spp.) specialist and a blue cohosh (Caulophyllum robustum) specialist, respectively, and their ranges differ accordingly. The two specialist species widely co‐occur but are reproductively isolated solely due to their high specialization to a particular host plant. Genomewide amplified fragment‐length polymorphism (AFLP) markers and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences demonstrated obvious genomewide divergence associated with both geographic distance and ecological divergence. However, a hybridization assessment for both AFLP loci and the mitochondrial sequences revealed a certain degree of unidirectional gene flow between the two sympatric specialist species. Principal coordinates analysis (PCoA) based on all of the variable AFLP loci demonstrated that there are genetic similarities between populations from adjacent localities irrespective of the species (i.e. host range). However, a further comparative genome scan identified a few fractions of loci representing approximately 1% of all loci as different host‐associated outliers. These results suggest that these three species had a complex origin, which could be obscured by current gene flow, and that ecological divergence can be maintained with only a small fraction of the genome is related to different host use even when there is a certain degree of gene flow between sympatric species pairs.     

Reproductive protein evolution in two cryptic species of marine chordate

Background  

Reproductive character displacement (RCD) is a common and taxonomically widespread pattern. In marine broadcast spawning organisms, behavioral and mechanical isolation are absent and prezygotic barriers between species often operate only during the fertilization process. Such barriers are usually a consequence of differences in the way in which sperm and egg proteins interact, so RCD can be manifest as faster evolution of these proteins between species in sympatry than allopatry. Rapid evolution of these proteins often appears to be a consequence of positive (directional) selection. Here, we identify a set of candidate gamete recognition proteins (GRPs) in the ascidian Ciona intestinalis and showed that these GRPs evolve more rapidly than control proteins (those not involved in gamete recognition). Choosing a subset of these gamete recognition proteins that show evidence of positive selection (CIPRO37.40.1, CIPRO60.5.1, CIPRO100.7.1), we then directly test the RCD hypothesis by comparing divergence (omega) and polymorphism (McDonald-Kreitman, Tajima's D, Fu and Li's D and F, Fay and Wu's H) statistics in sympatric and allopatric populations of two distinct forms of C. intestinalis (Types A and B) between which there are strong post-zygotic barriers.