- Research article
- Open Access
Characterisation of the enzyme transport path between shipworms and their bacterial symbionts
BMC Biology volume 19, Article number: 233 (2021)
Shipworms are marine xylophagus bivalve molluscs, which can live on a diet solely of wood due to their ability to produce plant cell wall-degrading enzymes. Bacterial carbohydrate-active enzymes (CAZymes), synthesised by endosymbionts living in specialised shipworm cells called bacteriocytes and located in the animal’s gills, play an important role in wood digestion in shipworms. However, the main site of lignocellulose digestion within these wood-boring molluscs, which contains both endogenous lignocellulolytic enzymes and prokaryotic enzymes, is the caecum, and the mechanism by which bacterial enzymes reach the distant caecum lumen has remained so far mysterious. Here, we provide a characterisation of the path through which bacterial CAZymes produced in the gills of the shipworm Lyrodus pedicellatus reach the distant caecum to contribute to the digestion of wood.
Through a combination of transcriptomics, proteomics, X-ray microtomography, electron microscopy studies and in vitro biochemical characterisation, we show that wood-digesting enzymes produced by symbiotic bacteria are localised not only in the gills, but also in the lumen of the food groove, a stream of mucus secreted by gill cells that carries food particles trapped by filter feeding to the mouth. Bacterial CAZymes are also present in the crystalline style and in the caecum of their shipworm host, suggesting a unique pathway by which enzymes involved in a symbiotic interaction are transported to their site of action. Finally, we characterise in vitro four new bacterial glycosyl hydrolases and a lytic polysaccharide monooxygenase identified in our transcriptomic and proteomic analyses as some of the major bacterial enzymes involved in this unusual biological system.
Based on our data, we propose that bacteria and their enzymes are transported from the gills along the food groove to the shipworm’s mouth and digestive tract, where they aid in wood digestion.
In nature, there are many organisms able to degrade lignocellulose, having evolved to utilise plant carbohydrates for their metabolism through the production of plant cell wall-degrading enzymes . Bacteria and fungi are well known for their ability to digest lignocellulosic materials, and many terrestrial and marine invertebrates (including species of insects, nematodes, molluscs and crustaceans) are also able to produce endogenous cellulases. In most of these invertebrates, the digestion of woody biomass is made more efficient by the presence of bacterial or fungal symbionts, with only the isopods Limnoria ssp. and the amphipods Chelura terebrans being shown able to digest wood using endogenous enzymes alone [2,3,4]. Microbial symbionts playing these roles are normally found in the host in the digestive system, where the biomass is processed.
Shipworms (Fig. 1) are marine xylophagous bivalve molluscs belonging to the Teredinidae family, which can live on a diet of wood [5, 6]. They use modified shells (Fig. 1) to burrow into the wood, which is ground into small particles that are ingested and digested . Shipworms have an important role in the marine ecosystem, accounting for 70% of wood turnover in mangrove systems . Shipworms are also considered pests, since their destructive action on man-made wooden structures such as piers, boats or navigation poles has substantial negative economic impacts [9, 10]. The ability of shipworms to feed on lignocellulose is dependent on the presence of endosymbiotic bacteria that live in the animal’s gills in specialised eukaryotic cells called bacteriocytes [11, 12]. These bacteria provide the animal with fixed nitrogen to supplement its diet [13,14,15], as well as hydrolytic enzymes to help wood digestion and secondary metabolites such as antibiotics [16,17,18,19,20]. Shipworms differ from other bivalve molluscs in possessing a caecum (Fig. 1 and Fig. 2A), packed with wood particles, which opens from the posterior stomach and in certain species occupies a considerable portion of the body cavity. The caecum is the main site of wood digestion and contains a large number of carbohydrate-active enzymes (CAZymes) of both endogenous and bacterial origin  despite harbouring few bacteria, belonging to ribotypes different from those found in the gills . The prokaryotic fraction of CAZymes found in the caecum has been shown to be produced by bacteria resident in the spatially distant gills [23, 24], but until now the mechanism by which these bacterial enzymes end up in the lumen of the distant caecum has remained elusive. A duct connecting the gills to the oesophagus, located within the afferent branchial veins, has been reported in anatomical drawings of the shipworms Bankia gouldi and Teredo furcifera . However, despite numerous later studies, no further evidence supporting the existence of this “duct of Deshayes” has been published , suggesting that a yet unidentified mechanism may be responsible for the translocation of enzymes from the gills to the caecum.
In filter feeding molluscs, food particles suspended in the water are captured by the action of cilia and mucus on the surface of the gills (Additional file 1: Fig. S1A, B). These food particles are then channelled into the food groove (Fig. 2A–C) in a stream of mucus propelled by the movement of surface cilia, which transports them to the animal’s mouth and digestive system . In bivalves, an invagination of the stomach hosts a rotating gelatinous structure called the crystalline style (Fig. 2D, E), which is reported to accumulate digestive enzymes that are released by abrasion of the style against a chitinous area of the stomach wall, referred to as the gastric shield [28, 29]. The crystalline style has received little attention in shipworms, and its possible role in the interaction between host and symbionts has not been investigated.
In this paper, we present an in-depth investigation into the anatomical and molecular basis underpinning enzyme transport in shipworms and their bacterial endosymbionts. Based on data from meta-transcriptomic, meta-proteomic, X-ray micro-computed tomographic (micro-CT) and electron microscopy analyses, we propose that the food groove has been co-opted in shipworms to translocate bacteria and their secreted enzymes from the gills to the digestive system, and that the rotating crystalline style disrupts the incoming bacteria to further release their digestive enzymes, facilitating wood digestion. We also present in vitro biochemical characterisation of five prokaryotic lignocellulolytic enzymes identified in this complex symbiosis.
Micro-CT analysis and scanning electron microscopy
In order to explore the internal anatomy of L. pedicellatus, a high-resolution three-dimensional rendering of an adult specimen was created using micro-CT (Fig. 3), of an animal extracted from the wood. This technique avoided disruption of the delicate anatomical structures, which inevitably occurs during manual dissection, and provided insight into the animal’s complex internal features. By analysing the 3D model generated, we were unable to detect a duct connecting the gills and digestive system via the route proposed by Sigerfoos . Direct comparison with the illustrations of Sigerfoos was not possible due to significant modification of the gill structure in our specimen due to the presence of numerous brooded larvae. Previous descriptions of the duct were based on species that either do not brood larvae as part of their life cycle , or were not currently brooding at the time of study . Future research could examine developmental stages prior to the onset of brooding. However, examination of sections from the anterior adductor to posterior adductor muscle, a region which encompasses the mouth, oesophagus and stomach, also failed to reveal any duct associated with the afferent branchial vein (Additional file 2: Fig. S2). Micro-CT analysis therefore reveals no evidence for a “duct” by which bacterial enzymes could move from the gills to the caecum. Further, scanning electron microscopy (Additional file 3: Fig. S3) shows how the food groove provides a direct connection between the gills and the digestive system in L. pedicellatus, and ciliated tracts leading from the crystalline style into the caecum have been mapped using coloured particles in live L. pedicellatus specimens .
Bacterial CAZymes are found in the gills, caecum and food groove of L. pedicellatus
To identify the mechanism by which the bacterial CAZymes move from the gills to the digestive system, sections of the gill lamellae and food groove of L. pedicellatus were examined by TEM and SEM (transmission and scanning electron microscopy) in freshly harvested animals. These analyses showed eukaryotic cells (ranging from 10 to 50 μm and featuring organelles such as nuclei, mitochondria and Golgi apparatus) in the gill lamella tissues, which contained rod-shaped bacteria (Additional file 4: Fig. S4A-C), as previously described for bacteriocytes in other shipworm species . These bacteria have the characteristic appearance of gram-negative species, with an inner and outer membrane separated by the periplasmic space (Additional file 4: Fig. S4C), as reported in other species of endosymbiont-bearing molluscs [32, 33].
To establish the route by which bacterial CAZymes produced in the gills reach the shipworm digestive system, immunogold labelling studies were performed (Fig. 4) using antibodies (Additional file 5: Fig. S5) raised against a bacterial glycoside hydrolase family 5 (GH5) identified as abundant both in our transcriptomic and proteomic analysis of shipworm organs (LpsGH5_8, see next section on characterisation of bacterial CAZymes). TEM analysis of sectioned animals showed specific nanogold labelling in bacteriocytes in sections of the gill lamellae (Fig. 4D), and in the lumen of the food groove (halfway between gills and mouth, Fig. 4B), and in the lumen of the caecum (Fig. 4C), indicating the presence of the secreted protein. Nanogold particles in the food groove were found to be in close proximity to structures resembling, in size (~ 200 nm diameter, ~ 2 μm length), and features (inner and outer membrane, periplasmic space, granular content), the prokaryotic cells identified in the gills. In the caecum, the labelling was observed on and around the wood fragments (Fig. 4C), and no putative prokaryotic cells were visible. No labelling was observed on negative controls (pre-immune serum) (Fig. Additional file 6: Fig. S6).
Bacterial CAZyme genes are transcribed in the gills and their proteins reach the crystalline style before entering the caecum
The crystalline style, found in bivalve molluscs and some gastropods, is an acellular rod rich in glycoproteins, located inside the crystalline style sac  (Fig. 2D,E). The crystalline style is thought to help digestion by mixing digestive enzymes, reducing the size of food particles and directing them to the appropriate compartment for digestion [35,36,37,38,39,40,41,42,43,44]. The proteins that form the structure of the style are thought to be produced by the style sac , although little is known about their identity and function. The possible role of the crystalline style in the symbiosis between gill bacteria and shipworms has not previously been investigated.
Here we carried out shotgun proteomics and transcriptome sequencing of the L. pedicellatus crystalline style and style sac respectively, and discovered that a significant portion of the proteome of the style (14.7%) is represented by CAZymes of both prokaryotic and eukaryotic origin (Table 1). Furthermore, analysis of relative transcript abundance in the gills, digestive glands and crystalline style sac revealed that all genes coding for eukaryotic CAZyme proteins found in the style are transcribed in the animal’s digestive glands, while the prokaryotic ones are exclusively transcribed in the gills (Table 2). Table S1 (Additional file 7) gives details of the annotation obtained with BlastX searches against the NCBI non-redundant database for the proteins detected in the crystalline style.
Bacterial CAZymes found in the crystalline style proteome are dominated by members of the GH5 (27.3% of all the CAZymes) and GH134 (16.1%) families, with a lesser contribution of proteins from the GH16 (5.2%), GH10 (5.0%) and carbohydrate esterase (CE) 3 families (3.9%). Most of the prokaryotic CAZymes found in the crystalline style possess carbohydrate binding modules (CBMs) belonging to family 10. The main substrates for which these CAZymes and associated modules have been reported to be active, are described in Table S2 (Additional file 8). The most abundant animal proteins identified by proteomic analysis in the crystalline style are structural proteins, mucin (22.5%), tubulin (15.0%) and collagen (11.4%). The main eukaryotic CAZy families present in the style are GH9 (12.7% of the CAZymes), GH7 (7.7%) and GH5 (6.8%), while GH20, GH16, GH35, GH45 and GH2 family proteins are found at 4.3%, 4.1%, 2.7%, 2.0% and 0.7%, respectively. Other abundant eukaryotic proteins in the crystalline style are annotated as complement component (2.5%), hemocytin-like proteins (4.1%) or containing a membrane attack complex/perforin (MACPF) domain (6.8%), all of which have a major role in innate and adaptive immunity [46,47,48].
In order to investigate whether the CAZymes found in the crystalline style are transported to the caecum and likely contribute to wood digestion, we repeated the proteomic analysis of the caecum contents of L. pedicellatus performed by Sabbadin and colleagues , using a more modern mass spectrometer, which allowed a greater dynamic range of detection and therefore a more complete identification of the proteome analysed. We found that all the CAZymes identified in the crystalline style are also present in the caecum (Additional file 9: Table S3). Only 17.3% of the CAZymes detected in the caecum are not found in the crystalline style and their coding genes are transcribed in the caecum itself, downstream of the crystalline style.
Characterisation of CAZymes from gill bacterial symbionts
The combined results from micro-CT analysis, electron microscopy, transcriptomics and proteomics indicate that bacterial CAZymes are transported by the food groove from the gills to the stomach, where they are found in the crystalline style together with eukaryotic CAZymes and are finally released into the caecum by ciliary currents. In order to define the activities of some of the main bacterial CAZymes found in L. pedicellatus proteome and transcriptome, recombinant versions were produced and characterised in vitro. Full-length coding sequences of five enzymes were successfully cloned from RNA extracted from bacteria located in the gills, heterologously expressed and purified. Selected targets belong to GH family 5 subfamily 8 (LpsGH5_8, where Lps stands for Lyrodus pedicellatus symbionts; accession LS999939; antibodies raised against this enzyme were used for all immunolabelling experiments described earlier), GH11 (LpsGH11, accession LS999940), GH134 (LpsGH134a and LpsGH134b, accession LS999941 and LS999941 respectively) and Auxiliary Activity 10 (LpsAA10A, accession LS999942). All native target proteins carried an N-terminal signal peptide for secretion, at least one CBM, and contained up to 10 disulphide bonds (Fig. 5A and Additional file 10 and 11: Table S4 and S5). BlastP searches against the NCBI non-redundant database (NCBI nr) found best matches with proteins from the gammaproteobacteria Teredinibacter turnerae or Alteromonadaceae Bs08 and Bs12 strains, two of the genetically distinct ribotypes previously identified within Bankia setacea endosymbiont populations .
The purified recombinant GH enzymes were tested in vitro on a range of polysaccharide substrates using carbohydrate gel electrophoresis (PACE) (Fig. 5B) and DNS-reducing sugar assays (Additional file 12: Fig. S7).
The assays revealed that LpsGH5_8 (either with or without CBM) is active on a range of mannans, releasing mannose, mannobiose and longer mannooligosaccharides from glucomannan, locust bean gum (LBG), mannan and Scots pine wood.
Consistent with its predicted function as a β-1,4-xylanase, LpsGH11 expressed with both its CBMs showed activity on arabinoxylan and xylan. PACE analysis of the products released from grass xylan (miscanthus stem alcohol-insoluble residues) indicated the release of xylose, xylobiose and decorated xylooligosaccharides.
LpsGH134a and LpsGH134b, which are putative endo-β-1,4-mannanases based on similarity to known CAZymes, were mainly active on glucomannan, mannan, LBG and galactan. PACE analysis showed the absence of mannose monosaccharides (and disaccharides for LpsGH134a), suggesting that these bacterial hydrolases act as endo-mannanases requiring mannosyl residues at the − 2, − 1, + 1 and + 2 subsites of the enzyme.
Activity assays with purified, copper-bound LpsAA10A [lytic polysaccharide monooxygenases (LPMO) domain only] were carried out on a panel of cellulosic, hemicellulosic and chitinous substrates (Methods) in the presence of the electron donor gallic acid. Samples were analysed by MALDI-TOF MS and peak masses of the reaction products compared to previously published data [49, 50]. This revealed a predominant C1-oxidation pattern and generation of C1-aldonic acids on both phosphoric acid swollen cellulose (PASC, Fig. 5C) and microcrystalline cellulose (Avicel, Additional file 8: Fig. S8D) by LpsAA10A in presence of an external electron donor. Oxidised products were not detected in the negative controls (Additional file 13: Fig. S8A-C and E-G).
We used a range of microscopic, transcriptomic, proteomic and biochemical analyses to shed light on the anatomical adaptations and the dynamics of digestive enzyme distribution in Lyrodus pedicellatus. We have confirmed that the genes for bacterial enzymes found in the digestive caecum are transcribed in the gills , and our microscopy and micro-CT analyses indicate that the food groove is the most likely direct connection between the gills and the mouth (which leads into the caecum). LpsGH5_8 was detected in the crystalline style through proteomics analysis, and localised in the gills, as well as in the lumen of the food groove and of the caecum, using immunogold labelling. The density of nanogold particles found in the food groove lumen was much higher than in the gill bacteriocytes, suggesting that bacterial secretion of CAZymes mostly occurs after the bacteria have left the gills. Interestingly, the immunogold signal in the food groove was mostly localised around bacterial cells, indicating that both free bacteria and their secreted enzymes may be translocated via the food groove. Finally, the presence of nanogold particles, but not bacteria-like structures, in the lumen of the caecum suggested that bacterial cells are likely disrupted upon leaving the food groove, either in the mouth or within the stomach.
Our proteomics analysis also showed that eukaryotic and prokaryotic CAZymes account for 15% of all the proteins in the crystalline style, a rotating rod-shaped structure found in the animal’s stomach. The apparent accumulation of bacterial CAZymes in the shipworm style, and the lack of visible bacterial cells downstream of it, suggests that this rotating structure is responsible for disrupting incoming bacteria (as well as food particles) from the food groove, potentially releasing more wood digestive enzymes. Our analysis of the style found a highly abundant protein belonging to the Membrane Attack Complex/Perforin (MACPF) perforin superfamily, involved in punching holes in the membranes of Gram-negative bacteria as part of the animal immune system . The possible role of the style as a bacterial “grinder” upstream of the caecum is also indirectly supported by previous FISH investigations that showed only very sparse bacterial populations in the caecum .
Based on our observations and available literature, we suggest the following scenario (Fig. 6). Bacterial endosymbionts residing in the gill bacteriocytes are either passively or actively expelled from the gill’s tissue and incorporated in the mucous stream of food particles originated from filter feeding. Bacteria (together with food particles) are channelled into the food groove, propelled by the movement of surface cilia and they start releasing abundant CAZymes, possibly as a result of cell death induced by the host’s immune system. Through the food groove, bacteria are transported to the animal’s mouth and stomach, where the rotating crystalline style, enriched with bacteria-lysing MACPF perforin, grinds the incoming food particles, wood fragments and bacteria against the gastric shield, effectively releasing more CAZymes and mixing them with endogenous enzymes coming from the animal’s glands. The homogenised mixture of enzymes and food is transported by tracts of cilia  from the region of the style to the caecum, where the combined action of bacterial and eukaryotic CAZymes degrades lignocellulose and releases simple sugars that are taken up by the abundant glucose transporters (solute carrier family 2 transporters and sodium-dependent glucose transporters) previously identified in the caecum tissues .
In the present work, we have also characterised some of the major bacterial CAZymes identified in the gills, food groove and crystalline style of L. pedicellatus (belonging to GH families 5, 134 and 11, and the AA10 family), and identified activities on cellulose and hemicelluloses, particularly galacto-glucomannan and xylans. Our microscopy and shotgun proteomics studies have shown that GH5_8, a highly active galacto- and glucomannanase, is abundant both in the food groove and in the crystalline style. Strong hemicellulolytic activities are indeed to be expected in animals grown in blocks of Scots pine, a gymnosperm wood particularly rich in mannans . The results are in line with the findings of O’Connor and colleagues , who identified numerous mannanases, mannosidases and xylanases in the caecum contents of the shipworm B. setacea. We also managed to obtain the recombinant version of a bacterial AA10 lytic polysaccharide monooxygenase (LPMO) from L. pedicellatus and characterised its ability to cleave crystalline and amorphous cellulose at the C1 position, similarly to a previously characterised bacterial AA10 from Teredinibacter turnerae . Our combined omics analysis suggests that, although L. pedicellatus carries endogenous AA15 LPMO sequences, both their gene expression and protein abundance are very low in the digestive system (digestive glands, crystalline style, caecum) and that bacterial LPMOs have been co-opted towards wood digestion. This strategy is in striking contrast with the solution evolved by the firebrat Thermobia domestica, a terrestrial insect that exploits a complex array of endogenous LPMOs to digest cellulose .
Our work provides new insight into the evolution of complex symbioses between animals and intracellular bacteria. Although intracellular symbionts involved in supplementing nutrients to host are common among insects such as aphids and weevils [53, 54] and marine invertebrates such as vestimentiferan tubeworms and bivalve molluscs [55, 56], shipworms offer the most striking example of symbionts that are resident in organs spatially distant from the site of food digestion and nutrient uptake. This mechanism allows shipworms to take full advantage of the abundant carbohydrates found in wood without competition from symbiotic bacteria; however, it creates the logistical problem of transporting necessary nutrients and enzymes to a distant organ. Similar arrangements are rare in nature. One example is the tortoise beetle Cassidia rubiginosa, which hosts pectin-degrading bacteria within special organs of the reproductive tract . Shipworms have elegantly solved the problem by adopting a transport system that was already in place as part of their filter feeding behaviour. This evolutionary stratagem enabled them to fill an ecological niche in the marine ecosystem that was only partially exploited by other organisms. Our observation that bacteria appear to be transported along the food groove towards the digestive system in L. pedicellatus suggests that shipworms could be using this organ not only to relocate bacterial CAZymes, but also to recycle bacteria as a source of nutrients, particularly nitrogen, which is scarce in wood. Indeed, marine bivalves have been shown to feed on bacteria as a food source when phytoplankton is in short supply [58, 59]. It is possible that the use of the food groove as a transport system for bacterial products is widespread among bivalves and other molluscs, where endosymbiotic bacteria are commonly hosted in the gills . Recycling bacterial products by transporting them via the food groove could supplement their nutrition and allow them to thrive in extreme environments, such as deep marine hydrothermal vent ecosystems .
We have identified a novel putative mechanism of translocation of bacterial enzymes across distant organs in shipworms, and provided evidence that the food groove and the crystalline style may play key roles in it. This mechanism originates from a unique combination of anatomical and molecular features evolved by shipworms for the difficult task of digesting wood lignocellulose.
Animal rearing and dissection
Shipworms of the species Lyrodus pedicellatus from the Atlantic lineage  were used for our work. Samples were collected from the pier of Portsmouth Harbour (50° 47′ 47″ N, 1° 01′48″ W). Larvae from the original wood were used to infest logs of Scots pine, which were kept in tanks in the laboratories of the Institute of Marine Science, University of Portsmouth. To rear the animals, the water was taken directly from the Langstone Harbour (34 PSU salinity) and kept aerated and at a temperature of 15–18 °C using a flow-through system. The wood logs were opened by splitting them with a hammer and screwdriver, and the animals were then extracted with tweezers, placed in sea water containing EDTA-free protease inhibitors (1% v/v, Thermo Scientific) and kept on ice until dissection to anesthetise them. Species identification was performed using the pallets as described in . The dissections were performed using a stereomicroscope (Leica MZ6) after removing the mantle to expose the organs.
Scanning and transmission electron microscopy
Ten shipworms ranging in size from 4 to 8 cm were freshly extracted for all the microscopy experiments. For transmission electron microscopy, caecum, gills and food groove samples were fixed for 1–2 h at room temperature in primary fixative (4% formaldehyde (w/v), 2.5% (w/v) glutaraldehyde in 100 mM sodium phosphate buffer pH 7.2), and then washed (3 × 10 min) in 100 mM sodium phosphate buffer pH 7.2. Samples were then incubated in secondary fixative (1% osmium tetroxide in 100 mM sodium phosphate buffer pH 7.2) for 1 h on ice and dehydrated through a graded ethanol series (15 min each), followed by two washes (5 min each) in epoxy propane. Samples were infiltrated with a series of epoxy propane/Epon araldite (25%, 50%, 75% Epon Araldite with a minimum of 1 h at each stage, all at 30 °C) plus at least two changes of Epon araldite resin over 24 h at 30 °C, and polymerised at 60 °C for 48 h in flat embedding moulds. Pale gold (70–90 nm) ultra-thin sections were cut with a Diatome diamond knife using a Leica Ultracut UCT microtome and mounted on hexagonal 200-mesh nickel grids. Sections were post-stained with 2% (w/v) aqueous uranyl acetate (10 min) followed by lead citrate (5 min) in a carbon dioxide-free chamber and viewed using a FEI Tecnai 12 BioTWIN G2 TEM operating at 120 kV. Images were captured using AnalySIS software and a Megaview III CCD camera.
For scanning electron microscopy, the samples were fixed in 4% (v/v) glutaraldehyde in a cacodylate buffer (0.2 M sodium cacodylate, 0.3 M sodium chloride, 2 mM calcium chloride) for 2 h at room temperature and then rinsed once in buffer for 30 min. Samples were then taken through an ethanol dehydration series (50-70-100% ethanol and twice in 100% acetone, each stage for 30 min), critical point dried and then mounted on aluminium stubs using adhesive carbon tabs. Sputter coating was carried out under an argon atmosphere using a gold and palladium target, at a voltage of 1.4 kV using a current of approximately 18 mA for 3 min. Specimens were examined using a Zeiss MA10 Scanning Electron Microscope with an accelerating voltage of 20 kV and the Zeiss Smart software.
Embedding for immunolabeling proved difficult and several attempts were made to allow resin infiltration and at the same time preserve antigenicity. Freshly dissected shipworm tissues (food groove, gills and caecum) were fixed with 4% paraformaldehyde, 0.2% glutaraldehyde in sodium cacodylate buffer (0.2 M sodium cacodylate, 0.3 M sodium chloride, 2 mM calcium chloride, pH 7.4) on ice in a vacuum chamber for 2 h, then on a rotator without vacuum for a further 12 h at 4 °C. Samples were washed in 0.2 M sodium cacodylate buffer (three washes of 20 min each) and dehydrated through a graded ethanol series initially on ice (50%) and subsequently at − 20 °C on a rotator (70%, 90%, 100%) with 20 min at each stage and two changes of 100% ethanol. Ethanol was gradually replaced with LR Gold resin (1:2, 1:1, 2:1 resin:ethanol) with 1 h at each stage, followed by three changes of 100% LR Gold resin, 12 h each, all at − 20 °C on a rotator. Tissues were embedded in closed gelatine capsules and polymerised with UV light at − 20 °C for 24 h, followed by 24 h at − 10 °C. Pale gold (70–80 nm) ultra-thin sections were cut with a Diatome diamond knife, using a Leica Ultracut UCT microtome, and mounted on hexagonal 200-mesh nickel grids. All immunolabeling steps were achieved by floating grids on droplets of reagent. Sections were incubated in blocker (3% BSA in phosphate-buffered saline—PBS, 137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4 and 1.8 mM KH2PO4, pH 7.0) for 30 min at ambient temperature before incubation with primary antibodies against the bacterial LpsGH5_8 (see further down for details on antibodies production) diluted 1:100 in 1% BSA in PBS at 30 °C for 1 h, followed by washing with PBS at ambient temperature (three brief washes followed by three washes of 10 min each). Sections were incubated in secondary antibody (goat anti-rabbit IgG conjugated to 10 nm gold) diluted 1:100 in 1% BSA in PBS for 1 h at 30 °C, followed by washes with PBS as before, and subsequently with ultrapure water. All immunolabeling procedures included negative controls treated exactly as the samples: both pre-immune controls (diluted 1:100 in 1% BSA in PBS) and buffer only (1% BSA in PBS). Sections were post-stained with 2% (w/v) aqueous uranyl acetate (10 min), then lead citrate (5 min) in a carbon dioxide-free chamber and viewed using a FEI Tecnai 12 BioTWIN G2 operating at 120 kV. Images were captured using AnalySIS software and a Megaview III CCD camera.
Antibody production and purification
Two milligrams of purified recombinant bacterial LpsGH5_8 (see further sections for cloning, expression and purification) were used to raise polyclonal antibodies in rabbits (ProteoGenix, France). To enrich the serum for antigen-specific antibodies, affinity columns were made using recombinant LpsGH5_8. Recombinant protein preparations were dialysed against coupling buffer (0.1 M NaHCO3, 0.5 M NaCl, pH 8.3) and bound to CNBr-activated Sepharose™ 4 Fast Flow resin (GE Healthcare Life Sciences) followed by affinity purification of an aliquot of the crude antibody serum according to the resin manufacturer’s instructions. The pre-immune serum was subject to the same purification procedure. Purified antibody and pre-immune serum fractions were characterised for their affinity by western blotting using both recombinant LpsGH5_8 and L. pedicellatus caecum fluids (Additional file 5: Fig. S5). Fractions showing the highest titre and no unspecific binding were selected for immunogold labelling.
One adult Lyrodus pedicellatus (Quatrefages, 1849) specimen, measuring 3.6 cm in length and 1.62 cm in width, was used for micro-CT scanning. The specimen was reared at the Institute of Marine Sciences, University of Portsmouth, UK, and extracted from wood in 2012. The sample was fixed in 4 % v/v glutaraldehyde in a cacodylate buffer (0.2 M sodium cacodylate, 0.3 M sodium chloride, 2 mM calcium chloride) for 1 h at room temperature, rinsed three times in buffer for 10 min each, post-fixed in 1 % w/v aqueous osmium tetroxide for 1 h and rinsed three times in seawater for 10 min each. Samples were then immediately ethanol dehydrated and dried with hexamethyldisilazane (HMDS).
The specimen was mounted onto the sample holder and secured using glue, and scanned at the Ghent University Centre for X-ray Tomography (UGCT), Woodlab-UGent, using a scanner developed at UGCT. The scanner consisted of two X-ray tubes and two X-ray detectors, specifically designed to obtain very high-resolution scans as well as scans of larger objects. Scans were carried out using a microfocus X-ray tube in combination with a Varian flat-panel detector with an exposure time of 1500 ms, a rotation angle of 0.25° resulting in an average scan time of 45–60 min and an approximate voxel pitch of 2.5 μm. Details of the scanner are outlined in Masschaele et al.  and Van den Bulcke et al. [64, 65]. Due to large specimen size, two stacked scans were performed. The dataset was reconstructed using the Octopus software package with beam hardening correction. The two reconstructed volumes were then loaded in VGStudio MAX and stitched into a single stack of cross-sections. All resulting image and video analysis was performed using visualisation software myVGL.
RNA extraction and sequencing
Gills, digestive glands and caecum from three healthy adult L. pedicellatus were dissected and prepared for the paper published by Sabbadin and colleagues , using Ribosomal RNA depletion with a RiboZero™ Magnetic Gold Kit (Epidemiology) (Epicentre) in order to isolate both eukaryotic and prokaryotic mRNA. RNA-Seq libraries were prepared from each mRNA sample according to the Ion Total RNA-Seq kit v2 (Thermo Fisher Scientific). Templates were synthesised from mRNA libraries using the Ion OneTouch 200 Template Kit v2 DL on a OneTouch system (Thermo Fisher Scientific) and sequenced on an Ion Torrent PGM™ using a 318 chip (IonPGM200Kit; Thermo Fisher Scientific). All raw sequence data are available in NCBI under BioProject PRJNA412369 (SRA files: SRR6106265, SRR6106266, SRR6106267, SRR6106268, SRR6106269, SRR6106270, SRR6106271, SRR6106272, SRR6106273).
Crystalline style sacs were freshly dissected from 38 animals (which were then pooled together), flash frozen in liquid nitrogen and stored at − 80 °C. Total RNA was extracted using TRIzol® Reagent (Thermo Fisher Scientific), DNase treatment was carried out with Turbo DNA-free (Ambion), RNA was cleaned with RNA Clean & Concentrator™-5 (Zymo Research) and then quantified with a Qubit 3.0 Fluorometer and Agilent TapeStation. RNA depletion for both eukaryotic and prokaryotic ribosomal RNA was performed with the Ribo-ZeroTM Magnetic Gold Kit Epidemiology (Epicentre) and mRNA was then concentrated using RNA Clean & Concentrator™-5 (Zymo Research). The sequencing of the crystalline style sacs was performed at the Next Generation Sequencing Facility at the University of Leeds with HiSeq3000 using Illumina Technology to generate the required 150 bp paired end data. After rRNA depletion, library construction was completed using Illumina’s TruSeq stranded mRNA library protocol, starting at the RNA fragmentation step as suggested by Illumina.
Transcriptome assembly and analysis
Two meta-transcriptomes were assembled, using raw EST sequencing reads from digestive gland, caecum and gills (meta-transcriptome 1)  and from crystalline style sacs (meta-transcriptome 2), respectively. The raw EST sequencing reads were trimmed using Trimmmatic (part of the Galaxy tool [66, 67]) and were assembled into contigs using the Trinity software v. 2.8.3 .. Raw reads were mapped back onto the contigs and gene expression levels were calculated as TPM values (transcripts per kilobase million) using Salmon as part of the online tool Galaxy [66, 69] using standard parameters. Annotation of the contigs was performed by BlastX searches against the non-redundant database of the NCBI. The online software dbCAN (DataBase for automated Carbohydrate-active enzyme ANnotation)  was used to search for carbohydrate-active domains. Results with an e-value < 1e−10 and those with a CBM (Carbohydrate Binging Module) but no annotation were excluded, as well as CAZymes (Carbohydrate-Active enZYmes) belonging to the class of glycosyl transferases.
The following protocol describes the preparation and analysis of caecum and crystalline style samples. The caeca of five animals grown on Scots pine were dissected in 50 mM sodium phosphate buffer pH 7 and the content (food particles and enzymes) was isolated and pooled together. Similarly, 21 crystalline styles were dissected and washed three times in PBS buffer. Caecum and crystalline styles samples were then boiled for 10 min in denaturing buffer (1% SDS, 2.5% beta-mercapto ethanol, 175 mM DTT), centrifuged, and the supernatant was run into a 10% polyacrylamide gel. The protein bands were excised and digested with trypsin, and the resulting peptides were analysed by label-free LC-MS/MS. Tandem mass spectra were searched against the combined meta-transcriptomes (digestive gland, caecum, gills and crystalline style sac) of L. pedicellatus (which includes both eukaryotic and prokaryotic sequences) using the Mascot search programme. emPAI values were converted into molar percentages, the identified proteins were ranked based on relative abundance and annotated using BlastX versus non-redundant NCBI databases. CAZy annotation was carried out using the online tool dbCAN .
CAZy annotation was carried out using the online software dbCAN and results with an e-value < 1e−10, with a CBM but no CAZy module, and CAZymes belonging to the class of glycosyl transferases were excluded. Signal peptides were identified using the server SignalP 4.1 for either eukaryotes or gram-negative bacteria (www.cbs.dtu.dk/services/SignalP/).
Bacterial CAZymes cloning, expression, purification and biochemical characterisation
RNA was extracted from the gill tissue using the TRIzol® method (Thermo Fisher Scientific), cleaned with RNA Clean & Concentrator™-5 (Zymo Research) and a polyA tail was added using the poly(A) polymerase and protocol from Takara . cDNA was produced using the SuperScript® II reverse transcriptase (Thermo Fisher Scientific) with a oligo-dT primer and was purified with the Clontech NucleoSpin PCR Clean-up and gel extraction Kit (Clontech). The DNA sequences encoding the bacterial proteins LpsGH5_8, LpsGH11, LpsGH134a and LpsGH134b were amplified from cDNA without their signal peptide using the primers and PCR setting listed in Table S6 (Additional file 14) and they were cloned with the StrataClone Blunt PCR Cloning Kit (Stratagene); the sequences were verified by Sanger sequencing. LpsAA10A was not successfully cloned from the cDNA and therefore a synthetic version of the gene was codon optimised for E. coli expression by GeneArt.
The In-Fusion HD cloning kit (Takara) was used for cloning LpsGH5_8, LpsGH134a and LpsGH134b into the vector pET52b+, which has N-terminal Strep-Tag II followed by the human rhinovirus (HRV) 3C protease cleavage site, and a C-terminal His-Tag. The vectors were then transformed into E. coli Rosetta-GamiTM2(DE3) competent cells by heat shock.
LpsGH11 was cloned into the vector pOPINS3C , which contains an N-terminal His-Tag, followed by a Halo-Tag for improved soluble expression and the HRV 3C protease cleavage site, and no signal peptide for secretion. The vector was transformed into Spodoptera frugiperda 9 (Sf9) insect cells by heat shock.
LpsAA10A was cloned, without its CBM, using the In-Fusion HD cloning kit (Takara) into a modified pET26 vector containing at the N-terminus the pelB leader sequence to direct protein production to the periplasm, and a C-terminal Strep-tag. The construct was transformed into RosettaTM2(DE3) competent cells by heat shock.
Expression and purification
The E. coli bacterial cells containing LpsGH5_8, LpsGH134a and LpsGH134b constructs were grown in LB broth supplemented with carbenicillin (50 μg/ml) and chloramphenicol (34 μg/ml) at 37 °C until OD600 = 0.7 and then induced with isopropyl β-D-1-thiogalactopyranoside (IPTG) 1 mM and grown overnight at 20 °C and 200 rpm. After harvesting the cells were pelleted, suspended in phosphate-buffered saline (PBS) with 0.01 mM 4-(2-aminoethyl) benzenesulfonyl fluoride hydrochloride (AEBSF) and were lysed by sonication. After addition of 5 mM MgCl2 and DNaseI (0.025 U/μl) and filtering through a 0.45-μm filter, the supernatant was run through a 5-ml StrepTrap column, washed with PBS and eluted with 2.5 mM desthiobiotin in PBS. The eluted fractions containing absorbance peaks were analysed by SDS/PAGE to confirm the presence of the recombinant protein, combined together and the strep-tag was removed with the HRV Turbo protease (ABnova) at a ratio of 1:100 overnight at 4 °C and gentle shaking. After removal of the HRV, gel filtration was performed with a HiLoad 16/600 Superdex 75 pg column (Ge Healthcare) and the relevant peaks were verified by SDS/PAGE.
Expression of LpsGH11 in the Sf9 insect cells was performed using the baculovirus expression system  with a virus dilution of 1:1000. Once harvested, the cells were pelleted, resuspended in the lysis buffer I-PER insect cell protein extraction reagent (Thermo Fisher Scientific) with 5 mM MgCl2 and DNaseI (0.025 U/μl) and incubated on ice for 10 min. They were then pelleted and the supernatant was affinity purified with a pre-equilibrated HisTrap 5-ml column. Ni affinity chromatography was run with an elution gradient of 30 to 500 mM imidazole. Fractions were collected and analysed by SDS PAGE. Fractions containing the protein were pooled and concentrated to 2.5 ml and were run on a HiLoad 16/600 Superdex 75 pg column (Ge Healthcare) and the relevant peaks were verified by SDS/PAGE. The fractions containing the protein were concentrated to 1 ml and the tags were removed with the HRV Turbo protease (ABnova) at a ratio of 1:100 overnight at 4 °C and gentle shaking. Purification was carried out manually using a 5-ml HisTrap column and a gradient of 20–500 mM imidazole. The fractions were run by SDS-PAGE to confirm tag cleavage.
E. coli bacterial cells expressing the LpsAA10A were grown in M9 Minimal Medium, containing 1% (w/v) glucose and the appropriate antibiotics at 37 °C until OD600 = 0.7. The culture was induced with IPTG (0.1 mM) and grown overnight at 20 °C. Cells were harvested by centrifugation, resuspended in 50 ml of 50 mM Tris-HCl/20% sucrose (pH 8.0) for each litre of original culture and kept on ice for 30 min. After centrifugation at 8000 rpm for 10 min the supernatant was discarded, the cells were resuspended in 50 ml of 5 mM MgSO4 for each litre of original culture and kept on ice for 30 min. After centrifugation the supernatant, containing the periplasmic fraction, was equilibrated with 0.2 M Na phosphate buffer pH 7.6 to a final concentration of 50 mM, applied to a 5-ml StrepTrap HP column, washed with binding buffer and eluted with 2.5 mM desthiobiotin. 5-fold excess copper was added as CuSO4, then unbound copper and desthiobiotin were removed by passing the protein in a HiLoad TM 16/600 Superdex 75 gel filtration column (Ge Healthcare) equilibrated with 10 mM sodium phosphate buffer pH 7.0. The protein was then concentrated using Microsep TM Advance Centrifugal Devices (Pall Corporation).
Substrates used for the DNS assay, or PACE: barley β-glucan (β-D-1,3-1,4-glucan), mannan (borohydride reduced), konjac glucomannan (β-D-1,4), larch arabinogalactan, wheat arabinoxylan, tamarind seed xyloglucan, potato galactan and galactan (Gal:Ara:Rha:Xyl:GalUA = 91:2:1.7:0.3:5), are all purchased from Megazyme; locust bean gum (LBG), carboxymethyl-cellulose (CMC), microcrystalline cellulose (Avicel) and beech wood xylan are purchased from Sigma-Aldrich. Phosphoric acid swollen cellulose (PASC) was prepared as in . Grass xylan (miscanthus stem alcohol-insoluble residues) was prepared as described in .
DNS-reducing sugar assays
The activity of LpsGH5_8, LpsGH11, LpsGH134a and LpsGH134b was determined by dinitrosalicylic acid (DNS)-reducing sugar assay on a range of polysaccharides (see paragraph “Substrates”). The 50-μl reactions were carried out in triplicates in 50 mM sodium phosphate buffer pH 6.0, 0.1% substrate and 3 μg of protein (0.5 μg for LpsGH11). They were incubated at 30 °C for 2 h with shaking at 320 rpm and then 9 μl of the reaction was added to 31 μl of DNS reagent and heated at 100 °C for 5 min. After cooling at room temperature and addition of 160 μl water, the 540 nm absorbance was measured in a micro-plate reader and the results were compared to a glucose standard curve. The A540 of the substrates was subtracted from that of the samples. The DNS reagent was prepared by mixing 0.75 g of dinitrosalycilic acid, 1.4 g NaOH, 21.6 g sodium potassium tartrate tetrahydrate, 0.53 mL phenol and 0.59 g sodium metabisulfite in 100 ml of distilled water, and it was filtered and kept in the dark before used.
Product analysis by mass spectrometry (MS)
Reactions with the purified LpsAA10A were carried out by mixing 4 mg mL − 1 substrate with purified copper-loaded enzyme (2 μM) and 4 mM electron donor (gallic acid), in 50 mM ammonium acetate buffer pH 6 in 2-mL plastic reaction tubes (reaction volume: 100 μL). The tubes were incubated for 24 h at 28 °C shaking at 1000 rpm, centrifuged at 14,000 rpm and the supernatant was collected for analysis through mass spectrometry. Briefly, 1 μl of supernatant was mixed with an equal volume of matrix solution (20 mg mL− 1 2,5-dihydroxybenzoic acid (DHB) in 50% acetonitrile plus 0.1% TFA), spotted on a SCOUT-MTP 384 target plate (Bruker) and analysed by positive-mode MALDI-TOF MS using an Ultraflex III matrix-assisted laser desorption ionisation-time of flight/time of flight (MALDI/TOF-TOF) instrument (Bruker).
Polysaccharide analysis by carbohydrate gel electrophoresis (PACE)
Purified enzyme at 20 μg/ml was mixed with 0.5% galactan, glucomannan, galactomannan, mannan or locust bean gum (LBG) or with 40 mg/ml of milled Scots pine wood (pre-treated in 0.5 N NaOH for 30 min at 90 °C and rinsed 5 times in 50 mM NaPO4 buffer) in 50 mM NaPO4 buffer pH 6.5 and incubated overnight at 30 °C with shaking. The samples were then centrifuged, supernatant was transferred to a new tube and undigested polysaccharides were removed by precipitation with 80% ethanol. Following centrifugation, supernatants were transferred to a new tube and dried.
Miscanthus stem AIR (alcohol-insoluble residues) was pre-treated in 4 M NaOH for 1 h at RT and neutralised with HCl. Resultant substrate at 1 mg/ml (of initial untreated AIR) was digested overnight at RT with various amounts of xylanase (3–40 μg/ml). All samples were purified on Nanosep 10 K and dried. Dried digestion products and manno-oligosaccharide and xylo-oligosaccharide standards and appropriate controls were labelled with 8-aminonaphthalene-1,3,6-trisulfonic acid (ANTS; Invitrogen, www.invitrogen.com) and separated by polyacrylamide gels, as described previously . PACE gels were visualised using a G-box (Syngene, www.syngene.com/). Experiments were carried out in triplicate, and the representative gels are shown.
Availability of data and materials
All data generated or analysed during this study are included in this published article, its supplementary information files and publicly available repositories or are otherwise available in the following repositories:
Proteomics data: Pesante, G., Sabbadin, F., Dowle, A-A. Analysis of the proteome of the crystalline style of the shipworm Lyrodus pedicellatus. MassIVE (Mass Spectrometry Interactive Virtual Environment).
Transcriptomics data: Pesante, G., Sabbadin, F., Li, Y. Lyrodus pedicellatus crystalline style sac RNA sequencing. NCBI (National Center for Biotechnology Information). BioProject ID PRJNA450378, SRA accession: SRP140494, www.ncbi.nlm.nih.gov/bioproject/450378.
Bacterial protein accessions: Pesante, G., Sabbadin, F., Elias, L. Project PRJEB28739. ENA (European Nucleotide Archive), www.ebi.ac.uk/ena/data/view/LS999939-LS999943.
MicroCT data: Shipway, J-R, Cragg, S-M. Proposed mechanism of enzyme transport between shipworms and their bacterial symbionts. Figshare.
Cragg SM, Beckham GT, Bruce NC, Bugg TDH, Distel DL, Dupree P, et al. Lignocellulose degradation mechanisms across the Tree of Life. Curr Opin Chem Biol. 2015;29:108–19. https://doi.org/10.1016/j.cbpa.2015.10.018.
Sleeter TD, Boyle PJ, Cundell AM, Mitchell R. Relationships between marine microorganisms and wood-boring isopod Limnoria-Tripunctata. Mar Biol. 1978;45(4):329–36. https://doi.org/10.1007/BF00391819.
Daniel G, Nilsson T, Cragg S. Limnoria-Lignorum ingest bacterial and fungal degraded wood. Holz Als Roh-Und Werkstoff. 1991;49(12):488–90. https://doi.org/10.1007/BF02619480.
King AJ, Cragg SM, Li Y, Dymond J, Guille MJ, Bowles DJ, et al. Molecular insight into lignocellulose digestion by a marine isopod in the absence of gut microbes. Proc Natl Acad Sci U S A. 2010;107(12):5345–50. https://doi.org/10.1073/pnas.0914228107.
Becker G. Biological investigations on marine borers in Berlin-Dahlem. In: Ray DL, editor. Marine Boring and Fouling Organisms. Seattle: University of Washington Press; 1959. p. 62–76.
Gallager SM, Turner RD, Berg CJ. Physiological aspects of wood consumption, growth, and reproduction in the shipworm Lyrodus pedicellatus Quatrefages (Bivalvia, Teredinidae). J Exp Mar Biol Ecol. 1981;52(1):63–77. https://doi.org/10.1016/0022-0981(81)90171-4.
Mann R. In: Costlow JD, Tipper RC, editors. Nutrition in the Teredinidae. Marine biodeterioration: an interdisciplinary study; 1984. p. 24–9.
Distel DL. The biology of marine wood boring bivalves and their bacterial endosymbionts. Acs Sym Ser. 2003;845:253–71. https://doi.org/10.1021/bk-2003-0845.ch014.
Turner RD. A survey and illustrated catalogue of the Teredinidae (Mollusca: Bivalvia). A survey and illustrated catalogue of the Teredinidae (Mollusca: Bivalvia); 1966.
Distel DL, Amin M, Burgoyne A, Linton E, Mamangkey G, Morrill W, et al. Molecular phylogeny of Pholadoidea Lamarck, 1809 supports a single origin for xylotrophy (wood feeding) and xylotrophic bacterial endosymbiosis in Bivalvia. Mol Phylogenet Evol. 2011;61(2):245–54. https://doi.org/10.1016/j.ympev.2011.05.019.
Popham JD, Dickson MR. Bacterial associations in the Teredo Bankia australis (Lamellibranchia: Mollusca). Mar Biol. 1973;19(4):338–40. https://doi.org/10.1007/BF00348904.
Distel DL, Beaudoin DJ, Morrill W. Coexistence of multiple proteobacterial endosymbionts in the gills of the wood-boring bivalve Lyrodus pedicellatus (Bivalvia : Teredinidae). Appl Environ Microbiol. 2002;68(12):6292–9. https://doi.org/10.1128/AEM.68.12.6292-6299.2002.
Carpenter EJ, Culliney JL. Nitrogen fixation in marine shipworms. Science. 1975;187(4176):551–2. https://doi.org/10.1126/science.187.4176.551.
Waterbury JB, Calloway CB, Turner RD. A cellulolytic nitrogen-fixing bacterium cultured from the gland of Deshayes in shipworms (Bivalvia: Teredinidae). Science. 1983;221(4618):1401–3. https://doi.org/10.1126/science.221.4618.1401.
Lechene CP, Luyten Y, McMahon G, Distel DL. Quantitative imaging of nitrogen fixation by individual bacteria within animal cells. Science. 2007;317(5844):1563–6. https://doi.org/10.1126/science.1145557.
Distel DL, Morrill W, MacLaren-Toussaint N, Franks D, Waterbury J. Teredinibacter turnerae gen. nov., sp. nov., a dinitrogen-fixing, cellulolytic, endosymbiotic gamma-proteobacterium isolated from the gills of wood-boring molluscs (Bivalvia: Teredinidae). Int J Syst Evol Microbiol. 2002;52(Pt 6):2261–9. https://doi.org/10.1099/00207713-52-6-2261.
Elshahawi SI, Trindade-Silva AE, Hanora A, Han AW, Flores MS, Vizzoni V, et al. Boronated tartrolon antibiotic produced by symbiotic cellulose-degrading bacteria in shipworm gills. Proc Natl Acad Sci U S A. 2013;110(4):E295–304. https://doi.org/10.1073/pnas.1213892110.
Han AW, Sandy M, Fishman B, Trindade-Silva AE, Soares CA, Distel DL, et al. Turnerbactin, a novel triscatecholate siderophore from the shipworm endosymbiont Teredinibacter turnerae T7901. PLoS ONE. 2013;8(10):e76151. https://doi.org/10.1371/journal.pone.0076151.
Haygood M, Altamia M, Elshahawi S, Han A, Lin Z, Betcher M, et al. Versatile bacterial symbionts of shipworms contribute to wood digestion, fix nitrogen and produce secondary metabolites. Planta Medica. 2015;81(11):IL43. https://doi.org/10.1055/s-0035-1556140.
Brito TL, Campos AB, Bastiaan von Meijenfeldt FA, Daniel JP, Ribeiro GB, Silva GGZ, et al. The gill-associated microbiome is the main source of wood plant polysaccharide hydrolases and secondary metabolite gene clusters in the mangrove shipworm Neoteredo reynei. PLoS ONE. 2018;13(11):e0200437. https://doi.org/10.1371/journal.pone.0200437.
Sabbadin F, Pesante G, Elias L, Besser K, Li Y, Steele-King C, et al. Uncovering the molecular mechanisms of lignocellulose digestion in shipworms. Biotechnol Biofuels. 2018;11(1):59. https://doi.org/10.1186/s13068-018-1058-3.
Betcher MA, Fung JM, Han AW, O'Connor R, Seronay R, Concepcion GP, et al. Microbial distribution and abundance in the digestive system of five shipworm species (Bivalvia: Teredinidae). PLoS ONE. 2012;7(9):e45309. https://doi.org/10.1371/journal.pone.0045309.
O'Connor RM, Fung JM, Sharp KH, Benner JS, McClung C, Cushing S, et al. Gill bacteria enable a novel digestive strategy in a wood-feeding mollusk. Proc Natl Acad Sci U S A. 2014;111(47):E5096–104. https://doi.org/10.1073/pnas.1413110111.
Ekborg NA, Morrill W, Burgoyne AM, Li L, Distel DL. CelAB, a multifunctional cellulase encoded by Teredinibacter turnerae T7902T, a culturable symbiont isolated from the wood-boring marine bivalve Lyrodus pedicellatus. Appl Environ Microbiol. 2007;73(23):7785–8. https://doi.org/10.1128/AEM.00876-07.
Sigerfoos CP. Natural history, organization, and late development of the Teredinidae, or ship-worms. Bull U S Bur Fisheries. 1908;27:191–231. https://doi.org/10.5962/bhl.title.24371.
De Moraes DT, Lopes SG. The functional morphology of Neoteredo reynei (Bartsch, 1920)(Bivalvia, Teredinidae). J Molluscan Stud. 2003;69(4):311–8. https://doi.org/10.1093/mollus/69.4.311.
Gosling E. Marine bivalve molluscs. Oxford: John Wiley & Sons; 2015. https://doi.org/10.1002/9781119045212.
McQuiston R. Fine structure of the gastric shield in the lamellibranch bivalve, Lasaea rubra (Montagu). J Molluscan Stud. 1970;39(1):69–75.
Sakamoto K, Uji S, Kurokawa T, Toyohara H. Immunohistochemical, in situ hybridization and biochemical studies on endogenous cellulase of Corbicula japonica. Comp Biochem Physiol B Biochem Mol Biol. 2008;150(2):216–21. https://doi.org/10.1016/j.cbpb.2008.03.004.
Nair NB, Saraswathy M. The biology of wood-boring teredinid molluscs. Adv Mar Biol. 1971;9:335–509. https://doi.org/10.1016/S0065-2881(08)60345-4.
Morton B. The functional anatomy of the organs of feeding and digestion of Teredo navalis Linnaeus and Lyrodus pedicellatus (Quatrefages). Proc MalacSoc Lond. 1970;39(2-3):151–68. https://doi.org/10.1093/oxfordjournals.mollus.a065089.
Distel D, Felbeck HJMB. Endosymbiosis in the lucinid clams Lucinoma aequizonata, Lucinoma annulata and Lucina floridana: a reexamination of the functional morphology of the gills as bacteria-bearing organs. Mar Biol. 1987;96(1):79–86.
Distel DL, Roberts SJ. Bacterial endosymbionts in the gills of the deep sea woodboring bivalves Xylophaga atlantica and Xylophaga washingtona. Biol Bull (Woods Hole). 1997;192(2):253–61. https://doi.org/10.2307/1542719.
Mackenzie AL, Marshall C. Proteins in the crystalline styles of the marine mussels Perna canaliculus Gmelin and Mytilus galloprovincialis Lamarck. J Shellfish Res. 2014;33(3):673–85. https://doi.org/10.2983/035.033.0302.
Nelson TC. On the origin, nature, and function of the crystalline style of lamellinranchs. J Morphol. 1917;31(1):53–111. https://doi.org/10.1002/jmor.1050310105.
Edmondson CH. The reformation of the crystalline style in Mya arenaria after extraction. J Exp Zool. 1920;30(3):259–91. https://doi.org/10.1002/jez.1400300302.
Yonge CM. Studies on the comparative physiology of digestion. Br J Exp Zool. 1931;1(1):15–63. https://doi.org/10.1242/jeb.1.1.15.
Lavine TF. A study of the enzymatic and other properties of the crystalline style of clams; evidence for the presence of a cellulase. J Cell Physiol. 1946;28(2):183–95. https://doi.org/10.1002/jcp.1030280203.
Morton J. The role of the crystalline style. J Molluscan Stud. 1952;29(2-3):85–92.
Bailey K, Worboys BD. The lamellibranch crystalline style. Biochem J. 1960;76(3):487–91. https://doi.org/10.1042/bj0760487.
Horiuchi S, Lane EC. Carbohydrases of the crystalline style and hepatopancreas of Strombus gigas Linné. Comp Biochem Physiol. 1966;17(4):1189–97. https://doi.org/10.1016/0010-406X(66)90294-5.
Kristensen JH. Structure and function of crystalline styles of bivalves. Ophelia. 1972;10(1):91–108. https://doi.org/10.1080/00785326.1972.10430106.
Judd W. Crystalline style proteins from bivalve mollusks. Comp Biochem Physiol, B Biochem Mol Biol. 1987;88(1):333–9. https://doi.org/10.1016/0305-0491(87)90124-6.
Alyakrinskaya IO. The dimentions, characteristics and functions of the crystalline style of molluscs. Biol Bull Acad Sci USSR. 2001;28(5):523–35. https://doi.org/10.1023/A:1016756629952.
Judd W. The secretions and fine structure of bivalve crystalline style sacs. Ophelia. 1979;18(2):205–33. https://doi.org/10.1080/00785326.1979.10425500.
Peitsch MC, Tschopp J. Assembly of macromolecular pores by immune defense systems. Curr Opin Cell Biol. 1991;3(4):710–6. https://doi.org/10.1016/0955-0674(91)90045-Z.
Rosado CJ, Kondos S, Bull TE, Kuiper MJ, Law RH, Buckle AM, et al. The MACPF/CDC family of pore-forming toxins. Cell Microbiol. 2008;10(9):1765–74. https://doi.org/10.1111/j.1462-5822.2008.01191.x.
Yamakawa M, Tanaka H. Immune proteins and their gene expression in the silkworm, Bombyx mori. Dev Comp Immunol. 1999;23(4-5):281–9. https://doi.org/10.1016/S0145-305X(99)00011-7.
Quinlan RJ, Sweeney MD, Leggio LL, Otten H, Poulsen J-CN, Johansen KS, et al. Insights into the oxidative degradation of cellulose by a copper metalloenzyme that exploits biomass components. Proc Natl Acad Sci. 2011;108(37):15079–84. https://doi.org/10.1073/pnas.1105776108.
Sabbadin F, Hemsworth GR, Ciano L, Henrissat B, Dupree P, Tryfona T, et al. An ancient family of lytic polysaccharide monooxygenases with roles in arthropod development and biomass digestion. Nat Commun. 2018;9(1):756. https://doi.org/10.1038/s41467-018-03142-x.
Fengel D, Wegener G. Wood: chemistry, ultrastructure, reactions. Berlin: Walter de Gruyter; 1989.
Fowler CA, Sabbadin F, Ciano L, Hemsworth GR, Elias L, Bruce N, et al. Discovery, activity and characterisation of an AA10 lytic polysaccharide oxygenase from the shipworm symbiont Teredinibacter turnerae. Biotechnol Biofuels. 2019;12(1):12(1). https://doi.org/10.1186/s13068-019-1573-x.
Buchner P. Endosymbiosis of animals with plant microorganims; 1965.
Douglas A. Mycetocyte symbiosis in insects. Biol Rev. 1989;64(4):409–34. https://doi.org/10.1111/j.1469-185X.1989.tb00682.x.
Savazzi E. A review of symbiosis in the Bivalvia, with special attention to macrosymbiosis. Paleontol Res. 2001;5(1):55–73.
Roeselers G, Newton IL. On the evolutionary ecology of symbioses between chemosynthetic bacteria and bivalves. Appl Microbiol Biotechnol. 2012;94(1):1–10. https://doi.org/10.1007/s00253-011-3819-9.
Salem H, Bauer E, Kirsch R, Berasategui A, Cripps M, Weiss B, et al. Drastic genome reduction in an herbivore’s pectinolytic symbiont. Cell. 2017;171(7):1520–+.
Mchenery JG, Birkbeck TH. Lysozyme of the mussel, Mytilus-Edulis (L). Mar Biol Lett. 1979;1(2):111–9.
Langdon CJ, Newell RIE. Utilization of detritus and bacteria as food sources by 2 bivalve suspension-feeders, the oyster Crassostrea-Virginica and the mussel Geukensia-Demissa. Mar Ecol Prog Ser. 1990;58(3):299–310. https://doi.org/10.3354/meps058299.
Distel DL, Baco AR, Chuang E, Morrill W, Cavanaugh C, Smith CR. Marine ecology - do mussels take wooden steps to deep-sea vents? Nature. 2000;403(6771):725–6. https://doi.org/10.1038/35001667.
Borges LMS, Merckelbach LM, Sampaio I, Cragg SM. Diversity, environmental requirements, and biogeography of bivalve wood-borers (Teredinidae) in European coastal waters. Front Zool. 2014;11(1):13. https://doi.org/10.1186/1742-9994-11-13.
Turner RD. Identification of marine wood-boring molluscs. In: EBG J, Eltringham SK, OECD, editors. Marine Borers and Fouling Organisms of Wood. Paris: Organisation for economic co-operation and development; 1971. p. 17–64.
Masschaele BC, Cnudde V, Dierick M, Jacobs P, Van Hoorebeke L, Vlassenbroeck J. UGCT: new X-ray radiography and tomography facility. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment. 2007;580(1):266–9. https://doi.org/10.1016/j.nima.2007.05.099.
Van den Bulcke J, De Windt I, Defoirdt N, De Smet J, Van Acker J. Moisture dynamics and fungal susceptibility of plywood. Int Biodeterioration Biodegradation. 2011;65(5):708–16. https://doi.org/10.1016/j.ibiod.2010.12.015.
Van den Bulcke J, De Windt I, Defoirdt N, Van Acker J. Non-destructive evaluation of wood decay. In: 42nd Annual meeting of the International Research Group on Wood Protection (IRG/WP); 2011. International Research Group on Wood Protection.
Afgan E, Baker D, van den Beek M, Blankenberg D, Bouvier D, Cech M, et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update. Nucleic Acids Res. 2016;44(W1):W3–W10. https://doi.org/10.1093/nar/gkw343.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat Biotechnol. 2011;29(7):644–52. https://doi.org/10.1038/nbt.1883.
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–+.
Yin YB, Mao XZ, Yang JC, Chen X, Mao FL, Xu Y. dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012;40(W1):W445–51. https://doi.org/10.1093/nar/gks479.
Sippel AE. Purification and characterization of adenosine triphosphate: ribonucleic acid adenyltransferase from Escherichia coli. Eur J Biochem. 1973;37(1):31–40. https://doi.org/10.1111/j.1432-1033.1973.tb02953.x.
Bird LE. High throughput construction and small scale expression screening of multi-tag vectors in Escherichia coli. Methods. 2011;55(1):29–37. https://doi.org/10.1016/j.ymeth.2011.08.002.
Berger I, Fitzgerald DJ, Richmond TJ. Baculovirus expression system for heterologous multiprotein complexes. Nat Biotechnol. 2004;22(12):1583–7. https://doi.org/10.1038/nbt1036.
Tryfona T, Sorieul M, Feijao C, Stott K, Rubtsov DV, Anders N, et al. Development of an oligosaccharide library to characterise the structural variation in glucuronoarabinoxylan in the cell walls of vegetative tissues in grasses. Biotechnol Biofuels. 2019;12(1):109. https://doi.org/10.1186/s13068-019-1451-6.
Bromley JR, Busse-Wicher M, Tryfona T, Mortimer JC, Zhang Z, Brown DM, et al. GUX1 and GUX2 glucuronyltransferases decorate distinct domains of glucuronoxylan with different substitution patterns. Plant J. 2013;74(3):423–34. https://doi.org/10.1111/tpj.12135.
The staff of the Institute of Marine Sciences (ISM) in Portsmouth gave invaluable help with the rearing, retrieval and dissection of the shipworms. The Oxford Protein Production Facility was very helpful in protein production. Thank you also to Jan Van den Bulcke, Ghent University, for performing the micro-CT scan on the shipworm specimen. Acknowledgements to Rachel Bates and Meg Stark from the Bioscience Technology Facility in the University of York for proteomics analysis and imaging; to Ian Carr, Sally Harrison, Carolina Lascelles and Ummey Hany from the Next Generation Sequencing Facility at the University of Leeds for RNA sequencing; to Rita Delgado Silva Marques from the Department of Biochemistry in Cambridge for PACE analysis. Thanks to Joe Bennett and Andrew Smith for their help with protein annotation.
This project (grant number BB/L001926/1) was funded by the UK Biotechnology and Biological Sciences Research Council. The York Centre of Excellence in Mass Spectrometry was created thanks to a major capital investment through Science City York, supported by Yorkshire Forward with funds from the Northern Way Initiative, and subsequent support from EPSRC (EP/K039660/1; EP/M028127/1).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Scanning electron microscopic images of the gills and food groove of L. pedicellatus. A) Close up of the gills lamellae to show the numerous cilia that capture food and draw it to the food groove. B) Close up of the food groove to show the numerous cilia (most of the mucus is lost during critical point drying). File format .DOCX.
Searching for the opening of the duct of Deshayes in the shipworm digestive system. Left, transverse sections through the shell valve, from the anterior to posterior adductor muscle (A-E). Right, 3D rendered model of the whole shipworm, with the green highlighted region showing the position of the transverse cross section displayed on the left. File format .DOCX.
Scanning electron microscopy of L. pedicellatus. A dissected specimen with the mantle removed, showing the food groove connecting the gills to the mouth. File format .DOCX.
Bacterial symbionts in L. pedicellatus. A) TEM image of a gill bacteriocyte. Arrows indicate the numerous rod-shaped bacteria in the cell. B) Close-up view of gill bacteria (TEM). C) Detailed TEM showing a cross section of a gill bacterium and some of its features (periplasmic space and membranes). File format .DOCX.
SDS-PAGE and nitrocellulose western blot obtained by loading 0.6 μg of the recombinant purified bacterial LpsGH5_8 (without the appended CBM). A) Coomassie Brilliant Blue stained SDS-PAGE gel of the purified protein. B) Western blot detection performed with the purified pre-immune serum. C) Western blot detection performed with the purified anti-protein serum. File format .DOCX.
Immunogold labelling of LpsGH5_8 in L. pedicellatus using pre-immune serum (negative control). A) SEM image illustrating the anatomical position of the gills, food groove, caecum and mouth in the shipworm body. B-D) Immunogold labelling of the lumen of the food groove (B), lumen of the caecum (C), and the gills (D) performed with pre-immune serum. No gold particles are observed in any of the samples. B = bacteria. File format .DOCX.
Details of the proteins identified through proteomic analysis of the crystalline style. The annotation was performed against the NCBI non-redundant database, the signal peptides were identified with the Signal P4.1 server and the CAZy domains were searched using dbCAN. Many of the contigs were out of frame or not full length (particularly the bacterial ones) and therefore the identification of the signal peptide was not possible. File format .DOCX.
CAZy families. List of the CAZy families mentioned in the results, with the description of the activities that have been recorded for enzymes listed in each family. The information has been gathered from the CAZy database (http://www.cazy.org). File format .DOCX.
CAZy families of the caecum. Table showing the relative abundance (calculated from the emPAI score) of eukaryotic and prokaryotic CAZy families identified in the proteomic analysis of the caecum content of L. pedicellatus. File format .DOCX.
Details of the proteins for which heterologous recombinant expression from symbionts cDNA was successful. The annotation was performed against the NCBI non-redundant database, the signal peptides were identified with the Signal P4.1 server and the CAZy domains were searched using dbCAN. Dis. bonds = disulfide bonds. File format .DOCX.
The amino acid sequences of the bacterial proteins LpsGH5_8, LpsGH11, LpsGH134a, LpsGH134b and LpsAA10A. Text in red represents the signal peptide. File format .DOCX.
Characterisation of the recombinant bacterial CAZymes encoded by the endosymbionts. DNS reducing sugars assays showing activities on a number of substrates for LpsGH5_8, LpsGH11, LpsGH134a and LpsGH134b. CMC = carboxymethyl cellulose, LBG = locust bean gum. The nanomoles of sugars released by the different proteins cannot be compared quantitatively as different amounts were used for the assay. File format .DOCX.
MALDI-TOF MS analysis of in vitro activity assays with purified LpsAA10A under the same experimental conditions as in Fig. 5 C. Panels from a to d show spectra of products obtained after incubation of the enzyme with 4 mg mL-1 Avicel (a), 4 mg mL-1 Avicel plus 4 mM gallic acid (b), 4 mg mL-1 Avicel plus 2 μM LPMO (c) and 4 mg mL-1 Avicel plus 2 μM LPMO and 4 mM gallic acid (d). In panels a to d, 100% relative intensity represents 1.3 × 104 arbitrary units (a.u.). The panels from e to g show spectra of products obtained after incubation of 4 mg mL-1 PASC (e), 4 mg mL-1 PASC plus 4 mM gallic acid (f) and 4 mg mL-1 PASC plus 2 μM LPMO (g). In panels e to g, 100% relative intensity represents 1.0 × 104 arbitrary units. File format .DOCX.
Primers used for the cloning of the bacterial proteins LpsGH5_8, LpsGH11, LpsGH134a and LpsGH134b. LpsAA10A could not be amplified from the cDNA and therefore a synthetic version of the gene was codon-optimised for E-coli expression. File format .DOCX.
About this article
Cite this article
Pesante, G., Sabbadin, F., Elias, L. et al. Characterisation of the enzyme transport path between shipworms and their bacterial symbionts. BMC Biol 19, 233 (2021). https://doi.org/10.1186/s12915-021-01162-6
- Lyrodus pedicellatus
- Crystalline style
- Food groove