MBE Advance Access originally published online on May 23, 2008
Molecular Biology and Evolution 2008 25(8):1705-1713; doi:10.1093/molbev/msn121
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Articles |
The Sperm Proteins from Amphioxus Mirror Its Basal Position among Chordates and Redefine the Origin of Vertebrate Protamines



* Department of Biochemistry and Microbiology, University of Victoria, Victoria, British Columbia, Canada
Departamento de Biología Celular y Molecular, Universidade da Coruña, Campus de A Zapateira s/n, A Coruña, Spain
Departament de Ciències Fisiològiques II, Facultat de Medicina, Campus de Bellvitge, Universitat de Barcelona, L'Hospitalet de Llobregat, Barcelona, Spain
Departament d'Enginyeria Química, ETSEIB, Universitat Politècnica de Catalunya, Diagonal, 647, Barcelona, Spain
E-mail: jausio{at}uvic.ca
| Abstract |
|---|
|
|
|---|
The sperm nuclear basic proteins (SNBPs) that participate in chromatin condensation in spermatozoa belong to 3 groups: histone (H), protamine-like (PL), and protamine (P) type. They share a common origin with histone H1 resulting from the segregation of PL components, corresponding to different regions of an H1 precursor molecule (N-terminal, winged-helix, C-terminal domains), becoming independent and following a subsequent process of parallel vertical evolution (H
PL
P). In the present work, we describe the sequence and primary structure of the main SNBP component in the sperm of the cephalochordate Branchiostoma floridae (amphioxus), revealing that it represents the deuterostome counterpart of the PL-III SNBP component from molluscs corresponding to the H1 N-terminal region. Until now, this has been a missing piece needed to complete the evolutionary history of SNBPs in metazoan genomes. The discovery of this PL lineage in deuterostomes definitively validates the parallel vertical evolution of SNBPs across metazoans, giving further support to the "basal" position of amphioxus among chordates, with respect to tunicates. Sequence analyses suggest that later on in evolution, the appearance of positively selected arginine-rich protamines, derived from the H1 C-terminal region, led to the extinction of this PL lineage in the genomes of early protostomes and deuterostomes. Given that tunicates are now viewed as a sister group of vertebrates, the lysine to arginine transition responsible for the origin of vertebrate protamines must be set a step back from tunicates.
Key Words: reproductive proteins protamines metazoans long-term evolution chromatin
| Introduction |
|---|
|
|
|---|
In somatic eukaryotic cells, the DNA is compacted within the nucleus by chromatin, a nucleoprotein complex that results from the association of DNA with basic chromosomal proteins known as histones (van Holde 1988
This observation represents an important starting point in assembling the pieces necessary to reconstruct the evolutionary histories of histone H1 and SNBPs. Piece one, we already evidenced that the long-term evolution of the histone H1 family has been subject to a birth-and-death process under strong purifying selection at the protein level, leading to the high levels of diversification and differentiation observed among the family members in higher eukaryotes (Eirín-López et al. 2004
; Nei and Rooney 2006
). Piece 2, the long-term evolution of SNBPs seems to be best described by a process of parallel vertical evolution (H
PL
P) subject to purifying selection at least in H-type and PL-type lysine-rich lineages in both protostomes and deuterostomes (Eirín-López, Lewis, et al. 2006
). Accordingly, only H- or PL-precursor SNBP types would be present in those taxa that arose early in metazoan evolution, whereas the more specialized PL and P types would represent a characteristic feature of those taxa located at the uppermost evolutionary branches of bilaterian evolution. Piece 3, the transition from lysine-rich PLs toward the highly specialized arginine-rich protamines was caused by a frameshift mutation in the PL component corresponding to the C-terminal region of H1 (Lewis et al. 2004
) and involved a shift from negative selection to an adaptive selection process (Eirín-López, Lewis, et al. 2006
).
However, an important piece is still needed in order to articulate and validate the whole evolutionary picture of these proteins. In one hand, functional PL components corresponding to the N-terminal region of an H1 precursor have been only identified in mollusc species (protostomes) (Rocchini et al. 1995
). This observation is in contrast to the expectations made by the model of parallel vertical evolution of SNBPs, which suggests that the differentiation of the 3 PL components occurred before the split between protostomes and deuterostomes. On the other hand, and more importantly, very recent analyses have traded the position of tunicates with that occupied by cephalochordates (also known as amphioxus or lancelets) in chordate evolution (Bourlat et al. 2006
), with the amphioxus now being viewed as the most "basal" chordate and the tunicates as the sister group of the vertebrates (Bourlat et al. 2006
; Holland 2007
).
In the present work, we found that the SNBP component in the sperm of the cephalochordate Branchiostoma floridae is exclusively composed of the PL fraction corresponding to the N-terminal region of H1. The discovery of this PL lineage in deuterostomes definitively validates the parallel vertical evolution of SNBPs across metazoans, giving further support to the basal position of the amphioxus among chordates, with respect to tunicates. The appearance of arginine-rich protamines later on in evolution led to the extinction of this PL lineage due to its inability to undergo a lysine to arginine transition.
| Materials and Methods |
|---|
|
|
|---|
SNBP Extraction and Fractionation
SNBPs were extracted from testes of B. floridae with 0.4 N HCl and precipitated with acetone as described in Wang and Ausió (2001)
In Silico Isolation of the Amphioxus PL Protein and Genomic DNA and cDNA Sequence Determination
Protein sequences corresponding to the PL component of the sperm of molluscs, tunicates and sea urchins were used in protein–protein basic local alignment and search tool (BLAST) searches against the complete genome sequence of B. floridae (unmasked v1.0, DOE Joint Genome Institute, http://www.jgi.doe.gov). The genomic region of the only SNBP candidate found (sequence ID 125082) was subsequently retrieved from the amphioxus database and used in the design of gene-specific primers. Total RNA was extracted from whole B. floridae animals using Trizol reagent (GibcoBRL, Burlington, ON), and mRNA was isolated using an mRNA purification kit (Amersham Bioscience, Piscataway, NJ). Genomic DNA was extracted using a DNeasy Tissue Kit (Qiagen, Mississauga, ON) following the manufacturers instructions. Gene-specific primers were designed using the genomic sequence retrieved from the amphioxus genome as a reference, and they were subsequently used for polymerase chain reaction with cDNA and genomic DNA templates. In addition, rapid amplification of cDNA ends (RACE) was employed to obtain the complete cDNA sequence using the First Choice RACE Kit (Ambion, Austin, TX). Sequencing was done by the DNA Sequencing Facility, Centre for Biomedical Research at the University of Victoria.
Molecular Evolutionary Analyses
A total of 206 amino acid sequences (see supplementary table, Supplementary Material online) were used in the analyses, including 91 nonredundant histone H1 somatic sequences (68 replication dependent [RD] and 23 replication independent [RI]), 1 oocyte-specific H1, 16 testis-specific H1, 97 SNBP sequences (3 H type, 18 PL type, and 76 P type), and the PL protein from amphioxus. Multiple alignments of the amino acid sequences were conducted using the BIOEDIT (Hall 1999
) program, and all molecular evolutionary analyses were conducted using the program MEGA version 3.1 (Kumar et al. 2004
). The extent of the amino acid divergence was estimated by means of the uncorrected differences (p distance), and the minimum evolution tree-building method (Saitou and Nei 1987
) was used to reconstruct the phylogenetic trees. In order to assess that our results were not dependent on this method, the simplified phylogeny was compared with a maximum parsimony tree reconstructed using the close-neighbor-interchange tree search method. The reliability of the topologies was tested by the bootstrap method (Felsenstein 1985
) and the interior branch test (confidence probability) method (Sitnikova 1996
), based on 1,000 replicates. Codon usage analyses on PL, P, and H1 genes were performed by using the DnaSP version 3 program (Rozas et al. 2003
), using the "effective number of codons" (ENC) to estimate the codon usage bias (Wright 1990
).
| Results and Discussion |
|---|
|
|
|---|
The Main Sperm Chromosomal Protein of Amphioxus Corresponds to the N-terminal Region of a Histone H1 Precursor
Figure 1A shows the SNBP composition of B. floridae (lane 5) in comparison to that of representative species from different classes of echinoderms and a tunicate. As seen in this figure (lanes 1–3), echinoderm SNBPs belong to the H type with the exception of the class Holoturoidea, where a short PL (protein Phi-0) (Subirana 1970
|
Partial information on the N-terminal residues of the main SNBP component of B. floridae was obtained (Saperas et al. 1994
|
The sequence of the gene encoding this putative SNBP was retrieved from the B. floridae genome database and had 3 exons in the coding region (fig. 2B). The full-length protamine cDNA was subsequently obtained, assessing the existence and functionality of this protein in vivo in amphioxus, consisting of a coding region which was 333 nt in length, and encoding a protein with an amino acid composition virtually identical to that previously obtained (Saperas et al. 1994
Molecular Evolutionary Genomics of SNBPs
The hypothesis for the parallel vertical evolution of SNBPs (Ausió 1999
) predicts that only H- or PL-precursor SNBP types would be found in those taxa that arose early in metazoan evolution. The phylogenetic relationships among H1 and SNBPs reconstructed in the present work (fig. 3) clearly support this notion, depicting an evolutionary diversification process that mirrors the evolution of triploblastic animals. In fact, histone H1 and protamines appear as 2 clearly differentiated groups, with PLs occupying an intermediate position. The occurrence of RD and RI H1 histones, as well as PL- and P-type SNBPs in protostome and deuterostome representatives, provides support that their evolutionary origin occurred prior to the split between these 2 groups of organisms.
|
The topology shown in figure 3 is in good agreement with the taxonomic relationships among species, indicating a monophyletic origin for the group of RI H1s and a functional clustering pattern for the different mammalian somatic H1 isoforms (Eirín-López, Lewis, et al. 2006
|
The differentiation of the amphioxus SNBP lineage reported in figure 3 agrees with the evidence provided by recent reports supporting the status of the amphioxus as the most basal chordate, in exchange of tunicates (Bourlat et al. 2006
The Amphioxus PL Lineage Is Subject to Different Evolutionary Constraints than the Protamines or the PL-I Lineage from Tunicates
The progressive specialization of SNBPs has followed a parallel vertical pattern of evolution defined as H
PL
P (Ausió 1999
; Eirín-López, Frehlick, and Ausió 2006
), encompassing a progressive replacement of lysines (K) by arginines (R) leading to a more efficient DNA packaging in the sperm as well as in the activation of enzymatic cascades after fertilization in vertebrates (Ausió et al. 1984
; Ohtsuki et al. 1996
; Rooney and Zhang 1999
). The gradual variation in the K/R ratio across H1/SNBP evolution is evident in the phylogeny shown in figure 3, ranging from around 8.7 in the case of H1 to 0.1 in the case of protamines. Furthermore, the increase in arginine is concomitant with the increase in the numbers of nonsynonymous substitutions among P lineages (see fig. 3), suggesting that this trait is positively selected.
It has been recently shown that a frameshift mutation involving only 2 nt would be enough to duplicate the percent of arginine in the PL component of the tunicate Ciona intestinalis, due to the high frequency of AAG and AGA triplets encoding polylysine and polyarginine tracts, respectively (Lewis et al. 2004
). By studying the relative codon usage for lysine and arginine residues (fig. 5A), it is evident that, although the SNBP from B. floridae presents similar amounts of lysine and arginine (as does PL-III from molluscs), all lysines are encoded by the same codon (AAA) with a higher proportion of arginine residues encoded by CGU and CGC. This codon preference is in sharp contrast to that shown by histone H1 and the PL components from molluscs, where arginines are mainly encoded by AGA and AGG triplets, and would ultimately represent an insurmountable barrier for a lysine to arginine transition based on a frameshift mutation.
|
Although molluscan PLs and the amphioxus SNBP are very biased in terms of their codon usage, somatic histone H1s and tunicate PLs exhibit intermediate and very low codon bias, respectively. This may be explained by the fact that although PL proteins rely equally on lysines and arginines for DNA packaging, constraints acting on protamines are only dependent on the overall arginine content of the molecule, rather than on the specific position of the arginines or the nature of their encoding triplets (Rooney et al. 2000
The PL component from Amphioxus Provides a Footprint in the Process of Parallel Evolution Observed in SNBPs across Metazoan Genomes
There is now plenty of evidence indicating that the exclusion of specialized H1 genes from the main repetitive histone units to a solitary location in the genome resulted in the generation of an "orphon" group of H1 genes early in the evolution of metazoans (Eirín-López et al. 2002
, 2004
). The independent evolution of this group eventually led to the differentiation of RD and RI H1 variants, followed by the subsequent differentiation of SNBPs within the second lineage (Eirín-López, Lewis, et al. 2006
). During the differentiation and diversification of these chromosomal proteins, the transition toward arginine-rich protamines appears to have involved in the progressive specialization of a PL protein corresponding to the C-terminal region of H1 (Lewis et al. 2004
), which eventually suffered a shift from a purifying selection process (present in PLs) to adaptive selection favoring high arginine contents in protamines (fig. 5B) (Balhorn 2007
). However, the presence and fate of the SNBP component corresponding to the N-terminal part of H1 have remained elusive until now.
| Conclusions |
|---|
|
|
|---|
The identification of a deuterostome counterpart (represented by the amphioxus PL reported here) to the protostome PL component corresponding to the N-terminal region of a precursor histone H1 (represented by the PL-III protein from molluscs) provides a footprint that definitively validates the parallel vertical evolution of SNBPs (summarized in fig. 6). This result is even more critically relevant in light of recent studies placing the amphioxus in the basal position in chordate evolution (Bourlat et al. 2006
|
The SNBP component corresponding to the H1 N-terminal region appears to be the initial choice for sperm DNA packaging in early metazoans. The extinction of this PL lineage in the genome was probably due to its inability to undergo a lysine to arginine transition based on the observed patterns of codon usage bias for lysine residues. In this regard, the PL-III protein from molluscs and the SNBP component identified in amphioxus would represent end points of the SNBP evolutionary process in deuterostome and protostome genomes that were promptly discarded by the fast specialization experienced by arginine-rich PL- and P-type SNBPs which arose through a process of positive selection.
| Supplementary Material |
|---|
|
|
|---|
Supplementary table and figure are available at Molecular Biology and Evolution online (http://www.mbe.oxfordjournals.org/). The sequence described in the present work has been deposited in the GenBank database with accession number EU271675.
| Acknowledgements |
|---|
|
|
|---|
We are very thankful to Linda Holland of the Marine Biology Research Division at the Scripps Institution of Oceanography (UCSD) for providing us with specimens of B. floridae during the early biochemical characterization of SNBPs and to Nancy Sherwood and her laboratory at the Department of Biology at UVic for providing us with material for the molecular biology analyses. This work was supported by Natural Science and Engineering Research Council of Canada (NSERC) Grant OGP 0046399-02 to J.A. and by Grants BFU 2005-00123 and BFU2006-04035/BMC from the Spanish MCYT PGC-FEDER. J.M.E.-L. is a recipient of a Postdoctoral Marie Curie International Fellowship within the sixth European Community Framework Program. L.J.F. holds an NSERC PGSD fellowship.
| Footnotes |
|---|
Billie Swalla, Associate Editor
| References |
|---|
|
|
|---|
Ausió J. Histone H1 and the evolution of the nuclear sperm-specific proteins. In: Advances in spermatozoal phylogeny and taxonomy—Jamieson BGM, Ausió J, Justine JL, eds. (1995) Memoires du Museum National d'Histoire Naturelle. Paris (France). p. 447–462.
Ausió J. An unusual cysteine-containing histone H1-like protein and two protamine-like proteins are the major nuclear proteins of the sperm of the bivalve mollusc Macoma nasuta. J Biol Chem (1988) 263:10141–10150.
Ausió J. Histone H1 and evolution of sperm nuclear basic proteins. J Biol Chem (1999) 274:31115–31118.
Ausió J, Eirín-López JM, Frehlick LJ. Evolution of vertebrate chromosomal sperm proteins: implications for fertility and sperm competition. Soc Reprod Fertil Suppl (2007) 65:63–79.[Medline]
Ausió J, Greulich KO, Watchel E. Characterization of the fluorescence of the protamine thynnine and studies of binding to double-stranded DNA. Biopolymers (1984) 23:2559–2571.[Web of Science][Medline]
Azorin F, Olivares C, Jordan A, Perez-Grau L, Cornudella L, Subirana JA. Heterogeneity of the histone-containing chromatin in sea cucumber spermatozoa. Distribution of the basic protein phi 0 and absence of non-histone proteins. Exp Cell Res (1983) 148:331–344.[CrossRef][Web of Science][Medline]
Balhorn R. The protamine family of sperm nuclear proteins. Genome Biol (2007) 8:227.[CrossRef][Medline]
Bourlat SJ, Juliusdottir T, Lowe CJ, et al, (14 co-authors). Deuterostome phylogeny reveals monophyletic chordates and the new phylum Xenoturbellida. Nature (2006) 444:85–88.[CrossRef][Medline]
Chiva M, Saperas N, Cáceres C, Ausió J. Nuclear basic proteins from the sperm of tunicates, cephalochordates, agnathans and fish. In: Advances in spermatazoal phylogeny and taxonomy—Jamieson BGM, Ausió J, Justine JL, eds. (1995) Memoires du Museum National d'Histoire Naturelle. Paris (France). p. 501–514.
Dunn CW, Hejnol A, Matus DQ, et al, (18 co-authors). Broad phylogenomics sampling improves resolution of the animal tree of life. Nature (2008) 452:745–749.[CrossRef][Medline]
Eirín-López JM, Frehlick LJ, Ausió J. Protamines, in the footsteps of linker histone evolution. J Biol Chem (2006) 281:1–4.
Eirín-López JM, González-Tizón AM, Martínez A, Méndez J. Molecular and evolutionary analysis of mussel histone genes (Mytilus spp.): possible evidence of an "orphon origin" for H1 histone genes. J Mol Evol (2002) 55:272–283.[CrossRef][Web of Science][Medline]
Eirín-López JM, González-Tizón AM, Martínez A, Méndez J. Birth-and-death evolution with strong purifying selection in the histone H1 multigene family and the origin of orphon H1 genes. Mol Biol Evol (2004) 21:1992–2003.
Eirín-López JM, Lewis JD, Howe L, Ausió J. Common phylogenetic origin of protamine-like (PL) proteins and histone H1: evidence from bivalve PL genes. Mol Biol Evol (2006) 23:1304–1317.
Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution Int J Org Evolution (1985) 39:783–791.[CrossRef][Web of Science]
Frehlick LJ, Eirín-López JM, Prado A, Su HW, Kasinsky HE, Ausió J. Sperm nuclear basic proteins of two closely related species of Scorpaeniform fish (Sebastes maliger, Sebastolobus sp.) with different sexual reproduction and the evolution of fish protamines. J Exp Zoolog A Comp Exp Biol (2006) 305:277–287.[Medline]
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser (1999) 41:95–98.
Holland LZ. A chordate with a difference. Nature (2007) 447:153–155.[CrossRef][Medline]
Jutglar L, Borrell JI, Ausió J. Primary, secondary, and tertiary structure of the core of a histone H1-like protein from the sperm of Mytilus. J Biol Chem (1991) 266:8184–8191.
Kasinsky HE, Lewis JD, Dacks JB, Ausió J. Origin of H1 histones. FASEB J (2001) 15:34–42.
Kumar S, Tamura K, Nei M. MEGA3: integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform (2004) 5:150–163.
Lewis JD, Saperas N, Song Y, Zamora MJ, Chiva M, Ausió J. Histone H1 and the origin of protamines. Proc Natl Acad Sci USA (2004) 101:4148–4152.
Lewis JD, Song Y, de Jong M, Bagha S, Ausió J. A walk though vertebrate and invertebrate protamines. Chromosoma (2003) 111:473–482.[Web of Science][Medline]
Nei M, Rooney AP. Concerted and birth-and-death evolution in multigene families. Annu Rev Genet (2006) 39:121–152.[CrossRef][Web of Science]
Ohtsuki K, Nishikawa Y, Saito H, Munakata H, Kato T. DNA-binding sperm proteins with oligo-arginine clusters function as potent activators for egg CK-II. FEBS Lett (1996) 378:115–120.[CrossRef][Web of Science][Medline]
Prats E, Cornudella L. The gene encoding the sperm-specific basic nuclear protein fo from sea cucumber. In: Advances in spermatozoal phylogeny and taxonomy. Memoires du Museum National d'Histoire Naturelle—Jamieson BGM, Ausió J, Justine JL, eds. (1995) Paris (France): NHBS Evironment Bookstore. 491–500.
Rocchini C, Rice P, Ausió J. Complete sequence and characterization of the major sperm nuclear basic protein from Mytilus trossulus. FEBS Lett (1995) 363:37–40.[CrossRef][Web of Science][Medline]
Rooney AP, Zhang J. Rapid evolution of a primate sperm protein: relaxation of functional constraint or positive Darwinian selection? Mol Biol Evol (1999) 16:706–710.[Abstract]
Rooney AP, Zhang J, Nei M. An unusual form of purifying selection in a sperm protein. Mol Biol Evol (2000) 17:278–283.
Rozas J, Sánchez-del Barrio JC, Messeguer X, Rozas P. DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics (2003) 19:2496–2497.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol (1987) 4:406–425.[Abstract]
Saperas N, Chiva M, Ausió J. Purification and characterization of the protamines and related proteins from the sperm of a tunicate, Styela plicata. Comp Biochem Physiol (1992) 103:969–974.[CrossRef]
Saperas N, Chiva M, Ribes E, Kasinsky HE, Rosenberg E, Youson JH, Ausió J. Chromosomal proteins of the sperm of a cephalochordate (Branchiostoma floridae) and an agnathan (Petromyzon marinus). Compositional variability of the nuclear sperm proteins of deuterostomes. Biol Bull (1994) 186:101–114.[Abstract]
Sitnikova T. Bootstrap method of interior-branch test for phylogenetic trees. Mol Biol Evol (1996) 13:605–611.[Abstract]
Subirana JA. Nuclear proteins from a somatic and a germinal tissue of the echinoderm Holothuria tubulosa. Exp Cell Res (1970) 63:253–260.[CrossRef][Medline]
Telford MJ. Resolving animal phylogeny: a sledgehammer for a tough nut? Dev Cell (2008) 14:457–459.[CrossRef][Web of Science][Medline]
van Holde KE. Chromatin (1988) New York: Springer.
Wang X, Ausió J. Histones are the major chromosomal protein components of the sperm of the nemerteans Cerebratulus californiensis and Cerebratulus lacteus. J Exp Zool (2001) 290:431–436.[CrossRef][Web of Science][Medline]
Wright F. The effective number of codons used in a gene. Gene (1990) 87:23–29.[CrossRef][Web of Science][Medline]
![]()
CiteULike
Connotea
Del.icio.us What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||





approximate).