2-D Wavelet Packet Analysis of Structural Self-Organization and Morphogenic Regulation in Filamentous Fungal Colonies
Cameron L. Jones
Email: CJONES@swin.edu.au
This preprint should be referenced as follows:
TITLE: "2-D Wavelet Packet Analysis of Structural Self-Organisation and Morphogenic Regulation in Filamentous Fungal Colonies"AUTHOR: Cameron L. Jones
SUBMITTED TO: Complex Systems Conference - From Local Interactions to Global Phenomena; July 14-17, 1996 Charles Sturt University, Albury, Wodonga, Australia.
DATE: 31 May 1996
URL REFERENCE: http://www.swin.edu.au/chem/bio/cs96/camjones.htm
Abstract
Introduction
Filamentous fungal colonies are eukaryotic microorganisms which grow by repeated branching of the hyphal tips to form a ramified mycelial network. Understanding how fungi grow is important since many strains are opportunistic medical pathogens, while others are fundamental for use in industrial fermentations and for the manufacture of antibiotics and bio-pharmaceuticals. Each fungal species displays a unique hyphal morphology and secretes exo-enzymes following nutrient assimilation. An objective of this work has been to investigate the link between colony structure and enzymatic function, with the aim of providing new methods to increase the rate of protein expression for industrial harvest. The branching patterns observed in fungal colonies are similar in appearence to Diffusion-Limited Aggregates (DLA), and new image-analysis software has been developed to quantify the global Fractal Dimension. By growing fungal colonies on different nutrient media formulations, we hoped to correlate branching morphology (indexed with the fractal dimension) with enzyme expression yield. This paper reports the results obtained using the new 2-D Wavelet Packet routine, which is an extension of previous work with the 1-D Wavelet Packet Transform (Jones et al., 1996a). Wavelet analysis is similar to, but an extension of the Fourier Transform (FT). The FT is useful for the analysis of linear functions but is less suited for nonlinear signal fluctuations. Wavelets transform a signal into a series of coefficients having discrete energies. The original signal is completely specified by these coefficients and the analyzing wavelet. In 2-D, wavelets offer a method to quantify how energy is spatially distributed at multiple scales and also how it is partitioned in the horizontal, vertical and diagonal directions at different resolutions. This method is superior to alternative imaging procedures principally due to the fact that multiple analysis of the same image is not required (c.f. the Sandbox method - Forest and Witten, 1979) and it does not suffer from edge-bias (c.f. dilation/erosion Box Counting - Donelly et al., 1995). Furthermore, all statistics are returned in a single window. 2-D WPA showed good accuracy when it was applied to 'test objects' of known dimension and returned a standard deviation of ±0.0239. Complete statistics regarding accuracy are available in a preprint referenced below. We also report that the addition of the paramorphogen, cellobiose results in an increased branching response (a higher fractal dimension) with a concomitant increase in the expression of the phenol-oxidase enzyme, laccase. Possible physiological reasons why cellobiose alters the self-organisational dynamics of fungal colony branching is also explored.
Modification of Protein Expression in Fungi:
This paper focuses on a new method to experimentally influence fungal protein expression mechanisms, with the aim to overproduce or hypersecrete extracellular oxidase enzymes. A paramorphogen is thought to stimulate or induce changes in fungal cell physiology by altering either the:- Cell wall biochemistry
- cytosol molar ratios of transport proteins which are involved with protein synthesis or,
- the cell-mycelium morphology by enhancing or repressing hyphal branching.
Our interest in paramorphogen stimulation developed from an earlier study which established a correlation between hyphal branching and the number density of active tip production per unit area and extracellular oxidase expression levels (Jones et al., 1996b). Image-analysis of mycelial morphology was used to quantitatively measure the macroscopic branching index of isolated whole fungal colonies growing on solid surface culture (SSC). SSC can be on a large industrial scale or 'bench scale' using agar filled Petri-plates. The Fractal Dimension, D was shown to be a sensitive numerical estimate of branching behaviour (Jones et al., 1993, 1994). This index estimates the amount of surface coverage by the fungal mycelium. For two dimensional (2-D) objects in the plane, 1<D<2. We consider fungal colonies which grow across the agar surface or on a membrane covering the agar, to be two dimensional since we view them as a projection against a flat 2-D plane using transmitted illumination. Straight lines have a D value=1 and a plane filling curve has a D value=2. Values of D within this range therefore characterize the extent of space filling. Highly branched fungal colonies have a 'high' fractal dimension compared against sparsely branched colonies. All measurements must be carried out at a standardized time, and all culture conditions must be well defined. Obviously, measuring various D values for discrete colonies growing on SSC is interesting, but of greater import is the possible relationship between growth in 2- D surface culture and growth in unrestricted 3-D liquid culture. Parallel experiments in liquid culture have established a positive correlation between the Fractal dimension of branching and extracellular oxidase secretion into the medium (Jones et al., 1996b). Over time, the fungal mycelium secretes exo-enzymes and the kinetic 'efficiency' of their production is of considerable industrial importance (Lonergan et al., 1995; Hatakka, 1985) . The activity of the crude enzyme (i.e. laccase in Pycnoporus cinnabarinus) for performing useful substrate transformation (oxidation) was determined by measuring the decolorization (degradation) of the aminoanthraquinone dye, Remazol Brilliant Blue R. The crude extracellular fluid was analyzed daily and shown to contain (i) oxidase enzyme activity and that (ii) oxidase enzyme could decolorize the dye. Of most importance was the fact that the fractal dimension measured in very young colonies (24hrs.) was positively correlated with future enzyme secretion potential.
Global Objective:
The aim was to identify compounds which modify the branching response (in a similar way to the dye example described above) and which could be used industrially to control protein synthesis and expression of fungal exo-enzymes.The chosen paramorphogen - cellobiose is a carbohydrate-linked glucose disaccharide. The literature documents that cellobiose generally causes a restriction in fungal growth rate and may induce profuse branching, compared against untreated or wild-type growth. Its mode of action is thought to act upon cell wall metabolism, and we review this evidence below and then report on our results with cellobiose induction in P. cinnabarinus.
Several studies have confirmed that cellobiose can induce both macroscopic changes to structural morphology, and microscopic changes at the level of intracellular metabolism (Wilson & Niedepruem,(1967a),1967b; Wilson, 1970). A substance called R-glucan (present in the cell wall) is thought to provide a structural shape preserving action to the fungal cell. Since cellobiose is known to effect cell wall metabolism, and the R-glucan fraction maintains hyphal cell shape (Wessells, 1965) a practical experiment carried out by Wilson (1970) sought to evaluate the changes in R-glucanase activity in response to cellobiose induction. Measurement of R-glucanase (the enzyme against R-glucan) therefore provides an index of potential hydrolytic action against this particular cell wall fraction. It was confirmed (Wilson, 1970) that when several basidiomycete fungi were grown on cellobiose containing media, the levels of several wall-hydrolysing enzymes (including R-glucanase) were all increased. It was concluded that elevated levels of cell wall hydrolysing enzymes provide multiple regions of wall softening, thereby resulting in increased branching per unit area. It has also been noted that hydrolytic enzymes can participate in transfer reactions which suggests the possibility that cell wall synthesis may also be a factor in morphogenesis (Wilson, 1970). They concluded that the elevation in R-glucanase activity could physiologically increase cell wall plasticity and result in an increased branching response.
Materials and Methods
Discrete young colonies (24hrs) of P. cinnabarinus were examined microscopically and characterized with image analysis to determine the global fractal dimension of branching. Similar evaluations were performed on mycelia germinated on media lacking cellobiose. The nutrient media consisted of non-defined powdered newspaper (200 mesh pass). 'The Age' newspaper was used for preparation of all media and contained no color pages. It was decided to use newspaper as a model nutrient substrate since one of the future applications of hypersecreting fungi (and their exo-enzymes) is for the bio-treatment of recycled pulps (Blanchette and Burnes, 1988; Jurasek and Paice, 1988; Kaneko et al., 1995). An additional benefit of this research was the fact that old newspaper in Australia is currently collected, but very little is recycled (approximately 2%) while the remainder is landfilled. This made newspaper an attractive waste substrate for a nutrient medium since it was readily available, and on a larger commercial scale would be reasonably cheap to acquire.Powdered newspaper at 0.2% w/v was solidified with 2% w/v Technical Grade Agar No. 3 (Oxoid). The sole carbon source was therefore newspaper only (N) for the control, while the cellobiose treated plates were prepared by adding 2% w/v cellobiose (Sigma). Media was poured into 9cm diameter Petri plates at 25-30mL final volume. Colonies of P. cinnabarinus were inoculated following the previously described protocol for the germination of discrete, single colonies on microporous polycarbonate membrane overlays (Jones et al.,1994). The colonies were also contrast stained with Coomassie Brilliant Blue and examined using an Olympus CH- 2 optical microscope set up for bright field transmission. The methods used to apply 2-D WPA are essentially the same as those used for 1-D WPA discussed in detail in Jones et al. (1996a). However we have used the 2-D symmetrical (s.12) wavelet filter with 12 vanishing methods, using input images of size: 512x512. More complete details will shortly be made available when the software is lodged with the StatLib software repository for S-Plus written software (Carnegie Mellon University Statistics Department). A brief preprint detailing step-by-step instructions for operating the software is available here (Jones, 1996c). Source code (in Microsoft Word ver 6.x for Windows format) is available individually for the following 2-D image sizes: 128x128, 192x192, 256x256, 384x384, 300x200, 512x512.
To examine the effect of cellobiose supplementation versus newspaper only controls, the following experiment was also conducted. Liquid media were prepared using the same newspaper and cellobiose concentrations as for solid media in Petri plates, however technical agar was not added. Static batch flasks (250mL volume) containing 25mL of liquid media were inoculated with 2x5mm diam. plugs of 5 day old stock cultures of P. cinnabarinus taken from the exponentially growing frontier. All flasks were incubated under humidified conditions at 37°C for five days, by incubating the flasks in trays covered with plastic wrap. Aliquots (700µl volume) were removed daily under aseptic conditions for spectrophotometric analysis of oxidase activity. Phenol oxidase activity of laccase was measured at 415nm over a period of 3 minutes at room temperature. Oxidase enzyme activity was measured as activity against ABTS (Niku-Paavola et al., 1990; Bourbonnais and Paice, 1992). This spectrophotometric assay used a reaction mixture containing 0.5mL of 100mM sodium succinate buffer at pH 4.5; 0.1mL of 2.92mM ABTS; 0.3mL of distilled water and 0.1mL of extracellular culture fluid. One unit (U) of enzyme was defined as the amount which caused a change in absorbance (415nm) of 1.0 per minute at room temperature. All activities were reported as U mL-1 (Table 1).
Although we have monitored only oxidase enzyme, the fate of other extracellular proteins have been established in two ways: Total protein was measured on Day 5 run-out samples and could be compared against oxidase activity to provide a specific activity index on each day. This led to further enzyme expression profiles following cellobiose induction (Table 2). Specifically, the Day 5 extracellular fluid was concentrated by Amicon Ultracentrifugation using Centriprep-10 and Centricon-10 filtration units (10,000 dalton M.W. cut-off) for three hours. This reduced the total volume from approximately 20mL to 40µL of concentrated sample. Proteins were separated by electrophoresis on native PAGE (12%) by loading 20µL of sample and 5µL of buffer. Activity staining using guaiacol established the band position of oxidase, and the gels could be further stained with Coomassie Brilliant Blue R or silverstain enhanced to reveal additional protein bands.
Extracellular protein was assayed at Day 5 for both media treatments using a Bio-Rad Microassay Kit for protein, using bovine serum albumin as standard. The specific enzyme activity (U/mg) was estimated by dividing the oxidase activity (U/mL) by the protein concentration (g/mL). This provides a semi-quantitative biochemical index of the amount of extracellular proteins which are oxidases.
Results and Discussion
Colonies of P. cinnabarinus grown on newspaper media supplemented with 2% cellobiose (N+Cell) displayed the following characteristics: (i) radial growth rate was significantly reduced at 24hrs compared against growth achieved on N only. This was obvious since an increased magnification of (x100) was required to visualise each discrete fungal colony at a comparable final size to compare against the N only control colonies which were recorded at x50 magnification. (ii) A significant increase in branching was observed throughout the colony both apically and subapically. (iii) Some colonies exhibited tip splitting and selected hyphae were wavy or 'gnarled' in appearance. (iv) The oldest regions of the colony exhibited density filling usually via straight extending hyphae with little or no branching. Control (N only) colonies were composed of longer, hyphae which displayed little sub-apical branching and were not wavy or 'knarled' in appearance.The mean fractal dimension measured with 2-D WPA for newspaper only is: D=1.406±0.124 and the mean value for N+Cell is: D=1.501±0.076. The global fractal dimension provides an unbiased estimate of the density scaling of hyphal branching and clearly confirms that the N+Cell media induced an increase in branching compared against the control. Although the colonies germinated on N+Cell appear significantly more branched than on the control, this was not reflected in the global fractal dimension (where the increase was approximately 0.1 scaling units). This was because more straight line linear branches were produced on N+Cell which did not bifurcate and there was little tip splitting. The fractal dimension estimates straight lines (or unbranched hyphae) as a D value = 1. A multifractal estimation algorithm would detect a larger difference between the two treatments compared against monofractal estimation. Software and the appropriate algorithms to do this are currently being developed.
The addition of cellobiose to the media clearly caused a
significant induction in extracellular oxidase enzyme
activity. A Table 1 shows the yield of oxidase produced on both
nutrient media formulations. Clearly cellobiose
supplementation induced the greatest amount of oxidase
synthesis and expression. The standard deviation between
the yields measured in triplicate flasks for each
treatment are provided as error bar estimates. There is
little variation in yield for the newspaper only, but for
the cellobiose supplemented media, the yield fluctuations
are significantly higher. Notably we see an increase in
the standard deviation between Day 1 and 3, with maximum
deviation on Day 3. Activity fluctuations then decrease
between Days 4 and 5 although the yield is still
increasing exponentially, while the newspaper only
control has entered the plateau growth phase. Table 1
provides a detailed listing of the oxidase activity
measured on each day and shows the percentage of oxidase
enzyme activity for cellobiose measured against newspaper
only - standardized as 100%. At Day 5, the yield of
oxidase on cellobiose was over 70 fold that seen on
newspaper only controls.
Table 1. Phenol oxidase (laccase) enzyme yield measured
against ABTS. The contol enzyme activity decreased over
time, and was benchmarked at 100% to compare against the
yield expression of laccase in cellobiose treated
cultures. Therefore for the cellobiose treated cultures,
the percentage shows the relative % increase against the
cotrol on that particular day.
By estimating the percentage increase in oxidase activity
versus the control, it was possible to quantify the
percentage induction caused by cellobiose (Table 1
above). The increase was exponential from day 1 to 5
with the following respective percentage trend: 825,
1693, 5100, 8051, 7070%. This is a very large increase
in enzyme activity and supports the hypothesis that
cellobiose acts as a paramorphogen to increase branching
which in turn promotes the expression and secretion of
exoenzyme - specifically the phenol-oxidase, laccase into
the extracellular medium.
These results detailing an increase in oxidase enzyme following
cellobiose induction is the first report of an
extracellular phenol-oxidase (laccase) being
exogeneously controlled by a paramorphogen. This
demonstrates that measurement of the fractal dimension in
2-D surface grown colonies to detect an increase in the
fractal dimension is a sensitive indicator of potential
secretion behavior of the same fungal microorganism grown
in liquid media. It should be stressed that the static
batch flask arrangement is sub-optimal for gas exchange
and nutrient diffusion due to non-stirring conditions.
It is anticipated that scale-up to more 'microbially optimal'
fermentation conditions should result in superior yield (work in progress).
We note that the result for induction of oxidase
following exposure to cellobiose is in agreement with the
finding of Wilson and Niederpreum (1967b) who established
that ß-glucosidase was induced by cellobiose. The present
results coupled with the review of the relevant
literature confirms the paramorphogenic enzyme induction
ability of cellobiose in several strains of basidiomycete
fungi, and has now been confirmed to also effect similar
changes in P. cinnabarinus. At Day 5, the protein concentration and specific enzyme activities are as follows:
Table 2. Comparison of total extracellular protein,
specific enzyme activity, and oxidase activity measured
in Day 5 run-out samples grown on the two treatment
medias. All values are the means of triplicate readings.
Table 2 confirms that although both treatments return
similar total protein concentrations, the very high
specific enzyme activity for the cellobiose media implies
that most of the induced protein was oxidase activity.
Polyacrylamide gel electrophoresis (PAGE) was also used to
provide a measure of the number and location of protein
bands which were expressed on each of the two media
formulations. We were especially interested in the
position and size of the laccase band. Activity
staining with guaicol was used to detect the phenol-
oxidase enzyme, laccase. [ These results imply that a 'non-linear dynamical' approach to biology can yield fascinating new results. First motivations were to develop an improved computer algorithm to quantify the global fractal dimension for digital images. Because wavelet analysis is especially sensitive to fluctuations in spatial pixel positions at multiple scales and resolutions, it proved extremely useful with comparable accuracy to established computer methods. Secondly, it was decided to apply this new imaging technique to measure morphology changes which occur when filamentous fungi grow on different nutrient formulations. At first, fungi were seen simply as 'nice' stochastic examples of self-similar fractals. However a positive correlation was found between the amount of branching in 2-D and the expression of laccase enzyme in 3-D liquid media. This implied that in this example, a simple morphology index was correlated with physiological function. This study would not have been realised with established experimental techniques. The theoretical framework of fractal geometry as well as the unifying mechanism of self-organised criticality underpinned the experimental approach taken here. In addition, these results are impressive (regardless of theory) from a practical industrial viewpoint. This paper is the result of a dual approach between practical application and fundamental, curiosity driven research. This study supports previous efforts to identify and measure fractal, percolation and critical phenomena in branching fungal colonies (Jones et al., 1995a,b). Morphogenic self-organisation in fungi has been shown to be an interconnected dynamic process which links spatial structure with biochemical function at the cellular and sub-cellular levels of scale.
Page counts provided by Web-Counter:
Day Oxidase activity on Newspaper Only (N Only) U/mL. % Oxidase activity on N + 2% Cellobiose U/mL. % 1 0.1230 +- 0.0184 100 1.0150 +- 0.5640 825 2 0.0810 +- 0.0050 100 1.3710 +- 0.7620 1693 3 0.0370 +- 0.0035 100 1.8867 +- 1.1460 5100 4 0.0257 +- 0.0055 100 2.0690 +- 0.7750 8051 5 0.0327 +- 0.0185 100 2.3120 +- 0.4720 7070
Biochemical Test Newspaper Only N + 2% Cellobiose Protein concentration (µ g/mL) 2.0258 2.2431 Specific enzyme activity (U/mg) 16.1418 1030.7164 Oxidase activity (U/mL) 0.0327 2.312 Acknowledgements
This work was supported by an ARC Collaborative Industry grant with Visy Paper RecyclingBibliography