Skip to Content
Web-Counter says you are visitor number

2-D Wavelet Packet Analysis of Structural Self-Organization and Morphogenic Regulation in Filamentous Fungal Colonies

Cameron L. Jones

Centre for Applied Colloid and BioColloid Science, Swinburne University of Technology, School of Chemical Sciences, GPO Box 218, Hawthorn, Victoria. 3122, Australia
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

Fungal colonies are composed of filaments called hyphae which self-organise into an interconnected, branching mycelium during growth. These organisms display species-specific branching behaviour and secrete exo-enzymes in response to nutrient assimilation. This paper details a new image-analysis software routine called 2-D Wavelet Packet Analysis (2-D WPA) which quantifies the global fractal dimension, and provides a unique index of fungal branching morphology. The 2-D WPA is superior to alternative analytical methods since it estimates scaling in the horizontal, vertical and diagonal directions at multiple spatial resolutions by mapping pixel points onto wavelet coefficients. The logarithmic decay of wavelet packet coefficient energy versus position provides an unbiased estimate of fractal scaling. Full source code is provided. The fungus Pycnoporus cinnabarinus secretes the phenol-oxidase exo-enzyme, laccase which is of use industrially for the bleaching of recycled paper/Kraft pulps and in other areas which use peroxidase chemistry. The fractal dimension was enumerated on basal Newspaper Only media and on media supplemented with the paramorphogen, Cellobiose. This compound caused profuse hyphal branching and a reduction in radial expansion. These changes in morphology were accompanied by an exponential increase in the expression of laccase into the extracellular medium in liquid culture. Hypersecretion was stimulated approx. 7070% over the control at 5 days. This study shows that the fractal dimension is a useful indicator of future enzyme yield and demonstrates a strong correlation between spatial pattern phenomena and physiological function.

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. Examples of both fungal morphologies are shown here.

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 graph of the relative oxidase activity in Units/mL produced over a five day induction/growth cycle has been prepared.

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.

DayOxidase activity on Newspaper Only (N Only) U/mL.%Oxidase activity on N + 2% Cellobiose U/mL.%
10.1230 +- 0.01841001.0150 +- 0.5640825
20.0810 +- 0.00501001.3710 +- 0.76201693
30.0370 +- 0.00351001.8867 +- 1.14605100
40.0257 +- 0.00551002.0690 +- 0.77508051
50.0327 +- 0.01851002.3120 +- 0.47207070

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:

Biochemical TestNewspaper OnlyN + 2% Cellobiose
Protein concentration (µ g/mL)2.02582.2431
Specific enzyme activity (U/mg)16.14181030.7164
Oxidase activity (U/mL)0.03272.312

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. [Click here to see the gel]. This stain revealed a single band for both newspaper only and N+2% cellobiose substrate media. However, the size of the two bands were clearly different with size an indication of relative activity concentration compared to volume. Further staining of the gel with silver express staining followed by Bio-Rad silver staining also demonstrated a single main band for laccase for both treatments at the same position, but revealed a significantly larger, more intense band for the N+2% cellobiose induction. This supports the other evidence (total protein, oxidase activity profiles, specific enzyme activity) that cellobiose is a potent paramorphogen which not only influences morphology (increasing the fractal dimension of branching), but alters the mechanism of protein synthesis and excretion.

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.

Acknowledgements

This work was supported by an ARC Collaborative Industry grant with Visy Paper Recycling

Bibliography

  1. Jones, C.L., Lonergan, G.T. and Mainwaring, D.E. 1996a. Wavelet packet computation of the Hurst exponent. Journal of Physics A: Mathematical and General. 29: 2509-2527.
  2. Forrest, S.R. and Witten Jr., T.A. 1979. Long-range correlations in smoke-particle aggregates. Journal of Physics A: Mathematical and General. 12: L109-117.
  3. Donelly, D.P., Wilkins, M.F. and Boddy, L. 1995. An intergrated image analysis approach for determining biomass, radial extent and box-count fractal dimension of macroscopic mycelial systems. Binary. 7: 19-28.
  4. Jones, C. L., Mainwaring, D.E. and Lonergan, G.T. 1996b. Prediction of fungal oxidase exo-enzyme expression using the fractal dimension. The 10th. International Biotechnology Symposium. August 25-30 (Sydney - Australia) Accepted.
  5. Jones, C.L., Lonergan, G.T. and Mainwaring, D.E. 1993. A rapid method for the fractal analysis of fungal colony growth using image processing. Binary. 5: 171-180.
  6. Jones, C.L., Lonergan, G.T., and Mainwaring, D.E. 1994. Specimen preparation for high contrast image processing and fractal analysis of branching fungal colonies. Biotechnology Techniques. 8:313-318.
  7. Lonergan, G.T., Panow, A., Jones, C.L., Schliephake, K., Ni, C.J. and Mainwaring, D.E. 1995. Physiological and biodegradative behaviour of the white-rot fungus, Pycnoporus cinnabarinus in a packed-bed bioreactor. Australasian Biotechnology. 5(2): 107-111.
  8. Hatakka, A.I. 1985. Degradation of veratric acid and other lignin-related aromatic compounds by the white-rot fungus Pycnoporus cinnabarinus. Archives of Microbiology. 141: 22-28.
  9. Wilson, R.W. and Niederpruem, D.J. 1967a. Cellobiose as a paramorphogen in Schizophyllum commune. Canadian Journal of Microbiology. 13:1663-1670.
  10. Wilson, R.W. and Niederpruem, D.J. 1967b. Control of ß- glucosidases in Schizophyllum commune. Canadian Journal of Microbiology. 13:1009-1020.
  11. Wilson, R.W. 1970. Influence of cellobiose on the colonial morphology of glycosidases of several fungi. Canadian Journal of Microbiology. 16:629-634.
  12. Wessels, J.G.H. 1965. Morphogenesis and Biochemical processes in Schizophyllum commune Fr. Wentia 13: 1-113.
  13. Blanchette, R.A. and Barnes, T.A. 1988. Selection of white-rot fungi for biopulping. Biomass. 15: 93-101.
  14. Jurasek, L. and Paice, M.G. 1988. Biological treatments of pulps. Biomass. 15: 103-108.
  15. Kaneko, R., Iimori, T., Miyawaki, S., Muchida, M. and Murakami, K. 1995. Biobleaching with manganese peroxidase purified from the pulp bleaching fungus SKB-1152. Bioscience, Biotechnology and Biochemistry. 59(8): 1584-1585.
  16. Jones, C.L. 1996c (May 28). Preprint of wavelet packet analysis fractal-software operating instructions. Centre for Applied Colloid and BioColloid Science, Swinburne University of Technology. URL: http://www.swin.edu.au/chem/bio/splus/splusfaq.htm
  17. Niku-Paavola, M.L. Raaska, L. and Itavaaru, M. 1990. Detection of white-rot fungi by a non-toxic stain. Mycological Research. 94: 27-31.
  18. Bourbonnais, R. and Paice, M.G. 1990. Oxidation of non phenolic substrates: An expanded role for laccase in lignin biodegradation. FEBS Letters. 267: 99-102.
  19. Jones, C.L., Lonergan, G.T. and Mainwaring, D.E. 1995a. Acid phosphatase positional correlations in solid surface fungal cultivation: a fractal interpretation of biochemical differentiation. Biochemical and Biophysical Research Communications. 208(3): 1159-1165.
  20. Jones, C.L., Lonergan, G.T. and Mainwaring, D.E. 1995b. Sierpinski fractal analysis and percolation threshold implications in fungal colonies. Bioimages. 3(2): 71-84.

Page counts provided by Web-Counter: