- Research article
- Open Access
Molt-dependent transcriptomic analysis of cement proteins in the barnacle Amphibalanus amphitrite
© Wang et al. 2015
- Received: 18 August 2015
- Accepted: 8 October 2015
- Published: 24 October 2015
A complete understanding of barnacle adhesion remains elusive as the process occurs within and beneath the confines of a rigid calcified shell. Barnacle cement is mainly proteinaceous and several individual proteins have been identified in the hardened cement at the barnacle-substrate interface. Little is known about the molt- and tissue-specific expression of cement protein genes but could offer valuable insight into the complex multi-step processes of barnacle growth and adhesion.
The main body and sub-mantle tissue of the barnacle Amphibalanus amphitrite (basionym Balanus amphitrite) were collected in pre- and post-molt stages. RNA-seq technology was used to analyze the transcriptome for differential gene expression at these two stages and liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) was used to analyze the protein content of barnacle secretions.
We report on the transcriptomic analysis of barnacle cement gland tissue in pre- and post-molt growth stages and proteomic investigation of barnacle secretions. While no significant difference was found in the expression of cement proteins genes at pre- and post-molting stages, expression levels were highly elevated in the sub-mantle tissue (where the cement glands are located) compared to the main barnacle body. We report the discovery of a novel 114kD cement protein, which is identified in material secreted onto various surfaces by adult barnacles and with the encoding gene highly expressed in the sub-mantle tissue. Further differential gene expression analysis of the sub-mantle tissue samples reveals a limited number of genes highly expressed in pre-molt samples with a range of functions including cuticular development, biominerialization, and proteolytic activity.
The expression of cement protein genes appears to remain constant through the molt cycle and is largely confined to the sub-mantle tissue. Our results reveal a novel and potentially prominent protein to the mix of cement-related components in A. amphitrite. Despite the lack of a complete genome, sample collection allowed for extended transcriptomic analysis of pre- and post-molt barnacle samples and identified a number of highly-expressed genes. Our results highlight the complexities of this sessile marine organism as it grows via molt cycles and increases the area over which it exhibits robust adhesion to its substrate.
- Cement protein
- Molt cycle
- Amphibalanus amphitrite
Adhesion related to biological processes is often important to maintain proper function, and this is especially true of sessile marine organisms where substrate adhesion is critical for survival. A prominent example is the barnacle, which first relies on temporary adhesion on a suitable substrate as a cyprid and then transitions to permanent attachment as the animal undergoes metamorphosis to a sessile juvenile and matures to an adult with either a calcareous or membranous base attached to the underlying surface . Barnacles are arthropods in the Crustacea subphylum and are partly characterized by their growth via molt cycles. Molting in adult acorn barnacles is evidenced by shedding of the exuvia as well as the formation of radial ecdysal lines underneath the animal. Much effort has been dedicated to understanding the intertwined processes of growth and adhesion as well as the composition and the molecular mechanisms that lead to a strongly adhered organism. These efforts are complicated by the presence of a rigid outer calcareous shell, under which cement components are delivered and deposited by the barnacle under the leading edge of the barnacle/substrate interface and sheltered from the external marine environment. As with growth, evidence suggests the process of adhesion occurs in a cyclic, multi-step manner synchronized with the molt cycle [2, 3]. The means by which cement precursor materials are delivered to the leading edge of the barnacle/substrate interface are not completely understood, yet a network of capillaries originating from cement glands and terminating at ducts at the substrate interface are thought to play a prominent role [2–4].
Barnacle cement itself is composed predominantly of protein [5–7], though other materials, such as lipids, polysaccharides and biomineralization components, have also been found at the barnacle-substrate interface . Several putative cement proteins have been identified and studied, mainly through solubilization and analysis of cement plaques that remained on substrates following removal of the barnacle shell and body [7–12]. These include cp100k, cp68k, cp52k, cp20k, cp19k, where the number corresponds with the protein molecular weight in daltons. Detailed structural information about each protein remains unknown, though amino acid composition and conserved and/or repeat subunits have shed light on potential function [9, 12, 13]. Additionally, advances in gene sequencing technology have aided in the development of cDNA libraries as well as the discovery of genes that encode cement protein homologues in different barnacle species [14–16]. This knowledge continues to fill in the picture of how barnacles are able to create a robust adhesive interface in a relatively harsh marine environment.
With this backdrop, we set out to test the hypothesis that cement proteins are actively regulated at different stages of the molt cycle and in different tissues of adult barnacles. Molting of several adult Amphibalanus amphitrite (basionym Balanus amphitrite) barnacles was tracked based on shedding of the exuvia under controlled laboratory conditions. Protein and mRNA were collected from the sub-mantle tissue (containing the cement glands) of barnacles in pre- and post-molt stages. Material was also collected from the main body of adult barnacles for comparison. Using RNA-Seq technology, the expression levels of cement protein genes in the sub-mantle tissue of pre- and post-molt samples were found to be the same, suggesting these proteins are not regulated at the transcript level when synchronized to the molt cycle. Comparison of cement protein gene expression in the sub-mantle tissue versus the main barnacle body revealed significantly elevated levels in sub-mantle tissue. These results fit with the current understanding of barnacle anatomy, where cement pre-cursor material is thought to be produced and stored in cement glands located in the sub-mantle tissue. We also report the discovery of a novel 114kD protein, named Aacp114k, and found its transcript to be highly expressed in the sub-mantle tissue and the protein identified in barnacle cement secretions collected at the barnacle/substrate interface. This protein shares significant homology to a previously identified 100kD cement protein in A. amphitrite as well as others identified in the barnacle species Megabalanus rosa and Teraclita japonica formosana [7, 15, 16].
In the absence of a complete genome, the information collected from the comparison of pre- and post-molt samples allowed for a more general analysis of differential gene expression in these two stages. A number of genes were identified as being up-regulated in the sub-mantle tissue predominantly in pre-molt barnacle samples with a range of functions including cuticular development, biominerialization, and proteolytic activity. Collectively, these results add a novel and potentially prominent protein to the mix of cement-related components and further our understanding of the complex processes at play as barnacles grow via cyclic molting and extend their adhesive area to substrates.
Sample collection and molt stage confirmation
RNA-seq transcriptome analysis
To understand barnacle cement protein gene expression levels in pre- and post-molt growth stages, RNA-seq technology was used to compare transcriptomic profiles of pre- and post-molt barnacles, specifically examining the sub-mantle tissue containing cement gland tissue. We first provide a general summary of these results. Complementary DNA (cDNA) libraries were constructed from three independent samples of sub-mantle tissue collected in pre-molt or post-molt stages. Two additional cDNA libraries were made from two pooled main body samples as pilot experiments for a total of eight separate samples.
Summary of A. amphitrite transcriptome
Average length of each read (bp)
Total assembled sequence size (bp)
GC content (%)
Contig mean length (bp)
Contig N50 (bp)
Transcripts annotated by SwissPro
Transcripts with GO term
Cement protein gene expression: pre-molt versus post-molt
Comparison of cement proteins in A. amphitrite
GenBank Accession Number
KP863707 (this study)
KP863709 (this study)
KP863708 (this study)
Cement protein gene expression: main body versus sub-mantle tissue
In addition to examining cement protein gene expression levels, material secreted from barnacles at the substrate interface was analyzed using mass spectrometry. There is not a standard method to collect barnacle secretions, so different approaches were employed to identify the presence of cement proteins and attempt to correlate their presence in the sub-mantle tissue and the barnacle/substrate interface. See Methods for detailed explanation of the coverslip, medallion, mortar, and bead collection techniques.
Number of unique peptides identified for cp100k
Extended analysis of differential gene expression: Pre- and post-molt
Of the 135 differentially expressed transcripts in the pre- and post-molt samples (Additional file 6), proteins encoded by 27 of them (20 %) were identified in the secretion samples studied by proteomics, as highlighted in Table 4. Additional file 7 provides detailed information about protein identifications from these samples. Interestingly, one of the encoding genes highly expressed in pre-molt samples, comp97185_c0 (Additional file 6) shares 44 % sequence homology to the settlement inducing protein complex (SIPC) of both A. amphitrite  and M. rosa . Peptides of this protein were recovered and identified in all four barnacle secretion collections.
In this study, we compare gene expression levels of pre- and post-molt barnacles while also examining the main body versus the sub-mantle tissue. The first comparison was enabled by tracking the shed exuvia of several barnacles and collecting sub-mantle tissue at time intervals preceding or following molt. The molt stage was confirmed by examining cirri anatomy via light microscopy and correlated with established stages of barnacle molting . The two comparisons presented in this study have resulted in a thorough analysis of principal genes differentially expressed in pre- and post-molt stages and also a more complete picture of the abundance and presence of putative barnacle cement proteins involved in adhesion to the underlying surface of an adult barnacle.
The assembled transcriptome size (~114 Mbp) in this study is more than four times greater than the number previously reported in A. amphitrite by Chen et al. . In RNA-Seq analysis of differential gene expression in different conditions, it is very common that de novo assembling of eukaryotic transcriptome often reports fragmented transcripts of one gene or repeated contigs that only differ by a single-nucleotide polymorphism (SNP) or a base insertion or deletion (indel). Such expanded transcriptome size is usually evidence of isoforms of the same gene and artifacts such as sequencing errors, repeats, and variations in coverage and genetics within diploid individuals and pooled samples . In our study, RNA of the pre- and post-molt samples was from six individuals and RNA of the whole bodies was collected from two pooled samples, each containing ten individual barnacles. The transcriptome generated from those eight samples resulted in many paths in the graph construction during the assembly process, effectively fragmenting the assembly which ultimately produced 163,929 assembled contigs. The main contribution to the larger transcriptome size might be heterogeneities within input samples, which was indicated by following evidence: (1) when sequence reads from the pre- and post-molt conditions were assembled separately, the contig numbers were greatly reduced to 73,973 and 70,840, respectively, and the total transcriptome sizes were reduced to approximately 42 Mbp (see Additional file 8). (2) Further clustering these individual assemblies using a homology approach (CD-HIT, 95 % homology) resulted in an overall reduced assembly of approximately 54.8 Mbp in 95,873 contig clusters (35 % reduction from combined). Clustering the individual assemblies results in ca. 60 k contig clusters in each condition, which might be closer to the expected gene number of the barnacle. The fact that the contig sequences from the pre- and post-molting samples do not cluster to a greater degree at 95 % suggested that variation (splicing, SNP, sequencing error) played a greater role in the transcriptome assembly. (3) Finding a six nucleotide difference between the assembled and the cloned Aacp114 gene sequences (data not shown) also indicated the heterogeneity phenomenon present within individuals. Further inspecting the reads mapping to the Aacp114k gene, we observed many differences between the consensus sequences not shared across individual samples (data not shown).
The Illumina sequencing approach employed in the current study produced sequence coverage at ~114x, which is an order of magnitude higher than 8.6x coverage using the 454 method reported by Chen et al. . Higher depth of sequencing would reveal more real alternative splicing transcripts and splicing variants, which contribute to greater numbers of assembled sequences and contigs. On the other hand, higher coverage may also amplify issues with genomic contamination, partially spliced products, and other artifacts of the method. In this study, we were only interested in genes changing within the conditions queried, negating many of these issues in the final differential analysis. As such, reducing the putative transcript/gene (cluster) tally was not necessary. Here, only about 23 % of the ~160,000 predicted transcripts had annotated functions based on available database sequences, meaning the vast majority are unidentified. This trend was consistent even among the differentially expressed proteins in the pre- and post-molt tissue samples, indicating the need for a much larger effort of creating a fully annotated barnacle genome.
The main purpose of targeting the sub-mantle tissue for temporal comparative analysis is that, anatomically, this region is where cement-related tissue of adult barnacles resides. While it was straightforward to separate main body from sub-mantle tissue with little cross-contamination, collection of only cement gland tissue within the sub-mantle cavity proved more challenging. Ovarian tissue, identified by its distinct color and texture, was often intimately intertwined with the surrounding tissue, which included the cement glands [25, 26]. During sample collection, care was taken to remove the vast majority of ovarian tissue, though its complete removal was beyond the scope of this investigation given the limits of dissection. Therefore, potential contamination of a minimal amount of ovarian tissue in the sub-mantle samples is a possibility. Despite this, the aim was to collect as much of the non-ovarian sub-mantle tissue as possible in order to have an optimal sampling of material thought to contain cement precursor materials. Judging by the gene expression data comparison between sub-mantle tissues and the main body, the sub-mantle tissues clearly demonstrates elevated expression of cement proteins and resulted in the identification of a new cement-related protein.
In regard to temporal expression of cement protein genes, there have been previous claims of differential cement protein gene expression as a function of barnacle molting [9, 13], however, these are unsubstantiated and cite unpublished work. In the current study, several barnacle samples with regular molt cycles were chosen in either a pre- or post-molt state accompanied by collection of cirri for molt cycle confirmation based on established criteria . For the detected cement protein genes (Fig. 3), there was no significant difference in the pre- and post-molt gene expression levels despite that fact that several other genes were found to be differentially expressed (Table 4). It may be possible to increase the resolution of molt-dependent gene expression using qRT-PCR according to eight (sub-)stages described by Davis et al. . However, differential gene expression was intentionally limited to pre-molt and post-molt due to the considerable number of other parameters that, collectively, are difficult to control and monitor. These include barnacle size, age, and season, time variability within each molt stage and transcript heterogeneity. Despite these variables, two phenotypes were observed in the present study and were consistent among themselves.
The presence of cement proteins Aacp20-1 and −2 has been reported previously and were shown to localize and express in different types of secretory cells in the cyprid cement gland . In adult barnacles, He et al. reported that protein Aacp20k-1 was present in soft tissue, whereas protein Aacp20k-2 was dominant in the basal shell of adult barnacles, which agrees with our gene expression results (Fig. 4). Interestingly, Aacp20k-1 is the only cement protein with a low pI of ~5 (whereas the pI values of all other cement proteins are ~9 or greater, see Table 2) and its gene expression was detected at very low levels in the sub-mantle tissue and not at all in the body. It is possible that transcripts of Aacp20k-1 are not stable and were degraded, making it difficult to assess if these results are coincident. Nevertheless, the questions of how and why this transcript is regulated are of interest.
We report high levels of gene expression of a novel paralogue cp100k protein termed Aacp114k. Irrespective of the pre- or post-molt stage, this newly identified protein had the highest gene expression compared to all other putative cement proteins located in the sub-mantle tissue based on transcript numbers (Fig. 3). These results and the relative change in expression level in the sub-mantle tissue versus body (Fig. 4) suggest that this protein is localized in sub-mantle tissue. Proteomic analysis confirms this protein is found in material secreted by the barnacle to a greater extent than the Aacp100k paralogue, as observed in coverslip, medallion, microsphere and mortar samples (Table 3). Based upon our gene expression and mass spectrometry results, we suggest the Aacp114k paralogue identified in the current study is a prominent protein component of barnacle cement in A. amphitrite. Interestingly, phylogenetic analysis of the various cp100 homologues show a closer relation between the cp100 homologue in A. amphitrite first identified by Chen et al.  and in T. j. formosana  than between the two paralogues found in A. amphitrite (Additional file 3).
A previous report on M. rosa claimed cp52k, cp68k, and cp100k as being the major protein components of barnacle cement . However, we were not able to identify a cp68k-related gene from the transcriptome in A. amphitrite mainly because the cp68k sequence is unavailable in public domains (such as Genbank, SwissPro) despite references to it having been identified and its sequence [9, 13]. Using N-terminal sequences of eight cyanogen bromide (CB) peptides determined from the insolubilized cement fraction of M. rosa  as queries, we found that only five of those peptides were homologous to either Aacp52K or Aacp100K in our transcriptome database. Therefore, more studies are required to verify whether cp68k exists as single protein in A. amphitrite and, if so, whether the expression level is differentiated throughout molting.
In comparing the main body to the sub-mantle tissue using RT-PCR, cement protein genes were predominantly expressed in sub-mantle tissue, which is the expected location of barnacle cement precursor materials. Proteomic analysis of barnacle secretions collected using different techniques indicate that Aacp114k and Aacp100k are secreted by the animal at the barnacle/substrate interface. Of course, the inability to detect other putative cement proteins from the secreted materials does not necessarily mean they are not present. Two factors could contribute to this finding. First, the amount of secreted material is expected to be relatively low since the collection methods used in this study are targeted at newly secreted material over a limited time period with a maximum of 2 weeks and minimum of 1 day. (This is in comparison to collection of an entire, fully-formed cement plaque). In addition, MS analysis may lie outside the limits of detection or may indicate the proteins have been cross-linked prior to enzymatic digestion, escaping detection. Combined, these factors are likely to contribute to the limited detection of cement proteins in the various barnacle secretion collections. Nevertheless, the newly described variations to the collection methods of cement material over relatively short time frames using four different approaches (the coverslip, mortar, medallion, and bead techniques) allowed for proteomic comparison with sub-mantle tissue and offer a fresh means of analyzing recently secreted material and newly formed barnacle cement.
We expand upon the insights gained from cement protein (gene) detection and analysis to the more general findings of differential gene expression in the pre- and post-molt samples of the sub-mantle tissue. Most highly expressed genes were found in pre-molt samples and, of those identified, generally encoded proteins in four categories: cuticular proteins, metalloproteases, C-type lectins, and antioxidants. These categories can be associated with molting, biomineralization, morphogenesis, and cell adhesion. Here, we provide a more in-depth discussion of some of the highly expressed protein genes identified in the sub-mantle tissue of pre-molt barnacle samples.
The increased expression of cuticular protein genes was not unexpected since molting requires significant production of cuticular material and tissue scaffolding to replace the soon-to-be-shed exuvia. Several cuticle-related protein genes in the pre-molt samples were identified and four were found to have the chitin-binding Rebers-Riddiford consensus motif . Multiple transcripts of RR-containing cuticle proteins were reported to be temporally and spatially expressed in the decapod crustacean Callinectes sapidus and involved in biomineralization . More than 60 homologs of this type of cuticle protein gene were identified from the transcriptome data in this study, indicating a variety of functions at different growth stages or conditions. RR-containing cuticle protein genes highly expressed in the pre-molt sub-mantle tissue of barnacles strongly suggest an involvement with chitin and possibly a role in biominerialization, both of which have a role in barnacle growth .
Two genes encoding C-type lectins, BRA-3 and nattectin, were also found to be highly expressed in the pre-molt condition. C-type lectins usually bind carbohydrate structures in a Ca2+-dependent manner and are well-recognized immune response initiator proteins . They have been found to participate in mineralization  and inhibit the growth of calcium carbonate crystals . High expression of these two genes in the pre-molt condition suggests that they have significant roles in biomineralization during growth/molting and coordinating post-molt responses.
Coincident differential expression of genes encoding peroxiredoxin-6 and thioredoxin domain-containing protein 5 in pre-molt samples suggests that the barnacle experiences oxidative stress during or following molting, since peroxiredoxin uses thioredoxin to recharge after reducing hydrogen peroxide (H2O2) so as to restore its catalytic activity . Antioxidant proteins like peroxiredoxin are often coupled with others (like thioredoxin) to limit damage of oxidative stress to host tissue .
Increased expression of genes encoding proteases (including zinc metalloproteinases) in the pre-molt stage points toward proteolytic activation of proteins, including contributors to barnacle cement, having an essential role for effective cement cross-linking and co-adhesion within the interface. Zinc metalloproteinases have a significant role in inflammation as enzymes that rework extracellular matrices and remodeling of structural material . In addition, it has previously been suggested that barnacle cement cross-linking involves proteolytic activation of structural proteins in order to maximize the potential for bonding interactions with other proteins and with the surface .
A glycoprotein, settlement inducing protein complex (SIPC) has been reported as the species-specific cue to settlement of barnacles [22, 23] and expressed in the cuticle of adults . It has been reported as a complex of three proteins (98, 88 and 76 kDa) but retains some biological activity as separate subunits [37, 38]. In the transcriptomic data, we found SIPC that is 100 % identical to the previously reported one  and its gene expressed at similar levels in both pre- and post-molt conditions. Interestingly, we also found a homologous SIPC gene, comp97185_c0 (Additional file 6), highly expressed in the pre-molt condition and its encoded peptides were identified in all secretion samples analyzed with MS. This finding suggests that the barnacle may secrete a new and different type or combination of SIPC molecule with alternative (non-SIPC) ligands as a molting signal, and its function is worth further investigation.
This work advances the understanding of the complex process of barnacle molting/growth by defining the temporal and tissue distribution of various proteins, including those reported to be major components of barnacle cement. Our analysis suggests there is not differential expression of cement protein genes in pre- and post-molt barnacles and most appear to be highly expressed in the sub-mantle tissue compared to the main body. In addition, we identify a novel cement protein homologue, Aacp114k, with high gene expression levels in the sub-mantle tissue and confirmed its presence in secretions at the barnacle/substrate interface. In the absence of a fully annotated genome, our transcriptomic analysis reveals highly expressed genes identified in pre-molt barnacle samples and paints a picture of extensive cuticular development, biominerialization, and proteolytic activity that may be associated with the precursor elements of barnacle cement. Our results also highlight the fact that the majority of differentially expressed genes in the pre- and post-molt stage are unidentified and potentially unique to barnacles.
The description of several new cement collection techniques opens opportunity to build a profile of the macromolecules present and active at the barnacle/substrate interface with the goal of elucidating the mechanism of the cementing process. Future directions include developing approaches to collect and profile pre-secreted material as well as establishing protocols to identify the location of various cement components prior to secretion. We note that in situ hybridization studies have been carried out in A. amphitrite cyprids  and offer great potential for identifying the tissue (or cell) specific expression of particular genes relevant to the current study. Given the increased complexity of the adult barnacle anatomy and the fact that all tissue is surrounded by an opaque calcified shell, there are challenges with in situ hybridization that are beyond the scope of the current work yet are a promising future direction.
Ultimately, the research goals are to identify the various precursor components of barnacle cement, understand the mechanism of delivery, and how the components interact at the substrate interface. A complete understanding of this process can lead to targeted strategies to mitigate barnacle adhesion as well as mimic its functionality in underwater environments.
A. amphitrite cyprids were settled and raised on silicone coated panels as previously described  until they were mature enough to be fed Artemia spp. nauplii. Samples were then shipped to the Naval Research Laboratory where they were maintained in an incubator at 23 °C with 12 h daylight cycles in 32 ppt artificial seawater and fed Artemia spp. nauplii every other day. In addition, the water for barnacles was changed weekly during which the algal growth was removed. Barnacles were reattached from the silicone panels to alternative substrates for the experiments .
Molt cycle analysis and cirri imaging
A subset of over ten barnacles were reattached to glass slides, kept individually in separate containers to track the molt cycles over the course of several weeks (minimum of 2 weeks). Molting was noted by the presence of the exuvia in an individual barnacle container and, depending on the age or size of the barnacles, occurred anywhere from every 3 to 7 days. Barnacles molting in a regular pattern were selected for this study. For pre- and post-molt comparison, three barnacles were sacrificed the day prior and three the day after a molt, for a total of six individual samples. The cirri were removed for light microscopy observation (see below). In other barnacles undergoing regular molt cycles, only the main body was collected for subsequent protein and mRNA analysis and the cirri were removed for light microscopy analysis.
Cirri were examined using differential interference contrast (DIC) light microscopy on a Nikon TE2000 inverted microscope to confirm the molting stage of each sample. Following the procedure of Davis et al. , setae from the upper portion of cirri (rami) were imaged under 20× and 40× magnification. Several segments were imaged for consistency and structural features identifying the molt stage were noted. Davis et al. categorize four main stages, A through D, of barnacle molting proceeding through a short time period following molting (stage A) to the period prior to molting (stage D). These main stages are discriminated further into eight sub-stages. For the purposes of this study we have generalized these stages to two main categories: i) pre-molt (proecdysis, stage D) and ii) post-molt (post- or interecdysis, stages B or C).
Total RNA and protein collection
Isolation of total RNA from whole body tissue without the sub-mantle tissue was performed in the following manner. First, the main body was removed from a barnacle shell that had been cracked open with sharp tweezers followed by the addition of TRIzol Reagent (Life Technologies Inc., Grand Island, NY) according to the manufacturer’s protocol. The soft body tissues in TRIzol were vigorously vortexed and resuspended several times to break them apart. Chloroform was then added and mixed evenly by vortexing. The mixture was allowed to sit at room temperature for 2–3 min and centrifuged. Subsequently, the liquid phase was removed without disturbing the interphase precipitate and mixed with isopropanol to precipitate total RNA, followed by centrifugation and washing the pellet with 75 % ethanol. The rest of liquid phase including the interphase precipitate was then subjected to DNA and protein isolation according to the manufacturer’s protocol. The total RNA was further treated with DNAse I to get rid of DNA contaminants using on-column digestion (Qiagen Inc., Valencia, CA).
A similar procedure was used for isolating total RNA from sub-mantle tissue. The main body of the animal was removed as well as the underlying inner mantle and the sub-mantle tissue carefully collected with several gentle pipette suctions for processing. A schematic of the barnacle anatomy highlighting the main body and sub-mantle tissue is shown in Fig. 1. Ovarian tissue, which was clearly evident based on its color and texture with respect to surrounding tissue, was carefully removed with tweezers to minimize/eliminate its inclusion for subsequent processing/analysis.
For all samples, isolated total RNA were examined by agarose gel electrophoresis. A quality check was performed using a Bioanalyzer (Agilent Technologies) and RNA samples were found to be of high quality with no DNA contamination.
mRNA sequencing and data analysis
Three separate samples of the sub-mantle tissue from the pre- and post-molting conditions and two pooled samples from the whole body (10 samples in each; all without the sub-mantle tissue) were collected for RNA sequencing (RNA-seq) experiments. Enrichment of poly(A) RNA, synthesis of cDNA and construction of cDNA libraries were processed using protocols provided by Illumina, Inc. (San Diego, CA). Libraries were sequenced on an Illumina HisSeq2000 Sequencer in 100 paired-end configuration. Following sequencing, data was trimmed for both adaptor and quality using a combination of EA-utils and Btrim [41, 42]. Trimmed-paired reads were then de novo assembled using Trinity software using the default parameters found in the online protocol . Assembled transcripts were also further clustered into putative gene clusters using CD-HIT-EST . Annotation of transcripts was performed using the Transdecoder/Trinotate package [http://www.vcru.wisc.edu/simonlab/bioinformatics/programs/trinity/docs/annotation/Trinotate.html] (blast based annotation using Swissprot) with the default E-value as 1e-6. Read alignment and estimation of transcript abundance was performed using Bowtie2  and RSEM  as described in the Trinity abundance estimation protocol. Transcripts differentially expressed between two conditions were identified by comparing absolute number of transcripts across biological replicates with FDR < 0.1 as the cutoff value using the EdgeR software package , again as described in the Trinity protocol. Raw sequence reads were deposited in Genbank as BioProject PRJNA271096.
The primers used for the quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) assays are listed in Additional file 4. The cytochrome b gene (cyb) of A. amphitrite was chosen as an internal control. The qRT-PCR was conducted on an iCycler (BioRad Laboratories) using the iScript OneStep RT-PCR Kit with SYBR Green (BioRad Laboratories). qRT-PCR reaction mixtures consisted of 10 ng of total RNA from each extracted sample, 200 nM of primers and were subjected to the following cycling conditions: one cycle at 50 °C for 30 min and 95 °C for 15 min, followed by 40 cycles of 94 °C for 15 s, 52 °C for 30 s and 72 °C for 30 s. Triplicates for each RNA sample were collected for statistical analysis. The average cycle threshold (Ct) value for each tested RNA sample was obtained and the gene expression level was then used for comparison among tissue samples. Relative quantities of transcript were determined using the 2-ΔΔCt formula. ΔCt is the difference in Ct of the gene of interest and Ct of the normalizer gene, and ΔΔCt is the difference in ΔCt from the sub-mantle sample and ΔCt from the body sample. Cytochrome b gene (cyb) was used to normalize expression levels of genes of interest because its transcript was constant in all conditions in this study.
Cloning and sequencing of Aacp114k encoding gene
Total RNA isolated from sub-mantle tissues were reverse transcribed to cDNA according to the manufacturer’s protocol (Life Technologies). The first strand cDNA was then used as a template along with the primers, 114KNcoIF (5′AAAAACCATGGGCATGCTGCGGCTCTCGCTAGC3′) and 114KXhoIR (5′GGGGGCTCGAGGCACTTGAAGTAGTCGTAC3′), to amplify the Aacp114k transcript by conducting PCR. Resulting 3 kb DNA fragments were digested with NcoI and XhoI, followed by cloning into a pET22b vector. The cloned plasmids were propagated in E. coli Top10 in LB culture supplemented with ampicillin. The insert of isolated plasmid was sequenced using the Sanger sequencing method (Eurofins Genomics).
Barnacle secretion collection
There is not a standard procedure to collect barnacle cement and secretions, so to maximize the opportunity to identify potential cement proteins collection of barnacle secretions was carried out in three different ways to produce four different samples for subsequent analysis. From each collection method, the majority of protein was solubilized. Any insoluble fraction was minimal and remained on the respective surface. A description of each technique is provided below.
In the first, adult barnacles were gently removed from silicone substrates and all shell plates, including the base plate, were rinsed and cleaned with a cotton swab in deionized water. The barnacles were then inverted and settled onto a Plexiglas® platform engineered with circular openings, which allowed the apical portion of the organism to be partially submerged in sea water to allow for continuous feeding and growth during the cement collection process. A No. 2, 12 mm diameter borosilicate cover glass (Fisher Scientific #22-293-232) was placed on top of the exposed base plate as a substrate to collect secreted material from the barnacle (see Additional file 9 for schematic of setup). After 24–48 h at room temperature, the cover glass substratum was gently removed using a dissecting needle. In order to maximize sample collection, the entire cover glass was placed into a 1.65 mL microfuge tube and crushed with a glass rod. Next, 100 μL 1× Laemmli sample buffer was directly added to the sample and heated for three minutes at 95 °C. Finally, 50 μL of the solubilized sample was separated by SDS-PAGE (4-20 % TGX precast gel; BioRad) and used for further mass spectrometry analysis. These samples are referred to as “coverslip” samples.
Mortar and medallions
In the second, adult barnacles were transferred from silicone panels and settled on sodium aluminoborate (Na2O • Al2O3 • 3B2O3) glass substrates, which form a hydrated reaction layer (<25 μm thick) in aqueous environments that is resistant to barnacle adhesion . After 2 weeks, barnacle bodies and side plates were carefully removed, leaving the base of the barnacles attached to the substrates which were then cleaned with a cotton swab in deionized water to remove loosely bound organic matter (see Additional file 9). Base plates were then de-minerialized by immersing substrates in 0.1 M ethylenediaminetetraacetic acid (EDTA) at room temperature for 48–72 h. The organic remnants of the base plates were collected and placed in 50 μL of 0.1 % (w/v) sodium dodecyl sulfate (SDS) solution in deionized water for subsequent analysis; these are referred to as “mortar” samples. Likewise, barnacle “medallions,” consisting of the cuticular layer and underlying barnacle secretions, were gently rinsed with deionized water then peeled off the aluminoborate glass substrates and placed in 50 μL of the SDS solution. Next, 50 μL of Laemmli sample buffer containing 300 mM dithiothreitol (DTT) was added to both samples and the resulting suspensions were incubated for 15 min at 95 °C to solubilize proteins in the samples. After centrifugation (9000 rpm for 1 min), supernatants were separated by SDS-PAGE for further analysis. This procedure resulted in the cuticle being pelleted to the bottom of the centrifuge tube.
In the third type of preparation, borosilicate glass microspheres (<50 μm in diameter) were placed in the bottom of dishes filled with ASW to form a bed of microspheres 2–3 mm in depth. Adult barnacles (n = 5) were then transferred from silicone panels onto the bed of microspheres. After 24 h, barnacles were lifted out, and those with microspheres attached to their underside (three out of the five) were gently scraped to remove the majority of the microspheres without damaging the barnacles (see Additional file 9). The collected microspheres were placed in 50 μL of an SDS solution, vortexed, and centrifuged at 9000 rpm for 1 min. The resulting supernatant was collected, mixed in a 1:1 ratio with Laemmli sample buffer, and then separated by SDS-PAGE for further analysis and referred to as “bead” samples.
Individual bands from protein extracts of each sample separated on SDS-PAGE gels were excised and digested in gel by trypsin. Peptides were extracted by 2 % formic acid in 50/50 acetonitrile/water, followed by 100 % acetonitrile. Digests were analyzed by liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) using Tempo-MDLC coupled to a QStar Elite mass spectrometer (AB Sciex, Foster City, CA). Tandem mass spectra were extracted by AB Sciex MS data convertor version 2. Charge state deconvolution and deisotoping were not performed. All MS/MS samples were analyzed using Mascot (Matrix Science, London, UK; version 2.4.1) and X! Tandem (The GPM, thegpm.org; version CYCLONE (2010.12.01.1)). Mascot was set up to search the BALAM_adhesproteins_003 database (1492 entries) assuming the digestion enzyme trypsin. X! Tandem was set up to search a subset of the BALAM database also assuming trypsin. Mascot and X! Tandem were searched with a fragment ion mass tolerance of 0.80 Da and a parent ion tolerance of 0.80 Da. Deamidation of asparagine and glutamine, oxidation of methionine, acetyl of the N-terminus and carbamidomethyl of cysteine were specified in Mascot as variable modifications. Glu- > pyro-Glu of the N-terminus, ammonia-loss of the N-terminus, gln- > pyro-Glu of the N-terminus, deamidated of asparagine and glutamine, oxidation of methionine, acetyl of the N-terminus and carbamidomethyl of cysteine were specified in X! Tandem as variable modifications. Scaffold (version Scaffold_4.2.1, Proteome Software Inc., Portland, OR) was used to validate MS/MS based peptide and protein identifications. Peptide identifications were accepted if they could be established at greater than 70 % probability by the Peptide Prophet algorithm  with Scaffold delta-mass correction. Protein identifications were accepted if they could be established at greater than 90 % probability and contained at least 2 identified peptides. Protein probabilities were assigned by the Protein Prophet algorithm . Proteins that contained similar peptides and could not be differentiated based on MS/MS analysis alone were grouped to satisfy the principles of parsimony. Proteins sharing significant peptide evidence were grouped into clusters.
Availability of supporting data
The data sets supporting the results of this article are available in the Proteomics IDEntifications (PRIDE) repository, http://www.ebi.ac.uk/pride/archive/projects/PXD001982. Other supporting data are presented in Additional files 1, 2, 3, 4, 5, 6, 7, 8 and 9 and also found in the same repository.
We gratefully acknowledge Prof. Dan Rittschof and Ms. Beatriz Orihuela from Duke University Marine Lab for supplying the juvenile and adult barnacle samples. T.E.B. and S.E.H were supported through the Naval Research Enterprise Internship Program and A.M. through an NRL summer internship. Funding for this work was provided by Naval Research Laboratory base funding and the Office of Naval Research.
Open AccessThis 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.
- Walker G. Adhesion of Barnacles. J Adhesion. 1981;12:51–8.View ArticleGoogle Scholar
- Saroyan JR, Lindner E, Dooley CA. Repair and Reattachment in Balanidae as Related to Their Cementing Mechanism. Biol Bull. 1970;139:333–50.View ArticleGoogle Scholar
- Burden DK, Barlow DE, Spillmann CM, Orihuela B, Rittschof D, Everett RK, et al. Barnacle Balanus amphitrite Adheres by a Stepwise Cementing Process. Langmuir. 2012;28:13364–72.View ArticlePubMedGoogle Scholar
- Burden DK, Spillmann CM, Everett RK, Barlow DE, Orihuela B, Deschamps JR, et al. Growth and development of the barnacle Amphibalanus amphitrite: time and spatially resolved structure and chemistry of the base plate. Biofouling. 2014;30:799–812.PubMed CentralView ArticlePubMedGoogle Scholar
- Walker G. Biochemical Composition of Cement of 2 Barnacle Species, Balanus-Hameri and Balanus-Crenatus. J Mar Biol Assoc UK. 1972;52:429–35.View ArticleGoogle Scholar
- Naldrett MJ. The Importance of Sulfur Cross-Links and Hydrophobic Interactions in the Polymerization of Barnacle Cement. J Mar Biol Assoc UK. 1993;73:689–702.View ArticleGoogle Scholar
- Kamino K, Inoue K, Maruyama T, Takamatsu N, Harayama S, Shizuri Y. Barnacle cement proteins - Importance of disulfide bonds in their insolubility. J Biol Chem. 2000;275:27360–5.PubMedGoogle Scholar
- Kamino K, Odo S, Maruyama T. Cement proteins of the acorn barnacle, Megabalanus rosa. Biol Bull. 1996;190:403–9.View ArticlePubMedGoogle Scholar
- Kamino K. Underwater adhesive of marine organisms as the vital link between biological science and material science. Mar Biotechnol. 2008;10:111–21.View ArticlePubMedGoogle Scholar
- Jonker JL, von Byern J, Flammang P, Klepal W, Power AM. Unusual adhesive production system in the barnacle Lepas anatifera: An ultrastructural and histochemical investigation. J Morphol. 2012;273:1377–91.View ArticlePubMedGoogle Scholar
- Jonker JL, Abram F, Pires E, Coelho AV, Grunwald I, Power AM. Adhesive Proteins of Stalked and Acorn Barnacles Display Homology with Low Sequence Similarities. Plos One. 2014;9:e108902.PubMed CentralView ArticlePubMedGoogle Scholar
- Nakano M, Kamino K. Amyloid-like Conformation and Interaction for the Self-Assembly in Barnacle Underwater Cement. Biochem. 2015;54:826–35.View ArticleGoogle Scholar
- Kamino K. Barnacle Underwater Cement. In: Smith AM, Callow JA, editors. Biological Adhesives. Berlin: Springer; 2006. p. 145–66.View ArticleGoogle Scholar
- De Gregoris TB, Borra M, Biffali E, Bekel T, Burgess JG, Kirby RR, Clare AS. Construction of an adult barnacle (Balanus amphitrite) cDNA library and selection of reference genes for quantitative RT-PCR studies. BMC Mol Biol. 2009; 10:doi:10.1186/1471-2199-1110-1162.
- Chen ZF, Matsumura K, Wang H, Arellano SM, Yan XC, Alam I, et al. Toward an Understanding of the Molecular Mechanisms of Barnacle Larval Settlement: A Comparative Transcriptomic Approach. Plos One. 2011;6:e22913.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin HC, Wong YH, Tsang LM, Chu KH, Qian PY, Chan BKK. First study on gene expression of cement proteins and potential adhesion-related genes of a membranous-based barnacle as revealed from Next-Generation Sequencing technology. Biofouling. 2014;30:169–81.View ArticlePubMedGoogle Scholar
- Davis CW, Fyhn UEH, Fyhn HJ. Intermolt Cycle of Cirripeds - Criteria for Its Stages and Its Duration in Balanus-Amphitrite. Biol Bull. 1973;145:310–22.View ArticleGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.PubMed CentralView ArticlePubMedGoogle Scholar
- Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.View ArticlePubMedGoogle Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.PubMed CentralView ArticlePubMedGoogle Scholar
- Rebers JE, Riddiford LM. Structure and Expression of a Manduca-Sexta Larval Cuticle Gene Homologous to Drosophila Cuticle Genes. J Mol Biol. 1988;203:411–23.View ArticlePubMedGoogle Scholar
- Dreanno C, Matsumura K, Dohmae N, Takio K, Hirota H, Kirby RR, et al. An alpha(2)-macroglobulin-like protein is the cue to gregarious settlement of the barnacle Balanus amphitrite. P Natl Acad Sci USA. 2006;103:14396–401.View ArticleGoogle Scholar
- Yorisue T, Matsumura K, Hirota H, Dohmae N, Kojima S. Possible molecular mechanisms of species recognition by barnacle larvae inferred from multi-specific sequencing analysis of proteinaceous settlement-inducing pheromone. Biofouling. 2012;28:605–11.View ArticlePubMedGoogle Scholar
- Davidson NM, Oshlack A. Corset: enabling differential gene expression analysis for de novo assembled transcriptomes. Genome Biol. 2014;15:410.Google Scholar
- Lacombe D. A Comparative Study of Cement Glands in Some Balanid Barnacles (Cirripedia-Balanidae). Biol Bull. 1970;139:164–79.View ArticleGoogle Scholar
- Walker G. Histochemistry and Ultrastructure of Cement Apparatus of 3 Adult Sessile Barnacles, Elminius-Modestus, Balanus-Balanoides and Balanus-Hameri. Mar Biol. 1970;7:239–48.View ArticleGoogle Scholar
- He LS, Zhang G, Qian PY. Characterization of Two 20 kDa-Cement Protein (cp20k) Homologues in Amphibalanus amphitrite. Plos One. 2013;8:e64130.PubMed CentralView ArticlePubMedGoogle Scholar
- Wynn A, Shafer TH. Four differentially expressed cDNAs in Callinectes sapidus containing the Rebers-Riddiford consensus sequence. Comp Biochem Phys B. 2005;141:294–306.View ArticleGoogle Scholar
- Robinson MJ, Sancho D, Slack EC, LeibundGut-Landmann S, Sousa CRE. Myeloid C-type lectins in innate immunity. Nat Immunol. 2006;7:1258–65.View ArticlePubMedGoogle Scholar
- Kamiya H, Jimbo M, Yako H, Muramoto K, Nakamura O, Kado R, et al. Participation of the C-type hemolymph lectin in mineralization of the acorn barnacle Megabalanus rosa. Mar Biol. 2002;140:1235–40.View ArticleGoogle Scholar
- Muramoto K, Yako H, Murakami K, Odo S, Kamiya H. Inhibition of the Growth of Calcium-Carbonate Crystals by Multiple Lectins in the Celomic Fluid of the Acorn Barnacle Megabalanus-Rosa. Comp Biochem Phys B. 1994;107:401–9.Google Scholar
- Rhee SG, Kang SW, Chang TS, Jeong W, Kim K. Peroxiredoxin, a novel family of peroxidases. IUBMB Life. 2001;52:35–41.View ArticlePubMedGoogle Scholar
- Robinson MW, Hutchinson AT, Dalton JP, Donnelly S. Peroxiredoxin: a central player in immune modulation. Parasite Immunol. 2010;32:305–13.View ArticlePubMedGoogle Scholar
- Parks WC, Wilson CL, Lopez-Boado YS. Matrix metalloproteinases as modulators of inflammation and innate immunity. Nat Rev Immunol. 2004;4:617–29.View ArticlePubMedGoogle Scholar
- Dickinson GH, Vega IE, Wahl KJ, Orihuela B, Beyley V, Rodriguez EN, et al. Barnacle cement: a polymerization model based on evolutionary concepts. J Exp Biol. 2009;212:3499–510.PubMed CentralView ArticlePubMedGoogle Scholar
- Dreanno C, Kirby RR, Clare AS. Locating the barnacle settlement pheromone: spatial and ontogenetic expression of the settlement-inducing protein complex of Balanus amphitrite. P R Soc B. 2006;273:2721–8.View ArticleGoogle Scholar
- Matsumura K, Nagano M, Fusetani N. Purification of a larval settlement-inducing protein complex (SIPC) of the barnacle, Balanus amphitrite. J Exp Zool. 1998;281:12–20.View ArticleGoogle Scholar
- Clare AS, Matsumura K. Nature and perception of barnacle settlement pheromones. Biofouling. 2000;15:57–71.View ArticlePubMedGoogle Scholar
- Holm ER, Orihuela B, Kavanagh CJ, Rittschof D. Variation among families for characteristics of the adhesive plaque in the barnacle Balanus amphitrite. Biofouling. 2005;21:121–6.View ArticlePubMedGoogle Scholar
- Rittschof D, Orihuela B, Stafslien S, Daniels J, Christianson D, Chisholm B, et al. Barnacle reattachment: a tool for studying barnacle adhesion. Biofouling. 2008;24:1–9.View ArticlePubMedGoogle Scholar
- Aronesty E. EA-utils: Command-line tools for processing biological sequencing data, version 1.1.2-537. [http://code.google.com/p/ea-utils]
- Kong Y. Btrim: A fast, lightweight adapter and quality trimming program for next-generation sequencing technologies. Genomics. 2011;98:152–3.View ArticlePubMedGoogle Scholar
- Li WZ, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9.View ArticlePubMedGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011; 12:doi 10.1186/1471-2105-1112-1323.
- Fears K. Surface-active glasses as regenerative anti-fouling materials. 2014, USPTO# 14/504,657.Google Scholar
- Keller A, Nesvizhskii AI, Kolker E, Aebersold R. Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem. 2002;74:5383–92.View ArticlePubMedGoogle Scholar
- Nesvizhskii AI, Keller A, Kolker E, Aebersold R. A statistical model for identifying proteins by tandem mass spectrometry. Anal Chem. 2003;75:4646–58.View ArticlePubMedGoogle Scholar