Skip to main content

Proteome analysis of shell matrix proteins in the brachiopod Laqueus rubellus



The calcitic brachipod shells contain proteins that play pivotal roles in shell formation and are important in understanding the evolution of biomineralization. Here, we performed a large-scale exploration of shell matrix proteins in the brachiopod Laqueus rubellus.


A total of 40 proteins from the shell were identified. Apart from five proteins, i.e., ICP-1, MSP130, a cysteine protease, a superoxide dismutase, and actin, all other proteins identified had no homologues in public databases. Among these unknown proteins, one shell matrix protein was identified with a domain architecture that includes a NAD(P) binding domain, an ABC-type transport system, a transmembrane region, and an aspartic acid rich region, which has not been detected in other biominerals. We also identified pectin lyase-like, trypsin inhibitor, and saposin B functional domains in the amino acid sequences of the shell matrix proteins. The repertoire of brachiopod shell matrix proteins also contains two basic amino acid-rich proteins and proteins that have a variety of repeat sequences.


Our study suggests an independent origin and unique mechanisms for brachiopod shell formation.


Many organisms have a diversity of biominerals that have a large number of ecologically important functions, including body support and defense against predators. Most metazoan biominerals are likely acquired independently in different phyla in the early Cambrian [1]. On the other hand, presence of some biomineralization genes shared by several phyla suggests that at least a part of the mechanisms was derived from a common ancestor [2]. A crucial factor in understanding the mechanisms and evolution of biomineralization is an in-depth analysis of the organic matrix present in hard tissues. A large number of skeletal proteins have been identified and their roles in the biomineralization process have been analyzed [3]. Recently, it became possible to identify skeletal proteins comprehensively using transcriptome analysis combined with proteome analysis across a wide range of phyla [49]. Brachiopods are marine invertebrates that appeared in the Cambrian and have two valves composed of calcium carbonate or calcium phosphate. Organic shell matrices in brachiopods have been studied in some detail [10, 11]. Previous SDS-PAGE analyses showed several major bands stained by Coomassie Brilliant Blue (CBB) [10, 12], and an amino acid composition analysis revealed that glycine and alanine account for a large proportion of bulk shell extracts [13, 14]. In addition to these studies, shell matrix proteins of brachiopods have been studied for the characterization of peptides and amino acids preserved in fossil brachiopod shells and immunological assays showed that some peptide sequences in the shells are preserved for about 1 million years [15, 16]. However, the amino acid sequences of the shell matrix proteins in brachiopods have not been identified, except for the partial amino acid sequencing of a chromoprotein, named ICP-1 [1719]. In this study, we performed a large-scale analysis of shell matrix proteins of the brachiopod Laqueus rubellus using proteomics combined with transcriptomics. The analysis identified 40 shell matrix proteins. These new proteome data contained 35 protein sequences that have no database-related homologues and indicate that the mechanism and evolution of brachiopod shell formation are unique.


SDS-PAGE of the soluble organic matrix extracted from the secondary layer showed two major bands (< 6.5 and 12 kDa) and a minor band of 62 kDa when stained with CBB (Fig. 1a). Silver staining showed a major band at 35 kDa and many minor bands in addition to the three bands stained with CBB (Fig. 1b). Transcriptome sequencing in mantle tissue generated a total of 125,437 reads. To identify shell matrix proteins from Laqueus rubellus, we performed a proteome analysis of the material from these different sources, namely (1) the soluble organic matrix from the whole shell, (2) insoluble organic matrix from the whole shell and (3) soluble organic matrix from the shell secondary layer. The generated MS/MS data were subjected to SEQUEST searches against the protein sequence database obtained from the transcriptome analysis. As a result, we identified 40 shell matrix proteins (Tables 1, 2, 3, 4 and Fig. 2). Among these proteins, 18 proteins were identified from the soluble organic matrix of the shell secondary layer, and the soluble and insoluble organic matrix from the whole shell (Table 1). While five proteins were identified from both the soluble and insoluble organic matrix of the whole shell, 17 proteins were only identified from either the soluble or insoluble organic matrices of the whole shell (Tables 2, 3, 4). A total of 22 proteins were deduced to be possibly complete sequences, because these sequences have a stop codon and either an in-frame start codon just after an in-frame stop codon in the 5′ region of the sequence or a potential signal peptide (Tables 1, 2, 3, 4). A blast search against the GenBank non-redundant database showed that 36 out of the 40 shell matrix proteins have no annotated homologous sequences. Among the proteins that have sequence homology, isotig 00281 showed relatively high sequence similarity with MSP130, which is a skeletal protein identified from sea urchins (Fig. 3) [20, 21]. Isotig 01587, isotig 00949, and isotig 00959 showed high sequence similarity with extracellular copper/zinc superoxide dismutase, actin I, and cathepsin L cysteine proteinase, respectively (Tables 2 and 4). Predicted molecular mass and isoelectric point are shown in Tables 1, 2, 3, 4, and domains predicted by the SMART program and a National Center for Biotechnology Information (NCBI) conserved domain search are schematically shown in Fig. 2. Using the known amino acid sequences of ICP-1 [19], we performed local blast searches against sequence data obtained from this study. As a result, isotig 00046 was found to have a high sequence similarity with ICP-1 (Fig. 4). Gene expression analysis showed that ICP-1 gene is expressed in lophophore as well as in mantle tissues (Additional file 1). To estimate the abundance of each protein, we calculated the relative copy number based on the identified spectral counting [2225]. The number of identified spectra of each protein was divided by the number of theoretically observable tryptic peptide ions, which have a mass-to-charge ratio of 400 to 1500 at two or three charge states, to generate an abundance index as the relative copy number. The result showed that isotig 00046 (ICP-1) is the most abundant protein in the shell extracts among the proteins identified in this study (Fig. 5). The percentage of spectra that were matched to peptide contained in protein sequences translated from transcriptome data was about 10 % of the total MS/MS spectra acquired. The proportion in this case was lower than those of other mass spectrometric analyses for organisms of which genome has been completely sequenced (S. cerevisiae, 40-50 %), suggesting presence of many more proteins in shell matrix samples than proteins identified in this study. Messenger RNA sequences in the transcriptome data may not include complete list for proteins in shell matrix or unknown chemical modifications that were not considered in our database search may occur in a large proportion of shell matrix proteins. The depth of isotigs is shown in Additional file 2.

Fig. 1

EDTA-soluble extracts from shell secondary layer of Laqueus rubellus were fractionated by SDS-PAGE. a: CBB b: silver staining M: marker

Table 1 Shell matrix proteins identified from the whole shell and shell secondary layer
Table 2 Shell matrix proteins identified from the whole shell (the soluble and insoluble organic matrix)
Table 3 Shell matrix proteins identified from the whole shell (the soluble organic matrix)
Table 4 Shell matrix proteins identified from the whole shell (the insoluble organic matrix)
Fig. 2

Schematic of the domains in shell matrix proteins identified in this study. CC: coiled coil; PL: pectin lyase-like; Vir28: variable surface protein Vir28; TI: trypsin inhibitor-like cysteine-rich domain; NAD(P): NAD(P)-binding Rossmann-fold domains; ABC-TS: ABC-type transport system; TR: transmembrane region; PAMG: Pneumovirinae attachment membrane glycoprotein G; SOD: copper/zinc superoxide dismutase; SapB: saposin B domains; Ac: Actin; I29: cathepsin propeptide inhibitor domain; and Papain: papain family cysteine protease

Fig. 3

Alignment of the amino acid sequences of MSP-130 and isotig 00281. Sk: Saccoglossus kowalevskii (NCBI Acc. No. XP_002739468.1); Sp: Strongylocentrotus purpuratus (NCBI Acc. No. NP_001116986.1); He: Heliocidaris erythrogramma (NCBI Acc. No. CAC20358.1); and Cg: Crassostrea gigas (NCBI Acc. No. EKC20477.1)

Fig. 4

Alignment of the amino acid sequences of ICP-1 and isotig 00046. Nl: Neothyris lenticularis; Ts: Terebratella sanguinea; and Ci: Calloria inconspicua

Fig. 5

Abundance index of shell matrix proteins identified in this study


Repertoire of matrix proteins found in brachiopod shells


ICP-1 (IntraCrystalline Protein-1) is a shell matrix protein extract from the calcitic shell of three brachiopod species: Neothyris lenticularis, Calloria inconspicua, and Terebratella sanguinea. The partial N-terminal amino acid sequences of ICP-1 from these species have been determined by Edman degradation [1719]. The local blast searches showed that the partial sequences of ICP-1 have a sequence similarity with a part of isotig 00046 (Fig. 4). Isotig 00046 showed the highest abundance in the shell among the shell matrix proteins identified in this study (Fig. 5). This is consistent with the results of a previous study [19]. ICP-1 was originally identified as a 6.5-kDa major band revealed by SDS-PAGE analysis of shell extracts from Neothyris lenticularis, Calloria inconspicua, and Terebratella sanguinea [19]. The predicted molecular mass of isotig 00046 is just over 11.9 kDa, while the thickest band in the SDS-PAGE analysis of Laqueus rubellus is < 6.5 kDa (Fig. 1). One possibility to explain this size discrepancy is that ICP-1 undergoes proteolytic cleavage after translation. We have identified cathepsin L cysteine proteinase from the shell extracts. This protein could be involved in post-translational modification. HPLC analysis in a previous study indicated that carotenoids are bound to ICP-1 [19]. The observation that ICP-1 shows the highest abundance in the shell suggests that ICP-1 plays key roles in the biomineralization processes in brachiopods. RT-PCR analysis showed that ICP-1 gene is also expressed in the lophophore tissues, suggesting that ICP-1 is also involved in the formation of the calcareous loop structure embraced by the lophophore.


Isotig 00281 showed relatively high sequence similarity with MSP130 (The e-value against MSP130 from Saccoglossus kowalevskii was 1e−10). MSP130 was originally identified from primary mesenchyme cells in the sea urchin [20, 21] and was subsequently detected in the hard tissues of sea urchins [26, 27]. Homologues and closely related proteins of MSP130 were also reported to be present in molluscan shells [28]. In addition, MSP130 have been found in genomes of hemichordate, cephalochordate, bacteria, and green algae [29]. MSP130 is predicted to have been acquired by independent horizontal gene transfer in Cambrian, because this gene exists in bacteria and has an extremely wide phylogenetic distribution [29]. However, it appears possible that MSP130 gene was already present in the metazoan or bilaterian last common ancestor because the phylogenetic tree of MSP130 constructed in the previous study was divided into two clusters of the bacterial/green algae clade and the metazoan (bilaterian) clade [29]. If the MSP130 gene was transferred to animals horizontally many times independently, the animal MSP130 genes would not form a monophyletic cluster. The functions of MSP130 have not been determined, but this protein is predicted to function at the cell surface [30]. Identification of MSP130 from the brachiopod shells in addition to the skeletons of sea urchin and molluscs suggests that this protein plays an important role in biomineralization processes.

Digestive enzymes and inhibitors

Many digestive enzymes have been identified from the shell of Laqueus rubellus. Isotig 00227 has a pectin lyase-like domain, which enzymatically breaks down pectin. Pectin is usually found in plant cells, but some neutral sugars have been detected in terebratulide brachiopod shells [12]. Thus, isotig 00227 could be involved in the breakdown of these neutral sugars. Two lysosomal proteins were identified from the shell: (i) isotig 00213 is a saposin protein that is involved in lipid degradation and (ii) isotig 00959 is cathepsin L, which is a cysteine protease that plays a major role in intracellular protein catabolism. However, a trypsin inhibitor domain was found in the amino acid sequence of isotig 01016. A number of protease inhibitors have been identified from molluscan shells and a possible function of these proteins could be to protect the shell matrix proteins against particular proteases [31, 32]. In brachiopod shell formation, a protease inhibitor could play a similar role, but since secreted proteases and protease inhibitors are involved in modification of extracellular and membrane bound proteins in many systems, their roles may not be specific to shell matrix modifications.

Membrane protein

Isotig 01176 has a transmembrane region, indicating that this protein binds to the cell membranes of the mantle tissue. Interestingly, several skeletal proteins of a coral also have a transmembrane region [8]. If this protein is involved in shell formation, the mantle epithelium would be closely attached to the inner shell surface when mineralization occurs. Isotig 01176 also has a NAD(P) binding domain. NADP is a hydrogen carrier that is used in metabolic pathways such as the photosynthetic pathway and glycolysis as a reducing agent. Although the function of this domain in isotig 01176 is unknown, the high importance of H+ in CaCO3 synthesis suggests that the NAD(P) binding domain of isotig 01176 could function to sequester H+ from the crystallization milieu. Isotig 01176 may also have an ABC-type transport system, which is involved in the transportation of many types of substrates. This suggests that isotig 01176 could also control ion concentrations in the space where crystallization occurs. As isotig 01176 has all these functional domains, this protein is likely to have key roles in brachiopod shell formation.

Secreted protein

Among the 40 shell matrix proteins, 18 have a signal peptide, indicating that many shell matrix proteins are secreted from cells. Generating the calcium carbonate in the extracellular space is consistent with brachiopod shells being an exoskeleton. However, seven proteins do not have a signal peptide even though they have complete sequences. One possibility is that these proteins bind with other proteins that have a signal peptide. Actin, which is one of the major intracellular proteins, was identified. However, there is a possibility that this protein is occluded from trapping of proteins involved in the secretory processes.

Amino acid composition and isoelectric points

Shell matrix proteins of molluscs contain aspartic acid-rich proteins that are postulated to interact with Ca2+ [33]. Isotig 01176 has a repeat sequence comprising aspartic acid, suggesting that this region also binds Ca2+ (Fig. 6). However, the brachiopod shell matrix proteins identified in this study do not have unusually low isoelectric points seen for aspartic acid-rich proteins in the molluscan shell (Fig. 7). Among the shell matrix proteins identified in this study, isotig 01521 and isotig 02158 have a relatively high concentration of basic amino acids (Fig. 6). Although it is possible that unusually acidic shell matrix proteins have not been identified in this study, there may exist a general difference in isoelectric points between mollusc and brachiopod shell matrix proteins.

Fig. 6

Amino acid sequences of isotig 01176, isotig 01521 and isotig 02158. Acidic amino acids are highlighted in red and basic amino acids are highlighted in blue

Fig. 7

Molecular mass and isoelectric points. Red symbols represent shell matrix proteins from brachiopod and gray symbols represent shell matrix proteins from mollusca (Marin et al. 2008)

Repeat sequences

Many of the shell matrix proteins identified in this study have repeat sequences, which are generally represented as a low complexity domain in the SMART prediction (Fig. 2). Repeat sequences exist in many skeletal proteins identified from other phyla and are thought to have important roles in biomineralization processes. The GXN (Glycine-X-Asparagine) repeat that is present in the amino acid sequence of isotig 02671 also exists in Nacrein, the carbonic anhydrase in the molluscan shell (Fig. 8). This repeat sequence has been proposed to function by inhibiting the precipitation of calcium carbonate [34]. Therefore, the GXN repeat in isotig 02671 could also be involved in the control of CaCO3 growth. Besides this repeat sequence, the repeat motifs PPRG, GGX, and GGQNTGX are also present in sequences of the shell matrix proteins (Fig. 8). Although the exact function of these repeat sequences is not clear, the existence of a variety of repeat sequences in skeletal proteins suggests that this sequence structure has fundamental roles in biomineralization.

Fig. 8

Repeat motifs found in shell matrix proteins identified in this study

Possible mechanisms of brachiopod biomineralization

Many shell matrix proteins identified in this study have a signal peptide, suggesting that crystallization occurs in the outer cell space. This is consistent with shell matrix proteins in mollusca. Isotig 01176, which has an asparatic acid-rich region, NAD(P) binding domain, and ABC-type transport system, is a possible candidate that controls Ca2+ and H+ levels, whereas isotig 01521, which contains many basic amino acids, may interact with HCO3 . These proteins do not exist in skeletal proteins of other phylum, suggesting that the brachiopod has a unique ion control system for shell formation. In addition to these proteins, MSP130 and ICP-1 are thought to have important roles in shell formation.

Evolution of brachiopod biomineralization

Skeletal proteins from brachiopods and other phyla share common structural features. The existence of MSP130 in brachiopod shell suggests that this protein have significant functions for biomineralization, and the several phyla share a part of biomineralization mechanisms. This protein could have been acquired by horizontal gene transfer from bacteria, but as discussed above, the data shown by Ettensohn (2014) are not conclusive, and MSP130 could also have been present in the last common ancestor of bilaterians or metazoans. On the other hand, proteinase and proteinase inhibitors did not show significant sequence similarities to the skeletal proteins in other phyla, suggesting that these proteins have been recruited as shell matrix proteins independently. In addition, almost all the brachiopod shell matrix proteins identified in this study had no homologous sequences in the public data banks, suggesting unique origins for those proteins. These results support the hypothesis that brachiopod biominerals were acquired independently from other phyla in the Cambrian explosion. On the other hand, there is also a possibility that brachiopod biominerals share the same ancestral biomineralization system with some other phyla, but after acquisition of biominerals, the matrix proteins evolved rapidly and/or novel lineage-specific proteins were added to the ancestral biomineralization system. The shell proteome of the mollusc Lottia gigantea showed that 25 out of 39 shell matrix proteins had no homologues, and only a few proteins showed high sequence similarity to shell matrix proteins in bivalves and skeletal proteins in other phyla [7]. This observation may well be reflecting the rapid nature of the evolution of shell matrix proteins, but can also be considered as reflecting the independent origin of the skeletal matrix proteins of molluscs from those of other phyla. Although more data are needed to address this problem, extremely low numbers of homologous shell matrix proteins between brachiopods and molluscs, combined with the presence of a possible unique ion control system involving basic shell matrix proteins, tend to support an independent origin for the brachiopod shells as expected from the phylogenetic relationships and the fossil record [1, 35].


Our results identified two interesting shell matrix proteins, ICP-1 and MSP130. ICP-1 is a brachiopod shell matrix protein sequenced partially in previous studies, and MSP130 is a skeletal protein identified originally in sea urchins and oysters. Our data also showed novel shell proteins containing unique structures, including NAD(P) binding domains, suggesting the involvement of a hitherto unknown ion control system for shell formation. In addition, most other shell matrix proteins of Laqueus rubellus do not have a homologue in skeletal proteins of other phyla, suggesting an independent origin of the brachiopod shell. To further address the mechanisms and evolution of shell formation in brachiopods, additional studies, such as gene expression analysis and functional analysis, using this large-scale sequence information are necessary.

Materials and methods

RNA extraction and amplification

Live individuals of Laqueus rubellus were collected by dredging operations in Sagami Bay, Japan (2 km off Jogashima, 90 m water depth). The mantle tissues from a single individual were separated from the shell using tweezers and homogenized in 500 μL of Isogen (Nippon Gene, Tokyo, Japan) in a 1.5-mL microcentrifuge tube. Total RNA was extracted following the manufacturer’s protocol and purified using the RNeasy Mini Kit (Qiagen, Hilden, Germany). Amplification of mRNA was performed using MessageAmp II aRNA Amplification Kit (Life Technologies, Carlsbad, CA, USA) to obtain sufficient quantities of mRNA (200 ng) for transcriptome analysis. The amount of the initial total RNA subjected to the amplification was 1 μg and the incubation time for in vitro transcription was 8 h.


The nucleotide sequences of the cDNA expressed in the mantle tissue were determined using 454 GS Junior (Roche, Basel, Switzerland). Template DNA was prepared according to the supplier’s protocol. A total of 125,437 reads of the GS Junior sequences were generated. The obtained reads were assembled using Newbler v.2.8. software with default settings, resulting in 2,342 isotigs.

Extraction of shell matrix proteins

The shells were incubated overnight in a 5 % bleach solution with gentle shaking at room temperature to remove surface contaminants. After thorough washing with ultrapure water, the shells were crushed in water and dried. Organic materials were extracted by dissolution of the calcium carbonate using 0.5 M EDTA (pH 8.0) at a ratio of 23 mL to 1 g shell with shaking at 4 °C. After the solution was centrifuged at 20,000 g for 1 h, the insoluble organic materials separated by the centrifugation step were dissolved in an aqueous solution containing 9 M urea and 2 % (v/v) Triton X-100. The supernatant was concentrated and the EDTA removed using an Amicon Ultra-15 centrifugal filter unit with an Ultracel-3 membrane (Millipore, Billerica, CA, USA). The solution was then centrifuged at 15,000 g for 15 min and the supernatant was used for proteome analysis to identify the shell proteins in the soluble organic matrix. In addition to the matrix proteins of the whole shell, we also extracted proteins contained specifically in the shell secondary layer. The fibers of the shell secondary layer were collected by decantation after crushing the shells in water, and the soluble organic matrix was collected by the same method described above.

SDS-PAGE analysis

To check the concentration and heterogeneity of the extracted proteins, the soluble organic matrix from the shell secondary layer was subjected to SDS-PAGE analysis using 4–20 % Mini-PROTEAN TGX Precast Gels (Bio-Rad, Hercules, CA, USA) with the standard method [36], and stained with CBB and silver staining.

Peptide sample preparation by trypsin digestion

The soluble matrix proteins extracted with 0.5 M EDTA and the proteins that was insoluble in the first extraction step with EDTA and re-dissolved in solution containing 9 M urea and 2 % (v/v) triton X-100 were diluted in 0.1 M Tris–HCl (pH 8.5) to a final volume of 200 μL. Methanol (600 μL), chloroform (150 μL), and distilled water (450 μL) were added one-by-one and mixed thoroughly. After centrifugation at 12,000 rpm for 5 min at 4 °C, the upper aqueous layer was removed and 500 μL of methanol was added. After centrifugation at 15,000 rpm for 5 min at 4 °C, the supernatant was removed and the resultant precipitated proteins were dried using a Speed Vac for 2 min. Protein samples were re-dissolved in 10 μL of 8 M urea, 0.1 M Tris–HCl (pH 8.5), and mixed for 1 h. Subsequently, 0.5 μL of 0.1 M DTT was added and incubated at 37 °C for 1 h. Then, 0.5 μL of 208 mM iodoacetamide was added and incubated for 1 h in the dark. After adding 30 μL of 0.1 M Tris–HCl and 60 μL of ultra-pure water, sequencing grade modified trypsin (Promega, Fitchburg, WI, USA) was added, and the solution was incubated at 37 °C for over 15 h.

LC-MS/MS analysis

The tryptic peptides were analyzed with a LTQ Orbitrap mass spectrometer (Thermo Fisher Scientific, Waltham, MA, USA) coupled with a DiNa nanoLC system (KYA Technologies, Tokyo, Japan). Precursor ions were detected over a range of 400–1,500 m/z, and the top four high-intensity ions were selected for MS/MS analyses in a data-dependent mode. Acquired MS/MS spectra were subjected to a database search against the protein sequence database translated from the transcriptome data from the mantle tissues of Laqueus rubellus with the SEQUEST program using Proteome Discoverer software version 1.2 (Thermo Fisher Scientific). The parameters was set as below: the charge state of the precursor ions: automatically recognized, the mass range of tryptic peptides: 800 to 4,500, mass tolerances for precursor ions: 10 ppm mass tolerances for fragment ions: 1 Da. Up to two missed cleavages and modifications of carbamidomethylation (+57.021) of cysteine and oxidation (+15.995) of methionine were considered for calculation of the theoretical masses. False discovery rates (FDRs) was calculated based on a decoy database and using the Proteome Discoverer software. A list of the identified peptides that include a <1 % false discovery rate was obtained after filtering low confidence identification.

Sequence analyses

We performed BlastP similarity searches using the non-redundant protein sequence data stored in GenBank ( using the default settings. The e-value cutoff was set at 1e-10. The domains in the protein sequences were predicted using the SMART program [37, 38] (, including the optional searches for outlier homologues and homologues of known structure, Pfam domains and signal peptides, and the NCBI conserved domain search with the default settings [3941] ( Isoelectric points and molecular masses of the predicted proteins were calculated using Genetyx version 6 (Genetyx, Tokyo, Japan). Sequence alignments were performed using Genetyx version 6.

RT-PCR for gene expression analysis

The mantle and lophophore tissues from a single individual were homogenized in 500 μL of Isogen (Nippon Gene, Tokyo, Japan) in a 1.5-mL microcentrifuge tube. Total RNA was extracted following the manufacturer’s protocol, and treated with RQ1 RNase-Free DNase (Promega, Fitchburg, WI, USA). Reverse transcription was catalyzed by ReverTra Ace (Toyobo, Osaka, Japan), primed with a random primer. The amount of resultant cDNA was quantified by Qubit 2.0 Fluorometer (Life Technologies, Carlsbad, CA, USA), and an amount of 1 ng each of cDNA was used as template for polymerase chain reaction (PCR). PCR was catalyzed by Ex Taq Hot Start Version (Takara, Otsu, Japan). Partial sequence of ICP-1was amplified using the primer pair of Lr46S-1 (5’-GGC CAC ACC TCT GAT GGA TCA T) and Lr46A-1 (5’-TAC ACA CTT AAT GGA GAC CAG GC), and the annealing temperature was set at 58 °C. The following primer pairs were used for amplification of EF-1α; EF-B (5’-CCN CCD ATY TTR TAN ACR TCY TG) and EF-3 (5’-GGN CAY MGN GAY TTY RTN AAR AAY ATG AT), and the annealing temperature was set at 50 °C. Size and amount of RT-PCR products were verified by 1.5 % agarose gel electrophoresis.


  1. 1.

    Murdock DJ, Donoghue PC. Evolutionary origins of animal skeletal biomineralization. Cells Tissues Organs. 2011;194:98–102.

    Article  PubMed  Google Scholar 

  2. 2.

    Jackson DJ, Macis L, Reitner J, Degnan BM, Wörheide G. Sponge paleogenomics reveals an ancient role for carbonic anhydrase inskeletogenesis. Science. 2007;316:1893–5.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Marin F, Luquet G, Marie B, Medakovic D. Molluscan shell proteins: primary structure, origin, and evolution. Curr Top Dev Biol. 2008;80:209–76.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Mann K, Wilt FH, Poustka AJ. Proteomic analysis of sea urchin (Strongylocentrotus purpuratus) spicule matrix. Proteome Sci. 2010;8:33.

    PubMed Central  Article  PubMed  Google Scholar 

  5. 5.

    Joubert C, Piquemal D, Marie B, Manchon L, Pierrat F, Zanella-Cleon I, et al. Transcriptome and proteome analysis of Pinctada margaritifera calcifying mantle and shell: focus on biomineralization. BMC Genomics. 2010;11:613.

    PubMed Central  Article  PubMed  Google Scholar 

  6. 6.

    Marie B, Joubert C, Tayale A, Zanella-Cleon I, Belliard C, Piquemal D, et al. Different secretory repertoires control the biomineralization processes of prism and nacre deposition of the pearl oyster shell. Proc Natl Acad Sci U S A. 2012;109:20986–91.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  7. 7.

    Marie B, Jackson DJ, Ramos-Silva P, Zanella-Cleon I, Guichard N, Marin F. The shell-forming proteome of Lottia gigantea reveals both deep conservations and lineage-specific novelties. Febs j. 2013;280:214–32.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Ramos-Silva P, Kaandorp J, Huisman L, Marie B, Zanella-Cleon I, Guichard N, et al. The skeletal proteome of the coral Acropora millepora: the evolution of calcification by co-option and domain shuffling. Mol Biol Evol. 2013;30:2099–112.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  9. 9.

    Drake JL, Mass T, Haramaty L, Zelzion E, Bhattacharya D, Falkowski PG. Proteomic analysis of skeletal organic matrix from the stony coral Stylophora pistillata. Proc Natl Acad Sci U S A. 2013;110:3788–93.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  10. 10.

    Cusack M, Walton D, Curry G. Shell biochemistry. In: Kaesler R, editor. Treatise on Invertebrate Paleontology, Part H Brachiopoda (Revised) 1. Boulder: Geological Society of America Inc., and The University of Kansas; 1997. p. 243–66.

  11. 11.

    Cusack M, Williams A. Biochemistry and Diversity of Brachiopod Shells. In: Selden P, editor. Treatise on Invertebrate Paleontology, Part H Brachiopoda (Revised) 6. Boulder: Geological Society of America Inc., and The University of Kansas; 2007. p. 2373–95.

  12. 12.

    Collins MJ, Muyzer G, Curry GB, Sandberg P, Westbroek P. Macromolecules in brachiopod shells - characterization and diagenesis. Lethaia. 1991;24:387–97.

    Article  Google Scholar 

  13. 13.

    Jope M. The protein of brachiopod shell—I. Amino acid composition and implied protein taxonomy. Comp Biochem Physiol. 1967;20:593–600.

    CAS  Article  Google Scholar 

  14. 14.

    Walton D, Cusack M, Curry GB. Implications of the amino-acid-composition of recent new-zealand brachiopods. Palaeontology. 1993;36:883–96.

    Google Scholar 

  15. 15.

    Collins MJ, Muyzer G, Westbroek P, Curry GB, Sandberg PA, Xu SJ, et al. Preservation of fossil biopolymeric structures - conclusive immunological evidence. Geochim Cosmochim Acta. 1991;55:2253–7.

    CAS  Article  Google Scholar 

  16. 16.

    Endo K, Walton D, Reyment RA, Curry GB. Fossil intra-crystalline biomolecules of brachiopod shells - diagenesis and preserved geo-biological information. Org Geochem. 1995;23:661–73.

    CAS  Article  Google Scholar 

  17. 17.

    Curry G, Cusack M, Endo K, Walton D, Quinn R. Intracrystalline molecules from brachiopod shells. In: Suga S, Nakahara H, editors. Mechanisms and phylogeny of mineralization in biological systems. Tokyo: Springer; 1991. p. 35–9.

    Google Scholar 

  18. 18.

    Curry GB, Cusack M, Walton D, Endo K, Clegg H, Abbott G, et al. Biogeochemistry of brachiopod intracrystalline molecules. Philos Trans R Soc Lond Ser B Biol Sci. 1991;333:359–66.

    CAS  Article  Google Scholar 

  19. 19.

    Cusack M, Curry G, Clegg H, Abbott G. An intracrystalline chromoprotein from red brachiopod shells: implications for the process of biomineralization. Comp Biochem Physiol B. 1992;102:93–5.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Anstrom JA, Chin JE, Leaf DS, Parks AL, Raff RA. Localization and expression of msp130, a primary mesenchyme lineage-specific cell surface protein in the sea urchin embryo. Development. 1987;101:255–65.

    CAS  PubMed  Google Scholar 

  21. 21.

    Leaf DS, Anstrom JA, Chin JE, Harkey MA, Showman RM, Raff RA. Antibodies to a fusion protein identify a cDNA clone encoding msp130, a primary mesenchyme-specific cell surface protein of the sea urchin embryo. Dev Biol. 1987;121:29–40.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Kito K, Ito T. Mass spectrometry-based approaches toward absolute quantitative proteomics. Curr Genomics. 2008;9:263–74.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  23. 23.

    Ishihama Y, Oda Y, Tabata T, Sato T, Nagasu T, Rappsilber J, et al. Exponentially modified protein abundance index (emPAI) for estimation of absolute protein amount in proteomics by the number of sequenced peptides per protein. Mol Cell Proteomics. 2005;4:1265–72.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Zybailov B, Mosley AL, Sardiu ME, Coleman MK, Florens L, Washburn MP. Statistical analysis of membrane proteome expression changes in Saccharomyces cerevisiae. J Proteome Res. 2006;5:2339–47.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Lu P, Vogel C, Wang R, Yao X, Marcotte EM. Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat Biotechnol. 2007;25:117–24.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Mann K, Poustka AJ, Mann M. In-depth, high-accuracy proteomics of sea urchin tooth organic matrix. Proteome Sci. 2008;6:33.

  27. 27.

    Mann K, Poustka AJ, Mann M. The sea urchin (Strongylocentrotus purpuratus) test and spine proteomes. Proteome Sci. 2008;6:22.

  28. 28.

    Marie B, Zanella-Cleon I, Guichard N, Becchi M, Marin F. Novel proteins from the calcifying shell matrix of the Pacific oyster Crassostrea gigas. Mar Biotechnol (N Y). 2011;13:1159–68.

    CAS  Article  Google Scholar 

  29. 29.

    Ettensohn CA. Horizontal transfer of the msp130 gene supported the evolution of metazoan biomineralization. Evol Dev. 2014;16:139–48.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Parr BA, Parks AL, Raff RA. Promoter structure and protein sequence of msp130, a lipid-anchored sea urchin glycoprotein. J Biol Chem. 1990;265:1408–13.

    CAS  PubMed  Google Scholar 

  31. 31.

    Bedouet L, Duplat D, Marie A, Dubost L, Berland S, Rousseau M, et al. Heterogeneity of proteinase inhibitors in the water-soluble organic matrix from the oyster nacre. Mar Biotechnol (N Y). 2007;9:437–49.

    CAS  Article  Google Scholar 

  32. 32.

    Marie B, Zanella-Cleon I, Le Roy N, Becchi M, Luquet G, Marin F. Proteomic analysis of the acid-soluble nacre matrix of the bivalve Unio pictorum: detection of novel carbonic anhydrase and putative protease inhibitor proteins. Chembiochem. 2010;11:2138–47.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Weiner S, Hood L. Soluble protein of the organic matrix of mollusk shells: a potential template for shell formation. Science. 1975;190:987–9.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Miyamoto H, Miyoshi F, Kohno J. The carbonic anhydrase domain protein nacrein is expressed in the epithelial cells of the mantle and acts as a negative regulator in calcification in the mollusc Pinctada fucata. Zoolog Sci. 2005;22:311–5.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Skovsted CB, Balthasar U, Brock GA, Paterson JR. The tommotiid Camenella reticulosa from the early Cambrian of South Australia: morphology, scleritome reconstruction, and phylogeny. Acta Palaeontol Pol. 2009;54:525–40.

    Article  Google Scholar 

  36. 36.

    Laemmli UK. Cleavage of structural proteins during the assembly of the head of bacteriophage T4. Nature. 1970;227:680–5.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Letunic I, Doerks T, Bork P. SMART 7: recent updates to the protein domain annotation resource. Nucleic Acids Res. 2012;40:D302–5.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  38. 38.

    Schultz J, Milpetz F, Bork P, Ponting CP. SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci U S A. 1998;95:5857–64.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  39. 39.

    Marchler-Bauer A, Bryant SH. CD-Search: protein domain annotations on the fly. Nucleic Acids Res. 2004;32:W327–31.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  40. 40.

    Marchler-Bauer A, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, et al. CDD: specific functional annotation with the Conserved Domain Database. Nucleic Acids Res. 2009;37:D205–10.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  41. 41.

    Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, et al. CDD: a Conserved Domain Database for the functional annotation of proteins. Nucleic Acids Res. 2011;39:D225–9.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

Download references


We thank the staff of the Misaki Marine Biological Station (The University of Tokyo) for collecting brachiopod samples. We also thank Elika Iioka and Misa Kiuchi (The University of Tokyo) for DNA sequencing. This study was supported by the Japan Society for the Promotion of Science Grant-in-Aid for Scientific Research No. 15104009.

Author information



Corresponding author

Correspondence to Kazuyoshi Endo.

Additional information

Competing interests

The authors declared that they have no competing interests.

Authors’ contributions

YI collected Laqueus rubellus samples, and curried out RNA extraction, shell matrix protein extraction, SDS-PAGE analysis, sequence analysis, and drafted the manuscript. IS participated in extraction of shell matrix proteins and SDS-PAGE analysis. KK carried out LC-MS/MS analysis. KO and MH carried out transcriptome analysis. KE conceived of the study, participated in its design and coordination, and helped to draft the manuscript. All authors read and approved the final manuscript.

Additional files

Additional file 1:

RT-PCR analysis of ICP-1 and EF-1α. M: mantle; L: lophophore. (TIFF 897 kb)

Additional file 2:

The depth of the isotigs of shell matrix proteins in this study. (TIFF 597 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Isowa, Y., Sarashina, I., Oshima, K. et al. Proteome analysis of shell matrix proteins in the brachiopod Laqueus rubellus . Proteome Sci 13, 21 (2015).

Download citation


  • Brachiopoda
  • Biomineralization
  • Shell matrix protein
  • Proteome
  • Transcriptome