DA confusion matrices.
Epizootic shell disease (ESD) is a persistent malady that affects American lobsters (Homarus americanus) in the southern extent of the commercial fishery. Emerging at the turn of the 21st century, ESD presented as bacterial ulcerations on the carapace of affected lobsters. The research presented here examined the bacterial community of the lobster carapace and represented the first such attempt to characterize the lobster surface microbiome. Culture‐independent techniques, such as amplicon length heterogeneity and pyrosequencing, yielded sequence data of hypervariable regions of the genes for ribosomal RNA that upon comparison revealed the likely identities of the taxa present on the lobster carapace. Although some researchers have identified a novel chitinolytic bacterium of the genus Aquimarina (A. homaria) as consistently appearing on lobsters with shell disease, this research found no evidence of a correlation of this species with the disease. Instead, analysis revealed that the genus Aquimarina was ubiquitous and correlated only weakly with the diseased state. The data suggest that this disease is not caused by a single pathogen but by a state of dysbiosis where normally occurring microflora emerge as potential opportunistic pathogens when there is some apparent environmental stressor that alters the interaction of the surface biofilm of the lobster.
- epizootic shell disease
- lobster microbiome
- next‐generation sequencing
- multitag pyrosequencing
- amplicon length heterogeneity
Since the end of the 20th century, American lobsters (
One such challenge to lobster health that is reducing the quality, and possibly the quantity, of commercial landings is epizootic shell disease (ESD), a condition that is manifesting itself in geographically isolated portions of the commercial lobster fishing grounds .
Documentation of diseases that cause lesions and ulcerations in the carapace of lobsters can be found at least as early as 1937 . Initial reports of shell disease in other crustaceans are contemporary with Sindermann . In many cases, the pathology appeared to be linked to degraded environments, such as the proximity to pollution sources, for example, the municipal waste dump site in an offshore canyon used by New York City where lobsters and crabs presented with shell disease [6, 7]. Another common place to find crustacean shell diseases is in commercial impoundment facilities . In this situation, chitinolytic bacteria appeared to be eroding the carapace faster than it could be replaced by molting, leading to a pitting of the cuticle. However, there has only been one report of an infectious process after the lobster's carapace was abraded and swabbed with
In the epicuticle, calcification consists of spherulitic calcite surrounded by a lipid‐protein matrix; in the exocuticle and endocuticle, calcite crystals are dispersed throughout chitin‐protein fibers that are referred to as lamellae. In addition, the exocuticle contains trabeculae that are composed of apatite (calcium phosphate) that resembles spongy bone tissue in vertebrates . The endocuticle is the most calcified layer of the exoskeleton .
According to Smolowitz et al. , ESD is identifiable by erosions of the carapace that have a unique histology. The presence of lesions on the carapace of an infected lobster is the most visible sign of the disease. Histological examination revealed bacteria as the most common organisms in the lesions, whereas some protist constituents were found in more advanced cases. The lesions were grouped into three categories based on the depth of the bacterial incursion, presuming that depth reflects a progressive erosion of the cuticle.
Category 1 is the least severe erosion, with shallow lesions extending into the epicuticle and exocuticle. The margins of the lesions often exhibited evidence of melanization, but inflammation in the underlying connective tissue or other evidence of an immune response is rarely observed at this stage. Bacteria are found in the leading edges of the lesion and in the crystalline chitin lattice.
Category 2 lesions are moderately deep, penetrating the calcified endocuticle. The crystalline lattice structure of the chitin takes on a “pillar‐like” appearance as the bacteria degrade the protein structure between the lattice crystals. The endocuticle exhibits melanization, especially in the vertical areas of bacterial incursion. Evidence of immune response to the infection includes inflammation of the underlying cuticular epithelium and “moderate numbers” of hemocytes in the tissues. Secondary invasion by small protists occurs during this stage. These organisms are apparently responsible for the degradation of the crystal lattice structure of the chitin. In some category 2 lesions, an “inflammatory cuticle” forms between the uncalcified endocuticle and the cuticular epithelium. The latter has some areas of hyperplasia and hypertrophy.
Deeper erosions into the uncalcified endocuticle are characterized as category 3 lesions. At this stage, the overlying structures of the carapace are absent, and the exposed areas are melanized. The cuticular epithelium is hyperplastic and hypertrophic and exhibits an intense inflammatory response accompanied by greater hemocyte infiltration. The underlying connective tissue also exhibited signs of inflammation and immune response. The most extreme types of category 3 lesions had no more tissue than the inflammatory cuticle overlying the cuticular epithelium. In some cases, the lesions progressed to ulcerations, which were characterized by a complete absence of cuticular tissue and cuticular epithelium. Degranulated hemocytes developed a multilayered pseudomembrane to cover the connective tissue. The outer layer of the pseudomembrane was necrotic and melanized.
Certain groups of bacteria have been found as common constituents on the carapaces of moribund lobsters, and overall populations of bacteria are higher than on the carapaces of healthy lobsters . Using a combination of culture‐dependent techniques, polymerase chain reaction (PCR), and denaturing gradient gel electrophoresis, the group of Chisosterdov isolated three novel bacteria that were ubiquitous in ESD lesions:
The characterization of microbial communities using terminal restriction fragment length polymorphism revealed that, although no grossly significant difference in communities on healthy lobsters was present compared to diseased lobsters, there were some minor differences . More anaerobic bacteria and greater numbers of α‐ and β‐proteobacteria were found in the lesions. Bell et al. used lobsters from eastern and western Long Island Sound. Lobsters from a coastal Maine site were used as reference specimens. This study also elucidated the activities of four bacterial ectohydrolases on shell samples from healthy and diseased specimens and found that, whereas chitinase activity was high in all samples, cellulase and proteinase activities were significantly higher in diseased specimens. Lipase activity was higher in Long Island Sound lobsters compared to those from Maine but was apparently similar between healthy and diseased lobsters. They concluded that chitin degradation was not as important to the progression of the disease as the degradation of lipids and protein in the carapace. Bell et al.  interpreted the spatial difference in lipase activity in bacteria to indicate that the degradation of the lipid moieties in the epicuticle may be an important initial step in disease progression and hypothesized that it should correlate positively with geographic areas of ESD prevalence.
The works of Smolowitz et al.  and Bell et al.  make clear that ESD differs from earlier types of shell disease. ESD is ostensibly bacterial in origin, but other research indicates that the bacterial infection may only be a proximate cause.
Shields et al.  found that some idiopathic conditions were more prevalent in lobsters from Rhode Island compared to lobsters from Maine. Rhode Island lobsters presented with vibriosis, hepatopancreatitis, and eye lesions at a significantly higher rate than did Maine lobsters. When Rhode Island lobsters with and without ESD were compared, these and other idiopathic conditions were present but did not occur more frequently in lobsters with the disease.
Laufer et al.  found that ESD lobsters had higher levels of alkylphenols than did unaffected lobsters and that bottom sediments had higher than normal levels of these compounds in areas where diseased lobsters resided. Using gas chromatography‐mass spectrometry, Laufer et al. identified four alkylphenol species: 2‐
Tarrant et al. found that lobsters with ESD had gene expression patterns that were consistent with exposure to xenobiotics. Specifically, they found an elevated expression of ecdysteroid receptors and elevated expression of CYP450 (a cytochrome
In response to the ESD crisis, a group of scientists, fisheries managers, and lobstermen formed the New England Lobster Health Initiative. It was this ad hoc committee that ultimately received funding that was awarded to several university and institutional researchers, including the award that funded the study reported here. The findings of each group were published as a special edition of the
A systems approach to ESD was employed in this study to elucidate the surface biofilms of lobsters in three states of health: Healthy, Diseased, and Healthy‐on‐Diseased (areas of the carapace on a lobster presenting with ESD that do not yet have the characteristic lesions). Using culture‐independent molecular techniques coupled with multivariate statistical and network analysis models, this study surveyed the biofilms of the carapaces of lobsters from the geographic area that is most affected by ESD.
This study represents the first survey of the surface microbiome of the American lobster. Although the sequencing techniques can provide a comprehensive overview of bacteria in a sample, this study cannot be regarded as a global assessment of the lobster microbiome. The lobsters surveyed here are from a one geographic location and represent a statistically small sample of lobsters. Elucidating the complete lobster microbiome should make use of lobsters from as much of the natural range as possible and should include a statistically large sample of specimens. Bell et al. , for example, have demonstrated that there is a biochemical activity in the lobster surface microbiome that varies geographically.
This survey does, however, offer a glimpse into the subtle microbiome differences between lobsters that present with ESD (“Diseased”) and those that do not (“Healthy”). In addition, a third class, “Healthy‐on‐Diseased,” has been defined as a region of the shell of a lobster with ESD that does not present with lesions. This third class may provide some insight into the disease progression.
As described previously, Quinn et al.  were successful in initiating lesions similar to ESD with three bacteria that were isolated from ESD lesions. Some of the infections were initiated with only one bacterium (
The polymicrobial nature of ESD, and the desire to illuminate something of its etiology, is the focus of this research. To that end, culture‐independent techniques were employed to elucidate the bacterial community, and multivariate statistics were used to interpret the results. Specifically, length heterogeneity‐PCR (LH‐PCR) was used to survey a typical lobster. This provided an overview of the distribution of bacteria on different regions of the lobster carapace. LH‐PCR data do not provide sequence information but identify PCR amplicons by size (basepair length) that indicate different taxa. The number of representative fragments found in the cephalothorax region of the carapace, combined with the knowledge that this region had the highest incidence of ESD lesions, led to the decision to focus the next phase of molecular interrogation on the cephalothorax. This phase used multitag pyrosequencing (MTPS) to generate raw sequence data of each PCR amplicon, which could be identified by comparing these sequences to known sequences in a standard database.
Discriminant analysis (DA) was chosen as a means of determining if there is a statistically significant difference between the bacterial communities found on the three classes of data (Healthy, Diseased, and Healthy‐on‐Diseased). In addition, DA assigns each variable (bacterial taxon) with coefficients that indicate its contribution to the discriminant function. Correlational network analysis, which uses techniques similar to social network analysis to create a visual image of correlations between the bacterial taxa, was employed to identify which taxa are positively and negatively correlated in each of the three states examined. Using correlational difference analysis to compare the states, it may be possible to more clearly identify the changes in bacterial population on the lobster carapace that are correlated with a decline in health.
We approached the problem using the following hypothesis:
H: The bacterial communities on lobsters with ESD are significantly different in quality and quantity than unaffected lobsters.
HA: The bacterial communities on the carapaces of healthy and shell‐diseased lobsters are similar.
2. Materials and methods
We initially performed a survey of the microbiome on normal lobsters using a culture‐independent method to characterize the taxa composition of the biofilm on various locations of the lobster shell. Characterizing a microbial population using molecular techniques is advantageous over culture‐dependent techniques in that the latter relies on the growth of bacteria on artificial media usually at a controlled temperature. Many bacteria have particular nutritional and environmental conditions that render them difficult to grow in laboratory cultures and would not be present in cultures taken from the environment. Existence in culture medium, as discussed in the previous chapter, requires that the bacteria assume a planktonic phase rather than embedding in a biofilm, which may be an abnormal state for most species . These molecular techniques are sometimes referred to as “culture independent” in that they identify bacteria from a sample without first attempting to grow them under artificial conditions. Identifying bacteria by amplifying and sequencing the entire genome of each individual would be impractical and inordinately time‐consuming. Several techniques exist that enable the identification of bacteria by taking advantage of variable regions within genes that are ubiquitous in all prokaryotes. This research used the first two variable regions of the genes encoding for the small (16S) subunit of ribosomal RNA (rRNA) to discriminate between different operational taxonomic units (OTUs) that are surrogates for taxa identification. There are nine such variable regions within this gene.
The technique employed in this research is called LH‐PCR . The variability, or heterogeneity, of the basepair lengths of the variable regions can be used to define the OTUs that represent individual taxa, but the technique does not assign a taxon name to them.
An LH‐PCR survey of the diversity of the bacteria in the shell microbiome of a store‐bought lobster was initially performed to determine the compositional patterns throughout the surface of the lobster. In this procedure, bacterial DNA was extracted from various regions of the lobster carapace (claw, abdomen, cephalothorax, and telson). The carapace samples were dissolved in EDTA and proteinase K to recover all microorganisms from the surface and subsurface, and DNA was extracted using the FastDNA Spin Kit (MP Biomedicals, Solon, OH). PCR was employed to amplify the bacterial genes for the first two hypervariable regions (V1 and V2) of the 16S rRNA subunit using universal primers 27F (5'‐AGA GTT TGA TCM TGG CTC AG‐3') and 355R (5'‐GCT GCC TCC CGT AGG AGT‐3'; Invitrogen Corporation, Carlsbad, CA) .
The LH‐PCR products were diluted according to their intensity on agarose gel electrophoresis and mixed with ILS‐600 size standards (Promega, Madison, WI) and HiDi Formamide (Applied Biosystems, Foster City, CA). The diluted samples were then separated on an ABI 3130xl fluorescent capillary sequencer (Applied Biosystems) and processed using the Genemapper™ software package (Applied Biosystems). Normalized peak areas were calculated using a custom PERL script, and OTUs constituting less than 1% of the total community from each sample were eliminated from the analysis to remove the variable low abundance components within the communities.
The technique described above employed a culture‐independent technique that is capable of resolving amplified fragments of DNA by size (OTUs). This can be useful in making observations on gross changes of microbial communities over time, such as in a site that is contaminated by petroleum or some other organic material. In this case, LH‐PCR was used to demonstrate that the biofilm on the cephalothorax was the best representative sample of the shell microbiome. The next phase of interrogation employed MTPS  to further elucidate the nature of the microbial community. The OTUs from the lobster samples were sequenced using next‐generation sequencing technology, and the sequence reads were matched as closely as possible to known bacterial taxa. This provided us with data that could be used in a variety of ways to demonstrate the correlation between bacterial taxa and various states of health in the lobster.
Cuticle samples (0.5 cm2) were harvested as part of the “100 Lobsters” Project  and in a manner described therein. The lobsters were collected from within Narragansett Bay, Rhode Island. A total of 55 lobsters had ESD, whereas 47 were categorized as apparently healthy (i.e., the lobsters did not have visible ESD lesions). The samples were shipped in plastic vials on dry ice and stored at ‐80°C until processed. Three types of carapace samples were taken: those from lesions of diseased lobsters (Diseased), those from lesion‐free areas of diseased lobsters (Healthy‐on‐Diseased), and those from apparently healthy lobsters (Healthy). The carapace samples were dissolved in EDTA and proteinase K to recover all microorganisms from the surface and subsurface, and DNA was extracted using the FastDNA Spin Kit (MP Biomedicals).
PCR was employed to amplify the bacterial genes for two hypervariable regions of the rRNA subunit using modified primers 27F (5'‐AGA GTT TGA TCM TGG CTC AG‐3') and 355R (5'‐GCT GCC TCC CGT AGG AGT‐3'; Invitrogen Corporation). The forward primer listed above was modified with a “barcode” sequence that was unique to each DNA sample extracted from individual lobsters, allowing us to track the various bacterial sequence reads to their specific sample origin (i.e., the lobster and carapace region from which the bacteria were sampled). MTPS was employed to characterize the microbiome from a subset of the carapace samples that were used in the LH‐PCR analysis. We generated a set of 96 fusion primers that contained the emulsion PCR linkers (454 Life Sciences) on the 27F and 355R universal 16S rRNA primers along with different eight‐base “barcodes” on the 27F primer. Each lobster DNA sample was then amplified with a unique set of tagged 16S rRNA primers, pooled, subjected to emulsion PCR, and pyrosequenced using a GS‐FLX pyrosequencer as per manufacturer's instructions (Roche, Branchburg, NJ). Data from the pooled sample were “deconvoluted” by sorting the sequences into sample bins based on the barcodes using custom PERL scripts. This technique allows the rapid sequencing of multiple samples at one time, yielding thousands of sequence reads per sample. The sequence reads were identified using the Bayesian analysis from the Ribosomal Database Project . We used a custom PERL script to normalize the abundances of the taxa in a sample based on the total reads in that sample.
The sequences were aligned that were identified as
The software program Quantitative Insights into Microbial Ecology (QIIME)  was used to compare the microbial communities on the three types of carapace samples and to construct a UniFrac neighbor joining tree that graphically displays the similarities in microbiome structure between each lobster sample in the study. The technique was validated using jackknifing to estimate the dispersion of the data . According to Lozupone and Knight , the UniFrac metric “measures the phylogenetic distance between sets of taxa in a phylogenetic tree as a fraction of the branch length that leads to descendents from either one environment or the other, but not both.” In the process, UniFrac was able to depict the evolutionary differences between environments. In this application, UniFrac was used to compare phylogenetic differences between bacterial communities and provide a graphic representation of the differences. Each terminal node, therefore, represents a unique phylogenetic tree of that sample. The distance between the nodes is the distance between these trees.
DA, sometimes referred to as discriminant function analysis, is a multivariate statistical method that identifies variables of different cases that are useful in the discrimination of those cases into different classes. The analysis of the data results in the construction of a linear equation in which the variables are factors of both unstandardized and standardized coefficients. The products of these variable coefficient pairs are summed to produce a score for each case, and that score determines the membership of that case in a specific class . The summation of the unstandardized coefficients includes a constant. Standardizing the coefficients eliminates the constant. These coefficients are partial coefficients that compare the relative importance of the independent variables . We define the cases as individual lobsters in the study and the variables as the bacterial taxa that were identified by MTPS.
In addition to producing unstandardized and standardized linear coefficients that can be used to discriminate classes, DA produces full coefficients that are referred to as structure coefficients. These are measurements of the importance of the independent variables and the discriminant function itself [34, 35]. Structure coefficients are a more definitive coefficient for determining relative importance of the variables . In this research, the structure coefficients are useful in providing us with a ranking of the variables, which reflects their relative importance to ESD.
The abundance data for the 170 variables (taxa) from the 102 lobster samples were analyzed by DA with PASW version 18 (IBM, Chicago, IL). In the first analysis, we compared the taxa on the cuticle from Diseased versus Healthy animals, the second analysis compared the taxa on the cuticle from Diseased animals versus Healthy‐on‐Diseased animals, and the third analysis compared the taxa on the cuticle from Healthy animals versus Healthy‐on‐Diseased animals.
2.4. Correlational network analysis
The computations described above are useful in identifying variables (bacterial taxa) that are likely candidates for determining the class of each sample, in this case, the state of health of the lobster being sampled. DA, however, does not account for the effects of one variable upon another. The members of a microbial community interact not only with their host or source of food but also with each other .
Network analysis is a statistical‐graphical technique that has wide applications in fields as diverse as the social sciences and epidemiology and has gained some popularity among people working in “homeland security.” Arguably the most common of these types of applications is social network analysis, which can be defined as a “study of human relationships by means of graph theory” .
Microbes occupying a biofilm are exposed to the biochemical output of their neighboring microorganisms. Some of these molecules are metabolic wastes of one microbe but serve as nutritive material for another, which can be viewed as metabolic cooperativity . Members of a community also secrete antimicrobial agents in an effort to reduce competition, and still other bacteria engage in quorum sensing, the detection of chemical signals as a means of monitoring population density . It would be incorrect to equate human interactions with those of microorganisms, but it is arguably valid to recognize that biofilms exhibit a social structure that is mediated via biochemical signaling. Therefore, network analysis has the potential to elucidate microbial interactions within a biofilm. Within this study, the chemical interaction between various members of the biofilm was not directly investigated. Instead, we calculated the statistical correlations between bacterial taxa as a means of inferring their interactions and relationships.
The pyrosequencing abundance data were analyzed using Spearman rank correlation using a custom R module . Spearman rank correlation is a nonparametric statistical function that allows the comparison of nonlinear data . All possible correlation coefficients between variables are calculated, and correlations between variables are ranked in an output table. The parameters for the ranking were set as follows:
The correlation tables were generated by the R module into the networking software Cytoscape  version 2.8.2 to construct correlational network diagrams of the three states of health of the lobster microbiome: Diseased, Healthy, and Healthy‐on‐Diseased. In addition, the data were used to construct correlation difference networks that plot edges between features whose correlations have significantly changed between the two states being compared. One can interpret these correlation difference networks as a depiction of what correlations were altered between the two states. The author constructed two such maps: Healthy versus Diseased and Healthy‐on‐Diseased versus Diseased.
Initially, the lobster carapace was surveyed to determine the bacterial community distribution on various regions of the carapace. Figure 1 is a histogram of the normalized abundance of OTUs that were present on the various regions of the carapace (pink and red: cephalothorax, yellow: claws, green: abdomen, and tail: blue). The LH‐PCR profiles revealed that, with a few exceptions, the OTUs present on the cephalothorax are representative of the microbiome present on other regions, such as the claw and the abdomen. Of the 29 OTUs identified, 22 are present on the cephalothorax, 17 are present on the claws, 17 are present on the tail (telson), and 13 are present on the abdomen. Using the Shannon index of diversity , the cephalothorax region had a index of 2.81 and an eveness score of 0.85. The claw region had an index value of 2.80 and a eveness score of 0.833. The tail had an index value of 2.37 and an eveness score of 0.75. The abdomen had an index value of 1.95 and a eveness score of 0.61.
The peaks in the fingerprints are described as OTUs, because this method alone does not identify the bacteria. In fact, a peak may contain more than one species, and strains of the same species may show up in different peaks. The normalized OTU abundance for each sample was plotted as a stacked histogram to demonstrate the distribution of the microbiome on the lobster carapace.
MTPS was performed on 102 lobster samples and the resulting sequence reads were sorted based on their tags or barcodes. The analysis yielded 212,019 reads with an average of 1594 reads per sample. We identified 170 bacteria present on the cuticle of lobsters from Narragansett Bay, Rhode Island, having culled taxa that were less than 1% of the total community under the
Figure 2 is a histogram of the average abundance of the taxa found in each sample class (i.e., Diseased, Healthy‐on‐Diseased, and Healthy) rank ordered by the average abundance of the Disease class (red bars). The entire histogram for all 170 taxa is depicted in the insert. The most abundant taxa found in the carapace microbiome include the genera
Figure 3 is a neighbor joining tree that displays all of the members of the genus
The weighted UniFrac tree is displayed in Figure 4. Although a cursory inspection of the figure may lead one to conclude that there is indeed a clustering of the three disease classes (Diseased, Healthy, and Healthy‐on‐Diseased), significance tests reveal that the differences between the three classes are not significant (Diseased vs. Healthy‐on‐Diseased:
Figure 5 is a principal coordinate analysis (PCO) of the normalized abundances of the genera identified in each disease class using a Bray‐Curtis distance metric. The PCO performs an Eigen analysis that clusters the data based on the variance of all features (genera) of the data matrix . One can see a clustering of the Diseased samples (red dots) that is somewhat distinct from the Healthy (yellow) and Healthy‐on‐Diseased (blue) samples.
The DA for the first comparison, the bacterial taxa on cuticle from Diseased versus apparently Healthy lobsters, included a tolerance test that eliminated 112 variables because they lacked variance between groups and thus did not significantly contribute to the discriminant function. Thus, the analysis between the two groups relied on 58 variables. The functions at the group centroids were 1.463 for disease class 1 (Diseased) and ‐1.1712 for disease class 2 (Healthy). The group centroids are a central measure of a set of multivariate data similar to the mean in univariate analyses. These scores demonstrate that the centroids were well separated. Such a separation between centroids demonstrates that the function is able to discriminate one class from another. The canonical correlation coefficient (
The classification results (confusion matrix) of all three analyses are displayed in Table 1. The confusion matrix is run as a test of the ability of the discriminant function to predict the membership of the cases accurately. In the first analysis (Diseased vs. Healthy), the function predicted 94.1% of the original grouped cases correctly. The second analysis (Diseased vs. Healthy‐on‐Diseased) predicted 98.9% and the third analysis (Healthy vs. Healthy‐on‐Diseased) predicted 100% of original grouped cases correctly. This is further evidence that the respective functions for each of these comparisons are reliably discriminating between cases.
|Classification results: Diseased vs. Healthy|
|Disease class||Predicted group membership|
|Disease class||Predicted group membership||Total|
|Disease class||Predicted group membership||Total|
Table 2 displays the standardized canonical discriminant coefficients for all three classes. These coefficients are semipartial coefficients that compare the relative importance of the variable to the independent variable. They are analogous to the β weights of regression analysis . In this table, the variables with positive correlation coefficients are listed in descending order, whereas those with negative coefficients are listed in ascending order. This is done to reflect that the magnitude conveys the quantity of the correlation, and the sign indicates the quality (positive or negative). An inspection of this table reveals that only a small number of taxa have a number higher than 1, indicating that their correlation with the disease is not strong.
|Diseased vs. Healthy||Diseased vs. Healthy‐on‐Diseased||Healthy vs. Healthy‐on‐Diseased|
Table 3 displays the structure coefficient table. The structure coefficients are full coefficients, meaning that they are pooled as within‐groups correlations between the independent variable and the standardized canonical discriminant coefficients . According to Klecka , the structure coefficient is the authoritative coefficient for determining the importance of the independent variable to the discriminant function. As with the standardized canonical discriminant coefficient, the absolute value represents the quantity and the sign indicates quality of the correlation. In the case of structure coefficients, however, the largest absolute is 1.0. If a taxon had a correlation close to +1, then it could be identified as a cause of the disease. An examination of the table reveals that
|Diseased vs. Healthy||Diseased vs. Healthy‐on‐Diseased||Healthy vs. Healthy‐on‐Diseased|
For the second comparison of the bacterial taxa on the cuticle of Diseased lobsters versus Healthy cuticle on Diseased lobsters, the functions at the group centroids were well separated with a value of 2.060 for disease class 1 (Diseased) and ‐3.541 for disease class 3 (Healthy‐on‐Diseased;
Although these are the taxa with the highest positive structure coefficients, they are considerably less than 1. This indicates that, although they do discriminate for the disease compared to Healthy‐on‐Diseased samples, they do so weakly.
For the third comparison of the bacterial taxa on the cuticle of Healthy lobsters versus Healthy cuticle on Diseased lobsters, the discriminant functions at the group centroids demonstrated a good separation with values of 2.439 for disease class 2 (Healthy) and ‐3.582 for disease class 3 (Healthy‐on‐Diseased;
3.4. Correlational network analysis
Figure 6 is the network diagram of bacteria from the lobsters identified as Diseased. The edges (connecting links) between the nodes (taxa) are either red (negative correlation) or blue (positive). The width of the edge corresponds to the magnitude of the correlation coefficient. An inspection of the map reveals that the genera
Figure 7 is the network diagram of bacteria from the lobsters identified as Healthy. In this diagram,
Figure 8 represents the Healthy‐on‐Diseased microbiome correlations. In this map, there are five unconnected clusters. The largest consists of four taxa in which all are positively correlated. The connection is an unbranched chain with
Comparing the three maps, the average clustering coefficient of the Diseased microbiome is 0.289, the Healthy microbiome has an average clustering coefficient of 0.337, and the Healthy‐on‐Diseased microbiome has an average clustering coefficient of 0.0. Clustering coefficients are the ratios of the actual number of connections in comparison to the total possible connections and range from 0 to 1 . In other words, the clustering coefficient is the probability that one node is connected to another.
Figure 9 is a correlational difference network diagram that compares Healthy and Diseased microbiomes. This map displays taxa whose correlation coefficient changes significantly from Healthy to Diseased class. They are color coded to reflect the nature of the change. Red edge color indicates that the taxa correlate positively with the Healthy state and negatively with the Diseased state. Blue indicates that they correlate with both states positively but have a larger positive correlation coefficient in the Healthy state. Green indicates a negative correlation with both states but a larger negative correlation coefficient with the Healthy state.
The network consists of three disconnected groups of taxa.
This network is more complex, with a large group consisting of 10 taxa, 3 dyads, and 1 triad, suggesting that a number of potentially pathogenic processes or interactions are different in the Healthy‐on‐Diseased biofilm that do not occur in the Healthy biofilm. Within the large group,
The triad, similar to the triad in the Healthy versus Diseased correlational difference map, forms an incomplete clique , with
4. Discussion and conclusion
Smolowitz et al.  found that, although more advanced cases of ESD presented with extensive lesions throughout the dorsal surface, most cases were restricted to the cephalothorax.
The bacterial diversity as measured by the Shannon index appears to be roughly similar on all regions, with the claw and cephalothorax regions having the most similar indices (see Figure 1). Smolowitz et al.'s observations that ESD lesions were more prevalent on the upper cephalothorax, however, influenced the decision that a further molecular study would be focused on this region of the carapace. Furthermore, as the taxa distribution is relatively similar in all regions, we would be unlikely to miss key taxa by focusing on the cephalothorax.
As mentioned previously, the inset histogram (Figure 2) represents 170 taxa that were identified and at a normalized abundance greater than 1%. The most abundant taxa include the genera
These taxa do conform, however, to what has been observed about ESD. There is some, but not extensive, degradation of chitin. In these six genera, only one genus (
Bell et al.  found that ectohydrolase activity differed within the biofilms of lobsters from the geographic areas with high ESD incidence compared to lobsters inhabiting areas of little to no ESD. This could indicate that there is a shift in metabolic activity of the bacteria that are present or could indicate a shift in population.
The neighbor joining tree (Figure 3) of
Researchers have implicated a novel species (
An inspection of the UniFrac data that are depicted by the tree (Figure 4) reveals that, although there is some clustering, the lobsters identified as Healthy and Diseased show enough similarity that they exist as neighbors on all of the major branches of the tree. There is no statistically significant difference between the three classes as demonstrated by the enumeration of their
The PCO reveals much the same pattern; although there are identifiable regions, there is an overlap between the classes.
Our original simple hypothesis is that there is a significant difference between the bacterial communities in the Healthy and Disease states, indicating that this difference is associated with the etiology of the disease. The null hypothesis asserts that the bacterial communities are similar on the Healthy and Diseased lobsters. The phylogenetic metrics employed to compare the microbiomes of subjects in this study indicate that there is no significant difference between the three classes, which leads us to accept the null hypothesis and reject the original hypothesis. However, these data suggest that, although there may not be a gross difference between the microbiomes of Healthy and Diseased lobsters, there is a subtle shift that requires further investigation and the definition of a more complex hypothesis. A major avenue of inquiry would be to determine whether this shift is one of the metabolic activities among essentially the same bacteria, if there is an emergence of some members of the population that displace other extant members, or if there are new recruits to the community from outside the biofilm. Wahl et al.  suggested the use of emergent techniques, such as desorption electrospray ionization‐mass spectrometry, as a means of surveying the surface of the biofilm to elucidate the metabolic compounds that are present in biofilms. Coupled with confocal microscopy techniques, as suggested by Costerton , researchers could begin to develop a spatial representation of the lobster surface microbiome that would include both the identity of the bacterial cells and the metabolic activity that was occurring in their vicinity.
From the perspective of disease etiology, this kind of shift may point to a cause other than bacterial. If the same microflora, albeit in different abundances, are found on lobsters that have the disease as well as on those that are apparently healthy, then perhaps this dysbiosis represents a sign of some kind of disorder that is caused by something other than bacteria. Laufer et al.  established a correlation between lobsters with ESD and high tissue concentrations of alkylphenols. As noted earlier, these compounds retard the cross‐linking of tyrosine moieties in the carapace, which makes shell hardening take longer. Tarrant et al.  found evidence to suggest that gene expression in lobsters affected by ESD might indicate increased exposure to xenobiotics. Homerding et al.  found that lobsters living in the geographic area with a high incidence of ESD presented with reduced immunocompetence. This evidence supports an assertion that ESD is not primarily a bacterial disease but that bacterial lesions are a manifestation of a systemic dysfunction.
Kunkel et al.  discovered that lobsters with ESD had a loss of calcification in and around the disease lesions, particularly in the trabecular structures that are composed of apatite. One question that needs to be investigated further is whether the decalcification is related to the presence of alkylphenols. Could exposure of the underlying structure to the environment lead to decalcification? Alternatively, could acidic metabolic products of bacterial action produce the same results?
To use an analogy, suppose there is a house whose roof is defective. Imagine the nails were weak and could not hold the protective components of the roof in place, similar to the protein‐chitin structures in the lobster carapace in which the tyrosine cross‐linking was not complete. Rainwater would leak into the underlying structure, analogous to seawater getting into the trabecular apatite of the lobster carapace. In both cases, this would physically weaken the underlying structure. In both cases, normal microflora would have access to parts of the structure from which they would normally be excluded and may start to use these energy and carbon sources. This would further weaken and disfigure the structure and increase the likelihood of a breach in the outer protective structure. In this analogy, it is not the microflora that causes the problem in the roof. Rather, the rot that would ensue from its invasion is the result of the ineffectiveness of barriers whose function is to keep these organisms out of the underlying structure. In the case of ESD, perhaps the bacterial lesions are not the cause of the disease but merely the inevitable result of structural weaknesses in the lobster carapace. Additionally, the trabecular apatite and the other calcium moieties present in the lobster cuticle may have antimicrobial properties themselves . Their physical dissolution could be the result of the intrusion of acidified seawater, as there is now evidence that ocean acidification is affecting organisms that produce calcareous shells . The absence calcium moieties from the site of lesions could represent an additional reduction in the ability to exclude microbiota from the exoskeleton.
Although there have been other studies of biofilms related to crustaceans [50, 51], this represents the first survey of the surface microbiome of the lobster that uses high‐throughput, culture‐independent molecular techniques. As mentioned above, it is by no means an exhaustive analysis of the lobster microbiome, as the subjects of this study are geographically restricted and are not a large enough size to yield a global generalization. It is comprehensive enough, however, to demonstrate that, although there is no gross difference between the microbiomes of lobsters with and without ESD in Rhode Island samples, there may be a subtle shift in the microbial population that correlates with dysfunction. As mentioned previously, Bell et al.'s work indicates a variance in the metabolic activity of bacteria in eastern Long Island Sound, an area of high ESD occurrence, compared to those found on the surface of lobsters in regions that have little to no ESD incidence . This shift may be more closely related to a shift in activity rather than of population.
Considering the first DA analysis (Diseased vs. Healthy), the canonical correlation coefficient (
This analysis supports the affirmative hypothesis that there is a difference between the microbial taxa in shell‐diseased lobsters compared to healthy lobsters. The study identifies 58 bacteria that are significantly different between the two classes and rejects 112 that do not contribute to discrimination between them. The relatively low correlation coefficients for these 58 taxa, however, do indicate that they are present on all three classes but in slightly different abundances.
This analysis also provides additional evidence that ESD in the American lobster correlates with dysbiosis rather than the presence of a single overt pathogen. In a disease caused by a single pathogen, one would expect a structure coefficient of one of the variables to approach 1. An inspection of the structure coefficients reveals that none of the bacterial taxa have a coefficient greater than 0.325.
The role of
This analysis supports the observations that were articulated above. The weak correlations of some of the taxa with the three states of health reinforce the assertion that, although there are differences in the abundances of these taxa when comparing these three states, none of them are exclusive to one or the other state. The weak correlations also demonstrate that this is not a disease caused by a discrete pathogen, which would have a much stronger correlation. The works of the author's colleagues within the New England Lobster Health Initiative strongly suggest that there are underlying causes that render some lobsters, particularly within the waters of eastern Long Island Sound and Rhode Island, susceptible to an opportunistic infection [10, 14–16, 21]. Exposure to abnormally high levels of endocrine‐disrupting compounds such as alkylphenols appears to fit the profile of intervening agents, which is described by Laufer et al. and Tarrant et al.; not only are they capable of eliciting a response consistent with an endocrine‐disrupting compound, but they also appear to retard proper shell development. Hitherto, there is a correlation, but no direct evidence, linking alkylphenol exposure to ESD. Obtaining such evidence would require an experiment in which a group of lobsters was exposed to alkylphenols, whereas a control group would be kept in an alkylphenol‐free environment. As discussed previously, however, captive lobsters are more susceptible to impoundment shell diseases . Designing an experiment that would mitigate such confounding results is likely to prove challenging.
The abundance tables, as well as the DA computations reported in the previous chapter, reveal that
The correlational difference networks may shed some light on the role of
The only data used in these correlational networks were those derived from MTPS of sample lobsters from the “100 Lobsters” Project . As described in Section 1, samples from these lobsters were used in other investigations, and the data obtained have been recorded [14, 16, 20]. Integrating these data by means of correlational network analysis with alkylphenol concentrations in the tissues and ectohydrolases on the surface of the same lobster, for example, may produce a more definitive correlation between ESD and alkylphenols and may allow researchers to correlate which bacterial taxa are involved in what metabolic activity.
When comparing the data obtained using previous methods, a slightly different picture emerges. Taxa such as
As stated above, the hypotheses addressed are as follows:
H: The bacterial communities on lobsters with ESD are significantly different in quality and quantity than unaffected lobsters.
HA: The bacterial communities on the carapaces of healthy and shell‐diseased lobsters are similar.
There are no major shifts in the microbiome between the Healthy and Diseased states, so we accept the alternative hypothesis and reject the original simple hypothesis. However, there are minor differences in the microbiomes of Healthy and Diseased lobsters. Furthermore, the microbiomes of lesion‐free carapace on Diseased lobsters exhibited different microbiome compositions. However, the microbiomes of all three classes of the Diseased state that were identified here were similar enough that they occupied the same branches on a weighted UniFrac tree, and the difference between the microbiomes was determined to be statistically insignificant. Although there appears from these data to be a subtle difference between the microbiomes of samples from these three sample classes, the difference is not enough to fully support the original simple hypothesis. We must move forward with a more complicated hypothesis where environmental factors play a major role in the etiology of the disease. As such, we can start to define ESD as a complex environmental disease.
Arguably one of the most important advances in health and disease prevention has been the recognition that multiple pathogens are the cause of many diseases and it is now well accepted that these “polymicrobial disease” are in fact quite common . As the body of knowledge has increased and we have extended our limits of detection and sequencing throughput, a growing number of diseases and syndromes have emerged that have been shown to be polymicrobial in nature. ESD is an example of one such polymicrobial disease.
This research has employed culture‐independent techniques coupled with multivariate statistical treatment of the resultant data that present a shift in the lobster microbiome that correlates with the disease. This phenomenon is defined here as a dysbiosis. Although this research has elucidated this subtle shift in the microbiome of the lobster, it has not addressed the etiology of the disease. Indeed, it is not clear that the bacterial manifestations of this disease are anything more than a proximal cause. Researchers have successfully induced a condition that resembles ESD in captive lobsters but under extraordinary conditions. The evidence of involvement of alkylphenol contamination, for example, is circumstantial at this point. A direct evidence might be obtained through controlled experimentation as long as the confounding effects of captivity can eliminated.
The bulk method of extracting bacterial DNA from the lobster samples is itself confounding to the process of understanding the disease. Removing the biofilm and extracting the microbiomic DNA in bulk fashion precludes the interrogation of the spatial aspect of the biofilm. If microbes are cooperatively harvesting material and energy from a site, understanding their positional relationship may be useful. In addition, the nature of bulk extraction is that a relatively large amount of microbiota is captured. Such an approach can obscure members of the microbial community whose importance to the ecological function is disproportional to their relative abundance. It may be that minority members of the community, when observed in their spatial context within the biofilm, play a pivotal role in metabolic function. In the research described here, these bacteria may have been discarded due to insufficient abundance.
As noted previously, this research did not subdivide the disease lesions based on severity. Doing so would have yielded a more complete picture of the disease process and may have provided us with more insight into the roles of
In response to the ESD crisis, a group of scientists, fisheries managers, and lobstermen formed the New England Lobster Health Initiative. It was this ad hoc committee that ultimately received funding that was awarded to several university and institutional researchers, including the award that funded what is reported here. The findings of each group were published as a special edition of the