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

Proteomic analysis of quail calcified eggshell matrix: a comparison to chicken and turkey eggshell proteomes

Abstract

Background

Eggshell mineralization in commercially important species such as chicken, turkey or quail is of interest as a general model of calcium carbonate biomineralization. Knowledge of proteins and molecular mechanisms in eggshell assembly may also pave the way to manipulation of thickness of the calcified layer or other features. Comparison of eggshell matrix proteomes of different species may contribute to a better understanding of the mineralization process. The recent publication of the quail genome sequence now enables the proteomic analysis of the quail shell matrix and this comparison with those of chicken and turkey.

Results

The quail eggshell proteome comprised 622 identified proteins, 311 of which were shared with chicken and turkey eggshell proteomes. Forty-eight major proteins (iBAQ-derived abundance higher than 0.1 % of total identified proteome) together covered 94 % of total proteome mass. Fifteen of these are also among the most abundant proteins in chicken and turkey eggshell matrix. Only three proteins with a percentage higher than 1.0 % of the total had not previously been identified as eggshell matrix proteins. These were an uncharacterized member of the latexin family, an uncharacterized protease inhibitor containing a Kunitz domain, and gastric intrinsic factor. The most abundant proteins were ovocleidin-116, ovalbumin and ovocalyxin-36 representing approximately 31, 13 and 8 % of the total identified proteome, respectively. The major phosphoproteins were ovocleidin-116 and osteopontin. While osteopontin phosphorylation sites were predominantly conserved between chicken and quail sequences, conservation was less in ovocleidin-116.

Conclusions

Ovocleidin-116 and ovocalyxin-36 are among the most abundant eggshell matrix proteins in all three species of the family Phasianidae analyzed so far, indicating that their presently unknown function is essential for eggshell mineralization. Evidence for other chicken eggshell-specific proteins in quail was inconclusive. Therefore measurement of additional eggshell proteomes, especially from species of different families and preferentially from outside the order Galliformes, will be necessary.

Background

The biological function of the avian eggshell is to provide protection against mechanical impact, to form a first line of defense against microbial infection, to regulate water and gas exchange and to provide calcium for the developing embryo [1]. The shell is formed in the eggshell gland (uterus) and consists to approximately 95 % of calcium carbonate in the form of calcite. The calcified layer is pervaded by an organic matrix consisting of proteins and proteoglycans, which is thought to provide a 3-dimensional network guiding and controlling the mineralization process [2]. The most studied avian eggshell is that of the chicken, probably because of commercial importance and easy accessibility. Among the first identified chicken eggshell matrix proteins were major egg white proteins, such as ovalbumin [3], lysozyme [4], and ovotransferrin [5]. These proteins were shown by immunohistochemical methods to be part of the calcified matrix rather than surface contaminants. Egg white proteins are produced and secreted predominantly in the magnum section of the oviduct [2]. Soluble remnants of egg the white assembly may migrate with the egg into the eggshell gland, where they are eventually incorporated into the mineralizing matrix. In addition, messages for lysozyme and ovotransferrin are also present at much lower concentration in white isthmus [4, 5]. These messages can even be detected in red isthmus and uterus after extensive amplification [4, 5]. At present it is unknown what percentage of these egg white proteins found in the eggshell matrix may be contributed by these tissues. Other proteins are produced by eggshell gland epithelial cells but apparently not in other sections of the oviduct and not in other selected tissues analyzed. The first of these so-called eggshell-specific proteins was ovocleidin-17 (OC17) [6], subsequently shown to belong to the C-type lectin-like family of proteins [7]. This was followed by cloning and characterization of ovocleidin-116 (OC116) [810], which was first identified by its N-terminal sequence in a dermatan sulfate proteoglycan preparation from chicken eggshell [11]. It was subsequently also detected in chicken bone [12, 13] and is thus not strictly eggshell-specific. Other members of this group with possible eggshell-specific distribution are ovocalyxin-32 (OCX32), a member of the latexin family of carboxypeptidase inhibitors [14], and ovocalyxin-36 (OCX36) [15], a member of the BPI/LBP/PLUNC family of anti-microbial proteins [16]. Other ovocalyxins occasionally mentioned in publications but poorly characterized are ovocalyxins-25 and −21. OCX-25 contains protease inhibitor domains and OCX-21 is apparently identical to gastrokine-2, a constituent of the gastric secretome [17]. Two other eggshell matrix proteins of widespread distribution, osteopontin and glypican-4, are induced in eggshell gland epithelia by the mechanical strain exerted upon entry of the egg into the gland [18, 19]. The exact role of these proteins remains unknown at present, but in vitro studies have shown that uterine fluid as well as isolated eggshell components may influence calcite crystallization in vitro and thus may also control eggshell mineralization in vivo [2024]. Other matrix proteins may participate in the egg’s anti-microbial defense during and after egg production [25, 26].

The almost complete sequencing of the chicken genome [27] provided the possibility to identify more eggshell matrix proteins using high-throughput mass-spectrometry-based proteomics. The first compartment of the egg to be analyzed by such methods was the acid-soluble fraction of the chicken eggshell calcified layer [28, 29], leading to identification of 528 proteins in a wide abundance range. Measurement of the acid-insoluble matrix [30, 31], the cuticle [32, 33], and the soluble fraction of the eggshell membranes in conjunction with the innermost eggshell calcified layer (mammillary cones) [34] yielded several additional protein identifications. Soluble eggshell membrane proteins were also analyzed at different stages of chicken embryo development [35]. In addition, the proteome of the uterine fluid bathing the egg during shell mineralization was compared to the proteome of the calcified shell [36]. In that study a total of 577 proteins were identified in uterine fluid and 466 in eggshell, with an overlap of 244. The proteomes of uterine fluid at different stages of eggshell mineralization has also been compared [37]. In sum, these studies identified a total of 675 eggshell proteins [37] and yielded important information on the distribution of shell proteins in different compartments and on the temporal sequence of their appearance in the uterine fluid. Proteomic studies were complemented by transcriptomic studies aiming at identifying genes expressed in uterus with possible importance for eggshell production [3842].

Compared to chicken, the eggshell proteomes of other species were less well explored. This was due to the lack of comprehensive sequence databases, still a prerequisite for high-throughput proteomics. The publication of the almost complete turkey genome sequence [43] provided the possibility to compare the eggshell proteomes of chicken and turkey [44]. The turkey eggshell matrix yielded 697 proteins. The overlap with the chicken eggshell proteome was 52 %. However, if only turkey proteins with an abundance of >0.01 % of the total were compared, the overlap increased to 95 %. This indicated that most of the major proteins were conserved between species, but that there were also a few potentially important differences.

In the present report, we use the recently published genome-derived sequence database of Japanese quail [45] to compare another eggshell proteome of the same avian family, the Phasianidae, to the preceding ones. Previously Western blotting analysis of eggshell matrices of several avian species using antibodies against known hen eggshell matrix proteins identified ovotransferrin, osteopontin and ovalbumin as components of the quail eggshell [46]. In addition, ovomucoid and lysozyme were identified by N-terminal sequence analysis of electrophoretically separated matrix proteins [47]. More recently, an unknown 32 kDa protein of the cuticle was characterized by a short N-terminal sequence [48], and the calcified matrix was shown to contain a protein with sequence similarity to ovocleidin-116 [49]. Ovocleidin-116, ovocleidin-17, ovocalyxin-32, clusterin, cystatin, lysozyme, osteopontin, ovalbumin, ovoinhibitor, ovomucoid and ovotransferrin were also identified in quail eggshell matrix in a previous proteomic study using the chicken sequences for protein identification [50]. Here we show that the number of quail eggshell matrix proteins is similar to that of chicken and turkey and explore similarities and differences among the major matrix proteins.

Materials and methods

Matrix and peptide preparation

Fresh quail eggs were bought at a local market. The broken shells were emptied, cleaned under a jet of water, and washed with 5 % EDTA at 6 °C for 30 min to facilitate mechanical removal of the cuticle and the membranes. The cuticles were then removed by brushing under a jet of de-ionized water, and pieces of calcified shell were stripped off the wet membranes. The dried pieces of calcified eggshell were demineralized in 50 % acetic acid (20 ml/g of shell) at 4–8 °C for 15 h with constant stirring. The turbid mixture was dialyzed (Spectra Por 6, cut off 2000; Spectrum Europe/Carl Roth GmbH, Karlsruhe, Germany) against 2 × 10 vol. of 10 % acetic acid and 2 × 10 vol. of 5 % acetic acid, and freeze dried.

SDS-PAGE was done using pre-cast 4–12 % Novex Bis-Tris gels in the MOPS buffer system using reagents and protocols supplied by the manufacturer (Invitrogen, Carlsbad, CA). The kit sample buffer was modified by adding β-mercaptoethanol to a final concentration of 2 %, and the sample was suspended in 30 μl sample buffer/100 μg of organic matrix and heated to 70 °C for 10 min. Samples were centrifuged to sediment insoluble material and gels were loaded with the dissolved proteins of 100 μg of matrix per lane and stained with colloidal Coomassie (Invitrogen) after electrophoresis. Gels were cut into 20 slices for in-gel reduction, carbamidomethylation and digestion with trypsin [51]. Peptides were cleaned with C18 Stage-Tips [52] before mass spectrometric analysis.

LC-MS and data analysis

Peptide mixtures were analyzed by in-line nanoflow liquid chromatography using the EASY-nLC system (Proxeon Biosystems, Odense, Denmark; now part of Thermo Fisher Scientific) with 15 cm capillary columns of an internal diameter of 75 μm filled with 3 μm Reprosil-Pur C18-AQ resin (Dr. Maisch GmbH, Ammerbuch-Entringen, Germany). The gradient consisted of 5–30 % acetonitrile in 0.5 % acetic acid at a flow rate of 250 nl/min for 85 min, 30–60 % acetonitrile in 0.5 % acetic acid at a flow rate of 250 nl/min for 5 min and 60-80 % acetonitrile in 0.5 % acetic acid at a flow rate of 250 nl/min for 7 min. The eluate was electrosprayed into an LTQ Orbitrap Velos (Thermo Fisher Scientific, Bremen, Germany) through a Proxeon nanoelectrospray ion source. The Orbitrap Velos was operated in an HCD top 10 mode essentially as described [53] at a resolution of 60,000 for full scans and of 7500 for fragment measurement (both specified at m/z 400). The dynamic exclusion time was 120 s.

Raw files were processed using the Andromeda search engine-based version 1.5.1.6 of MaxQuant (http://www.maxquant.org/) [5456] with enabled second peptide identification, iBAQ quantification and match between runs (match time window 0.5 min; alignment time window 20 min) options. We used the predicted gene database of Coturnix japonica [45] (http://www.nodai-genome.org/japanese_quail.html? lang = en; 30810 entries; downloaded November 2014) combined with a Coturnix subset of the UniProtKB database (Release 2014_9, 813 entries). In some trial searches, the sequences of chicken ovocleidins and ovocalyxins were added. The corresponding reversed databases and the sequences of common contaminants possibly introduced during sample preparation and handling were appended to the database. Carbamidomethylation was set as a fixed modification. Variable modifications were oxidation (M), N-acetyl (protein), pyro-Glu/Gln (N-term), and phosphorylation (Phospho (STY)). In some trial searches, hydroxyproline was added as a variable modification. Initial peptide mass tolerance and allowed MS/MS mass deviation were 20 ppm. Two missed cleavages were allowed and the minimal length required for a peptide identification was seven amino acids. The protein false discovery rate (FDR) and PSM (peptide spectral match) FDR were set to 0.01. The minimal peptide score for unmodified and modified peptides was set to 60. Identifications with less than three sequence-unique peptides were in each case validated using the MaxQuant Expert System software [57] considering the assignment of major peaks, occurrence of uninterrupted y- or b-ion series of at least four consecutive amino acids, preferred cleavages N-terminal to proline bonds, the possible presence of a2/b2 ion pairs and immonium ions, and mass accuracy. We only accepted protein identifications with at least two sequence-unique peptides occurring in at least two replicates with a total of three peptides. Exceptions were protein identifications sharing peptides with very similar proteins and fragments of proteins obviously belonging together. The iBAQ (intensity-based absolute quantification) [58] option of MaxQuant was used to calculate, based on the sum of peak intensities, the approximate share of each protein in the total identified proteome mass, and thus discerning major from minor proteins.

Sequence homology searches were performed with FASTA (http://www.ebi.ac.uk/Tools/sss/fasta/) [59] against current releases of Uniprot Knowledgebase (UniProtKB). We also searched the quail sequence databases for some proteins of particular interest using the Local Blast function [60] of BioEdit Sequence Alignment Editor version 7.2.5 from http://www.mbio.ncsu.edu/bioedit/bioedit.html and with the TBLASTN function of the BLAST site provided by the NODAI Genome Research Center (http://www.nodai-genome.org/blast/Japanesequail_scaf_1000/blast/). We also used Kalign (http://www.ebi.ac.uk/Tools/msa/kalign/) [59] and Clustal Omega for sequence alignments (http://www.ebi.ac.uk/Tools/msa/clustalo/) [61], InterPro (http://www.ebi.ac.uk/interpro/) [62] for domain predictions, and the Venn diagram plotter (http://omics.pnl.gov/software/VennDiagramPlotter.php) for preparing Venn diagrams. Kinase motifs apart from FAM20C consensus phosphorylation motifs were predicted using NetPhosK (http://www.cbs.dtu.dk/services/NetPhosK/) [63].

Results and discussion

In this study, we used the total PAGE sample buffer-soluble matrix of quail eggshell without separation of acid-soluble and acid-insoluble fractions because only approximately 10 % of the quail eggshell matrix is acid-soluble. Comparison of this total matrix fraction to the quail acid-soluble fraction alone showed that the band pattern remained the same but that some higher molecular weight bands were enriched in the total matrix (Fig. 1, lanes 2 and 3). In addition lane 1 of Fig. 1 shows the chicken eggshell acid-soluble matrix for comparison with the quail acid-soluble quail matrix (Fig. 1, lane 2). Both matrices showed similar complexity but low-molecular weight bands below 21 kDa were less prominent in quail than in chicken. In chicken this mobility range includes major matrix proteins, such as ovocleidin-17 and lysozyme.

Fig. 1
figure 1

PAGE separation of eggshell matrices. Lane 1, acid-soluble chicken eggshell matrix. Lane 2, acid-soluble quail eggshell matrix. Lane 3, total quail eggshell matrix. 100 μg of matrix were separated per lane. Molecular weight of markers is shown on the left in kDa. Gel sections used for in-gel digestion and analysis are indicated on the right

To ensure we obtained a representative average shell proteome, we pooled the calcified layers of eighteen washed eggshells (12.85 g) and performed demineralization and proteome analysis in triplicate. The yield of quail eggshell matrix was 37 mg/g of air-dried eggshell calcified layer. Proteins were separated by SDS-PAGE and the gels were sliced into 20 sections (Fig. 1) for in-gel reduction, carbamidomethylation and digestion with trypsin.

Processing the resulting 60 raw-files with MaxQuant yielded 1173 identified protein groups (Additional file 1: Table S1). Protein groups collect protein entries that cannot be distinguished by the obtained MS evidence (however, both proteins and protein groups are used interchangeably in the manuscript). The corresponding peptide data are shown in Additional file 2: Table S2. The amino acid sequences of identified accessions were compared to sequences contained in the UniProt Knowledgebase using FASTA. The results indicated that the NODAI quail sequence database obviously contained singles entries comprising several unrelated proteins as well as different entries likely containing fragments of the same protein. The first possibility is illustrated by entry 713 (Fig. 2). The sequence of its first 230 amino acids is nearly identical to the sequence of bone sialoprotein 2. This sequence region comprised three identified sequence-unique peptides. The sequence of amino acids 231–991 was very similar to chicken ovocleidin-116 and contained 45 identified sequence-unique peptides. The remainder of the sequence of 713 was very similar to osteopontin and the four peptides identified in this part of the entry (Fig. 2) were also contained in the Coturnix japonica osteopontin entry Q9I832_COTJA. An example of different entries containing fragments of the same protein is illustrated using accessions 14241 and 14979, which both matched to different regions of pentraxin (Q5UMH8) without overlap in matching regions (Fig. 3). Finally, entries 15278 and part of 3793 (3793b) are an example of entries already indicating their origin from one and the same protein by partially shared sets of identified peptides (Fig. 4). Both quail database entries were highly similar to cathepsin D (CATD_CHICK).

Fig. 2
figure 2

Analysis of accession 713 amino acid sequence. The sequence of quail accession no. 713 is aligned to chicken bone sialoprotein 2 (SIAL_CHICK), chicken ovocleidin-116 (OC116_CHICK) and chicken osteopontin (F1NSM8_CHICK). Identical amino acids are shaded yellow. Sequence regions covered by identified peptides are shown in bold red letters

Fig. 3
figure 3

Alignment of pentraxin fragments. Quail entries 14241 and 14979 matched to different sequence regions of pentraxin (Q5UMH8_CHICK) strongly indicating that both were fragments of the same protein. Identical amino acids are shaded yellow. Sequence regions covered by identified peptides are shown in bold red letters

Fig. 4
figure 4

Alignment of cathepsin D fragments. Quail database entry 15278 and part of 3793 (3793b) contained overlapping sequences and were both highly similar to chicken cathepsin D. Sequences confirmed by MS/MS sequencing are shown in red. Peptides identified only in one of the entries and not in the other one are underlined. Peptides matching only entry 3793b align to a gap in the alignment of this entry to chicken cathepsin D, and peptides matching only entry 15278 align to gaps in 3793b. Because of the very high sequence identity of both entries to chicken cathepsin D and the occurrence of several sequenced peptides shared by both entries we suggest that both entries are incomplete and likely belong to the same protein. Identical amino acid positions are shaded yellow

After elimination of identifications not matching the criteria detailed above (Materials and methods, LC-MS and data analysis) and tentatively grouping together protein fragments matching identical UniProt Knowledgebase accessions, we accepted 622 identified protein groups (Additional file 3: Table S3). However, this number remained tentative, because fragment assembly and grouping may contain errors and must be confirmed by determination of complete gene or mRNA sequences. Furthermore, proteins appear to be missing from the genome-derived database as indicated by identification of quail Uniprot protein entries not contained in the genome-derived database. The number of identified proteins is also lower than that of chicken (675; [37]) and turkey (697; [44]. Furthermore, identities between quail sequences and similar chicken or turkey sequences were frequently much lower than expected (Additional file 3: Table S3) and in almost all instances this was due to gaps in the respective sequence alignments caused by shorter quail sequences. The overlap of the quail shell proteome with turkey and chicken shell proteomes was approximately 63 and 60 %. Approximately 50 % of the proteins/protein groups were identified in all three species (Fig. 5).

Fig. 5
figure 5

Quantitative comparison of eggshell proteomes. The number of chicken proteins derived from different reports and fractions was taken from [37]. The overlap between turkey and chicken proteomes [44] was updated to include new data [34, 36]

As before [44, 64] we used intensity-based absolute quantification (iBAQ) as implemented in MaxQuant to discern major from minor proteins (Additional file 3: Table S3). To obtain more correct abundances for entries obviously containing sequences from different proteins, these were dissected into different database entries distinguished by small letters following the original accession number (Additional file 4: Figure S1). We also joined entries tentatively identified as fragments of the same protein to form a single accession according to the differentially shaded groups in Additional file 3: Table S3. This table contains abundances before and after database modification. Using threshold of ≥ 0.1 % of the total identified proteome to differentiate between minor and major proteins, we obtained 48 major proteins that together covered 94 % of the proteome (Table 1).

Table 1 Major proteins (≥0.1 %) of the quail eggshell calcified layer

Eggshell-specific proteins

A group of major eggshell proteins generally referred to as eggshell-specific proteins was once thought to occur specifically in the uterus because they were neither found in other sections of the oviduct nor in a few selected other tissues. These are the ovocleidins and ovocalyxins and they will be discussed in this section.

The first eggshell-specific protein to be detected in the chicken eggshell as a major matrix component was ovocleidin-17 (OC17) [6]. Sequence information of this C-type lectin-like protein was not contained in the published genome sequence of chicken but was known from direct sequence analysis of the isolated protein [7]. A full length clone containing the OC17 sequence was isolated and characterized only very recently and its mRNA expression level was found to correlate negatively with eggshell strength [65]. Furthermore, the message was found to be expressed in uterus and, at a lower level, in the preceding oviduct section, the isthmus, but not in 13 other tissues tested. However, its possible function in eggshell mineralization remains ill-defined at present and may include antimicrobial activity [25] and direct interaction with calcite surfaces or carbonate ions [22, 23, 66]. We could not identify turkey OC17 in our previous study [44], however, identification of a few peptides after addition of the chicken sequence to the searched database and previous detection of a protein in turkey eggshell matrix in Western blotting analysis using anti-chicken OC17 antiserum [46] indicated its presence in that eggshell matrix, but not in the genome sequence database. Homology searching of the quail sequence database with the chicken sequence did not yield a significant match and addition of the chicken sequence to the database searched by MaxQuant did not identify any matching peptide. This contradicts the previous identification of OC17 in quail eggshell using such a trans-species approach [50], but agrees with Western blotting results [46]. Very thick eggshells, such as the shells of ostrich [67], emu and rhea [68] all have two major proteins of this family in their shell matrix. Chicken and goose with their much thinner shell both only have one of these proteins [28, 69]. The apparent absence of an OC17 homolog in the even thinner shell of quail may thus indicate a connection between number or concentration of these proteins and eggshell thickness.

Ovoclein-116 was by far the most abundant protein of the quail and turkey eggshell matrix (Table 1), and also one of the most abundant chicken eggshell matrix components [28]. Most of the quail protein sequence was contained in the middle region of accession 713 (713b; Fig. 2). However, a second entry, accession 2298, was even more similar to chicken OC-116 than accession 713. However, we did not identify any peptides shared between these sequences. Closer examination showed that two overlapping, very abundant, peptides contained in 2298 (168 times identified; Additional file 3: Table S3) were located in the first 40 amino acids, preceded by a predicted secretion signal peptide. The N-terminus of the predicted mature protein was similar, but not identical, to a short N-terminal sequence published previously [49] (accession 2298 in Additional file 4: Figure S1). The remainder of the amino acid sequence of accession 2298 contained only two identified low-abundance peptides identified just four times altogether. In summary, the evidence indicated that accession 2298 contained the N-terminus of OC-116 missing in accession 713. We consequently combined this part of the sequence (aa1-40; 2298a) into one entry termed 713b to obtain a better estimate of OC-116 abundance. OC116 was mainly localized to the palisade layer [8], in agreement with its predominant secretion during the growth phase of eggshell mineralization [37]. The protein was also identified in chicken bone [12, 13], thus establishing a link between these two different biomineralization systems in this species. OC116 was suggested to have a mammalian counterpart, matrix extracellular phosphoprotein (MEPE), belonging to the secretory calcium-binding phosphoprotein (SCPP) group of proteins [70] or small integrin-binding ligand N-linked glycoproteins (SIBLING), a group of proteins including dentin matrix protein, osteopontin, dentin sialoprotein and other proteins important for skeletal and dental mineralization and remodeling [71, 72]. In fact, chicken OC116 is N-glycosylated [10] and phosphorylated [29]. While the phosphorylation status of the turkey protein was not determined, quail OC116 was also a major phosphoprotein of the eggshell matrix (see below).

Ovocalyxin-36 (OCX36) is related to lipopolysaccharide-binding proteins (LPS), bactericidal permeability-increasing proteins (BPIP) and Plunc family proteins [15]. Immunofluorescence staining indicated a distribution of the protein throughout the eggshell and the membrane with highest intensity at the interface between membranes and mammillary cones [15]. Its expression was greatest in the uterus during eggshell mineralization and its distribution was restricted to isthmus and uterus [15]. Its similarity to the LPS/BPI/PLUNC family of proteins and its ability to bind to bacterial lipopolysaccharide suggested a function in egg antimicrobial defense [15, 16, 73, 74]. A sequence with approximately 70 % identity to chicken OCX36 was among the most abundant proteins of the quail eggshell proteome (Additional file 3: Table S3). Thus OCX36 is a major component of the eggshell proteomes of all three species analyzed (Table 1), indicating an essential role in eggshell production.

Ovocalyxin-32 (OCX32) is a member of the latexin family of carboxypeptidase inhibitors [14], with unknown function in the eggshell matrix. It was most abundant in uterine fluid during the initial phase of mineralization [37]. We found OCX32 neither in turkey eggshell matrix [44] nor in quail matrix, despite the addition of the chicken sequence to databases for possible cross-species identification with peptides of identical sequence. This is in contrast to the previously reported presence of OCX32 in quail eggshell using the chicken sequence [50]. However, we identified a highly abundant quail database entry (accession 11366; Additional file 3: Table S3) with high similarity to latexin. Closer examination of the sequence and the distribution of identified peptides indicated that the highly abundant peptides all matched to the N-terminal 71 amino acids of this entry (11366a), which was not very similar to latexin but showed a weak similarity to chicken OCX32 (C7G541_CHICK; Additional file 3: Table S3; accession 11366 in Additional file 4: Figure S1) and contained a predicted proteinase inhibitor I47 (latexin) domain. The remainder of the sequence of this entry with high similarity to latexin (11366b) was represented by only one peptide identified only once. In conclusion, evidence for the presence of OCX32 in the quail eggshell matrix was inconclusive, but we unequivocally identified a highly abundant fragment of a protein belonging to the latexin family (Table 1) that may have a function similar to that of OCX32.

Two less characterized so-called eggshell-specific proteins occasionally appearing in the literature are ovocalyxin-25 (OCX25) and ovocalyxin-21 (OCX21). OCX25 was described as a protease inhibitor with WAP and Kunitz domains [37]. A Kunitz protease inhibitor-type domain was also predicted for the very abundant uncharacterized protein in quail accession 23208 (Table 1). The gene id given for chicken OCX25 [37], LOC771972, was linked to five chicken protein sequences of different length that all showed blocks of almost complete sequence identity separated by gaps in the alignments to the uncharacterized chicken protein F1NPR2. The N-terminal 50 amino acids of accession 23208 (length 120aa) also showed a weak similarity to the same chicken protein. This was also the accession retrieved by searching the quail sequence database with the sequence of LOC771972 and accession 23208 was 45 % identical to aa178-311 of LOC771972 (isoform X3, XP_004947248.1; Fig. 6). However, we are not sure whether this is sufficient to establish the protein encoded in entry 23208 as a homolog of OCX25.

Fig. 6
figure 6

Alignment of quail accession 23208 to LOC771972 (OCX25). The isoform if LOC771972 shown is X3 (XP_004947248.1). Identical amino acids are shaded yellow. Sequence regions covered by identified peptides are shown in bold red

OCX21 was among the most abundant proteins of the chicken eggshell matrix (IPI00574331.1) [28]. This protein is identical to gastrokine-2 (E1C2G7_CHICK) [32, 37], a protein of the gastric mucosal secretome and is therefore not truly eggshell-specific. We did not identify this protein in quail eggshell, although the quail sequence database contained two accessions (29506, 30653) with more than 90 % sequence identity to chicken and turkey OCX21/gastrokine-2 sequences. Therefore, this major chicken eggshell protein is not part of the quail eggshell proteome. Remarkably, however, all three species contain among their major shell matrix proteins a protein of the gastric secretome (see below).

Other major proteins

Egg white proteins, such as ovalbumin, lysozyme C and ovotransferrin were among the first eggshell proteins identified [35] and were invariably found among the major proteins of eggshell matrices [28, 44] (Table 1). Among its major components quail eggshell matrix contained ovalbumin, avidin, ovostatin-like protein, ovalbumin-related protein X, ovotransferrin, lysozyme C, cystatin and ovomucoid. All of these were also identified in chicken and turkey eggshell matrices at different abundances (Table 1), but ovalbumin was always among the most abundant eggshell matrix components, in agreement with its very high abundance in egg white. No specific function in the mineralization process at the molecular level has been reported for these proteins so far. Some of these proteins show antimicrobial activity mediated by different mechanisms, such as degradation of bacterial cell walls, inhibition of microbial proteases or sequestration of components essential for microbial growth [75]. Therefore, these proteins are thought to contribute to general egg defense against microbial contamination, which can also happen during shell growth in the uterus. However, because several of these proteins were also shown to influence calcium carbonate mineralization in vitro, a direct effect on eggshell mineralization cannot be excluded at present. Many of these egg white proteins are very abundant and therefore even weak binding of calcite could affect mineralization. Egg white proteins are predominantly produced in the magnum section of the oviduct and not in uterus [38, 42] and may migrate together with the unfinished egg from the site of their synthesis to the uterus. However, the messages for lysozyme and ovotransferrin [4, 5] were also detected at much lower levels in white isthmus, red isthmus and uterus tissues. At present it is not known whether these messages are translated into protein, and, if they are, what percentage of these abundant shell matrix proteins is derived from these alternative sources.

Major proteins that we identified only in quail eggshell matrix were an uncharacterized member of the latexin family (accession 11366a) and an uncharacterized proteinase inhibitor containing a predicted Kunitz domain (accession 23208). While the first may have a similar function as OCX32, the latter may be functionally related to OCX25. Another major new eggshell matrix protein was gastric intrinsic factor (GIF). GIF is part of the gastric mucosal secretome and binds cobalamin (vitamin B12). Interestingly each of the three proteomes considered in this report contained among its most abundant proteins a protein of gastric origin that was not identified in the respective other species’ eggshell proteome. In chicken this was gastrokine-2/OCX21 [28] and in turkey this was trefoil factor 2 (TFF2) [44]. The role of these gastric proteins in eggshell production and the origin of the eggshell matrix forms remain unknown at present.

EGF-like repeats and discoidin I-like domains 3 protein (EDIL3/Del1) and tsukushin were among the most abundant proteins with an iBAQ of more than 1 % (Table 1) in all three species analyzed. EDIL3 is a widespread extracellular calcium-binding protein that was most abundant in chicken uterine fluid during the initial phase of eggshell mineralization [37]. No connection of this protein to mineralization or eggshell production is known at present. This is also true for tsukushin, a member of the small leucine-rich proteoglycan family [76]. This protein was highly abundant in quail and turkey eggshell matrix (Table 1) but comparatively less abundant in chicken eggshell matrix. In chicken uterus fluid it reached the highest abundance during growth and terminal phase [37]. The expression of both genes was up-regulated when an unfinished egg was present in the uterus [42] indicating a specific, but presently unknown, function.

The list of major proteins (Table 1) contains two proteins that were previously shown to be induced and secreted in the uterus upon entry of the egg before the advent of large scale transcriptomic studies, and may be supposed to have some specific function in shell production. These are osteopontin and glypican-4 [18, 19]. Osteopontin is a member of the SIBLING family of mineralization-related secreted phosphoproteins [71, 72]. It was identified as a major protein of the shell matrix of all three phasianid species [28, 29, 44] (Table 1). The involvement of this multifunctional protein in mineralization processes is reviewed in [72, 77, 78]. Its activity is generally inhibitory and depends on phosphorylation. The localization of osteopontin predominantly at the surface of mammillary cones and eggshell pores and at the margins of calcite columns of the palisade layer supports an inhibitory function in eggshell mineralization by binding to selected crystal surfaces [24]. Glypican-4 is one of several related cell surface heparane sulfate proteoglycans with a GPI anchor and may be released from the surface by Notum protein (present in the quail eggshell matrix as a minor protein; Additional file 3: Table S3). Glypicans have been implicated in various regulatory processes at the cell surface [79], but no connection to mineralization events has been reported so far. Glypican-4 was identified in all three species analyzed, but was a major protein only in turkey and quail eggshell matrix (Table 1).

Also among the most abundant eggshell proteins in all three species were serum albumin, lactadherin (milk fat globule membrane protein 8), extracellular fatty acid-binding protein, the extracellular chaperone clusterin, extracellular serpin E2, the extracellular matrix protein fibronectin, and ubiquitin. Nothing is known about their specific role, if any, in eggshell production. Lactadherin and clusterin were most abundant in chicken uterus fluid during the initial phase of mineralization, while ubiquitin was most abundant in the terminal phase and albumin in the growth phase [37].

Phosphoproteins

We showed previously that higher energy collisional dissociation (HCD) fragmentation, also used in the present report, is well suited to determine peptide phosphorylation sites [80]. Applied to low-complexity proteomes such as those of biomineral matrices, this technique can yield useful information about major phosphoproteins and their phosphorylation sites without prior enrichment of modified peptides [64]. Because phosphorylation was reported previously to potentially affect biomineralization processes (reviewed in [8183]) we included it among the variable modifications used for MaxQuant search of our raw files.

In this way, we identified 21 phosphoproteins with a total of 56 different phosphorylation sites each with a localization probability in the sequence of >0.75 [84] (Table 2; only highest probability shown). This was less than in a previous study using phosphopeptide enrichment prior to analysis of the chicken eggshell phosphoproteome, which enables the identification of even traces of phosphorylation [29]. However, such sites may be less important for protein function, at least in proteins without catalytic or regulatory function. We compared the number of repeatedly identified phosphorylated versions of peptides to the number of non-phosphorylated forms of these peptides. This indicated that most sites were phosphorylated only partially (Table 2). The highest site occupancy, close to 100 %, was found for Ser346 of ovalbumin. This site and a second one (S69), with an estimated site occupancy of approximately 25 %, were already detected previously in several avian species including chicken [85]. In addition, we detected three previously unknown phosphorylation sites at Thr76, Ser77 and Thr266 of ovalbumin with very low site occupancy (Table 2). Ser80 in cystatin [86], corresponding to Ser103 in accession 14227 is a previously known phosphorylation site in a quail protein. However, we could not confirm the complete site occupancy reported previously (Table 2). As in chicken [29], the major phosphoproteins of the quail eggshell were OC116 with 10 different phosphosites (Fig. 7) and site occupancies up to approximately 80 %, and osteopontin with 14 different phosphosites and occupancy of up to 75 %.

Table 2 Phosphoproteins and phosphosites
Fig. 7
figure 7

Representative spectra of ovocleidin-116 phosphopeptides. HCD spectra of two selected ovocleidin-116 phosphopeptides (compare Table 2). Y-ions are shown in red, b-ions in blue, and water or ammonia losses in orange. *indicates loss of H3PO4 from a phosphorylated amino acid. Such losses can occur only C-terminal to the phosphorylated amino acid. In both spectra y7 is the most intense ion, due to the frequently observed preferential cleavage N-terminal to a proline. A tryptophane immonium ion (m/z 159.0922) in the lower spectrum is labeled Wimm

While nothing is known about the specific function of OC116 phosphorylation, osteopontin is one of the few matrix proteins better characterized in this regard. This multifunctional glycosylated phosphoprotein generally inhibits mineralization in a phosphorylation-dependent way. This was shown by comparison of phosphorylated osteopontin to enzymatically de-phosphorylated or un-phosphorylated recombinant osteopontin in in vitro crystallization assays using different minerals, such as calcium phosphate [87], calcium oxalate [88], or calcium carbonate [89]. Calcification of cultured human smooth muscle cells, used as a model of vascular calcification, was also inhibited by phosphorylated osteopontin, but not by an enzymatically de-phosphorylated form [90]. The extent of phosphorylation may also play a role. Moderately phosphorylated osteopontin, such as bone osteopontin, inhibited mineral formation in a hydroxyapatite crystallization assay, while the more heavily phosphorylated milk osteopontin rather promoted mineralization [91]. Eleven of fourteen quail osteopontin phosphorylation sites were previously identified in chicken eggshell osteopontin by mass spectrometric analysis after phosphopeptide enrichment [29] or by Edman degradation of peptides isolated from metabolically 32P-labeled cultured chicken osteoblast osteopontin [92] (Fig. 8). Differences may be species specific, tissue specific, or due to differences in isolation and detection methods.

Fig. 8
figure 8

Comparison of phosphorylation sites in quail and chicken osteopontin. Phosphorylated amino acid residues of eggshell osteopontin are highlighted by yellow shading. Phosphorylation sites identified in metabolically 32P-labeled chicken osteoblast osteopontin [92] are underlined

Comparison between quail and chicken ovocleidin-116 phosphorylation sites provided a different result (Fig. 9). Only four of the twelve quail OC116 phosphorylation sites were shared between the species, although eight of the amino acids phosphorylated in the quail protein were conserved in the chicken sequence. The higher number of phosphorylation sites in chicken OC116 may be due to phosphopeptide enrichment prior to analysis, enabling the identification of traces of phosphopeptides, possibly with no functional importance [93]. Overall, phosphorylation sites in OC116 were less well conserved than osteopontin phosphorylatin sites. It is possible that the overall extent of phosphorylation is more important in OC116 than conservation of exact positions.

Fig. 9
figure 9

Comparison of phosphorylation sites in quail and chicken ovocleidin-116. Phosphorylation sites are highlighted by yellow shading. Only phosphorylation sites with a site localization probability of >0.75 [84] are shown

A total of 12 of the phosphoproteins were previously identified as such in chicken [29], and three in mammals (Table 2). More than half of the phosphorylation sites (29 sites) agree with the consensus motif for phosphorylation by the secreted kinase FAM20C (S/T-X-[D,E,pS]), a kinase that appears to be identical with Golgi casein kinase and is known to modify many biomineralization-related extracellular proteins, such as the members of the SIBLING cluster [94, 95]. We also identified this kinase among the major quail eggshell matrix proteins (Table 1). Other kinase motifs detected in phosphopeptides were those recognized by casein kinase 2 (CH-II, four sites), protein kinase A (PKA, three sites), protein kinase C (PKC, two sites), pyruvate kinase (PK1, one site), ribosomal S6 kinase (RSK9, one site), and casein kinase 1 (CK1, one site). The phosphorylation site data, such as best score, mass error, and site probability as determined by MaxQuant are in Additional file 5: Table S4, which also contains data of peptides with lower site probability than 0.75 that were not included into Table 2.

Conclusions

The quail eggshell matrix proteome shares 50 % of its 622 identified proteins with the shell proteomes of chicken and turkey. Nevertheless we did not find the quail homologs of several major proteins of chicken and turkey eggshell matrix, including members of the so-called eggshell-specific group of proteins. Only ovocleidin-116 and ovocalyxin-36 were unequivocally identified as major eggshell proteins in all three species analyzed and these may therefore be essential for eggshell mineralization in general (Table 1). Osteopontin, like OC116 a member of the SIBLING cluster of biomineralization-related proteins, and a major phosphoprotein in all three eggshell matrices (Table 1), may also be essential. For almost all major proteins no specific function in eggshell mineralization is known. However, there are many suggestions for possible functions based on known enzymatic or binding activities of isolated proteins, results of in vitro crystallization experiments, immunohistochemical localization, transcriptomic studies relating message expression level to different shell mineralization stages or eggshell properties, and most recently, proteomic analysis of the shell itself at different mineralization stages [96]. It remains to be seen how this accumulated knowledge can help to elucidate specific functions of particular proteins at the molecular level. Proteins of general importance for eggshell matrix assembly and mineralization should be conserved in other species. Therefore proteomic analysis of other eggshell matrices, especially those of species not belonging to the order Galliformes, may provide further clues concerning distribution and importance of particular matrix proteins. However, our study also emphasizes the need for better, more comprehensive, and less redundant sequence databases to facilitate such comparative studies.

Abbreviations

Aa:

Amino acid

FDR:

False discovery rate

HCD:

Higher-energy collision-induced dissociation

iBAQ:

Intensity-based absolute quantification

OC:

Ovocleidin

OCX:

Ovocalyxin

PAGE:

Polyacrylamide gel electrophoresis

SIBLING:

Small integrin-binding ligand N-linked glycoproteins

References

  1. Solomon SE. The eggshell: strength, structure and function. Brit Poult Sci. 2010;51 Suppl 1:52–9.

    Article  Google Scholar 

  2. Hincke MT, Nys Y, Gautron J, Mann K, Rodriguez-Navarro AB, McKee MD. The eggshell: structure, composition and mineralization. Front Biosci. 2012;17:1266–80.

    Article  CAS  Google Scholar 

  3. Hincke MT. Ovalbumin is a component of the chicken eggshell matrix. Connect Tissue Res. 1995;31:227–33.

    Article  CAS  PubMed  Google Scholar 

  4. Hincke MT, Gautron J, Panhéleux M, Garcia-Ruiz J, McKee MD, Nys Y. Identification of lysozyme as a component of eggshell membranes and eggshell matrix. Matrix Biol. 2000;19:443–53.

    Article  CAS  PubMed  Google Scholar 

  5. Gautron J, Hincke MT, Panhéleux M, Garcia-Ruiz JM, Boldicke T, Nys Y. Ovotransferrin is a matrix protein of the hen eggshell membranes and basal calcified layer. Connect Tissue Res. 2001;42:255–67.

    Article  CAS  PubMed  Google Scholar 

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

    Article  CAS  Google Scholar 

  7. Mann K, Siedler F. The amino acid sequence of ovocleidin-17, a major protein of the avian eggshell calcified layer. Biochem Mol Biol Int. 1999;47:997–1007.

    CAS  PubMed  Google Scholar 

  8. Hincke MT, Gautron J, Tsang CPW, McKee MD, Nys Y. Molecular cloning and ultrastructural localization of the core protein of an eggshell matrix proteoglycan, ovocleidin-116. J Biol Chem. 1999;46:32915–23.

    Article  Google Scholar 

  9. Mann K, Hincke MT, Nys Y. Isolation of ovocleidin-116 from chicken eggshells, correction of its amino acid sequence and identification of disulfide bonds and glycosylated Asn. Matrix Biol. 2002;21:383–7.

    Article  CAS  PubMed  Google Scholar 

  10. Nimtz M, Conradt HS, Mann K. LacdiNAc (GalNAcβ1-4GlcNAc) is a major motif in N-glycan structures of the chicken eggshell protein ovocleidin-116. Biochim Biophys Acta. 2004;1675:71–80.

    Article  CAS  PubMed  Google Scholar 

  11. Carrino DA, Rodriguez JP, Caplan AI. Dermatan sulfate proteoglycans from the mineralized matrix of the avian eggshell. Connect Tissue Res. 1997;36:175–93.

    Article  CAS  PubMed  Google Scholar 

  12. Horvat-Gordon M, Yu F, Burns D, Leach RM. Ovocleidin (OC-116) is present in avian skeletal tissues. Poult Sci. 2008;87:1618–23.

    Article  CAS  PubMed  Google Scholar 

  13. Bardet C, Vincent C, Lajarille M-C, Jaffredo T, Sire J-Y. OC-116, the chicken ortholog of mammalian MEPE found in eggshell, is also expressed in bone cells. J Exp Zool. 2010;314B:653–62.

    Article  CAS  Google Scholar 

  14. Gautron J, Hincke MT, Mann K, Panhéleux M, Bain M, McKee MD, et al. Ovocalyxin-32, a novel chicken eggshell matrix protein. Isolation, amino acid sequencing, cloning and immunohistochemical localization. J Biol Chem. 2001;42:39243–52.

    Article  Google Scholar 

  15. Gautron J, Murayama E, Vignal A, Morisson M, McKee MD, Réhault S, et al. Cloning of ovocalyxin-36, a novel chicken eggshell protein related to lipopolysaccharide-binding proteins, bactericidal permeability-increasing proteins and Plunc family proteins. J Biol Chem. 2007;282:5273–86.

    Article  CAS  PubMed  Google Scholar 

  16. Gautron J, Réhault-Godbert S, Pascal G, Nys Y, Hincke MT. Ovocalyxin-36 and other LBP/BPI/Plunc-like proteins as molecular actors of the mechanisms of the avian egg natural defences. Biochem Soc Trans. 2011;39:971–6.

    Article  CAS  PubMed  Google Scholar 

  17. Nys Y, Gautron J, Garcia-Ruiz JM, Hincke MT. Avian eggshell mineralization: biochemical and functional characterization of matrix proteins. Comptes Rendus Palevol. 2004;3:549–62.

    Article  Google Scholar 

  18. Lavelin I, Meiri N, Pines M. New insight in eggshell formation. Poult Sci. 2000;79:1014–7.

    Article  CAS  PubMed  Google Scholar 

  19. Lavelin I, Meiri N, Einat M, Genina O, Pines M. Mechanical strain regulation of the chicken glypican-4 gene expression in the avian eggshell gland. Am J Physiol Integr Comp Physiol. 2002;283:R853–61.

    Article  Google Scholar 

  20. Gautron J, Bain M, Solomon SE, Nys Y. Soluble matrix of hen’s eggshell extracts changes in vitro the rate of calcium carbonate precipitation and crystal morphology. Brit Poult Sci. 1996;37:853–66.

    Article  CAS  Google Scholar 

  21. Dominguez-Vera JM, Gautron J, Garcia_Ruiz JM, Nys Y. The effect of avian uterine fluid on the growth behavior of calcite crystals. Poult Sci. 2000;79:901–7.

    Article  CAS  PubMed  Google Scholar 

  22. Marin-Garcia L, Frontana-Uribe BA, Reyes-Grajeda JP, Stojanoff V, Serrano-Posada HJ, Moreno A. Chemical recognition of carbonate anions by proteins involved in biomineralization processes and their influence on calcite crystal growth. Crystal Growth Design. 2008;8:1340–5.

    Article  CAS  Google Scholar 

  23. Hernández-Hernández A, Gómez-Morales J, Rodríguez-Navarro AB, Gautron J, Nys Y, García-Ruiz JM. Identification of some active proteins in the process of hen eggshell formation. Crystal Growth Design. 2008;8:4330–9.

    Article  Google Scholar 

  24. Chien Y-C, Hincke MT, Vali H, McKee MD. Ultrastructural matrix-mineral relationships in avian eggshell, and effects of osteopontin on calcite growth in vitro. J Struct Biol. 2008;163:84–99.

    Article  CAS  PubMed  Google Scholar 

  25. Wellman-Labadie O, Lakshminarayanan R, Hincke MT. Antimicrobial properties of avian eggshell-specific C-type lectin-like proteins. FEBS Lett. 2008;582:699–704.

    Article  CAS  PubMed  Google Scholar 

  26. Wellman-Labadie O, Picman J, Hincke MT. Antimicrobial activity of cuticle and outer eggshell protein extracts from three species of domestic birds. Brit Poultry Sci. 2008;49:133–43.

    Article  CAS  Google Scholar 

  27. International Chicken Genome Sequencing Consortium. Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Science. 2004;432:695–716.

    Google Scholar 

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

  29. Mann K, Olsen JV, Maček B, Gnad F, Mann M. Phosphoproteins of the chicken eggshell calcified layer. Proteomics. 2007;7:106–15.

    Article  CAS  PubMed  Google Scholar 

  30. Mikšik I, Eckhardt A, Sedláková P, Mikulíková K. Proteins of insoluble matrix of avian (Gallus gallus) eggshell. Connect Tissue Res. 2007;48:1–8.

    Article  PubMed  Google Scholar 

  31. Mikšik I, Sedláková P, Lacinová K, Pataridis S, Eckhardt A. Determination of insoluble avian eggshell matrix proteins. Anal Bioanal Chem. 2010;397:205–14.

    Article  PubMed  Google Scholar 

  32. Rose-Martel M, Du J, Hincke MT. Proteomic analysis provides new insight into the eggshell cuticle. J Proteom. 2012;75:2697–706.

    Article  CAS  Google Scholar 

  33. Mikšík I, Ergang P, Pácha J. Proteomic analysis of chicken eggshell cuticle membrane layer. Anal Bioanal Chem. 2014;406:7633–40.

    Article  PubMed  Google Scholar 

  34. Rose-Martel M, Smiley S, Hincke MT. Novel identification of matrix proteins involved in calcitic biomineralization. J Proteomics. 2015;116:81–96.

    Article  CAS  PubMed  Google Scholar 

  35. Kaweewong K, Garnjanagoonchorn W, Jirapakkul W, Roytrakul S. Solubilization and identification of hen eggshell membrane proteins during different times of chicken embryo development using the proteomic approach. Protein J. 2013;32:297–308.

    Article  CAS  PubMed  Google Scholar 

  36. Sun C, Xu G, Yang N. Differential label-free quantitative proteomic analysis of avian eggshell matrix and uterine fluid proteins associated with eggshell mechanical property. Proteomics. 2013;13:3523–36.

    Article  CAS  PubMed  Google Scholar 

  37. Marie P, Labas V, Brionne A, Harichaux G, Hennequet-Antier C, Nys Y, et al. Quantitative proteomics and bioinformatics analysis provide new insight into protein function during avian eggshell biomineralization. J Proteomics. 2015;113:178–93.

    Article  CAS  PubMed  Google Scholar 

  38. Dunn IC, Wilson PW, Lu Z, Bain MM, Crossan CL, Talbot RT, et al. New hypotheses on the function of the avian shell gland derived from microarray analysis comparing tissue from juvenile and sexually mature hens. Gen Comp Endocr. 2009;163:225–32.

    Article  CAS  PubMed  Google Scholar 

  39. Jonchère V, Réhault-Godbert S, Hennequet-Antier C, Canau C, Sibut V, Cogburn LA, et al. Gene expression profiling to identify eggshell proteins involved in physical defense of the chicken egg. BMC Genomics. 2010;11:55.

    Article  Google Scholar 

  40. Jonchère V, Brionne A, Gautron Y, Nys Y. Identification of uterine ion transporters for mineralization precursors of the avian eggshell. BMC Physiol. 2012;12:10.

    Article  PubMed Central  PubMed  Google Scholar 

  41. Fan Y-F, Hou Z-C, Yi G-Q, Xu G-Y, Yang N. The sodium channel gene family is specifically expressed in hen uterus and associated with eggshell quality traits. BMC Genet. 2013;14:90.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Brionne A, Nys Y, Hennequet-Antier C, Gautron J. Hen uterine gene expression profiling during eggshell formation reveals putative proteins involved in the supply of minnerals or in the shell mineralization process. BMC Genomics. 2014;15:220.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Dalloul RA, Long JA, Zimin AV, Aslam L, Beal K, Blomberg LE, et al. Multi-Platform next-generation sequencing of the domestic turkey (Meleagris gallopavo): Genome Assembly and analysis. PLoS Biol. 2010;8, e1000475.

    Article  PubMed Central  PubMed  Google Scholar 

  44. Mann K, Mann M. The proteome of the calcified layer organic matrix of turkey (Meleagris gallopavo) eggshell. Proteome Sci. 2013;11:40.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Kawahara-Miki R, Sano S, Nunome M, Shimmura T, Kuwayama T, Takahashi S, et al. Next-generation sequencing reveals genomic features of Japanese quail. Genomics. 2013;101:345–53.

    Article  CAS  PubMed  Google Scholar 

  46. Panhéleux M, Bain M, Fernandez MS, Morales I, Gautron J, Arias JL, et al. Organic matrix composition and ultrastructure of eggshell: a comparative study. Brit Poultry Sci. 1999;40:240–52.

    Article  Google Scholar 

  47. Lakshminarayanan R, Loh XJ, Gayathri S, Sindhu S, Banerjee Y, Kini RM, et al. Formation of transient amorphous calcium carbonate precursor in quail eggshell mineralization: an in vitro study. Biomacromolecules. 2006;7:3202–9.

    Article  CAS  PubMed  Google Scholar 

  48. Rahman MA, Moriyama A, Iwasawa A, Yoshizaki N. Cuticle formation in quail eggs. Zool Sci. 2009;26:496–9.

    Article  CAS  PubMed  Google Scholar 

  49. Iwasawa A, Rahman MA, Roy TK, Moriyama A, Yoshizaki N. Morphological and histochemical changes in the uterus epithelium during eggshell formation in quail. J Poult Sci. 2010;47:183–9.

    Article  CAS  Google Scholar 

  50. Portugal SJ, Cooper HJ, Zampronio CG, Wallace LL, Cassey P. Can museum egg specimens be used for proteomic analysis? Proteome Sci. 2010;8:40.

    Article  PubMed Central  PubMed  Google Scholar 

  51. Shevchenko A, Tomas H, Havlis J, Olsen JV, Mann M. In-gel digestion for mass spectrometric characterization of proteins and proteomes. Nat Protoc. 2006;1:2856–60.

    Article  CAS  PubMed  Google Scholar 

  52. Rappsilber J, Mann M, Ishihama Y. Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using StageTips. Nat Protoc. 2007;2:1896–906.

    Article  CAS  PubMed  Google Scholar 

  53. Olsen JV, Schwartz JC, Griep-Raming J, Nielsen ML, Damoc E, Denisov E, et al. A dual pressure linear ion trap-Orbitrap instrument with very high sequencing speed. Mol Cell Proteomics. 2009;8:2759–69.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized ppb-range mass accuracies and proteome-wide protein quantification. Nature Biotechnol. 2009;26:1367–72.

    Article  Google Scholar 

  55. Cox J, Neuhauser N, Michalski A, Scheltema RA, Olsen JV, Mann M. Andromeda – a peptide search engine integrated into the MaxQuant environment. J Proteome Res. 2011;10:1794–805.

    Article  CAS  PubMed  Google Scholar 

  56. Tyanova S, Temu T, Carlson A, Sinitcyn P, Mann M, Cox J. Visualization of LC-MS/MS proteomics data in MaxQuant. Proteomics. 2015;15:1453–6.

    Article  CAS  PubMed  Google Scholar 

  57. Neuhauser N, Michalski A, Cox J, Mann M. Expert system for computer-assisted annotation of MS/MS spectra. Mol Cell Proteom. 2012;11:1500–9.

    Article  Google Scholar 

  58. Schwanhäusser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, et al. Global quantification of mammalian gene expression control. Nature. 2011;473:337–42.

    Article  PubMed  Google Scholar 

  59. Goujon M, McWilliam H, Li W, Valentin F, Squizzato S, Paern J, et al. A new bioinformatics analysis tools framework at EMBL-EBI (2010). Nucleic Acids Res. 2010;38(Suppl):W695–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Altschul SF, Thomas TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1977;25:3389–402.

    Article  Google Scholar 

  61. Sievers F, Wilm A, Dineen DG, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539.

    Article  PubMed Central  PubMed  Google Scholar 

  62. Hunter S, Jones P, Mitchell A, Apweiler R, Attwood TK, Bateman A, et al. InterPro in 2011: new developments in the family and domain prediction database. Nucleic Acids Res. 2011;40:D306–12.

    Article  PubMed Central  PubMed  Google Scholar 

  63. Blom N, Sicheritz-Ponten T, Gupta R, Gammeltoft S, Brunak S. Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics. 2004;4:1633–49.

    Article  CAS  PubMed  Google Scholar 

  64. Mann M, Edsinger E. The Lottia gigantea shell matrix proteome: re-analysis including MaxQuant iBAQ quantitation and phosphoproteome analysis. Proteome Sci. 2014;12:28.

    Article  PubMed Central  PubMed  Google Scholar 

  65. Zhang Q, Liu L, Zhu F, Ning Z, Hincke MT, Yang N, et al. Integrating de novo transcriptome assembly and cloning to obtain chicken ovocleidin-17 full-length cDNA. Plos One. 2014;9, e93452.

    Article  PubMed Central  PubMed  Google Scholar 

  66. Freeman CL, Harding JH, Quigley D, Rodger PM. Simulations of ovocleidin17 binding to calcite surfaces and its implications for eggshell formation. J Phys Chem. 2011;115:8175–83.

    CAS  Google Scholar 

  67. Mann K, Siedler F. Ostrich (Struthio camelus) eggshell contains two different C-type lectin-like proteins. Isolation, amino acid sequence, and post-translational modifications. Biochim Biophys Acta. 2004;1696:41–50.

    Article  CAS  PubMed  Google Scholar 

  68. Mann K, Siedler F. Amino acid sequences and phosphophorylation sites of emu and rhea eggshell C-type lectin-like proteins. Comp Biochem Physiol B. 2006;143:160–70.

    Article  PubMed  Google Scholar 

  69. Lakshminarayanan R, Valiyaveetil S, Rao VS, Kini RM. Purification, characterization, and in vitro mineralization studies of a novel goose eggshell matrix protein, ansocalcin. J Biol Chem. 2003;278:2928–36.

    Article  CAS  PubMed  Google Scholar 

  70. Kawasaki K, Weiss KM. Evolutionary genetics of vertebrate tissue mineralization: the origin and evolution of the decretory calcium-binding phosphoprotein family. J Exp Zool (Mol Dev Evol). 2006;306B:295–316.

    Article  CAS  Google Scholar 

  71. Rowe PSN. The chicken or the egg: PHEX, FGF23 and SIBLINGs unscrambled. Cell Biochem Funct. 2012;30:355–75.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  72. Staines KA, MacRae VE, Farquharson C. The importance of the SIBLING family of proteins on skeletal mineralization and bone remodeling. J Endocrinol. 2012;214:241–55.

    Article  CAS  PubMed  Google Scholar 

  73. Cordeiro CMM, Esmaili H, Ansah G, Hincke MT. Ovocalyxin-36 is a pattern recognition protein in chicken eggshell membranes. PLoS One. 2013;8, e84112.

    Article  PubMed Central  PubMed  Google Scholar 

  74. Kovacs-Nolan J, Cordeiro C, Young D, Mine Y, Hincke MT. Ovocalyxin-36 is an effector protein modulating the production of proinflammatory mediators. Vet Immunol Immunop. 2014;160:1–11.

    Article  CAS  Google Scholar 

  75. Abdou AM, Sato K, Kim M. Functional proteins and peptides of hen’s egg origin. INTECH Open Access Publisher;2013. http://www.intechopen.com/books/bioactive-food-peptides-in-health-and-disease/functional-proteins-and-peptides-of-hen-s-egg-origin

  76. Dellet M, Hu W, Papadaki V, Ohnuma S. Small leucine rich proteoglycan family regulates multiple signaling pathways in neural development and maintenance. Dev Growth Differ. 2012;54:327–40.

    Article  Google Scholar 

  77. Giachelli CM, Steitz S. Osteopontin: a versatile regulator of inflammation and biomineralization. Matrix Biol. 2000;19:615–22.

    Article  CAS  PubMed  Google Scholar 

  78. Sodek J, Ganss B, McKee MD. Osteopontin. Crit Rev Oral Biol Med. 2000;11:279–303.

    Article  CAS  PubMed  Google Scholar 

  79. Filmus J, Capurro M, Rast J. Glypicans. Genome Biol. 2008;9, e224.

    Article  Google Scholar 

  80. Nagaraj N, D’Sousa RCJ, Cox J, Olsen JV, Mann M. Feasibility of large-scale phosphoproteomics with higher energy collisional dissociation fragmentation. J Proteome Res. 2010;9:6786–94.

    Article  CAS  PubMed  Google Scholar 

  81. Veis A, Sfeir C, Wu CB. Phosphorylation of the proteins of the extracellular matrix of mineralized tissues by casein kinase-like activity. Crit Rev Oral Biol Med. 1997;8:360–79.

    Article  CAS  PubMed  Google Scholar 

  82. George A, Veis A. Phosphorylated proteins and control over apatite nucleation, crystal growth, and inhibition. Chem Rev. 2008;108:4670–93.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  83. Alvares K. The role of acidic phosphoproteins in biomineralization. Connect Tissue Res. 2014;55:34–40.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  84. Olsen JV, Blogoev B, Gnad F, Maček B, Kumar C, Mortensen P, et al. Global, in vivo, and site-specific phosphorylation dynamics in signaling networks. Cell. 2006;127:635–48.

    Article  CAS  PubMed  Google Scholar 

  85. Henderson JY, Moir AJ, Fothergill LA, Fothergill JE. Sequences of sixteen phosphoserine peptides from ovalbumins of eight species. Eur J Biochem. 1981;114:439–50.

    Article  CAS  PubMed  Google Scholar 

  86. Gerhartz B, Engh RA, Mentele R, Eckerskorn C, Torquato R, Wittmann J, et al. Quail cystatin: isolation and characterization of a new member of the cystatin family and its hypothetical interaction with cathepsin B. FEBS Lett. 1997;412:551–8.

    Article  CAS  PubMed  Google Scholar 

  87. Hunter GK, Kyle CL, Goldberg HA. Modulation of crystal formation by bone phosphoproteins: structural specificity of the osteopontin-mediated inhibition of hydroxyapatite formation. Biochem J. 1994;300:723–8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  88. Hunter GK, Grohe B, Jeffrey S, O’Young J, Sørensen ES, Goldberg HA. Role of phosphate groups in inhibition of calcium oxalate crystal growth by osteopontin. Cells Tissues Organs. 2009;189:44–50.

    Article  CAS  PubMed  Google Scholar 

  89. Hincke MT, St.Maurice M. Phosphorylation-dependent modulation of calcium carbonate precipitation by chicken eggshell matrix proteins. In: Goldberg M, Boskey A, Robinson C, editors. Chemistry and biology of mineralized tissues. Rosemont: American Academy of Orthopaedic Surgeons; 2000. p. 13–7.

    Google Scholar 

  90. Jono S, Peinado C, Giachelli CM. Phosphorylation of osteopontin is required for inhibition of vascular smooth muscle cell calcification. J Biol Chem. 2000;275:20197–203.

    Article  CAS  PubMed  Google Scholar 

  91. Gericke A, Quin C, Spevak L, Fujimoto Y, Butler WT, Sørensen ES, et al. Importance of phosphorylation for osteopontin regulation of biomineralization. Calcif Tissue Int. 2005;77:45–54.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  92. Salih E, Ashkar S, Gerstenfeld LC, Glimcher MJ. Identification of the phosphorylated sites of metabolically 32P-labeled osteopontin from cultured chicken osteoblasts. J Biol Chem. 1997;272:13966–73.

    Article  CAS  PubMed  Google Scholar 

  93. Lienhard G. Non-functional phosphorylations? Trends Biochem Sci. 2008;33:351–2.

    Article  CAS  PubMed  Google Scholar 

  94. Ishikawa HO, Xu A, Ogura E, Manning G, Irvine KD. The Raine syndrome protein FAM20C is a Golgi kinase that phosphorylates bio-mineralization proteins. PLoS One. 2012;7, e42988.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  95. Tagliabracci VS, Engel JL, Wen J, Wiley SE, Worby CA, Kinch LN, et al. Secreted kinase phosphorylates proteins that regulate biomineralization. Science. 2012;336:1150–3.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Gaby Sowa for preparing the capillary columns for LC-MS/MS, Korbinian Mayr and Igor Paron for keeping the mass spectrometers in excellent condition, Markus Räschle for help with setting up Bioedit programs, and Ryouka Kawahara-Miki (Genome Research Center of the Nodai Research Institute) for his response to our questions.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Karlheinz Mann.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KM conceived the study, performed sample preparation and data acquisition. MM supplied methodological expertise. Both authors took part in the design of the study and were critically involved in manuscript drafting. Both authors read and approved the final manuscript.

Additional files

Additional file 1: Table S1.

This table is derived from the MaxQuant output table ProteinGroups and contains the complete list of identified proteins/protein groups including the ones that we did not accept after application of criteria described in Materials and Methods. This table also contains additional data such as the complete set of accession numbers forming one group, the distribution of peptides among the 20 PAGE sections analyzed separately, the calculated molecular weight of each entry, the iBAQ intensity and the percentages calculated from it. (XLSX 265 kb)

Additional file 2: Table S2.

Table S2, derived from MaxQuant output table Peptides, contains all identified peptide sequences belonging to the proteins of Table S1 in alphabetic order starting with the first amino acid: It also contains other relevant data such as scores, charge states, and mass accuracy. (XLSX 2472 kb)

Additional file 3: Table S3.

All accepted identifications with annotations derived from FASTA searches. Proteins and protein groups predicted to belong together are indicated by alternate white and grey background. The table also contains UniProt accession numbers of similar chicken and turkey proteins with FASTA-derived percentages of sequence identity and e-values, number of peptides and iBAQ-derived abundances. (DOCX 419 kb)

Additional file 4: Figure S1.

Alignments of quail accessions predicted to contain sequences (or partial sequences) of different proteins with the most similar chicken or turkey protein sequences. Also included is the location of identified peptides. (DOCX 143 kb)

Additional file 5: Table S4.

Curated version of the MaxQuant output file Phospho (STY) sites with phosphopeptide sequences, localization probabilities, scores, mass errors and other relevant data. Also included are phosphopeptides with site probabilities lower than 0.75, which were not included into Table 2. (XLS 99 kb)

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mann, K., Mann, M. Proteomic analysis of quail calcified eggshell matrix: a comparison to chicken and turkey eggshell proteomes. Proteome Sci 13, 22 (2015). https://doi.org/10.1186/s12953-015-0078-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12953-015-0078-1

Keywords