- Research article
- Open Access
Molecular characterization of Pseudomonas from Agaricus bisporus caps reveal novel blotch pathogens in Western Europe
BMC Genomics volume 21, Article number: 505 (2020)
Bacterial blotch is a group of economically important diseases affecting the cultivation of common button mushroom, Agaricus bisporus. Despite being studied for more than a century, the identity and nomenclature of blotch-causing Pseudomonas species is still unclear. This study aims to molecularly characterize the phylogenetic and phenotypic diversity of blotch pathogens in Western Europe.
In this study, blotched mushrooms were sampled from farms across the Netherlands, United Kingdom and Belgium. Bacteria were isolated from symptomatic cap tissue and tested in pathogenicity assays on fresh caps and in pots. Whole genome sequences of pathogenic and non-pathogenic isolates were used to establish phylogeny via multi-locus sequence alignment (MLSA), average nucleotide identity (ANI) and in-silico DNA:DNA hybridization (DDH) analyses.
The known pathogens “Pseudomonas gingeri”, P. tolaasii, “P. reactans” and P. costantinii were recovered from blotched mushroom caps. Seven novel pathogens were also identified, namely, P. yamanorum, P. edaphica, P. salomonii and strains that clustered with Pseudomonas sp. NC02 in one genomic species, and three non-pseudomonads, i.e. Serratia liquefaciens, S. proteamaculans and a Pantoea sp. Insights on the pathogenicity and symptom severity of these blotch pathogens were also generated.
A detailed overview of genetic and regional diversity and the virulence of blotch pathogens in Western Europe, was obtained via the phylogenetic and phenotypic analyses. This information has implications in the study of symptomatic disease expression, development of diagnostic tools and design of localized strategies for disease management.
Commercial button mushroom cultivation relies heavily on the dynamic interactions between Agaricus bisporus and the casing soil microflora . The transformation of vegetative mycelium into a fruiting body is initiated by beneficial microbes in the casing soil [2,3,4]. However, the casing soil also introduces pathogenic microbes into mushroom farms, including blotch causing Pseudomonas species [5,6,7]. The humid and mesophilic conditions required for mushroom production are highly conducive to the enrichment and spread of such pathogens. Reliable identification and early detection is thus essential to avoid disease outbreaks.
The genus Pseudomonas is one of the most complex genera of Gram negative bacteria due to its large size of 114 species . They form a major proportion (~ 40%) of the total culturable bacteria obtained from casing soil in mushroom farms . While some of these are essential for stimulating the pinning of button mushrooms (e.g. P. putida) [1, 10], others are detrimental to crop health (e.g. P. tolaasii) [11, 12].
Bacterial blotch is a group of diseases that result in discolouration and disfiguration of mushroom caps in A. bisporus, due to fungal production of phenols and tyrosinases . This reduces the total marketable crop due to compromised aesthetic value, lowers the shelf-life post-harvest, and lessens the overall yields due to pin death. These aspects of bacterial blotch jointly lead to severe economic losses [14,15,16]. Various Pseudomonas species are the main causative agents of blotch diseases on mushroom caps .
P. tolaasii causes small sunken dark brown spots or lesions on the mushroom cap that are referred to as “brown blotch” [11, 12]. “P. reactans” is known to cause varying discoloration from dark to light, accompanied by a surface depression  and “P. gingeri” produces ginger coloured discolorations that are more spread out on the cap surface, called “ginger blotch” . Both of these species have not been formally described. P. agarici is the causative agent of “drippy gill” on A. bisporus and “yellow blotch” on oyster mushrooms (Pleurotus spp.), where it leads to relatively pale discolorations . Global reports also indicate the role of other Pseudomonas, such as P. costantinii, P. fluorescens and P. marginalis in bacterial blotch diseases, with large phenotypic variation within and across species .
Blotch pathogens can be considered as endemic to the casing soil, an artificially prepared growth media composed of peat and lime, that is added on top of the compost [5, 21]. They have been found on healthy crops at similar densities to that of diseased crops . It has thus been suggested that not just the pathogen density, but the composition of Pseudomonas species in the casing soil, especially the relative abundance of beneficial and disease-causing species, can be an important indicator for disease outbreaks . A deeper understanding of the beneficials and pathogens within the genus is hence necessary.
Bacterial blotch has been studied for over a century [22, 23], despite which the identity and nomenclature of blotch-causing Pseudomonas is still unclear. Recent molecular investigations clarify the taxonomy of some blotch pathogens [24,25,26,27]. However, knowledge on the identity, diversity and pathogenicity of mushroom-associated Pseudomonas species at the regional scale is still lacking. This information is instrumental for the development of localized strategies for diagnostics, disease control and breeding of varieties.
In this study, we isolated Pseudomonas from blotched mushrooms on farms in the Netherlands, United Kingdom and Belgium. We performed whole genome sequence analyses of pathogenic isolates to develop a deeper understanding of the genetic diversity among the pathogens in Western Europe. Some non-pathogenic isolates were also included in the study. The molecular characterization of these isolates provides insights into the phylogenetic relationships between beneficial and blotch-causing Pseudomonas species commonly associated with the button mushroom, A. bisporus.
Blotched mushrooms were sampled from commercial farms in the Netherlands, United Kingdom and Belgium for isolation of blotch-causing Pseudomonas species. Biopsies from symptomatic tissue of the cap surface (2 cm2 area) were made in sterilized Ringer’s solution , and homogenized in a polyethylene bag (Bioreba, Switzerland). The extract from each biopsy was dilution plated on King’s B medium . After incubation at 25 °C for 48 h, single colonies were picked and re-plated. In total, 161 single colonies of suspected Pseudomonas spp., that were fluorescent under UV light (365 nm), were plated to pure cultures by re-streaking on King’s B medium. One isolation was also made from a healthy mushroom that did not display visible blotch symptoms. A list of bacterial isolates is presented as Additional file 1.
All isolates were tested in an in vitro assay to check their pathogenicity. Bacterial strains were cultured in King’s B medium  at 25 °C for 24 h, and tested in the bioassay. Similarly sized cap surfaces (4–5 cm in diameter) of healthy mushrooms were placed on damp filter paper and inoculated with 20 μl of aqueous bacterial suspension of 106 colony-forming units (cfu) per ml from the isolate, and tested in replicates of three. The mushrooms were incubated under high humidity conditions for 72 h at 20 °C. The development of blotch symptoms on the cap surface was observed visually and photographed. The isolates were scored, between 0 and 3, with the ascending numbers referring to non-pathogenic, mild, moderate, and severe symptoms for bacterial blotch . Negative controls consisted of uninoculated mushroom caps and sterile water inoculated mushroom caps.
A selection of isolates were re-tested in pot assays. Mushrooms were grown in plastic pots (230 mm diameter × 220 mm depth) containing 4 kg of Phase III compost, spawn-run with the most commonly cultivated mushroom strain, Sylvan A15. The pots were cased with 1.3 L of casing soil (moist mixture of peat and sugar beet lime). The pots were watered with sterile water and incubated at 25 °C for 7 days. The room was then ventilated and the air temperature reduced to 18 °C and the relative humidity was maintained at 91–93%, until the end of the cultivation cycle. After 5 days, the casing soil in each pot was inoculated with 50 ml of aqueous bacterial suspension of 107 cfu/ml. The development of blotch symptoms on the mushrooms was recorded over two flushes and scored as above. The type of blotch symptoms (brown, ginger or others) was recorded and photographed. Negative controls consisted of casing soil inoculated with sterile water.
DNA extraction and sequencing
For NL and BE isolates, 250 mg of bacterial slime was picked from a pure culture on an agar plate, and used as starting material. Total DNA was extracted using Wizard Magnetic DNA Purification System for Food (Promega, United States) according to the manufacturer’s protocol, including the DNase-free RNAse treatment. Library construction was performed using Illumina Truseq Nano (Illumina, United States) with 1 μg of bacterial DNA. 125 bp paired-end sequencing of the DNA libraries was done using HiSeq2500 (Illumina, United States).
For UK isolates, a single colony was picked from each agar plate and extracted using the Qiagen DNeasy Blood and tissue kit following the manufacturer’s protocol. The DNA was quantified fluorometrically using a Quant-iT PicoGreen dsDNA Assay Kit (Thermo Fisher Scientific, United States) on the Infinite M200 PRO (Tecon, Switzerland) and then stored at − 80 °C for downstream processing. Library construction was performed using Illumina Nextera XT library preparation kit (Illumina, United States) with 0.8 ng of bacterial DNA. Sequencing of the DNA libraries was performed on the MiSeq (Illumina, United States) using the V3 Reagent Kit, generating 300 bp paired-end sequences.
The combined dataset included 68 newly generated genome sequences from bacteria isolated from symptomatic cap tissue, 30 sequences of mushroom-associated Pseudomonas species from a previous sampling  and 15 sequences of related strains obtained from NCBI (https://www.ncbi.nlm.nih.gov/). Quality control was performed on the raw reads prior to read mapping using CLC Genomics Workbench (QIAGEN, Germany). Adapter sequences were removed from the raw reads. Bases with Phred quality scores less than 20 based on a modified-Mott algorithm were trimmed. Raw reads greater than 1000 bp and less than 45 bp were discarded. Reads were trimmed to a final length of 125 bp (NL and BE isolates) and 300 bp (UK isolates). Trimmed reads were mapped to the reference genomes without masking. Non-specific matches were randomly mapped.
Determination of prokaryotic taxonomy
Multi-Locus Sequence Alignment (MLSA) with trimmed coding sequences of eleven barcoding genes from 13 reference strains were used to establish phylogeny between the isolates . Housekeeping genes were chosen as phylogenetic molecular markers based on several criteria. The genes had a single copy number, they code universally for ubiquitous proteins with housekeeping functions, are likely recalcitrant to the effects of horizontal gene transfer, are long enough (> 900 bp) to contain sufficient information, and can predict whole-genome relationships . Trimmed reads were mapped to the concatenated sequences of individual barcoding genes from multiple reference strains, using Map Reads to Reference 1.6 with a similarity and length fraction of 0.9 (CLC Genomics Workbench 11.0.2). Consensus DNA sequences were extracted from the mapping, and used for making phylogenetic trees with maximum likelihood and maximum parsimony methods . Graphics from phylogenetic trees were made in RStudio  using package ggtree .
Genome assemblies were performed on the trimmed reads using De Novo Assembly 1.4 with a minimum contig length of 200 bp (CLC Genomics Workbench 11.0.2). Legacy BLAST  based Average Nucleotide Identity (ANI) analysis  was performed on the contig sequences from the assembled genomes using pyani 0.2.9 . Similarity value of 95% was used as cut-off threshold for identification of a unique taxonomic group. The similarity values were used for phylogenetic analyses and to create graphics in RStudio . To clarify the taxonomy of isolates that did not cluster together with any of the reference strains in the ANI or the MLSA, a digital DNA:DNA hybridization  was performed using the Genome-genome distance calculator (GGDC) . A threshold of 70% for digital DNA:DNA hybridization and 1% for difference in percentage guanine-cytosine content were used for determination of species and subspecies boundaries via the Type (Strain) Genome Server . Phylogenetic trees were constructed from the alignment of the whole-genomes and their corresponding 16S rRNA sequences, using a GreedyWithTrimming algorithm on FastME 2.0 .
Pathogenicity of isolates
102 bacterial isolates were tested for their ability to cause bacterial blotch symptoms on fresh mushroom caps (Fig. 1). 6 out of the 17 strains that belonged to international culture collections could cause blotch symptoms. Out of the 85 bacterial isolates recovered from blotched mushroom tissue, 55 isolates caused mild to severe symptoms in the pathogenicity cap test. From this panel, the pathogenicity of 30 bacterial isolates and strains were also validated in the pot test, by inoculation of the pathogen in the casing soil (Fig. 2). The pot test and cap test gave similar results (Additional file 1). The pathogenicity of the isolates is further described in this text with reference to bacterial blotch only.
Whole genome sequences
In total, whole genome sequences from 113 bacteria were analysed. The consortium of sequenced bacteria contained 85 isolates from symptomatic mushroom tissue and 28 reference strains from international culture collections at LMG (Laboratory of Microbiology, Belgium), ATCC (American Type Culture Collection, United States) and NCPPB (National Collection of Plant Pathogenic Bacteria, United Kingdom). The total number of reads per isolate, averaged across the dataset, were 16,075,321 indicating good sequencing depth, with mean Phred score of 38 which suggests high sequence quality and mean GC content of 63%. Additional file 1 describes the sequence identifiers and their metadata.
Multi-locus sequence analysis
Eleven taxon-specific sequences that are stable with regard to rapid genetic mutations were selected as barcoding genes from known literature reports [8, 25, 42], namely, atpD, fusA, glnS, groeL, gyrB, ileS, recA, recN, rpoB, rpoD and uvrC (Table 1). Concatenated sequences of the barcoding genes from 13 reference strains of well-known mushroom-associated Pseudomonas species (Table 2) were used for MLSA, such that the variability within the barcoding genes predicts the overall whole genome relatedness . Two of these strains from the P. fluorescens group were added to the reference list based on preliminary data exploration. P. yamanorum is a psychrotolerant soil bacterium from the P. fluorescens group , that contains a paralog of the tolaasin gene fragment . Pseudomonas sp. NC02 is a recently isolated soil bacterium, which is closely related to P. yamanorum.
The largest clusters in the phylogenetic tree are comprised of “P. gingeri”, Pseudomonas sp. NC02 “P. reactans” and P. tolaasii, in decreasing order of size. “P. gingeri” and “P. reactans” isolates formed multiple clusters within the species (Fig. 3). A few isolates also mapped to barcoding genes from P. putida, P. agarici, P. veronii, P. costantinii and P. yamanorum. With the exception of “P. gingeri”, non-pathogenic isolates were also found in species clusters that contained pathogenic isolates.
None of the isolates mapped to reference strains of P. poae, P. protegens, P. fluorescens or P. syringae. 15 out of 113 isolates mapped non-specifically, with relatively low percentage identity, to multiple reference strains indicating the presence of other Pseudomonas species. 5 isolates did not map to any of the references and could potentially be non-pseudomonads. Within species clusters of the phylogenetic tree, the individual isolates had low numbers of substitutions per sequence site, indicating short evolutionary distances within species (Fig. 3).
Average nucleotide identity
Average nucleotide identity analysis recognized 32 unique bacterial phylotypes associated with the cap tissue of blotched button mushrooms in Western Europe (Fig. 4). These are phylogenetically distinct whole genome sequences differing at species level, based on a < 95% similarity cut-off for delineation. Similar to the MLSA, the largest numbers of blotch pathogens were identified as close relatives of Pseudomonas sp. NC02 or “P. gingeri”, although pathogenic isolates of P. tolaasii, P. costantinii, “P. reactans” and P. yamanorum were also discovered. Isolates that were unable to cause blotch were primarily identified as “P. reactans”, P. agarici, P. veronii and strains belonging to the same species as Pseudomonas sp. NC02. Within phylotypes, isolates did not cluster according to geographic region, year of outbreak or pathogenicity scores. Twelve unidentified phylotypes consisted of isolates with varying levels of symptom severity on fresh mushroom caps and did not contain any reference or type strains.
In-silico DNA:DNA hybridization
Five non-pseudomonad isolates were identified from in-silico DDH of the whole genome sequences (Table 3). Two isolates that cause moderate blotch on fresh caps, C2002 (phylotype 15) and C7002 (phylotype 17) were identified as Serratia liquefaciens and Serratia proteamaculans respectively (Fig. 5). Non-pathogenic isolate P7753 (phylotype 30) was closely related to Brevundimonas bullata. Non-pathogenic isolate P7760 (phylotype 31) was closely related to Cedecea neteri. Moderate blotch-causing isolate B9002 (phylotype 13) was related to reference genomes of multiple Pantoea species.
Among the Pseudomonas, three isolates that cause blotch on fresh caps (phylotype 12) belonged to P. edaphica and severe blotch-causing isolate IPO3765 (phylotype 25) was identified as P. salomonii (Fig. 6). Several Pseudomonas spp. did not hybridize sufficiently (> 70%) with any of the type strains or share a similar %GC content (< 1%). Seven of these isolates (phylotypes 4 and 26) were closely related to P. yamanorum LMG 27247T. Isolates P7548 (phylotype 9) and B6002 (phylotype 11), clustered together as the same genomic species, and were closely related to P. fluorescens DSM 50090T (> 60% dDDH). Isolate D2002 (phylotype 21) was also closely related to P. fluorescens. Isolate B3002 (phylotype 16) did not map sufficiently with the type strain of “P. reactans”, in contrast to the MLSA and ANI results.
The taxonomy of several reference strains from culture collections have been corrected based on the combined results from MLSA, ANI and dDDH analyses. ATCC 51312, ATCC 51311, LMG 2343 were formerly identified as P. tolaasii, but they share less than 95% genome similarity with P. tolaasii genome cluster (phylotype 6). They also do not map to barcoding genes of P. tolaasii LMG 2342T. Non-pathogenic strains ATCC 51312 and ATCC 51311 were instead identified close relatives of Pseudomonas sp. NC02 (phylotype 2) in the MLSA and ANI, and LMG 2343 (phylotype 12) was identified as P. edaphica in the dDDH. Non-pathogenic isolates P7774, P7753 and P7760 were formerly described as P. veronii, P. syringae and P. agarici, respectively, in internal culture collections. Based on ANI and dDDH results from this study they were re-identified as “P. reactans” (phylotype 16), close relative of Brevundimonas bullata (phylotype 30) and close relative of Cedecea neteri (phylotype 31) respectively.
Known blotch pathogens
All strains of P. tolaasii (phylotype 6) isolated from blotch outbreaks in the Netherlands were found to cause moderate or severe brown blotch symptoms, with the exception of IPO 3746. They also map to barcoding genes of P. tolaasii LMG 2342T in the MLSA. Two pathogenic isolates from the United Kingdom, 21,815,971 and 21,815,972 (phylotype 3) that cause severe blotch and pitting in both cap and pot tests, clustered together with the blotch-causing reference strain of P. costantinii, LMG 22119T in the MLSA and ANI. Three isolates from the “P. reactans” clusters, B3002, C5002, IPO 3751 (phylotypes 9 and 16), were also pathogenic. All 26 isolates of “P. gingeri” obtained from blotch outbreaks in the Netherlands, United Kingdom and Belgium were able to cause blotch, with varying symptom severity in the pathogenicity tests performed. All of them mapped to barcoding genes of “P. gingeri” LMG 5327T. They form 3 MSLA clusters, but 5 distinct ANI clusters (phylotypes 1, 5, 14, 22 and 24), indicating that multiple closely-related species can cause ginger blotch.
New blotch pathogens
Several isolates presented in this study are reported as blotch-pathogens for the first time. 19 isolates (phylotype 2) with varying severity of blotch symptoms were identified as belonging to the same species as Pseudomonas sp. NC02. Two blotch-causing isolates (phylotype 10) from the Netherlands were identified as P. yamanorum. They were also closely related to the Pseudomonas sp. NC02 isolates. Three pathogenic isolates (phylotype 12) from Dutch blotch outbreaks were identified as P. edaphica. A severely pathogenic isolate (IPO3765) from Netherlands was recognized as P. salomonii.
Four ANI clusters (phylotypes 4, 9, 21 and 26) contained pathogenic Pseudomonas isolates without a reference or type strain. Phylotype 4 contained only one severe blotch-causing isolate, A5002. Phylotypes 9 and 21 contained only a single isolate each, both of which (B6002 and D2002) cause brown blotch symptoms. Phylotype 26 contained two unidentified isolates, IPO 3778 and IPO 3779 that cause mild-blotch on fresh caps.
Unidentified Pseudomonas species
10 isolates demonstrating a range of blotch severity, from non-pathogenic to severe symptoms, remained unidentified. They mapped to barcoding genes from multiple reference strains in the MLSA, formed 7 distinct clusters (phylotypes 4, 9, 11, 21, 26, 30 and 31) in the ANI analysis (Fig. 4), and were closely but insufficiently (< 70%) related to any other species in the dDDH (Table 3). Their contigs mapped in-specifically to multiple reference strains, when compared via BLAST against the NCBI database. Contamination of these isolates during laboratory work was ruled out by confirming the presence of single-copies of the barcoding genes, and low genome duplication in the range of 1.50–2.84% (Additional file 2). These could potentially be new species that have yet to be discovered.
Five non-pathogenic isolates, P7787, P8021, P774, IPO 3748 and IPO3750, mapped to barcoding genes from “P. reactans” LMG 5329T in the MLSA, but were split over phylotype 8 and 16 in the ANI analysis. All isolates from phylotype 19 did not cause blotch and mapped to barcoding genes of P. agarici LMG 2112T, better known as the causative agent of drippy gill. Fluorescent Pseudomonas spp. such as P. putida (phylotype 23, 29, 30), P. veronii (phylotype 27) and P. poae (phylotype 20) that are usually considered as beneficial organisms, were also confirmed to be non-pathogenic.
Genetic diversity of blotch-causing Pseudomonas
The established phylogeny of bacterial isolates associated to symptomatic cap tissue is largely consistent between the MLSA, ANI and dDDH analyses. This substantiates the conclusions made about the identity of blotch pathogens, and the related taxonomic corrections. “P. gingeri”, P. tolaasii, “P. reactans” and P. costantinii are known causative agents of bacterial blotch [4, 11, 17, 24]. In addition, several novel blotch pathogens were identified in these analyses, namely P. yamanorum, P. edaphica, P. salomoni and strains clustering together with Pseudomonas sp. NC02, as the same genomic species.
P. yamanorum is a recently described psychrotolerant member of the P. fluorescens group that was found in Argentinian soil under cold and humid climatic conditions . In a recent study, the tolaasin gene fragment of a few brown blotch-causing isolates in Turkey were also found to share homology with that of P. yamanorum . Pseudomonas sp. NC02 is a recently isolated soil bacterium about which not much is known . P. edaphica is another newly identified fluorescent Pseudomonas that was isolated from rhizosphere soil, and is closely related to P. brenneri . P. salomonii was first characterized as a pathogen of ‘Café au lait’ disease on garlic , and was later also recovered from ready-to-eat vegetables . This is the first report of the blotch-causing ability of this pathogen on mushrooms. Its pathogenicity in the cap test was also verified by pot tests.
“P. gingeri” and “P. reactans” are both invalidly named species known to cause bacterial blotch [17, 18, 50]. In this study, all isolates previously classified as “P. gingeri” and “P. reactans” map to barcoding genes of their respective reference strain, although they form multiple distinct clusters in the MLSA as observed from the phylogenetic tree (Fig. 3). They also split into multiple phylotypes in the ANI analysis instead of clustering together (Fig. 4). Delineation of these phylotypes on a genome level (< 95% similarity) despite mapping to barcoding references from the same species, perhaps indicates the existence of a species complex with multiple taxonomically-related species, as opposed to other blotch pathogens which cluster as single phylotypes.
On the basis of genomic and phylogenetic analyses such as MLSA, ANI and dDDH, this study also reports four blotch-causing Pseudomonas spp. (phylotypes 4, 9, 21 and 26) which do not belong to any of the existing species within the genus Pseudomonas. These are as yet unidentified, and could potentially be new species. The taxonomic assignment of these novel blotch-causing phylotypes as new species requires a polyphasic approach involving additional phenotypic and chemotaxonomic data such as fatty acid methyl esters, polar lipids and respiratory quinones, etc. .
Pathogenicity and phylogeny
Koch’s postulates were evaluated to identify causative agents of bacterial blotch . Bacteria were isolated from symptomatic mushroom tissue and successfully grown to pure culture. In every case, blotch symptoms could be reproduced when isolates were inoculated on mushroom caps. The 4th postulate of Koch, i.e. re-isolation of pathogen from the experimental host and its identification as originally inoculated pathogen, was not fulfilled. The mushroom cap surfaces also harbour endemic Pseudomonas spp. that are closely related to the isolates, and have similar characteristics and colony morphology [53, 54]. Due to lack of diagnostic methods, the inoculated isolates cannot be easily differentiated from the endemic microflora. However, the disease developement can be attributed to the inoculated pathogens with certainty, due to the presence of multiple mock-inoculated negative controls which did not develop any blotch symptoms.
24, 8, 39 and 28 isolates were confirmed to cause non-pathogenic, mild, moderate and severe blotch symptoms respectively, on fresh mushroom caps in an in-vitro pathogenicity test. Bacterial isolates did not cluster together according to their pathogenicity scores in the MLSA or ANI analysis. Often blotch-causing and non-pathogenic isolates were part of the same phylotype or branch (Figs. 3 and 4). This elaborates the need for a pan-genomic analysis to identify the pathogenicity determinants of bacterial blotch diseases. Fluorescent Pseudomonas share a rather small core genome of 1344-protein coding genes . The size and diversity of the pan-genome is thus determined by secondary metabolites biosynthesis clusters that can be strain specific and responsible for pathogenicity . Pan-genomic elements shared exclusively by pathogenic isolates can be ideal targets for the further development of DNA-based diagnostic tools.
Regional diversity of blotch-causing Pseudomonas
Within species, bacterial isolates did not cluster together according to their geographic region or year of outbreak. The two Belgian isolates (IPO 3781 and IPO 3782), both caused mild ginger blotch and clustered together with those of the Netherlands. Two Dutch isolates of “P. gingeri” (IPO 3757 and IPO 3756) that caused severe ginger blotch, clustered differently from each other in phylotypes 14 and 24 respectively. The phylogenetic distribution of the common pathogens independently from their region of outbreak indicates that Western Europe shares a well-mixed pathogen pool.
However, a few potentially region-specific blotch pathogens were also recovered in this study. P. yamanorum, P. edaphica and P. salomonii isolates have so far only been found in the Netherlands. P. costantinii isolates in this study, 21,815,971 and 21,815,972, which cause brown blotch and pitting, have been found only in the United Kingdom, although other blotch-causing P. costantinii have been reported in Finland too . It is also important to consider that very few isolates from these species were found in this study. This regional exclusivity could also be an artefact of low sampling due to the low presence of these species.
Blotch causing non-pseudomonads
Two Serratia species, S. liquefaciens, (C2002, phylotype 15) and S. proteamaculans, (C7002, phylotype 17) were recovered from blotched caps in this study. Serratia species are present in casing soil, compost, and on mushroom caps [53, 58,59,60]. Serratia liquefaciens is a known agent of yellow blotch on the caps of A. bisporus . They cause browning of the A. bisporus sporophore  by secreting chitinolytic enzymes and chitin-binding proteins . They also produce a surfactant, serwettin, which enables the bacterium to colonise fresh areas of the mushroom cap, as tolaasin does [64, 65]. P. fluorescens also uses a similar mechanisms to spread on broccoli heads .
A moderately pathogenic Pantoea isolate (B9002, phylotype 13) was also identified in this study. Pantoea species have been previously isolated from caps of A. bisporus , and are frequently recovered from soil and water environments too [67, 68]. This gram negative bacterium is an agent of soft rot disease on the stipes and pileus of Pleurotus eryngii , although it has not yet been reported to be a pathogen of A. bisporus. Pantoea species from waste mushroom beds are also used for production of plant auxin, indole-acetic acid, for phosphate solubilization .
A non-pathogenic Cedecea isolate (P7760, phylotype 31), closely related to C. neteri, was also found on a blotched mushroom cap. Cedecea were one of the most common genera recovered from blotched tissue of A. brazilensis  and is often present on A. bisporus also . Although Cedecea neteri is an agent of soft rot in Pholiota nameko  and yellow sticky disease in Flammulina velutipes , only Cedecea davisae strains have been reported to cause blotch in A. bisporus .
Beneficial bacteria on symptomatic cap tissue
Non-pathogenic and potentially beneficial bacteria were also recovered from cap tissue of blotched mushrooms. P. putida (P7771, P7772, phylotype 27) and P. veronii (P7765, phylotype 32) isolates were confirmed to be non-pathogenic on fresh mushroom caps. P. putida, P. poae and P. veronii are known to stimulate primordia formation of A. bisporus by the removal of volatile C8 compounds [75, 76]. The relative abundance of these beneficial bacteria on blotched caps is yet to be explored, in relation to disease expression by blotch causing bacteria. Brevundimonas isolate, P7753 (phylotype 30), closely related to B. bullata was also discovered on symptomatic cap tissue. During the cultivation cycle, increased Brevundimonas populations have been observed in the casing soils from A. bisporus and P. ostreatus farms [58, 77, 78], and Brevundimonas sp. can be thus expected to also colonize mushroom caps. It has also been reported as a symbiont on the apothecia of the cup fungus, Scutellinia scutellata .
This study provides a detailed overview of the regional, genetic and phenotypic diversity of bacterial blotch pathogens in Western Europe. It describes the presence of known pathogens inhabiting blotched mushroom caps: “P. gingeri”, P. tolaasii, P. costantinii and “P. reactans”. It also reports seven new pathogens of bacterial blotch: four Pseudomonas species, P. yamanorum, P. edaphica, P. salomonii and strains belonging to the same species as Pseudomonas sp. NC02, and three non-pseudomonads, S. liquefaciens, S. proteamaculans and a Pantoea sp. Their epidemiology and aetiology deserve further attention. This molecular investigation of bacterial species that colonize blotched mushroom caps allows further study of the mechanisms and microbial ecology of symptomatic disease expression on the A. bisporus sporophore. It also highlights the need for development of diagnostic assays against these newly discovered pathogens. Several new blotch-causing phylotypes discovered in this study are as yet unidentified and require additional chemotaxic data for their taxonomic assignment as new species.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the NCBI repository under BioProject number PRJNA607442.
Hayes WA, Randle PE, Last FT. The nature of the microbial stimulus affecting sporophore formation in Agaricus bisporus (Lange) sing. Ann Appl Biol. 1969;64(1):177–87.
Eger G. The" Halbschalentest"-a simple method for testing-casing materials. MGA Bull Grow Assoc. 1962;148:159–68.
Thomas DS, Mullins JT, Block SS. Involvement of volatiles in mushroom sporophore initiation. Mushroom Grow Assoc Bull. 1964;177:429–33.
Preece TF, Wong WC. Quantitative and scanning electron microscope observations on the attachment of Pseudomonas tolaasii and other bacteria to the surface of Agaricus bisporus. Physiol Plant Pathol. 1982;21(2):251–7.
Wong WC, Preece TF. Pseudomonas tolaasii in mushroom crops: a note on primary and secondary sources of the bacterium on a commercial farm in England. J Appl Bacteriol. 1980;49(2):305–14.
Fletcher JT, White PF, Gaze RH. Mushrooms: pest and disease control 2nd edition. Andover: UK Intercept; 1989. p. 174.
Fletcher JT, Gaze RH. Mushroom pest and disease control: a colour handbook: CRC Press; 2007.
Gomila M, Peña A, Mulet M, Lalucat J, García-Valdés E. Phylogenomics and systematics in Pseudomonas. Front Microbiol. 2015;6:214.
Godfrey SAC. Molecular investigation of pseudomonads causative of Agaricus bisporus blotch disease in New Zealand mushroom farms; 2003.
Hume DP, Hayes WA. The production of fruit-body primordia in Agaricus bisporus (Lange) sing. On agar media. Mushroom Sci. 1972;7:527–32.
Tolaas AG. A bacterial disease of cultivated mushrooms. Phytopathology. 1915;5(1):U51–5.
Paine SG. A brown blotch disease of cultivated mushrooms. Ann Appl Biol. 1919;5:206–19.
Brennan M, Le Port G, Gormley R. Post-harvest treatment with citric acid or hydrogen peroxide to extend the shelf life of fresh sliced mushrooms. LWT-Food Sci Technol. 2000;33(4):285–9.
Burton KS. The effects of pre-and post-harvest development on mushroom tyrosinase. J Hortic Sci. 1988;63(2):255–60.
Burton KS, Noble R. The influence of flush number, bruising and storage temperature on mushroom quality. Postharvest Biol Technol. 1993;3(1):39–47.
Soler-Rivas C, Jolivet S, Arpin N, Olivier JM, Wichers HJ. Biochemical and physiological aspects of brown blotch disease of Agaricus bisporus. FEMS Microbiol Rev. 1999;23(5):591–614.
Wells JM, Sapers GM, Fett WF, Butterfield JE, Jones JB, Bouzar H, et al. Postharvest discoloration of the cultivated mushroom Agaricus bisporus caused by Pseudomonas tolaasii, P.’reactans’, and P.’gingeri’. Phytopathology. 1996;86(10):1098–104.
Wong WC, Fletcher JT, Unsworth BA, Preece TF. A note on ginger blotch, a new bacterial disease of the cultivated mushroom, Agaricus bisporus. J Appl Bacteriol. 1982;52(1):43–8.
Young JM. Drippy gill: a bacterial disease of cultivated mushrooms caused by Pseudomonas agarici n. sp. New Zeal J Agric Res. 1970;13(4):977–90.
Abou-Zeid MA. Pathogenic variation in isolates of Pseudomonas causing the brown blotch of cultivated mushroom, Agaricus bisporus. Brazilian J Microbiol. 2012;43(3):1137–46.
Nair NG, Bradley JK. Mushroom blotch bacterium during cultivation. Mushroom J. 1980;90:201–3.
Elphinstone J, Noble R. Mushrooms: Identification, detection and control of different Pseudomonas species causing bacterial blotch symptoms. Kenilworth: UK Agriculture and Horticulture Development Board (AHDB) Horticulture, Report M063; 2018.
Osdaghi E, Martins S, Ramos-Sepulveda L, Vieira FR, Pecchia J, Beyer DM, et al. 100 Years since Tolaas: Bacterial Blotch of Mushrooms in the 21st Century. Plant Dis. 2019;103(11):2714–32.
Munsch P, Alatossava T. The white-line-in-agar test is not specific for the two cultivated mushroom associated pseudomonads, Pseudomonas tolaasii and Pseudomonas “reactans”. Microbiol Res. 2002;157(1):7–11.
van der Wolf JM, Kastelein P, Krijger MC, Hendriks MJA, Baars JJP, Amsing JGM, et al. Characterization of Pseudomonas species causing brown blotch of Agaricus bisporus. Int Soc Mushroom Sci. 2016;19:104–8.
Godfrey SAC, Marshall JW, Klena JD. Genetic characterization of Pseudomonas ‘NZI7’–a novel pathogen that results in a brown blotch disease of Agaricus bisporus. J Appl Microbiol. 2001;91(3):412–20.
Godfrey SAC, Harrow SA, Marshall JW, Klena JD. Characterization by 16S rRNA sequence analysis of pseudomonads causing blotch disease of cultivated Agaricus bisporus. Appl Environ Microbiol. 2001;67(9):4316–23.
Marino PL, Sutin KM, del Sol Jaquotot MJ. El Libro de la UCI 3rd Edition; 2008.
King EO, Ward MK, Raney DE. Two simple media for the demonstration of pyocyanin and fluorescin. J Lab Clin Med. 1954;44(2):301–7 Available from: https://www.translationalres.com/article/0022-2143(54)90222-X/abstract.
Zeigler DR. Gene sequences useful for predicting relatedness of whole genomes in bacteria. Int J Syst Evol Microbiol. 2003;53(6):1893–900.
Adékambi T, Butler RW, Hanrahan F, Delcher AL, Drancourt M, Shinnick TM. Core gene set as the basis of multilocus sequence analysis of the subclass Actinobacteridae. PLoS One. 2011;6(3):e14792.
Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci. 2004;101(30):11030–5.
Team RC. R: A language and environment for statistical computing; 2013.
Yu G, Smith DK, Zhu H, Guan Y, Lam TT. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol Evol. 2017;8(1):28–36.
Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. DNA–DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol. 2007;57(1):81–91.
Richter M, Rosselló-Móra R. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci. 2009;106(45):19126–31.
Pritchard L, Cock P, Esen Ö. pyani v0. 2.8: average nucleotide identity (ANI) and related measures for whole genome comparisons; 2019.
Auch AF, von Jan M, Klenk H-P, Göker M. Digital DNA-DNA hybridization for microbial species delineation by means of genome-to-genome sequence comparison. Stand Genomic Sci. 2010;2(1):117.
Meier-Kolthoff JP, Auch AF, Klenk H-P, Göker M. Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinformatics. 2013;14(1):60.
Meier-Kolthoff JP, Göker M. TYGS is an automated high-throughput platform for state-of-the-art genome-based taxonomy. Nat Commun. 2019;10(1):2182.
Lefort V, Desper R, Gascuel O. FastME 2.0: a comprehensive, accurate, and fast distance-based phylogeny inference program. Mol Biol Evol. 2015;32:2798–800. https://doi.org/10.1093/molbev/msv150.
Glaeser SP, Kämpfer P. Multilocus sequence analysis (MLSA) in prokaryotic taxonomy. Syst Appl Microbiol. 2015;38(4):237–45.
Arnau VG, Sánchez LA, Delgado OD. Pseudomonas yamanorum sp. nov., a psychrotolerant bacterium isolated from a subantarctic environment. Int J Syst Evol Microbiol. 2015;65(2):424–31.
Basim E, Basim H. Characterization of Pseudomonas tolaasii isolates, causing Brown blotch disease of mushroom (Agaricus bisporus). Int J Agric Biol. 2018;20(12):2721–30.
Holland BR, Huber KT, Dress A, Moulton V. δ plots: a tool for analyzing phylogenetic distance data. Mol Biol Evol. 2002;19(12):2051–9.
Cerra J, Donohue H, Kral A, Oser M, Rostkowski L, Zappia L, et al. Complete genome sequence of Pseudomonas sp. strain NC02, isolated from soil. Genome Announc. 2018;6(7):e00033–18.
Ramírez-Bahena MH, Salazar S, Santín PJ, Sánchez-Rodríguez JA, Fernández-Pascual M, Igual JM, et al. Pseudomonas edaphica sp. nov., isolated from rhizospheric soil of Cistus ladanifer L. in Spain. Int J Syst Evol Microbiol. 2019;69(10):3141–7.
Gardan L, Bella P, Meyer J-M, Christen R, Rott P, Achouak W, et al. Pseudomonas salomonii sp. nov., pathogenic on garlic, and Pseudomonas palleroniana sp. nov., isolated from rice. Int J Syst Evol Microbiol. 2002;52(6):2065–74.
Caldera L, Franzetti LV, Van Coillie E, De Vos P, Stragier P, De Block J, et al. Identification, enzymatic spoilage characterization and proteolytic activity quantification of Pseudomonas spp. isolated from different foods. Food Microbiol. 2016;54:142–53.
Iacobellis NS, Cantore PL. Pseudomonas “reactans” a new pathogen of cultivated mushrooms. In: Pseudomonas syringae and related pathogens. Dordrecht: Springer; 2003. p. 595–605.
Kämpfer P, Glaeser SP. Prokaryotic taxonomy in the sequencing era–the polyphasic approach revisited. Environ Microbiol. 2012;14(2):291–317.
Fulton MM. A simple lab exercise demonstrating Koch’s postulates. Am Biol Teach. 1981;43(7):394–6.
Doores S, Kramer M, Beelman R. Evaluation and bacterial populations associated with fresh mushrooms (AgaricusBisporus). In: Developments in crop science: Elsevier; 1987. p. 283–94.
Cutri SS, Macauley BJ, Roberts WP. Characteristics of pathogenic non-fluorescent (smooth) and non-pathogenic fluorescent (rough) forms of Pseudomonas tolaasii and Pseudomonas gingeri. J Appl Bacteriol. 1984;57(2):291–8.
Garrido-Sanz D, Meier-Kolthoff JP, Göker M, Martin M, Rivilla R, Redondo-Nieto M. Genomic and genetic diversity within the Pseudomonas fluorescens complex. PLoS One. 2016;11(2):e0150183.
Loper JE, Hassan KA, Mavrodi DV, Davis EW II, Lim CK, Shaffer BT, et al. Comparative genomics of plant-associated Pseudomonas spp.: insights into diversity and inheritance of traits involved in multitrophic interactions. PLoS Genet. 2012;8(7):e1002784.
Munsch P, Alatossava T, Marttinen N, Meyer J-M, Christen R, Gardan L. Pseudomonas costantinii sp. nov., another causal agent of brown blotch disease, isolated from cultivated mushroom sporophores in Finland. Int J Syst Evol Microbiol. 2002;52(6):1973–83.
Mcgee CF, 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(11):751–61.
Jarial RS, Shandilya TR, Jarial K. Casing in mushroom beds-a review. Agric Rev Res Commun Cent INDIA. 2005;26(4):261.
Reyes JE, Venturini ME, Oria R, Blanco D. Prevalence of Ewingella americana in retail fresh cultivated mushrooms (Agaricus bisporus, Lentinula edodes and Pleurotus ostreatus) in Zaragoza (Spain). FEMS Microbiol Ecol. 2004;47(3):291–6.
Sivanesan D. Diversity among bacteria causing blotch disease on the commercial mushroom, Agaricus bisporus; 2003.
Grewal PS. Influence of bacteria and temperature on the reproduction of Caenorhabditis elegans (Nematoda: Rhabditidae) infesting mushrooms (Agaricus bisporus). Nematologica. 1991;37(1–4):72–82.
Watanabe T, Kimura K, Sumiya T, Nikaidou N, Suzuki K, Suzuki M, et al. Genetic analysis of the chitinase system of Serratia marcescens 2170. J Bacteriol. 1997;179(22):7111–7.
Matsuyama T, Kaneda K, Nakagawa Y, Isa K, Hara-Hotta H, Yano I. A novel extracellular cyclic lipopeptide which promotes flagellum-dependent and-independent spreading growth of Serratia marcescens. J Bacteriol. 1992;174(6):1769–76.
Matsuyama T, Murakami T, Fujita M, Fujita S, Yano I. Extracellular vesicle formation and biosurfactant production by Serratia marcescens. Microbiology. 1986;132(4):865–75.
Hildebrand PD. Surfactant-like characteristics and identity of bacteria associated with broccoli head rot in Atlantic Canada. Can J Plant Pathol. 1989;11(3):205–14.
Brenner DJ, Family I. Enterobacteriaceae. Bergey’s Man Syst Bacteriol. 1984;1:408–516.
Holt JG, Krieg NR, Sneath PHA. Bergey’s manual of determinative bacterology; 1994.
Kim MK, Ryu JS, Lee YH, Yun HD. First report of Pantoea sp. induced soft rot disease of Pleurotus eryngii in Korea. Plant Dis. 2007;91(1):109.
Walpola BC, Noh J-G, Kim CK, Kyung K-C, Kong W-S, Yoon M-H. Optimization of indole-3-acetic production by phosphate solubilization bacteria isolated from waste mushroom bed of Agaricus bisporus. 한국버섯학회지. 2013;11(2):53–62.
Jhune C-S, Yun H-S, Park Y-J, Weon H-Y, Yoo Y-B, Lee C-J, et al. Studies on the outbreak cause of inner tissue browning and malformed symptoms on fruiting body of Agaricus brazilensis. J Mushroom. 2009;7(4):168–72.
Rossouw W, Korsten L. Cultivable microbiome of fresh white button mushrooms. Lett Appl Microbiol. 2017;64(2):164–70.
Yan JJ, Chen BZ, Wang RQ, Liu F, Wei Q, Xie BG, et al. First report of Cedecea neteri-induced soft rot disease of Pholiota nameko in China. Plant Dis. 2018;102(10):2025.
Yan JJ, Liu YY, Wang RQ, Mukhtar I, Liu F, Lin ZY, et al. First report of Cedecea neteri causing yellow sticky disease in Flammulina velutipes in China. Plant Dis. 2019;103(5):1014.
Fermor T, Lincoln S, Noble R, Dobrovin-Pennington A, Colauto N. Microbiological properties of casing. Mushroom Sci. 2000;15(1):447–54.
Noble R, Dobrovin-Pennington A, Hobbs PJ, Pederby J, Rodger A. Volatile C8 compounds and pseudomonads influence primordium formation of Agaricus bisporus. Mycologia. 2009;101(5):583–91.
Nazareth SA, Surridge-Talbot AKJ, Korsten L. Bacterial profiling of casing materials for white button mushrooms (Agaricus bisporus) using denaturing gradient gel electrophoresis; 2010.
Cho YS, Weon HY, Joh JH, Lim JH, Kim KY, Son ES, et al. Effect of casing layer on growth promotion of the edible mushroom Pleurotus ostreatus. Mycobiology. 2008;36(1):40–4.
Giordano R, Nardi JB, Bee CM, Miller LA, Raja H. Bacterial symbionts that inhabit apothecia of the cup fungus Scutellinia scutellata. Nov Hedwigia. 2013;97(1–2):1–18.
We would like to thank P Kastelein, B Brankovics, J Baars, H Pufal, B Carter, A Dobrovin-Pennington and J Vicente for their support.
The project in Netherlands was granted by the Topsector ‘Horticulture and propagation material’ of the Dutch Ministry of Economic Affairs. It was conducted under grant number KV1509–050 in collaboration with members of the Dutch and Belgian mushroom supply chain. Work in the United Kingdom was funded by the Agriculture and Horticulture Development Board under project number M063. The funders had no role in research design, data generation, analysis, and preparation of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Isolates and reference strains that were sequenced, with metadata on initial identity, pathogenicity, source, region of outbreak, year of outbreak, MLSA, ANI and dDDH.
Barcoding genes and genome duplication in Pseudomonas isolates that map in specifically and remain unidentified. It describes the number of copies of nine barcoding genes found in the assembled genome, and the overall percentage of genome duplication in the assembly.
About this article
Cite this article
Taparia, T., Krijger, M., Haynes, E. et al. Molecular characterization of Pseudomonas from Agaricus bisporus caps reveal novel blotch pathogens in Western Europe. BMC Genomics 21, 505 (2020). https://doi.org/10.1186/s12864-020-06905-3
- Multilocus sequence alignment
- Average nucleotide identity
- In-silico DNA DNA hybridization
- “Pseudomonas gingeri”
- Pseudomonas sp. NC02
- Serratia spp.
- Pseudomonas salomonii
- Pseudomonas yamanorum
- Pseudomonas edaphica
- Bacterial blotch
- Pot test
- Cap test