Skip to main content

Advertisement

Transcriptional characterisation of the Exaiptasia pallida pedal disc

Article metrics

Abstract

Background

Biological adhesion (bioadhesion), enables organisms to attach to surfaces as well as to a range of other targets. Bioadhesion evolved numerous times independently and is ubiquitous throughout the kingdoms of life. To date, investigations have focussed on various taxa of animals, plants and bacteria, but the fundamental processes underlying bioadhesion and the degree of conservation in different biological systems remain poorly understood. This study had two aims: 1) To characterise tissue-specific gene regulation in the pedal disc of the model cnidarian Exaiptasia pallida, and 2) to elucidate putative genes involved in pedal disc adhesion.

Results

Five hundred and forty-seven genes were differentially expressed in the pedal disc compared to the rest of the animal. Four hundred and twenty-seven genes were significantly upregulated and 120 genes were significantly downregulated. Forty-one condensed gene ontology terms and 19 protein superfamily classifications were enriched in the pedal disc. Eight condensed gene ontology terms and 11 protein superfamily classifications were depleted. Enriched superfamilies were consistent with classifications identified previously as important for the bioadhesion of unrelated marine invertebrates. A host of genes involved in regulation of extracellular matrix generation and degradation were identified, as well as others related to development and immunity. Ab initio prediction identified 173 upregulated genes that putatively code for extracellularly secreted proteins.

Conclusion

The analytical workflow facilitated identification of genes putatively involved in adhesion, immunity, defence and development of the E. pallida pedal disc. When defence, immunity and development-related genes were identified, those remaining corresponded most closely to formation of the extracellular matrix (ECM), implicating ECM in the adhesion of anemones to surfaces. This study therefore provides a valuable high-throughput resource for the bioadhesion community and lays a foundation for further targeted research to elucidate bioadhesion in the Cnidaria.

Background

Adhesion is ubiquitous in nature, having evolved in all kingdoms ranging in length-scale from macromolecule to whole organisms, and for applications as diverse as cell-cell adhesion, biofilm formation, locomotion, surface attachment, prey-capture, reproduction and the formation of protective casings [9, 31, 35, 86]. The ability of some organisms to attach themselves to surfaces, either permanently or temporarily, has received particular attention. Marine organisms have been highly represented in such studies due to their impressive ability to circumvent environmental conditions that challenge synthetic adhesives [78]. Having evolved independently on numerous occasions [9, 44, 67], bioadhesion has provided organisms from single-cell prokaryotes to complex vertebrates with the adaptability to specialise and survive in niche environments.

Most organisms adhering to surfaces do so via a complex interplay of morphology, behaviour and chemistry. A liquid secretion is often involved, the composition of which can be highly variable and may or may not harden post-secretion. Adhesion to surfaces has been studied, to varying degrees, in sponges [41], freshwater hydrozoans (e.g. Hydra spp.; [80, 82]), polychaete worms (e.g. Phragmatopoma californica; [48]), gastropod molluscs (e.g. Arion spp.; [87]), bivalve molluscs (e.g. Mytilus spp.; [85, 99]), barnacles (e.g. Pollicipes & Balanus spp.; [37, 79, 88]), caddisfly larvae (Stenopsychie marmorata; [91]), echinoderms [32, 33], and tunicates [1, 108] to name several. The diversity of metazoan taxa in this restricted list demonstrates not only the importance of bioadhesion for survival in aquatic environments, but also the interest of the research community in identifying convergent themes that run through the adhesives of unrelated species, and which may provide inspiration for future biotechnological applications.

To identify and characterise ancestral mechanisms of adhesion it is necessary to examine adhesion in basal metazoan species. The phylum Cnidaria includes corals, anemones and jellyfish, among others, and is a sister group of the bilateria that diverged > 500 million years ago [90]. Unlike the Porifera (sponges), cnidarians possess body-axis symmetry, a nervous system [51] and other characteristics of complex eukaryotes. To date, studies of cnidarian species have increased our understanding of animal development [27], neural networks [101], immunity [8] and symbiosis [47]. Although the field of bioadhesion is in its infancy, much can therefore be learned from established cnidarian models (e.g. Exaiptasia pallida.; Hydra magnipapillita) that have toolkits of molecular techniques available for them.

A recent study of the freshwater hydrozoan H. magnipapillita [80] used a combination of transcriptomics, proteomics and in-situ hybridisation to provide foundational knowledge of the molecular regulation of bioadhesion in that species. However, general understanding of cnidarian bioadhesion is lacking and the evolutionary distance between marine cnidarians, such as anemones, and freshwater hydrozoans invites more detailed investigation of marine species. Surprisingly, literature on the adhesion of anemones to surfaces is almost completely absent. Early observations of the swimming anemone, Stomphia coccinea, suggested that nematocysts were involved in surface adhesion, whilst locomotion and detachment were controlled by muscle contraction [29]. In Actinia equina, it was later suggested that bioadhesion could rely on protein-protein interactions rather than mucopolysaccharides or nematocysts [107]. Only one, small-scale, study exists on the screening and identification of adhesion proteins in the cadherin-catenin complex (CCC) associated with cell adhesion in Nematostella vectensis [21]. There is, however, no evidence that this moiety relates directly to surface adhesion and N. vectensis is, in any case, a sediment-dweller that does not adhere to solid substrates.

The glass anemone, Exaiptasia pallida [38], is a fast-growing, symbiotic species native to shallow waters of the western Atlantic, Caribbean Sea and Gulf of Mexico. Growing to approximately three centimetres when extended, the body plan consists of a pedal disc attached to the substratum, a slender peduncle and an oral disc (with mouth) surrounded by up to 96 tentacles (Fig. 1). Two rows of slits on the peduncle contain nematocyst-armed acontia that extend for defence. Interest in this anemone has grown in recent years for three reasons: 1) Its taxonomic position within the Anthozoa places it close to hard corals and, like many corals, it exhibits a symbiotic relationship with photosynthetic zooxanthellae. It has thus been developed into a model for studying coral-reef, climate-change dynamics; 2) E. pallida is easy to culture in the laboratory, making it ideal for manipulation experiments; 3) Publication of the E. pallida (strain CC7) genome [13] has provided a suite of DNA sequences, expression-supported mRNAs and predicted proteins for on-going [47, 55, 68] and future high-throughput-omics research.

Fig. 1
figure1

The anatomical body plan of Exiptasia pallida. The pedal disc is located at the aboral end of the oral-aboral axis. Image taken by authors of this study

mRNA-sequencing (mRNA-seq) is a powerful technology for understanding how organisms function, develop and respond to their environments [23], allowing for the detection of novel genes, metabolic pathways, diseases and bioactive compounds [73, 84, 100]. mRNA-seq has been used previously in the field of bioadhesion (e.g. [42, 54]), allowing for interpretation of the genes, molecular processes and regulatory networks associated with adhesion in a number of species [19, 24, 58, 63, 80]. To our knowledge, however, this approach has not been applied to understanding the transcriptional processes involved in adhesion of the anemone pedal disc to surfaces. In this study, a genome reference-assisted transcriptomics approach was taken to identify the regulatory processes associated with the pedal disc of E. pallida. The findings highlight the importance of the basal tissue not only for adhesion to surfaces, but also for development, defence and immunity.

Results

Transcriptome assembly and differential expression analysis

After discarding low quality base pairs, rRNA contamination and ambiguous reads, 300,294,692 paired-end reads were uniquely mapped to the E. pallida genome [13]. For individual libraries, the percentage of reads that uniquely mapped to the genome ranged from 78.70 to 85.21% (Table 1). Twenty-four thousand, eight hundred and forty-four genes (84.88%) from the genome were found to be expressed in the transcriptome (supported by one or more counts in a library). Of these, 547 (2.20%) genes were differentially expressed, 427 (1.72%) were significantly upregulated and 120 (0.48%) were significantly downregulated in the pedal disc. A principle component analysis (PCA) represents the within and between group similarity of normalised mRNA libraries, with principle components 1 and 2 representing 94% of variability (Fig. 2).

Table 1 Number of reads used for the alignment of RNA libraries to the E. pallida genome after rRNA removal. Number and percentage of reads uniquely mapped to the E. pallida genome. Ambiguous multi-mapped reads were discarded. WA = whole animal, AM = amputated animal (without pedal disc)
Fig. 2
figure2

PCA plot showing 94% of the variance between and within groups of mRNA libraries. WA = Whole animal (with pedal disc); AM = Amputated animal (without pedal disc). One of the AM biological replicates was an outlier and as a result it was removed from the experiment (not shown)

Gene ontologies and protein superfamilies of the pedal disc

A total of 141 GO terms were enriched in the pedal disc (FDR: 0.05; ≥ 5 genes represented by GO term; Additional file 1): 93 biological processes (BP), 24 molecular functions (MF) and 24 cellular components (CC). Collapsing of terms with ReviGO resulted in a total of 41 terms (20 BP; 12 MF; 9 CC; Fig. 3). Of the enriched BP terms (Fig. 3a): modification of morphology or physiology of other organism (GO:0035821), cytolysis (GO:0019835) and cell killing (GO:0001906) were among the most enriched. Further queries of the genes representative of these terms identified genes encoding Matrix metalloproteinase-9, Leukotoxin, toxins (types AvTX, CaTX, CfTX, CrTX, PsTX) and Venom prothrombin activator oscutarin-C non-catalytic subunits. The TX (Toxin) genes identified in this study received a functional annotation to Actineria villosa (AvTX), Carybdea alata (CaTX), Chironex fleckeri (CfTX), Carybdea rastonii (CrTX) and Phyllodiscus semoni (PsTX) homologs. Terms including: regeneration (GO:0031099), positive regulation of developmental process (GO:0051094), developmental process (GO:0032502), response to steroid hormone (GO:0048545), proteolysis (GO:0006508), cellular modified amino acid biosynthetic process (GO:0042398), biological adhesion (GO:0022610) and cell adhesion (GO:0007155) were also enriched.

Fig. 3
figure3

Enriched (FDR: 0.05; ≥ 5 genes per term) gene ontology (GO) terms identified in the pedal disc. Terms are representative of collapsed ReviGO semantic SimRel terms. a Biological Processes; b Molecular Functions; c Cellular Components. Significant enrichment P-values (Log10) are represented by the X-axis. GO Term Log size (representative gene grouping size) is represented by the Y-axis

Among the enriched MF terms (Fig. 3b): Calcium ion binding (GO:0005509), peptidase activity, acting on L-amino acid peptides (GO:0070011), serine hydrolase activity (GO:0017171) and iron ion binding were identified (GO:0005506). Of the nine CC terms identified (Fig. 3c), the terms, extracellular region (GO:0005576) and extracellular region part (GO:0044421) were most enriched. Vesicle (GO:0031982), collagen trimer (GO:0005581), endoplasmic reticulum lumen (GO:0005788), membrane-enclosed lumen (GO:0031974), extracellular membrane-bounded organelle (GO:0065010), other organism cell membrane (GO:0044218) and nematocyst (GO:0042151) were also identified.

Eight condensed GO terms were identified as depleted in the pedal disc (FDR: 0.05; ≥ 5 genes represented by GO term; Additional file 1). These included three BP terms: Protein autoubiquitination (GO:0051865), extracellular matrix organisation (GO:0030198) and extracellular structure organisation (GO:0043062). Two MF terms were depleted: Calcium ion binding (GO:0005509) and extracellular matrix structural constituent (GO:0005201), as well as three CC terms: Extracellular region (GO:0005576), axon (GO:0030424) and extracellular matrix component (GO:0044420). Extracellular matrix structural constituent (GO:0005201) was represented by one Fibrillin-1, two Fibrillin-2, a Mucin-4, a Papillin and four Collagen alpha-1 related genes.

Nineteen protein superfamily domain classifications (Fig. 4; Additional file 1) were enriched in the pedal disc (FDR: 0.05; ≥ 5 genes represented by superfamily classification). Top superfamilies included: EGF/Laminin (SCOP: 57196), beta-Roll (SCOP: 51120), Thrombospondin-1 (TSP-1) type 1 repeats (SCOP: 82895), Invertebrate chitin binding proteins (SCOP: 57625) and PR-1-like (SCOP: 55797). Other superfamilies were associated with oxidoreductase, binding and enzymatic functions including Trypsin-like serine proteases (SCOP: 50494; Fig. 4). Eleven superfamilies were depleted in the downregulated gene set (FDR: 0.05; ≥ 5 genes represented by superfamily classification; Additional file 1). Superfamilies depleted in the pedal disc included: Spectrin repeat (SCOP: 46966), Thyroglobulin type-1 domain (SCOP: 57610), Cadherin-like repeats (SCOP: 49313), and Immunoglobulin (SCOP: 48726).

Fig. 4
figure4

Superfamily protein domain classifications significantly enriched (FDR: 0.05; ≥ 5 genes per classification) in the E. pallida pedal disc

Upregulated genes of the pedal disc

There were 427 significantly upregulated genes in the pedal disc, expression differences ranged from 2.00 to 88.65 fold greater (Log2 Fold Change (FC) 1.00–6.47). Of those, 250 genes (58.55%) had meaningful functional annotations (e.g not putative, hypothetical, predicted or uncharacterised; [13]). Furthermore, 173 genes (40.52%) were predicted to encode for extracellularly secreted proteins (ESPs). 145 (33.96%) upregulated genes were predicted as containing a signal peptide (classically secreted) and 28 (6.56%) genes were predicted to be secreted without a signal peptide (non-classically secreted). Table 2 and Fig. 5 depict the top 35 genes that received a functional annotation. Eighty-five (49.13%) of the 173 genes had no meaningful functional annotation (as defined above). For a full repertoire of all 173 genes (including those without meaningful functional annotations), see Additional file 1.

Table 2 Top 35 ab initio predicted extracellular secreted proteins which received a meaningful functional annotation. Mode = Secretory mode: C = Classical secretory pathway (Predicted by SignalP5.0); NC = Non-classical secretory pathway (Predicted by SecretomeP2.0). FC = Log2 fold change of gene expression
Fig. 5
figure5

Heatmap depicting the top 35 functionally annotated extracellularly secreted protein (ESP) encoding genes upregulated (Orange) in the pedal disc of E. pallida. Z-scores are representative of centred gene expression values (Log2). WA_1 – WA_4 = Whole Animal; AM_1 – AM_3 = Amputated animals (pedal disc amputated). * = Non-classically secreted; ** = Classically secreted

The most significantly upregulated gene in the pedal disc (FC 6.47; Aipgene2358; Deleted in malignant brain tumors 1 protein; DMBT1), encoded an ESP containing a single spermadhesion-type CUB domain and two scavenger receptor cysteine-rich (SRCR) domains (Fig. 6). Aipgene7095 (FC 3.04) and Aipgene595 (FC 1.55; Fig. 6) also encoded DMBT1 homologs.

Fig. 6
figure6

Protein domain architecture of select predicted extracellularly secreted proteins. Length of protein is highlighted as AA (Amino Acids). CUB = Complement C1r/C1s, Uegf, Bmp1; EGF/LAM = Epidermal Growth Factor/ Laminin; SRCR-like = Scavenger receptor cysteine-rich-like; ICBP = Invertebrate chitin binding protein domain; KSPI = Kazal-type Serine Protease Inhibitor; TSP-1 type 1 = Thrombospondin-1 type 1; YWTD = YWTD AA repeat; Colour of domain classifications are consistent with labelling. Aipgene8196 contained overlapping domains as shown

Fifteen upregulated genes had collagen-related functional annotations. Seven of these encoded for ESPs, including: Aipgene20335 (Collagen alpha-1(XXVII) chain; FC 2.13), Aipgene20425 (Collagen alpha-1(XII) chain; FC 1.77), Aipgene14299 (Collagen alpha-1(XII) chain; FC 1.61), Aipgene17088 (Short-chain collagen C4 (Fragment); FC 2.24), Aipgene24564 (Short-chain collagen C4 (Fragment); FC 3.68); Aipgene5529 and Aipgene13236 (Collagen triple helix repeat-containing protein 1, CTHRC1; FC 4.20; FC 2.22). Aipgene 20,425 and Aipgene14299 contained von Willebrand Type A (VWA) domains.

Nine metalloproteinase genes were identified encoding for ESPs. These included: Aipgene11959 (Matrix Metalloproteinase-24; FC 3.51), Aipgene16517 (Disintegrin and metalloproteinase domain-containing protein 10; FC 2.13), Aipgene9574 (Zinc metalloproteinase nas-13; FC 2.07), Aipgene1297 (Matrix metalloproteinase-9; FC 1.62), Aipgene25980 (MAM and LDL-receptor class A domain-containing protein 1; FC 2.01), Aipgene23879 and Aipgene23885 (Nematocyte expressed protein 6 (NEP-6); FC 1.71, 3.36), Aipgene28234 and Aipgene28816 (Blastula Protease 10; FC 3.35, 2.31).

Extracellularly secreted glycoproteins and proteoglycans were found to be encoded by upregulated genes. These included: Aipgene8196 (Uromodulin; FC 1.59; Fig. 6), Aipgene25531 (Fibronectin; FC 2.13), Aipgene23599 and Aipgene23611 (Rhamnose-binding lectin; FC 1.76, 1.37), Aipgene16897 (Fibropellin-1; FC 1.74; Fig. 6), Aipgene14232 (Neurogenic locus notch homolog protein 3; FC 1.93), Aipgene27786 (Agrin; FC 2.50; Fig. 6), Aipgene4050 (Cubilin; FC 2.75) and Aipgene10560 (Chondroitin Proteoglycan 2, FC 2.17; Fig. 6). Aipgene16897 (Fibropellin-1; Fig. 6) contained 22 EGF (Epidermal Growth Factor) repeats, while Aipgene14232 (Neurogenic locus notch protein) contained 7. Although not predicted as being extracellularly secreted in this study, Aipgene16924 (Neurogenic locus notch protein) was differentially expressed and contained 19 EGF repeats. Aipgene10560 (Chondroitin Proteoglycan 2) contained 11 Invertebrate chitin-binding (peritriphin-A) domains. Additionally, five genes encoding extracellularly secreted Semaphorin-5B (one gene; FC 2.68) and Hemicentin-1 (four genes) contained more than one TSP-1 type 1 repeat (Fig. 6; Additional file 1).

In conjunction with the enrichment of GO terms associated with proteases, a host of trypsin-like serine protease genes were identified as encoding ESPs. These included Aipgene7203 (Coagulation factor XI, FC 4.67), Aipgene5888 (Chymotrypsinogen B, FC 4.59), Aipgene7224 (Plasma kallikrein, FC 2.23), Aipgene7188 and Aipgene7180 (Transmembrane protease serine 3; FC 1.91, 1.85), Aipgene7205 (Transmembane protease serine 5; FC 1.88) and Aipgene9301 (chymotrypsin-like protease CTRL-1; FC 1.53). On the contrary, the ESP encoded by Aipgene27786 (Agrin; FC, 2.50) contained 10 kazal-type serine protease inhibitors (KSPI). Aipgene8196 (Uromodulin; FC 1.59) also contained one KSPI.

Discussion

Transcriptome assembly

In this study a transcriptome was assembled to enable profiling and characterisation of genes expressed in the pedal disc of E. pallida. The high percentage of reads (78.70 to 85.21%) aligning to the genome after quality checks and screening provided ample sequencing depth for a comparative transcriptomic analysis [66]. Analysis of the transcriptional processes up-regulated in the pedal disc provided significant insight into the processes of extracellular matrix (ECM) secretion, adhesion, defence, immunity, and the development of the tissue in contact with the substratum.

Gene ontology enrichment

The enriched biological process (BP) terms indicated a prominent role for immunity and pathogen defence in the pedal disc. Further queries of the genes associated with toxin production suggested conservation of toxin-encoding sequences throughout species of the Anthozoan and Cubozoan lineages. The majority of toxin genes identified are involved with the formation of the membrane-attack complex/perforin in anemone species [71, 83]. Although prothrombin activators are found in venomous animals, these genes promote thrombin synthesis required for coagulation, cellular aggregation and inflammatory responses in eukaryotes [12]. These genes may play a dual role in defence and cellular aggregation within the pedal disc. Increased defence and immunity have often been associated with adhesive tissues in order to prevent disease and degradation [26, 52]. Until now, most immune studies in the Cnidaria have focussed on H. magnipapillata, Aurelia aurita [18] and corals [74]. These findings are therefore valuable.

Development, regeneration, cellular homeostasis and degradation processes were prominent in the pedal disc. Enriched terms in the pedal disc indicated that cellular turnover rates may be higher than in the rest of the animal, an observation that aligns well with the asexual fission performed at the pedal disc of E. pallida. Maintenance of the ECM throughout pedal disc adhesion would also require constant remodelling and protein turnover [81]. This is especially important given that E. pallida is somewhat motile, with the ability to reverse its adhesion to a surface facilitating locomotion. The enrichment of terms, biological adhesion (GO:0022610) and cell adhesion (GO:0007155) imply strongly that the ECM is involved in pedal disc adhesion. Furthermore, depletion of ECM-associated terms in the down-regulated genes imply significant changes and ECM remodelling through the regulation of gene expression.

Given that concentrations of calcium in the extracellular space are four to five orders of magnitude greater (typically 1.2 mM) than intercellular concentrations [50], enrichment of calcium ion binding (GO:0005509) is logical, especially when many ECM protein domains, including EGF/laminins, TSP-1, and C-type lectins bind calcium in order to structurally stabilise the ECM [43]. Enrichment of terms associated with peptidase activity, iron ion binding and peptidase regulation in conjunction with enrichment of KSPIs (SCOP: 100895) and Trypsin-like serine proteases (SCOP: 50494) highlight the importance of serine proteases and protease inhibitors in the pedal disc. Indeed, serine-rich proteins are often found within the adhesives of marine invertebrates [30, 32, 37, 99]. Serine-rich proteins confer adhesive and cohesive ability in marine organisms [109]. Furthermore, use of a serine protease resulted in degradation of the cyprid adhesive of Balanus amphitrite [2]. Hypothetically, serine proteases would be required within an ECM based attachment/detachment system to disrupt cross-linking of phosphoserine residues [42] and, as discussed, have a recognised role in degrading and remodelling the ECM [39].

With respect to enriched CC GO terms, the enrichment of collagen trimer (GO:0005581) is congruent with the literature in which collagen, secreted by fibroblast cells, is one of the major constituents of the ECM [89]. The enriched CC GO terms also imply that the secretory pathway is enhanced in the pedal disc. In theory, for significant extracellular secretion of a proteinaceous matrix, an enriched (larger) endoplasmic reticulum lumen would be required (where post-translational modifications and protein folding occur; [62]), along with increased numbers of vesicles ([97]; Fig. 7). Results of this study are therefore in agreement with those of Young et al.[107] who suggested that adhesion is strongly associated with protein-protein interactions. Overall in comparison to other anemone studies, common GO groups are represented across many different tissue types and not just one type [11, 61, 93]. For example, cell adhesion (GO:0007155) and calcium ion binding (GO:0005509) are represented in the nematosomes, mesentries and tentacles of Nematostella vectensis [11]. Statistical enrichment does; however, suggest GO groupings identified in this study are important for the functioning of the pedal disc.

Fig. 7
figure7

Schematic illustration depicting, a Cellular machinery involved in extracellular matrix protein synthesis and trafficking of extracellularly secreted proteins to the extracellular space. b ECM protein components suspected to form the basis of pedal disc adhesion in E. pallida. Image taken and illustrations created by authors of this study

Protein superfamily enrichment

For metazoans to evolve from single-celled organisms, cell adhesion was obligatory. New protein domains including EGF, TSP-1, C-type lectins, collagen triple-helix domains, laminins, cadherins and integrins evolved to facilitate adhesion [46]. Enrichment of EGF domains in the pedal disc, with high numbers of this domain appearing in extracellularly secreted Aipgene16897 (Fibropellin-1, Fig. 6) and Aipgene14232 (Neurogenic locus notch homolog protein3) proteins is in line with results of other molecular adhesion studies. In the flatworm, Macrostomum lignano, the 17 EGF domain-containing protein, Mlig-ap1, exhibits a cohesive function in the adhesive [105]. EGF repeats have also been identified in echinoderm and mussel adhesion [45]. Whilst the large glycoprotein, Fibropellin, forms an ECM layer known as the apical lamina in the sea urchin, Strongylocentrotus purpuratus [16] and sea cucumber Apostichopus japonicus [10]. Neurogenic locus notch protein is a receptor involved in cell-cell interaction, cell fate and differentiation [49].

Eleven Invertebrate chitin-binding (peritrophin-A) domains were present within Aipgene10560 (Chondroitin proteoglycan 2, FC 2.17, Fig. 6). In H. magnipapillata, chitin-binding (peritrophin-A) domains were suggested to play a role in pedal disc adhesion [80]. Proteoglycans are a major constituent of the ECM and basal membrane [75, 94], consisting of a core peptide with heavily glycosylated sidechains. Chitin-binding proteins have been identified in the gastrolith ECM of crayfish where they cross-link to harden the ECM via the oxidation of phenols or catechols [36].

TSP-1 type 1-containing semaphorin glycoproteins are known guidance molecules in axon development, associated with the ECM [4]. Hemicentin-1 (Fig. 6) is a large extracellular glycoprotein of the immunoglobulin family, conserved in the eukaryotic lineage [98]. In H. magnipapillata, N. vectensis [96], oysters [34] and echinoderm adhesomes, hemicentin homologs have also been identified [102]. This protein supports architectural and structural integrity of animal tissues, including the ECM [106]. Enriched PR-1 (Pathogenesis related 1)-like domains appeared in 13 genes (Golgi-associated plant pathogenesis related protein 1; GAPR1) of this study. The majority of these proteins are located in the Golgi apparatus, however some are known to be secreted extracellularly [70]. Eight were predicted to be extracellularly secreted in this study. These proteins contain cysteine-rich regions associated with innate immunity and allergenic effects. In the Cnidaria they are commonly associated with nematocysts [72]. Closer examination of depleted superfamilies revealed that genes representing some of these classifications were associated with cytoskeletal development.

Upregulation of specific genes in the pedal disc

Although 58.55% of the upregulated genes in this study had a functional annotation, the remaining proportion had no functional annotation or a meaningful annotation. The lack of annotation calls for further functional characterisation of the E. pallida genome. Generating longer reads via a sequencing platform such as PacBio or Oxford Nanopore and improvement of annotation databases has potential to improve the completeness and accuracy of the genome assembly [13]. Like all reference genomes, the E. pallida genome assembly and databases will most likely continue to evolve over time [40].

The most differentially expressed gene in the pedal disc, Aipgene2358 (DMBT1), contained one spermadhesion type CUB domain and two SRCR protein domains (Fig. 6). Another DMBT1 homolog identified in this study, Aipgene595, contained additional domains (Fig. 6). The spermadhesion, CUB domain possesses the ability to bind to ligands e.g. carbohydrates, glycosaminoglycans, phospholipids and protease inhibitors [95]. CUB domains are typically found in proteins involved in developmental processes [17]. SRCR domains, consisting of serine-threonine-rich amino acid motifs, are conserved across the metazoan lineage. DMBT1 has been found to encode for three types of glycoprotein: Deleted in malignant brain tumours 1 protein, salivary agglutinin (DMBTSAG) and lung glycoprotein-340 (DMBT1GP340; [57]). This protein is a versatile mucin-like molecule that has involvement in epithelial differentiation, agglutination and defence [64]. Given the absence of organs in Cnidaria, it is likely that these glycoproteins serve a role in innate immunity and agglutination of other proteins in the ECM of the pedal disc. DMBT1 was found to be upregulated in the corals Acropora millepora and Orbicella faveolata in response to bacterial challenge [104].

Seven genes were predicted to code for extracellularly secreted proteins associated with collagen metabolism: Aipgene20335 (Collagen alpha-1(XXVII) chain), Aipgene20425 and Aipgene14299 (Collagen alpha-1(XII) chain), Aipgene17088 and Aipgene24564 (short-chain collagen C4 (Fragment)), Aipgene5529 and Aipgene13236 (CTHRC1). Collagen alpha-1 is the main constituent of type 1 collagen. Short-chain collagen C4 (a component of Type IV collagen) is associated with forming sheet-like structures underlying the basal membranes of epithelial and endothelial tissues, surrounding muscle cells, peripheral nerves and adipocytes [7]. The two genes encoding collagen alpha-1(XII) chains contained von Willebrand Type A domains involved in adhesion [32, 87, 103]. In contrast, the glycoprotein Aipgene5529 (CTHRC1) inhibits collagen synthesis in ECMs and could be involved in a regulatory feedback loop controlling collagen deposition in the pedal disc. This protein is also required for epithelial-mesenchymal transition and cellular migration [59].

Matrix metalloproteinases (MMP) are implicated in the regulation of growth factors and their receptors, cytokines, chemokines, adhesion receptors and cell surface proteoglycans in order to alter cellular responses to the environment [89]. MMPs are required for normal developmental processes, regeneration of tissue and degradation of ECM proteins [6]. Blastula Protease 10 is a member of the astacin family of zinc-dependant endopeptidases and in sea urchins this MMP is involved in development, contains a tyrosine switch and is influenced by calcium binding [22]. NEP-6 metalloproteinases are associated with nematocysts and are thought to have a dual-role, acting as a potassium-channel toxin and degrading ECM proteins [61]. In conjunction with the enrichment of nematocyst (GO:0042151), the identification of NEP-6 genes and GAPR1 genes, an atypical gland-like nematocyst may play a role in pedal disc adhesion or defence [65].

Uromodulin (Tamm-Horsefall protein) was classified as an upregulated ESP. This protein undergoes heavy glycosylation and promotes protein-protein interactions, resulting in gelification [77]. It contained two EGF/laminin and a C-lectin domain (dependant on calcium for binding; Fig. 6). In agreement with the current study, uromodulin mRNA was found to be localised in the distal portion of the aboral pore of N. vectensis, where it may play a role in defence. If tissue damage occurs, mRNA accumulates at the wound site and is considered to be involved in regeneration [28]. The upregulated ESP, fibronectin, is another high-molecular weight glycoprotein that binds large numbers of cell adhesion receptors including proteoglycans, growth factors, collagen, integrins, and other ECM proteins to strengthen the structure of the ECM [110]. Additionally, active fibronectin fibrillogenesis is a prerequisite for the deposition of collagen type I to the ECM [53]. Cubilin, also an upregulated ESP in this study, is a receptor found on the cell surface binding galectin-3, a lectin that promotes cell-matrix interactions [69]. Two rhamnose-binding lectins were also classified as upregulated ESPs in this study. In H. magnipapillata, six copies of rhamnose-binding lectins were found to be potentially associated with adhesion [80]. This number may of course be exaggerated by the use of de novo transcriptomics in that study, where difficulties can be encountered defining unigenes and splice variants. As discussed by Rodrigues et al. [80], glycan cross-binding has the potential to facilitate non-covalent cross-linking, increasing cohesion and adhesion strength [54]. Agrin was also an upregulated ESP in this study (Fig. 6). It is a large heparan sulphate proteoglycan known to bind laminin and integrins in the ECM, and is involved in postsynaptic clustering [102]. A schematic diagram details some protein classes believed, on the basis of these data, to be involved in the ECM and possibly adhesion of E. pallida to surfaces (Fig. 7).

Conclusion

The findings of this study, in conjunction with past observations, suggest that adhesion of the E. pallida pedal disc may be facilitated through the secretion of an ECM-like proteinaceous matrix, containing collagen, glycoproteins, proteoglycans and lectins. Conversely, the metalloproteinases and serine proteinases identified here may play roles in immunity, degradation and remodelling of the ECM to facilitate and prime the detachment of E. pallida from surfaces. Cross-linking via oxidoreductase reactions may also occur to strengthen adhesion. The methods used resulted in a list of high-confidence ab initio predicted extracellularly secreted proteins. Functional characterisation of proteins, morphological analyses and localisation of mRNA are now necessary to validate the predicted role of these elements. This study and the datasets supplied thus provide a foundation for future research investigating the regenerative capability, fission, immunity and, in particular, the bioadhesion capability of E. pallida.

Methods

Culturing and sampling

Symbiotic Exaiptasia pallida (Strain CC7) were cultured and maintained in polycarbonate 5 L containers within a temperature (26 °C) and light-controlled (~ 60 μmol m− 2 s− 1 12 h light: 12 h dark) incubator at the School of Natural and Environmental Sciences, Newcastle University. Salinity was maintained at 35 parts per thousand (ppt) using artificial seawater (TropicMarin™). Anemones were fed with stage 1 Artemia sp. nauplii three times per week. Artificial seawater was replenished after feeding. Anemones were starved 48 h prior to sampling in order to minimise contaminating molecular artefacts. For RNA samples, three anemones were pooled to form one biological replicate. In total, four biological replicates were used in this study per tissue type. Two tissue types were sampled and snap frozen in liquid nitrogen: Whole animal (WA) and amputated animal (AM), the latter consisting of the entire animal without the pedal disc. Pedal discs were surgically removed using a sterile razor. Amputated animals (minus the pedal discs) were used instead of pedal disc tissue. This method was previously used on H. magnipapillata [80] to minimise tissue damage.

Library preparation, sequencing, quality checks and assembly

Total RNA was obtained from RNA samples by homogenisation with TRIzol Reagent® and the use of a Direct-zol™ RNA Miniprep Plus kit (per manufacturer’s guidelines; Zymo Research). In-house quality checks were performed at the Leeds University NGS facility to ensure all samples achieved a RIN (RNA Integrity number) value of 7 or above for high-quality sequencing. mRNA was enriched using an Illumina TruSeq kit and sequenced with an Illumina Nextseq 500 sequencer (Paired-end: 76 bp × 2). Quality of reads was visualised using FastQC version 0.11.8 software [5]. Reads were quality trimmed and any adaptor sequences were removed with BBDuk of the BBTools software package [20]. The following parameters were used in conjunction with the BBDuk adapters.fa reference file: ordered = t; ktrim = r; k = 23; mink = 11; hdist = 1; qtrim = rl; trimq = 10; minlength = 35; tpe; tbo. Although mRNA enrichment was completed prior to sequencing (Illumina TruSeq Kit), a second rRNA screening and removal step was conducted using BBDuk (parameters: k = 31; hdist = 1) in conjunction with the associated custom ribokmers.fa file [20]. The decontaminated reads were aligned to the E. pallida genome (Version 1.0; [13]; http://aiptasia.reefgenomics.org/download/) using STAR ultrafast aligner software (Version 2.7; [25]). A genome index was created with the settings: sjdbGTFtagExonParentTranscript Parent and sjdbOverhang = 75. The largest intronic region of the genome was calculated in conjunction with the mRNA .fasta and .gff3 files of the version 1.0 genome and, as a result, the option alignIntronMax = 70,000 was used as an additional alignment parameter. Subread Featurecounts software (version 1.6.4; [56]) was used to obtain a summary of gene counts for all samples against the genome. The count matrix was imported into R Studio (Version 3.5.3; “Great Truth”). Reads with no evidential support were discarded (0 counts in all libraries). DESeq2 software (Version 1.22.2; [60]) was used to perform normalisation and differential gene expression analysis. Genes were identified as differentially expressed according to the following criteria: P-value = 0.05; alpha = 0.05; LogThreshold = 1(log2 based value). Additional file 2 depicts the bioinformatics workflow utilised in this study. Raw sequencing reads were deposited under the NCBI Sequence Read Archive (SRA) accession number: PRJNA540572.

Gene ontology and protein superfamily enrichment analysis

The sub-sets of significantly upregulated and downregulated genes were defined as the foreground datasets and the entire expressed transcriptome of this study was defined as the background dataset for enrichment analyses. Gene ontology and superfamily protein domain classifications were obtained from the Reef Genomics repository and the authors of the genome publication [13]. Fisher’s exact statistical test was conducted within R using the fisher.test function with the option, alternative = “greater”. Multiple comparison correction was performed with an FDR threshold of 0.05 [15]. Additionally GO terms or protein superfamily domain classifications not represented by at least 5 or more differentially expressed genes were discarded to improve confidence. Those enriched GO terms that met the criterion were condensed and visualised using the Revigo GO visualisation web portal (http://revigo.irb.hr/revigo.jsp; [92]). A ‘small’ (0.5) threshold was chosen for SimRel semantic similarity measure, all other parameters were kept as default. Resulting .csv tables were exported and used for graphically plotting the condensed GO terms.

Identification of genes encoding for putative extracellularly secreted proteins

Classically secreted proteins (CSPs) contain a signal peptide domain located in the N-terminus of the protein, whereas non-classically secreted proteins (NCSPs) do not contain a signal peptide. Rather, they are transported to their end-destination by carrier proteins and lipids [14]. In order to form a list of high-confidence ab initio extracellularly secreted proteins encoded by the significantly upregulated genes, the predicted protein sequences [13] of the upregulated gene set were first scanned with SignalP5.0 [76] with ‘Eukarya’ selected as organism group. Protein sequences which did not possess a signal peptide according to SignalP5.0 were scanned using SecretomeP 2.0 for NCSPs with an applied neural network (NN) threshold of 0.6. Resulting data of identified CSPs and NCSPs was concatenated and scanned using DeepLoc-1.0 which uses deep-learning and neural networks to classify the final location of proteins [3]. ‘Profiles matrix’ was chosen instead of a Blossum62 matrix for searches to increase accuracy. Proteins classified as ‘extracellular’ by DeepLoc-1.0 provided a final set of genes which putatively encoded for extracellular secreted proteins.

Availability of data and materials

Availability and sources of data used in this study are detailed throughout the article.

Abbreviations

AA:

Amino acids

AM:

Amputated animal

AvTX:

Actineria villosa toxin

BP:

Biological processes

CaTX:

Carybdea alata toxin

CC:

Cellular components

CCC:

Cadherin-catenin complex

CfTX:

Chironex fleckeri toxin

CrTx:

Carybdea rastonii toxin

CSPs:

Classically secreted proteins

CTHRC1:

Collagen triple helix repeat-containing protein 1

CUB:

Complement C1r/C1s, Uegf, Bmp1

DMBT1:

Deleted in malignant brain tumors 1 protein

ECM:

Extracellular matrix

EGF:

Epidermal growth factor

ESPs:

Extracellularly secreted proteins

FC:

Fold Change

FDR:

False discovery rate

GAPR1:

Golgi-associated plant pathogenesis related protein 1

GO:

Gene ontology

KSPI:

Kazal-type serine protease inhibitor(s)

MF:

Molecular functions

MMP:

Matrix metalloproteinases

mRNA:

messenger RNA

mRNA-seq:

mRNA-sequencing

NCSPs:

Non-classically secreted proteins

NEP-6:

Nematocyte expressed protein 6

NN:

Neural network

PCA:

Principle component analysis

ppt:

Parts per thousand

PR-1:

Pathogenesis related 1

PsTX:

Phyllodiscus semoni toxin

RIN:

RNA integrity number

rRNA:

ribsomonal RNA

SCOP:

Structural classification of proteins

SRA:

Sequence Read Archive

SRCR:

Scavenger receptor cysteine-rich

TSP-1:

Thrombospondin-1

TX:

Toxin

VWA:

von Willebrand type A

WA:

Whole animal

References

  1. 1.

    Aldred N, Clare AS. Mini-review: impact and dynamics of surface fouling by solitary and compound ascidians. Biofouling. 2014;30(3):259–70.

  2. 2.

    Aldred N, Phang IY, Conlan SL, Clare AS, Vancso GJ. The effects of a serine protease, Alcalase®, on the adhesives of barnacle cyprids ( Balanus amphitrite ). Biofouling. 2008;24(2):97–107.

  3. 3.

    Almagro Armenteros JJ, Sønderby CK, Sønderby SK, Nielsen H, Winther O. DeepLoc: prediction of protein subcellular localization using deep learning. Hancock J, editor. Bioinformatics. 2017;33(21):3387–95.

  4. 4.

    Alto LT, Terman JR. Semaphorins and their Signaling Mechanisms. New York: Humana Press; 2017. p. 1–25.

  5. 5.

    Andrews S. FastQC: a quality control tool for high throughput sequence data; 2010.

  6. 6.

    Angerer L, Hussain S, Wei Z, Livingston BT. Sea urchin metalloproteases: a genomic survey of the BMP-1/tolloid-like, MMP and ADAM families. Dev Biol. 2006;300(1):267–81.

  7. 7.

    Aouacheria A, Geourjon C, Aghajari N, Navratil V, Deleage G, Lethias C, et al. Insights into early extracellular matrix evolution: Spongin short chain collagen-related proteins are homologous to basement membrane type IV collagens and form a novel family widely distributed in invertebrates. Mol Biol Evol. 2006;23(12):2288–302.

  8. 8.

    Augustin R, Bosch TCG. Cnidarian Immunity: A Tale of Two Barriers. Boston: Springer; 2010. p. 1–16.

  9. 9.

    Autumn K, Peattie AM. Mechanisms of adhesion in geckos. Integr Comp Biol. 2002;42(6):1081–90.

  10. 10.

    Ba H, Yao F, Yang L, Qin T, Luan H, Li Z, et al. Identification and expression patterns of extracellular matrix-associated genes fibropellin-ia and tenascin involved in regeneration of sea cucumber Apostichopus japonicus. Gene. 2015;565(1):96–105.

  11. 11.

    Babonis LS, Martindale MQ, Ryan JF. Do novel genes drive morphological novelty? An investigation of the nematosomes in the sea anemone Nematostella vectensis. BMC Evol Biol. 2016;16(1):114.

  12. 12.

    Banfield DK, Irwin DM, Walz DA, MacGillivray RTA. Evolution of prothrombin: isolation and characterization of the cDNAs encoding chicken and hagfish prothrombin. J Mol Evol. 1994;38(2):177–87.

  13. 13.

    Baumgarten S, Simakov O, Esherick LY, Liew YJ, Lehnert EM, Michell CT, et al. The genome of Aiptasia, a sea anemone model for coral symbiosis. Proc Natl Acad Sci U S A. 2015;112(38):11893–8.

  14. 14.

    Bendtsen JD, Jensen LJ, Blom N, von Heijne G, Brunak S. Feature-based prediction of non-classical and leaderless protein secretion. Protein Eng Des Sel. 2004;17(4):349–56.

  15. 15.

    Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J Royal Stat Soc Ser B. 1995;57:289–300.

  16. 16.

    Bisgrove B, Andrews M, Raff R. Evolution of the fibropellin gene family and patterns of fibropellin gene expression in sea urchin phylogeny. J Mol Evol. 1995;41(1):34–45.

  17. 17.

    Bork P, Beckmann G. The CUB domain: a widespread module in developmentally regulated proteins. J Mol Biol. 1993;231(2):539–45.

  18. 18.

    Bosch TCG. The Path Less Explored: Innate Immune Reactions in Cnidarians. Berlin, Heidelberg: Springer; 2008. p. 27–42.

  19. 19.

    Buffet J-P, Corre E, Duvernois-Berthet E, Fournier J, Lopez PJ. Adhesive gland transcriptomics uncovers a diversity of genes involved in glue formation in marine tube-building polychaetes. Acta Biomater. 2018;72:316–28.

  20. 20.

    Bushnell B. BBMap – sourceforge.net/projects/bbmap/.; 2014.

  21. 21.

    Clarke DN, Miller PW, Lowe CJ, Weis WI, Nelson WJ. Characterization of the cadherin–catenin complex of the sea Anemone Nematostella vectensis and implications for the evolution of metazoan cell–cell adhesion. Mol Biol Evol. 2016;33(8):2016–29.

  22. 22.

    da Silva GFZ, Reuille RL, Ming L-J, Livingston BT. Overexpression and mechanistic characterization of blastula protease 10, a metalloprotease involved in sea urchin embryogenesis and development. J Biol Chem. 2006;281(16):10737–44.

  23. 23.

    Davey PA, Pernice M, Sablok G, Larkum A, Lee HT, Golicz A, et al. The emergence of molecular profiling and omics techniques in seagrass biology; furthering our understanding of seagrasses. Funct Integr Genomics. 2016;16(5):465–80.

  24. 24.

    DeMartini DG, Errico JM, Sjoestroem S, Fenster A, Waite JH. A cohort of new adhesive proteins identified from transcriptomic analysis of mussel foot glands. J R Soc Interface. 2017;14(131):20170151.

  25. 25.

    Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.

  26. 26.

    Dreanno C, Kirby RR, Clare AS. Smelly feet are not always a bad thing: the relationship between cyprid footprint protein and the barnacle settlement pheromone. Biol Lett. 2006;2(3):423–5.

  27. 27.

    DuBuc TQ, Ryan JF, Shinzato C, Satoh N, Martindale MQ. Coral comparative genomics reveal expanded Hox cluster in the cnidarian-Bilaterian ancestor. Integr Comp Biol. 2012;52(6):835–41.

  28. 28.

    DuBuc TQ, Traylor-Knowles N, Martindale MQ. Initiating a regenerative response; cellular and molecular features of wound healing in the cnidarian Nematostella vectensis. BMC Biol. 2014;12(1):24.

  29. 29.

    Ellis VL, Ross DM, Sutton L. The pedal disc of the swimming sea anemone Stomphia coccinea during detachment, swimming, and resettlement. Can J Zool. 1969;47(3):333–42.

  30. 30.

    Endrizzi BJ, Stewart RJ. Glueomics: an expression survey of the adhesive gland of the sandcastle worm. J Adhes. 2009;85(8):546–59.

  31. 31.

    Favi PM, Yi S, Lenaghan SC, Xia L, Zhang M. Inspiration from the natural world: from bio-adhesives to bio-inspired adhesives. J Adhes Sci Technol. 2014;28(3–4):290–319.

  32. 32.

    Flammang P, Demeuldre M, Hennebert E, Santos R. Adhesive secretions in echinoderms: a review. In: Biological adhesives. Cham: Springer International Publishing; 2016. p. 193–222.

  33. 33.

    Flammang R, Michel S, Cauwenberge V, Alexandre P, Jangoux M. A study of the temporary adhesion of the podia in the sea star Asterias rubens (Echinodermata, asteroidea) through their footprints. J Exp Biol. 1998;201(Pt 16):2383–95.

  34. 34.

    Foulon V, Boudry P, Artigaud S, Guérard F, Hellio C. In silico analysis of Pacific oyster (Crassostrea gigas) transcriptome over developmental stages reveals candidate genes for larval settlement. Int J Mol Sci. 2019;20(1):197.

  35. 35.

    Garrett TR, Bhakoo M, Zhang Z. Bacterial adhesion and biofilms on surfaces. Prog Nat Sci. 2008;18(9):1049–56.

  36. 36.

    Glazer L, Sagi A. On the involvement of proteins in the assembly of the crayfish gastrolith extracellular matrix. Invertebr Reprod Dev. 2012;56(1):57–65.

  37. 37.

    Gohad NV, Aldred N, Hartshorn CM, Jong Lee Y, Cicerone MT, Orihuela B, et al. Synergistic roles for lipids and proteins in the permanent adhesive of barnacle larvae. Nat Commun. 2014;5(1):4414.

  38. 38.

    Grajales A, Rodríguez E. Morphological revision of the genus Aiptasia and the family Aiptasiidae (Cnidaria, Actiniaria, Metridioidea); 2014.

  39. 39.

    Green KA, Lund LR. ECM degrading proteases and tissue remodelling in the mammary gland. BioEssays. 2005;27(9):894–903.

  40. 40.

    Haas BJ, Delcher AL, Mount SM, Wortman JR, Smith RK Jr, Hannick LI, et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assemblies. Nucleic Acids Res. 2003;31(19):5654–66.

  41. 41.

    Haseley SR, Vermeer HJ, Kamerling JP, Vliegenthart JF. Carbohydrate self-recognition mediates marine sponge cellular adhesion. Proc Natl Acad Sci U S A. 2001;98:9419–24.

  42. 42.

    Hennebert E, Maldonado B, Ladurner P, Flammang P, Santos R. Experimental strategies for the identification and characterization of adhesive proteins in animals: a review. Interface Focus. 2015;5(1):–20140064.

  43. 43.

    Hohenester E, Yurchenco PD. Laminins in basement membrane assembly. Cell Adhes Migr. 2013;7(1):56–63.

  44. 44.

    Hutter H, Vogel BE, Plenefisch JD, Norris CR, Proenca RB, Spieth J, et al. Conservation and novelty in the evolution of cell adhesion and extracellular matrix genes. Science. 2000;287(5455):989–94.

  45. 45.

    Hwang DS, Zeng H, Masic A, Harrington MJ, Israelachvili JN, Waite JH. Protein-and metal-dependent interactions of a prominent protein in mussel adhesive plaques. Biol Chem. 2010;285(33):25880–58.

  46. 46.

    Hynes RO, Zhao Q. The evolution of cell adhesion. J Cell Biol. 2000;150(2):F89–96.

  47. 47.

    Jones VAS, Bucher M, Hambleton EA, Guse A. Microinjection to deliver protein, mRNA, and DNA into zygotes of the cnidarian endosymbiosis model Aiptasia sp. Sci Rep. 2018;8(1):16437.

  48. 48.

    Kaur S, Weerasekare GM, Stewart RJ. Multiphase adhesive coacervates inspired by the sandcastle worm. ACS Appl Mater Interfaces. 2011;3(4):941–4.

  49. 49.

    Kitagawa M. Notch signalling in the nucleus: roles of Mastermind-like (MAML) transcriptional coactivators. J Biochem. 2015;159(3):mvv123.

  50. 50.

    Koch PJ, de Viragh PA, Scharer E, Bundman D, Longley MA, Bickenbach J, et al. Lessons from Loricrin-deficient mice. J Cell Biol. 2000;151(2):389–400.

  51. 51.

    Koizumi O. Origin and evolution of the nervous system considered from the diffuse nervous system of cnidarians. In: The Cnidaria, past, present and future. Cham: Springer International Publishing; 2016. p. 73–91.

  52. 52.

    Koutsogiannaki S, Kaloyianni M. Pollutants effects on ROS production and DNA damage of Mytilus galloprovincialis hemocytes: role of NHE, PKC and cAMP. Comp Biochem Physiol Part A Mol Integr Physiol. 2010;157(157):S38.

  53. 53.

    Kubow KE, Vukmirovic R, Zhe L, Klotzsch E, Smith ML, Gourdon D, et al. Mechanical forces regulate the interactions of fibronectin and collagen I in extracellular matrix. Nat Commun. 2015;6(1):8026.

  54. 54.

    Lengerer B, Hennebert E, Flammang P, Salvenmoser W, Ladurner P. Adhesive organ regeneration in Macrostomum lignano. BMC Dev Biol. 2016;16(1):20.

  55. 55.

    Li Y, Liew YJ, Cui G, Cziesielski MJ, Zahran N, Michell CT, et al. DNA methylation regulates transcriptional homeostasis of algal endosymbiosis in the coral model Aiptasia. Sci Adv. 2018;4(8):eaat2142.

  56. 56.

    Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30.

  57. 57.

    Ligtenberg AJM, Karlsson NG, Veerman ECI, Ligtenberg AJM, Karlsson NG, Veerman ECI. Deleted in malignant brain Tumors-1 protein (DMBT1): a pattern recognition receptor with multiple binding sites. Int J Mol Sci. 2010;11(12):5212–33.

  58. 58.

    Lin H-C, Wong YH, Tsang LM, Chu KH, Qian P-Y, 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(2):169–81.

  59. 59.

    Liu J, Li W, Liu S, Zheng X, Shi L, Zhang W, et al. Knockdown of collagen triple Helix repeat containing 1 (CTHRC1) inhibits epithelial-mesenchymal transition and cellular migration in glioblastoma cells. Oncol Res Featur Preclin Clin Cancer Ther. 2017;25(2):225–32.

  60. 60.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

  61. 61.

    Macrander J, Broe M, Daly M. Tissue-specific venom composition and differential gene expression in sea anemones. Genome Biol Evol. 2016;8(8):2358–75.

  62. 62.

    Malhotra JD, Kaufman RJ. The endoplasmic reticulum and the unfolded protein response. Semin Cell Dev Biol. 2007;18(6):716–31.

  63. 63.

    Miao Y, Zhang L, Sun Y, Jiao W, Li Y, Sun J, et al. Integration of transcriptomic and proteomic approaches provides a Core set of genes for understanding of scallop attachment. Mar Biotechnol. 2015;17(5):523–32.

  64. 64.

    Mollenhauer J, Herbertz S, Holmskov U, Tolnay M, Krebs I, Merlo A, et al. DMBT1 Encodes a Protein Involved in the Immune Defense and in Epithelial Differentiation and Is Highly Unstable in Cancer 1. Cancer Res. 2000;60:1704–10.

  65. 65.

    Moran Y, Praher D, Schlesinger A, Ayalon A, Tal Y, Technau U. Analysis of soluble protein contents from the nematocysts of a Model Sea Anemone sheds light on venom evolution. Mar Biotechnol. 2013;15(3):329–39.

  66. 66.

    Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.

  67. 67.

    Nichols SA, Dirks W, Pearse JS, King N. Early evolution of animal cell signaling and adhesion genes. Proc Natl Acad Sci U S A. 2006;103(33):12451–6.

  68. 68.

    Oakley CA, Ameismeier MF, Peng L, Weis VM, Grossman AR, Davy SK. Symbiosis induces widespread changes in the proteome of the model cnidarian Aiptasia. Cell Microbiol. 2016;18(7):1009–23.

  69. 69.

    Ochieng J, Furtak V, Lukyanov P. Extracellular functions of galectin-3. Glycoconj J. 2002;19(7–9):527–35.

  70. 70.

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

  71. 71.

    Oshiro N, Kobayashi C, Iwanaga S, Nozaki M, Namikoshi M, Spring J, et al. A new membrane-attack complex/perforin (MACPF) domain lethal toxin from the nematocyst venom of the Okinawan Sea anemone Actineria villosa. Toxicon. 2004;43(2):225–8.

  72. 72.

    Özbek S. The cnidarian nematocyst: a miniature extracellular matrix within a secretory vesicle. Protoplasma. 2011;248(4):635–40.

  73. 73.

    Ozsolak F, Milos PM. RNA sequencing: advances, challenges and opportunities. Nat Rev Genet. 2011;12(2):87–98.

  74. 74.

    Palmer CV, Traylor-Knowles N. Towards an integrated network of coral immune mechanisms. Proc R Soc B Biol Sci. 2012;279(1745):4106–14.

  75. 75.

    Paulsson M. Basement membrane proteins: structure, assembly, and cellular interactions. Crit Rev Biochem Mol Biol. 1992;27(1–2):93–127.

  76. 76.

    Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8(10):785–6.

  77. 77.

    Rampoldi L, Scolari F, Amoroso A, Ghiggeri G, Devuyst O. The rediscovery of uromodulin (Tamm–Horsfall protein): from tubulointerstitial nephropathy to chronic kidney disease. Kidney Int. 2011;80(4):338–47.

  78. 78.

    Richter K, Grunwald I, von Byern J. In: Adams RD, Öchsner A, da Silva LFM, editors. Handbook of Adhesion Technology. Berlin Heidelberg: Springer-Verlag; 2011. p. 1597–640.

  79. 79.

    Rocha M, Antas P, Castro LFC, Campos A, Vasconcelos V, Pereira F, et al. Comparative analysis of the adhesive proteins of the adult stalked goose barnacle Pollicipes pollicipes (Cirripedia: Pedunculata). Mar Biotechnol. 2019;21(1):38–51.

  80. 80.

    Rodrigues M, Ostermann T, Kremeser L, Lindner H, Beisel C, Berezikov E, et al. Profiling of adhesive-related genes in the freshwater cnidarian Hydra magnipapillata by transcriptomics and proteomics. Biofouling. 2016;32(9):1115–29.

  81. 81.

    Sarras MP Jr. Components, structure, biogenesis and function of the Hydra extracellular matrix in regeneration, pattern formation and cell differentiation. Int J Dev Biol. 2012;56(6–7–8):567–76.

  82. 82.

    Sarras MP, Yan L, Grens A, Zhang X, Agbas A, Huff JK, et al. Cloning and biological function of laminin in Hydra vulgaris. Dev Biol. 1994;164(1):312–24.

  83. 83.

    Satoh H, Oshiro N, Iwanaga S, Namikoshi M, Nagai H. Characterization of PsTX-60B, a new membrane-attack complex/perforin (MACPF) family toxin, from the venomous sea anemone Phyllodiscus semoni. Toxicon. 2007;49(8):1208–10.

  84. 84.

    Shi C-Y, Yang H, Wei C-L, Yu O, Zhang Z-Z, Jiang C-J, et al. Deep sequencing of the Camellia sinensis transcriptome revealed candidate genes for major metabolic pathways of tea-specific compounds. BMC Genomics. 2011;12(1):131.

  85. 85.

    Silverman HG, Roberto FF. Byssus formation in Mytilus. In: Biological adhesive systems. Vienna: Springer Vienna; 2010. p. 273–83.

  86. 86.

    Smith AM. The Role of Suction in the Adhesion of Limpets. J Exp Biol. 1991;161:1.

  87. 87.

    Smith AM, Papaleo C, Reid CW, Bliss JM. RNA-Seq reveals a central role for lectin, C1q and von Willebrand factor a domains in the defensive glue of a terrestrial slug. Biofouling. 2017;33(9):741–54.

  88. 88.

    So CR, Fears KP, Leary DH, Scancella JM, Wang Z, Liu JL, et al. Sequence basis of barnacle cement nanostructure is defined by proteins with silk homology. Sci Rep. 2016;6(1):36219.

  89. 89.

    Stamenkovic I. Extracellular matrix remodelling: the role of matrix metalloproteinases. J Pathol. 2003;200(4):448–64.

  90. 90.

    Steele RE, David CN, Technau U. A genomic view of 500 million years of cnidarian evolution. Trends Genet. 2011;27(1):7–13.

  91. 91.

    Stewart RJ, Wang CS. Adaptation of caddisfly larval silks to aquatic habitats by phosphorylation of H-fibroin Serines. Biomacromolecules. 2010;11(4):969–74.

  92. 92.

    Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms. PLoS One. 2011;6(7):e21800.

  93. 93.

    Surm JM, Smith HL, Madio B, Undheim EAB, King GF, Hamilton BR, et al. A process of convergent amplification and tissue-specific expression dominates the evolution of toxin and toxin-like genes in sea anemones. Mol Ecol. 2019;28(9):2272–89.

  94. 94.

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

  95. 95.

    Töpfer-Petersen E, Romero A, Varela PF, Ekhlasi-Hundrieser M, Dostàlovà Z, Sanz L, et al. Spermadhesins: a new protein family. Facts, hypotheses and perspectives. Andrologia. 2009;30(4–5):217–24.

  96. 96.

    Tucker RP, Adams JC. Adhesion networks of cnidarians: a postgenomic view. Int Rev Cell Mol Biol. 2014;308:323–77.

  97. 97.

    van Niel G, D’Angelo G, Raposo G. Shedding light on the cell biology of extracellular vesicles. Nat Rev Mol Cell Biol. 2018;19(4):213–28.

  98. 98.

    Vogel BE, Hedgecock EM. Hemicentin, a conserved extracellular member of the immunoglobulin superfamily, organizes epithelial and other cell attachments into oriented line-shaped junctions. Development. 2001;128:883–94.

  99. 99.

    Waite JH. Mussel adhesion-essential footwork; 2017.

  100. 100.

    Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.

  101. 101.

    Watanabe H, Fujisawa T, Holstein TW. Cnidarians and the evolutionary origin of the nervous system. Develop Growth Differ. 2009;51(3):167–83.

  102. 102.

    Whittaker CA, Bergeron K-F, Whittle J, Brandhorst BP, Burke RD, Hynes RO. The echinoderm adhesome. Dev Biol. 2006;300(1):252–66.

  103. 103.

    Whittaker CA, Hynes RO. Distribution and evolution of von Willebrand/integrin a domains: widely dispersed domains with roles in cell adhesion and elsewhere. Pollard TD, editor. Mol Biol Cell. 2002;13(10):3369–87.

  104. 104.

    Wright RM, Kenkel CD, Dunn CE, Shilling EN, Bay LK, Matz MV. Higher stress and immunity responses are associated with higher mortality in reef-building coral exposed to a bacterial challenge. bioRxiv. 2016:059444.

  105. 105.

    Wunderer J, Lengerer B, Pjeta R, Bertemes P, Kremser L, Lindner H, et al. A mechanism for temporary bioadhesion. Proc Natl Acad Sci U S A. 2019;116(10):4297–306.

  106. 106.

    Xu X, Xu M, Zhou X, Jones OB, Moharomd E, Pan Y, et al. Specific structure and unique function define the hemicentin. Cell Biosci. 2013;3(1):27.

  107. 107.

    Young GA, Yule AB, Walker G. Adhesion in the sea anemones actinia equina L. and metridium senile (L.). Biofouling. 1988;1(2):137–46.

  108. 108.

    Zeng F, Wunderer J, Salvenmoser W, Hess MW, Ladurner P, Rothbächer U. Papillae revisited and the nature of the adhesive secreting collocytes. Dev Biol. 2018.

  109. 109.

    Zhao H, Waite JH. Linking adhesive and structural proteins in the attachment plaque of Mytilus californianus. J Biol Chem. 2006;281(36):26150–8.

  110. 110.

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

Download references

Acknowledgements

Leeds University NGS Facility provided sequencing services and discussion. Access to the Newcastle University Rocket High-Performance Computing Cluster facilitated the bioinformatic analyses. We are grateful to Dr. Annika Guse of Heidelberg University for providing original broodstock of the CC7 Exaiptasia clonal line used in the study, and for her support throughout.

Funding

Newcastle University received financial support for the research from The Leverhulme Trust (Research Project Grant RPG-2017-017 to NA). NA was also supported by a Newcastle University Faculty Research Fellowship. Apart from providing financial support, funding bodies were not involved in the study design, data analyses, interpretation of results or manuscript preparation.

Author information

PD conducted the bioinformatic analyses, provided schematic graphics and wrote the manuscript. MR, NA and PD designed the research. MR extracted and prepared RNA for sequencing. JC maintained Exaiptasia pallida strain CC7 cultures. All authors contributed to and approved manuscript revisions.

Correspondence to Peter A. Davey.

Ethics declarations

Ethics approval and consent to participate

Ethics approval was not required for this study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Additional files

Additional file 1:

Bioinformatic analyses additional data. (XLSX 195 kb)

Additional file 2:

Schematic diagram detailing the bioinformatics workflow used in this study. (DOCX 44 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

Verify currency and authenticity via CrossMark

Keywords

  • Anemone
  • Adhesion
  • Exaiptasia
  • Bioadhesion
  • Cnidarian
  • mRNA-Seq
  • Transcriptomics