Skip to main content
  • Research article
  • Open access
  • Published:

The skeletome of the red coral Corallium rubrum indicates an independent evolution of biomineralization process in octocorals

Abstract

Background

The process of calcium carbonate biomineralization has arisen multiple times during metazoan evolution. In the phylum Cnidaria, biomineralization has mostly been studied in the subclass Hexacorallia (i.e. stony corals) in comparison to the subclass Octocorallia (i.e. red corals); the two diverged approximately 600 million years ago. The precious Mediterranean red coral, Corallium rubrum, is an octocorallian species, which produces two distinct high-magnesium calcite biominerals, the axial skeleton and the sclerites. In order to gain insight into the red coral biomineralization process and cnidarian biomineralization evolution, we studied the protein repertoire forming the organic matrix (OM) of its two biominerals.

Results

We combined High-Resolution Mass Spectrometry and transcriptome analysis to study the OM composition of the axial skeleton and the sclerites. We identified a total of 102 OM proteins, 52 are found in the two red coral biominerals with scleritin being the most abundant protein in each fraction. Contrary to reef building corals, the red coral organic matrix possesses a large number of collagen-like proteins. Agrin-like glycoproteins and proteins with sugar-binding domains are also predominant. Twenty-seven and 23 proteins were uniquely assigned to the axial skeleton and the sclerites, respectively. The inferred regulatory function of these OM proteins suggests that the difference between the two biominerals is due to the modeling of the matrix network, rather than the presence of specific structural components. At least one OM component could have been horizontally transferred from prokaryotes early during Octocorallia evolution.

Conclusion

Our results suggest that calcification of the red coral axial skeleton likely represents a secondary calcification of an ancestral gorgonian horny axis. In addition, the comparison with stony coral skeletomes highlighted the low proportion of similar proteins between the biomineral OMs of hexacorallian and octocorallian corals, suggesting an independent acquisition of calcification in anthozoans.

Background

Biomineralization is a widespread process in metazoans, and biomineral structures provide a multitude of functions (protection, maintenance, etc.) [1]. Among marine calcifiers, corals represent one major group. They produce a CaCO3 exoskeleton to sustain the vertical and horizontal growth of the coral colony. The vernacular term “coral” refers to calcifying organisms within the phylum Cnidaria. Anthozoa, a class of Cnidaria, encompasses two major subclasses (Fig. 1): Hexacorallia and Octocorallia. Among hexacorallians only stony corals (order Scleractinia) are calcifiers, they are the basis of the edification of coral reefs [2]. On the other hand, octocorallians, with very few exceptions, are all calcifiers in the sense that they produce sclerites, i.e. micrometric calcified skeletal elements found throughout the tissues [3]. Many octocorallians, such as gorgonians (order Alcyonacea), also have a central horny axis composed of sclerotized collagen [4, 5]. In red corals (order Alcyonacea), the central axis is a calcified axis [6, 7]. Hexacorallians and octocorallians diverged approximately 600 Mya [8, 9]. Comparison of the mechanisms that govern the production of calcium carbonate (CaCO3) biominerals in hexacorallians (scleractinians) and octocorallians is thus expected to shed light on the conservation of the protein toolkit used by these organisms for making their biocalcification. Moreover, studying corals that can produce distinct skeletal structures holds the potential to provide novel insight into the evolutionary aspect of biomineralization.

Fig. 1
figure 1

Phylogenetic relationship of anthozoan species used in the present study. Anthozoans are divided into two subclasses: Octocorallia and Hexacorallia. In Octocorallia, 3 species are represented: the red coral Corallium rubrum (order Alcyonacea), which is the model of the present study (red frame), and 2 other species, the soft coral Dendronephthya gigantea (order Alcyonacea) and the blue coral Heliopora coerulea (order Helioporacea). In Hexacorallia, 4 species are represented: 2 stony corals (order Scleractinia) and 2 sea anemones (order Actinaria). Phylogeny is based on www.tolweb.fr. The presence of axial skeleton and/or sclerites and the CaCO3 polymorph are indicated by colored stars

In this respect, the red coral Corallium rubrum is a model of interest for the study of biomineralization. This species, endemic to the Mediterranean Sea, is noteworthy for its intense red color but also because it produces two distinct biominerals, the axial skeleton and the sclerites. The sclerites are formed by specialized cells named scleroblasts and are spread within the mesoglea [10, 11]. Axial skeleton formation is the result of two separated processes that produce two distinct structures, the medullar and the annular regions. First, the axial skeleton extends at the apex of the branches by an aggregation of the sclerites (longitudinal branch extension), forming the medullar part [7, 12,13,14]. Then, calcite is deposited concentrically around the medulla by the axial skeleton epithelium. While the axial skeleton has a clear role in supporting the coral colony, the sclerites have suggested roles in protection against predators, waste storage, and as initiation sites for axial skeleton formation [14]. Both biominerals are composed of a high-magnesium calcite polymorph as opposed to the aragonite skeletons of stony corals [7, 15,16,17,18]. In C. rubrum, the organic matrix (OM; 1–2% of the biomineral), which can be dissociated as a water-soluble organic matrix (SOM) and water-insoluble organic matrix (IOM), essentially contains proteins, glycoproteins and polysaccharides [19]. The term “organic matrix” refers here to the components of the matrix forming the scaffold involved in the process of biomineralization, i.e. the structural components displaying extracellular matrix (ECM) properties (scaffold), as well as the components required for the precipitation of CaCO3 [20, 21]. In cnidarians, the process of OM assembly remains to be elucidated.

Our understanding of the role of OM components in biomineralization is derived from many in vitro and a few in vivo experiments. These molecules can stabilize amorphous calcium carbonate and control the nucleation, the orientation and the polymorph selection of the crystal [20, 22,23,24,25]. So far organic matrices of C. rubrum have been electrophoretically characterized and hitherto only one OM protein, scleritin, has been sequenced [13, 19, 26,27,28]. Recently, the first investigation of the OM proteome of the two biominerals was performed in the octocorallian species Corallium konojoi (Japanese red and pink coral species). The authors showed a low similarity of proteins between the SOM of the axial skeleton and the sclerites, with only 9 out of the 147 identified proteins being common to both biominerals. However, at the time of this study, neither octocorallian transcriptomes nor genomes were available, and the authors had to rely on phylogenetically distant hexacorallian transcriptomes to identify their octocorallian OM proteins [29]. Hence, identification of the proteome components may have been artefactual in parts.

In the present study, we performed high-resolution mass spectrometry analysis to characterize the proteomes of the SOM and the IOM of the axial skeleton, and the sclerites of C. rubrum. By analyzing these results in combination with the C. rubrum transcriptome, we have identified a total of 102 OM proteins that were further manually curated and categorized according to their possible function. Also, we analyzed their conservation with other sequenced anthozoans: (i) 3 octocorals, (ii) 4 scleractinian species known to have a fast skeletal growth rate, and (iii) 3 non-calcifying sea anemones.

Results

The proteome of Corallium rubrum biominerals

For each biomineral, the extraction process produced 2 OM fractions, the water-soluble organic matrix (SOM) and the water-insoluble organic matrix (IOM). After 1D SDS-PAGE electrophoresis and trypsin digestion of the SOM and IOM of the two biominerals, spectra from tandem mass spectrometry were analyzed and interrogated against our in-house transcriptome assembly (CrubmRNA.fasta; data available at http://data.centrescientifique.mc/CSMdata-crubrum_data.html) using Mascot. We identified and retained for further analysis 107 proteins with a minimum of 2 unique peptides per protein sequence in at least one of the 4 fractions.

The 107 protein sequences were then annotated using the Blast2GO pipeline as well as other prediction tools (see Methods). They were named CR_n, in the descending order of their unique peptide count. According to our BLAST analysis, out of the 107 proteins, 34 were distributed into 9 groups of similar proteins (Additional file 1). Two groups were particularly represented within the red coral skeletome; the collagen-like proteins (10 proteins), and the agrin-like proteins (8 proteins). Both families are known as extracellular matrix (ECM) proteins [30, 31]. Also of interest, 20 proteins did not match with known proteins (no BLAST hit). Finally, when analyzed by a BLAST search against NCBI Metazoa and Prokaryota RefSeq databases (April 2020), one protein (CR_22) exhibited high similarity with prokaryote proteins, and with uncharacterized proteins from only one metazoan species: Dendronephthya gigantea (Octocorallia/Alcyonacea). Concurrent to the proteomic and transcriptomic project of the red coral, we are currently assembling the C. rubrum genome (PG, MA, DA, ST; personal communication). Although at a draft stage, we were able to map all the aforementioned 21 (20 + 1) proteins’ mRNA sequences to our C. rubrum genome assembly. Beside the fact that all proteins were encoded in the genome (Table 1), they also contained introns (except CR_91), suggesting that they were genuine C. rubrum proteins.

Table 1 Corallium rubrum proteins identified in the axial skeleton and the sclerites

To avoid redundancy, we manually cross-examined the sequences contained within each group/family against our C. rubrum genome draft assembly (PG, MA, DA, ST; personal communication) and another C. rubrum published transcriptome [32]. Several of the initially assembled transcripts turned out to be partial and were in fact portions of the same larger sequence. Indeed, CR_3 and CR_4 are 2 moieties of the same CAP sequence (cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 protein) and the three CR_10, CR_28, CR_60 and the three CR_62, CR_76, CR_82 were portions of two large agrin-like proteins (Additional file 1, Additional file 2). The non-redundant protein set of the C. rubrum OM is thus 102 proteins (Table 1, Additional file 3). The rest of the analysis was performed on these 102 proteins.

Sclerites vs axial skeleton organic matrix composition

The axial skeleton is formed by medullar part (occluding sclerites) and annular part (without sclerites). Relative contribution of medullar and annular parts in the axis is dependent on the colony age, since the medulla is of constant diameter whereas the annular part increases throughout time. Using 7 published axial skeleton cross-sections (Additional file 4) [7, 33], we estimated the proportion of the sclerites into the axial skeleton. The surface ratio between the medulla and the annula (easily discernible after staining, see references) showed that the medulla surface represented 5.37% (SEM ± 1.05%) of the annula (Additional file 4). In medulla, sclerites make up about half of the medullar area [14], which corresponds to a proportion of about 3% of the axis surface. Thus, in our study, the axial skeleton also contained proteins from the sclerite OM, although the overall proportion is low.

We explored the distribution of OM proteins within the axial skeleton and the sclerites (Fig. 2a, Table 1). Among the 102 proteins composing the red coral proteome, the OM of the axial skeleton and the sclerites was composed of 79 and 75 proteins, respectively. Twenty-seven were specific to the axial skeleton whereas 23 were specific to the sclerites; 52 were common to both. These numbers correspond to proteins identified in individual fractions with a cut-off of minimum 2 peptides per protein. With a cut-off decreased to minimum of 1 peptide per protein, the 102 OM proteins would correspond to 92 and 84 proteins present in the axial skeleton and the sclerites, respectively, with 74 proteins common to both biominerals (Additional file 3 and Additional file 5). Nevertheless, we here favored the stringent cut-off (2 peptides/protein) in the rest of the analysis in order to avoid false positive protein identification in our samples [34].

Fig. 2
figure 2

Distribution of the proteins identified in the organic matrix of the red coral axial skeleton and sclerites. a Among the total 102 proteins identified in the C. rubrum biominerals, 79 are present in the axial skeleton (left) and 75 are present in the sclerites (right), 52 being shared between the two. The Venn diagram shows the number of shared and non-shared proteins between the 4 extracted organic matrices: SOMax: soluble organic matrix of the axial skeleton, IOMax: insoluble matrix of the axial skeleton, SOMsc: soluble organic matrix of the sclerites, and IOMsc: insoluble organic matrix of the sclerites. b Protein domains identified in the proteome of the red coral biominerals and their distribution using circle representation in SOMax (dark blue), IOMax (blue), SOMsc (yellow) and IOMsc (red). ADAM10 desintegrin and metalloproteinase, AMOP adhesion-associated domain present in MUC4 and other proteins, BPI bactericidal permeability-increasing protein/lipopolysaccharide-binding protein/cholesteryl ester transfer protein N-terminal domain, CA carbonic anhydrase, EGF-like epidermal growth factor, EF-hand a calcium-binding domain, F5/8 type C coagulation factor 5/8 C-terminal domain, FAS1 fasciclin-like, HYR hyalin repeat, IG-like immunoglobulin-like, KAZAL serine protease inhibitor, LamG laminin G, LCD low complexity domain, NIDO extracellular domain of unknown function in nidogen, PAM peptidylglycine α-amidating monooxygenase, PHM peptidylglycine α-hydroxylating monooxygenase, ShK toxic stichodactyla toxin, SCP Cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 proteins (CAP) superfamily proteins, TIMP tissue inhibitor of metalloproteinase, vWA von Willebrand type A, vWD von Willebrand type D, UBA52 Ubiquitin A-52, UOMPs uncharacterized organic matrix proteins, WAP whey acidic protein ‘four-disulfide core’

The protein distribution between the SOM and IOM fractions of the organic matrix was somewhat equivalent in both the axial skeleton (65 SOM and 64 IOM proteins) and the sclerites (62 SOM and 58 IOM proteins, Fig. 2a). In the axial skeleton and the sclerites, only minorities of the proteins were specific to the IOM (10 and 9, respectively) or the SOM (10 and 10, respectively), and 34 were common to the SOM and IOM of both C. rubrum biominerals.

Using the prediction tools InterProScan and SMART, we analyzed the domain composition of proteins of the C. rubrum biominerals. We found domains such as vWA, vWD, EGF-like, thrombospondin, sushi, reeler, fibrinogen, fibronectin and protocadherin (Fig. 2b). We also noticed the presence of sugar-binding domains in 10 proteins (CR-21, 22, 45, 47, 57, 58, 65, 70, 85 and 97; Table 1). Interestingly, the fibronectin (CR_84), protocadherin (CR_24 and CR_34), and the peroxidase (CR_92) domains were found in proteins present in the axial skeleton but not in the sclerites. Finally, for 18 of the OM proteins, no domains were identified and they were classified as Uncharacterized Organic Matrix Proteins (UOMPs; Table 1).

We quantified the relative abundance of proteins (Additional file 3 sheet#1) in the four OM fractions using two label-free methods: exponentially modified protein abundance index (emPAI) and weighted spectral count (WS). According to emPAI values, the 4 proteins CR_1, CR_3-4, CR_5 and CR_10-28–60 were listed within the top 10 most abundant proteins in all 4 samples (i.e. SOM and IOM in both the axial skeleton and the sclerites; Additional file 6). We noticed that CR_1 (scleritin), previously identified in the OM of the sclerites of C. rubrum [13], was the highest abundant protein in SOM and IOM of the axial skeleton (emPAI: 32.011 and 30.225, respectively) and the sclerites (emPAI: 71.871 and 66.78, respectively), which is supported by WS results (Additional file 6). Although of unknown function and with no functional domain annotation, its high abundance in both biominerals emphases its importance in the biomineral formation of C. rubrum. CR_1 was followed distantly by CR_3-4 (CAP family protein, a prototype extracellular protein) that was more abundant in the axial skeleton OM than in the sclerites OM.

Major protein families of the red coral OM

Among the 102 proteins found in the OM of C. rubrum biominerals, we choose to limit in this section those with potentially evolutionary interest.

Collagen-like proteins

In the red coral, out of 10 collagen-like proteins, 8 were identified in the two biominerals (CR_5, 9, 17, 35, 38, 39, 50, 54, 56 and 68), one in the axial skeleton (CR_35) and one in the sclerites (CR_68) (Fig. 3a, Additional file 7). Collagens are ECM proteins that constitute a large family involved in the ultrastructural organization of the organic matrix network [35]. Most collagens contains a repeat motif of 3 amino acids (G–X–Y) and form fibrils after post-translation modification [36]. Initially, 12 proteins were similar to collagens according to the BLAST search results. However, 2 of these candidates, CR_6 and CR_98, were arguably assigned as collagens: they did not possess triplet repeats G–X–Y, but harbored fibrinogen-related and von Willebrand A domains, which are also present within collagen sequences of other species. As a consequence, they were automatically annotated as “collagen” (Additional file 8), but here, we did not consider these as such. The same holds true for two other proteins (P14 and P18) identified as collagen in the OM of the scleractinian coral S. pistillata [18]. However, with the completion of the S. pistillata genome [37], reassessment of these sequences shows that P14 is a coadhesin, and P18 a vWA domain-containing protein (Additional file 8). Finally, only one collagen-like sequence was identified in the biomineral of another scleractinian coral, Acropora millepora, which corresponds to a collagen containing G–X–Y motif (Acc. B8V7R6; [16]). In C. rubrum, the remaining 10 collagen-like proteins (CR_5, 9, 17, 35, 38, 39, 50, 54, 56 and 68) had typical collagen domains (Fig. 3a, Additional file 7). Five of them also had a C-terminal fibrillar collagen domain. The 5 without (CR_5, 38, 39 and CR_9, 35) were encoded on two gene clusters (Fig. 3b). Taken individually, none of these 10 collagen-like proteins accounted for the major protein of a biomineral fraction. However, when summed up (Fig. 3c), “collagen” was the most abundant protein of the axial skeleton OM (~ 14%). In the sclerites, collagen-like proteins were present at similar quantities (~ 10%), though scleritin was by far the most abundant OM protein (~ 40%). In vertebrates, collagens are classified in types according to their physico-chemical properties and tissue specific expression. Based on BLAST searches and domain signatures, the C. rubrum collagen-like proteins were α-type fibrillar assigned as type I, II, III, V, VI, and XI (Table 1). Of note, collagen type IV, which is characterized by its NC1 domain, is present in cnidarians [38] including C. rubrum (P. Ganot, personal communication), but was not detected in the biominerals’ OM.

Fig. 3
figure 3

Schematic representation of the 10 collagen-like protein sequences identified in the proteome of the red coral biominerals. The collagen sequence identified in the Acropora millepora (UniProt accession number B8V7R6) organic matrix is included. a Schematic protein domain (described in the box below) representation of the different collagen proteins with their respective ID (CR_n). Colored circles indicate presence of the protein in different fractions as described in Fig. 2b. EGF-like epidermal growth factor-like, SP signal peptide, TSPN Thrombospondin N-terminal-like domain, vWA von Willebrand factor type A. Scale bars = 100 amino acids. b Genomic synteny of the CR_5, CR_39, and CR_38 genes and of the CR_9 and CR_35 genes; exons are in green boxes, introns are in blue lines, and intergenic sequences are in dotted lines. Scale bar = 1000 base pairs. c The overall proportion of collagen-like proteins and the most abundant protein, scleritin (CR_1), in each sample in percentage of weighted spectra

Agrin-like proteins

The agrin-like protein family is the second most diversified family, in terms of number of genes, in the C. rubrum OM. In vertebrates, agrin is a single copy gene that produces different proteins through alternative splicing [39, 40]. Agrins are long proteoglycans (heavily glycosylated proteins) found in many extracellular environments. They are composed of two main parts, the N-terminal moiety, which is composed of a series of Foln/Kazal domains (Fig. 4a) where glycosaminoglycan (GAG) side chains are attached, and the C-terminal moiety, which contains binding sites for various ligands including LRP4 (Low-density Lipoprotein Receptor-related Protein 4), and dystroglycan, among others [39, 40]. Agrins are highly repetitive proteins, especially in their N-terminal moiety, and as such, their nucleotide sequence is difficult to assemble, especially with high-throughput RNAseq technologies using short reads. In our C. rubrum transcriptome assembly, we found different transcripts corresponding to pieces of agrin homologs (CR_10, 19, 28, 53, 60, 62, 76 and 82). However, when mapped to our draft genome and to the transcriptome of Pratlong [32], we found that CR_10, CR_60, and CR_28, and CR_76, CR_82 and CR_62 were parts of two agrin proteins, agrin-1 and agrin-2, respectively (Fig. 4a, Additional file 1 and Additional file 2). In all, we were able to identify 4 different agrin homologs in the OM of the C. rubrum biominerals. CR_10-28-60, CR_62-76-82 and CR_19 were present in both biominerals and CR_53 was only identified in the sclerite OM (Fig. 4a, Additional file 3 sheet#2). Although their protein sequence may yet be partial for some of them, all 4 corresponded to homologs of the N-terminal moiety of vertebrate agrins, i.e. the moiety carrying the sugar chains.

Fig. 4
figure 4

Schematic representation of the 7 proteases and the 4 agrin-like proteins identified in the proteome of the red coral biominerals. a The peptidase S9 (CR_90) and the matrix metallopeptidase 13 (MMP13, CR_61) are specific to the axial skeleton OM. The cathepsin Z (CR_81) and one MMP12A (CR_64) are specific to the sclerites. The peptidase S8 (CR_30), the other MMP12A (CR_16) and the MMP12B (CR_20) are found in both biominerals. b CR_10, CR_60 and CR_28 are part of the same agrin-1 protein. CR_76, CR_82 and CR_62 are part of the same agrin-2 protein. Proteome ID are indicated, circle representation in the different biominerals fractions is as in Fig. 2b, and protein domains are described in the box. C345C netrin C-terminal domain, DISIN disintegrin domain, EGF-like epidermal growth factor-like, F5/8 type C coagulation factor 5/8 C-terminal domain, FOLN follistatin-N-terminal domain-like, LDLa low-density lipoprotein receptor domain class A, ShK toxin stichodactyla toxin, SP signal peptide, TM transmembrane, TY thyroglobulin type I repeats, WAP four-disulfide core domain

Matrix peptidases

We identified 7 peptidases in the C. rubrum OM: 4 matrix metallopeptidases (MMPs), 2 serine peptidases and 1 cysteine peptidase (Table 1, Fig. 4b). These proteases are not equally distributed between the 2 biominerals: the serine peptidases S9 (CR_90) and MMP13 (CR_61) were specific to the axial skeleton OM whereas the cysteine protease (CR_81) and one MMP12 (CR_64) were specific to the sclerite OM (the remaining 3 being present in both biominerals; Additional file 3 sheet#2).

Carbohydrate-binding proteins

Among the striking features of the C. rubrum OM composition is the large number of proteins containing sugar-binding domains, most related to lectins (Fig. 2b; Additional file 3 sheet#2). In C. rubrum, 10 extracellular proteins contained carbohydrate-binding domains. CR_21, 22, 45, 47, 57 and 58 were found in both sclerites and axial skeleton OM whereas CR_65, 70, 85 and 97 were only found in the axial skeleton OM. BLAST search against NCBI_nr database revealed that CR_22 was very similar to proteins in the two octocorallians Dendronephthya gigantea and Heliopora coerulea, but not to other metazoan species (Fig. 5, Additional file 9). CR_22 was also similar to proteins in prokaryotes, with first hits in Myxococcus and Actinoplanes (Additional file 9). In order to conduct a phylogenetic analysis of the putative CR_22 homologs, we recurrently retrieved the first hits (regardless of the e-value) from different BLAST searches conducted against specific subdivisions of the metazoan tree (Additional file 9). PhyML analysis using the protein sequences from these BLAST results demonstrated that CR_22 had homologs in Octocorallia and Prokaryota, but not in other metazoan lineages, including Hexacorallia (Additional file 9). Importantly, CR_22 and its D. gigantea homolog are both encoded in their respective octocorallian genome; their genes contain introns with similar exon/intron organization (Additional file 9). Domain analyses indicated that CR_22 (and its homologs) encompasses a duplicate carbohydrate-binding domain family CBM25 [41], although with a low Pfam e-value. Thus, a plausible explanation is that the sugar-binding protein CR_22 could have been horizontally acquired through a prokaryote gene transfer early during Octocorallia evolution.

Fig. 5
figure 5

Conservation of the C. rubrum biomineral proteomes with other anthozoan species. Venn diagram depicts the conservation of the 102 C. rubrum OM proteins within anthozoan species. Among the 102 proteins, 77 share sequence similarity with hexacorallian (Scleractinia and Actiniaria) and octocorallian proteins, 1 (CR_72) is shared between Actiniaria and Octocorallia, and 24 are only found in octocorallians. Among the latter 24, 18 are shared in octocorallians and 6 are specific to C. rubrum. The 18 octocorallian proteins are listed and the 6 C. rubrum specific proteins are detailed as in Fig. 3b (pI: isoelectric point). Sequence names in blue correspond to proteins remarkable for their acidic pI and number of negative charges at pH = 8; names in red correspond to protein remarkable for their alkaline pI and number of positive charges at pH = 8 (see Additional file 3 sheet#2)

Carbonic anhydrases

Carbonic anhydrases (CAs) are metallo-enzymes that catalyze the reversible hydration of carbon dioxide into bicarbonate, the source of inorganic carbon for the precipitation of CaCO3. Here, CR_29, which corresponds to CruCA4 previously identified in the C. rubrum tissues [42], was identified with equivalent abundance in the proteome of SOM and IOM of the axial skeleton and the sclerites (emPAI: 1.15, 1.24, 0.64 and 0.71, respectively; Additional file 3 sheet#1).

The C. rubrum OM compared with other cnidarians

In order to better understand the conservation of the C. rubrum OM proteins in regard to other anthozoans, the 102 C. rubrum OM proteins were BLAST searched against a selection of nucleotide databases from 3 anthozoan taxa (3 octocorallians, 3 actiniarians and 4 scleractinians; Additional file 10). Hits were retained for TBLASTN e-values below 1e-15 (no seg filter) and sequences were classified according to their taxon affiliation. Of the 102 OM proteins, 77 were found in at least one database of both subclasses (Octocorallia and Hexacorallia); 1 (CR_72) was found in octocorallians and actiniarians (not in scleractinians); 18 were shared only within the 4 octocorallian species and 6 were specific to C. rubrum (Fig. 5). Among the latter 6, structural analysis showed that Low Complexity Domains (LCDs) were dominant. LCDs have been described in the secreted OM of biominerals in different metazoan taxa [43,44,45,46]. We identified D-rich (aspartate-rich) LCD in CR_37 and CR_48, E-rich (glutamate-rich) and K-rich (lysine-rich) in CR_106. These D/E-rich domains involve a very acidic pI (Fig. 5, Additional file 3 sheet#2). Other LCDs such as G-rich have been identified in the OM of mollusk shell [47], but these domains were absent from the C. rubrum biomineral proteomes. Among the 24 proteins only identified in octocorallians, 9 were strongly alkaline (with pI > 9; CR_1, CR_21, CR_31, CR_47, CR_75, CR_84, CR_88, CR_91 and CR_97; Additional file 3 sheet#2), including the highly abundant scleritin (CR_1, conserved in the 3 tested octocorallians). Seven out of the 24 proteins (CR_75, 84, 85, 91, 97, 103 and 106; Additional file 3 sheet#2) were detected in the axial skeleton biomineral only.

With regards to the conservation of OM proteins between C. rubrum OM and stony corals, only 11 proteins (CR_2, 3–4, 7, 14, 18, 24, 27, 29, 34, 58 and 68) had homologs with constituents of the OM of the scleractinian biominerals (Fig. 6, Additional file 3 sheet#2 and Additional file 11). Within these 11 proteins, 2 were galaxin-like proteins, 1 was LRP and 2 were protocadherins. The two galaxin-like proteins (CR_14 and CR_27) exhibited 10 and 11 positions with a double di-cysteine motif that we identified using MEME discovery motif (Additional file 12). This motif was similar to the galaxin motif previously identified in the scleractinian coral, A. millepora [48]. In C. rubrum, these two galaxin-like proteins are present in both the axial skeleton and the sclerites. CR_7 (identified in both sclerites and axial skeleton) contains five low-density lipoprotein-receptor YWTD domains, as found in the mammalian LRP [49]. The mucin-like (CR_2) is one of the most abundant proteins. CR_3-4 is a prototype extracellular protein composed of the GAPR-1 (Golgi associated pathogenesis related-1) domain from the CAP superfamily followed by 2 TSP1 domains. Of note, we arguably accounted the cytoplasmic actin (CR_49) and the ubiquitin-60S ribosomal protein L40 (CR_67) as cellular contaminants, although they were also found in the OM of scleractinian corals.

Fig. 6
figure 6

Schematic representation of eleven putatively conserved proteins identified in OM of the red coral and the scleractinian coral biominerals and their percentage of identity and similarity. Names, Proteome_ID and sequence description of the 11 C. rubrum proteins are indicated (left). Their presence in the different biomineral fractions is represented with circle as in Fig. 2b, and protein domains are described in the box. The 11 putative homologs were present in the published skeletome of 3 scleractinian species [16, 18, 21] and the OM of C. rubrum biominerals of our study (see Additional file 3 sheet#2). Percentage of identity and similarity between proteins of C. rubrum and proteins identified in the skeletome of the 3 other scleractinian species (Additional file 11) were calculated using https://www.bioinformatics.org/sms2/ident_sim.html and are indicated in the table (right; N: no putative homolog). For each scleractinian protein, GenBank accession number is indicated and published sequence name is given in brackets. Scale bars = 100 amino acids. AMOP adhesion-associated domain present in MUC4 and other proteins, Cadh cadherin repeats, Cadh_C cadherin_C, CAP cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 protein, CA carbonic anhydrase, EGF epidermal growth factor, HH_Nterm hedgehog, IG-like immunoglobulin-like, IGc2 immunoglobulin C-2 type, LamG laminin G, LCD low complexity domain, LDLa low-density lipoprotein receptor domain class A, LY low-density lipoprotein receptor YWTD domain, NIDO extracellular domain of unknown function in nidogen, SH2 Src homology 2, SP signal peptide, TM transmembrane, Tryp_SPc trypsin-like serine protease, TSP1 thrombospondin-1, vWA von Willebrand type A, vWD von Willebrand type D

Discussion

Organic matrices (OMs) from biominerals are specialized ECMs that support the extracellular process of calcification. Among anthozoans, the OM of coral biominerals is widely studied in scleractinian (Hexacorallia) species [16, 18, 50,51,52,53], but poorly known in Octocorallia. Here, we studied the OM composition of the two high-magnesium calcite biominerals (the axial skeleton and the sclerites) produced by the octorallian species Corallium rubrum. This study aims to compare the protein repertoire (i) in both biominerals of C. rubrum, and (ii) with skeletomes of other anthozoan species. This new proteome provides novel understanding into the evolution of skeletomes in anthozoan biomineralization.

The red coral axial skeleton and sclerite proteomes

The LC–MS/MS analysis of the axial skeleton and the sclerite OMs extracted from C. rubrum allowed us to identify 102 non-redundant proteins (Table 1). Organic matrices of both biominerals contain 52 common proteins, which may be due to the presence of the sclerites in the axial skeleton. In our study, scleritin, which was previously identified only in sclerites [13], is predominant in sclerites and in the axial skeleton. According to our estimates, sclerites only make up a minor component of the axial skeleton (3%). Given this, and the very high abundance of scleritin in the axial skeleton proteome, we suspect that scleritin is also secreted by the axial skeleton epithelium. The two biominerals also present specificities with 27 proteins only identified in the axial skeleton and 23 in the sclerites. Regarding SOM and IOM, the difference in solubility is unknown but may be explained by different degrees of post-translational modifications and/or increasing matrices complexity due to cross-linking with other types of molecules (e.g. polysaccharides). The present study is limiting by the unique peptide threshold (2 unique peptides per protein) and the fact that one replicate per sample was used for mass spectrometry analysis.

Organic matrix formation requires synthesis (secretion), assembly (e.g. cross-linking) and remodeling (e.g. digestion) of the various proteins forming the OM network [54, 55]. In the red coral proteome, we identified characteristic domains of proteins normally found in various ECMs (vWA, vWD, EGF-like, thrombospondin, sushi, reeler, fibrinogen, fibronectin and protocadherin) [56]. Several have matrix-remodeling function such as peptidases and peptidase inhibitors, which constitute a significant amount of the proteins fraction in various metazoan taxa [16, 57, 58]. The red coral skeletome possesses 7 peptidases including 4 matrix metallopeptidases (MMPs). MMPs are calcium-dependent zinc-containing endopeptidases with different substrate specificities central to the modeling of the OMs of mineralized tissues [59]. In mammals, MMP13 is known as the collagenase processing triple helical collagens (type I, II, and III among others) involved in bone remodeling and mineralization processes, whereas MMP12 digests elastin and a number of ECM molecules, and is also found in bone tissues [60,61,62]. The fact that each biomineral contains a different set of peptidases (serine peptidase and MMP13 in the axial skeleton, and cysteine peptidase and MMP12 in the sclerites) implies that their respective OM will be remodeled in different (specific) manners. Likewise, a prolyl 4-hydroxylase (P4H) homolog (CR_93) is found in the OM of the axial skeleton but not in the OM of the sclerites (Additional file 3 sheet#2). P4H is responsible for proline residues hydroxylation [63, 64], and its presence in the axial skeleton may reflect differences in collagen fibers post-translational modifications. The presence of proteins containing fibronectin (CR_84), protocadherin (CR_24 and CR_34) and peroxidase (CR_92) domains in the axial skeleton but not in the sclerites suggests functional differences between the set of OM proteins of each biomineral. In bilaterians, fibronectin binds to various other extracellular proteins [65, 66], protocadherin are trans-membrane adhesion proteins [67] and extracellular haem peroxidases have been involved in extracellular protein cross-linking [68, 69]. Altogether, these results suggest different matrix assembly properties of the axial skeleton OM compared with the OM of the sclerites, despite a similar composition in structural matrix proteins (such as collagen-like and agrin-like proteins).

Evolutionary perspectives on coral calcification

The proteome composition of the red coral OMs exhibits numerous differences compared with stony coral skeletomes [16, 18, 21]. The most striking differences are the abundance of the structural collagen-like and agrin-like proteins, the cohort of protein carrying sugar-binding domains and the octocorallian specific OM proteins such as the predominant scleritin.

In C. rubrum, we identified 8 different collagen-like proteins in both biominerals (CR_5, 9, 17, 38, 39, 50, 54 and 56), 1 in the axial skeleton (CR_35) and 1 in the sclerites (CR_68). After synthesis of fibrillar collagens, the chains are extensively modified (e.g. in conserved triplet repeats G–X–Y, where Y is often proline, many proline residues are hydroxylated), an essential step for glycosylation and tertiary structuration of the collagen fibers [36]. The C. rubrum collagen-like sequences contain G–X–Y repeat motifs suggesting that they can assemble into fibers. Quantitatively, they represent a total 10–15% of the total amount of the OM proteins (Fig. 3c). This is highly contrasting with reef-building corals, where 1 collagen sequence was identified in A. millepora [16]. In gorgonian species, collagens have been found to be a major component of the sclerotized axis [70]. Thus, collagen-rich axis appears as a shared trait with gorgonians.

Then, we also identified 4 agrin-like proteins among which 3 are present in both C. rubrum biominerals. They represent candidate core-proteins for carrying glycoaminoglycans (GAGs), such as heparan sulfate or chondroitin sulfate. The presence of GAGs in the OM of C. rubrum has already been detected [19], but they were not identified in the OM of stony corals [16, 18]. In addition, the presence of an exostosin-3 homolog (CR_41), a protein responsible for the extension of heparan sulfate chains, corroborates the presence of GAGs in the red coral biominerals. The finding of CR_10-28-60 in the top 10 of highest abundant proteins and with equivalent relative abundance in the SOMax, IOMax, SOMsc and IOMsc (Additional file 6), points at this agrin-1 proteoglycan as an important structural component of the OM of both biominerals.

Lectins (carbohydrate-binding proteins) are widespread in animals [71,72,73], and have been described in mollusk shells, echinoderms skeletons and bird eggshells [74,75,76]. Polysaccharides make an important part of the skeletome of hexacorallian and octocorallian corals [16, 19]. The presence of at least 10 proteins containing sugar-binding domains suggests manifold possible interactions between proteins and polysaccharides. Among these, 6 (CR_21, 22, 45, 47, 57 and 58) are found in both biominerals of the red coral and 4 proteins (CR_65, 70, 85 and 97) are only identified in the axial skeleton OM, suggesting that the two matrices may also differ on the basis of specific sugar–protein interactions. Out of the 10 identified carbohydrate-binding proteins, 5 are only found in octocorallian species (CR_21, 22, 47, 85 and 97). Interestingly, CR_22 may have arisen through horizontal gene transfer (HGT) from a prokaryote to an octocoral ancestor (Fig. 5, Additional file 9). In cnidarians, HGTs have already been demonstrated for other protein families [77,78,79].

Based on sequence similarity (BLAST search with evalue < 1e-15), 24 C. rubrum OM proteins appear to be restricted to octocorallians. Among them, 18 are found in all 4 octocoral databases investigated, whereas 6 are only found in C. rubrum. Out of the 6 red coral-specific proteins, 3 proteins are LCD-rich. LCDs have been described in the secreted OM of biominerals in different metazoan taxa [43,44,45,46]. In previous proteomic studies on stony coral skeletons, these types of proteins were called Secreted Aspartic Acid-Rich Proteins (SAARPs [16]) and Coral Acidic Rich-Proteins (CARPs [18]). Such highly acidic macromolecules have a high affinity to positively charged ions such as calcium and may thus play a role in CaCO3 formation [52, 80, 81]. On the other hand, 9 of the 24 octocorallian-specific proteins are strongly alkaline, including scleritin (CR_1). Scleritin was first identified as a sclerite specific protein [13].

Stony coral skeletons are composed of aragonite, whereas red coral skeletons are composed of calcite [7, 15,16,17,18]. Our proteomic study demonstrates that the underlying skeletomes of stony and red corals also differ, which extend a recent proteomic study realized on the sclerites of 3 octocoral and 1 scleractinian species [82]. Additionally, our study reveals similarities in the axis composition of alcyonaceans in the form of a collagen-rich matrix. Thus, we propose that the red coral axial skeleton biomineralization is the evolutionary result of the calcification of the gorgonian axis. This hypothesis is also supported by other evidence of axis mineralization in alcyonaceans. For instance, in the bamboo coral (Keratoisis genus) the axis is made of alternating structures of gorgonin (nodes) and calcite (internodes) [83]. In the sea rod Plexaurella nutans, calcite loculi precipitate in the axis [84]. However, only red corals produce a fully calcified axis. According to estimated phylogenetic divergence between Octocorallia and Hexacorallia (Precambrian) [53, 85, 86], calcification would have been acquired independently in these two subclasses during anthozoan evolution.

Calcification toolkit in anthozoans

Among the 77 proteins shared between octocorallians, actiniarians and hexacorallians, only 11 are present in the OM of both stony and red corals biominerals (Fig. 6, Additional file 3 sheet#2). The 11 proteins are distributed in 9 protein families (1 mucin, 1 CAP-like protein, 2 galaxin-like proteins, 2 protocadherins, 1 LRP, 1 CA, 1 collagen-like protein, 1 pikachurin-like protein and 1 unknown protein).

Mucins are a family of high molecular weight, heavily glycosylated proteins (glycol-conjugates) produced by epithelial tissues in most animals. In addition to their role in vertebrate mucus formation [87], they also participate in the calcification of mollusk shells, echinoderm skeletons and vertebrate bones [88,89,90]. Their conservation into the process of anthozoan calcification pinpoints them as major players into the mechanism of calcification in metazoans. In the CAP-like protein (CR_3-4), two domains were identified, GAPR-1 and TSP1. The GAPR-1 domain found has been shown to form amyloid-like fibrils in the presence of acidic phospholipids [91]. TSP1 repeats can physically interact with a variety of ligands, including structural components of the ECM, other matricellular proteins, growth factors, and proteases [92, 93]. Although a structural role may be envisaged, the function of CR_3-4 remains to be clarified. LRP protein was also previously identified in the skeletome of scleractinian corals [21] and echinoderms [94]. LRP4 is expressed in chondrocytes and promotes cell-autonomous cartilage growth in association with agrins [95]. However, LRPs are cell-surface (co)receptors involved in intracellular signal transduction and virtually nothing is known on the signaling control of calcification in anthozoans. Likewise, protocadherins (CR_24 and CR_34) are also known as trans-membrane proteins involved in cell–cell adhesion. The conserved function of these cell surface macromolecules in the calcification process remains to be deciphered.

Two galaxin-like proteins (CR_14 and CR_27) are identified in both the axial skeleton and the sclerites of the red coral. Galaxin was first identified in the exoskeleton of the scleractinian coral Galaxea fascicularis and was described as a tandem repeat structure with a di-cysteine motif fixed at nine positions [96]. Since this discovery, galaxin homologs have been observed in the exoskeleton of other scleractinian species [16, 21, 97]. Galaxin is associated with the developmental onset of calcification after larval stage in A. millepora [48]. Galaxin-like proteins were first suggested to be cnidarian-specific proteins [48, 96] until their discovery in mollusks and annelids [98,99,100]. In the squid Euprymna scolopes, a galaxin-like protein called EsGal1, was identified as an antimicrobial protein involved in the selection and modulation of growth of its symbiont, the Gram-negative bacterium Vibrio fischeri [98]. The high conservation of this motif in the different taxa suggests that this di-cysteine motif might play an essential role.

In metazoans, carbonic anhydrases (CAs) belong to multigenic family and are widely known to be involved in biomineralization in diverse metazoans such as sponge spicules [101, 102], mollusk shells [103, 104], sea urchin skeleton [105] and bird eggshells [106, 107], as well as in stony corals [16, 18, 50, 51, 108, 109]. In C. rubrum, 6 CAs were previously identified (CruCA1 to 6) and based on tissue expression profiles, CruCA4 was proposed to be the CA involved in the biomineralization process [25]. We now identified CruCA4 (CR_29) in both the axial skeleton and the sclerites, which supports its possible role in the biomineralization of the red coral. With regards to phylogeny, CruCA4 likely represents a duplicated copy of CruCA6, which is not phylogenetically related to other coral CAs involved in calcification [42]. Previous studies emphasize that the involvement of CAs in calcification likely represents a functional convergence of evolution of the calcification process [42, 45]. Even with some similar proteins in the calcification toolkit of stony and red corals, the independent acquisition of calcification is the most credible scenario. The presence of similar/homologous domains/proteins in the coral skeletomes could be the result of independent gene recruitment in anthozoan taxa.

Conclusions

This work is the first proteomic study of the organic matrix of the precious red coral (C. rubrum) biominerals. We report a total of 102 identified proteins, scleritin being the most abundant protein in the OM of both the sclerites and the axial skeleton. We noticed the abundant presence of collagen-like proteins (9 in each biomineral on 10 identified in our study), and highlighted 4 agrin-like proteins (3 in both biominerals and 1 more in the sclerites), which are ECM glycoproteins. These reuslts support the view that the calcified axial skeleton of the red coral evolved from the horny skeleton of gorgonians. Comparative analysis of the protein composition between the two biominerals pointed to differences in the processing (assembly/remodeling) of the OM rather than its structural composition. Indeed, specific enzymes (e.g. endopeptidases) and proteins harboring distinct sugar-binding domains are particular to each biomineral. Finally, we identified 18 octocorallian-specific proteins (of which 1 is a putative HGT from prokaryotes) and 6 unique C. rubrum novelties. Compared with stony corals, only few proteins were conserved (e.g. mucin, carbonic anhydrase, galaxins and protocadherins). The presence of these proteins in both octocorallian and scleractinian biomineralizations potentially represents a convergent co-option of calcification toolkit in corals.

Methods

Biological material

Corallium rubrum colonies were collected by the IMBE of Marseille in the Mediterranean Sea near the Marseille coast (Plane Island, Gulf of Lion: 43° 11.190′ N, 5° 23.470′ E), at a depth of 10 m. The collection was authorized by the Maritime Prefect of Bouches-du-Rhône (France) and the sampling did not involve permanent injuries, i.e. colonies recovered within two weeks after cutting. Colonies were transferred to the Centre Scientifique de Monaco to be cultured in open-circuit seawater aquaria supplied with Mediterranean seawater at a temperature of 18 °C ± 2 °C and were fed five times per week with frozen red plankton, rotifers and living artemia.

Extraction procedure of organic matrices

Axial skeletons and sclerites of red corals were isolated from living tissue by incubating branches in 10% hypochlorite under slow stirring and then rinsed five times in milli-Q water. Samples were subsequently dried at 60 °C. The axial skeleton and the sclerites were separated and washed 3 times in 10% hypochlorite for 24 h under agitation. Sclerites present in the axis cannot be separated as they are cemented in the medulla of the axial skeleton. Axis skeleton samples were pooled and crushed in liquid nitrogen for 5 min to obtain a homogeneous powder. The powder of axis and the non-crushed sclerites were cleaned in 10% hypochlorite overnight at 4 °C under stirring, then rinsed 8 times with milli-Q water and dried according to a published protocol [27]. The powder from pooled axis (10 g) and the pooled sclerites (5 g) were decalcified using 4% acetic acid for 24 h at 4 °C under stirring. Decalcified samples were centrifuged at 3500 g for 10 min at 4 °C to separate the supernatant containing the acidic-soluble fraction (soluble organic matrix—SOM) from the pellet containing the acidic-insoluble fraction (insoluble organic matrix—IOM). To remove acidic solutions, each SOM sample was diafiltered in a 3 kDa ultrafiltration cell (Amicon, Millipore) by centrifugation at 3500 g for 20 min at 4 °C. After the first passage of the decalcification solution, milli-Q water was added and centrifuged again. In total eight cycles of centrifugation-resuspension were performed. The pellets containing the IOMs were resuspended in milli-Q water, centrifuged and the supernatants discarded.

Electrophoresis and trypsin treatment of organic matrices

A total of 100 µg of SOM, and non-quantified IOM, extracted from pooled axial skeletons and sclerites (called SOMax, IOMax and SOMsc, IOMsc), were denatured in Laemmli buffer for 5 min at 95 °C. Proteins were loaded on a 12% polyacrylamide SDS-PAGE electrophoretic gel for rapid run (15 min) in order to obtain one single band. Gels were colored using Coomassie blue staining and rinsed several times with milli-Q water. Gels were then packed in sterile bags with 1% acetic acid until the trypsin treatment. The single bands were excised from the gels and cut into pieces. Reduction was performed using 10 mM dithiothreitol (DTT) in 100 mM NH4HCO3, 30 min at 37 °C. Gels were then cooled at room temperature and the excess DTT was removed. Alkylation was done with 55 mM Iodoacetamide in 100 mM NH4HCO3 for 1 h at room temperature, with occasional vortexing. Gels were washed with 100 mM NH4HCO3, dehydrated in 100 μl of acetonitrile and re-swelled 100 mM NH4HCO3 for 5 min. Gels were dehydrated once again with 100 μl of acetonitrile and dried in vacuum centrifuge (Savant Speedvac). For digestion, gels were firstly incubated in digestion solution containing 12.5 ng/μl of trypsin (Sequencing Grade Modified Trypsin, Promega, ref. V5111) in 50 mM NH4HCO3 during 30 min at 4 °C. Then, the excess of trypsin was removed and 50 mM NH4HCO3 was added on gels and trypsin in gel digested proteins overnight at 37 °C. After digestion, samples were cooled at room temperature and centrifuged 10 min at 6000 rpm. The supernatant was saved in a new microtube. Gels were then successively treated with 20 mM NH4HCO3 and 5% formic acid in 50% aqueous acetonitrile, centrifuge 5 to 10 min at 6000 rpm and supernatants were saved. The three supernatants/samples were finally pooled and dried in vacuum centrifuge.

NanoLC-MS/MS analysis

NanoLC-MS/MS analysis was performed once on each sample (SOMax, IOMax, SOMsc and IOMsc) using an on-line system consisting of a nano-pump UltiMate™ 3000 UHPLC binary HPLC system (Dionex, ThermoFisher) coupled to a Q-Exactive HF mass spectrometer (ThermoFisher, Germany). Peptides were resuspended in 20 μl of sample buffer (3% acetonitrile, 0.1% formic acid) and 2 μl was injected into a pre-column 300 µm × 5 mm (Acclaim PepMap, 5 µm particle size). After loading, peptides were eluted to an Acclaim PepMap100 C18 capillary column (75 µm × 15 cm, 100 Å, 3 μm particle sizes). Peptides were eluted into the MS, at a flow rate of 300 nl/min using a 40 min gradient from 5 to 40% mobile phase B. Mobile phase A was 0.1% formic acid in H2O and mobile phase B was 80% acetonitrile and 0.1% formic acid. The mass spectrometer was operated in positive and data-dependent mode, with a single MS scan (m/z 350–1400 at 60,000 resolution (at m/z 200) in a profile mode) followed by MS/MS scans on the 10 most intense ions at 15,000 resolution. Ions selected for MS/MS scan were fragmented using higher energy collision dissociation (HCD) at normalized collision energy of 28% and using an isolation window of m/z 1.8.

The raw mass spectrometric data have been deposited to PRIDE via the ProteomeXchange Consortium (https://www.ebi.ac.uk/pride/archive/) with the project accession PXD020332.

Protein identification and in silico characterization

The RAW files from Q-Exactive HF were converted into Mascot generic format (mgf) files using Proteome Discoverer v1.4 (Thermo Scientific). These files were submitted to Mascot v2.3 (Matrix Sciences Ltd, United Kingdom) for database search against the red coral transcriptome database of the Centre Scientifique de Monaco (CrubmRNA.fasta; data available at http://data.centrescientifique.mc/CSMdata-crubrum_data.html). The mass tolerance was set to 20 ppm for precursors, and 0.5 Da for the MS/MS fragment ion. The fixed modifications were set to carbamidomethyl and variable modifications were set to oxidation at methionine. The Mascot result files were processed using Scaffold v4.8.3 (Proteome Software Inc. USA; Additional file 13) software using the Prophet algorithm for validation of peptide and protein identifications with a threshold of 95% and 99%, respectively, and a minimum of two peptides in one of the 4 fractions for the identification of one protein. The relative abundance of identified proteins in each organic matrix fraction was estimated by calculating the exponentially modified protein abundance index (emPAI) using Scaffold v4.8.3. In order to compare the abundance of proteins in each condition, the weighted spectra (WS) value was also calculated as the number of MS/MS spectra for each protein (Additional file 3 sheet#1).

Identified proteins using Mascot were cross-checked against the published C. rubrum transcriptome [32], which used a different assembly strategy. This allowed us to merge several sequences from the CrubmRNA.fasta database into longer cDNA sequences (Additional file 2). Cellular localization was determined from the prediction of signal peptide (SignalP; http://www.cbs.dtu.dk/), transmembrane domain (TMHMM; http://www.cbs.dtu.dk/) and GPI-anchor (assigned GPI if both in predGPI and in GPISOM; http://gpcr.biocomp.unibo.it/predgpi/pred.htm: http://gpi.unibe.ch/). Sequence annotation was carried out using Blast2GO pipeline (B2G: https://www.blast2go.com/) using BLASTP against refseq (Eval > 1E-15) for homology search. ProtParam (http://web.expasy.org/protparam/) was used to determine in silico molecular weight and pI of proteins, net charge pH = 8 was calculated with PROTEIN CALCULATOR v3.4 (http://protcalc.sourceforge.net/), Low Complexity Region were determined using segmasker (SEG: http://mendel.imp.ac.at/METHODS/seg.server.html), and InterProScan (B2G) and SMART (http://smart.embl-heidelberg.de) were used for domain predictions (Additional file 3 sheet #2).

Multiple sequence alignment and phylogenetic analysis

Multiple sequence alignments for collagen-like and galaxin-like proteins were performed using Mafft (https://www.ebi.ac.uk/Tools/msa/mafft).

For CR_22 phylogenetic analysis, we did independent BLAST search against different taxonomic subdivisions of the NCBI_nr protein database: Cnidaria/Octocorallia, Cnidaria/Hexacorallia, Cnidaria/Hydrozoa, Deuterostomia, Protostomia, Placozoa, Porifira), as well as the Prokaryota subdivision. The first hit of each search as retrieved and subjected to clustalΩ MSA. Best-fit models of amino-acid replacement were calculated using ProtTest v3.4 [110]. Phylogenetic tree was reconstructed using PhyML v3.0 [111]. The tree was build using VT model with a gamma shape of 3.48, and likelihood topology was estimated by SPR (subtree pruning and regrafting) moves.

For BLAST homology search within anthozoan, local BLAST (2.2.27 +) was used. Species-specific databases were merged and interrogated as detailed in main text.

Scleractinian OM homologs to C. rubrum OM proteins were identified based on the published domain compositions [16, 18, 21] in addition to BLAST homology searches described above. We next retrieved the full-length sequences corresponding to the published sequences for A. millepora, A. digitifera and S. pistillata (which some were partial). Then, each set of C. rubrum protein and their respective putative homologs were aligned using Mafft. Finally, the percentages of amino-acid identity and similarity for each MSA were calculated using “Ident and Sim” (https://www.bioinformatics.org/sms2/ident_sim).

Availability of data and materials

The transcriptome of Corallium rubrum (CrubmRNA.fasta) that is used for analyses in the present study is available in the public database of the Centre Scientifique de Monaco (http://data.centrescientifique.mc/CSMdata-crubrum_data.html). The dataset supporting the conclusions of this article is included within the article and its Additional files 2 and 13. The raw mass spectrometric data have been deposited to PRIDE via the ProteomeXchange Consortium (https://www.ebi.ac.uk/pride/archive/) with the project accession PXD020332.

Abbreviations

CaCO3 :

Calcium carbonate

DTT:

Dithiothreitol

ECM:

Extracellular matrix

emPAI:

Exponentially modified protein abundance index

GAG:

Glycosaminoglycan

HGT:

Horizontal gene transfer

IOM:

Insoluble organic matrix

LRP:

Low-density lipoprotein receptor-related protein

MMP:

Matrix metalloproteinase

MSA:

Multiple sequence alignment

NanoLC-MS/MS:

Nano-liquid chromatography-tandem mass spectrometry

SOM:

Soluble organic matrix

OM:

Organic matrix

vWA:

Von Willebrand factor type A

WS:

Weighted spectral count

References

  1. Williams RJP. An introduction to biominerals and the role of organic molecules in their formation. Phil Trans R Soc Lond B. 1984;304:411–24.

    Article  CAS  Google Scholar 

  2. Muir PR, Wallace CC, Pichon M, Bongaerts P. High species richness and lineage diversity of reef corals in the mesophotic zone. Proc R Soc B. 2018;285:20181987.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Watabe N, Kingsley RJ. Calcification in octocorals. In: Suga S, Watabe N, editors. Hard tissue mineralization and demineralization. Tokyo: Springer; 1992. p. 127–47.

    Chapter  Google Scholar 

  4. Kingsley RJ, Watabe N. Ultrastructure of the axial region in Leptogorgia virgulata (Cnidaria: Gorgonaceae). Trans Am Microsc Soc. 1982;101:325.

    Article  Google Scholar 

  5. Goldberg WM. Evidence of a sclerotized collagen from the skeleton of a gorgonian coral. Comp Biochem Physiol Part B Comp Biochem. 1974;49:525–6.

    Article  CAS  Google Scholar 

  6. Dauphin Y. Mineralizing matrices in the skeletal axes of two Corallium species (Alcyonacea). Comp Biochem Physiol Part A Mol Integr Physiol. 2006;145:54–64.

    Article  CAS  Google Scholar 

  7. Vielzeuf D, Garrabou J, Baronnet A, Grauby O, Marschal C. Nano to macroscale biomineral architecture of red coral (Corallium rubrum). Am Mineral. 2008;93:1799–815.

    Article  CAS  Google Scholar 

  8. Guzman C, Shinzato C, Lu TM, Conaco C. Transcriptome analysis of the reef-building octocoral, Heliopora coerulea. Sci Rep. 2018;8:1–11.

    Article  CAS  Google Scholar 

  9. Jeon Y, Park SG, Lee N, Weber JA, Kim HS, Hwang SJ, et al. The draft genome of an octocoral Dendronephthya gigantea. Genome Biol Evol. 2019;11:949–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Lacaze-Duthiers H. Histoire naturelle du corail. Paris: J. B. Bail; 1864.

    Google Scholar 

  11. Grillo MC, Goldberg WM, Allemand D. Skeleton and sclerite formation in the precious red coral Corallium rubrum. Mar Biol. 1993;117:119–28.

    Article  Google Scholar 

  12. Marschal C, Garrabou J, Harmelin JG, Pichon M. A new method for measuring growth and age in the precious red coral Corallium rubrum (L.). Coral Reefs. 2004;23:423–32.

    Article  Google Scholar 

  13. Debreuil J, Tambutté É, Zoccola D, Deleury E, Guigonis JM, Samson M, et al. Molecular cloning and characterization of first organic matrix protein from sclerites of red coral, Corallium rubrum. J Biol Chem. 2012;287:19367–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Perrin J, Vielzeuf D, Ricolleau A, Dallaporta H, Valton S, Floquet N. Block-by-block and layer-by-layer growth modes in coral skeletons. Am Mineral. 2015;100:681–95.

    Article  Google Scholar 

  15. Rahman MA, Oomori T. Structure, crystallization and mineral composition of sclerites in the alcynarian coral. J Cryst Growth. 2008;310:3528–34.

    Article  CAS  Google Scholar 

  16. Ramos-Silva P, Kaandorp J, Huisman L, Marie B, Zanella-Cléon 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Goffredo S, Vergni P, Reggi M, Caroselli E, Sparla F, Levy O, et al. The skeletal organic matrix from Mediterranean coral Balanophyllia europaea influences calcium carbonate precipitation. PLoS ONE. 2011;6:e22338.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. 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 USA. 2013;110:3788–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Allemand D, Cuif J, Watabe N, Oishi M, Kawaguchi T. The organic matrix of skeletal structures of the Mediterranean red coral, Corallium rubrum. Bulletin de l’Institut océnaographique de Monaco. 1994;14:129–39.

    Google Scholar 

  20. Rahman MA, Oomori T, Wörheide G. Calcite formation in soft coral sclerites is determined by a single reactive extracellular protein. J Biol Chem. 2011;286:31638–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Takeuchi T, Yamada L, Shinzato C, Sawada H, Satoh N. Stepwise evolution of coral biomineralization revealed with genome-wide proteomics and transcriptomics. PLoS ONE. 2016. https://doi.org/10.1371/journal.pone.0156424.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Bentov S, Weil S, Glazer L, Sagi A, Berman A. Stabilization of amorphous calcium carbonate by phosphate rich organic matrix proteins and by single phosphoamino acids. J Struct Biol. 2010;171:207–15.

    Article  CAS  PubMed  Google Scholar 

  23. Söllner C, Burghammer M, Busch-Nentwich E, Berger J, Schwarz H, Riekel C, et al. Control of crystal size and lattice formation by starmaker in otolith biomineratization. Science. 2003;302:282–6.

    Article  PubMed  CAS  Google Scholar 

  24. Takeuchi T, Sarashina I, Iijima M, Endo K. In vitro regulation of CaCO3 crystal polymorphism by the highly acidic molluscan shell protein Aspein. FEBS Lett. 2008;582:591–6.

    Article  CAS  PubMed  Google Scholar 

  25. Falini G, Fermani S, Goffredo S. Coral biomineralization: a focus on intra-skeletal organic matrix and calcification. Semin Cell Dev Biol. 2015;46:17–26.

    Article  PubMed  Google Scholar 

  26. Allemand D. The biology and skeletogenesis of the Mediterranean red coral. Precious Corals Octocoral Res. 1993;2:19–39.

    Google Scholar 

  27. Debreuil J, Tambutté S, Zoccola D, Segonds N, Techer N, Allemand D, et al. Comparative analysis of the soluble organic matrix of axial skeleton and sclerites of Corallium rubrum: Insights for biomineralization. Comp Biochem Physiol Part B Biochem Mol Biol. 2011;159:40–8.

    Article  CAS  Google Scholar 

  28. Debreuil J, Tambutté S, Zoccola D, Segonds N, Techer N, Marschal C, et al. Specific organic matrix characteristics in skeletons of Corallium species. Mar Biol. 2011;158:2765–74.

    Article  CAS  Google Scholar 

  29. Rahman MA, Karl K, Nonaka M, Fujimura H, Shinjo R, Oomori T, et al. Characterization of the proteinaceous skeletal organic matrix from the precious coral Corallium konojoi. Proteomics. 2014;14:2600–6.

    Article  CAS  PubMed  Google Scholar 

  30. Kielty CM, Grant ME. The collagen family: structure, assembly, and organization in the extracellular matrix. In: Connective tissue and its heritable disorders. Hoboken: Wiley; 2002. p. 159–221.

    Chapter  Google Scholar 

  31. Cotman SL, Halfter W, Cole GJ. Identification of extracellular matrix ligands for the heparan sulfate proteoglycan agrin. Exp Cell Res. 1999;249:54–64.

    Article  CAS  PubMed  Google Scholar 

  32. Pratlong M, Haguenauer A, Chabrol O, Klopp C, Pontarotti P, Aurelle D. The red coral (Corallium rubrum) transcriptome: a new resource for population genetics and local adaptation studies. Mol Ecol Resour. 2015;15:1205–15.

    Article  CAS  PubMed  Google Scholar 

  33. Chaabane S, López Correa M, Ziveri P, Trotter J, Kalle N, Douville E, McCulloch M, Taviani M, Linares C, Montagna P. Elemental systematics of the calcite skeleton of Corallium rubrum and implications for the Mg/Ca temperature proxy. Chem Geol. 2019;524:237–58.

    Article  CAS  Google Scholar 

  34. Higdon R, Kolker E. A predictive model for identifying proteins by a single peptide match. Bioinformatics. 2007;23:277–80.

    Article  CAS  PubMed  Google Scholar 

  35. Exposito J-Y, Cluzel C, Garrone R, Lethias C. Evolution of collagens. Anat Rec. 2002;268:302–16.

    Article  CAS  PubMed  Google Scholar 

  36. Rodriguez-Pascual F, Slatter DA. Collagen cross-linking: insights on the evolution of metazoan extracellular matrix. Sci Rep. 2016;6:37374.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Voolstra CR, Li Y, Liew YJ, Baumgarten S, Zoccola D, Flot JF, et al. Comparative analysis of the genomes of Stylophora pistillata and Acropora digitifera provides evidence for extensive differences between species of corals. Sci Rep. 2017;7:17583.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Fidler AL, Vanacore RM, Chetyrkin SV, Pedchenko VK, Bhave G, Yin VP, et al. A unique covalent bond in basement membrane is a primordial innovation for tissue evolution. Proc Natl Acad Sci USA. 2014;111:331–6.

    Article  CAS  PubMed  Google Scholar 

  39. Zong Y, Zhang B, Gu S, Lee K, Zhou J, Yao G, et al. Structural basis of agrin-LRP4-MuSK signaling. Genes Dev. 2012;26:247–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Winzen U, Cole GJ, Halfter W. Agrin is a chimeric proteoglycan with the attachment sites for heparan sulfate/chondroitin sulfate located in two multiple serine-glycine clusters. J Biol Chem. 2003;278:30106–14.

    Article  CAS  PubMed  Google Scholar 

  41. Boraston AB, Healey M, Klassen J, Ficko-Blean E, Van Bueren AL, Law V. A structural and functional analysis of α-glucan recognition by family 25 and 26 carbohydrate-binding modules reveals a conserved mode of starch recognition. J Biol Chem. 2006;281:587–98.

    Article  CAS  PubMed  Google Scholar 

  42. Le Goff C, Ganot P, Zoccola D, Caminiti-Segonds N, Allemand D, Tambutté S. Carbonic anhydrases in cnidarians: Novel perspectives from the octocorallian Corallium rubrum. PLoS ONE. 2016;11:1–22.

    Google Scholar 

  43. Marie B, Joubert C, Tayalé A, Zanella-Cléon 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 USA. 2012;109:20986–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Jackson DJ, McDougall C, Woodcroft B, Moase P, Rose RA, Kube M, et al. Parallel evolution of nacre building gene sets in molluscs. Mol Biol Evol. 2010;27:591–608.

    Article  CAS  PubMed  Google Scholar 

  45. Le Roy N, Jackson DJ, Marie B, Ramos-Silva P, Marin F. The evolution of metazoan α-carbonic anhydrases and their roles in calcium carbonate biomineralization. Front Zool. 2014;11:75.

    Article  CAS  Google Scholar 

  46. Kocot KM, Aguilera F, Mcdougall C, Jackson DJ, Degnan BM. Sea shell diversity and rapidly evolving secretomes: insights into the evolution of biomineralization. Front Zool. 2016;13:23.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Yano M, Nagai K, Morimoto K, Miyamoto H. Shematrin: a family of glycine-rich structural proteins in the shell of the pearl oyster Pinctada fucata. Comp Biochem Physiol B Biochem Mol Biol. 2006;144:254–62.

    Article  PubMed  CAS  Google Scholar 

  48. Reyes-Bermudez A, Lin Z, Hayward DC, Miller DJ, Ball EE. Differential expression of three galaxin-related genes during settlement and metamorphosis in the scleractinian coral Acropora millepora. BMC Evol Biol. 2009;9:178.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  49. May P, Woldt E, Matz RL, Boucher P. The LDL receptor-related protein (LRP) family: An old family of proteins with new physiological functions. Ann Med. 2007;39:219–28.

    Article  CAS  PubMed  Google Scholar 

  50. Moya A, Tambutté S, Bertucci A, Tambutté E, Lotto S, Vullo D, et al. Carbonic anhydrase in the scleractinian coral Stylophora pistillata: Characterization, localization, and role in biomineralization. J Biol Chem. 2008;283:25475–84.

    Article  CAS  PubMed  Google Scholar 

  51. Bertucci A, Tambutté S, Supuran CT, Zoccola D. A new coral carbonic anhydrase in Stylophora pistillata. Mar Biotechnol. 2011;13:992–1002.

    Article  CAS  Google Scholar 

  52. Puverel S, Tambutté E, Pereira-Mouriès L, Zoccola D, Allemand D, Tambutté S. Soluble organic matrix of two scleractinian corals: Partial and comparative analysis. Comp Biochem Physiol Part B Biochem Mol Biol. 2005;141:480–7.

    Article  CAS  Google Scholar 

  53. Lin MF, Moya A, Ying H, Chen CA, Cooke I, Ball EE, et al. Analyses of corallimorpharian transcriptomes provide new perspectives on the evolution of calcification in the Scleractinia (corals). Genome Biol Evol. 2017;9:150–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Theocharis AD, Skandalis SS, Gialeli C, Karamanos NK. Extracellular matrix structure. Adv Drug Deliv Rev. 2016;97:4–27.

    Article  CAS  PubMed  Google Scholar 

  55. Young MF. Bone matrix proteins: their function, regulation, and relationship to osteoporosis. Osteoporos int. 2003;14(Suppl):3.

    Google Scholar 

  56. Halper J, Kjaer M. Basic components of connective tissues and extracellular matrix: elastin, fibrillin, fibulins, fibrinogen, fibronectin, laminin, tenascins and thrombospondins. In: Halper J, editor. Advances in experimental medicine biology. Dordrecht: Springer; 2014. p. 31–47.

    Google Scholar 

  57. Marie B, Arivalagan J, Mathéron L, Bolbach G, Berland S, Marie A, et al. Deep conservation of bivalve nacre proteins highlighted by shell matrix proteomics of the deep conservation of bivalve nacre proteins highlighted by shell matrix proteomics of the Unionoida Elliptio complanata and Villosa lienosa. J R Soc Interface. 2017;14:20160846.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  59. Tokuhara CK, Santesso MR, de Oliveira GSN, da Ventura TMS, Doyama JT, Zambuzzi WF, de Oliveira RC. Updating the role of matrix metalloproteinases in mineralized tissue and related diseases. J Appl Oral Sci. 2019. https://doi.org/10.1590/1678-7757-2018-0596.

    Article  PubMed  Google Scholar 

  60. Tang SY, Herber RP, Ho SP, Alliston T. Matrix metalloproteinase-13 is required for osteocytic perilacunar remodeling and maintains bone fracture resistance. J Bone Miner Res. 2012;27:1936–50.

    Article  CAS  PubMed  Google Scholar 

  61. Hou P, Troen T, Ovejero MC, Kirkegaard T, Andersen TL, Byrjalsen I, et al. Matrix metalloproteinase-12 (MMP-12) in osteoclasts: new lesson on the involvement of MMPs in bone resorption. Bone. 2004;34:37–47.

    Article  CAS  PubMed  Google Scholar 

  62. Nagase H, Visse R, Murphy G. Structure and function of matrix metalloproteinases and TIMPs. Cardiovasc Res. 2006;69:562–73.

    Article  CAS  PubMed  Google Scholar 

  63. Pihlajaniemi T, Myllylä R, Kivirikko KI. Prolyl 4-hydroxylase and its role in collagen synthesis. J Hepatol. 1991;13(Suppl 3):S2-7.

    Article  CAS  PubMed  Google Scholar 

  64. Walmsley AR, Batten MR, Lad U, Bulleid NJ. Intracellular retention of procollagen within the endoplasmic reticulum is mediated by prolyl 4-hydroxylase. J Biol Chem. 1999;274:14884–92.

    Article  CAS  PubMed  Google Scholar 

  65. Singh P, Carraher C, Schwarzbauer JE. Assembly of fibronectin extracellular matrix. Annu Rev Cell Dev Biol. 2010;26:397–419.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  66. Zollinger AJ, Smith ML. Fibronectin, the extracellular glue. Matrix Biol. 2017;60–61:27–37.

    Article  PubMed  CAS  Google Scholar 

  67. Chen X, Gumbiner BM. Paraxial protocadherin mediates cell sorting and tissue morphogenesis by regulating C-cadherin adhesion activity. J Cell Biol. 2006;174:301–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Bathish B, Turner R, Paumann-Page M, Kettle AJ, Winterbourn CC. Characterisation of peroxidasin activity in isolated extracellular matrix and direct detection of hypobromous acid formation. Arch Biochem Biophys. 2018;646:120–7.

    Article  CAS  PubMed  Google Scholar 

  69. Bailey D, Basar MA, Nag S, Bondhu N, Teng S, Duttaroy A. The essential requirement of an animal heme peroxidase protein during the wing maturation process in Drosophila. BMC Dev Biol. 2017;17:1.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  70. Kingsley RJ, Tsuzaki M, Watabe N, Mechanic GL. Collagen in the spicule organic matrix of the gorgonian Leptogorgia virgulata. Biol Bull. 1990;179:207–13.

    Article  CAS  PubMed  Google Scholar 

  71. Drickamer K, Taylor ME. Biology of animal lectins. Annu Rev Cell Biol. 1993;9:237–64.

    Article  CAS  PubMed  Google Scholar 

  72. Drickamer K. C-type lectin-like domains. Curr Opin Struct Biol. 1999;9:585–90.

    Article  CAS  PubMed  Google Scholar 

  73. Kilpatrick DC. Animal lectins: a historical introduction and overview. Biochim Biophys Acta Gen Subj. 2002;1572:187–97.

    Article  CAS  Google Scholar 

  74. Weiss IM, Kaufmann S, Mann K, Fritz M. Purification and characterization of perlucin and perlustrin, two new proteins from the shell of the mollusc Haliotis laevigata. Biochem Biophys Res Commun. 2000;267:17–21.

    Article  CAS  PubMed  Google Scholar 

  75. Flores RL, Gonzales K, Seaver RW, Livingston TB. The skeletal proteome of the brittle star Ophiothrix spiculata identifies C-type lectins and other proteins conserved in echinoderm skeleton formation. AIMS Mol Sci. 2016;3:357–67.

    Article  CAS  Google Scholar 

  76. Hincke MT, Tsang CPW, Courtney M, Hill V, Narbaitz R. Purification and immunochemistry of a soluble matrix protein of the chicken eggshell (ovocleidin 17). Calcif Tissue Int. 1995;56:578–83.

    Article  CAS  PubMed  Google Scholar 

  77. Denker E, Bapteste E, Le Guyader H, Manuel M, Rabet N. Horizontal gene transfer and the evolution of cnidarian stinging cells. Curr Biol. 2008;18:R858–9.

    Article  CAS  PubMed  Google Scholar 

  78. Chapman JA, Kirkness EF, Simakov O, Hampson SE, Mitros T, Weinmaier T, et al. The dynamic genome of Hydra. Nature. 2010;464:592–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Leclère L, Horin C, Chevalier S, Lapébie P, Dru P, Peron S, et al. The genome of the jellyfish Clytia hemispherica and the evolution of the cnidarian life-cycle. Nat Ecol Evol. 2019;3:801–10.

    Article  PubMed  Google Scholar 

  80. Addadi L, Weiner S. Interactions between acidic proteins and crystals: Stereochemical requirements in biomineralization. Proc Natl Acad Sci USA. 1985;82:4110–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Marin F, Le Roy N, Marie B. The formation and mineralization of mollusk shell. Front Biosci Sch. 2012;4:1099–125.

    Article  Google Scholar 

  82. Conci N, Lehmann M, Vargas S, Wörheide G. Comparative proteomics of octocoral and scleractinian skeletomes and the evolution of coral calcification. Genome Biol Evol. 2020. https://doi.org/10.1093/gbe/evaa162.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Noé SU, Dullo W-C. Skeletal morphogenesis and growth mode of modern and fossil deep-water isidid gorgonians (Octocorallia) in the West Pacific (New Zeland and Sea of Okhotsk). Coral Reefs. 2006;25:303–20.

    Article  Google Scholar 

  84. Lewis JC, Barnowski F, Telesnicki GJ. Characteristics of carbonates of gorgonian axes (Coelenterata, Octocorallia). Biol Bull. 1992;183:278–96.

    Article  CAS  PubMed  Google Scholar 

  85. Miglietta MP, McNally L, Cunningham CW. Evolution of calcium carbonate skeletons in the Hydractiniidae. Integr Comp Biol. 2010;50:428–35.

    Article  CAS  PubMed  Google Scholar 

  86. Conci N, Wörheide G, Vargas S. New non-bilaterian transcriptomes provide novel insights into the evolution of coral skeletomes. Genome Biol Evol. 2019;11:3068–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Meldrum OW, Yakubov GE, Bonilla MR, Deshmukh O, McGuckin MA, Gidley MJ. Mucin gel assembly is controlled by a collective action of non-mucin proteins, disulfide bridges, Ca2+ -mediated links, and hydrogen bonding. Sci Rep. 2018;8:5802.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

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

    Article  CAS  Google Scholar 

  89. Midura RJ, Hascall VC. Bone sialoprotein—a mucin in disguise? Glycobiology. 1996;6:677–81.

    Article  CAS  PubMed  Google Scholar 

  90. Boskey AL. Biomineralization: an overview. Connect Tissue Res. 2003;44(Suppl 1):5–9.

    Article  CAS  PubMed  Google Scholar 

  91. Sheng J, Olrichs NK, Geerts WJ, Kaloyanova DV, Helms JB. Metal ions and redox balance regulate distinct amyloid-like aggregation pathways of GAPR-1. Sci Reports. 2019;9:15048.

    Google Scholar 

  92. Olrichs NK, Mahalka AK, Kaloyanova D, Kinnunen PK, Bernd HJ. Golgi-associated plant pathogenesis related protein 1 (GAPR-1) forms amyloid-like fibrils by interaction with acidic phospholipids and inhibits Aβ aggregation. Amyloid. 2014;21:88–96.

    Article  CAS  PubMed  Google Scholar 

  93. Resovi A, Pinessi D, Chiorino G, Taraboletti G. Current understanding of the thrombospondin-1 interactome. Matrix Biol. 2014;37:83–91.

    Article  CAS  PubMed  Google Scholar 

  94. Flores RL, Livingston BT. The skeletal proteome of the sea star Patiria miniata and evolution of biomineralization in echinoderms. BMC Evol Biol. 2017. https://doi.org/10.1186/s12862-017-0978-z.

    Article  PubMed  PubMed Central  Google Scholar 

  95. Eldridge S, Nalesso G, Ismail H, Vicente-Greco K, Kabouridis P, Ramachandran M, et al. Agrin mediates chondrocyte homeostasis and requires both LRP4 and α-dystroglycan to enhance cartilage formation in vitro and in vivo. Ann Rheum Dis. 2016;75:1228–35.

    Article  CAS  PubMed  Google Scholar 

  96. Fukuda I, Ooki S, Fujita T, Murayama E, Nagasawa H, Isa Y, et al. Molecular cloning of a cDNA encoding a soluble protein in the coral exoskeleton. Biochem Biophys Res Commun. 2003;304:11–7.

    Article  CAS  PubMed  Google Scholar 

  97. Karako-Lampert S, Zoccola D, Salmon-Divon M, Katzenellenbogen M, Tambutté S, Bertucci A, et al. Transcriptome analysis of the scleractinian coral Stylophora pistillata. PLoS ONE. 2014. https://doi.org/10.1371/journal.pone.0088615.

    Article  PubMed  PubMed Central  Google Scholar 

  98. Heath-Heckman EAC, Gillette AA, Augustin R, Gillette MX, Goldman WE, McFall-Ngai MJ. Shaping the microenvironment: evidence for the influence of a host galaxin on symbiont acquisition and maintenance in the squid-vibrio symbiosis. Environ Microbiol. 2014;16:3669–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  99. Allam B, Pales Espinosa E, Tanguy A, Jeffroy F, Le Bris C, Paillard C. Transcriptional changes in Manila clam (Ruditapes philippinarum) in response to brown ring disease. Fish Shellfish Immunol. 2014;41:2–11.

    Article  CAS  PubMed  Google Scholar 

  100. Sanchez S, Hourdez S, Lallier FH. Identification of proteins involved in the functioning of Riftia pachyptila symbiosis by subtractive suppression hybridization. BMC Genomics. 2007;8:1–16.

    Article  CAS  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  102. Voigt O, Adamski M, Sluzek K, Adamska M. Calcareous sponge genomes reveal complex evolution of α-carbonic anhydrases and two key biomineralization enzymes. BMC Evol Biol. 2014. https://doi.org/10.1186/s12862-014-0230-z.

    Article  PubMed  PubMed Central  Google Scholar 

  103. Miyamoto H, Miyashita T, Okushima M, Nakano S, Morita T, Matsushiro A. A carbonic anhydrase from the nacreous layer in oyster pearls. Proc Natl Acad Sci USA. 1996;93:9657–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  104. Marie B, Le Roy N, Zanella-Cléon I, Becchi M, Marin F. Molecular evolution of mollusc shell proteins: Insights from proteomic analysis of the edible mussel Mytilus. J Mol Evol. 2011;72:531–46.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  Google Scholar 

  106. Mann K, Maček B, Olsen JV. Proteomic analysis of the acid-soluble organic matrix of the chicken calcified eggshell layer. Proteomics. 2006;6:3801–10.

    Article  CAS  PubMed  Google Scholar 

  107. Marie P, Labas V, Brionne A, Harichaux G, Hennequet-Antier C, Rodríguez-Navarro AB, et al. Quantitative proteomics provides new insights into chicken eggshell matrix protein functions during the primary events of mineralisation and the active calcification phase. J Proteomics. 2015;126:140–54.

    Article  CAS  PubMed  Google Scholar 

  108. Goreau TF. The physiology of skeleton formation in corals. I. A method for measuring the rate of calcium deposition by corals under different conditions. Biol Bull. 1959;116:59–75.

    Article  CAS  Google Scholar 

  109. Tambutté S, Tambutté E, Zoccola D, Caminiti N, Lotto S, Moya A, et al. Characterization and role of carbonic anhydrase in the calcification process of the azooxanthellate coral Tubastrea aurea. Mar Biol. 2007;151:71–83.

    Article  CAS  Google Scholar 

  110. Darriba D, Taboada GL, Doallo R, Posada D. ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011;27:1164–5.

    Article  CAS  PubMed  Google Scholar 

  111. Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors gratefully thank Anne Haguenauer and Dr. Didier Aurelle from the IMBE/Marseille for the C. rubrum sampling. We would like to thank M. François Rougaignon, President of the Fondation Paul Hamel for his confiance and support. We gratefully thank the reviewers for their constructive analysis of our work that significantly helped to improve this manuscript.

Funding

This work was supported by Postdoctoral funding of the Centre Scientifique de Monaco (n°500/685229; Government of Principality of Monaco). Proteomic analysis was supported by the Centre Scientifique de Monaco research program, funded by the Government of Principality of Monaco, and the C. rubrum genome and transcriptome analyses were supported by the Fondation Paul HAMEL.

Author information

Authors and Affiliations

Authors

Contributions

The study was designed by NLR and ST. NLR performed the preparation of samples, the organic matrix extractions and electrophoresis for proteomic analysis. MA generated mass spectrometry data. PG elaborated the transcriptome dataset used for de novo analysis. NLR, PG and MA analyzed and interpreted the proteomic dataset. NLR and PG wrote the original draft of the manuscript and NLR, PG, DA, ST, and MA reviewed and edited it. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Nathalie Le Roy.

Ethics declarations

Ethics approval and consent to participate

No ethical approval was required for any of the experimental research described here.

Consent for publication

Not applicable.

Competing interests

The authors declare that there have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

List of the 34 identified proteins in Corallium rubrum biominerals that are distributed into 9 protein families.

Additional file 2:

Genomic and transcriptomic sequences corresponding to the 107 proteins initially identified in Corallium rubrum biominerals. The final sequences for CR3_4, CR_10_28_60, and CR_62_76_82 are also listed.

Additional file 3:

Table of details for the 102 identified proteins in C. rubrum biominerals. For each identified protein in each OM fraction, unique peptide count, spectrum count, percentage of spectra, protein probability, percent protein coverage, exponentially modified protein abundance index (emPAI) and weighted spectra (WS) values are reported in sheet#1. Amino acid length, pI, charge at pH = 8, molecular weight (MW in Dalton Da), low complexity region, InterProscan protein signatures, description, subcellular (SubCell) localization prediction, BLAST (NCBI_ref-seq) result with first hit species and taxonomy, conservation in other anthozoan (BLAST e-val > 1e-15) and in scleractinian OM are listed in sheet#2. SOMax: Soluble Organic Matrix of the axial skeleton (dark blue), IOMax: Insoluble Organic Matrix of the axial skeleton (blue), SOMsc: Soluble Organic Matrix of the sclerites (yellow), IOMsc: Insoluble Organic Matrix of the sclerites (red). Empty circles: 1 unique peptide count.

Additional file 4:

Relative estimation of the proportion of sclerites into the axial skeleton using 7 published axial skeleton cross-sections [7, 33]. a: annular part; m: medullar part. Medullar and annular areas were drawn by hand after the published pictures with Photoshop®. The resulting png files were analyzed using ImageJ. The relative calculated areas are shown in the table.

Additional file 5:

Venn diagram distribution of the 102 identified proteins using 1 unique peptide threshold, in the IOM and SOM of the axial skeleton and the sclerites of C. rubrum.

Additional file 6:

List of the 10 most abundant proteins in each C. rubrum sample. a. Table of the 10 most abundant proteins in each OM biomineral fraction according to their exponentially modified protein abundance index (emPAI) value. Proteins highlighted in grey are in the top 10 of each sample, proteins highlighted in red are in the top 10 of SOM and IOM of the axial skeleton and proteins highlighted in blue are in the top 10 of SOM and IOM of the sclerites. b. Relative abundance in weighted spectra of the ten most abundant proteins in the red coral biominerals. SOMax: Soluble Organic Matrix of the axial skeleton (dark blue bars), IOMax: Insoluble Organic Matrix of the axial skeleton (blue bars), SOMsc: Soluble Organic Matrix of the sclerites (yellow bars), IOMsc: Insoluble Organic Matrix of the sclerites (red bars).

Additional file 7:

Multiple sequence alignment of the 10 collagen-like proteins identified in the proteome of the C. rubrum biominerals. Multiple alignment was carried out using Multalin (http://multalin.toulouse.inra.fr/multalin/multalin.html) and edited using JalView. Note the G-X–Y triplet repeats (glycine in blue, proline in red) characteristic of collagens.

Additional file 8:

Wrongly assigned collagen sequences in C. rubrum (top) and S. pistillata (bottom). Both in our analysis and in Drake and coworkers [18], the C. rubrum CR_98 and the S. pistillata XP_022783415 (P14) and AGG36343 (P18) sequences were wrongly assigned as “collagen”, probably because they are containing a von Willebrand factor type A domain (vWA) very homologous to another vWA domain sequence from another organism, which is there associated to a collagen protein sequence. Hence, by transitivity, the coral sequences were automatically annotated as “collagen”, although not being a collagen. Likewise, the C. rubrum CR_6 was also automatically annotated as “collagen”, as it contains a ColF1 domain. However, since the rest of the sequence is fibrinogen-related and does not contain G-X-Y triplets, CR_6 is assigned as “fibrillar collagen”, but not considered in the set of the collagens. COLFI: fibrillar collagen C-terminal domain; FA58C: coagulation factor 5/8 C-terminal domain; FBG: fibrinogen-related domain; TSP1: thrombonspondin-1.

Additional file 9:

CR_22 as a gene horizontally transferred from unicellular organisms. a. Results of the BLAST search in the different taxa in NCBI database (Porifera, Placozoa, Cnidaria/Octocorallia, Cnidaria/Hexacorallia, Cnidaria/Hydrozoa, Protostomia, Deuterostomia and Bacteria). b. Multiple alignment of CR_22 with its 2 octocorallian homologs (HeCoe: Heliopora coerulea; DeGig: Dendronephthya gigantea) and 2 prokaryotes (MyXan: Myxococcus xanthus; AcXin: Actinoplanes xinjiangensis). The two carbohydrate-binding domains CBM25 (Pfam) of CR_22 are framed in red box. c. Genomic exon (green squares)/intron (dashed blue lines) structure of the CR_22 and DeGig_CBM25 gene coding sequences. Numbers correspond to the introns’ sizes in base pair. d. Phylogenetic tree (PhyML) of CR_22 homologs and proteins corresponding to first BLAST hit from a.

Additional file 10:

List of the six anthozoan species used for the comparison with Corallium rubrum biominerals proteome. Nucleotide databases are indicated for each species.

Additional file 11:

List of putative homologs between Corallium rubrum biomineral proteomes and skeletomes of three scleractinian species Acropora millepora, Acropora digitifera and Styllophora pistilata [16, 18, 21].

Additional file 12:

Galaxin-like proteins. The published protein sequences for galaxin homologs from Acropora millepora, Galaxea fascicularis, Stylophora pistillata and Euprymna scolopes [16, 48, 96,97,98], as well as the present Corallium rubrum CR_14 and CR_27 were submitted to the MEME discovery motif and the MAST motif scanning programs (http://meme-suite.org/). The motifs is a double di-cysteine motif similar to galaxin motif found by Reyes-Bermundez and colleagues [48], and is shown as logo on top. On the listed fasta sequences, this double di-cysteine motif, or galaxin motif, is highlighted in yellow in each sequence. The signal peptide sequence (http://www.cbs.dtu.dk/services/SignalP-5.0/) is in italic and underlined.

Additional file 13:

List of proteins and peptides identified in the SOMs and IOMs of the axial skeleton and the sclerites of Corallium rubrum. For each identified protein, identified peptides are listed. For each peptide, amino acid sequence, spectra and Mascot information are listed. SOMax: soluble organic matrix of the axial skeleton (sheet#1); IOMax: insoluble organic matrix of the axial skeleton (sheet#2); SOMsc: soluble organic matrix of the sclerites (sheet#3); IOMsc: insoluble organic matrix of the sclerites (sheet#4).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Le Roy, N., Ganot, P., Aranda, M. et al. The skeletome of the red coral Corallium rubrum indicates an independent evolution of biomineralization process in octocorals. BMC Ecol Evo 21, 1 (2021). https://doi.org/10.1186/s12862-020-01734-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-020-01734-0

Keywords