Introduction
Speciation in fertilization is critical in the beginning of development. While mechanisms of fertilization have diverged from invertebrates to humans, several strictly defined steps have been governed by specific molecules from eggs and sperm for the recognition of one another (Vacquier, Reference Vacquier1998). These steps include sperm activation and/or attraction to the egg, as well as sperm acrosome reaction upon egg contact. Understanding the mechanisms of these molecular interactions is important, as any disruptions of these processes could cause infertility and destroy species barriers at the molecular level.
In the sea star, Asterias amurensis, the asteroidal sperm-activating peptide (asterosap) has been purified and identified from the extracellular jelly coat of the egg. Asterosap is composed of 34 amino acids and has an intramolecular disulphide bond that is essential for its physiological activity (Nishigaki et al., Reference Nishigaki, Chiba, Miki and Hoshi1996). The asterosap precursor mRNA encodes 10 asterosap molecules that are tandemly connected by spacer sequences (Matsumoto et al., Reference Matsumoto, Briones, Nishigaki and Hoshi1999). Asterosap binds to the extracellular domain of its receptor, a membrane-bound guanylate cyclase, in the sperm flagellum to activate its intracellular catalytic domain for the production of cGMP (Nishigaki et al., Reference Nishigaki, Chiba and Hoshi2000; Matsumoto et al., Reference Matsumoto, Solzin, Helbig, Hagen, Ueno, Kawase, Maruyama, Ogiso, Godde, Minakata, Kaupp, Hoshi and Weyand2003). The cooperation of three egg jelly coat components, ARIS (acrosome reaction-inducing substance), Co-ARIS and asterosap, induces the acrosome reaction (Hoshi et al., Reference Hoshi, Nishigaki, Ushiyama, Okinaga, Chiba and Matsumoto1994). Asterosap is essential for the induction of the acrosome reaction by elevating the internal pH of sperm (Kawase et al., Reference Kawase, Minakata, Hoshi and Matsumoto2005). In addition, asterosap induces hyperactivation-like flagellar movement in sperm that are firmly attached to a solid surface (Shiba et al., Reference Shiba, Tagata, Ohmuro, Mogami, Matsumoto, Hoshi and Baba2006), potentially allowing an effective acrosome reaction on the surface of the jelly coat. Furthermore, asterosap also has sperm-attracting activity (Bohmer et al., Reference Bohmer, Van, Weyand, Hagen, Beyermann, Matsumoto, Hoshi, Hildebrand and Kaupp2005). The specificity of asterosap activity is conserved within the superorder level (Nakachi et al., Reference Nakachi, Moriyama, Hoshi and Matsumoto2006).
In sea urchin, resact (SAP-IIA) from Arbacia punctulata is the only peptide known to have sperm-attracting activity among more than 70 sperm-activating peptides identified from the egg jelly coat (Ward et al., Reference Ward, Brokaw, Garbers and Vacquier1985). Resact binds directly to a membrane-bound guanylate cyclase receptor in sperm and activates it (Singh et al., Reference Singh, Lowe, Thorpe, Rodriguez, Kuang, Dangott, Chinkers, Goeddel and Garbers1988; Kaupp et al., Reference Kaupp, Solzin, Hildebrand, Brown, Helbig, Hagen, Beyermann, Pampaloni and Weyand2003).
In order to investigate the molecular details of the interaction between asterosap and its receptor, we cloned the cDNA encoding these proteins from five closely related sea star species. We report that although the amino acid sequences are highly conserved over these species, small differences in the sequences can explain their speciation. The functional mechanism of the asterosap receptor is also discussed by analogy to another membrane-bound guanylate cyclase, the atrial natriuretic peptide receptor.
Materials and Methods
Animals and gametes
Sea stars were collected locally during their annual breeding seasons as described previously (Nakachi et al., Reference Nakachi, Moriyama, Hoshi and Matsumoto2006). The sea stars used in this study were Asterias amurensis (Aa), A. rubens (Ar), A. forbesi (Af), Aphelasterias japonica (Apj) and Distolasterias nipon (Dn).
Cloning of genes
Total RNA was prepared from fresh testes or ovaries by the acid guanidium–phenol–chloroform (AGPC) method (Sambrook et al., Reference Sambrook, Fritsch and Maniatis1989). First-strand cDNA was synthesized using reverse transcriptase with either a random 9-mer using Rever Tra Ace (TOYOBO) or the ImProm-2 reverse transcription system (Promega) and was subsequently used as a template for PCR. DNA fragments of asterosap and spacer sequences were amplified by PCR from ovary cDNA using primers TN1 5′-ATCGCTCGAGTAAGACAA-3′ and TN2 5′-TCGTTTGCTCTGACATTG-3′. Extracellular domain sequences of the asterosap receptors were amplified by PCR from testes cDNA using primers hani1 5′-AGGTGTTTGATGCTGTCTGTCG-3′, hani2 5′-GCGTAGTCTCCGTTATTGAGACC-3′ MN1′ 5′-AAGGAGACAACAAGGATTTAACGTG-3′ and MN2 5′-GACCGGCATATCTAGAGGTGG-3′. Guanylate cyclase domain sequences of the asterosap receptors were also amplified from testes cDNA using primers P2 5′-GACATCGTCGGGTTTACC-3′ and P6 5′-AGTGTTGACGGTGTCACC-3′. PCR products were separated by size via agarose gel electrophoresis, extracted and cloned into the pGEM-T vector (Promega) and sequenced with the M13 forward and reverse primers using BigDye chemistry on an ABI PRISM 3100 DNA sequencer (Applied Biosystems). We deposited the sequences into GenBank and the accession numbers are listed in the Supplemental materials section.
Sequence analysis
Multiple alignments of nucleotide and amino acid sequences were done using Clustal X (Thompson et al., Reference Thompson, Gibson, Plewniak, Jeanmougin and Higgins1997). Phylogenetic reconstruction of nucleotide sequences was performed using PHYLIP version 3.66 (Felsenstein, Reference Felsenstein2005) and the phylogenetic trees were viewed using njplot (Perriere & Gouy, Reference Perriere and Gouy1996). Secondary structures and β–turn structures of the amino acid sequences were predicted using PORTER (Pollastri & McLysaght, Reference Pollastri and McLysaght2005) and BetaTPred2 (Kaur & Raghava, Reference Kaur and Raghava2003), respectively. Hydropathy index and pI of amino acid sequences were calculated by sequence analysis (Gilbert, Reference Gilbert2006) using Kyte & Doolittle parameters (Kyte & Doolittle, Reference Kyte and Doolittle1982) in a running average of five amino acids. The coiled-coil structures in the amino acid sequences were predicted by COILS (Lupas et al., Reference Lupas, Van Dyke and Stock1991). Accession numbers for the guanylate cyclase genes were X12874 from sea urchin (Arbacia punctulata; Ap) and AB070354 from Aa. Accession numbers for the mitochondrial cytochrome c oxidase 1 sequences other than ours were AF498647 (Af), AF031028 (Ap) and AF525453 (Hemicentrotus pulcherrimus, Hp, sea urchin). Sequence alignments referring to their 3D structures were performed by the 3D-PSSM (Kelley et al., Reference Kelley, MacCallum and Sternberg2000) and Swiss model (Schwede et al., Reference Schwede, Kopp, Guex and Peitsch2003) servers. The hormone-bound, atrial natriuretic peptide receptor extracellular domain (Ogawa et al., Reference Ogawa, Qiu, Ogata and Misono2004) was accessed from the Protein Data Bank (ID code 1t34).
Results
In order to investigate the molecular details of the interaction between asterosap and its receptor, we cloned the cDNAs encoding these proteins from five closely related sea star species.
Sequence variation among asterosap
Sea star asterosap sequences exist as 10 consecutive repeats interspersed with short spacers (Matsumoto et al., Reference Matsumoto, Briones, Nishigaki and Hoshi1999). These 10 asterosap sequences from the same genome are not identical but vary among them. The cDNAs of randomly chosen seven asterosap isoforms were each cloned including their spacer regions from five sea star species: Asterias amurensis (Aa), A. rubens (Ar), A. forbesi (Af) and Aphelasterias japonica (Apj), from the Asteriinae subfamily and Distolasterias nipon (Dn) from the Coscinasteriinae subfamily. All these five species belong to the family Asteriidae.
The deduced amino acid sequences of the asterosaps were highly conserved among the species with more variation found in the N-terminal region (Fig. 1). Of the 35 cloned sequences, the asterosaps had 10 completely conserved amino acid sites out of 34, covering 29% of the peptide length. While all asterosaps were of the same length, the spacer regions varied ranging from 12 to 22 amino acids. The spacer regions had four completely conserved amino acids sites out of 12 (excluding gap sites), corresponding to 33% identity.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160705004222-01028-mediumThumb-S0967199408004759_fig1g.jpg?pub-status=live)
Figure 1 Comparison of asterosap and spacer amino acid sequences from five sea star species. Plus (+) and minus (−) position numbers at the top indicate positions within the spacers from the N- and C-termini, respectively. A dot (·) shows the residue identical to Aa01; a minus sign (−) indicates a gap. The slash (/) in the ruler shows the potential cleavage site. Sequences used for PCR primers are underlined by dotted lines in Aa01 and primer names are indicated in square brackets ([]) above with arrows indicating the amplification directions. Two conserved Cys residues are marked by number signs (#) above Aa01. Positively and negatively charged conserved residues in asterosaps are shaded with light and dark grey, respectively. The ‘Common’ row shows the symbols based on the sequence similarity for each position after CLUSTAL X: an asterisk (*) for a completely conserved position; a colon (:) for the position conserved among the ‘strong’ amino acid groups; a period (.) for the position conserved among the ‘weaker’ amino acid groups. The ‘Motif’ row shows the predicted secondary structures for Aa04: a tilde (∼) for the coil structure; and an equal (=) for the helical structure. The ‘Phobic’ row shows Kyte & Doolittle hydropathy index from a running average of five amino acids: a plus (+) for hydrophobicity; a vertical bar (|) for hydrophilicity.
In A. amurensis (Aa), asterosap has an intramolecular disulphide bond between Cys8 and Cys32, which is known to be essential for its activity (Nishigaki et al., Reference Nishigaki, Chiba, Miki and Hoshi1996). These Cys residues were among the completely conserved amino acids from all asterosap sequences examined. As shown in Fig. 1, some amino acid substitutions were unique in Asterias, in Apj or in Dn. All the spacer regions had paired basic residues, Lys–Arg or Lys–Lys, at both ends, which are commonly found at the cleavage sites of pro-hormone precursors (Hook et al., Reference Hook, Yasothornsrikul, Greenbaum, Medzihradszky, Troutner, Toneff, Bundey, Logrinova, Reinheckel, Peters and Bogyo2004). The results suggest that the asterosap precursor undergoes endoproteolytic processing at either the C- or N-termini of paired basic residues followed by the removal of the remaining basic residues as shown in other pro-hormone processing pathways (Fricker, Reference Fricker1988; Steiner, Reference Steiner1998; Hook et al., Reference Hook, Yasothornsrikul, Greenbaum, Medzihradszky, Troutner, Toneff, Bundey, Logrinova, Reinheckel, Peters and Bogyo2004). The pI represents the favourable environmental pH range of biologically active molecules such as asterosap. While the distribution of pIs for asterosap was between 4.5 and 10.0 in Aa, Af, Ar and Apj, pIs for asterosaps from Dn range from 3.8 to 6.3.
Every asterosap was predicted to have two amphipathic helices between two conserved Cys residues (Fig. 1). Positive charges (Arg11 and 13) in the first helix and negative charges (Asp20 and Glu21) in the second helix were highly conserved. Then, we assumed that these helical structures are important for asterosap binding to its receptor. Hydropathy index of the asterosaps indicated the overall hydrophilicity, explaining its diffusive characteristic in seawater. Conserved hydrophobic residues in each helix may also be involved in either binding to its receptor or peptide folding. Spacer regions had turn structures and were entirely hydrophilic.
Asterosap receptor
A. amurensis asterosap receptor has four domains: extracellular, transmembrane, kinase-like and guanylate cyclase (Matsumoto et al., Reference Matsumoto, Solzin, Helbig, Hagen, Ueno, Kawase, Maruyama, Ogiso, Godde, Minakata, Kaupp, Hoshi and Weyand2003). We cloned the cDNAs of the extracellular and guanylate cyclase domains of the asterosap receptor from four additional sea star species and compared these sequences. The deduced amino acid sequences (Fig. 2) shared high similarities in these domains: 86% in the extracellular domain and 94% in the guanylate cyclase domain. Most amino acid differences were found between Asterias (Aa, Ar and Af) compared with Apj and Dn: more than 10% of the extracellular domain sequences differed. In the extracellular domain, all four Cys residues and three potential glycosylation sites were conserved among the five species. In the guanylate cyclase domain, all putative catalytic sites (Liu et al., Reference Liu, Force and Bloch1997) were conserved (Fig. 2).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160705004222-48914-mediumThumb-S0967199408004759_fig2g.jpg?pub-status=live)
Figure 2 Comparison of asterosap receptor amino acid sequences from five sea star species. A dot (·) indicates the amino acid identical to Aa at the corresponding position. Mature proteins start from the residue +1 after the N-terminal signal peptides. Each domain starts at |→ and ends at ←|. Sequences used for PCR primers are underlined with dotted lines in Aa and primer names are given in square brackets ([]) with arrows indicating the amplification directions. An asterisk (*), a colon (:) and a period (.) at the bottom show the levels of sequence similarities for each position (see the Fig. 1 legend). Extracellular Cys residues are marked by a number signs (#). Potential N-glycosylation sites (Asn–X–Ser/Thr) are marked by a plus (+) in parentheses. Parts of the Aa sequence subjected to the homology modelling are underlined. A potential ligand-binding region is double underlined with positively and negatively charged residues shaded with light and dark grey, respectively. Potential phosphorylation sites in the region equivalent to the rat atrial natriuretic peptide receptor are shaded with light grey in Aa (from positions 500 to 517). A potential coiled-coil region is underlined by short- and long-dotted lines. Potential GTP binding sites in the guanylate cyclase domain are marked by ‘g’. Except for Aa, only extracellular and guanylate cyclase domains were cloned and sequenced.
Molecular evolution of asterosap and its receptor
In order to compare the evolutionary patterns between asterosap and its receptor sequences, we performed phylogenetic reconstruction at the nucleotide level (Fig. 3). The phylogenetic tree reconstructed using the 35 asterosap sequences showed three major groups, each corresponding to the three subfamilies: Dn, Apj and Asterias, although the bootstrap analysis did not support these three clusters (Fig. 3a). Interestingly, clustering among Asterias asterosap sequences was not consistent with species clustering. For example, the clusters from a to d in Fig. 3a included more than two Asterias species. Moreover, identical asterosap sequences were found between Af02 and Aa13 and between Apj07 and Af05. More phylogenetics analysis may be required to establish taxonomical relationships among these sea star species, especially the phylogenetic location of Apj. The occurrence of close relationships among Asterias species found in phylogenies agrees with our previous results on the gamete-compatibility grouping in sea stars (Nakachi et al., Reference Nakachi, Moriyama, Hoshi and Matsumoto2006).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160705004222-11462-mediumThumb-S0967199408004759_fig3g.jpg?pub-status=live)
Figure 3 Phylogenetic relationships reconstructed based on the nucleotide sequences of asterosap (a) and its receptor (b and c). The receptor phylogenies were reconstructed using the extracellular domain (B, 1299 nucleotides) and the guanylate cyclase domain (C, 326 nucleotides) separately. Neighbour-joining phylogenetic method (Saitou and Nei, 1987) was used on nucleotide distances estimated by F84 model (Kishino and Hasegawa, 1989). Bootstrap analysis was done with 1000 pseudoreplicates and supporting values greater than 80% are shown on each consensus phylogeny. For asterosap receptor analyses, resact receptor sequences from sea urchin were added as outgroups: Ap, Arbacia punctulata; Hp, Hemicentrotus pulcherrimus.
Discussion
Interaction between asterosap and its receptor
The structure of the asterosap receptor and the binding of the asterosap to its receptor were modelled based on the structure of the rat atrial natriuretic peptide and its receptor (ANP and NPR-A; PDB entry 1t34). The sequence identity between the asterosap receptor and NPR-A sequences was 27%, which is considered sufficient for structural modelling. NPR-A is known to form a dimer (Wilson & Chinkers, Reference Wilson and Chinkers1995). Our native PAGE result showed that the Aa asterosap receptor in sperm appears to be present as a dimer and oligomer (data not shown). In addition, the Aa asterosap receptor has a potential coiled-coil region in the area similar to that of NPR-A, which is important for dimerization (Wilson & Chinkers, Reference Wilson and Chinkers1995). A sequence alignment between the Aa asterosap receptor and NPR-A, taking their 3-D structures into consideration, indicated that they share key amino acid residues, Trp71 and Asp59 in the asterosap receptor and Trp74 and Asp62 in NPR-A, as well as Arg 13 in asterosap and Arg14 in ANP (Fig. 4). These amino acids have been strictly conserved in all examined sea star sequences.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20160705004222-49906-mediumThumb-S0967199408004759_fig4g.jpg?pub-status=live)
Figure 4 A schematic model of hormone and receptor binding. The rat atrial natriuretic peptide (ANP) and its receptor (NPR-A) are shown on the left. Asterosap and its receptor are shown on the right. Two identical receptor subunits are connected by ionic interactions at the top. A large hydrophobic tryptophan (W) residue from each subunit acts as a spring to secure the hinge. Interactions between the negatively charged aspartic acid (D) and the positively charged arginine (R) play important roles in the binding between ligands and the receptors. The molecular size of asterosap peptide (3700 Da, 34 amino acids) is larger than that of ANP (3100 Da, 28 amino acids).
Resact is a sperm-activating peptide from Arbacia punctulata (Ap) and is composed of 14 residues with an intramolecular disulphide bond. The primary structure of the resact receptor was also compared with the asterosap receptor to predict their functional similarity. Their extracellular domains have high sequence similarity. While the resact receptor lacks the large C-terminal portion of the guanylate cyclase domain including almost all putative key catalytic residues, the asterosap receptor and resact receptor may have similar modes of ligand binding but different catalytic regulatory mechanisms.
Survival strategies of sea stars
The acrosome reaction is an important fertilization process in sea stars, which takes place upon the arrival of sperm onto the egg jelly surface. ARIS, a huge proteoglycan, is the major component of egg jelly in sea stars. A sulphate moiety in ARIS plays an important role in the acrosome reaction. The pH of the egg jelly and the neighbouring seawater can be affected by the content of the sulphate moiety caused by its acidic properties. Each species has several varieties of asterosaps at different pIs. This variation is advantageous for sea stars because a wide distribution of pIs allows with the organism to adjust to any environmental changes. Alternatively, differences in pI can reflect the different sugar components of egg jelly in sea stars. For example, Dn has fewer sulphate moieties in its egg jelly, which results in a different pH in the environment.
Aa, Apj and Dn live in close proximity off the Japanese shore in the North Pacific Ocean. There is, therefore, a certain risk of cross-fertilization between them. The fertilization season of Apj (fall) is different from that of Aa and Dn (early spring). While Aa and Apj are able to induce cross-acrosome reactions in vitro using each other's egg jelly (Nakachi et al., Reference Nakachi, Moriyama, Hoshi and Matsumoto2006), natural cross-fertilization is prohibited by the season. The cross-acrosome reaction (i.e. fertilization) is, therefore, unlikely to occur between Aa, Apj and Dn. These facts are evident in the phylogenic relationship presented in Fig. 3.
We have identified homologous guanylate cyclase signalling systems from five different sea star species, presumably the first attempt in marine invertebrates from the North and South Pacific Ocean as well as the Atlantic Ocean. The cloned asterosap and guanylate cyclase sequences both exhibited speciation. Although the number of clones was limited caused by specimen availability, the analysed codon utilization did not have a strong bias. As we accumulate more clones, we will have the capacity to test evolutionary hypotheses, such as positive selection and co-evolution of asterosap and its receptor.
Acknowledgements
We thank the directors and staff of the following universities and centres for their assistance in collecting sea stars: the Research Center for Marine Biology (Tohoku University, Asamushi, Japan), the International Coastal Research Center (University of Tokyo, Otsuchi, Japan), the Ushimado Marine Laboratory (Okayama University, Japan), the Marine Research Laboratories (University of Tasmania, Australia), the University Marine Biological Station (Millport, Scotland) and the Marine Biological Laboratory (Woods Hole, USA). We thank Ms Mamta Bajaj (H.M. Lab) for her initial phylogenetic analysis.
Supplemental material
GenBank accession numbers
[asterosap]
Aa01 DQ985961
Aa03 DQ985962
Aa04 DQ985963
Aa06 DQ985964
Aa08 DQ985965
Aa11 DQ985966
Aa13 DQ985967
Af02 DQ985968
Af03 DQ985969
Af04 DQ985970
Af05 DQ985971
Af08 DQ985972
Af09 DQ985973
Af11 DQ985974
Ar01 DQ985975
Ar02 DQ985976
Ar03 DQ985977
Ar04 DQ985978
Ar05 DQ985979
Ar06 DQ985980
Ar07 DQ985981
Apj01 DQ985982
Apj02 DQ985983
Apj03 DQ985984
Apj06 DQ985985
Apj07 DQ985986
Apj10 DQ985987
Apj11 DQ985988
Dn03 DQ985989
Dn04 DQ985990
Dn05 DQ985991
Dn06 DQ985992
Dn07 DQ985993
Dn08 DQ985994
Dn09 DQ985995
[asterosap receptor extracellular domain]
ArGC_ec DQ992550
AfGC_ec DQ992551
ApjGC_ec DQ992552
DnGC_ec DQ992553
[asterosap receptor guanylate cyclase domain]
ArGC_gc DQ992554
AfGC_gc DQ992555
ApjGC_gc DQ992556
DnGC_gc DQ992557
[mitochondrial cytochrome c oxidase 1]
Aa DQ992558
Ar DQ992559
Apj DQ992560
Dn DQ992561