Next Article in Journal
The Small Ras Superfamily GTPase Rho4 of the Maize Anthracnose Fungus Colletotrichum graminicola Is Required for β-1,3-glucan Synthesis, Cell Wall Integrity, and Full Virulence
Next Article in Special Issue
Identification of Key Regulatory Pathways of Basidiocarp Formation in Pleurotus spp. Using Modeling, Simulation and System Biology Studies
Previous Article in Journal
Morphological and Phylogenetic Characterizations Reveal Five New Species of Astrothelium (Trypetheliales, Ascomycota) from China
Previous Article in Special Issue
Research Progress on Elements of Wild Edible Mushrooms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metaprofiling of the Bacterial Community in Colonized Compost Extracts by Agaricus subrufescens

by
Matheus Rodrigo Iossi
1,
Isabela Arruda Palú
2,
Douglas Moraes Soares
3,
Wagner G. Vieira, Jr.
1,
Lucas Silva Alves
1,
Cassius V. Stevani
3,
Cinthia E. C. Caitano
1,
Samir V. F. Atum
3,4,
Renato S. Freire
3,
Eustáquio S. Dias
5 and
Diego Cunha Zied
2,*
1
Programa de Pós-Graduação em Microbiologia Agropecuária, Faculdade de Ciências Agrárias e Veterinárias (FCAV), Universidade Estadual Paulista (UNESP), São Paulo 14884-900, Brazil
2
Faculdade de Ciências Agrárias e Tecnológicas (FCAT), Universidade Estadual Paulista (UNESP), São Paulo 17900-000, Brazil
3
Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo (USP), São Paulo 05508-220, Brazil
4
Departamento de Bioquímica, Instituto de Química, Universidade de São Paulo, São Paulo 05508-220, Brazil
5
Departamento de Biologia, Universidade Federal de Lavras (UFLA), Lavras 37200-900, Brazil
*
Author to whom correspondence should be addressed.
J. Fungi 2022, 8(10), 995; https://doi.org/10.3390/jof8100995
Submission received: 23 August 2022 / Revised: 15 September 2022 / Accepted: 19 September 2022 / Published: 22 September 2022
(This article belongs to the Special Issue Edible and Medicinal Macrofungi)

Abstract

:
It is well-known that bacteria and fungi play important roles in the relationships between mycelium growth and the formation of fruiting bodies. The sun mushroom, Agaricus subrufescens, was discovered in Brazil ca. 1960 and it has become known worldwide due to its medicinal and nutritional properties. This work evaluated the bacterial community present in mushroom-colonized compost extract (MCCE) prepared from cultivation of A. subrufescens, its dynamics with two different soaking times and the influence of the application of those extracts on the casing layer of a new compost block for A. subrufescens cultivation. MCCEs were prepared through initial submersion of the colonized compost for 1 h or 24 h in water followed by application on casing under semi-controlled conditions. Full-length 16S rRNA genes of 1 h and 24 h soaked MCCE were amplified and sequenced using nanopore technology. Proteobacteria, followed by Firmicutes and Planctomycetes, were found to be the most abundant phyla in both the 1 h and 24 h soaked MCCE. A total of 275 different bacterial species were classified from 1 h soaked MCCE samples and 166 species from 24 h soaked MCCE, indicating a decrease in the bacterial diversity with longer soaking time during the preparation of MCCE. The application of 24 h soaked MCCE provided increases of 25% in biological efficiency, 16% in precociousness, 53% in the number of mushrooms and 40% in mushroom weight compared to control. Further investigation is required to determine strategies to enhance the yield and quality of the agronomic traits in commercial mushroom cultivation.

1. Introduction

Agaricus subrufescens is a Brazilian Agaricales fungus, the mushrooms of which are considered functional foods with medicinal and nutraceutical properties, resulting in its popularization across the world in recent decades [1]. As mushroom cultivation does not require special conditions, this species is cultivated in a rustic way without controlled conditions (e.g., temperature, relative humidity, light period) by small growers, contributing to increases in family incomes [2,3,4]. In fact, this type of cultivation can, in some cases, favor productivity when linked to the adaptability of the strain.
Mushrooms of the genus Agaricus require a casing, a peat layer over the compost colonized by the mycelium, to fructify [5]. Aside from providing physical support for the formation of mushrooms [6], this creates the physicochemical and biological conditions that promote fruiting [7]. The material used must meet some requirements, including a medium that is close to neutral (pH 7), high water-retention and -release capacities, porosity that allows the exchange of gases between the substrate and the medium and resistance to stabilization after frequent irrigation [6,7,8,9,10].
Despite the existence of some reports on the casing irrigation in Agaricus bisporus cultivation [9,11], there are few specific studies for A. subrufescens, which opens the possibility of exploring this technique with the aim of improving production performance. Furthermore, the concomitant use of chemical and/or biological agents with the irrigation can provide an opportunity to modify the growth medium in the casing layer and contribute to increasing the yield and quality of mushrooms, as successfully demonstrated by the application of fungicides at low concentrations [12] and by a study that used a spent mushroom substrate extract as a biological control for mushroom diseases [13].
In addition to the physicochemical parameters, bacteria–fungi interactions also have a strong influence on mushroom production, ranging from antagonism to mutualism. Recent reports have attributed an important role to the compost and casing microbiome in the growth and fructification of cultivated mushrooms, mainly by influencing the quality and maturation of the compost, mycelium running, primordia induction and formation of fruit bodies [5,14,15].
In addition to providing water supply and maintaining a humid microclimate for the casing, the use of mushroom-colonized compost extract (MCCE) may modify the microbiome of the casing. Thus, the characterization of microbial communities occurring in MCCEs and the effects of their application on the cultivation of commercial mushrooms may reveal important bacteria–fungi associations related to higher yields and superior quality for agronomic traits [5]. Current long-read DNA sequencing technology, including portable nanopore devices, have been successfully applied for the analysis of full-length 16S rRNA gene amplicons, providing a reliable, rapid and cost-effective strategy to assess bacterial communities [16,17,18].
In this scenario, we evaluated the bacterial diversity of MCCE obtained from the cultivation of A. subrufescens and the effects of its application on the casing layer of a new compost block for the mushroom’s cultivation. This is the first report on 16S rRNA gene-based nanopore metagenomics applied to the microbiome analysis of substrates used for A. subrufescens production. Together with data from agronomic traits of mushrooms, these new findings may help to bring new insights about the influences of the bacterial communities, occurring in both the compost and casing on A. subrufescens mycelium growth.

2. Materials and Methods

2.1. Experimental Design

Experiments were conducted under a completely randomized design and carried out in a double factorial scheme (5 application dates × 3 extraction times, including control) including 15 treatments with 4 replications, resulting in 60 experimental units.

2.2. Compost Preparation

The substrate was prepared using the traditional method of composting for the cultivation of A. bisporus. Briefly, phase 1 lasted 22 days and phase 2 lasted 7 days, totaling 29 days. For the substrate formulation, 900 kg of wheat straw, 1400 kg of sugarcane bagasse, 45 kg of soybean meal, 4.5 kg of urea, 5 kg of ammonium sulfate, 10 kg of simple superphosphate and 40 kg of limestone were used. Bulky materials (wheat straw and sugarcane bagasse) were moistened for 9 days and turned over every 2 days. During phase 1, the concentrated materials (soybean meal, urea, ammonium sulfate, simple superphosphate and limestone) were added over the substrate windrow after each turn [19]. In phase 2, the substrate was maintained for 18 h at 59 ± 1 °C for pasteurization, followed by conditioning for 5 days at 47 ± 2 °C. All materials were purchased from Compobras® (Castro, PR, Brazil), and the mineral analysis is shown in Table S1.

2.3. Spawn Production

The strain used in the experiments was A. subrufescens ABL 04/49, which is deposited in the public mycological collection of the FCAT at the Centro de Estudos em Cogumelos (Dracena, SP, Brazil). The inoculum was produced from sorghum grains (Sorghum bicolor) cooked for approximately 40 min in boiling water. After cooking, excess water was drained and 2% calcite was added to the grains based on their wet weight to provide calcium and raise the pH [20]. After homogenization with calcite, the grains were autoclaved for a period of 4 h at 121 °C. After cooling, the grains were inoculated under aseptic conditions in a laminar flow chamber. The inoculum was incubated at a temperature of 28 °C. Substrate inoculation was performed using 1% inoculum for each experimental unit, which contained 2.8 kg wet weight. The experiments were conducted in a semi-controlled chamber with only a cooling system, which resulted in thermal oscillations between day and night periods.

2.4. Casing Layer and Fruiting Induction

Peat moss (+terra®, Castro, PR, Brazil) was used as casing layer, deposited approx. 3.0 cm above the leveled colonized compost [3]. The ruffling was performed before the mycelium reached the surface of the casing layer, 6 days after the addition of the casing [21]. After 24 h of ruffling, when the mycelium resumed its metabolic activity, the fruiting induction process began, and the temperature was decreased by 2 °C per day until it reached 20 °C, then increased by 2 °C per day until it reached 26 °C [21]. In the harvest phase, the CO2 content ranged from 600 to 1000 ppm.

2.5. Mushroom-Colonized Compost Extract (MCCE) Preparation and Addition to the Casing Layer

For the preparation of the extract, the compost colonized with the A. subrufescens ABL 04/49 strain 16 days after spawning was submerged in tap water in a 1:5 (substrate: water) ratio (Figure 1). Two MCCEs were prepared with different soaking times: 1 h and 24 h. Subsequently, the substrate was removed and the solution was filtered through a 10 mesh (2.0 mm) sieve. MCCEs were produced from the dilution of the filtered solution to a concentration of 12.5%; i.e., 87.5 L of water was added to 12.5 L of compost extract. Each experimental unit—boxes 15 cm × 30 cm × 40 cm (h × w × l) containing 2.8 kg of compost—was watered with 100 mL of either 1 h or 24 h soaked MCCE on the 3rd, 6th, 9th, 12th and 15th day after the addition of the casing layer to the compost. This volume was equivalent to 2.212 L of irrigation per square meter. The control treatment consisted of the water irrigation over the casing layer without MCCE.

2.6. Harvest Phase

The fruiting period started 17 days after the addition of the casing layer. Harvesting was performed manually up to three times a day and mushrooms were harvested with the closed cap [22]. Post-harvest procedures included the cutting of the base of the stipe and the washing of mushrooms to clean the cap and remove casing layer residues.

2.7. Agronomic Trait Evaluation

Data from agronomic traits collected daily were used to estimate the biological efficiency, precociousness, number and weight of the mushrooms. All parameters were calculated following Atila [23]. Biological efficiency was expressed as the ratio of the fresh mass of mushrooms and the dry mass of the initial substrate. Precociousness was calculated as the production of the first half of the harvest time divided by the total production. The number of mushrooms was calculated as the sum of all collected basidiomes. The weight of mushrooms was obtained from the ratio of the total mass of fresh mushrooms to the number of basidiomes collected.

2.8. Statiscal Analysis

Agronomic traits were measured in four replicates for the control and the 1 h and 24 h soaked MCCEs for each application date (Table S4). Standard deviations were calculated from results obtained under different conditions, with the presentation data of the media, accompanied by the range of variation verified in each repetition within the treatment using SAS statistical software (SAS Institute Inc., Cary, NC, USA) (Table S5).

2.9. Estimation of Bacterial Diversity by 16S rRNA Gene Sequencing

In order to characterize the 16S rRNA gene-based microbiome, 1 h and 24 h soaked MCCEs were filtered and the filter surface was used for DNA isolation with the DNeasy Power Water kit (QIAGEN, Germany), according to the manufacturer’s instructions. PCR amplification of full-length 16S rRNA genes was carried out in 50 µL reactions using Taq DNA Polymerase (Cellco Biotech, Sao Carlos, Brazil) and the universal primers 27F and 1492R [24]. PCR reactions were incubated in the SimpliAmp thermal cycler (Applied Biosystems, Waltham, MA, USA) under the following settings: initial denaturation at 95 °C for 2 min, 35 cycles of denaturation at 95 °C for 30 s, annealing at 52 °C for 30 s, extension at 72 °C for 90 s and a final extension at 72 °C for 120 s. DNA amplicons (~1.5 kb) were separated by agarose gel electrophoresis and purified using the QIAquick PCR purification kit (QIAGEN, Hilden, Germany). DNA libraries were prepared from 2 µg of the purified DNA amplicons. Double-strand DNA molecules were end-repaired using the NEBNext FFPE DNA Repair Mix and NEBNext Ultra II End Repair/dA-tailing Module reagents (New England Biolabs, Ipswich, MA, USA). A total of 10 unique DNA barcodes from the Native Barcoding Expansion kit (EXP-NBD104) (Oxford Nanopore Technologies, Oxford, UK) were attached to the end-prepped DNA (1–5 for 1 h and 6–10 for 24 h soaked MCCE samples). Equimolar amounts (350 ng) of each barcoded sample were combined to produce a pooled sample, used for the adapter ligation with the Ligation Sequencing kit (SQK-LSK109) (Oxford Nanopore Technologies, UK). The DNA library was cleaned up using the Agencourt AMPure XP beads (Beckman Coulter, Brea, CA, USA) and quantified using the Qubit fluorometer (Applied Biosystems, USA). A sequencing library was loaded into the minION Mk 1B flow cell R9.4.1 and the nanopore sequencing was performed over 48 h.

2.10. Bioinformatic Analysis

Basecalling, demultiplexing and quality control were undertaken in guppy v. 5.0.16 with the high quality r9.4.1_450bps_hac model. Screening against possible contaminants was conducted using BLAST v. 2.12 [25] against the UniVec database v. 10.0 (http://www.ncbi.nlm.nih.gov/tools/vecscreen/univec/, accessed on 14 August 2022). Nanofilt v. 2.8.0 [26] was used to filter the reads by size, selecting only those between 1400 and 1600 bp. After filtering, kraken2 v.2.1.1 [27] was used with the Greengenes database (v. 13.5 downloaded at 24 March 2022) for classification. To support the Greengenes results, minimap2 v. 2.22 [28] was also used to align the reads to the NCBI 16S RefSeq database (https://www.ncbi.nlm.nih.gov/refseq/targetedloci/16S_process/—downloaded on 24 March 2022). For visualization of the results, Pavian v. 1.0 [29] was used to generate Sankey graphs and classification tables of phyla and genera ranked by the number of reads. Data from barcode A, corresponding to the 16S-based microbiome for 1 h soaked MCCE, were combined, as were those from barcode B, representing the 24 h soaked MCCE samples (Tables S2 and S3).

3. Results

3.1. Composition and Dynamics of MCCE Microbiome

To investigate the bacterial diversity and the microbiome dynamics in 1 h and 24 h soaked MCCEs, these solutions were initially filtered, and the membrane was used for DNA isolation. The full-length 16S rRNA gene was successfully amplified, and the amplicons were properly barcoded and prepared for nanopore sequencing. After bioinformatics processing, including basecalling, demultiplexing and quality control steps, 1,269,157 reads were obtained and 99.85% of them were taxonomically classified.
Proteobacteria, followed by Firmicutes and Planctomycetes, were found to be the most abundant phyla in both 1 h and 24 h soaked MCCE (Figure 2 and Figure 3A). Proteobacteria was the dominant phylum in both MCCEs and its abundance was noticeably increased from 53.8% in 1 h soaked MCCE to 97.1% in 24 h soaked MCCE, which suggests that these changes in microbiome dynamics were influenced by the duration of the soaking step during MCCE preparation. Most Proteobacteria were assigned to the Gammaproteobacter class, which includes several Gram-negative microbes with scientific and health relevance, such as the families Enterobacteriaceae, Vibrionaceae and Pseudomonadaceae.
More than 92% of Firmicutes, the second most abundant phylum, were classified as Bacilli, a Gram-positive, aerobic and spore-forming bacteria class. The abundance of Firmicutes was dramatically reduced from 23% in 1 h soaked MCCE to 1.7% in 24 h soaked MCCE, which may have been related to lower oxygen availability during the longer soaking times, causing a reduction in the growth rate of this class of microorganisms.
Serratia was the most abundant genus in both MCCEs, increasing from 19.3% in 1 h soaked MCCE to 39.7% in 24 h soaked MCCE (Figure 3B). Likewise, the abundance of Klebsiella genus was slightly increased in 24 h soaked MCCE (from 7.2% to 10.3%). In addition to these two common genera, Bacillus, Paenibacillus and Planctomyces were the most abundant in 1 h soaked MCCE and Pseudomonas, Citrobacter and Erwinia in 24 h soaked MCCE (Figure 3B).
A total of 275 different bacterial species were classified from the 1 h soaked MCCE samples (Table S2) and 166 species from the 24 h soaked MCCE (Table S3), indicating a reduction in the bacterial diversity with longer soaking times during the preparation of MCCE. In both groups of samples, Serratia marcescens was the dominant species, with abundance significantly higher in 24 h soaked MCCE when compared to that obtained after 1 h of soaking (Tables S2 and S3). In a more discreet way, the abundance of Klebsiella oxytoca was positively correlated to S. marcescens, being slightly higher in 24 h soaked MCCE (Tables S2 and S3). A longer soaking time also contributed to increasing the abundance of bacteria from the Pseudomonas genus, which accounted for up to 15.9% of the total classified genera in 24 h soaked MCCE samples (Figure 3B).
A strong reduction in the percentage of genera classified as “others” was observed in the microbiome of the 24 h soaked MCCE samples (25.9%) in comparison to the 1 h soaked MCCE (56.7%), as well as a significant increment in the representativeness of the top five most abundant genera in the 24 h soaked MCCE microbiome (corresponding to 74.1% of the total genera classified) compared to the 1 h soaked MCCE microbiome (representing only 43.3% of the total genera), indicating a scenario of dramatic changes in the dynamics of the microbial diversity, and its reduction in MCCE prepared with a longer soaking time (Figure 3B).

3.2. Evaluation of the Application of MCCE on the Agronomic Traits of Mushrooms

Longer soaking time also contributed to improving all agronomic parameters. Figure 4 presents the average of the application times (3rd, 6th, 9th, 12th and 15th days after casing) with the standard deviation for each analyzed variable. 24 h soaking of the MCCE increased the biological efficiency by 25%, the precociousness by 16%, the number of mushrooms by 53% and the weight by 40% compared to control (Tables S4 and S5).
Regarding the 1 h soaked MCCE, it stood out for its shorter preparation time and lower precociousness. Despite the fact that the biological efficiency did not increase significantly, the precociousness decreased, which allowed a shorter crop cycle and a substantial increase of 46% in the weight of mushrooms. Evaluating the best duration for the application of MCCE in the casing layer, and considering the significant increase in biological efficiency, it is recommended to apply the 1 h MCCE between the 9th and 12th days and the 24 h MCCE on the 12th day after casing (Figure 5). The red band indicates the value of the control treatment, with biological efficiency of 20 ± 5%. The results obtained for each agronomic parameter are described in Table S4.

4. Discussion

Mushroom production technology has advanced with the optimization of agronomic parameters related to chemical, physical-chemical, microbiological and environmental aspects [30]. Several studies have demonstrated the role of microorganisms in mushroom production [31,32,33], and analysis of the microbiome in basidiomes revealed that it is highly influenced by the bacterial diversity of the casing [34]. It is evident that, from the inoculation phase to the fructification and harvesting of the mushroom, there is a microbial succession in the compost and casing. Actually, there are several reports concerning the changes in the microbiome dynamics in the casing and compost during mushroom cultivation [34,35,36].
Compared to the classical approaches used to assess microbial diversity, next-generation sequencing (NGS) offers a more powerful and holistic tool for the taxonomic classification of microorganisms, and it is able to detect them even in complex environmental samples and does not require time-consuming steps for enrichment, isolation and cultivation. Most of metagenomics studies characterizing microbiomes from substrates for mushroom cultivation rely on DNA sequencing of the V3-V4 portion of the 16S rRNA gene, a hypervariable region commonly used for bacterial taxonomy, using the short-read technology of Illumina platforms [34,36]. In contrast, we present here the first successful report of a full length 16S rRNA gene-based metagenomics analysis of A. subrufescens mushroom-colonized compost extracts (MCCEs) using the long-read sequencing technology of a portable nanopore sequencer (minION, Oxford Nanopore Technologies, UK), providing more complete genetic data for a proper taxonomic classification of reads.
Our sequencing data revealed the Proteobacteria to be the dominant phylum, with a significant increment in 24 h soaked MCCE. This finding is in accordance with previous microbiome analysis in the casing throughout the crop cycle of A. bisporus, a mushroom-forming species of the same genus of A. subrufescens [34,37]. Other abundant phyla, including Actinobacteria, Bacteroides, Firmicutes, Proteobacteria and Planctomycetes, were also found during the A. bisporus mushroom cultivation [15,34,38], and may be related to the beneficial relationship between these bacterial phyla and Agaricus fungi.
Similarly, Firmicutes and Actinobacteria were also found to be dominant phyla in the microbiomes of different developmental stages during Hypsizygus marmoreus mushroom cultivation [36]. Firmicutes bacteria mostly belong to the class Bacilli, represented by the most abundant genera Bacillus and Paenibacillus. Several studies have reported the use of Bacillus strains, naturally occurring in casings, as biological agents to suppress the growth of fungal pathogens and increase the yield of commercial mushrooms [5,39,40]. Despite being the third most abundant genus in the microbiome of the 1 h soaked MCCE, Bacillus abundance was strongly reduced with longer soaking time, being replaced with Pseudomonas in 24 h soaked MCCE.
Likewise, the relative abundance of Pseudomonas in the casing increased during the crop cycle of A. bisporus, becoming the second most abundant genus by the end of the second flush [34]. The beneficial effects of the proliferation of this genus in the casing are well-documented and mainly involve assisting the hyphae in decomposing complex substances present in the substrate and facilitating the efficient absorption of nutrients [36]. In addition, some species, such as Pseudomonas putida can break down the 1-octen-3-ol bond and induce fruiting body formation in A. bisporus [41,42]. Similarly, Pseudomonas fluorescens is able to remove A. bisporus deposits of crystalline calcium oxalate, a prerequisite for fruiting initiation [43], as well as promote the formation of the primordium and enhance the development of fruiting bodies of Agaricus bitorquis [42], Pleurotus eryngii [44] and Pleurotus ostreatus [45].
In both MCCE microbiomes, the Gram-negative and facultative anaerobic bacteria K. oxytoca and S. marcescens were found, S. marcescens being the most abundant bacterial species. It is known that S. marcescens exhibits high efficiency in the degradation of chitin through the activity of the chitinase enzyme [46]. Fungi have chitin in the cell wall [47], which contributes to the development of S. marcescens in this environment. The dominance of the Serratia genus was also reported during the mature period of H. marmoreus cultivation. In fact, there are several reports on the benefits of Serratia on mushroom cultivation, including its role in biocontrol [48,49], phosphorus solubilization [50] and growth promotion [36].
A recent report showed that the Serratia odorifera HZSO-1 strain, a symbiotic bacterium that resides in the H. marmoreus hyphae, can shorten the fruiting cycle by 3–4 days and increase the fruiting body yield by 12% [36]. Among secondary metabolites produced by S. odorifera HZSO-1, the secretion of quorum-sensing molecules, such as N-acyl homoserine lactones, could stimulate H. marmoreus to secrete some lignin-degrading enzymes, promoting growth and development and increasing the yield of mushrooms [36,51]. Furthermore, it has been demonstrated that the Serratia plymuthica PRI-2C strain produces the metabollite sodorifen in response to volatile organic compounds (VOCs) emitted by the fungal pathogen Fusarium culmorum, suggesting that terpenes may be important molecules in the communication between fungi and bacteria [52]. In this context, it seems plausible that the application of S. marcescens-containing MCCE could be explored towards obtaining better agronomic traits in mushroom-forming fungi.
The optimal date for application of MCCE to the casing layer for maximal biological efficiency was found to be the 12th day, when the formation of clumping occurred before the pin phases. During the first week after casing, there was a large increase in the bacterial population, which tended to stabilize or decrease after the 10th day. Precautions should be taken regarding the use of MCCE, such as using the same strain/species of the fungus as the one to be cultivated, checking for the absence of pests and diseases in the colonized compost and, finally, using newly colonized compost with high mycelial growth, which denotes an active metabolic state. The precautions should be the same as for the “compost add casing—CAC’ing” technique, described in detail by Pardo-Giménez et al. [53]. Thus, the use of mushroom growth-promoting microorganisms can involve two technological application levels: one with high investment, involving specific culture medium (solid or liquid) and known and improved species for such application; and another with low investment, involving a simple extraction process (known as “on farn”), with the selection of natural species that represent a high diversity of microorganisms.
Overall, the recent report on H. marmoreus suggests that the cultivation of edible mushrooms by adding beneficial microbes to the substrate is a feasible practice [36]. Several microorganisms present in these environments may strongly influence, and in some cases, are required for, the growth and fructification of cultivated mushrooms [5]. Understanding the dynamics of the microbiome in the compost and casing may determine the composition of the mushroom holobiont, the association of fungus and microorganisms [5]. Therefore, the elucidation of the microbiota structure and the symbiotic bacterial communities in mushroom substrates is a key step that can be explored for isolation of growth-promoting strains related to higher yield and better agronomic traits in mushrooms and for the development of synthetic communities for application to mushroom production and disease control [34,36].

5. Conclusions

Metagenomics of extracts prepared from A. subrufescens compost with different soaking times were used as a powerful tool to elucidate the dynamics of bacterial diversity and bacterial interactions with the hyphae of basidiomycetes during its cultivation. In this context, the long-read technology of nanopore sequencing is a great choice for full-length 16S rRNA gene-based metagenomics, providing the genetic data to explore the taxonomic classification in microbiome analysis. The 24 h soaked MCCE increased the biological efficiency, precociousness, number and weight of the mushrooms. Further investigation is required to determine strategies to enhance the yield and quality of agronomic traits in commercial mushroom cultivation.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jof8100995/s1. Table S1: Mineral analysis of the compost used in the experiment. Table S2: Abundance of bacteria present in the mushroom colonized compost extract (MCCE) obtained by immersion of the compost for 1 h. Table S3: Abundance of bacteria present in the mushroom colonized compost extract (MCCE) obtained by immersion of the compost for 24 h. Table S4: Agronomic traits of A. subrufescens mushrooms. Table S5: Mean values and standard deviation of agronomic traits of A. subrufescens mushrooms from the 3rd application until the 15th (see Table S2). Figure S1. Chamber and compost temperature, mushroom harvested and mushroom flushes, and oscillation temperature during crop period.

Author Contributions

Conceptualization, D.C.Z. and E.S.D.; methodology, M.R.I., D.M.S. and C.V.S.; investigation, M.R.I., I.A.P., W.G.V.J. and C.E.C.C.; data curation, L.S.A. and R.S.F.; writing—original draft preparation, M.R.I., I.A.P. and S.V.F.A.; writing—review and editing, D.M.S., C.V.S. and D.C.Z.; project administration, D.C.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundação de Amparo à Pesquisa do Estado de São Paulo (process numbers 2018/21492–2 and 2017/22501-2) and CAPES (grant no. 001) to Matheus Rodrigo Iossi.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Largeteau, M.L.; Llarena-Hernández, R.C.; Regnault-Roger, C.; Savoie, J.M. The medicinal Agaricus mushroom cultivated in Brazil: Biology, cultivation and non-medicinal valorisation. Appl. Microbiol. Biotechnol. 2011, 92, 897–907. [Google Scholar] [CrossRef] [PubMed]
  2. Tomizawa, M.M.; Dias, E.S.; Assis, L.J.D.; Gomide, P.H.O.; Santos, J.B.D. Genetic variability of mushroom isolates Agaricus blazei using markers RAPD. Ciência E Agrotecnologia 2007, 31, 1242–1249. [Google Scholar] [CrossRef]
  3. Wisitrassameewong, K.; Karunarathna, S.C.; Thongklang, N.; Zhao, R.; Callac, P.; Moukha, S.; Férandon, C.; Chukeatirote, E.; Hyde, K.D. Agaricus subrufescens: A review. Saudi J. Biol. Sci. 2012, 19, 131–146. [Google Scholar] [CrossRef] [PubMed]
  4. Colauto, N.B.; Linde, G.A. Avances sobre el cultivo del “cogumelo do sol” en Brasil. In El Cultivo de los Hongos Comestibles en Iberoamérica; Sánchez, J.E., Mata, G., Eds.; Colegio de la Frontera Sur: Tapachula, Mexico, 2012; pp. 103–128. [Google Scholar]
  5. Carrasco, J.; Preston, G.M. Growing edible mushrooms: A conversation between bacteria and fungi. Environ. Microbiol. 2020, 22, 858–872. [Google Scholar] [CrossRef]
  6. Ghasemi, K.; Emadi, M.; Bagheri, A.; Mohammadi, M. Casing material and thickness effects on the yield and nutrient concentration of Agaricus bisporus. Sarhad. J. Agric. 2020, 36, 958–965. [Google Scholar] [CrossRef]
  7. Murmu, R.; Maurya, A.K.; John, V. Mycoflora of certain casing materials used in the production of white button mushroom (Agaricus bisporus (Lange) Imbach). Int. J. Chem. Stud. 2020, 8, 2863–2868. [Google Scholar] [CrossRef]
  8. Levanon, D.; Danai, O.; Raz, D. Method of Indoor Mushroom Cultivation. European Patent Application EP 3106025A1, 2016. [Google Scholar]
  9. Navarro, M.J.; Gea, F.J.; Pardo-Giménez, A.; Martínez, A.; Raz, D.; Levanon, D.; Danay, O. Agronomical valuation of a drip irrigation system in a commercial mushroom farm. Sci. Hort. 2020, 265, 109234. [Google Scholar] [CrossRef]
  10. Kumar, P.; Kumar, V.; Adelodun, B.; Bedeković, D.; Kos, I.; Širić, I.; Silva, L.F. Sustainable Use of Sewage Sludge as a Casing Material for Button Mushroom (Agaricus bisporus) Cultivation: Experimental and Prediction Modeling Studies for Uptake of Metal Elements. J. Fungi 2022, 8, 112. [Google Scholar] [CrossRef] [PubMed]
  11. Noble, R.; Goodchild, M.; Jenkins, D.; Burton, K.; Walker, D.; Moorhouse, E. Precision Irrigation of Mushrooms Using a Novel Dielectric Tensiometer. Mushroom News 2017, 65, 8–17. [Google Scholar]
  12. Zied, D.C.; Caitano, C.E.C.; Pardo-Giménez, A.; Dias, E.; Zeraik, M.L.; Pardo, J.E. Using of appropriated strains in the practice of compost supplementation for Agaricus subrufescens production. Front. Sustain. Food Syst. 2018, 2, 26. [Google Scholar] [CrossRef]
  13. Gea, F.J.; Carrasco, J.; Diánez, F.; Santos, M.; Navarro, M.J. Control of dry bubble disease (Lecanicillium fungicola) in button mushroom (Agaricus bisporus) by spent mushroom substrate tea. Eur. J. Plant Pathol. 2014, 138, 711–720. [Google Scholar] [CrossRef]
  14. Suwannarach, N.; Kumla, J.; Zhao, Y.; Kakumyan, P. Impact of Cultivation Substrate and Microbial Community on Improving Mushroom Productivity: A Review. Biology 2022, 11, 569. [Google Scholar] [CrossRef] [PubMed]
  15. Vieira, F.R.; Pecchia, J.A. Bacterial community patterns in the Agaricus bisporus cultivation system, from compost raw materials to mushroom caps. Microb. Ecol. 2021, 84, 20–32. [Google Scholar] [CrossRef] [PubMed]
  16. Fujiyoshi, S.; Muto-Fujita, A.; Maruyama, F. Evaluation of PCR conditions for characterizing bacterial communities with full-length 16S rRNA genes using a portable nanopore sequencer. Sci. Rep. 2020, 10, 12580. [Google Scholar] [CrossRef] [PubMed]
  17. Ciuffreda, L.; Rodríguez-Pérez, H.; Flores, C. Nanopore sequencing and its application to the study of microbial communities. Comput. Struct. Biotechnol. J. 2021, 19, 1497–1511. [Google Scholar] [CrossRef]
  18. Matsuo, Y.; Komiya, S.; Yasumizu, Y.; Yasuoka, Y.; Mizushima, K.; Takagi, T.; Kryukov, K.; Fukuda, A.; Morimoto, Y.; Naito, Y.; et al. Full-length 16S rRNA gene amplicon analysis of human gut microbiota using MinION™ nanopore sequencing confers species-level resolution. BMC Microbiol. 2021, 21, 35. [Google Scholar] [CrossRef]
  19. Andrade, M.C.N.; Kopytowski Filho, J.; Minhoni, M.T.D.; Coutinho, L.N.; Figueiredo, M.B. Productivity, biological efficiency, and number of Agaricus blazei mushrooms grown in compost in the presence of Trichoderma sp. and Chaetomium olivacearum contaminants. Braz. J. Microbiol. 2007, 38, 243–247. [Google Scholar] [CrossRef]
  20. Aguiar, L.V.B.D.; Sales-Campos, C.; Gouvêa, P.R.D.S.; Vianez, B.F.; Dias, E.S.; Chevreuil, L.R. Substrate disinfection methods on the production and nutritional composition of a wild oyster mushroom from the Amazon. Ciência E Agrotecnologia 2021, 45, e010321. [Google Scholar] [CrossRef]
  21. Pardo-Giménez, A.; Pardo, J.E.; Dias, E.S.; Rinker, D.L.; Caitano, C.E.C.; Zied, D.C. Optimization of cultivation techniques improves the agronomic behavior of Agaricus subrufescens. Sci. Rep. 2020, 10, 1–9. [Google Scholar] [CrossRef]
  22. Stoknes, K.; Scholwin, F.; Jasinska, A.; Wojciechowska, E.; Mleczek, M.; Hanc, A.; Niedzielski, P. Cadmium mobility in a circular food-to-waste-to-food system and the use of a cultivated mushroom (Agaricus subrufescens) as a remediation agent. J. Environ. Manag. 2019, 245, 48–54. [Google Scholar] [CrossRef]
  23. Atila, F. Compositional changes in lignocellulosic content of some agro-wastes during the production cycle of shiitake mushroom. Sci. Hort. 2019, 245, 263–268. [Google Scholar] [CrossRef]
  24. Lane, D.J. 16S/23S rRNA sequencing. In Nucleic Acid Techniques in Bacterial Systematics; Stackebrandt, E., Goodfellow, M., Eds.; John Wiley & Sons: New York, NY, USA, 1991; pp. 115–148. [Google Scholar]
  25. Boratyn, G.M.; Schäffer, A.A.; Agarwala, R.; Altschul, S.F.; Lipman, D.J.; Madden, T.L. Domain enhanced lookup time accelerated BLAST. Biol. Direct. 2012, 17, 7–12. [Google Scholar] [CrossRef] [PubMed]
  26. Coster, W.D.; D’Hert, S.; Schultz, D.T.; Cruts, M.; Broeckhoven, C.V. NanoPack: Visualizing and processing long-read sequencing data. Bioinformatics 2018, 34, 2666–2669. [Google Scholar] [CrossRef] [PubMed]
  27. Wood, D.E.; Lu, J.; Langmead, B. Improved metagenomic analysis with Kraken 2. Genome. Biol. 2019, 20, 257. [Google Scholar] [CrossRef]
  28. Li, H. Minimap2: Pairwise alignment for nucleotide sequences. Bioinformatics 2018, 15, 3094–3100. [Google Scholar] [CrossRef]
  29. Breitwieser, F.P.; Salzberg, S.L. Pavian: Interactive analysis of metagenomics data for microbiome studies and pathogen identification. Bioinformatics 2020, 36, 1303–1304. [Google Scholar] [CrossRef] [PubMed]
  30. Bellettini, M.B.; Fiorda, F.A.; Maieves, H.A.; Teixeira, G.L.; Ávila, S.; Hornung, P.S.; Ribani, R.H. Factors affecting mushroom Pleurotus spp. Saudi. J. Biol. Sci. 2019, 26, 633–646. [Google Scholar] [CrossRef]
  31. Mcgee, C.F.; Byrne, H.; Irvine, A.; Wilson, J. Diversity and dynamics of the DNA and cDNA-derived bacterial compost communities throughout the Agaricus bisporus mushroom cropping process. Ann. Microbiol. 2017, 67, 751–761. [Google Scholar] [CrossRef]
  32. An, Q.; Wu, X.J.; Dai, Y.C. Comparative genomics of 40 edible and medicinal mushrooms provide an insight into the evolution of lignocellulose decomposition mechanisms. 3 Biotech 2019, 9, 1–15. [Google Scholar] [CrossRef]
  33. Zhang, H.L.; Wei, J.K.; Wang, Q.H.; Yang, R.; Gao, X.J.; Sang, Y.X.; Cai, P.P.; Zhang, G.Q.; Chen, Q.J. Lignocellulose utilization and bacterial communities of millet straw based mushroom (Agaricus bisporus) production. Sci. Rep. 2019, 9, 1–13. [Google Scholar] [CrossRef]
  34. Carrasco, J.; Tello, M.L.; Toro, M.; Tkacz, A.; Poole, P.; Pérez-Clavijo, M.; Preston, G. Casing microbiome dynamics during button mushroom cultivation: Implications for dry and wet bubble diseases. Microbiology 2019, 165, 611–624. [Google Scholar] [CrossRef] [PubMed]
  35. Carrasco, J.; García-Delgado, C.; Lavega, R.; Tello, M.L.; Toro, M.; Barba-Vicente, V.; Rodríguez-Cruz, M.S.; Sánchez-Martín, M.J.; Pérez, M.; Preston, G. Holistic assessment of the microbiome dynamics in the substrates used for commercial champignon (Agaricus bisporus) cultivation. Microb. Biotechnol. 2020, 13, 1933–1974. [Google Scholar] [CrossRef]
  36. Sun, S.; Li, F.; Xu, X.; Liu, Y.; Kong, X.; Chen, J.; Liu, T.; Chen, L. Study on the community structure and function of symbiotic bacteria from different growth and developmental stages of Hypsizygus marmoreus. BMC Microbiol. 2020, 20, 311. [Google Scholar] [CrossRef] [PubMed]
  37. Pecchia, J.; Cortese, R.; Albert, I. Investigation into the microbial community changes that occur in the casing layer during cropping of the white button mushroom, Agaricus bisporus. In Proceedings of the 8th International Conference on Mushroom Biology and Mushroom Products, New Delhi, India, 19–22 November 2014; ICAR-Directorate of Mushroom Research: Himachal Pradesh, India, 2014; pp. 309–313. [Google Scholar]
  38. Vos, A.M.; Heijboer, A.; Boschker, H.T.; Bonnet, B.; Lugones, L.G.; Wösten, H.A. Microbial biomass in compost during colonization of Agaricus bisporus. AMB Express 2017, 7, 1–7. [Google Scholar] [CrossRef] [PubMed]
  39. Liu, C.; Sheng, J.P.; Chen, L.; Zheng, Y.Y.; David, Y.W.L.; Yang, Y.; Xu, M.S.; Shen, L. Biocontrol activity of Bacillus subtilis isolated from Agaricus bisporus mushroom compost against pathogenic fungi. J. Agric. Food Chem. 2015, 63, 6009–6018. [Google Scholar] [CrossRef]
  40. Pandin, C.; Védie, R.; Rousseau, T.; Le Coq, D.; Aymerich, S.; Briandet, R. Dynamics of compost microbiota during the cultivation of Agaricus bisporus in the presence of Bacillus velezensis QST713 as biocontrol agent against Trichoderma aggressivum. Biol. Control 2018, 127, 39–54. [Google Scholar] [CrossRef]
  41. Ralph, N.; Andreja, D.P.; Philip, J.H.; Jemma, P.; Alison, R. Volatile C 8 compounds and pseudomonads influence primordium formation of Agaricus bisporus. Mycologia 2009, 101, 583–591. [Google Scholar] [CrossRef]
  42. Colauto, N.B.; Fermor, T.R.; Eira, A.F.; Linde, G.A. Pseudomonas putida stimulates primordia on Agaricus bitorquis. Curr. Microbiol. 2016, 72, 482–488. [Google Scholar] [CrossRef]
  43. Jarial, R.S.; Shandilya, T.R.; Kumad, J. Casing in Mushroom Beds—A review. Agric. Rev. 2005, 26, 261–271. [Google Scholar]
  44. Kim, M.K.; Math, R.K.; Cho, K.M.; Shin, K.J.; Kim, J.O.; Ryu, J.S.; Lee, Y.H.; Yun, H.D. Effect of Pseudomonas sp. P7014 on the growth of edible mushroom Pleurotus eryngii in bottle culture for commercial production. Bioresour Technol. 2008, 99, 3306–3308. [Google Scholar] [CrossRef]
  45. Cho, Y.S.; Kim, J.S.; Crowley, D.E.; Cho, B.G. Growth promotion of the edible fungus Pleurotus ostreatus by fluorescent pseudomonads. FEMS Microbiol. Lett. 2003, 218, 271–276. [Google Scholar] [CrossRef]
  46. Brurberg, M.B.; Eijsink, V.G.; Haandrikman, A.J.; Venema, G.; Nes, I.F. Chitinase B from Serratia marcescens BJL200 is exported to the periplasm without processing. Microbiology 1995, 141, 123–131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Kaomek, M.; Mizuno, K.; Fujinira, T.; Sriyotha, P.; Cairna, J.R.K. Cloning, expression and characterization of an antifungal chitinase from Leucaena leucocephla De Wit. Biosci. Biotechnol. Biochem. 2003, 67, 667–676. [Google Scholar] [CrossRef]
  48. Van Houdt, R.; Givskov, M.; Michiels, C.W. Quorum sensing in Serratia. FEMS Microbiol. Rev. 2007, 31, 407–424. [Google Scholar] [CrossRef] [PubMed]
  49. Matteoli, F.P.; Passarelli-Araujo, H.; Regis, J.A.R.; da Rocha, L.O.; de Souza, E.M.; Aravind, L.; Olivares, F.L.; Venancio, T.M. Genome sequencing and assessment of plant growth-promoting properties of a Serratia marcescens strain isolated from vermicompost. BMC Genom. 2018, 39, 750. [Google Scholar] [CrossRef]
  50. Das, M. Chitinase produced by Serratia marcescens SMG isolated from decomposed Volvariella volvacea. Afr. J. Microbiol. Res. 2011, 5, 3220–3222. [Google Scholar] [CrossRef]
  51. Sun, S.J.; Liu, Y.C.; Weng, C.H.; Sun, S.W.; Li, F.; Li, H.; Zhu, H. Cyclic dipeptides mediating quorum sensing and their biological effects in Hypsizygus marmoreus. Biomolecules 2020, 10, 298. [Google Scholar] [CrossRef]
  52. Schmidt, R.; Jager, V.D.; Zühlke, D.; Wolff, C.; Bernhardt, J.; Cankar, K.; Garbeva, P. Fungal volatile compounds induce production of the secondary metabolite Sodorifen in Serratia plymuthica PRI-2C. Sci. Rep. 2017, 7, 1–14. [Google Scholar] [CrossRef]
  53. Pardo-Giménez, A.; Pardo-González, J.E.; Zied, D.C. Casing materials and techniques in Agaricus bisporus cultivation. In Edible and Medicinal Mushrooms: Technology and Applications; Zied, D.C., Pardo-Giménez, A., Eds.; Wiley Blackwell: West Sussex, UK, 2017; pp. 149–174. [Google Scholar]
Figure 1. Steps involved in the preparation of mushroom-colonized compost extract (MCCE) and its addition over the casing of compost blocks for the production of Agaricus subrufescens. The fully colonized compost block with the mycelium was ground by hand and soaked in water for 1 h or 24 h (step 1), followed by filtration (step 2), dilution and storage of the MCCE in the refrigerator (steps 3 to 5). Each bottle obtained with 1 h or 24 h soaking time was then added to another colonized compost block in the casing on the 3rd, 6th, 9th and 15th days after casing (step 6). Full-length 16S rRNA genes of the 1 h and 24 h soaked MCCEs were amplified and sequenced using nanopore technology (step 7). A. subrufescens mushrooms were grown in a chamber under semi-controlled conditions (steps 8 and 9).
Figure 1. Steps involved in the preparation of mushroom-colonized compost extract (MCCE) and its addition over the casing of compost blocks for the production of Agaricus subrufescens. The fully colonized compost block with the mycelium was ground by hand and soaked in water for 1 h or 24 h (step 1), followed by filtration (step 2), dilution and storage of the MCCE in the refrigerator (steps 3 to 5). Each bottle obtained with 1 h or 24 h soaking time was then added to another colonized compost block in the casing on the 3rd, 6th, 9th and 15th days after casing (step 6). Full-length 16S rRNA genes of the 1 h and 24 h soaked MCCEs were amplified and sequenced using nanopore technology (step 7). A. subrufescens mushrooms were grown in a chamber under semi-controlled conditions (steps 8 and 9).
Jof 08 00995 g001
Figure 2. Sankey visualization of taxonomic classification at domain (D), phylum (P), class (C), order (O), family (F), genus (G) and species (S) levels based on the 16S rRNA gene amplicons from samples of mushroom-colonized compost extracts (MCCEs) soaked in water for either 1 h (A) or 24 h (B).
Figure 2. Sankey visualization of taxonomic classification at domain (D), phylum (P), class (C), order (O), family (F), genus (G) and species (S) levels based on the 16S rRNA gene amplicons from samples of mushroom-colonized compost extracts (MCCEs) soaked in water for either 1 h (A) or 24 h (B).
Jof 08 00995 g002
Figure 3. Most abundant phyla (A) and genera (B) of bacteria present in the Agaricus subrufescens mushroom-colonized compost extracts (MCCE) soaked in water for either 1 h or 24 h. Others figures refer to Bacteroidetes, BRC 1, Chloroflexi, Cyanocateria, OD 1 (genera) and Verrucomicrobia (phyla).
Figure 3. Most abundant phyla (A) and genera (B) of bacteria present in the Agaricus subrufescens mushroom-colonized compost extracts (MCCE) soaked in water for either 1 h or 24 h. Others figures refer to Bacteroidetes, BRC 1, Chloroflexi, Cyanocateria, OD 1 (genera) and Verrucomicrobia (phyla).
Jof 08 00995 g003
Figure 4. Agronomic parameters of mushrooms cultivated with the control and application of 1 h- and 24 h soaked MCCE over the casing, which (A) indicate biological efficiency, (B) precociousness, (C) number of mushroom, and (D) weight of mushroom.
Figure 4. Agronomic parameters of mushrooms cultivated with the control and application of 1 h- and 24 h soaked MCCE over the casing, which (A) indicate biological efficiency, (B) precociousness, (C) number of mushroom, and (D) weight of mushroom.
Jof 08 00995 g004
Figure 5. Biological efficiency of the Agaricus subrufescens mushroom cultivation after the application of 1 h or 24 h soaked MCCE on the 3rd, 6th, 9th, 12th and 15th days after casing. The red line corresponds to the control treatment. The highlighted red area includes the control mean value and standard deviation (20 ± 5%, n = 4).
Figure 5. Biological efficiency of the Agaricus subrufescens mushroom cultivation after the application of 1 h or 24 h soaked MCCE on the 3rd, 6th, 9th, 12th and 15th days after casing. The red line corresponds to the control treatment. The highlighted red area includes the control mean value and standard deviation (20 ± 5%, n = 4).
Jof 08 00995 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Iossi, M.R.; Palú, I.A.; Soares, D.M.; Vieira, W.G., Jr.; Alves, L.S.; Stevani, C.V.; Caitano, C.E.C.; Atum, S.V.F.; Freire, R.S.; Dias, E.S.; et al. Metaprofiling of the Bacterial Community in Colonized Compost Extracts by Agaricus subrufescens. J. Fungi 2022, 8, 995. https://doi.org/10.3390/jof8100995

AMA Style

Iossi MR, Palú IA, Soares DM, Vieira WG Jr., Alves LS, Stevani CV, Caitano CEC, Atum SVF, Freire RS, Dias ES, et al. Metaprofiling of the Bacterial Community in Colonized Compost Extracts by Agaricus subrufescens. Journal of Fungi. 2022; 8(10):995. https://doi.org/10.3390/jof8100995

Chicago/Turabian Style

Iossi, Matheus Rodrigo, Isabela Arruda Palú, Douglas Moraes Soares, Wagner G. Vieira, Jr., Lucas Silva Alves, Cassius V. Stevani, Cinthia E. C. Caitano, Samir V. F. Atum, Renato S. Freire, Eustáquio S. Dias, and et al. 2022. "Metaprofiling of the Bacterial Community in Colonized Compost Extracts by Agaricus subrufescens" Journal of Fungi 8, no. 10: 995. https://doi.org/10.3390/jof8100995

APA Style

Iossi, M. R., Palú, I. A., Soares, D. M., Vieira, W. G., Jr., Alves, L. S., Stevani, C. V., Caitano, C. E. C., Atum, S. V. F., Freire, R. S., Dias, E. S., & Zied, D. C. (2022). Metaprofiling of the Bacterial Community in Colonized Compost Extracts by Agaricus subrufescens. Journal of Fungi, 8(10), 995. https://doi.org/10.3390/jof8100995

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