Hostname: page-component-745bb68f8f-g4j75 Total loading time: 0 Render date: 2025-02-11T01:43:15.611Z Has data issue: false hasContentIssue false

Conserved sequences of sperm-activating peptide and its receptor throughout evolution, despite speciation in the sea star Asterias amurensis and closely related species

Published online by Cambridge University Press:  01 August 2008

Mia Nakachi
Affiliation:
Department of Biosciences and Informatics, Keio University, Hiyoshi, Kouhoku-ku, Yokohama 223–8522, Japan. Department of Chemistry, Environmental Toxicology and Biotechnology, University of Nebraska – Lincoln, 28 Hamilton Hall, Lincoln, Nebraska 68588–0304, USA.
Motonori Hoshi
Affiliation:
Department of Biosciences and Informatics, Keio University, Hiyoshi, Kouhoku-ku, Yokohama 223–8522, Japan.
Midori Matsumoto
Affiliation:
Department of Biosciences and Informatics, Keio University, Hiyoshi, Kouhoku-ku, Yokohama 223–8522, Japan.
Hideaki Moriyama*
Affiliation:
Department of Chemistry, University of Nebraska – Lincoln, 28 Hamilton Hall, Lincoln, Nebraska 68588–0304, USA. Department of Chemistry, Environmental Toxicology and Biotechnology, University of Nebraska – Lincoln, 28 Hamilton Hall, Lincoln, Nebraska 68588–0304, USA.
*
All correspondence to: Hideaki Moriyama. Department of Chemistry, University of Nebraska – Lincoln, 28 Hamilton Hall, Lincoln, Nebraska 68588–0304, USA. Tel: +1 402 472 5367. Fax: +1 402 472 9402. e-mail: hmoriyama2@unl.edu
Rights & Permissions [Opens in a new window]

Summary

The asteroidal sperm-activating peptides (asterosaps) from the egg jelly bind to their sperm receptor, a membrane-bound guanylate cyclase, on the tail to activate sperm in sea stars. Asterosaps are produced as single peptides and then cleaved into shorter peptides. Sperm activation is followed by the acrosome reaction, which is subfamily specific. In order to investigate the molecular details of the asterosap–receptor interaction, corresponding cDNAs have been cloned, sequenced and analysed from the Asteriinae subfamily including Asterias amurensis, A. rubens, A. forbesi and Aphelasterias japonica, as well as Distolasterias nipon from the Coscinasteriinae subfamily. Averages of 29% and 86% identity were found from the deduced amino acid sequences in asterosap and its receptor extracellular domains, respectively, across all species examined. The phylogenic tree topology for asterosap and its receptor was similar to that of the mitochondrial cytochrome c oxidase subunit I. In spite of a certain homology, the amino acid sequences exhibited speciation. Conservation was found in the asterosap residues involved in disulphide bonding and proteinase-cleaving sites. Conversely, similarities were detected between potential asterosap-binding sites and the structure of the atrial natriuretic peptide receptor. Although the sperm-activating peptide and its receptor share certain common sequences, they may serve as barriers that ensure speciation in the sea star A. amurensis and closely related species.

Type
Research Article
Copyright
Copyright © Cambridge University Press 2008

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.

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

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

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.

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

References

Bohmer, M., Van, Q., Weyand, I., Hagen, V., Beyermann, M., Matsumoto, M., Hoshi, M., Hildebrand, E. & Kaupp, U.B. (2005). Ca2+ spikes in the flagellum control chemotactic behavior of sperm. EMBO J. 24, 2741–52.CrossRefGoogle ScholarPubMed
Felsenstein, J. (2005). PHYLIP (Phylogeny Inference Package) version 3.6. Distributed by the author. Department of Genome Sciences, University of Washington, Seattle.Google Scholar
Fricker, L. D. (1988). Activation and membrane binding of carboxypeptidase E. J. Cell Biochem. 38, 279–89.CrossRefGoogle ScholarPubMed
Gilbert, W. (2006). Sequence Analysis 1.6.0. Informagen.Google Scholar
Hook, V., Yasothornsrikul, S., Greenbaum, D., Medzihradszky, K.F., Troutner, K., Toneff, T., Bundey, R., Logrinova, A., Reinheckel, T., Peters, C. & Bogyo, M. (2004). Cathepsin L and Arg/Lys aminopeptidase: a distinct prohormone processing pathway for the biosynthesis of peptide neurotransmitters and hormones. Biol. Chem. 385, 473–80.CrossRefGoogle ScholarPubMed
Hoshi, M., Nishigaki, T., Ushiyama, A., Okinaga, T., Chiba, K. & Matsumoto, M. (1994). Egg-jelly signal molecules for triggering the acrosome reaction in starfish spermatozoa. Int. J. Dev. Biol. 38, 167–74.Google ScholarPubMed
Kaupp, U.B., Solzin, J., Hildebrand, E., Brown, J.E., Helbig, A., Hagen, V., Beyermann, M., Pampaloni, F. & Weyand, I. (2003). The signal flow and motor response controling chemotaxis of sea urchin sperm. Nat. Cell Biol. 5, 109–17.CrossRefGoogle ScholarPubMed
Kaur, H. & Raghava, G.P. (2003). Prediction of beta-turns in proteins from multiple alignment using neural network. Protein Sci. 12, 627–34.CrossRefGoogle ScholarPubMed
Kawase, O., Minakata, H., Hoshi, M. & Matsumoto, M. (2005). Asterosap-induced elevation in intracellular pH is indispensable for ARIS-induced sustained increase in intracellular Ca2+ and following acrosome reaction in starfish spermatozoa. Zygote 13, 6371.CrossRefGoogle ScholarPubMed
Kelley, L. A., MacCallum, R.M. & Sternberg, M.J. (2000). Enhanced genome annotation using structural profiles in the program 3D-PSSM. J. Mol. Biol. 299, 499520.CrossRefGoogle ScholarPubMed
Kyte, J. & Doolittle, R.F. (1982). A simple method for displaying the hydropathic character of a protein. J. Mol. Biol. 157, 105–32.CrossRefGoogle ScholarPubMed
Liu, H., Force, T. & Bloch, K.D. (1997). Nerve growth factor decreases soluble guanylate cyclase in rat pheochromocytoma PC12 cells. J. Biol. Chem. 272, 6038–43.CrossRefGoogle ScholarPubMed
Lupas, A., Van Dyke, M. & Stock, J. (1991). Predicting coiled coils from protein sequences. Science 252, 1162–4.CrossRefGoogle ScholarPubMed
Matsumoto, M., Briones, A.V., Nishigaki, T. & Hoshi, M. (1999). Sequence analysis of cDNAs encoding precursors of starfish asterosaps. Dev. Genet. 25, 130–6.3.0.CO;2-8>CrossRefGoogle ScholarPubMed
Matsumoto, M., Solzin, J., Helbig, A., Hagen, V., Ueno, S., Kawase, O., Maruyama, Y., Ogiso, M., Godde, M., Minakata, H., Kaupp, U.B., Hoshi, M. & Weyand, I. (2003). A sperm-activating peptide controls a cGMP-signaling pathway in starfish sperm. Dev. Biol. 260, 314–24.CrossRefGoogle ScholarPubMed
Nakachi, M., Moriyama, H., Hoshi, M. & Matsumoto, M. (2006). Acrosome reaction is subfamily specific in sea star fertilization. Dev. Biol. 298, 597604.CrossRefGoogle ScholarPubMed
Nishigaki, T., Chiba, K. & Hoshi, M. (2000). A 130-kDa membrane protein of sperm flagella is the receptor for asterosaps, sperm-activating peptides of starfish Asterias amurensis. Dev. Biol. 219, 154–62.CrossRefGoogle ScholarPubMed
Nishigaki, T., Chiba, K., Miki, W. & Hoshi, M. (1996). Structure and function of asterosaps, sperm-activating peptides from the jelly coat of starfish eggs. Zygote 4, 237–45.CrossRefGoogle ScholarPubMed
Ogawa, H., Qiu, Y., Ogata, C.M. & Misono, K.S. (2004). Crystal structure of hormone-bound atrial natriuretic peptide receptor extracellular domain: rotation mechanism for transmembrane signal transduction. J. Biol. Chem. 279, 28625–31.CrossRefGoogle ScholarPubMed
Perriere, G. & Gouy, M. (1996). WWW-query: an on-line retrieval system for biological sequence banks. Biochimie 78, 364–9.CrossRefGoogle ScholarPubMed
Pollastri, G. & McLysaght, A. (2005). Porter: a new, accurate server for protein secondary structure prediction. Bioinformatics 21, 1719–20.CrossRefGoogle Scholar
Qiu, Y., Ogawa, H., Miyagi, M. & Misono, K.S. (2004). Constitutive activation and uncoupling of the atrial natriuretic peptide receptor by mutations at the dimer interface. Role of the dimer structure in signalling. J. Biol. Chem. 279, 6115–23.CrossRefGoogle ScholarPubMed
Sambrook, J., Fritsch, E.F. & Maniatis, T. (1989). Molecular Cloning: A Laboratory Manual. Plainview, NewYork: Cold Spring Harbor Laboratory Press.Google Scholar
Schwede, T., Kopp, J., Guex, N. & Peitsch, M.C. (2003). SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res. 31, 3381–5.CrossRefGoogle ScholarPubMed
Shiba, K., Tagata, T., Ohmuro, J., Mogami, Y., Matsumoto, M., Hoshi, M. & Baba, S.A. (2006). Peptide-induced hyperactivation-like vigorous flagellar movement in starfish sperm. Zygote 14, 2332.CrossRefGoogle ScholarPubMed
Singh, S., Lowe, D.G., Thorpe, D.S., Rodriguez, H., Kuang, W.J., Dangott, L.J., Chinkers, M., Goeddel, D.V. & Garbers, D.L. (1988). Membrane guanylate cyclase is a cell-surface receptor with homology to protein kinases. Nature 334, 708–12.CrossRefGoogle ScholarPubMed
Steiner, D.F. (1998). The proprotein convertases. Curr. Opin. Chem. Biol. 2, 3139.CrossRefGoogle ScholarPubMed
Thompson, J.D., Gibson, T.J., Plewniak, F., Jeanmougin, F. & Higgins, D.G. (1997). The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25, 4876–82.CrossRefGoogle ScholarPubMed
Vacquier, V.D. (1998). Evolution of gamete recognition proteins. Science 281, 1995–8.CrossRefGoogle ScholarPubMed
Ward, G.E., Brokaw, C.J., Garbers, D.L. & Vacquier, V.D. (1985). Chemotaxis of Arbacia punctulata spermatozoa to resact, a peptide from the egg jelly layer. J. Cell Biol. 101, 2324–9.CrossRefGoogle ScholarPubMed
Wilson, E.M. & Chinkers, M. (1995). Identification of sequences mediating guanylyl cyclase dimerization. Biochemistry 34, 4696–701.CrossRefGoogle ScholarPubMed
Figure 0

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.

Figure 1

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.

Figure 2

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.

Figure 3

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