Viral Interactions and Pathogenesis during Multiple Viral Infections in Agaricus bisporus.

Viral interactions during multiple viral infections were examined in Agaricus bisporus cultures harboring 9 viruses (comprising 18 distinct viral RNAs) by statistically analyzing their relative abundance in fruitbodies. Four clusters of viral RNA were identified that suggested synergism and coreplication. Pairwise correlations revealed negative and positive correlations between clusters, indicating further synergisms and an antagonism involving a group containing a putative hypovirus and four nonhost ORFan RNAs (RNAs with no similarity to known sequences) possibly acting as defective interfering RNAs. The disease phenotype was observed in 10 to 15% of the fruitbodies apparently randomly located among asymptomatic fruitbodies. The degree of symptom expression consistently correlated with the levels of the multipartite virus AbV16. Diseased fruitbodies contained very high levels of AbV16 and AbV6 RNA2; these levels were orders of magnitude higher than those in asymptomatic tissues and were shown statistically to be discretely higher populations of abundance, indicating an exponential shift in the replicative capacity of the virus. High levels of AbV16 replication were specific to the fruitbody and not found in the underlying mycelium. There appeared to be a stochastic element occurring in these viral interactions, as observed in the distribution of diseased symptoms across a culture, differences in variance between experiments, and a number of additional viruses undergoing the step-jump in levels between experiments. Possible mechanisms for these multiple and simultaneous viral interactions in single culture are discussed in relation to known virus-host regulatory mechanisms for viral replication and whether additional factors could be considered to account for the 1,000-fold increase in AbV16 and AbV6 RNA2 levels.IMPORTANCE How viruses interact in a multiple-virus infection was examined by quantifying the levels of 18 viral RNAs in fruiting cultures of the agriculturally cultivated fungus Agaricus bisporus and statistically analyzing and modeling their abundance. Synergistic, antagonistic, and neutral interactions occurred simultaneously in cultures. The viral RNAs were grouped into four clusters, each displaying similar relative abundance, and between clusters, further interactions were found with positive, negative, or no correlations. Mushroom fruitbodies showing disease symptoms were distributed apparently randomly across the culture. These symptoms were associated with the presence of viral RNAs from two different clusters at very high levels, 1,000-fold higher than asymptomatic fruitbodies. The role of viral interaction together with stochastic factors and the regulation of host antiviral defenses in pathogenesis are discussed.

Multiple viral infections in Agaricus bisporus - Characterisation of 18 unique RNA viruses and 8 ORFans identified by deep sequencing.

Thirty unique non-host RNAs were sequenced in the cultivated fungus, Agaricus bisporus, comprising 18 viruses each encoding an RdRp domain with an additional 8 ORFans (non-host RNAs with no similarity to known sequences). Two viruses were multipartite with component RNAs showing correlative abundances and common 3' motifs. The viruses, all positive sense single-stranded, were classified into diverse orders/families. Multiple infections of Agaricus may represent a diverse, dynamic and interactive viral ecosystem with sequence variability ranging over 2 orders of magnitude and evidence of recombination, horizontal gene transfer and variable fragment numbers. Large numbers of viral RNAs were detected in multiple Agaricus samples; up to 24 in samples symptomatic for disease and 8-17 in asymptomatic samples, suggesting adaptive strategies for co-existence. The viral composition of growing cultures was dynamic, with evidence of gains and losses depending on the environment and included new hypothetical viruses when compared with the current transcriptome and EST databases. As the non-cellular transmission of mycoviruses is rare, the founding infections may be ancient, preserved in wild Agaricus populations, which act as reservoirs for subsequent cell-to-cell infection when host populations are expanded massively through fungiculture.

Functional analysis of Agaricus bisporus serine proteinase 1 reveals roles in utilization of humic rich substrates and adaptation to the leaf-litter ecological niche.

Agaricus bisporus is a secondary decomposer fungus and an excellent model for the adaptation, persistence and growth of fungi in humic-rich environments such as soils of temperate woodland and pastures. The A. bisporus serine proteinase SPR1 is induced by humic acids and is highly expressed during growth on compost. Three Spr1 gene silencing cassettes were constructed around sense, antisense and non-translatable-stop strategies (pGRsensehph, pGRantihph and pGRstophph). Transformation of A. bisporus with these cassettes generated cultures showing a reduction in extracellular proteinase activity as demonstrated by the reduction, or abolition, of a clearing zone on plate-based bioassays. These lines were then assessed by detailed enzyme assay, RT-qPCR and fruiting. Serine proteinase activity in liquid cultures was reduced in 83% of transformants. RT-qPCR showed reduced Spr1 mRNA levels in all transformants analysed, and these correlated with reduced enzyme activity. When fruiting was induced, highly-silenced transformant AS5 failed to colonize the compost, whilst for those that did colonize the compost, 60% gave a reduction in mushroom yield. Transcriptional, biochemical and developmental observations, demonstrate that SPR1 has an important role in nutrient acquisition in compost and that SPR1 is a key enzyme in the adaptation of Agaricus to the humic-rich ecological niche formed during biomass degradation.

Viral Agents Causing Brown Cap Mushroom Disease of Agaricus bisporus.

The symptoms of viral infections of fungi range from cryptic to severe, but there is little knowledge of the factors involved in this transition of fungal/viral interactions. Brown cap mushroom disease of the cultivated Agaricus bisporus is economically important and represents a model system to describe this transition. Differentially expressed transcript fragments between mushrooms showing the symptoms of brown cap mushroom disease and control white noninfected mushrooms have been identified and sequenced. Ten of these RNA fragments have been found to be upregulated over 1,000-fold between diseased and nondiseased tissue but are absent from the Agaricus bisporus genome sequence and hybridize to double-stranded RNAs extracted from diseased tissue. We hypothesize that these transcript fragments are viral and represent components of the disease-causing agent, a bipartite virus with similarities to the family Partitiviridae. The virus fragments were found at two distinct levels within infected mushrooms, at raised levels in infected, nonsymptomatic, white mushrooms and at much greater levels (3,500 to 87,000 times greater) in infected mushrooms exhibiting brown coloration. In addition, differential screening revealed 9 upregulated and 32 downregulated host Agaricus bisporus transcripts. Chromametric analysis was able to distinguish color differences between noninfected white mushrooms and white infected mushrooms at an early stage of mushroom growth. This method may be the basis for an "on-farm" disease detection assay.

Carbohydrate utilization and metabolism is highly differentiated in Agaricus bisporus.

BACKGROUND: Agaricus bisporus is commercially grown on compost, in which the available carbon sources consist mainly of plant-derived polysaccharides that are built out of various different constituent monosaccharides. The major constituent monosaccharides of these polysaccharides are glucose, xylose, and arabinose, while smaller amounts of galactose, glucuronic acid, rhamnose and mannose are also present. RESULTS: In this study, genes encoding putative enzymes from carbon metabolism were identified and their expression was studied in different growth stages of A. bisporus. We correlated the expression of genes encoding plant and fungal polysaccharide modifying enzymes identified in the A. bisporus genome to the soluble carbohydrates and the composition of mycelium grown compost, casing layer and fruiting bodies. CONCLUSIONS: The compost grown vegetative mycelium of A. bisporus consumes a wide variety of monosaccharides. However, in fruiting bodies only hexose catabolism occurs, and no accumulation of other sugars was observed. This suggests that only hexoses or their conversion products are transported from the vegetative mycelium to the fruiting body, while the other sugars likely provide energy for growth and maintenance of the vegetative mycelium. Clear correlations were found between expression of the genes and composition of carbohydrates. Genes encoding plant cell wall polysaccharide degrading enzymes were mainly expressed in compost-grown mycelium, and largely absent in fruiting bodies. In contrast, genes encoding fungal cell wall polysaccharide modifying enzymes were expressed in both fruiting bodies and vegetative mycelium, but different gene sets were expressed in these samples.

Agaricus bisporus genome sequence: a commentary.

The genomes of two isolates of Agaricus bisporus have been sequenced recently. This soil-inhabiting fungus has a wide geographical distribution in nature and it is also cultivated in an industrialized indoor process ($4.7bn annual worldwide value) to produce edible mushrooms. Previously this lignocellulosic fungus has resisted precise econutritional classification, i.e. into white- or brown-rot decomposers. The generation of the genome sequence and transcriptomic analyses has revealed a new classification, 'humicolous', for species adapted to grow in humic-rich, partially decomposed leaf material. The Agaricus biporus genomes contain a collection of polysaccharide and lignin-degrading genes and more interestingly an expanded number of genes (relative to other lignocellulosic fungi) that enhance degradation of lignin derivatives, i.e. heme-thiolate peroxidases and ß-etherases. A motif that is hypothesized to be a promoter element in the humicolous adaptation suite is present in a large number of genes specifically up-regulated when the mycelium is grown on humic-rich substrate. The genome sequence of A. bisporus offers a platform to explore fungal biology in carbon-rich soil environments and terrestrial cycling of carbon, nitrogen, phosphorus and potassium.

Environmental regulation of reproductive phase change in Agaricus bisporus by 1-octen-3-ol, temperature and CO2.

Reproductive phase change from vegetative mycelium to the initiation of fruiting in Agaricus bisporus is regulated in large part by the sensing of environmental conditions. A model is proposed in which three separate environmental factors exert control at different stages of the reproductive developmental process change. The eight carbon volatile 1-octen-3-ol controls the early differentiation from vegetative hyphae to multicellular knots; temperature reduction is essential for the later differentiation of primodia; and carbon dioxide level exerts quantitative control on the number of fruiting bodies developed. Analysis of transcriptomic changes during the reproductive phase change was carried out with initiation-specific microarrays, and the newly published A. bisporus genome was used to analyse the promoter regions of differentially regulated genes. Our studies have shown there to be both early and late initiation responses relating to sensing of eight carbon volatiles and temperature respectively. A subset of 45 genes was transcriptionally regulated during the reproductive phase change which exhibited a range of functions including cell structure, nitrogen and carbon metabolism, and sensing and signalling. Three gene clusters linking increased transcription with developmental stage were identified. Analysis of promoter regions revealed cluster-specific conserved motifs indicative of co-ordinated regulation of transcription.

Genome sequence of the button mushroom Agaricus bisporus reveals mechanisms governing adaptation to a humic-rich ecological niche.

Agaricus bisporus is the model fungus for the adaptation, persistence, and growth in the humic-rich leaf-litter environment. Aside from its ecological role, A. bisporus has been an important component of the human diet for over 200 y and worldwide cultivation of the "button mushroom" forms a multibillion dollar industry. We present two A. bisporus genomes, their gene repertoires and transcript profiles on compost and during mushroom formation. The genomes encode a full repertoire of polysaccharide-degrading enzymes similar to that of wood-decayers. Comparative transcriptomics of mycelium grown on defined medium, casing-soil, and compost revealed genes encoding enzymes involved in xylan, cellulose, pectin, and protein degradation are more highly expressed in compost. The striking expansion of heme-thiolate peroxidases and ß-etherases is distinctive from Agaricomycotina wood-decayers and suggests a broad attack on decaying lignin and related metabolites found in humic acid-rich environment. Similarly, up-regulation of these genes together with a lignolytic manganese peroxidase, multiple copper radical oxidases, and cytochrome P450s is consistent with challenges posed by complex humic-rich substrates. The gene repertoire and expression of hydrolytic enzymes in A. bisporus is substantially different from the taxonomically related ectomycorrhizal symbiont Laccaria bicolor. A common promoter motif was also identified in genes very highly expressed in humic-rich substrates. These observations reveal genetic and enzymatic mechanisms governing adaptation to the humic-rich ecological niche formed during plant degradation, further defining the critical role such fungi contribute to soil structure and carbon sequestration in terrestrial ecosystems. Genome sequence will expedite mushroom breeding for improved agronomic characteristics.

Genome organization and transcription response to harvest of two metallothionein-like genes in Agaricus bisporus fruiting bodies.

Metallothioneins are a class of small cysteine-rich proteins that have been associated with increased tolerance to metal and oxidative stresses in animals, plants, and fungi. We investigated a metallothionein-like (mt-like) gene shown previously to be upregulated in fruiting bodies of the fungus Agaricus bisporus in response to post-harvest storage. Analysis of an A. bisporus genomic DNA cosmid library identified two similar mt-like genes (met1 and met2) arranged as a bidirectional gene pair transcribed from the same promoter region. The promoter contained regulatory elements including 9 metal responsive elements and a CAAT box region 220 bp downstream of met1 that showed striking similarity to a feature in Coprinopsis cinerea mt-like gene promoters. Transcriptional analysis showed that both met genes are significantly and rapidly (within 3 hours) upregulated during post-harvest storage and expression is significantly greater in stipe and cap tissues compared with the gills. However, a strong directionality of the promoter was demonstrated, as transcript levels of met1 were at least two orders of magnitude greater than those of met2 in all samples tested.

Development of novel assays for lignin degradation: comparative analysis of bacterial and fungal lignin degraders.

Two spectrophotometric assays have been developed to monitor breakdown of the lignin component of plant lignocellulose: a continuous fluorescent assay involving fluorescently modified lignin, and a UV-vis assay involving chemically nitrated lignin. These assays have been used to analyse lignin degradation activity in bacterial and fungal lignin degraders, and to identify additional soil bacteria that show activity for lignin degradation. Two soil bacteria known to act as aromatic degraders, Pseudomonas putida and Rhodococcus sp. RHA1, consistently showed activity in these assays, and these strains were shown in a small scale experiment to breakdown lignocellulose, producing a number of monocyclic phenolic products. Using milled wood lignin prepared from wheat straw, pine, and miscanthus, some bacterial lignin degraders were found to show specificity for lignin type. These assays could be used to identify novel lignin degraders for breakdown of plant lignocellulose.

Influence of sporophore development, damage, storage, and tissue specificity on the enzymic formation of volatiles in mushrooms (Agaricus bisporus).

The enzymic oxidation of the polyunsaturated fatty acid-linoleic acid leads, in fungi, to the formation of a unique class of nonconjugated hydroperoxides, which are cleaved to form eight-carbon volatiles characteristic of mushroom and fungal flavor. However, the enzymes involved in this biosynthetic pathway, the bioavailability of the fatty acid substrate, and the occurrence of the reaction products (hydroperoxides and eight-carbon volatiles) are not fully understood. This study investigated the lipids, fatty acids, and hydroperoxide levels, as well as eight-carbon volatile variations in the fungal model Agaricus bisporus, according to four parameters: sporophore development, postharvest storage, tissue type, and damage. Eight-carbon volatiles were measured using solid phase microextraction and gas chromatography-mass spectrometry. Tissue disruption had a major impact on the volatile profile, both qualitatively and quantitatively; 3-octanone was identified as the main eight-carbon volatile in whole and sliced sporophore, an observation overlooked in previous studies due to the use of tissue disruption and solvent extraction for analysis. Fatty acid oxidation and eight-carbon volatile emissions decreased with sporophore development and storage, and differed according to tissue type. The release of 1-octen-3-ol and 3-octanone by incubation of sporophore tissue homogenate with free linoleic acid was inhibited by acetylsalicylic acid, providing evidence for the involvement of a heme-dioxygenase in eight-carbon volatile production.

Characterization of serine proteinase expression in Agaricus bisporus and Coprinopsis cinerea by using green fluorescent protein and the A. bisporus SPR1 promoter.

The Agaricus bisporus serine proteinase 1 (SPR1) appears to be significant in both mycelial nutrition and senescence of the fruiting body. We report on the construction of an SPR promoter::green fluorescent protein (GFP) fusion cassette, pGreen_hph1_SPR_GFP, for the investigation of temporal and developmental expression of SPR1 in homobasidiomycetes and to determine how expression is linked to physiological and environmental stimuli. Monitoring of A. bisporus pGreen_hph1_SPR_GFP transformants on media rich in ammonia or containing different nitrogen sources demonstrated that SPR1 is produced in response to available nitrogen. In A. bisporus fruiting bodies, GFP activity was localized to the stipe of postharvest senescing sporophores. pGreen_hph1_SPR_GFP was also transformed into the model basidiomycete Coprinopsis cinerea. Endogenous C. cinerea proteinase activity was profiled during liquid culture and fruiting body development. Maximum activity was observed in the mature cap, while activity dropped during autolysis. Analysis of the C. cinerea genome revealed seven genes showing significant homology to the A. bisporus SPR1 and SPR2 genes. These genes contain the aspartic acid, histidine, and serine residues common to serine proteinases. Analysis of the promoter regions revealed at least one CreA and several AreA regulatory motifs in all sequences. Fruiting was induced in C. cinerea dikaryons, and fluorescence was determined in different developmental stages. GFP expression was observed throughout the life cycle, demonstrating that serine proteinase can be active in all stages of C. cinerea fruiting body development. Serine proteinase expression (GFP fluorescence) was most concentrated during development of young tissue, which may be indicative of high protein turnover during cell differentiation.

Statistical modelling of transcript profiles of differentially regulated genes.

BACKGROUND: The vast quantities of gene expression profiling data produced in microarray studies, and the more precise quantitative PCR, are often not statistically analysed to their full potential. Previous studies have summarised gene expression profiles using simple descriptive statistics, basic analysis of variance (ANOVA) and the clustering of genes based on simple models fitted to their expression profiles over time. We report the novel application of statistical non-linear regression modelling techniques to describe the shapes of expression profiles for the fungus Agaricus bisporus, quantified by PCR, and for E. coli and Rattus norvegicus, using microarray technology. The use of parametric non-linear regression models provides a more precise description of expression profiles, reducing the "noise" of the raw data to produce a clear "signal" given by the fitted curve, and describing each profile with a small number of biologically interpretable parameters. This approach then allows the direct comparison and clustering of the shapes of response patterns between genes and potentially enables a greater exploration and interpretation of the biological processes driving gene expression. RESULTS: Quantitative reverse transcriptase PCR-derived time-course data of genes were modelled. "Split-line" or "broken-stick" regression identified the initial time of gene up-regulation, enabling the classification of genes into those with primary and secondary responses. Five-day profiles were modelled using the biologically-oriented, critical exponential curve, y(t) = A + (B + Ct)Rt + epsilon. This non-linear regression approach allowed the expression patterns for different genes to be compared in terms of curve shape, time of maximal transcript level and the decline and asymptotic response levels. Three distinct regulatory patterns were identified for the five genes studied. Applying the regression modelling approach to microarray-derived time course data allowed 11% of the Escherichia coli features to be fitted by an exponential function, and 25% of the Rattus norvegicus features could be described by the critical exponential model, all with statistical significance of p < 0.05. CONCLUSION: The statistical non-linear regression approaches presented in this study provide detailed biologically oriented descriptions of individual gene expression profiles, using biologically variable data to generate a set of defining parameters. These approaches have application to the modelling and greater interpretation of profiles obtained across a wide range of platforms, such as microarrays. Through careful choice of appropriate model forms, such statistical regression approaches allow an improved comparison of gene expression profiles, and may provide an approach for the greater understanding of common regulatory mechanisms between genes.

Hairpin-mediated down-regulation of the urea cycle enzyme argininosuccinate lyase in Agaricus bisporus.

A double-stranded (ds) RNA hairpin-mediated down-regulation system was developed for the cultivated mushroom Agaricus bisporus, and the role of the urea cycle enzyme argininosuccinate lyase (asl) in mushroom post-harvest development was investigated. Hairpin expression vectors were constructed to initiate down-regulation of asl and introduced into A. bisporus by Agrobacterium tumefaciens-mediated transformation. Transcripts of asl were significantly reduced (93.1 and 99.9%) in two transformants and hairpin vector transgene sequences were maintained throughout sporophore development. Single and multiple hairpin integration events were observed in Southern analysis. Transformants with down-regulated asl exhibited reduced yield and cap expansion during post-harvest sporophore development. There were no detectable differences in urea levels between the hairpin-transformed and control strains. This is the first report of reduced gene expression resulting from the introduction of dsRNA hairpins in A. bisporus and the applications of this technology will facilitate functional studies in the mushroom.

The role of ornithine aminotransferase in fruiting body formation of the mushroom Agaricus bisporus.

The complete oat gene and cDNA from the commercial mushroom, Agaricus bisporus, encoding ornithine aminotransferase (OAT) was characterized. The gene encodes a 466 amino acid protein and provides the first fully reported homobasidiomycete OAT protein sequence. The gene is interrupted by ten introns, and no mitochondrial targeting motif was present pointing to a cytoplasmic localization. The function of the gene was demonstrated by complementation of a Saccharomyces cerevisiae mutant unable to utilize ornithine as a sole source of nitrogen with an A. bisporus oat cDNA construct. Northern analysis of the oat gene together with the pruA gene (encoding Delta(1)-pyrroline-5-carboxylate dehydrogenase) showed that transcripts of both genes were lower during the first stages of fruiting body development. The higher expression of the oat gene in later stages of development, suggests the importance of ornithine metabolism for the redistribution of metabolites in the developing mushroom. Hplc analysis of all amino acids revealed that ornithine levels increased during fruiting body development whereas proline levels fell.

Argininosuccinate synthetase and argininosuccinate lyase: two ornithine cycle enzymes from Agaricus bisporus.

Accumulation of high quantities of urea in fruiting bodies is a known feature of larger basidiomycetes. Argininosuccinate synthetase (ASS) and argininosuccinate lyase (ASL) are two ornithine cycle enzymes catalysing the last two steps in the arginine biosynthetic pathway. Arginine is the main precursor for urea formation. In this work the nucleotide sequences of the genes and corresponding cDNAs encoding argininosuccinate synthetase (ass) and argininosuccinate lyase (asl) from Agaricus bisporus were determined. Eight and six introns were present in the ass and asl gene, respectively. The location of four introns in the asl gene were conserved among vertebrate asl genes. Deduced amino acid sequences, representing the first homobasidiomycete ASS and ASL protein sequences, were analysed and compared with their counterparts in other organisms. The ass ORF encoded for a protein of 425 amino acids with a calculated molecular mass of 47266Da. An alignment with ASS proteins from other organisms revealed high similarity with fungal and mammalian ASS proteins, 61-63% and 51-55% identity, respectively. The asl open reading frame (ORF) encoded a protein of 464 amino acids with an calculated mass of 52337Da and similar to ASS shared the highest similarity with fungal ASL proteins, 59-60% identity. Northern analyses of ass and asl during fruiting body formation and post-harvest development revealed that expression was significantly up-regulated from developmental stage 3 on for all the tissues studied. The expression reached a maximum at the later stages of fruiting body growth, stages 6 and 7. Both ass and asl genes were up-regulated within 3h after harvest showing that the induction mechanism is very sensitive to the harvest event and emphasizes the importance of the arginine biosynthetic pathway/ornithine cycle in post-harvest physiology.

Expression of the urease gene of Agaricus bisporus: a tool for studying fruit body formation and post-harvest development.

Fruit body initials of Agaricus bisporus contain high levels of urea, which decrease in the following developmental stages until stage 4 (harvest) when urea levels increase again. At storage, the high urea content may affect the quality of the mushroom, i.e. by the formation of ammonia from urea through the action of urease (EC 3.5.1.5). Despite the abundance of urea in the edible mushroom A. bisporus, little is known about its physiological role. The urease gene of A. bisporus and its promoter region were identified and cloned. The coding part of the genomic DNA was interrupted by nine introns as confirmed by cDNA analysis. The first full homobasidiomycete urease protein sequence obtained comprised 838 amino acids (molecular mass 90,694 Da, pI 5.8). An alignment with fungal, plant and bacterial ureases revealed a high conservation. The expression of the urease gene, measured by Northern analyses, was studied both during normal development of fruit bodies and during post-harvest senescence. Expression in normal development was significantly up-regulated in developmental stages 5 and 6. During post-harvest senescence, the expression of urease was mainly observed in the stipe tissue; expression decreased on the first day and remained at a basal level through the remaining sampling period.

L-malic acid biosensor for field-based evaluation of apple, potato and tomato horticultural produce.

A screen-printed three-electrode amperometric biosensor incorporating malic enzyme for the measurement of L-malic acid in apple, potato and tomato horticultural samples has been developed. The working electrode contained 0.38 mU of immobilised enzyme and was fabricated using rhodinised carbon to facilitate NADPH oxidation at an operating potential of +300 mV vs. Ag/AgCl compared with > +600 mV for bare carbon. The linear range of the sensor was 0.028-0.7 mM L-malic acid with relative standard deviations of 3.3-13.3%. When testing with real apple, potato and tomato samples, the sensor accuracy was within 13.7% of a standard commercially available photometric test kit. The sensor approach is cheap, simple to perform and rapid (6 min), requiring only buffer-electrolyte and a small sample volume.