Microbial Consortium Associated with Crustacean Shells Composting

Soil microbes play an essential role in the biodegradation of crustacean shells, which is the process of sustainable bioconversion to chitin derivatives ultimately resulting in the promotion of plant growth properties. While a number of microorganisms with chitinolytic properties have been characterized, little is known about the microbial taxa that participate in this process either by active chitin degradation or by facilitation of this activity through nutritional cooperation and composting with the chitinolytic microorganisms. In this study, we evaluated the transformation of the soil microbiome triggered by close approximation to the green crab shell surface. Our data indicate that the microbial community associated with green crab shell matter undergoes significant specialized changes, which was reflected in a decreased fungal and bacterial Shannon diversity and evenness and in a dramatic alteration in the community composition. The relative abundance of several bacterial and fungal genera including bacteria Flavobacterium, Clostridium, Pseudomonas, and Sanguibacter and fungi Mortierella, Mycochlamys, and Talaromyces were increased with approximation to the shell surface. Association with the shell triggered significant changes in microbial cooperation that incorporate microorganisms that were previously reported to be involved in chitin degradation as well as ones with no reported chitinolytic activity. Our study indicates that the biodegradation of crab shells in soil incorporates a consortium of microorganisms that might provide a more efficient way for bioconversion.


Introduction
Globally, around 600,000 tons of chitin waste are generated by the seafood industry annually [1], which brings substantial challenges for waste disposal; the disposal of one ton of shellfish waste through disposal outlets can cost up to USD 150 [2]. The routes on the crustacea waste disposal require a third party or company to take the material away from the possessing site for disposal by methods including aerobic and aerobic digestion, composting, and land spreading. Nevertheless, the reuse of waste from the industry is not a common practice, and a significant proportion of crustacean waste is deposited into the environment, generating potential environmental hazards due to fish tissue deterioration. Importantly, crustacean shells contain chitin, proteins, calcium carbonate, and carotenoids, which are valuable resources for pharmaceutical, agricultural, construction, and paper industries [1][2][3][4][5]. For example, treatments with chitin and shell derivatives have been shown to increase soil suppressiveness against plant fungal and nematode infections [6][7][8][9][10] and therefore attract significant attention as potential environmentally friendly approaches to control plant pathogens. Additionally, shell chitin derivatives are applied in agriculture to improve plant growth by influencing plant growth-promoting microorganisms, stimulating plant metabolism [11,12], and the induction of defense-related genes [7,13]. collected in the separate plastic bags. The spatula was cleaned between different sample collections to avoid any cross-contamination. A quantity of 15 control soil (control) samples from the adjacent location of the collection site were also collected at the same time. All the collected samples were then placed on ice and transported to the laboratory for processing. After transportation to the lab, 5 g of the soil samples was sieved (2 mm) and immediately stored at −86 • C until processing for DNA isolation. The shells were removed from the bags, vigorously shaken, frozen in liquid nitrogen, ground into a fine powder using a sterile pestle, and stored at −86 • C.

DNA Extraction and Sequencing
Approximately 250 µg of soil or shell tissue per sample was used to isolate fungal DNA using the QIAGEN Power Soil DNA extraction kit (Cat No.12888-100) following the manufacturer's protocols. At least 50 ng (10 µL) of DNA sample were sent to the Dalhousie University CGEB-IMR (https://imr.bio (accessed on 2 March 2022)) for V6-V8 16S rRNA gene (16S; forward: ACGCGHNRAACCTTACC; reverse: ACGGGCRGTGWGTRCAA) and fungal ITS2 region (ITS; forward: GTGAATCATCGAATCTTTGAA; reverse: TCCTC-CGCTTATTGATATGC) library preparation and sequencing. Samples were multiplexed using a dual-indexing approach and sequenced using an Illumina MiSeq with paired-end 300 + 300 bp reads. All PCR procedures and Illumina sequencing details were as previously described [29,30]. All sequences generated in this study are available in the NCBI sequence read archive under the accession numbers PRJNA835461 and PRJNA835468.

Sequence Processing
The overlapping paired-end forward and reverse reads were stitched together using PEAR [31] and exported into QIIME2 [32]. The sequences were trimmed of their primers using QIIME2's Cutadept plug-in [33,34]. Low-quality sequences were filtered from the dataset using QIIME2's q-score-joined function. Using QIIME2's Deblur plug-in, the sequences were organized into amplicon sequence variants (ASVs)-high resolution genomic groupings [33,35,36]. In order to account for the potential MiSeq bleed-through between runs (estimated by Illumina to be less than 0.1%), ASVs that accounted for less than 0.1% of the total sequences were removed. Taxonomic classifications were assigned to the ASV using QIIME2's naïve-Bayes scikit-learn function, referencing SILVA databases [37,38]. Additionally, ASV assigned to mitochondria and chloroplasts were filtered out [33]. Sequencing of one 16S soil sample failed. After filtering unclassified and plant-derived ASVs, a total of 413,496 and 732,021 high-quality 16S and ITS reads were obtained from 44 and 45 samples, respectively. These reads were distributed across 4915 and 1003 16S and ITS ASVs, respectively. To assess microbial diversity, the sets were normalized to the depth of 3123 16S and 3758 ITS reads per sample, resulting in the identification of 4901 bacterial and 1001 fungal ASVs. During the normalization process, 2 ITS control and 4 16S soil samples were removed from analysis.

Bioinformatics and Statistical Analysis
QIIME2's diversity function was used to calculate Shannon indices (alpha diversity) as well as UniFrac matrices (beta diversity) [39,40]. These UniFrac matrices were then subjected to an ADONIS test to determine what proportion of variance in the community structure could be attributed to treatment. Nonmetric multidimensional scaling (NMDS) of bacterial communities was performed on Bray-Curtis matrices using the Vegan R package [41]. Differential abundances in bacterial taxa were determined using ALDEx2 [42] with the Benjamini-Hochberg-corrected p-value of the Kruskal-Wallace test (p < 0.05). The graphics were produced using ggplot2 [43]. Core ASVs were identified using the QI-IME2 core-features plugin. The co-occurrence analysis was performed using the CCREPE (Compositionality Corrected by REnormalization and PErmutation) R package [44] with 1000 bootstrap iterations and default settings. To obtain comparable datasets from each treatment, 10 replicate samples from the Control, Soil, and Shell dataset were randomly selected. The co-occurrence and co-exclusion patterns in the samples were scored. The results were filtered to remove nonstatistically significant relationships. We generated the network based on correlations with p-values < 0.05. The networks were analyzed with Cytoscape [45]. Indicator species (p-values < 0.001, computed using 5000 permutations) were determined by multi-level pattern analysis using R package "Indicspecies" [46].

The Effect of Shell Proximity on the Diversity of Microbial Communities
We observed profound changes in microbial alpha-and beta-diversity with the increase in shell proximity. The niche, Control vs. Soil vs. Shell, explained 70% and 66% of bacterial and fungal communities' variation, respectively (Table S1; Figure 1). Substantial variation was also detected between bacterial and fungal communities from Control vs. Soil groups (75% and 58%, respectively), Control vs. Shell samples (68% and 69%, respectively), and Soil vs. Shell samples (42% and 50%, respectively) (Table S1; Figure 1). The Shell microbiome exhibited a significant decrease in both bacterial and fungal Shannon diversity and Evenness, compared to the Control and the Soil microbiomes and bacterial Shannon diversity, which had a significant decrease in the Soil compared to the Control microbiome ( Figure 2). graphics were produced using ggplot2 [43]. Core ASVs were identified using the QIIME2 core-features plugin. The co-occurrence analysis was performed using the CCREPE (Compositionality Corrected by REnormalization and PErmutation) R package [44] with 1000 bootstrap iterations and default settings. To obtain comparable datasets from each treatment, 10 replicate samples from the Control, Soil, and Shell dataset were randomly selected. The co-occurrence and co-exclusion patterns in the samples were scored. The results were filtered to remove nonstatistically significant relationships. We generated the network based on correlations with p-values < 0.05. The networks were analyzed with Cytoscape [45]. Indicator species (p-values < 0.001, computed using 5000 permutations) were determined by multi-level pattern analysis using R package "Indicspecies" [46].

The Effect of Shell Proximity on the Diversity of Microbial Communities
We observed profound changes in microbial alpha-and beta-diversity with the increase in shell proximity. The niche, Control vs. Soil vs. Shell, explained 70% and 66% of bacterial and fungal communities' variation, respectively (Table S1; Figure 1). Substantial variation was also detected between bacterial and fungal communities from Control vs. Soil groups (75% and 58%, respectively), Control vs. Shell samples (68% and 69%, respectively), and Soil vs. Shell samples (42% and 50%, respectively) (Table S1; Figure 1). The Shell microbiome exhibited a significant decrease in both bacterial and fungal Shannon diversity and Evenness, compared to the Control and the Soil microbiomes and bacterial Shannon diversity, which had a significant decrease in the Soil compared to the Control microbiome ( Figure 2). The difference between communities based on Bray-Curtis distance Adonis tests was used to assess whether beta-diversity is related to sample groupings, 999 permutations, R2, *** p < 0.001. The difference between communities based on Bray-Curtis distance Adonis tests was used to assess whether beta-diversity is related to sample groupings, 999 permutations, R2, *** p < 0.001.

The Effect of Shell Proximity on the Composition of Microbial Communities
The proximity of the shell had a dramatic effect on microbiome composition. A total of 55 bacterial and 17 fungal classes and 401 bacterial and 131 fungal genera were differentially represented between niches. Highly abundant fungal classes Mortierellomycetes, and Pezizomycetes and bacterial classes Bacteroidia, Gammaproteobacteria, Actinobacteria, Clostridia, and Bacilli were significantly enriched in the Shell microbiome. Considering highly abundant genera, bacteria Flavobacterium, Clostridium, Pseudomonas, Sanguibacter, Tissierella, and Streptomyces and fungi Mortierella, Scutellinia, and Paracremonium were enriched in the Shell microbiome ( Figure 3).
The fungal community showed a greater overlap (>21%) in the ASVs between all three niches (Control, Soil, and Shell) compared with the bacterial community (<10%)

The Effect of Shell Proximity on the Composition of Microbial Communities
The proximity of the shell had a dramatic effect on microbiome composition. A total of 55 bacterial and 17 fungal classes and 401 bacterial and 131 fungal genera were differentially represented between niches. Highly abundant fungal classes Mortierellomycetes, and Pezizomycetes and bacterial classes Bacteroidia, Gammaproteobacteria, Actinobacteria, Clostridia, and Bacilli were significantly enriched in the Shell microbiome. Considering highly abundant genera, bacteria Flavobacterium, Clostridium, Pseudomonas, Sanguibacter, Tissierella, and Streptomyces and fungi Mortierella, Scutellinia, and Paracremonium were enriched in the Shell microbiome ( Figure 3).
Microorganisms 2022, 10, x FOR PEER REVIEW 6 of 16 ( Figure 4A,B). The Control and the Shell microbiomes had minor overlap (<2%). The Soil microbiome shared more ASVs with the Shell compared to the Control microbiome-<10% of bacteria and <7% of fungal ASVs were shared between the Soil and Control, while >29% of both fungal and bacterial ASVs were shared between the Soil and Shell microbiomes. For both bacterial and fungal communities, the greatest proportion of specialists was found in the Control microbiome (41% and 39%, respectively). The shell microbiome specialists were represented by 6% and 3% of bacterial and fungal ASVs, respectively, and Soil microbiome specialists were represented by 3% and 10% of bacterial and fungal ASVs, respectively ( Figure 4A,B).

Core Microbiomes
Next, we looked at core microbiomes in the Control, Soil, and Shell microbiomes, which was defined as ASVs found in all but one sample from each niche, which represented at least 92.5% of all samples in each niche (Supplementary Tables S4-S9). In general, both bacterial and fungal core microbiomes from the Control microbiome contained less ASVs compared to those from the Shell and Soil microbiomes. The bacterial Shell, Soil, and Control core microbiomes were composed of 69, 117, and 48 ASVs, respectively, and fungal core microbiomes contained 58, 87, and 29 ASVs, respectively ( Figure 4C,D). The Shell and Soil core microbiomes shared the largest proportion of ASVs (17% of bacterial and 38% of fungal ASVs, respectively) compared to those shared between Soil and Control core microbiomes (3% of bacterial and 2% of fungal ASVs, respectively) ( Figure 4C,D). In addition, 5% of fungal ASVs were shared between all three niches and no bacterial ASVs were common for these niches. Additionally, no ASVs were common between bacterial and fungal Control and Shell core microbiomes.
When the relative abundances of core ASVs were combined together at the family taxonomic level, we found that the Shell core was dominated by a few bacterial and fungal families ( Figure 5) such as bacteria Clostridiaceae, Sanguibacteraceae, Pseudomonadaceae, Streptomycetaceae, Peptostreptococcales, and Flavobacteriaceae, and fungi Mortierellaceae, Trichocomaceae, Microascaceae, Pseudeurotiaceae, Piskurozymaceae, and Nectriaceae ( Figure 5). This families contained several taxa overrepresented in the overall Shell microbiome ( Figure 3). Together, these families comprised 61% and 82% of the total 16S and ITS Shell core reads, respectively. These families were also found in Soil core microbiomes and were represented by 22% and 63% of total 16S and ITS Soil core reads, respectively. On the other hand, only two of these bacterial (Clostridiaceae and Streptomycetaceae) and fungal (Mortierellaceae and Nectriaceae) taxa were found in the Control core microbiome, and they were represented by 2% and 10% of total 16S and ITS Control core reads, respectively.

Core Microbiomes
Next, we looked at core microbiomes in the Control, Soil, and Shell microbiomes, which was defined as ASVs found in all but one sample from each niche, which represented at least 92.5% of all samples in each niche (Supplementary Tables S4-S9). In general, both bacterial and fungal core microbiomes from the Control microbiome contained less ASVs compared to those from the Shell and Soil microbiomes. The bacterial Shell, Soil, and Control core microbiomes were composed of 69, 117, and 48 ASVs, respectively, and fungal core microbiomes contained 58, 87, and 29 ASVs, respectively ( Figure 4C,D). The Shell and Soil core microbiomes shared the largest proportion of ASVs (17% of bacterial and 38% of fungal ASVs, respectively) compared to those shared between Soil and Control core microbiomes (3% of bacterial and 2% of fungal ASVs, respectively) ( Figure  4C,D). In addition, 5% of fungal ASVs were shared between all three niches and no bacterial ASVs were common for these niches. Additionally, no ASVs were common between bacterial and fungal Control and Shell core microbiomes.
When the relative abundances of core ASVs were combined together at the family taxonomic level, we found that the Shell core was dominated by a few bacterial and fungal families ( Figure 5) such as bacteria Clostridiaceae, Sanguibacteraceae, Pseudomonadaceae, Streptomycetaceae, Peptostreptococcales, and Flavobacteriaceae, and fungi Mortierellaceae, Trichocomaceae, Microascaceae, Pseudeurotiaceae, Piskurozymaceae, and Nectriaceae ( Figure 5). This families contained several taxa overrepresented in the overall Shell microbiome (Figure 3). Together, these families comprised 61% and 82% of the total 16S and ITS Shell core reads, respectively. These families were also found in Soil core microbiomes and were represented by 22% and 63% of total 16S and ITS Soil core reads, respectively. On the other hand, only two of these bacterial (Clostridiaceae and Streptomycetaceae) and fungal (Mortierellaceae and Nectriaceae) taxa were found in the Control core microbiome, and they were represented by 2% and 10% of total 16S and ITS Control core reads, respectively.

Microbiome Cooperation
We generated a co-occurrence network by correlating relative abundances between bacterial and fungal ASVs from the Control, Soil, and Shell microbiomes. The list of ASVs found in each network and their characteristics can be found in Supplementary Tables S10-S12. The Shell and Soil network exhibited less cooperation compared to the Control network. The Shell co-occurrence network contained 214 ASVs, 884 interactions, and an average of 8.262 neighbor ASVs, and the Soil co-occurrence network contained 348 ASVs, 1606 interactions, and an average of 9.230 neighbors, while the Control network interaction contained 307 ASVs, 2689 interactions, and an average of 17.518 neighbors. We found a significant overlap between the Shell and Soil network-more than half of the ASVs in the Shell network were also part of the Soil network (116 ASVs) and 12 ASVs were common to all three networks. Several most connected in the Shell network ASVs also found in the Soil network included Dongia (30 and 12 interaction), Rhizobiales (29 and 7 interaction), 2 Tissierella ASVs (29/22 and 7/4 interaction), Mortierella hypsicladia (28 and 9 interaction); Scutellinia vitreola (27 and 10 interaction), and Mortierella (24 and 7 interaction) in the Shell and Soil network, respectively. On the other hand, the Control network did not have much commonality with the Shell and Soil networks: it had no unique overlap with the Shell and had only 28 ASVs of unique overlap with the Soil network ( Figure 4E).

16S and ITS2 ASVs with Strong Association with Shell Environment
To identify the most influential taxonomic groups associated with the shell environments, we selected the Shell core ASVs with increased relative abundances in the Shell compared to the Control microbiome and filtered out the ASVs that were not part of the Shell network interaction (Supplementary Tables S13 and S14). These selected ASVs did not belong to the Control core microbiome and network interaction. We also analyzed the strength and statistical significance of the relationship between the ASVs occurrence/abundance and their association with the specific niche and verified that all selected ASVs were indicator species of the Shell/Soil microbiome. As a result, 60 bacterial and 44 fungal ASVs comprising 38% and 56% of Shell microbiome 16S and ITS reads, respectively, were selected. These ASVs were collapsed to 45 bacterial and 27 fungal genera or taxa with lower annotation levels above the genus ( Table 1). The most abundant bacterial taxa included Actinobacteria, Gamma-proteobacteria Pseudomonas, Bacteroidia Flavobacterium and Chryseolinea, Clostridia Tissierella and Ruminococcaceae, Alpha-proteobacteria Dongia, and Bacilli Vagococcus. The most abundant fungal taxa included Mortierellomycetes Mortierella, Eurotiomycetes Talaromyces, Sordariomycetes Mycochlamys, Fusicolla, Sordariales and Zopfiella, Leotiomycetes Pseudeurotium and Pseudogymnoascus, and Tremellomycetes Apiotrichum. Together, they were represented by 25% and 50% of 16S and ITS reads in the Shell microbiome, respectively (Table 1). Fungi Mortierella were represented by 31% of Shell microbiome ITS reads and contained four species, M. hyaline (19%), M. zonata (3%), M. hypsicladia (0.9%), and M. humilis (0.6%) and annotated Mortierella (7%) (Supplementary Tables S13 and S14).

Discussion
A deep understanding of the process underlying the biological degradation of crustacean's shell waste is an important step to maximize the utilization of its valuable components into high-value-added products and to facilitate environmentally friendly waste disposal. In this study, we explored microbial communities tightly associated with green crab shells composted in the soil for over a one-year period to identify the microbial complex of potential chitin degraders. To gain an overall understanding of the microbial community involved in shell composting in soil, we did not separate between the microbiomes located in microporosities, inside shell matter, in the shell surface, and in the biofilm state. We detected profound changes in microbial structure and composition with the increase in shell proximity. Both bacterial and fungal Shannon diversity and Evenness were significantly decreased in the Shell microbiome compared to adjacent (Soil) and dis-tantly located (Control) soils. We detected a small but statistically significant decrease in bacterial Shannon diversity in the Soil compared to the Control. The decrease in the Shell microbiome alpha-diversity might reflect its strong specialization and enrichment with a relatively small set of taxa capable of obtaining nutrients from chitin matter or feed on chitin degradation products. Our results are consistent with a previous study that reported a decrease in diversity indices in the soils treated with chitin-rich amendments after 3 weeks post-incubation [8]. On the other hand, another ecological study reported an increase in the bacterial diversity after chitin enrichment after 35 days post-inoculation [21], indicating that our understanding of the effects of chitin on soil biology is far from complete.
The changes in microbiome composition with the increase in shell proximity were reflected in the differential representation of microbial taxa between niches. Our analysis identified a large number of ASVs differentially represented between the Control and Shell microbiomes. Interestingly, 47% and 63% of total 16S and ITS reads, respectively, represented ASVs with increased relative abundance in the Shell compared to Control microbiome, indicating that the adaptation of the microbiome to shell environment involved major microbial taxa. For example, highly relatively abundant bacterial genera Flavobacterium, Clostridium, Pseudomonas, Sanguibacter, Tissierella, Streptomyces, Mortierella, Scutellinia, and Paracremonium were enriched in the Shell microbiome. Previously, it was shown that Flavobacterium, Pseudomonas, Sanguibacter, Streptomyces, and fungi Mortierella and Aspergillus were representative genera of chitinolytic microbes [13,[20][21][22]47,48]. On the other hand, to the best of our knowledge, there are no reports indicating the effect of chitin amendments on the abundances of bacteria Tissierella and fungi Scutellinia and Paracremonium in soil. It is possible that these microorganisms might not possess any chitinolytic activities but could take an advantage of chitin degradation products available in the shell environment. Another possibility is that the previous studies were looking at microorganisms associated with chitin-amended soil, while our study was focused on the shell matter with very different physical and chemical properties, which might have provided additional selective pressure for enrichment of chitinolytic microorganisms. The enrichment in Clostridium spp. might be linked to the activation of fermentation processes in the Shell microbiome. We also observed a significant overlap between the Shell and Soil microbiomes, and about 75% of 16S and ITS Soil ASVs were shared with the Shell microbiome. On the other hand, only 18% of 16S and 37% of ITS ASVs were shared between the Control and Shell microbiomes, indicating a significant influence of shell matter on the composition of the adjacent soil microbiome.
The above conclusion was also supported by the fact that 48% of 16S and 76% of ITS Shell core ASVs were shared with the Soil core, while only 1% of ITS Shell core ASVs were shared with the Control core microbiome and none of the 16S Shell core ASVs were found in the Control core. Additionally, the six most abundant bacterial and fungal classes found in the Shell core microbiome were represented by 61% and 82% of total 16S and ITS Shell core reads, respectively, reiterating our previous conclusion of the involvement of the major taxa in microbial adaptation to the shell environment.
Microbial cooperation was also affected by shell proximity, which was reflected in the decrease in microbial cooperation in the Shell and Soil compared to the Control network. This suggested a stronger specialization within the shell-associated microbiome involving a limited number of microbial taxa capable of directly feeding on the shell matter or to benefit from chitin degradation products available in the shell environment. We also detected a strong overlap between the Soil and the Shell networks. Around 60% of microbial ASVs from the Shell network were part of the Soil network, which might be a result of a significant enrichment of adjacent soil with byproducts of shell degradation stimulating a similar type of cooperation between microorganisms. However, we also should consider a possibility that small shell particles could be present in the Soil samples and therefore influence the diversity and structure of the Soil microbiome.
In our study, we applied the analyses of taxonomic composition, core microbiome, and network cooperation to identify ASVs potentially involved in shell degradation in soil. These analyses produced 3 groups of ASVs, which could be involved in this process based on a single defined criterion. However, some of these ASVs were identified as important based on more than one criterion, while others were only found in one of the groups. We speculated that the ASVs found in the overlap between all three groups could be considered as the most active participants in the shell degradation consortium. This approach allowed us to select 60 bacterial and 44 fungal ASVs comprising 45 bacterial and 27 fungal genera. These ASVs represented a tiny proportion of ASVs identified in the Shell microbiome (3% and 8% of 16S and ITS ASVs, respectively) but comprised a large proportion of 16S and ITS Shell microbiome reads (38% and 56%, respectively), again emphasizing the involvement of major ASVs in shell degradation.
The most active fungi in the shell degradation consortium included Mortierella (30%), Talaromyces (5%), Mycochlamys (4%), Fusicolla (2%), Sordariales (2%), and Pseudeurotium (2%). It has been reported that some Mortierella spp. have a chitinolytic activity [47] and their abundance is positively correlated with soil chitin [9]. In our study, we also detected a substantial increase in relative abundance of Mortierella in the Shell compared to the Control microbiome. Two most dominant Mortierella spp. found in the consortiums, M. hyaline (19%) and M. zonata (3%), were reported to possess a strong plant growthpromoting capability [61,62]. Additionally, Mortierella spp. and Flavobacteriaceae spp. have been correlated with disease suppression [63], emphasizing the potential of shell-fish waste-based amendments in the manipulation of the soil mycobiome to promote disease suppression. There are several reports on the isolation of chitinases from Talaromyces spp., a highly abundant genus found in the shell degradation consortium [64][65][66]. However, as far as we know, the direct involvement of genera Fusicolla, Mycochlamys, and Pseudeurotium in chitin degradation has not been documented. Nevertheless, there was a report indicating an increase in relative abundances of Mycochlamys and Pseudeurotium in chitin-treated soil [13], and Fusicolla has been identified as part of the sea-food fermentation microbiome [67]. Given the active involvement of these fungi in the shell degradation consortium and their high abundances in the Shell and Soil microbiomes, the additional study of the role on these taxa in chitin degradation should be examined in greater detail.

Conclusions
In the present study, we demonstrated profound changes in microbial structure and composition with the increase in shell proximity, which was reflected in the enrichment of the Shell microbiome with the subset of taxa with potential to obtain nutrients from chitin matter or feed on chitin degradation products. Our analysis identified a large number of microbial taxa overrepresented in the Shell compared to the Control microbiome. We applied a combination of several bioinformatic tools to narrow down this list of microorganisms to focus on the most active participants in the shell degradation consortium. This approach allowed us to identify a group of microbes that might play a central role in the degradation of shell matter in soil and might form tightly cooperated consortiums to facilitate nutrient acquisition from shell matter. In addition to a number of microorganisms previously reported to have chitin degradation activities, several genera without a known role in this process were identified in the consortiums. These taxa might play an active role in chitin degradation or take advantage of degradation products available in the shell environment and therefore facilitate chitin degradation. These newly identified microorganisms should be further studied for their potential application in environmentally friendly chitin degradation.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/microorganisms10051033/s1. Figure S1: Microbial classes identified in the study. A, bacterial 16S rRNA; B, fungal ITS; Figure S2: Microbial families identified in the core microbiomes. A, bacterial 16S rRNA; B, fungal ITS; Table S1: Variation in sample groupings as explained by weighted UniFrac dissimilarity distances. Adonis tests were used to assess whether beta-diversity is related to sample groupings, 999 permutations, R2, *** p < 0.001; Table S2: Bacterial ASVs that showed significant differences in relative abundance between Shell, Soil, and Control microbiomes based on ALDEx2. kw.ep-Expected p-value of Kruskal-Wallace test; kw.eBH-Expected Benjamini-Hochberg-corrected p-value of Kruskal-Wallace test; glm.ep-Expected p-value of glm test; glm.eBH-Expected Benjamini-Hochberg-corrected p-value of glm test. In bold are the ASVs found significantly overrepresented based on kw.eBH; Table S3: Fungal ASVs that showed significant differences in relative abundance between Shell, Soil, and Control microbiomes based on ALDEx2. kw.ep-Expected p-value of Kruskal-Wallace test; kw.eBH-Expected Benjamini-Hochbergcorrected p-value of Kruskal-Wallace test; glm.ep-Expected p-value of glm test; glm.eBH-Expected Benjamini-Hochberg-corrected p-value of glm test. In bold are the ASVs found significantly overrepresented based on kw.eBH; Table S4: Bacterial ASVs identified in core Shell microbiome; Table S5. Fungal ASVs identified in core Shell microbiome; Table S6: Bacterial ASVs identified in core Soil microbiome; Table S7: Fungal ASVs identified in core Soil microbiome; Table S8: Bacterial ASVs identified in core Control microbiome; Table S9: Fungal ASVs identified in core Control microbiome; Table S10: Node statistics for network representing cooperation in Shell microbiome; Table S11: Node statistics for network representing cooperation in Soil microbiome; Table S12: Node statistics for network representing cooperation in Control microbiome; Table S13: Bacterial ASVs with potential to be directly involved in shell degradation; Table S14: Fungal ASVs with potential to be directly involved in shell degradation.