Skip to main content

Spotting disease disrupts the microbiome of infected purple sea urchins, Strongylocentrotus purpuratus

Abstract

Background

Spotting disease infects a variety of sea urchin species across many different marine locations. The disease is characterized by discrete lesions on the body surface composed of discolored necrotic tissue that cause the loss of all surface appendages within the lesioned area. A similar, but separate disease of sea urchins called bald sea urchin disease (BSUD) has overlapping symptoms with spotting disease, resulting in confusions in distinguishing the two diseases. Previous studies have focus on identifying the underlying causative agent of spotting disease, which has resulted in the identification of a wide array of pathogenic bacteria that vary based on location and sea urchin species. Our aim was to investigate the spotting disease infection by characterizing the microbiomes of the animal surface and various tissues.

Results

We collected samples of the global body surface, the lesion surface, lesioned and non-lesioned body wall, and coelomic fluid, in addition to samples from healthy sea urchins. 16S rRNA gene was amplified and sequenced from the genomic DNA. Results show that the lesions are composed mainly of Cyclobacteriaceae, Cryomorphaceae, and a few other taxa, and that the microbial composition of lesions is the same for all infected sea urchins. Spotting disease also alters the microbial composition of the non-lesioned body wall and coelomic fluid of infected sea urchins. In our closed aquarium systems, sea urchins contracted spotting disease and BSUD separately and therefore direct comparisons could be made between the microbiomes from diseased and healthy sea urchins.

Conclusion

Results show that spotting disease and BSUD are separate diseases with distinct symptoms and distinct microbial compositions.

Graphical abstract

Peer Review reports

Background

Spotting disease, or red spotting disease, has been described for many sea urchin species and shows a range of symptoms based on both the species that are impacted and the severity of the disease [1,2,3,4,5]. Spotting disease is characterized by discrete viscous surface lesions composed of blackened necrotic tissue on the body wall that typically occur at the ‘equator’ of these spheroid animals and that result in a loss of all surface appendages within the lesion. Spotting disease requires injury to initiate the infection [6, 7], which can happen when the spines of neighboring sea urchins abrade or pierce the surface tissues. This is consistent with the location of lesions on the equator of the body. Lesions range in size from a few mm2 to more than a third of the body surface [8] and the infected tissue can be discolored, which includes green [8], blue-green/olive [9], dark red/purple [1], or black [2, 4, 5]. The lesions expand to infect more of the body surface and deepen to expose and subsequently degrade the test. Spotting disease lesions have been reproduced in a laboratory setting by mechanical abrasion, which has also been used to identify causative agents of spotting disease by infecting wounds with specific bacteria following the abrasion [2, 4, 6, 7]. Many bacterial species reproduce the discrete lesions and therefore there are multiple causative agents underlying spotting disease [2, 4, 5, 7]. Spotting disease is a major problem for sea urchins housed in aquaculture facilities because of the high density of the animals, as well as frequent handling and substandard culturing conditions that increase the frequency of disease [10], which results in significant economic loss [11]. Separation of sea urchins into individual spaces in aquaculture facilities significantly reduces the impact from diseases [12].

A separate disease of sea urchins, bald sea urchin disease (BSUD), also affects sea urchin species across many locations. BSUD has variable descriptions in the literature, but is generally characterized by a loss of appendages over most or all of the surface [13,14,15]. However, surface lesions may be present that degrade the epidermal tissue and can damage the muscle fibers at the base of the spines, in the tube feet and the associated ampullae [6, 8, 13, 14, 16,17,18,19,20,21,22]. Because BSUD is commonly described with necrotic lesions, which are characteristic of spotting disease, it is generally lethal because the lesions may perforate the test plate causing death [1, 8, 17, 18]. As a result, the descriptions of BSUD and spotting disease overlap, which results in significant confusion regarding the differences between these two diseases. Because we have reported changes in the microbiome on sea urchins infected with BSUD in our aquaria [15], and spotting disease occurred on animals in the same aquaria, our aim is to understand whether the microbiomes associated with these two diseases are different when the animals are housed under the same conditions.

Here, we describe the lesions and the microbiomes associated with spotting disease on purple sea urchins, Strongylocentrotus purpuratus, housed in closed, recirculating aquaria. Although previous studies have focused on identifying the causative agent underlying spotting disease, which suggests a wide array of pathogenic bacteria, few studies have reported an analysis of microbiomes of tissues from sea urchins infected with spotting disease [5]. Furthermore, none report the microbiome of tissues from sea urchins housed in a closed aquarium system, which have distinct properties and dynamics compared to open aquarium systems, leading to direct impacts on the host microbiome [23]. Because the microbiome on sea urchins shifts when transferred from the ocean to a closed aquarium [24], sampling sea urchins housed in the same aquarium system allows direct comparisons among the microbiome samples. We collected global surface samples, lesion surfaces, lesioned and non-lesioned body wall tissues, and coelomic fluid from diseased and healthy sea urchins, and the 16S rRNA gene was sequenced to characterize the microbiomes of each sample. Results indicate that the global surface microbiomes consisting of the microbes present on the external surface differ between diseased and healthy sea urchins, despite being housed in the same aquarium. Furthermore, the microbiomes differ between lesioned and non-lesioned body wall tissue sampled from diseased sea urchins, and between non-lesioned body wall tissue and coelomic fluid from diseased compared to healthy sea urchins. Results also show that the microbiomes associated with spotting disease are distinct from those associated with BSUD for sea urchins housed in the same aquarium. This suggests that spotting disease and BSUD are different diseases and that this difference correlates with discrete lesions that are characteristic of spotting disease but not BSUD.

Materials and methods

Sea urchin husbandry

Purple sea urchins, Strongylocentrotus purpuratus, were purchased from the Southern California Sea Urchin Company (Corona del Mar, CA) and shipped to George Washington University in Washington DC. Sea urchins were housed in a 125 gallon aquarium (aquarium B) with recirculating artificial seawater (Premium Marine Salt, OmegaSea), salinity of 32–35 ppt, 13–14℃, and outfitted with both physical and bio-filters, a UV light housing, and a protein skimmer. The central aquarium pump (Pond-Mag 9.5, Pondmaster) that was positioned in the aquarium sump, circulated 950 gallons/hour through the system. The water quality was maintained with weekly seawater changes of 5 gallons that also served in solid waste removal. All animals were fed weekly with rehydrated brown seaweed, Saccharina angustata (Kjellman) (WEL-PAC). Diseased sea urchins were maintained in individual plastic floating holding boxes to prevent interaction with healthy sea urchins and to monitor disease progression.

Treatment of sea urchins with penicillin and streptomycin

Based on our standard protocol, all shipments of sea urchins were treated upon arrival by immersion for 1–2 h at 14℃ in a tray (8 L) of freshly prepared artificial seawater with 12 mg/L penicillin and 50 mg/L streptomycin sulfate (pen/strep). After treatment, sea urchins were returned to their aquarium.

Sample collection

Surface microbiome

To collect samples from the global surface microbiomes, diseased sea urchins (n = 4; D1-D4, Table 1) and healthy sea urchins (n = 4; H1-4) from aquarium B were placed in a funnel connected to a microbial collection system as described previously [15]. Briefly, seawater (500 ml) from aquarium B was poured onto the ventral side of the sea urchins, so that the seawater washed over the entire body and the cellular material was collected on nylon filters (0.22 𝜇m, 47 mm diameter; GVS Filter Technology) held on a filter-holder assembly. Seawater samples (500 ml, n = 2) (fSW) were filtered in the same manner and served as controls. A sample of 500 ml of freshly mixed OmegaSea seawater (foSW) was also filtered to serve as the negative control. The nylon filters were stored in 50 ml falcon tubes at -80℃ for later use.

Table 1 Definitions of sample abbreviations

Swab collections

Diseased sea urchins (n = 4) were rinsed with 500 ml 0.22 𝜇m filtered artificial seawater (faSW) that was poured directly onto the ventral side of the animals such that the seawater washed over the body surface and removed microbes that were not tightly associated with the lesions. Sterile swabs were used to collect cells from the lesion surface (LS) with gentle rotation and twisting until the swab tip was coated with material. The coated swab was inserted into a 50 ml falcon tube filled with 35 ml faSW, and the swab stick end was removed using sterile scissors. A sterile swab was dipped into the seawater in aquarium B and inserted into a 50 ml falcon tube to collect sample microbes in the aquarium seawater (sSW). The falcon tubes containing swabs were vortexed vigorously to release the microbes, and the contents were poured into the vacuum filtration system to transfer the cells to the nylon 0.22 𝜇m filters, as described above. An additional 20 ml of faSW was added to the falcon tube with the swabs, vortexed vigorously, and poured into the vacuum filtration system. Filters were stored in 50 ml falcon tubes at -80℃.

Sea urchin sacrifice for tissue collection

Surviving diseased sea urchins (n = 3) and healthy sea urchins (n = 3) were rinsed with 500 ml faSW, as described above. Sea urchins were sacrificed by cutting the peristomial membrane around Aristotle’s lantern with sterile scissors and removing the lantern. The whole coelomic fluid (CF, cells and fluid) from each animal was poured into a 50 ml falcon tube (diseased CF, DCF; healthy CF, HCF). For diseased animals, tissue samples of lesioned body wall (LBW, n = 7) and non-lesioned body wall (DBW, n = 6) were collected by cutting through the body wall with sterile scissors and sterile forceps. Non-lesioned areas of the body wall (HBW, n = 6) were also collected from healthy animals (H1-H3) by the same method. Replicate samples were collected for each tissue except for D3LBW (the entire lesion was collected in the first sample), and the CF samples (all of the coelomic fluid was collected for each sea urchin). Tissue samples were placed in 50 ml falcon tubes filled with 35 ml of ice chilled faSW and ground with glass rods and vortexed vigorously, whereas the coelomic fluid was allowed to clot on ice. Tissue samples were centrifuged at 4300 x g for 5 min at 4℃, and the supernatant was poured into the vacuum system to collect the microbes on 0.22 𝜇m filters, as described above. Additional faSW (40 ml) was added to each pelleted tissue sample, vortexed vigorously, centrifuged, and the supernatant collected on filters. This process was repeated once. Filters were stored in 50 ml falcon tubes at -80℃.

Genomic DNA isolation from samples

The genomic (g)DNA isolation from nylon filters was carried out according to Turner et al. [25] with modifications described by Shaw et al. [15]. Briefly, filters were incubated with 1 ml cetyltrimethylammonium bromide (CTAB; 2% CTAB, 100 mM Tris pH 7.4, 4 M NaCl, 1% polyvinylpyrrolidone, 20 mM EDTA), to which was added 1 ml of chloroform:isoamyl alcohol (24:1), precipitated with 2.5 M NaCl and 50% isopropanol, and resuspended in Tris-EDTA buffer (TE; 10 mM Tris base pH 7.4, 1 mM EDTA). The gDNA concentration was evaluated on a spectrophotometer (NanoDrop 2000c, ThermoFisher). The gDNA size and level of degradation was evaluated with a 0.75% agarose gel with Tris-acetate-EDTA buffer (TAE; 40 mM Tris; 20 mM acetic acid, 1 mM EDTA) plus ethidium bromide and imaged with a UV imaging system (Kodak Molecular Imaging, Kodak Gel Logic 1500 Imaging System).

Polymerase chain reactions

To verify that bacterial gDNA was present in the samples prior to sequencing, the 16S rRNA gene was amplified with one of two sets of test primers; either 27 F (AGA GTT TGA TCC TGG CTC AG) and 1492R (ACG GTT ACC TTG TTA CGA CTT) [26], or 331 F (TCC TAC GGG AGG CAG CAG T) and 797R (GGA CTA CCA GGG TAT CTA ATC CTG TT) [27]. The PCR reaction of 20 μl contained 1X ExTaq Buffer, 200 μM dNTPs, 0.5 or 1 μM each primer, and 0.5 units Ex Taq DNA polymerase (Takara). The two sets of primers were used in different PCR programs. The step-down program for the 27 F/1492R primers was 94℃ for 30 s, followed by 4 cycles of 98℃ for 10 s, 57.4℃ for 30 s and 72℃ for 90 s, then 4 cycles of 98℃ for 10 s, 55.5℃ for 30 s and 72℃ for 90 s, followed by 17 cycles of 98℃ for 10 s, 53.7℃ for 30 s and 72℃ for 90 s, with a final extension of 72℃ for 2 min and a 4℃ hold. The step-down program for the 313 F/797R primers was 94℃ for 30 s followed by 4 cycles of 98℃ for 10 s, 58℃ for 30 s, 72℃ for 30 s, with a -1℃ change in annealing temperature per cycle to 55oC, then 21 cycles of 98℃ for 30 s, 54℃ for 30 s and 72℃ for 30 s, finally 72℃ for 2 min, and a 4℃ hold. The amplicons were analyzed on a 0.8% or 1% agarose gel with TAE buffer plus ethidium bromide and imaged on an UV system as described above.

16S rRNA amplicon sequencing

The gDNA samples were processed and sequenced using the ZymoBIOMICS targeted sequencing service at ZymoResearch (Irvine CA). The targeted sequencing of the bacterial 16S rRNA gene was carried out as described [15], using the Quick-16S NGS Library Prep Kit (ZymoResearch) with custom-designed primers to amplify the V3-V4 region of the 16S rRNA gene, real-time PCR and qPCR fluorescence readings, the Select-a-Size DNA Clean & Concentrator (ZymoResearch), Tapestation (Agilent Technologies) and Qubit (ThermoFisher Scientific). The positive control sample used for library preparation was the ZymoBIOMICS Microbial Community DNA Standard (ZymoResearch). The negative controls included the foSW sample as a control for the DNA extraction method, and a blank for library preparation that was provided by ZymoResearch. The completed library was sequenced using a V3 reagent kit (600 cycles) on Illumina MiSeq, which was calibrated by a 10% spike-in of PhiX DNA. The raw sequence reads were uploaded to the Sequence Read Archive database at NCBI under the BioProject ID PRJNA937707.

Amplicon sequence analysis

Amplicon sequence analysis was carried out as reported [15] using the DADA2 pipeline [28], SILVA release 138.1 for taxonomy [29] and the Phyloseq package (ver 1.42.0) [30]. Beta diversity was analyzed using weighted unique fraction (UniFrac) and visualized using non-metric multidimensional scaling (NMDS). Statistical analysis and Linear Discriminant Analysis Effect Size (LEfSe) were also carried out as previously reported, by linear discriminant analysis (LDA) cutoff set to 2, and all analyses were performed using R version 4.1.1 [31]. The R code with the complete pipeline can be found in the GitHub repository Spotting_disease_S_purpuratus.

Results

Spotting disease progression is similar for all infected sea urchins

While being housed in our aquaria, a small number of sea urchins contracted spotting disease, which allowed us to observe disease progression in a controlled environment. Diseased sea urchins showed one or two discrete regions of blackened necrotic tissue of varying sizes (Fig. 1). For all infected sea urchins, progression of the infection was similar. Initially the lesions were small and nearly undetectable; however, as the infection continued, the epidermal tissue was degraded, resulting in the loss of all appendages within the lesioned area (Fig. 1). When the lesions expanded and deepened, the underlying test was exposed, and sea urchins lost the primary spines in areas surrounding the lesion. Eventually, the test disintegrated, which signaled that a sea urchin would succumb soon afterwards, indicating that spotting disease was fatal. The duration of infection typically lasted months to years, however, the final phase of infection after the test was degraded tended to develop rather quickly as a sea urchin became moribund. In this late stage of the disease, sea urchins demonstrated altered behaviors, which, in addition to primary spine loss in non-lesioned regions, included cessation of eating and failing to hold on to anything with their tube feet (i.e., aquarium wall, kelp, holding box). In all cases, infected sea urchins never showed signs of recovery from spotting disease. Spotting disease was not communicable, as illustrated when sea urchins with lesions were housed in the same aquarium with healthy sea urchins, and the disease never appeared on the healthy animals.

Fig. 1
figure 1

Sea urchins with spotting disease show discrete lesions at the equatorial to ventral body regions. All sea urchins are positioned with the oral surface facing down. A Diseased sea urchin (D1) has a single large black necrotic lesion. This sea urchin displays unusual orientation of its spines, which point in various directions rather than uniformly perpendicular to the body surface. This is a behavioral indication of disease that has been noted previously [15]. B Diseased sea urchin (D2) has two large black necrotic lesions, labelled “a” and “b”. Parts of the test (white) are exposed around the outer region of lesion a (arrows). This sea urchin has lost all primary spines, including non-lesioned areas of the body surface, which indicates its moribund condition. C Diseased sea urchin (D3) has one large black necrotic lesion. This sea urchin shows typical spine orientation of perpendicular to the body surface, which indicates better health despite the lesion. Exposed test is also evident within the lesion (arrow). D Diseased sea urchin (D4) has a small black necrotic lesion of approximately 1 cm in diameter. The mouth is located on the ventral side (yellow arrow). This sea urchin also shows indications of better health, including primary spines generally pointing perpendicular to the body surface

The four infected sea urchins were sampled at different stages of spotting disease as inferred from the different sizes and numbers of lesions (Fig. 1). Diseased animal 1 (D1) had one large lesion and at the time of sacrifice, it was not eating consistently, which was an indication of impending mortality (Fig. 1A). Diseased animal 2 (D2) had two lesions of different sizes that were referred to as ‘a’ and ‘b’ for the larger and smaller lesions, respectively (Fig. 1B). At the time of sacrifice, D2 had lost all primary spines, was no longer feeding, and did not display tube foot function to attach to the box in which it was housed, all of which were indicators of its moribund state. Diseased animal 3 (D3) had one large lesion but exhibited normal tube foot function, all spines were present except for those surrounding the lesion, and it was eating normally (Fig. 1C). Sea urchins D1, D2, and D3 were all treated twice with pen/strep; once at half concentration, and a second time at the standard concentration, one week apart. The pen/strep treatments had no effect on resolving the discrete lesions. Diseased animal 4 (D4) had one very small lesion that was first noticed once it was approximately 1 cm in diameter (Fig. 1D). D4 did not show any other signs of altered behaviors associated with spotting disease. Unfortunately, this animal died due to unrelated factors and therefore was not used in the analysis of dissected tissues that were collected from the other sea urchins upon sacrifice. However, D4 samples were collected for the analysis of the global surface microbiome and LS.

Sequence data identifies ASVs in all samples

Preliminary analyses by PCR of the microbial community gDNA isolated from the samples of both diseased and healthy sea urchins confirmed that bacterial gDNA was present in in all samples. 16S rRNA gene sequencing resulted in a total of 81,536 non-chimeric reads (Table 2).

Table 2 Read processing summary

The global surface microbiomes of diseased sea urchins show taxonomic differences compared to the global surface microbiomes of healthy sea urchins

The global surface microbiomes on echinoids are poorly studied. However, an initial report shows that they are altered when variegated sea urchins, Lytechinus variegatus, are transferred from open water to a closed system [24]. The surface microbiomes also shift as sea urchins recover from BSUD, and it can differ on sea urchins housed in different aquaria [15]. Because there were healthy sea urchins in the same aquarium with spotting disease, and the disease was not communicable, we investigated the global surface microbiome to identify differences between diseased and healthy sea urchins. The alpha diversity of the global surface microbiomes on diseased and healthy sea urchins were analyzed using observed species, Chao1 [32], and ACE [33] (Fig. 2A-C). Results showed that there were no significant differences in alpha diversity between the global surface microbiomes on the diseased compared to the healthy sea urchins (ANOVA, p > 0.05), although the diseased group microbiome had decreased alpha diversity compared to the healthy group for all metrics. Beta diversity was analyzed by weighted UniFrac that measures microbiome composition, which showed overlap of the samples from diseased and healthy sea urchins (PERMANOVA, p > 0.05; Fig. 2D), indicating that the microbiome composition was not significantly different between the two groups.

Fig. 2
figure 2

The global surface microbiome diversity is not different between diseased and healthy sea urchins. Alpha diversity is analyzed by A Observed Species, B Chao1, and C ACE. The box plots show the mean and quartile values for each group, which are not significantly different (ANOVA, p > 0.05). D Beta diversity is analyzed at the level of ASV sequences using weighted UniFrac and visualized with NMDS. Ellipses around sample groups show 95% confidence intervals assuming a multivariate t-distribution (solid line) or a multivariate normal distribution (dashed line). The beta diversity of the groups is not significantly different (PERMANOVA, p > 0.05), indicating that microbial composition is similar for the global surface microbiomes on diseased and healthy sea urchins housed in the same aquarium

The composition of bacterial taxa collected from the two groups of sea urchins were compared to identify differences among the global surface microbiomes, in addition to comparisons with the two samples of filtered aquarium seawater (fSW). The most abundant taxa were selected (Additional File; Tables S1, S2) to identify the similarities and differences. Although beta diversity showed no differences among the samples (Fig. 2D), the phyla present in the global surface microbiomes were different between the diseased and healthy sea urchins (Fig. 3). Samples were generally dominated by Proteobacteria, and three of the diseased surface microbiome samples had an elevated abundance of the phylum Bacteroidota compared to the healthy surface microbiome samples (Fig. 3A). There were also major differences in the genera identified in the global surface microbiomes among the groups (Fig. 3B). The microbial composition of the diseased surface microbiome samples differed from the healthy surface microbiome samples, and there were also differences among the samples within the diseased group. The samples collected from D1 and D2 had many taxa in common with elevated abundances of a genus in the Cryomorphaceae family, a genus in the Cyclobacteriaceae family, Lutibacter, and a genus in the Cellvibionaceae family. This was not evident in the global surface samples from D3 and D4 nor from the corresponding samples from the healthy sea urchins. The D3 and D4 surface samples were more similar to the healthy surface samples, and were mainly composed of Psychromonas, Colwellia, and a family in the Bacteroidia class. Samples from the fSW were distinct from all samples collected from the sea urchins. It is notable that sea urchins D3 and D4, with less extensive lesions, had global surface microbiomes that were more similar to those on the healthy sea urchins, whereas D1 and D2, with more extensive lesions, had microbiomes that were very different from the healthy sea urchins. LEfSe analysis showed that Cyclobateriaceae, Lutibacter and Pseudoteredinibacter were significantly differentially abundant in the diseased global surface microbiome samples, and that Gastranaerophilales and Enterobacterales were significantly differentially abundant in the corresponding samples from the healthy sea urchins (Fig. 4A; Additional File, Table S3). These results showed that the microbial composition of the global surface microbiome differed in samples from sea urchins that were in different stages of spotting disease progression. Results also showed that taxonomic differences were evident for the genera identified in the global surface microbiomes of the diseased compared to the healthy groups of the sea urchins that were housed in the same aquarium, and that they were also different from the microbes in the fSW samples. These results indicated that spotting disease could be characterized based on the differences in the global surface microbiomes of diseased compared to healthy sea urchins irrespective of the influences of the bacterial composition of the seawater in the aquarium and did not require direct sampling of the lesions.

Fig. 3
figure 3

Global surface microbiomes of diseased vs. healthy sea urchins show taxonomic differences. A All identified phyla for the sea urchin surface samples and the fSW sample (Additional File, Table S1) illustrate the relative abundance of each taxon in each sample. B Genera with an average relative abundance of > 0.1% across all samples (Additional File, Table S2) are illustrated by the relative abundance per sample. fSW is the average of fSW1 and fSW2. Taxa in A and B that could not be assigned at the level of phylum or genus are listed as the most specific known taxonomic level. BD2-3 is in the order Victivallales, the Pir4 lineage is in the family Pirellulaceae, vadinHA49 is in the phylum Planctomycetota, JGI-0000069-P22 is in the class Gracilibacteria, MSBL3 is in the family Kiritimatiellaceae, and HOC36 is in the class Gammaproteobacteria. ASV sequences that could not be assigned to a phylum are grouped as Bacteria. Sample name abbreviations are defined in Table 1

Fig. 4
figure 4

Many taxa are differentially abundant in the microbiome samples. A heatmap shows all taxa identified by LEfSe that are significantly differentially abundant (p < 0.05) and have an LDA score of > 2 as their abundance per group for A the global surface microbiome samples, B the lesioned body wall (LBW) and lesion surface (LS) microbiome samples, and C the tissue microbiome samples (Additional File, Tables S3, S6, S9). Sample name abbreviations are defined in Table 1

The lesions have a unique microbial composition

One of the reported causative agents of spotting disease, Vibrio coralliilyticus, forms biofilms on the surface of tissues [4]. During the formation of a biofilm, bacteria aggregate, function as a single unit, and produce an extracellular matrix that protects the bacterial community from external factors such as antimicrobial agents [34]. Based on the appearance of the lesions on infected sea urchins and based on the lack of response to the pen/strep treatment by the diseased sea urchins, microbes associated with lesions may have formed a biofilm. To investigate this possibility, the LS was collected by the swab method to determine whether the microbial composition on the surface of the lesions differed from the microbes that were associated with the lesioned body wall tissues, as identified in the LBW group (Fig. 5). Results from Observed Species, Chao1, and ACE did not show significant differences between the microbiomes collected from LS compared to the LBW (Fig. 5A-C). Beta diversity showed that the LS samples clustered closely together compared to the LBW samples, which were less similar and more spread out with a large confidence interval (Fig. 5D). This indicated that the microbial composition among the LS samples was similar, and that there was more variation among the LBW samples. Furthermore, because the two clusters overlapped, the microbial compositions between the LS group and the LBW group likely had many shared taxa. The microbial composition in the sSW was very different from the samples collected from both the LS and the LBW. These results suggested that microbiome compositions were similar among the LS samples and less similar among the LBW samples, and were both different from the microbes in the aquarium seawater.

Fig. 5
figure 5

The microbiome compositions are similar between the lesion surface and the lesioned body wall. Alpha diversity of the microbiomes from the LS and LBW sample groups are evaluated by A Observed Species, B Chao1, and C ACE. The box plots show the mean and quartile values for each group, which are not different (ANOVA, p > 0.05). D Beta diversity is evaluated at the ASV level using weighted UniFrac and visualized with NMDS. Ellipses around sample groups show 95% confidence intervals assuming a multivariate t-distribution (solid line) and a multivariate normal distribution (dashed line). The microbiomes of the LS compared to the LBW sample groups are not significantly different (PERMANOVA, p > 0.05). The sSW sample is shown for comparison. Sample name abbreviations are defined in Table 1

Further analysis of the microbiomes of the LS, LBW, and seawater samples (fSW and sSW) was carried out to compare the taxa identified in those sample groups (Fig. 6). At the level of phylum, the microbial composition of the LS and LBW microbiome samples were similar and were composed of mainly Proteobacteria and Bacteroidota, which varied inconsistently among the samples (Fig. 6A; Additional File, Table S4). At the level of genus, there were also many similarities in the microbial compositions (Fig. 6B; Additional File, Table S5), which were consistent with the beta diversity results (Fig. 5D). The LS samples from D1-D3 were highly similar to each other and were mainly dominated by a genus in the Cyclobacteriaceae family, a genus in the Cryomorphaceae family, and a genus in the Cellvibrionaceae family. The LS sample from D4, which had a single small lesion, differed from the others and was mainly composed of a genus in the Cellvibrionaceae family, the HOC36 group, and Pseudophaeobacter. The LBW samples also had an elevated abundance of a genus in the Cyclobacteriaceae family, and a genus in the Cryomorphaceae family, however, two of the LBW samples (D2a and D3) had elevated abundances of Vibrio and Candidatus Photodesmus. LEfSe results identified Cellvibrionaceae, Pseudophaeobacter, the HOC36 group, Gammaproteobacteria, and Roseobacter as significantly differentially abundant in the LS microbiome samples, whereas Canditatus Photodesmus was significantly differentially abundant in the LBW microbiome samples (Fig. 4B; Additional File, Table S6). Overall, all LBW samples were composed of a similar set of microbes.

Fig. 6
figure 6

The microbiomes of the lesion surface and the lesioned body wall are highly similar. A All identified phyla are shown as the relative abundance in each sample (Additional File, Table S4). B Genera with an average relative abundance of > 0.1% across all groups (Additional File, Table S5) are illustrated by the relative abundance per sample. Abundance of taxa from replicated samples are averaged. Both types of seawater control samples are included for comparisons, which are the microbes collected from 500 ml of filtered seawater (fSW) and seawater collected with a swab (sSW). Taxa in A and B that could not be assigned at the level of phylum or genus are listed as the most specific known taxonomic level. Taxa that could not be assigned to a phylum are grouped under Bacteria. Sample name abbreviations are defined in Table 1. BD2-3 is in the order Victivallales, the Pir4 lineage is in the family Pirellulaceae, vadinHA49 is in the phylum Planctomycetota, JGI-0000069-P22 is in the class Gracilibacteria, and HOC36 is in the class Gammaproteobacteria

When the LS microbiome samples from the diseased sea urchins were compared to each other, as well as to the sSW sample, D1-D3 had similarly abundant taxa, however, the LS microbiome from D4 differed from the other LS microbiomes (Fig. 6). Although an initial aim was to compare the microbiomes of the LS samples with non-lesioned surface areas of sea urchins with spotting disease, attempts to collect microbial gDNA from the surfaces of non-lesioned tissue using the swab method was unsuccessful. The 16S rRNA gene could not be amplified from the DNA isolated from the non-lesion surface swab samples indicating that were likely too few microbes associated with surfaces of healthy sea urchins, in agreement with a previous report [7]. Although another, more destructive collection method has been reported with more success for isolating bacterial gDNA from sea urchin surfaces and spines [24], this approach would not have been comparable to the LS sampling that was collected with swabs. Consequently, comparisons to non-lesioned surface microbiomes collected by the swab method could not be carried out.

Comparisons of the abundant taxa in the microbial compositions of the seawater samples showed that sSW sample included a genus in the Cellvibrionaceae family, Microbacterium, and Psychromonas. The abundant taxa in the fSW samples consisted mainly of the HOC36 group, the Pir4 lineage, and Psychromonas. Differences in the microbial composition of these seawater samples collected from the same aquarium may be attributed to differences in the volume of seawater collected. The small water volumes collected by swabs for the sSW samples may not have acquired an evenly distributed sample of microbes. Nonetheless, when the LS microbiomes were compared to the seawater samples, there were many taxa in the LS microbiomes that were absent from the microbiomes collected from either of the aquarium seawater samples. Overall, results indicated that the LS microbiome was distinct from the microbes in the seawater, but that it was highly similar to the LBW microbiome, which suggested similar a microbial composition on the surface of the lesion and within the lesioned body wall.

Microbiomes are different between diseased and healthy sea urchin tissues

A comparison of the microbiomes of dissected tissues from diseased compared to healthy sea urchins housed in the same closed aquarium system eliminates the variability of the microbial environment and is therefore informative for characterizing the microbiomes associated with the spotting disease infection on tissues that are outside of the lesioned areas. The alpha diversity metrics of the microbiomes from dissected tissues from infected sea urchins were evaluated and compared to each other and to the microbiomes of dissected tissues from healthy sea urchins (Fig. 7). Results from Observed Species, Chao1, and ACE, identified large variations in the samples within the groups of sea urchin tissues and therefore did not show significant differences among the groups (ANOVA, p > 0.05) (Fig. 7A-C). Beta diversity showed that the sample groups largely overlapped, but that the LBW samples tended to cluster together, and the HBW samples also clustered (Fig. 7D). The DBW samples had large variations, and therefore had a large confidence interval. The DCF samples tended to cluster with the LBW samples, whereas the HCF samples clustered with the HBW samples, indicating similar microbiome compositions for these pairs of sample groups. The beta diversity analysis revealed differing microbial compositions that suggested unique microbiomes for the LBW, DBW and HBW tissues.

Fig. 7
figure 7

The microbiomes of dissected tissues from diseased and healthy sea urchins are different. Alpha diversity is analyzed by A Observed Species, B Chao1, and C ACE. The box plots show the mean and quartile values for each group, which are not significantly different (ANOVA, p > 0.05). D Beta diversity is analyzed at the ASV level using weighted UniFrac. Ellipses around sample groups show 95% confidence intervals assuming a multivariate t-distribution (solid line) and a multivariate normal distribution (dashed line). Groups are not significantly different (PERMANOVA, p > 0.05). Sample name abbreviations are defined in Table 1

To understand the taxa underlying the bacterial compositions of the dissected tissues, the phyla identified in the microbiomes were compared among the samples (Fig. 8). Differences in the microbial composition were evident among the different tissues based on the relative abundance of the phyla (Fig. 8A; Additional File, Table S7). All groups were dominated by Proteobacteria, however, the LBW group including D1, D2a, b, and D3 (Fig. 1), was also largely composed of Bacteroidia, which was generally only present in low abundances in the other samples. The DBW samples from D1, D2, and D3 were largely composed of Proteobacteria with a few other phyla of low abundance, which was similar to the DCF samples from the same animals. The HBW group and HCF group from the healthy sea urchins similarly had low abundances of Bacteroidia, as well as Actinobacteriota and Verrucomicrobia that were not detected in the LBW samples. The taxa were also compared at the level of genus to identify differences among the sample groups (Fig. 8B; Additional File, Table S8). All LBW samples had an elevated abundance of a genus of the Cryomorphaceae family and a genus of the Cyclobacteriaceae family. In addition, there was an elevated abundance of Candidatus Photodesmus and Vibrio in two of the four LBW samples. The DBW samples differed in their microbial compositions, which was consistent with the large variation shown by beta diversity. The DCF samples, similar to the DBW samples, had differing microbial compositions for different sea urchins. The HBW samples consisted of a few taxa in common, which included Psychromonas, a genus of the Cellvibrionaceae family, Colwellia, and Microbacterium, among others. The microbial composition of the HCF samples differed based on the sea urchin, however, there were many taxa that were also present in the HBW samples, which included BD2-3, Microbacterium, Psychromonas, and a family of the Bacteroidia class.

Fig. 8
figure 8

The microbiome of each dissected tissue has a distinct microbial composition. A All phyla that were identified are shown as the relative abundance of each phylum in each sample (Additional File, Table S7). Taxa that could not be assigned to a phylum are grouped under Bacteria. B Genera with an average relative abundance of > 0.1% across all groups are shown as their relative abundance per sample (Additional File, Table S8). Taxa that could not be assigned at the level of genus are listed as the most specific known taxonomic level. Sample name abbreviations are defined in Table 1. BD2-3 is in the order Victivallales, vadinHA49 is in the phylum Planctomycetota, JGI-0000069-P22 is in the class Gracilibacteria, and MSBL3 is in the family Kiritimatiellaceae

Significant differential abundance of taxa for individual sample groups was identified by LEfSe analysis (Fig. 4C; Additional File, Table S9). Results showed that these taxa in the LBW samples were Candidatus Photodesmus, Cyclobacteriaceae, Cryomorphaceae, Vibrio, Pseudoteredinibacter, and Lutibacter. Taxa identified for the DCF samples were Pseudoalteromonas and Pseudophaeobacter, and taxa identified for the HCF samples were Microbacterium, BD2-3, Bacteroidia, vadinHA49, and Brevibacterium. Overall, the microbiome of the LBW was composed of a few taxa which were generally absent from the samples of the other groups. The microbiomes in the tissues from diseased sea urchins were distinct from those from healthy sea urchins, as evidenced by the differences between the DBW compared to the HBW samples, and the DCF compared to the HCF samples.

Spotting disease and BSUD may have different etiologies

During this study and as reported previously [15], we observed differences between spotting disease and BSUD that both occurred in our aquaria. After receiving a shipment of sea urchins and placing them in aquarium B, all animals from that shipment contracted BSUD (see Fig. 1 in [15]). In comparison, there were four sea urchins (D1-4) that were housed in aquarium C that originated from a different shipment and became infected with spotting disease. BSUD was communicable based on the spread of the disease to all animals in aquarium B, including a number of sea urchins that were transferred into aquarium B during the outbreak [15]. D3, with a single discrete spotting disease lesion, was one of the sea urchins transferred into aquarium B with animals infected with BSUD. It began to show symptoms of BSUD about 60 days after the transfer including agitated spine movement, tissue swelling at the base of the spines, and primary spine loss, consistent with the other sea urchins with BSUD [15]. As a result, sea urchin D3 appeared to be infected with both BSUD and spotting disease simultaneously and the symptoms of each disease could be discerned separately. D3 did not succumb to BSUD but followed the same infection and recovery process as the other animals in the aquarium, including primary spine regrowth over the course of 2–3 months. However, the spotting disease lesion on D3 did not resolve. BSUD and spotting disease appeared as separate diseases based on their different symptoms and different recoveries, particularly for D3.

Discussion

The global surface microbiomes on diseased sea urchins are different from those on healthy sea urchins

The microbial composition of the global surface microbiomes for the healthy compared to the diseased sea urchins show differences. This is notable for sea urchins D1 and D2 that had the largest surface lesions, which is consistent with differences in microbial taxa between D1 and D2 compared to D3 and D4. For example, the lesion on D4 was much smaller than the lesions on the other infected sea urchins (Fig. 1), perhaps because it may have been more recently contracted, suggesting less impact of lesion-associated microbes on the global surface microbiome. The possibility of an early-stage infection on D4 may be a basis for the similarity of the bacterial composition with the surface microbiomes of the healthy sea urchins rather than the surface microbiomes for D1 and D2 (Fig. 3B). Furthermore, some of the key genera in the LS microbiomes for D1 and D2 are members of the Cryomorphaceae and Cyclobacteriacae families, which are also some of the key genera in their global surface microbiomes. This infers that the taxa present on the LS, which is the most proximal source of microbes to the entire animal surface, may influence or perhaps invade and become established in the microbial composition of the global surface microbiome of diseased sea urchins. There is a correlation between larger lesions and increased changes to the microbiome on the global surface of the animal. Taxonomic differences are evident between the diseased and healthy surface microbiomes despite factors that may reduce these differences, such as the microbes in the fSW sample and the varying stages of the lesions among the diseased sea urchins. These differences are especially noteworthy given that all sea urchins are housed in the same aquarium. These results indicate that the spotting disease infection alters the global surface microbiome of diseased sea urchins, which is different from the microbes in the seawater of the aquarium.

The microbiomes of the lesion surface and the lesioned body wall are highly similar

The LS microbiome and the LBW microbiome are similar (Figs. 5 and 6) indicating that the microbes collected from the lesion surface and those identified from the entirety of the lesion including the body wall are largely uniform. This suggests that there is no surface biofilm composed of different taxa or a different relative composition of taxa. Among the diseased sea urchins, LS microbiome samples have similar microbial compositions (Fig. 5D), which is noteworthy considering that the compositions of the LBW microbiome samples are less similar. However, the LS of D4 (Fig. 6) is composed of different taxa than the other samples, which may be the outcome of the small lesion size that resulted in sampling difficulty. The similarities in the microbial composition between the LS and LBW microbiomes indicates that the same taxa are present throughout the lesion, and that each lesion has the same microbial composition for each of the diseased sea urchins. The identification of similar taxa based on two different methods of sampling strengthens the conclusion that these taxa underlie the spotting disease lesions.

The microbiome of diseased sea urchin tissues differs from that of healthy sea urchin tissues

The microbiomes characterized from dissected tissues show that the samples from diseased sea urchins are different from samples from the healthy sea urchins. The HBW samples all have similar microbiomes, and therefore may be considered as an indicator of a standard microbial membership and composition for healthy sea urchins in our aquarium. However, it should be noted that a single “healthy” microbiome likely does not exist, as a major characteristic of a healthy microbiome is its ability to be dynamic [23]. Nonetheless, because the microbiomes of the DBW samples differ from those of the HBW samples (regions of non-lesioned tissue) from sea urchins in the same aquarium, this indicates that the spotting disease infection affects surface areas of the sea urchin that are outside of the observable lesions. This is also the case for the DCF for which the microbial composition differs from that of HCF. Although microbes have been shown to be present in the CF of sea urchins [35], another report shows that the coelomic fluid is sterile and lacks bacteria completely [36]. Differences may be based on methods of CF collection. For example, the body wall of sea urchins cannot be sterilized prior to CF collection by needle aspiration, therefore the presence of microbes in the coelomic fluid may be due to external contamination. Here, microbes are identified in the HCF samples, however, because of the sacrifice method, an unknown portion of the microbes in the CF samples are likely the result of contamination from other tissues. Sources may include the body wall, which must be cut through for collection, and/or from the esophagus - gut intersection that is broken when Aristotle’s lantern is removed prior to collecting the CF. The likelihood of contamination is supported by our results because the microbial composition of the HCF microbiome shares many similarities with the HBW microbiome samples (Fig. 8). Previous studies of the gut microbiome from S. purpuratus in various housing conditions have found that the typical composition is mainly Proteobacteria and Bacteroidetes, as well as Vibrio, Arcobacter, Sulfurimonas, Desulfotalea, Psychromonas, and Shewanella among others [37, 38]). Based on the microbial composition identified for the tissue microbiome samples reported here, there are many similarities between these taxa and the taxa reported previously for the gut microbiome in S. purpuratus [37, 38]. Because all sea urchins were housed in the same aquarium for years, this supports the notion that microbes from the gut that contaminated the CF samples would be similar for all dissected sea urchins. However, the microbial composition of the DCF is different from that of the HCF, suggesting that the spotting disease lesions introduce microbes into the coelomic fluid, which are different from gut contaminants, and that they are not cleared by coelomocytes or other attributes of the immune system. This is consistent with many shared microbes in the DCF samples with microbes in the LBW including Cryomorphaceae, Vibrio, and Candidatus Photodesmus. Consequently, the bacteria in the lesions that degrade and penetrate the test gain access to the coelomic fluid.

The taxa we identify that are present in the spotting disease lesions includes Vibrio that are consistently present for most or all samples of the LBW and are increased compared to healthy sea urchins, which is consistent with a previous report on spotting disease [5]. However, because culturing marine bacteria often yields only a small subset of the total microbes in environmental samples [39], we opted not to culture bacteria and therefore we did not perform experiments to confirm a causative agent based on Koch’s postulates. Nonetheless, the pathogenic bacteria that cause spotting disease may be normally associated with the tissues of sea urchins, and that they only become opportunistic pathogens upon colonizing wounds on the sea urchin surface. Based on these results and based on the wide array of bacteria that have been reported as causative agents of spotting disease, the disease may be the outcome of many different pathogenic bacteria, that may normally be associated with sea urchin tissues.

The microbiome composition differs among reports of spotting disease infections

The microbial characterization associated with spotting disease reported here differs from previous reports. The bacteria associated with spotting disease are either based on identification of culturable bacteria [2, 7] or are identified through analysis of 16S rRNA gene sequences ([5], this study), and sometimes both [4]. There are multiple reports of BSUD and the associated bacteria, but that describe discrete lesions that are characteristic of spotting disease [8, 13, 19, 40]. Consequently, reports that include descriptions of discrete lesions will be considered here as spotting disease infections irrespective of the disease named in the paper. A variety of pathogenic bacteria cause or are associated with the discrete lesions, which include Vibrio sp [5, 40], V. alginolyticus [19], V. coralliilyticus [4], Flexibacter sp [2], Acinetobacter sp [7], Tenacibaculum sp [41], Colwellia sp, Flexibacteraceae, Rhodobacterales, Stappia [40], Psychrobacter, Staphylococcus [5], Saprospiraceae, and Cohaesibacter gelatinilyticus [13]. In this study, Cryomorphaceae, Cyclobacteriaceae, Candidatus Photodesmus, and Vibrio are the major taxa in the microbiomes of the spotting disease lesions on sea urchins in our aquarium. These taxa belong to either Bacteroidota or Gammaproteobacteria, but otherwise, there is little in common amongst the various bacteria that underlie the spotting disease infection. Species in the Vibrio genus appear to be generally associated with sea urchin diseases, including spotting disease [4,5,6, 11, 42, 43]. However, despite the discrete lesions on the body wall of infected sea urchins that have similar appearances across different species of sea urchins and in different locations in the oceans, the bacteria associated with the lesions are relatively diverse, spanning many different orders. The majority of taxa identified in this study have not been identified as pathogens, and may be commensals to echinoids. Colwellia has been isolated from seawater as well as a variety of marine organisms, which includes sea urchins [13]. It has been characterized as an opportunistic pathogen of the seas urchin, Strongylocentrotus intermedius [44], and therefore may be acting as a pathogen in this study. Cyclobacteriaceae have been isolated from seawater and echinoids, likely acting as commensals [45], although they have also been associated with diseased corals [46]. Cyclobacteriaceae belong to the bacterial order Cytophagales that are known for their ability to degrade chitin and other biopolymers such as pectin, cellulose, and agar, and may be involved in carbon cycling or remineralization of organic carbon [47,48,49,50]. While it is possible that Cyclobacteriaceae are actively involved in spotting disease progression, it may be more likely that they act as decomposers to degrade complex carbon compounds into simple units in marine environments. However, Cytophagales (and Rhodobacterales) also proliferate in coral tissues infected with Vibrio coralliilyticus, suggesting that they are opportunistic and act as destabilizers of the microbiome in diseased corals [51]. Thus, while the Cytophagales probably are not responsible for initiating spotting disease in echinoids, they may be involved in its progression. Cryomorphaceae have been isolated from diseased echinoids, however they have not been shown to be pathogenic [40]. It is more likely that they are opportunists or secondary pathogens along with Rhodobacteraceae, as in the case for corals [52, 53]. Candidatus Photodesmus is not well documented, and its functions in marine systems are unclear. The possibility that this bacterium is a commensal, pathogen, or both cannot be ruled out. Although Vibrio is involved in a variety of marine diseases, it also acts as a commensal and is an integral part of the sea urchin microbiome [5, 13, 38]. The Vibrio splendidus clade is associated with a disease affecting the sea urchin, Paracentrotus lividus [21], and Vibrio coralliilyticus has been identified as the primary agent of red spotting disease in S. intermedius [4]. Consequently, it is likely that it has important functions in spotting disease reported here. In addition to the taxa described above, many of those identified in S. purpuratus have been isolated from other marine organisms, suggesting functions as commensals in the echinoid microbiome. The host-microbe interaction is complex and commensals may become pathogenic in response to environmental changes. While limited information exists regarding the function of these specific taxa within the microbiota of marine organisms, the integration of our findings with the current literature suggests that spotting disease may be the result of, or in part due to, an overgrowth of commensal bacteria following a stressor that effects the host immune system. Our results add to this wide array of bacteria that are associated with spotting disease, and suggest that whether the spotting disease infection occurs in a closed aquarium system, in aquaculture facilities [2, 4, 5], or in natural environments [8, 19, 40], the putative pathogens are diverse. The taxa that invade tissues are likely the taxa that are present in or on the sea urchin tissues at the time of injury, which is assumed to be required to initiate a lesion [6, 7]. Overall, the integration of our results with previous reports suggests that the characteristic spotting disease lesions are likely caused by many different bacterial species, forming what is known as the “pathobiome”, which is a subset of microbes that is associated with negative effects on host health [54].

The appearance and outcome of spotting disease is distinct from bald sea urchin disease

Spotting disease and BSUD are described in the literature interchangeably with a wide range of overlapping symptoms that lead to confusions in distinguishing the diseases [1, 11, 15]. Based on our observations of an outbreak of BSUD in aquarium B, and separate infections of spotting disease in aquaria B and C, as well as sea urchin D3 that was infected with both BSUD and spotting disease simultaneously, we conclude that spotting disease and BSUD are separate diseases based on different symptoms, infectivity, severity, and lethality. The spotting disease infection that we observed is consistent with other reports of the disease, based on the descriptions of discrete surface lesions [1, 5]. The BSUD infection that we report [15], however, is different from other reports of BSUD that include descriptions of discrete lesions [7, 8, 19], which our sea urchins do not show. Our observations suggest that whole-body surface infection vs. discrete lesions are defining characteristics of BSUD vs. spotting disease, respectively, because lesions are present in all cases of spotting disease but are not necessarily a symptom of BSUD. In contrast, an infection that encompasses the entire animal surface and causes partial to complete loss of surface appendages is the key symptom of BSUD, but not spotting disease. The confusions in the literature may be based on the possibility that sea urchins can be infected simultaneously with both diseases ([4, 21], this study). Furthermore, when sea urchins infected with both diseases in our aquaria were treated with pen/strep; those with BSUD recover [15], however, those with spotting disease do not. The different outcomes from the antibiotic treatment suggest distinct characteristics that may be due to differences in the microbes that underlie the diseases. Based on the descriptions in the literature, there are reports of sea urchin diseases that have been identified as BSUD but should be classified as spotting disease because they describe discrete lesions [6, 7, 19, 20, 55].

The microbial composition associated with spotting disease is distinct from BSUD

Evaluation of the global surface microbiome associated with spotting disease employed the same approaches and methods to evaluate sea urchins in the same aquarium that had BSUD [15]. Consequently, the global surface microbiomes associated with the two diseases can be compared. There is no single taxon that dominates the microbiome composition on sea urchins with BSUD, and consequently, a combination of many taxa likely underlie the infection. Although spotting disease and BSUD have distinct symptoms, the genus Lutibacter is associated with both, suggesting that when sea urchins become compromised, this may be a taxon that is generally prevalent and is therefore likely to invade. Lutibacter belongs to the Flavobacteriaceae family (class Flavobacteriia) that is related to the Cytophagia class in a group previously known as Cytophaga, Flavobacteria, Bacteroides or CFB [56]. These bacteria are commonly found in marine habitats and are associated with certain species of sea urchins [57, 58]. Some members of the Flavobacteriaceae family are known pathogens, such as Flavobacterium psychrophilum that infects salmonids [59, 60]. The Flavobacteriaceae frequently express genes encoding proteases and proteins involved in degrading complex biopolymers, as well as virulence factors for host invasion [61]. It is notable that bacteria closely related to Lutibacter litoralis have been isolated from an American lobster (Homarus americanus) with ulcerative enteritis [62]. Therefore, we speculate that although Lutibacter spp may be commensals associated with sea urchins, they may have a propensity to degrade biopolymers that are found in these animals. Although Colwellia is elevated in the global surface microbiomes of the sea urchins infected with BSUD, in this report it is associated with the healthy sea urchins, suggesting its involvement mainly in BSUD and less in spotting disease. The variations in the taxa associated with these two diseases also suggests that different combinations of bacteria may underpin variations in the symptoms. This may be based on which bacteria become the dominant taxa, and how the proliferative opportunity arises; wound infection for spotting disease vs. altered microbial dynamics for BSUD. It is noteworthy that when the sea urchins recover from BSUD their surface microbiomes change and closely resemble the microbial composition of the fSW [15]. However, the microbial composition of the surface samples from the sea urchins with spotting disease do not share many similarities with the microbes in the fSW samples.

Reproducing discrete spotting disease lesions experimentally requires both surface abrasion and the introduction of microbes into the damaged tissue [7, 8]. Lesions do not develop in the absence of surface abrasion suggesting that both abrasion and microbial infection are required for the development of lesions. Unlike naturally occurring spotting disease however, sea urchins used for experimental induction of spotting disease typically resolve the lesions and recover [7, 8]. Although the appearance of the experimentally induced lesions are similar to naturally derived lesions, the recovery from spotting disease in experimentally induced infections compared to infections that appear in nature suggests that there may be other parameters in play that are required to result in spotting disease in addition to tissue damage. These may include (i) physiological stresses impacting sea urchin defense functions, (ii) differences in the bacterial composition of the species introduced into the injuries, and/or (iii) differences in the microbiomes on the global surface of the experimental sea urchins in aquaria compared to those in nature, which may impact or compete with the microbes that may invade the damaged tissue. Sea urchins infected with spotting disease in natural settings may be exposed to a greater quantity, concentration, or variety of opportunistic taxa that may exacerbate the pathogenicity. Sea urchins in different ecological habitats or in different aquarium systems also show significant differences in their global surface microbiomes [15, 24]. However, in this study, unlike experimentally induced spotting disease, the sea urchins that contracted lesions naturally while in a closed system did not resolve their lesions and did not recover. This suggests that the sea urchins in our aquaria were stressed physiologically, which likely compromises the immune response thereby releasing the microbiomes associated with the sea urchins from immunological control. The outcome may have been the proliferation and invasion of pathogenic bacteria into tissues rather than being restricted to the surface, which led to the development of focal lesions that could not be resolved. Strongylocentrotus purpuratus is commonly found in large clusters in nature that are necessary for successful spawning events and for communal capture of large drift kelp (e.g., see Fig. 19.2c in [63]). Injuries from spines among closely neighboring animals that compete for food may result in ongoing minor injuries that may normally resolve. The basis for whether the injuries progress to lesions that do not resolve under natural settings is likely an indication of poor condition of the sea urchins perhaps resulting from environmental stressors.

Conclusions

Here, we characterize the microbiomes associated with various tissues from sea urchins infected with spotting disease and housed in a closed aquarium. We show that the tissues of diseased and healthy sea urchins have distinct microbiomes, and that the spotting disease infection affects areas of the body wall and surface tissue outside of the lesioned area as well as the coelomic fluid. The combination of surface injuries, microbial infections of the injuries, and a compromised defense system, perhaps as an outcome of physiological stress in sea urchins, combine to result in spotting disease. The characteristic lesions of spotting disease are dominated by Cyclobacteriaceae, Cryomorphaceae, and a few other taxa, which may be the underlying pathogens, which invade injuries on healthy sea urchin tissues that progress to tissue necrosis. The key symptom of spotting disease is discrete surface lesions with restricted surface appendage loss whereas BSUD is a distinct disease with the key symptom of a global surface infection resulting in general spine loss.

Data Availability

The raw sequence reads are available in the Sequence Read Archive database at NCBI under the BioProject ID PRJNA937707. The R code with the complete pipeline is available in the GitHub repository under Spotting_disease_S_purpuratus.

References

  1. Shimizu M, Takaya Y, Ohsaki S, Kawamata K. Gross and histopathological signs of the spotting disease in the sea urchin Strongylocentrotus intermedius. Fish Sci. 1995;61:608–13. https://doi.org/10.1016/S0167-9309(07)80073-1

  2. Tajima K, Hirano T, Shimizu M, Ezura Y. Isolation and pathogenicity of the causative bacterium of spotting disease of sea urchin Strongylocentrotus intermedius. Fish Sci. 1997;63(2):249–52. https://doi.org/10.2331/fishsci.63.249

  3. Zhang W, Lv Z, Li C, Sun Y, Jiang H, Zhao M, Zhao X, Shao Y, Chang Y. Transcriptome profiling reveals key roles of phagosome and NOD-like receptor pathway in spotting diseased Strongylocentrotus intermedius. Fish Shellfish Immunol. 2019;84:521–31. https://doi.org/10.1016/j.fsi.2018.10.042

  4. Li R, Dang H, Huang Y, Quan Z, Jiang H, Zhang W, Ding J. Vibrio coralliilyticus as an agent of red spotting disease in the sea urchin Strongylocentrotus intermedius. Aquac Rep. 2020;16:100244. https://doi.org/10.1016/j.aqrep.2019.100244

  5. Wang L, He B, Chang Y, Ding J. Characterization of the bacterial community associated with red spotting disease of the echinoid Strongylocentrotus intermedius. Aquaculture. 2020;529:735606. https://doi.org/10.1016/j.aquaculture.2020.735606

  6. Gilles KW, Pearce JS. Disease in sea urchin Strongylocentrotus purpuratus: experimental infection and bacterial virulence. Dis Aquat Org. 1986;1:105–14. https://doi.org/10.3354/DAO001105

    Article  Google Scholar 

  7. Roberts-Regan D, Scheibling R, Jellett J. Natural and experimentally induced lesions of the body wall of the sea urchin Strongylocentrotus droebachiensis. Dis Aquat Org. 1988;5:51–62. https://doi.org/10.3354/dao005051

    Article  Google Scholar 

  8. Maes P, Jangoux M. The bald-sea-urchin disease: a biopathological approach. Helgoländer Meeresuntersuchungen. 1984b;37:217–24. https://doi.org/10.1007/BF01989306

    Article  Google Scholar 

  9. Mortensen T, Rosenvinge LK. Sur une algue cyanophycée Dactylococcopsis echini n. sp. parasite dans un oursin. Biologiske Meddelelser. 1934;11(7):1–10.

  10. Lafferty KD. Fishing for lobsters indirectly increases epidemics in sea urchins. Ecol Appl 14:51566–1573. https://doi.org/10.1890/03-5088

  11. Wang YN, Chang YQ, Lawrence JM. Chapter 12: Disease in sea urchins. Dev Aquac Fish Sci. 2013;38:179–86. https://doi.org/10.1016/B978-0-12-3964915.00012-5

    Article  Google Scholar 

  12. Hu F, Zhao C, Ding P, Li Y, Tian R, Qiao Y, Chang Y. An effective facility decreases disease transmission and promotes resistance ability of small sea urchins Strongylocentrotus intermedius: a potential application in the longline culture. Aquaculture. 2022;547:737542. https://doi.org/10.1016/j.aquaculture.2021.737542

  13. Brink M, Rhode C, Macey BM, Christison KW, Roodt-Wilding R. Metagenomic assessment of body surface bacterial communities of the sea urchin Tripneustes gratilla. Mar Genomics. 2019;47:100675. https://doi.org/10.1016/j.margen.2019.03.010

  14. Clemente S, Lorenzo-Morales J, Mendoza JC, López C, Sangil C, Alves F, Kaufmann M, Hernandez JC. Sea urchin Diadema africanum mass mortality in the subtropical eastern Atlantic: role of waterborne bacteria in a warming ocean. Mar Ecol Prog Ser. 2014;506:1–14. https://doi.org/10.3354/meps10829

  15. Shaw CG, Pavloudi C, Barela Hudgell MA, Crow RS, Saw JH, Pyron RA, Smith LC. Bald sea urchin disease shifts the surface microbiome on purple sea urchins in an aquarium. Path Dis. 2023;81:ftad025. https://doi.org/10.1093/femspd/ftad025

    Article  CAS  Google Scholar 

  16. Pearse J, Costa D, Yellin M, Agagian C. Localized mass mortality of red sea urchin Strongylocentrotus franciscanus near Santa Cruz California. Bull US Fish Comm. 1977;75:645–8.

    Google Scholar 

  17. Jangoux M. Diseases of Echinodermata. 1. Agents microorganisms and protistans. Dis Aquat Org. 1987;2:147–62.

    Article  Google Scholar 

  18. Jangoux M. Diseases of Echinodermata. In: Kinne O, editor. Diseases of Marine Animals Biologische Anstalt Helgoland. Hamburg Germany; 1990. pp. 439–567.

  19. Bauer J, Young C. Epidermal lesions and mortality caused by vibriosis in deep-sea Bahamian echinoids: a laboratory study. Dis Aquat Org. 2000;39:193–9. https://doi.org/10.3354/dao039193

    Article  CAS  Google Scholar 

  20. Becker PT, Gillan D, Eeckhaut I. Microbiological study of the body wall lesions of the echinoid Tripneustes gratilla. Dis Aquat Org. 2007;77:73–82. https://doi.org/10.3354/dao01821

    Article  Google Scholar 

  21. Grech D, Mandas D, Farina S, Guala I, Brundu R, Cristo B, Panzalis PA, Falati F, Carella F. Vibrio splendidus clade associated with a disease affecting Paracentrotus lividus (Lamarck 1816) in Sardinia (Western Mediterranean). J Invertebr Pathol. 2022;192:107783. https://doi.org/10.1016/J.JIP.2022.107783

    Article  PubMed  Google Scholar 

  22. Work TM, Dagenais J, Rameyer B, Breeden R, Weatherby TM. Mass mortality of collector urchins Tripneustes gratilla in Hawai’i. Dis Aquat Org. 2023;153:17–29. https://doi.org/10.3354/dao03716

    Article  CAS  Google Scholar 

  23. Lorgen-Ritchie M, Uren Webster T, McMurtrie J, Bass D, Tyler CR, Rowley A, Martin SA. Microbiomes in the context of developing sustainable intensified aquaculture. Front Microbiol. 2023;14:1200997. https://doi.org/10.3389/fmicb.2023.1200997

    Article  PubMed  PubMed Central  Google Scholar 

  24. Wessel GM, Kiyomoto M, Reitzel AM, Carrier TJ. Pigmentation biosynthesis influences the microbiome in sea urchins. Proc R Soc B. 2022;289:20221088. https://doi.org/10.6084/m9.figshare.c.6135664.v1

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  25. Turner CR, Miller DJ, Coyne KJ, Corush J. Improved methods for capture extraction and quantitative assay of environmental DNA from Asian bigheaded carp (Hypophthalmichthys spp). PLoS ONE. 2014;9(12):e114329. https://doi.org/10.1371/journal.pone.0114329

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Weisburg WG, Barns SM, Pelletier DA, Lane DJ. 16S ribosomal DNA amplification for phylogenetic study. J Bacteriol. 1991;173(2):697–703.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Nadkarni MA, Martin FE, Jacques NA, Hunter N. Determination of bacterial load by real-time PCR using a broad-range (universal) probe and primers set. Microbiology. 2002;148:257–66.

    Article  PubMed  CAS  Google Scholar 

  28. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3. https://doi.org/10.1038/nmeth.3869

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glӧckner FO. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nuc Acids Res. 2013;41(Database issue):D590–6. https://doi.org/10.1093/nar/gks1219

    Article  CAS  Google Scholar 

  30. McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8(4):e61217. https://doi.org/10.1371/journal.pone.0061217

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing Vienna Austria. 2021. https://www.R-project.org/

  32. Chao A. Nonparametric estimation of the number of classes in a population. Scand J Stat. 1984;11(4):265–70. https://doi.org/10.2307/4615964

    Article  Google Scholar 

  33. Chao A, Lee SM. Estimating the number of classes via sample coverage. J Am Stat Assoc. 1992;87:210–7.

    Article  Google Scholar 

  34. Qian PY, Cheng A, Wang R, Zhang R. Marine biofilms: diversity interactions and biofouling. Nat Rev Microbiol. 2022;20:671–84. https://doi.org/10.1038/s41579-022-00744-7

    Article  PubMed  CAS  Google Scholar 

  35. Faddetta T, Ardizzone F, Faillaci F, Reina C, Palazzotto E, Strati F, De Filippo C, Spinelli G, Puglia AM, Gallo G, Cavalieri V. Composition and geographic variation of the bacterial microbiota associated with the coelomic fluid of the sea urchin Paracentrotus lividus. Sci Rep. 2020;10:21443. https://doi.org/10.1038/s41598-020-78534-5

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Offret C, Jégou C, Mounier J, Fleury Y, Le Chevalier P. New insights into the haemo-and coelo‐microbiota with antimicrobial activities from Echinodermata and Mollusca. J Appl Microbiol. 2019;126(4):1023–31.

    Article  PubMed  CAS  Google Scholar 

  37. Hakim JA, Schram JB, Galloway AW, Morrow CD, Crowley MR, Watts SA, Bej AK. The purple sea urchin Strongylocentrotus purpuratus demonstrates a compartmentalization of gut bacterial microbiota predictive functional attributes and taxonomic co-occurrence. Microorg. 2019;7:35. https://doi.org/10.3390/microorganisms7020035

    Article  CAS  Google Scholar 

  38. Hakim JA, Green GBH, Watts SA, Crowley MR, Morrow CD, Bej AK. Microbial composition and genes for key metabolic attributes in the gut digesta of sea urchins Lytechinus variegatus and Strongylocentrotus purpuratus using shotgun metagenomics. Curr Issues Mol Biol. 2021;43(2):978–95. https://doi.org/10.3390/cimb43020070

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  39. Staley JT, Konopka A. Measurement of in situ activities of nonphotosynthetic microorganisms in aquatic and terrestrial habitats. An Rev Microbiol. 1985;39:321–46.

    Article  CAS  Google Scholar 

  40. Becker PT, Egea E, Eeckhaut I. Characterization of the bacterial communities associated with the bald sea urchin disease of the echinoid Paracentrotus lividus. J Invertebr Pathol. 2008;98(2):136–47. https://doi.org/10.1016/j.jip.2007.12.002

    Article  PubMed  CAS  Google Scholar 

  41. Masuda Y, Tajima K, Ezura Y. Resuscitation of Tenacibaculum sp. the causative bacterium of spotting disease of sea urchin Strongylocentroutus intermedius from the viable but non-culturable state. Fish Sci. 2004;70(2):277–84. https://doi.org/10.1111/j.1444-2906.2003.00801.x

  42. Wang B, Li Y, Li X, Qu JY, Zhao XM. Pathogenic mechanism of causative Vibrio found in red spotting diseased sea urchin Strongylocentrotus intermedius. J Dalian Fisheries Univ. 2005;20:11–5.

  43. Hira J, Stensvåg K. Evidence for association of Vibrio echinoideorum with tissue necrosis on test of the green sea urchin. Strongylocentrotus droebachiensis Sci Rep. 2022;12(1):4859. https://doi.org/10.1038/s41598-022-08772-2

  44. Kiselev KV, Ageenko NV, Kurilenko VV. Involvement of the cell-specific pigment genes pks and sult in bacterial defense response of sea urchins Strongylocentrotus intermedius. Dis Aquat Org. 2013;103:121–32. https://doi.org/10.3354/dao02570

  45. Weigel B. Sea cucumber intestinal regeneration reveals deterministic assembly of the gut microbiome. Appl Environ Microbiol. 2020;86(14):e00489–20. https://doi.org/10.1128/AEM.00489-20

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Zanotti AA, Gregoracci GB, Kitahara MV. The microbial profile of a tissue necrosis affecting the Atlantic invasive coral Tubastraea tagusensis. Sci Rep. 2021;11:9828. https://doi.org/10.1038/s41598-021-89296-z

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  47. Cottrell MT, Kirchman DL. Natural assemblages of marine proteobacteria and members of the Cytophaga-Flavobacter cluster consuming low- and high-molecular-weight dissolved organic matter. Appl Environ Microbiol. 2000;66(4):1692–7. https://doi.org/10.1128/AEM.66.4.1692-1697.2000

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  48. Kirchman DL. The ecology of Cytophaga-Flavobacteria in aquatic environments. FEMS Microbiol Ecol. 2002;39(2):91–100. https://doi.org/10.1111/j.1574-6941.2002.tb00910.x

    Article  PubMed  CAS  Google Scholar 

  49. Xie G, Bruce DC, Challacombe JF, Chertkov O, Detter JC, Gilna P, Han CS, Lucas S, Misra M, Myers GL, Richardson P, Tapia R, Thayer N, Thompson LS, Brettin TS, Henrissat B, Wilson DB, McBride MJ. Genome sequence of the cellulolytic gliding bacterium Cytophaga hutchinsonii. Appl Environ Microbiol. 2007;73(11):3536–46. https://doi.org/10.1128/AEM.00225-07

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  50. Bai X, Wang X, Wang S, Ji X, Guan Z, Zhang W, Lu X. Functional studies of β-glucosidases of Cytophaga hutchinsonii and their effects on cellulose degradation. Front Microbiol. 2017;8:140. https://doi.org/10.3389/fmicb.2017.00140

    Article  PubMed  PubMed Central  Google Scholar 

  51. Welsh RM, Rosales SM, Zaneveld JR, Payet JP, McMinds R, Hubbs SL, Vega Thurber RL. Alien vs. predator: bacterial challenge alters coral microbiomes unless controlled by Halobacteriovorax predators. PeerJ. 2017;5:e3315. https://doi.org/10.7717/peerj.3315

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Rosales SM, Miller MW, Williams DE, Traylor-Knowles N, Young B, Serrano XM. Microbiome differences in disease-resistant vs. susceptible Acropora corals subjected to disease challenge assays. Sci Rep. 2019;9:18279. https://doi.org/10.1038/s41598-019-54855-y

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  53. Gignoux-Wolfsohn SA, Aronson FM, Vollmer SV. Complex interactions between potentially pathogenic, opportunistic, and resident bacteria emerge during infection on a reef-building coral. FEMS Microbiol Ecol. 2017;93(7):fix080. https://doi.org/10.1093/femsec/fix080

    Article  CAS  Google Scholar 

  54. Bass D, Stentiford GD, Wang H, Koskella B, Tyler CR. The pathobiome in animal and plant diseases. Trends Ecol Evol. 2019;34:996–1008. https://doi.org/10.1016/j.tree.2019.07.012

    Article  PubMed  PubMed Central  Google Scholar 

  55. Maes P, Jangoux M. The bald sea urchin disease: a bacterial infection. In: Keegan BF, O’Connor BDS, editors. Echinodermata. CRC press. Boca Raton FL; 1984a. pp. 313–4.

  56. García-López M, Meier-Kolthoff JP, Tindall BJ, Gronow S, Woyke T, Kyrpides NC, Hahnke RL, Göker M. Analysis of 1,000 type-strain genomes improves taxonomic classification of Bacteroidetes. Front Microbiol. 2019;10:2083. https://doi.org/10.3389/fmicb.2019.02083

    Article  PubMed  PubMed Central  Google Scholar 

  57. Choi DH, Cho BC. Lutibacter litoralis gen. nov., sp. nov., a marine bacterium of the family Flavobacteriaceae isolated from tidal flat sediment. Int J Syst Evol Microbiol. 2006;56(Pt 4):771–6. https://doi.org/10.1099/ijs.0.64146-0

    Article  PubMed  CAS  Google Scholar 

  58. Nedashkovskaya OI, Van Trappen S, Zhukova NV, De Vos P. Lutibacter holmesii sp. nov., a marine bacterium of the family Flavobacteriaceae isolated from the sea urchin Strongylocentrotus intermedius, and emended description of the genus Lutibacter. Int J Syst Evol Microbiol. 2015;65(11):3991–6. https://doi.org/10.1099/ijsem.0.000525

  59. Starliper CE. Bacterial coldwater disease of fishes caused by Flavobacterium psychrophilum. J Adv Res. 2011;2(2):97–108. https://doi.org/10.1016/j.jare.2010.04.001

    Article  Google Scholar 

  60. Park J, Roh H, Lee Y, Lee J-Y, Kang H-Y, Seong MJ, Kang Y, Heo YU, Kim BS, Park C-I, Kim D-H. Characterization and pathogenicity of Flavobacterium psychrophilum isolated from rainbow trout (Oncorhynchus mykiss) in Korea. Microorg. 2023;11:2546. https://doi.org/10.3390/microorganisms11102546

    Article  CAS  Google Scholar 

  61. Guérin C, Lee BH, Fradet B, van Dijk E, Mirauta B, Thermes C, Bernardet J-F, Repoila F, Duchaud E, Nicolas P, Rochat T. Transcriptome architecture and regulation at environmental transitions in flavobacteria: the case of an important fish pathogen. ISME Commun. 2021;1:33. https://doi.org/10.1038/s43705-021-00029-9

    Article  PubMed  PubMed Central  Google Scholar 

  62. Battison AL, Després BM, Greenwood SJ. Ulcerative enteritis in Homarus americanus: case report and molecular characterization of intestinal aerobic bacteria of apparently healthy lobsters in live storage. J Invertebr Pathol. 2008;99(2):129–35. https://doi.org/10.1016/j.jip.2008.06.013

    Article  PubMed  CAS  Google Scholar 

  63. Smith LC, Byrne M, Gedan KB, Lipscomb DL, Majeske AJ, Tafesh-Edwards G. Ecological outcomes of echinoderm disease, mass die-offs, and pandemics. In: Rowley AF, Coates CJ, Whitten MMA, editors. Invertebrate Pathology. Oxford UK: Oxford University Press; 2022. pp. 563–88.

    Chapter  Google Scholar 

Download references

Acknowledgements

This work was completed in part with resources provided by the High Performance Computing Cluster at George Washington University, Information Technology, Research Technology Services. The authors are grateful to Dr. Mark Reeves for providing space in his lab for collecting microbes with the vacuum filtering system.

Funding

This research was supported by funding from the US National Science Foundation (NSF IOS 1855747) to LCS, a Sigelman Undergraduate Research Enhancement award and a Luther Rice Undergraduate Research Fellowship to CGS.

Author information

Authors and Affiliations

Authors

Contributions

CGS collected the samples, evaluated the gDNA quality, analyzed the data, generated figures and tables, and wrote the paper. CP carried out the amplicon sequence matching pipeline, analyzed the data, and edited the paper. RSC optimized data visualization. JHS guided the microbiology analysis, helped with bioinformatics, and edited the paper. LCS isolated the BAC DNA for sequencing, acquired funding, and wrote and edited the paper.

Corresponding author

Correspondence to L. Courtney Smith.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1: Table S1.

Relative abundance of phyla on the global surface microbiome samples. Table S2. Genera with an average relative abundance > 0.1% shown as the relative abundance per sample for the global surface microbiome samples. Table S3. Significantly differentially abundant taxa as identified by LEfSe for diseased and healthy global surface microbiome samples. Table S4. Relative abundance of phyla for the LBW and LS samples. Table S5. Genera with an average relative abundance of > 0.1% in the LBW and LS microbiome samples. Table S6. Significantly differentially abundant taxa as identified by LEFSe for the LBW and LS microbiome samples. Table S7. Relative abundance of phyla in the tissue microbiome samples. Table S8. Genera with an average relative abundance of > 0.1% in the tissue microbiome samples. Table S9. Significantly differentially abundance taxa as identified by LEfSe in the tissue microbiome samples

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shaw, C.G., Pavloudi, C., Crow, R.S. et al. Spotting disease disrupts the microbiome of infected purple sea urchins, Strongylocentrotus purpuratus. BMC Microbiol 24, 11 (2024). https://doi.org/10.1186/s12866-023-03161-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-023-03161-9

Keywords