Next Article in Journal
Enhanced Immunogenicity of Inactivated Dengue Vaccines by Novel Polysaccharide-Based Adjuvants in Mice
Next Article in Special Issue
Coriander (Coriandrum sativum L.) in Combination with Organic Amendments and Arbuscular Mycorrhizal Inoculation: An Efficient Option for the Phytomanagement of Trace Elements-Polluted Soils
Previous Article in Journal
Heterologous Expression, Purification, and Immunomodulatory Effects of Recombinant Lipoprotein GUDIV-103 Isolated from Ureaplasma diversum
Previous Article in Special Issue
Genomic, Antimicrobial, and Aphicidal Traits of Bacillus velezensis ATR2, and Its Biocontrol Potential against Ginger Rhizome Rot Disease Caused by Bacillus pumilus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Microbial Consortium Associated with Crustacean Shells Composting

1
USDA-ARS, Grain Legume Genetics and Physiology Research Unit, Prosser, WA 99350, USA
2
School of Science and the Environment, Grenfell Campus, Memorial University of Newfoundland and Labrador, Corner Brook, NL A2H 5G4, Canada
*
Author to whom correspondence should be addressed.
Microorganisms 2022, 10(5), 1033; https://doi.org/10.3390/microorganisms10051033
Submission received: 21 April 2022 / Revised: 10 May 2022 / Accepted: 12 May 2022 / Published: 16 May 2022

Abstract

:
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.

1. 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].
Chitin (polymer of (1→4)-β-linked N-acetyl-D-glucosamine) constitutes 15% to 40% of crustacean exoskeletons, representing around 75% of the shell organic fraction [14]. Industrial chitin extraction through chemical means has a number of limitations, including high energy requirements and the involvement of large amounts of sodium hydroxide and acidic treatments [15]. The biological extraction of chitin is a promising process, which can provide an alternative solution for decreasing production and environmental costs [1]. This process involves several steps, such as deproteinization, fermentation, and chitin-to-chitosan conversion, and could be facilitated by a number of microorganisms [7,16,17,18].
Culture-dependent and -independent analysis of chitinolytic microorganisms indicated that bacteria affiliated with Actinomyces, Proteobacteria, Flavobacteria, and Firmicutes and fungi Aspergillus and Mortierella are often identified as active chitin degraders [7,19,20,21]. More specifically, labeling studies pointed to a few key chitinolytic taxa, such as Pseudomonas, Massilia, and several families of Bacteroidetes as the most active in chitin degradations [22]. It was also suggested that chitinolytic bacteria sometimes process more chitin polymers than they are able to use themselves [23,24], providing the excess of the product to “satellite microbes”, opening the possibility for interspecific cross-feeding [19,22], a complex microbial cooperation within the chitin-degrading community. This in turn can improve the efficiency of chitin hydrolysis. The understanding of this already complex process is further complicated by the fact that chitin-degrading communities can undergo a temporal fluctuation [20] over a relatively short period of time: five to 90 days. Moreover, most of the studies aiming to identify chitinolytic microorganisms were focused on the soil bulk microbiome, while the analysis of those more tightly associated with shell microbes might provide additional information about the microbial complex involved in shell biodegradation. This approach could provide additional tools to facilitate environmentally friendly shell waste degradation.
The green crabs (Carcinus aestuarii) were first found in the 1950s in the eastern province of New Brunswick and extended to the southern part of Halifax in 1952. Since then, the population continued growing, and extended all over Newfoundland and Prince Edward Island [25]. The researchers worked on chitin extracted from green crab biomass to generate biodegradable polymers [26], and specialized chitin-based growing media had shown a reduction in soil pathogen levels leading to improved vegetable production [26]. The green crabs also had shells rich in calcium that can be used to neutralize acidic soils in Atlantic Canada. Similarly, the extracted chitin from crabs had a wide range of applications in Biomedicine, Pharmaceuticals, Food, Agriculture, and Personal Care Products [27]. The chitin-based compost could be used directly as a growing horticultural media to enhance soil pathogen control and significantly improve crop growth or as an extraction source for chitin/chitosan to be used for high-value applications in pharmaceuticals and the agricultural sector [27]. Previously, the crustacean composts were used in the nitrogen fixation, making the nitrogen available directly to the plants when added to the roots. The pathogens and microorganisms present in the crustaceans led to the degradation of the pathogens’ cell walls using the chitinases by composting [28].
The goal of this research was to assess the microbial communities, to evaluate the bacterial and fungal community tightly associated with green crab shell composted in Newfoundland, Canada, soils for a 1-year period, and to identify the microbial complex of potential chitin degraders.

2. Materials and Methods

2.1. Sampling Site Description and Sample Collection

The samples were collected on 14 October 2019, from Black Duck Siding, Western Newfoundland, Canada (48°34′00.1″ N 58°22′36.4″ W). The shells were originated from green crab processing facilities in Stephenville NL and were buried in soils for 1 year as part of the routine shell composting. A quantity of 15 crab shells (shell samples) were unearthed by a spatula from the soil at different locations on the same study site and packed in the labeled plastic bags, whereas the adherent soil (soil samples) to the crab shells was also 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.

2.2. 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: TCCTCCGCTTATTGATATGC) 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.

2.3. 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.

2.4. 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 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].

3. Results

3.1. Taxonomic Compositing

Overall, 79 bacterial and 25 fungal classes were detected in the dataset (Figure S1). In the total microbiome, Bacteroidia, Alphaproteobacteria, Gammaproteobacteria, Actinobacteria, and Clostridia were the dominant bacterial classes (18%, 16%, 14%, 13%, and 8%, respectively) and Mortierellomycetes, Sordariomycetes, and Leotiomycetes were the dominant fungi (24%, 22%, and 14%, respectively). The Control, Soil, and Shell microbiomes showed very different taxonomic compositions. The Control microbiome was dominated by unclassified Ascomycota, followed by Leotiomycetes, Dothideomycetes, and Sordariomycetes (28%, 21%, 13%, and 12%, respectively), while the Soil and Shell microbiomes were dominated by Sordariomycetes (33% and 20%, respectively), Mortierellomycetes (16% and 50%, respectively), Leotiomycetes (14% and 7%, respectively), and Eurotiomycetes (12% and 7%, respectively) (Figure S1).

3.2. 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).

3.3. 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).
Additionally, 430 bacterial and 173 fungal ASVs were differentially represented between the niches (Supplementary Tables S2 and S3), which represented 9% of 16S and 17% of ITS ASVs. These ASVs contained 55% and 75% of total 16S and ITS reads, respectively. The ASVs enriched in the Shell microbiome were represented by 47% and 63% of total 16S and ITS reads, respectively. Those highly abundant and overrepresented in the Shell microbiome ASVs were annotated as Clostridium (1 ASV, 6%), Sanguibacter (6%), Streptomyces (4%), Flavobacterium (3 ASVs: 3%, 2%, and 2%), Pseudomonas (2 ASVs: 2% and 2%), Mortierella (7%), M. hyaline (19%), M. hypsicladia (10%), M. macrospora (3%), M. hypsicladia (2 ASVs: 5% and 2%), M. zonata (3%), Pseudeurotium hygrophilum (2%), and Talaromyces atroroseus (2%). In total, these ASVs were represented by 27% and 51% to total 16S and ITS reads, respectively.
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%) (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).

3.4. 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.

3.5. 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).

3.6. 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).

4. 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 distantly 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.
Major bacterial classes identified by this approach included Actinobacteria, Gamma-proteobacteria, Bacteroidia, and Clostridia, represented by 13%, 8%, 5%, and 5% of Shell microbiome reads. Well-known chitinolytic bacteria Sanguibacter (6%) and Streptomyces (4%) [7,20,49] were the most abundant genera from class Actinobacteria. Gamma-proteobacteria Pseudomonas caeni and Bacteroidia Flavobacterium were represented by 5% and 2% of Shell microbiome reads, respectively. It was reported that some Pseudomonas and Flavobacterium spp. possess chitinolytic activities [22,47,50,51,52,53,54] and have been used for chitin extraction [55,56]. Highly abundant Clostridia Tissierella (2%) and Lactobacillales Vagococcus are known fermenters [57,58,59,60]. Tissierella was also identified as one of the most connected taxa in network interaction, which underlined the importance of this bacteria 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 growth-promoting 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.

5. 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–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 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.

Author Contributions

Conceptualization, S.N.Y. and M.C.; methodology, S.N.Y. and M.N.; software, S.N.Y.; validation, S.N.Y.; formal analysis, S.N.Y.; resources, S.N.Y. and M.C.; data curation, S.N.Y. and M.N.; writing—original draft preparation, S.N.Y.; writing—review and editing, S.N.Y., M.C. and M.N.; visualization, S.N.Y. supervision, S.N.Y. and M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a USDA ARS Project 2090-21000-002-00D and Discovery Grant from the Natural Sciences and Engineering Research Council of Canada to S.N.Y. Additional support was provided by the New Frontiers in Research Fund (NFRF) Transformation (Canada) to M.C.

Data Availability Statement

Data can be found within article and in supplementary files.

Acknowledgments

Authors greatly appreciate the support and materials provided by “Growing for Life” Stephenville, Newfoundland and Labrador and the assistance with English editing provided by Tomas Halford, English program, Memorial University of Newfoundland and Labrador.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Santos, V.P.; Marques, N.S.S.; Maia, P.; Lima, M.A.B.; Franco, L.O.; Campos-Takaki, G.M. Seafood Waste as Attractive Source of Chitin and Chitosan Production and Their Applications. Int. J. Mol. Sci. 2020, 21, 4290. [Google Scholar] [CrossRef] [PubMed]
  2. Yan, N.; Chen, X. Sustainability: Don’t waste seafood waste. Nature 2015, 524, 155–157. [Google Scholar] [CrossRef] [PubMed]
  3. Nguyen, T.T.; Barber, A.R.; Corbin, K.; Zhang, W. Lobster processing by-products as valuable bioresource of marine functional ingredients, nutraceuticals, and pharmaceuticals. Bioresour. Bioprocess. 2017, 4, 27. [Google Scholar] [CrossRef] [PubMed]
  4. Le, B.; Yang, S.H. Microbial chitinases: Properties, current state and biotechnological applications. World J. Microbiol. Biotechnol. 2019, 35, 144. [Google Scholar] [CrossRef] [PubMed]
  5. Ahuja, V.; Bhatt, A.K.; Sharma, V.; Rathour, R.K.; Rana, N.; Bhatia, R.K.; Varjani, S.; Kumar, M.; Magdouli, S.; Yung, H.; et al. Advances in glucosamine production from waste biomass and microbial fermentation technology and its applications. Biomass Convers. Biorefin. 2021, 1–23. [Google Scholar] [CrossRef]
  6. Radwan, M.A.; Farrag, S.A.A.; Abu-Elamayem, M.M.; Ahmed, N.S. Extraction, characterization, and nematicidal activity of chitin and chitosan derived from shrimp shell wastes. Biol. Fertil. Soils 2021, 48, 463–468. [Google Scholar] [CrossRef]
  7. Ilangumaran, G.; Stratton, G.; Ravichandran, S.; Shukla, P.S.; Potin, P.; Asiedu, S.; Prithiviraj, B. Microbial Degradation of Lobster Shells to Extract Chitin Derivatives for Plant Disease Management. Front. Microbiol. 2017, 8, 781. [Google Scholar] [CrossRef]
  8. Andreo-Jimenez, B.; Schilder, M.T.; Nijhuis, E.H.; Te Beest, D.E.; Bloem, J.; Visser, J.H.M.; van Os, G.; Brolsma, K.; de Boer, W.; Postma, J. Chitin- and Keratin-Rich Soil Amendments Suppress Rhizoctonia solani Disease via Changes to the Soil Microbial Community. Appl. Environ. Microbiol. 2021, 87, e00318-21. [Google Scholar] [CrossRef]
  9. Randall, T.E.; Fernandez-Bayo, J.D.; Harrold, D.R.; Achmon, Y.; Hestmark, K.V.; Gordon, T.R.; Stapleton, J.J.; Simmons, C.W.; VanderGheynst, J.S. Changes of Fusarium oxysporum f.sp. lactucae levels and soil microbial community during soil biosolarization using chitin as soil amendment. PLoS ONE 2020, 15, e0232662. [Google Scholar] [CrossRef]
  10. Cretoiu, M.S.; Korthals, G.W.; Visser, J.H.; van Elsas, J.D. Chitin amendment increases soil suppressiveness toward plant pathogens and modulates the actinobacterial and oxalobacteraceal communities in an experimental agricultural field. Appl. Environ. Microbiol. 2013, 79, 5291–5301. [Google Scholar] [CrossRef] [Green Version]
  11. El Hadrami, A.; Adam, L.R.; El Hadrami, I.; Daayf, F. Chitosan in plant protection. Mar. Drugs 2010, 8, 968–987. [Google Scholar] [CrossRef] [PubMed]
  12. Ramírez, M.Á.; Rodríguez, A.T.; Alfonso, L.; Peniche, C. Chitin and its derivatives as biopolymers with potential agricultural applications. Biotecnol. Apl. 2010, 27, 270–276. [Google Scholar]
  13. Debode, J.; De Tender, C.; Soltaninejad, S.; Van Malderghem, C.; Haegeman, A.; Van der Linden, I.; Cottyn, B.; Heyndrickx, M.; Maes, M. Chitin Mixed in Potting Soil Alters Lettuce Growth, the Survival of Zoonotic Bacteria on the Leaves and Associated Rhizosphere Microbiology. Front. Microbiol. 2016, 7, 565. [Google Scholar] [CrossRef] [PubMed]
  14. Fabritius, H.O.; Sachs, C.; Triguero, P.R.; Roobe, D. Influence of structural principles on the mechanics of a biological fiber-based composite material with hierarchical organization: The exoskeleton of the lobster Homarus americanus. Adv. Mater. 2009, 21, 391–400. [Google Scholar] [CrossRef]
  15. Rinaudo, M. Chitin and chitosan: Properties and applications. Prog. Polym. Sci. 2006, 31, 603–632. [Google Scholar] [CrossRef]
  16. Castro, R.; Guerrero-Legarreta, I.; Borquez, R. Chitin extraction from Allopetrolisthes punctatus crab using lactic fermentation. Biotechnol. Rep. 2018, 20, e00287. [Google Scholar] [CrossRef]
  17. Gong, X.; Tian, W.; Bai, J.; Qiao, K.; Zhao, J.; Wang, L. Highly efficient deproteinization with an ammonifying bacteria Lysinibacillus fusiformis isolated from brewery spent diatomite. J. Biosci. Bioeng. 2019, 127, 326–332. [Google Scholar] [CrossRef]
  18. Abirami, S.; Nagarajan, D. Extraction of Chitin from Shrimp Shell Wastes by Using Bacillus licheniformis and Lactobacillus plantarum. Inter. J. Recent Res. Asp. 2018, 307–315. [Google Scholar]
  19. Beier, S.; Bertilsson, S. Bacterial chitin degradation-mechanisms and ecophysiological strategies. Front. Microbiol 2013, 4, 149. [Google Scholar] [CrossRef] [Green Version]
  20. Iwasaki, Y.; Ichino, T.; Saito, A. Transition of the Bacterial Community and Culturable Chitinolytic Bacteria in Chitin-treated Upland Soil: From Streptomyces to Methionine-auxotrophic Lysobacter and Other Genera. Microbes Environ. 2020, 35, ME19070. [Google Scholar] [CrossRef] [Green Version]
  21. Jacquiod, S.; Franqueville, L.; Cecillon, S.; Vogel, T.M.; Simonet, P. Soil bacterial community shifts after chitin enrichment: An integrative metagenomic approach. PLoS ONE 2013, 8, e79699. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Wieczorek, A.S.; Schmidt, O.; Chatzinotas, A.; von Bergen, M.; Gorissen, A.; Kolb, S. Ecological Functions of Agricultural Soil Bacteria and Microeukaryotes in Chitin Degradation: A Case Study. Front. Microbiol. 2019, 10, 1293. [Google Scholar] [CrossRef] [PubMed]
  23. Baty, A.M., 3rd; Eastburn, C.C.; Diwu, Z.; Techkarnjanaruk, S.; Goodman, A.E.; Geesey, G.G. Differentiation of chitinase-active and non-chitinase-active subpopulations of a marine bacterium during chitin degradation. Appl. Environ. Microbiol. 2000, 66, 3566–3573. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Baty, A.M., 3rd; Eastburn, C.C.; Techkarnjanaruk, S.; Goodman, A.E.; Geesey, G.G. Spatial and temporal variations in chitinolytic gene expression and bacterial biomass production during chitin degradation. Appl. Environ. Microbiol. 2000, 66, 3574–3585. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Klassen, G.; Locke, A. A biological synopsis of the European green crab, Carcinus maenas. Can. Manuscr. Rep. Fish. Aquat. Sci. 2007, 2818, 1–82. [Google Scholar] [CrossRef]
  26. Skonberg, D.I.; Perkins, B.L. Nutrient composition of green crab (Carcinus maenus) leg meat and claw meat. Food Chem. 2002, 77, 401–404. [Google Scholar] [CrossRef]
  27. Hirano, S. Chitin biotechnology applications. Biotechnol. Ann. Rev. 1996, 2, 237–258. [Google Scholar] [CrossRef]
  28. Maksimov, I.V.; Abizgil’dina, R.R.; Pusenkova, L.I. Plant growth promoting rhizobacteria as alternative to chemical crop protectors from pathogens (review). Appl. Biochem. Microbiol. 2011, 47, 333–345. [Google Scholar] [CrossRef]
  29. Comeau, A.M.; Vincent, W.F.; Bernier, L.; Lovejoy, C. Novel chytrid lineages dominate fungal sequences in diverse marine and freshwater habitats. Sci. Rep. 2016, 6, 30120. [Google Scholar] [CrossRef]
  30. Yurgel, S.N.; Douglas, G.M.; Comeau, A.M.; Mammoliti, M.; Dusault, A.; Percival, D.; Langille, M.G. Variation in Bacterial and Eukaryotic Communities Associated with Natural and Managed Wild Blueberry Habitats. Phytobioms J. 2017, 1, 102–113. [Google Scholar] [CrossRef] [Green Version]
  31. Zhang, J.; Kobert, K.; Flouri, T.; Stamatakis, A. PEAR: A fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics 2014, 30, 614–620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef] [PubMed]
  33. Comeau, A.M.; Douglas, G.M.; Langille, M.G. Microbiome Helper: A Custom and Streamlined Workflow for Microbiome Research. mSystems 2017, 2, e00127-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  35. Amir, A.; McDonald, D.; Navas-Molina, J.A.; Kopylova, E.; Morton, J.T.; Zech Xu, Z.; Kightley, E.P.; Thompson, L.R.; Hyde, E.R.; Gonzalez, A.; et al. Deblur Rapidly Resolves Single-Nucleotide Community Sequence Patterns. mSystems 2017, 2, e00191-16. [Google Scholar] [CrossRef] [Green Version]
  36. Callahan, B.J.; McMurdie, P.J.; Holmes, S.P. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 2017, 11, 2639–2643. [Google Scholar] [CrossRef] [Green Version]
  37. Quast, C.; Pruesse, E.; Yilmaz, P.; Gerken, J.; Schweer, T.; Yarza, P.; Peplies, J.; Glockner, F.O. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013, 41, D590–D596. [Google Scholar] [CrossRef]
  38. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2′s q2-feature-classifier plugin. Microbiome 2018, 6, 90. [Google Scholar] [CrossRef]
  39. Lozupone, C.; Lladser, M.E.; Knights, D.; Stombaugh, J.; Knight, R. UniFrac: An effective distance metric for microbial community comparison. ISME J. 2011, 5, 169–172. [Google Scholar] [CrossRef] [Green Version]
  40. Kim, B.R.; Shin, J.; Guevarra, R.; Lee, J.H.; Kim, D.W.; Seol, K.H.; Lee, J.H.; Kim, H.B.; Isaacson, R. Deciphering Diversity Indices for a Better Understanding of Microbial Communities. J. Microbiol. Biotechnol. 2017, 27, 2089–2093. [Google Scholar] [CrossRef] [Green Version]
  41. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. vegan: Community Ecology Package; R Package Version 2.4-0; R Foundation for Statistical Computing: Vienna, Austria, 2016; Available online: https://CRAN.R-project.org/package=vegan (accessed on 5 March 2022).
  42. Fernandes, A.D.; Reid, J.N.; Macklaim, J.M.; McMurrough, T.A.; Edgell, D.R.; Gloor, G.B. Unifying the analysis of high-throughput sequencing datasets: Characterizing RNA-seq, 16S rRNA gene sequencing and selective growth experiments by compositional data analysis. Microbiome 2014, 2, 15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  44. Schwager, E.; Weingart, G.; Bielski, C.; Huttenhower, C. CCREPE: Compositionality Corrected by PErmutation and REnormalization. 2020. Available online: bioconductor.org (accessed on 20 April 2022).
  45. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef] [PubMed]
  46. De Caceres, M.; Legendre, P. Associations between species and groups of sites: Indices and statistical inference. Ecology 2009, 90, 3566–3574. [Google Scholar] [CrossRef] [PubMed]
  47. Swiontek Brzezinska, M.; Jankiewicz, U.; Burkowska, A.; Walczak, M. Chitinolytic microorganisms and their possible application in environmental protection. Curr. Microbiol. 2014, 68, 71–81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Xia, J.L.; Xiong, J.; Xu, T.; Zhang, C.G.; Zhang, R.Y.; Zhang, Q.; Wu, S.; Qiu, G.Z. Purification and characterization of extracellular chitinase from a novel strain Aspergillus fumigatus CS-01. J. Cent. South Univ. Technol. 2009, 16, 552–557. [Google Scholar] [CrossRef]
  49. Han, S.J.; Park, H.; Lee, S.G.; Lee, H.K.; Yim, J.H. Optimization of cold-active chitinase production from the Antarctic bacterium, Sanguibacter antarcticus KOPRI 21702. Appl. Microbiol. Biotechnol. 2011, 89, 613–621. [Google Scholar] [CrossRef]
  50. Starke, R.; Morais, D.; Vetrovsky, T.; Lopez Mondejar, R.; Baldrian, P.; Brabcova, V. Feeding on fungi: Genomic and proteomic analysis of the enzymatic machinery of bacteria decomposing fungal biomass. Environ. Microbiol. 2020, 22, 4604–4619. [Google Scholar] [CrossRef]
  51. Someya, N.; Ikeda, S.; Morohoshi, T.; Noguchi Tsujimoto, M.; Yoshida, T.; Sawada, H.; Ikeda, T.; Tsuchiya, K. Diversity of culturable chitinolytic bacteria from rhizospheres of agronomic plants in Japan. Microbes Environ. 2011, 26, 7–14. [Google Scholar] [CrossRef] [Green Version]
  52. Larsbrink, J.; Zhu, Y.; Kharade, S.S.; Kwiatkowski, K.J.; Eijsink, V.G.; Koropatkin, N.M.; McBride, M.J.; Pope, P.B. A polysaccharide utilization locus from Flavobacterium johnsoniae enables conversion of recalcitrant chitin. Biotechnol. Biofuels 2016, 9, 260. [Google Scholar] [CrossRef] [Green Version]
  53. Kharade, S.S.; McBride, M.J. Flavobacterium johnsoniae chitinase ChiA is required for chitin utilization and is secreted by the type IX secretion system. J. Bacteriol. 2014, 196, 961–970. [Google Scholar] [CrossRef] [Green Version]
  54. Wegner, C.E.; Liesack, W. Microbial community dynamics during the early stages of plant polymer breakdown in paddy soil. Environ. Microbiol. 2016, 18, 2825–2842. [Google Scholar] [CrossRef] [PubMed]
  55. Harkin, C.; Bruck, W.M.; Lynch, C. Isolation & identification of bacteria for the treatment of brown crab (Cancer pagurus) waste to produce chitinous material. J. Appl. Microbiol. 2015, 118, 954–965. [Google Scholar] [CrossRef] [PubMed]
  56. Ghorbel-Bellaaj, O.; Hmidet, N.; Jellouli, K.; Younes, I.; Maalej, H.; Hachicha, R.; Nasri, M. Shrimp waste fermentation with Pseudomonas aeruginosa A2: Optimization of chitin extraction conditions through Plackett-Burman and response surface methodology approaches. Int. J. Biol. Macromol. 2011, 48, 596–602. [Google Scholar] [CrossRef] [PubMed]
  57. Ren, F.; He, R.; Zhou, X.; Gu, Q.; Xia, Z.; Liang, M.; Zhou, J.; Lin, B.; Zou, C. Dynamic changes in fermentation profiles and bacterial community composition during sugarcane top silage fermentation: A preliminary study. Bioresour Technol. 2019, 285, 121315. [Google Scholar] [CrossRef]
  58. Ziganshina, E.E.; Bagmanova, A.R.; Khilyas, I.V.; Ziganshin, A.M. Assessment of a biogas-generating microbial community in a pilot-scale anaerobic reactor. J. Biosci. Bioeng. 2014, 117, 730–736. [Google Scholar] [CrossRef]
  59. Pang, L.; He, Y.; Liu, X.; Li, J.; Yang, P. The role of a newly isolated strain Corynebacterium pollutisoli SPH6 in waste activated sludge alkaline fermentation. Chemosphere 2020, 241, 125072. [Google Scholar] [CrossRef]
  60. Stiles, M.E.; Holzapfel, W.H. Lactic acid bacteria of foods and their current taxonomy. Int. J. Food Microbiol. 1997, 36, 1–29. [Google Scholar] [CrossRef]
  61. Johnson, J.M.; Ludwig, A.; Furch, A.C.U.; Mithofer, A.; Scholz, S.; Reichelt, M.; Oelmuller, R. The Beneficial Root-Colonizing Fungus Mortierella hyalina Promotes the Aerial Growth of Arabidopsis and Activates Calcium-Dependent Responses That Restrict Alternaria brassicae-Induced Disease Development in Roots. Mol. Plant-Microbe Interact. 2019, 32, 351–363. [Google Scholar] [CrossRef] [Green Version]
  62. Ozimek, E. Mortierella Species as the Plant Growth-Promoting Fungi Present in the Agricultural Soils. Agriculture 2021, 11, 7. [Google Scholar] [CrossRef]
  63. Gomez Exposito, R.; de Bruijn, I.; Postma, J.; Raaijmakers, J.M. Current Insights into the Role of Rhizosphere Bacteria in Disease Suppressive Soils. Front. Microbiol. 2017, 8, 2529. [Google Scholar] [CrossRef]
  64. Duo-Chuan, L.I.; Chen, S.; Jing, L.U. Purification and partial characterization of two chitinases from the mycoparasitic fungus Talaromyces flavus. Mycopathologia 2005, 159, 223–229. [Google Scholar] [CrossRef] [PubMed]
  65. Liu, Z.H.; Yang, Q.; Hu, S.; Zhang, J.D.; Ma, J. Cloning and characterization of a novel chitinase gene (chi46) from Chaetomium globosum and identification of its biological activity. Appl. Microbiol. Biotechnol. 2008, 80, 241–252. [Google Scholar] [CrossRef] [PubMed]
  66. Dwyer, K.; Bentley, I.S.; Tighe, E.; McGleenan, E.; Gaffney, D.; Walsh, G. Recombinant production and characterisation of two chitinases from Rasamsonia emersonii, and assessment of their potential industrial applicability. Appl. Microbiol. Biotechnol. 2021, 105, 7769–7783. [Google Scholar] [CrossRef] [PubMed]
  67. Zang, J.; Xu, Y.; Xia, W.; Yu, D.; Gao, P.; Jiang, Q.; Yang, F. Dynamics and diversity of microbial community succession during fermentation of Suan yu, a Chinese traditional fermented fish, determined by high throughput sequencing. Food Res. Int. 2018, 111, 565–573. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Nonmetric multidimensional scaling (NMDS) of bacterial and fungal communities at ASV level. 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.
Figure 1. Nonmetric multidimensional scaling (NMDS) of bacterial and fungal communities at ASV level. 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.
Microorganisms 10 01033 g001
Figure 2. Estimated Shannon diversity and Evenness of bacterial and fungal communities. For each variable, data followed by different letters are significantly different according to Kruskal-Wallis pairwise test (p < 0.05).
Figure 2. Estimated Shannon diversity and Evenness of bacterial and fungal communities. For each variable, data followed by different letters are significantly different according to Kruskal-Wallis pairwise test (p < 0.05).
Microorganisms 10 01033 g002
Figure 3. Bacterial and fungal taxa that were differentially represented between Shell, Soil, and Control microbiomes. Based on ALDEx2 Benjamini–Hochberg-corrected p-value of Kruskal–Wallace test.
Figure 3. Bacterial and fungal taxa that were differentially represented between Shell, Soil, and Control microbiomes. Based on ALDEx2 Benjamini–Hochberg-corrected p-value of Kruskal–Wallace test.
Microorganisms 10 01033 g003
Figure 4. Venn diagrams showing the overlap between microbial communities from Shell, Soil, and Control microbiomes. (A) Bacterial and (B) fungal ASVs found in total microbiomes; (C) bacterial and (D) fungal ASVs found in core microbiomes; (E) bacterial and fungal ASVs found in network interaction.
Figure 4. Venn diagrams showing the overlap between microbial communities from Shell, Soil, and Control microbiomes. (A) Bacterial and (B) fungal ASVs found in total microbiomes; (C) bacterial and (D) fungal ASVs found in core microbiomes; (E) bacterial and fungal ASVs found in network interaction.
Microorganisms 10 01033 g004
Figure 5. Microbial families identified in the core microbiomes. (A) Bacterial 16S rRNA; (B) fungal ITS.
Figure 5. Microbial families identified in the core microbiomes. (A) Bacterial 16S rRNA; (B) fungal ITS.
Microorganisms 10 01033 g005
Table 1. 16S rRNA and ITS ASVs with potential to be directly involved in shell degradation.
Table 1. 16S rRNA and ITS ASVs with potential to be directly involved in shell degradation.
Genera/Lowest AnnotationTotal *ControlSoilShellTaxonomy
16S ASVs
Sanguibacter3.0950.0000.2725.820Actinobacteriota Actinobacteria Micrococcales Sanguibacteraceae
Pseudomonas, 4 **3.1400.0061.7385.456Proteobacteria Gammaproteobacteria Pseudomonadales Pseudomonadaceae
Streptomyces, 22.2470.0000.4284.156Actinobacteriota Actinobacteria Streptomycetales Streptomycetaceae
Flavobacterium1.3800.0000.5572.462Bacteroidota Bacteroidia Flavobacteriales Flavobacteriaceae
Tissierella, 41.1650.0000.6412.026Firmicutes Clostridia Peptostreptococcales-Tissierellales
Dongia, 21.1820.0001.9801.651Proteobacteria Alphaproteobacteria Dongiales Dongiaceae
Vagococcus0.7920.0000.4991.358Firmicutes Bacilli Lactobacillales Vagococcaceae
Ruminococcaceae0.6280.0000.2121.132Firmicutes Clostridia Oscillospirales
Chryseolinea0.9350.0002.2191.108Bacteroidota Bacteroidia Cytophagales Microscillaceae
Micrococcaceae0.5160.0000.3310.884Actinobacteriota Actinobacteria Micrococcales
Clostridium, 20.4750.0000.6550.706Firmicutes Clostridia Clostridiales Clostridiaceae
Erysipelothrix0.4100.0000.2750.697Firmicutes Bacilli Erysipelotrichales Erysipelotrichaceae
Stenotrophomonas0.3460.0000.0680.640Proteobacteria Gammaproteobacteria Xanthomonadales Xanthomonadaceae
Microbacteriaceae, 20.3570.0020.1440.635Actinobacteriota Actinobacteria Micrococcales
Terrimonas0.4060.0000.6720.570Bacteroidota Bacteroidia Chitinophagales Chitinophagaceae
Chitinophagaceae0.3340.0000.2250.569Bacteroidota Bacteroidia Chitinophagales
Longispora0.3000.0000.1560.525Actinobacteriota Actinobacteria Micromonosporales Micromonosporaceae
Rhizobiales, 20.3360.0000.4640.499Proteobacteria Alphaproteobacteria
Ignatzschineria, 20.2790.0000.1440.488Proteobacteria Gammaproteobacteria Cardiobacteriales Wohlfahrtiimonadaceae
Leifsonia, 20.2910.0000.3070.461Actinobacteriota Actinobacteria Micrococcales Microbacteriaceae
Dokdonella0.20.0000.0890.429Proteobacteria Gammaproteobacteria Xanthomonadales Rhodanobacteraceae
Peptostreptococcus0.2450.0000.1410.424Firmicutes Clostridia Peptostreptococcales-Tissierellales Peptostreptococcaceae
Rhodanobacteraceae0.1740.0000.0560.316Proteobacteria Gammaproteobacteria Xanthomonadales
Pseudarthrobacter0.1800.0000.1380.301Actinobacteriota Actinobacteria Micrococcales Micrococcaceae
Lysinibacillus0.1510.0000.0030.286Firmicutes Bacilli Bacillales Planococcaceae
Blastocatellaceae0.1590.0000.1030.273Acidobacteriota Blastocatellia Blastocatellales
Lysobacter0.1560.0000.0820.272Proteobacteria Gammaproteobacteria Xanthomonadales Xanthomonadaceae
Bryobacter0.1640.0000.1540.267Acidobacteriota Acidobacteriae Bryobacterales Bryobacteraceae
Subgroup 100.1950.0000.3720.259Acidobacteriota Thermoanaerobaculia Thermoanaerobaculales Thermoanaerobaculaceae
Psychrobacter0.1480.0000.0770.258Proteobacteria Gammaproteobacteria Pseudomonadales Moraxellaceae
Shinella0.1400.0000.0500.251Proteobacteria Alphaproteobacteria Rhizobiales Rhizobiaceae
Flavihumibacter0.2300.0000.6220.249Bacteroidota Bacteroidia Chitinophagales Chitinophagaceae
Phenylobacterium, 20.1320.0000.0500.237Proteobacteria Alphaproteobacteria Caulobacterales Caulobacteraceae
Paeniglutamicibacter0.1290.0000.0540.230Actinobacteriota Actinobacteria Micrococcales Micrococcaceae
Psychrobacillus0.1130.0000.0120.212Firmicutes Bacilli Bacillales Planococcaceae
Solirubrobacteraceae0.1560.0000.3280.198Actinobacteriota Thermoleophilia Solirubrobacterales
Sandaracinus0.1060.0000.0350.192Myxococcota Polyangia Polyangiales Sandaracinaceae
Methylomonadaceae0.1300.0000.2190.182Proteobacteria Gammaproteobacteria Methylococcales
Lachnospiraceae0.1020.0000.0680.173Firmicutes Clostridia Lachnospirales
Allocatelliglobosispora0.0790.0000.0030.149Actinobacteriota Actinobacteria Micromonosporales Micromonosporaceae
Romboutsia0.0990.0000.1410.147Firmicutes Clostridia Peptostreptococcales-Tissierellales Peptostreptococcaceae
Bosea0.0760.0000.0480.130Proteobacteria Alphaproteobacteria Rhizobiales Beijerinckiaceae
Flavitalea0.0720.0000.0410.125Bacteroidota Bacteroidia Chitinophagales Chitinophagaceae
Sphingomonadaceae0.0870.0000.1450.122Proteobacteria Alphaproteobacteria Sphingomonadales Sphingomonadaceae
Microscillaceae0.1070.0000.3160.108Bacteroidota Bacteroidia Cytophagales
Pedobacter0.0570.0000.0610.090Bacteroidota Bacteroidia Sphingobacteriales Sphingobacteriaceae
ITS2 ASVs
Mortierella, 516.2130.8049.62630.505Mortierellomycota Mortierellomycetes Mortierellales Mortierellaceae
Talaromyces, 84.8360.0066.9955.054Ascomycota Eurotiomycetes Eurotiales Trichocomaceae
Mycochlamys, 24.7840.0037.9233.990Ascomycota Sordariomycetes Microascales Microascaceae
Pseudeurotium2.9980.0285.2802.169Ascomycota Leotiomycetes Thelebolales Pseudeurotiaceae
Fusicolla, 21.2040.0201.0821.914Ascomycota Sordariomycetes Hypocreales Nectriaceae
Sordariales1.1570.0061.1951.691Ascomycota Sordariomycetes
Pseudogymnoascus1.0170.0280.9441.583Ascomycota Leotiomycetes Thelebolales Pseudeurotiaceae
Zopfiella2.6510.0005.0461.549Ascomycota Sordariomycetes Sordariales Chaetomiaceae
Apiotrichum0.7240.0060.7081.096Basidiomycota Tremellomycetes Trichosporonales Trichosporonaceae
Scedosporium0.4560.0000.2370.905Ascomycota Sordariomycetes Microascales Microascaceae
Trichosporonaceae0.3650.0030.1560.756Basidiomycota Tremellomycetes Trichosporonales
Coniochaeta, 21.1180.0262.1290.641Ascomycota Sordariomycetes Coniochaetales Coniochaetaceae
Aspergillus, 31.2500.0012.5690.538Ascomycota Eurotiomycetes Eurotiales Aspergillaceae
Coprinopsis0.5420.0000.8270.525Basidiomycota Agaricomycetes Agaricales Psathyrellaceae
Humicola0.6040.0001.0490.456Ascomycota Sordariomycetes Sordariales Chaetomiaceae
Rasamsonia0.4430.0000.6950.408Ascomycota Eurotiomycetes Eurotiales Trichocomaceae
Paracremonium0.4070.0000.6160.399Ascomycota Sordariomycetes Hypocreales Nectriaceae
Parascedosporium0.3080.0000.4670.300Ascomycota Sordariomycetes Microascales Microascaceae
Cephalotrichum, 20.2120.0000.2780.252Ascomycota Sordariomycetes Microascales Microascaceae
Thermomyces0.4210.0000.8100.237Ascomycota Eurotiomycetes Eurotiales Trichocomaceae
Pseudeurotiaceae0.2540.0000.4070.226Ascomycota Leotiomycetes Thelebolales
Anguillospora0.3290.0210.5840.224Ascomycota Leotiomycetes Helotiales Helotiaceae
Arthrographis0.1810.0010.2490.201Ascomycota Dothideomycetes Eremomycetaceae
Cystobasidiomycetes0.1910.0000.2950.182Basidiomycota
Candida0.2040.0080.3430.162Ascomycota Saccharomycetes Saccharomycetales
Leotiomycetes0.2750.0000.5290.155Ascomycota
Sugiyamaella0.3440.0070.7070.144Ascomycota Saccharomycetes Saccharomycetales Trichomonascaceae
* Relative abundance in percent (%) of 16S and or ITS2 reads in the specific niche and ** Number of ASVs with same taxonomic annotation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yurgel, S.N.; Nadeem, M.; Cheema, M. Microbial Consortium Associated with Crustacean Shells Composting. Microorganisms 2022, 10, 1033. https://doi.org/10.3390/microorganisms10051033

AMA Style

Yurgel SN, Nadeem M, Cheema M. Microbial Consortium Associated with Crustacean Shells Composting. Microorganisms. 2022; 10(5):1033. https://doi.org/10.3390/microorganisms10051033

Chicago/Turabian Style

Yurgel, Svetlana N., Muhammad Nadeem, and Mumtaz Cheema. 2022. "Microbial Consortium Associated with Crustacean Shells Composting" Microorganisms 10, no. 5: 1033. https://doi.org/10.3390/microorganisms10051033

APA Style

Yurgel, S. N., Nadeem, M., & Cheema, M. (2022). Microbial Consortium Associated with Crustacean Shells Composting. Microorganisms, 10(5), 1033. https://doi.org/10.3390/microorganisms10051033

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop