Search

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 + ε. 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.

Contains fulltext : 36611.pdf (Publisher’s version ) (Closed access) 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.

Cloning and characterisation of a putative riboflavin-aldehyde-forming enzyme gene (raf) from the cultivated mushroom Agaricus bisporus and its expression during morphogenesis are described. Three cDNA clones were isolated following differential screening of cDNA libraries from rapidly expanding sporophores and post-harvest stored sporophores. The cDNA sequence and predicted translation analysis revealed an open reading frame (ORF) of 348 nucleotides encoding a polypeptide of 115 amino acids, with three introns (56-66 bases) interrupting the genomic ORF. Blast X searches of the databases with the gene sequence showed homology (40% identity and 56% similarity) to the riboflavin-aldehyde-forming enzyme gene from Schizophyllum commune. In A.bisporus, the raf gene sequence upstream of the ORF contained a large CT-rich putative regulatory element (-64 to -24 bases) found in highly expressed genes in various mushrooms, and a 6-base motif present in the 3' end of the genomic sequence, but not in the corresponding 3' non-coding part of the cDNA, was identified. The raf gene transcripts increased abundantly in rapidly developing sporophores as well in post-harvest stored sporophores. Differential expression of the raf gene transcripts in different tissues of the sporophore was also observed, with higher levels in the stipe compared with the cap and gills. The temporal and spatial expression patterns observed suggest transcriptional regulation of the raf gene during A. bisporus morphogenesis.

Fungi play a key role in ecosystem nutrient cycles by scavenging, concentrating, translocating and redistributing nitrogen. To quantify and predict fungal nitrogen redistribution, and assess the importance of the integrity of fungal networks in soil for ecosystem function, we need better understanding of the structures and processes involved. Until recently nitrogen translocation has been experimentally intractable owing to the lack of a suitable radioisotope tracer for nitrogen, and the impossibility of observing nitrogen translocation in real time under realistic conditions. We have developed an imaging method for recording the magnitude and direction of amino acid flow through the whole mycelial network as it captures, assimilates and channels its carbon and nitrogen resources, while growing in realistically heterogeneous soil microcosms. Computer analysis and modeling, based on these digitized video records, can reveal patterns in transport that suggest experimentally testable hypotheses. Experimental approaches that we are developing include genomics and stable isotope NMR to investigate where in the system nitrogen compounds are being acquired and stored, and where they are mobilized for transport or broken down. The results are elucidating the interplay between environment, metabolism, and the development and function of transport networks as mycelium forages in soil. The highly adapted and selected foraging networks of fungi may illuminate fundamental principles applicable to other supply networks. © Cambridge University Press.

Increases in both the levels and the activity of serine proteinase have been previously described in the senescing mushroom Agaricus bisporus. cDNA encoding serine proteinase was amplified by reverse transcriptase-polymerase chain reaction using a degenerate primer based on the N-terminal sequence of a previously isolated A. bisporus serine proteinase and then cloned. The cDNA was sequenced and shown to be homologous to those of other fungal serine proteinases. Northern analysis showed that this serine proteinase gene (Spr1) was not expressed in freshly harvested sporophores but was strongly up-regulated postharvest and found almost entirely in the stipe of the sporophore (approximately 0.08% of mRNAs 2 days after harvest). Low-level expression was detectable in the flesh (pileus trama) and gill (lamellae) tissues of the cap, but none was detected in the skin (pilei pellis). In three of the cloned cDNAs, sequence analysis showed that the poly(A) tail starts at different positions. Expression of Spr1 in Escherichia coli caused restricted colony growth.

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.