Growth and function of fungal mycelia in heterogeneous environments

Flow modeling and structural characterization in fungal pellets

Fungal pellets are hierarchical systems that can be found in an ample variety of applications. Modeling transport phenomena in this type of systems is a challenging but necessary task to provide knowledge-based processes that improve the outcome of their biotechnological applications. In this work, an upscaled model for total mass and momentum transport in fungal pellets is implemented and analyzed, using elements of the volume averaging and adjoint homogenization methods departing from the governing equations at the microscale in the intracellular and extracellular phases. The biomass is assumed to be composed of a non-Newtonian fluid and the organelles impervious to momentum transport are modeled as a rigid solid phase. The upscaled equations contain effective-medium coefficients, which are predicted from the solution of adjoint closure problems in a three-dimensional periodic domains representative of the microstructure. The construction of these domains was performed for Laccaria trichodermophora based on observations of actual biological structures. The upscaled model was validated with direct numerical simulations in homogeneous portions of the pellets core. It is shown that no significant differences are observed when the dolipores are open or closed to fluid flow. By comparing the predictions of the average velocity in the extracellular phase resulting from the upscaled model with those from the classical Darcy equation (i.e., assuming that the biomass is a solid phase) the contribution of the intracellular fluid phase was evidenced. This work sets the foundations for further studies dedicated to transport phenomena in this type of systems.

Imaging with spatio-temporal modelling to characterize the dynamics of plant-pathogen lesions

Within-host spread of pathogens is an important process for the study of plant-pathogen interactions. However, the development of plant-pathogen lesions remains practically difficult to characterize beyond the common traits such as lesion area. Here, we address this question by combining image-based phenotyping with mathematical modelling. We consider the spread of Peyronellaea pinodes on pea stipules that were monitored daily with visible imaging. We assume that pathogen propagation on host-tissues can be described by the Fisher-KPP model where lesion spread depends on both a logistic growth and an homogeneous diffusion. Model parameters are estimated using a variational data assimilation approach on sets of registered images. This modelling framework is used to compare the spread of an aggressive isolate on two pea cultivars with contrasted levels of partial resistance. We show that the expected slower spread on the most resistant cultivar is actually due to a significantly lower diffusion coefficient. This study shows that combining imaging with spatial mechanistic models can offer a mean to disentangle some processes involved in host-pathogen interactions and further development may allow a better identification of quantitative traits thereafter used in genetics and ecological studies.

Biophysical Manipulation of the Extracellular Environment by Eurotium halophilicum

Eurotium halophilicum is psychrotolerant, halophilic, and one of the most-extreme xerophiles in Earth's biosphere. We already know that this ascomycete grows close to 0 °C, at high NaCl, and-under some conditions-down to 0.651 water-activity. However, there is a paucity of information about how it achieves this extreme stress tolerance given the dynamic water regimes of the surface habitats on which it commonly occurs. Here, against the backdrop of global climate change, we investigated the biophysical interactions of E. halophilicum with its extracellular environment using samples taken from the surfaces of library books. The specific aims were to examine its morphology and extracellular environment (using scanning electron microscopy for visualisation and energy-dispersive X-ray spectrometry to identify chemical elements) and investigate interactions with water, ions, and minerals (including analyses of temperature and relative humidity conditions and determinations of salt deliquescence and water activity of extracellular brine). We observed crystals identified as eugsterite (Na4Ca(SO4)3·2H2O) and mirabilite (Na2SO4·10H2O) embedded within extracellular polymeric substances and provide evidence that E. halophilicum uses salt deliquescence to maintain conditions consistent with its water-activity window for growth. In addition, it utilizes a covering of hair-like microfilaments that likely absorb water and maintain a layer of humid air adjacent to the hyphae. We believe that, along with compatible solutes used for osmotic adjustment, these adaptations allow the fungus to maintain hydration in both space and time. We discuss these findings in relation to the conservation of books and other artifacts within the built environment, spoilage of foods and feeds, the ecology of E. halophilicum in natural habitats, and the current episode of climate change.

Mathematical modelling of fungi-initiated siderophore-iron interactions

Nearly all life forms require iron to survive and function. Microorganisms utilize a number of mechanisms to acquire iron including the production of siderophores, which are organic compounds that combine with ferric iron into forms that are easily absorbed by the microorganism. There has been significant experimental investigation into the role, distribution and function of siderophores in fungi but until now no predictive tools have been developed to qualify or quantify fungi-initiated siderophore-iron interactions. In this investigation, we construct the first mathematical models of siderophore function related to fungi. Initially, a set of partial differential equations are calibrated and integrated numerically to generate quantitative predictions on the spatio-temporal distributions of siderophores and related populations. This model is then reduced to a simpler set of equations that are solved algebraically giving rise to solutions that predict the distributions of siderophores and resultant compounds. These algebraic results require the calculation of zeros of cross products of Bessel functions and thus new algebraic expansions are derived for a variety of different cases that are in agreement with numerically computed values. The results of the modelling are consistent with experimental data while the analysis provides new quantitative predictions on the time scales involved between siderophore production and iron uptake along with how the total amount of iron acquired by the fungus depends on its environment. The implications to bio-technological applications are briefly discussed.

Phase-field modeling of constrained interactive fungal networks

Fungi develop structures that interact with their surroundings and evolve adaptively in the presence of geometrical constraints, finding optimal solutions for complex combinatorial problems. The pathogenic fungus Ophiocordyceps constitutes a perfect model for the study of constrained interactive networks. Modeling these networks is challenging due to the highly coupled physics involved and their interaction with moving boundaries. In this work, we develop a computational phase-field model to elucidate the mechanics of the emerging properties observed in fungal networks. We use a variational approach to derive the equations governing the evolution in time of the mycelium biomass and the nutrients in the medium. We present an extensive testing of our model, reproduce growing and decaying phenomena, and capture spatial and temporal scales. We explore the variables interplay mechanism that leads to different colony morphologies, and explain abrupt changes of patterns observed in the laboratory. We apply our model to simulate analogous processes to the evolution of Ophiocordyceps as it grows through confined geometry and depletes available resources, demonstrating the suitability of the formulation to study this class of biological networks.

Fungal bioremediation of soil co-contaminated with petroleum hydrocarbons and toxic metals

Abstract

Much research has been carried out on the bacterial bioremediation of soil contaminated with petroleum hydrocarbons and toxic metals but much less is known about the potential of fungi in sites that are co-contaminated with both classes of pollutants. This article documents the roles of fungi in soil polluted with both petroleum hydrocarbons and toxic metals as well as the mechanisms involved in the biotransformation of such substances. Soil characteristics (e.g., structural components, pH, and temperature) and intracellular or excreted extracellular enzymes and metabolites are crucial factors which affect the efficiency of combined pollutant transformations. At present, bioremediation of soil co-contaminated with petroleum hydrocarbons and toxic metals is mostly focused on the removal, detoxification, or degradation efficiency of single or composite pollutants of each type. Little research has been carried out on the metabolism of fungi in response to complex pollutant stress. To overcome current bottlenecks in understanding fungal bioremediation, the potential of new approaches, e.g., gradient diffusion film technology (DGT) and metabolomics, is also discussed.

Key points

• Fungi play important roles in soil co-contaminated with TPH and toxic metals.

• Soil characteristics, enzymes, and metabolites are major factors in bioremediation.

• DGT and metabolomics can be applied to overcome current bottlenecks.

Carbon metabolism and transcriptome in developmental paths differentiation of a homokaryotic Coprinopsis cinerea strain

The balance and interplay between sexual and asexual reproduction is one of the most intriguing mysteries in the study of fungi. The choice of developmental strategy reflects the ability of fungi to adapt to the changing environment. However, the evolution of developmental paths and the metabolic regulation during differentiation and morphogenesis are poorly understood. Here, an analysis was performed of carbohydrate metabolism and gene expression regulation during the early differentiation process from the vegetative mycelium, to the differentiated structures, fruiting body, oidia and sclerotia, of a homokaryotic fruiting Coprinopsis cinerea strain A43mutB43mut pab1-1 #326. Changes during morphogenesis and the evolution of developmental strategies were followed. Conversion between glucose and glycogen and between glucose and beta-glucan were the main carbon flows in the differentiation processes. Genes related to carbohydrate transport and metabolism were significantly differentially expressed among paths. Sclerotia displayed a set of specifically up-regulated genes that were enriched in the carbon metabolism and energy production and conversion processes. Evolutionary transcriptomic analysis of four developmental paths showed that all transcriptomes were under the purifying selection, and the more stressful the environment, the younger the transcriptome age. Oidiation has the lowest value of transcriptome age index (TAI) and transcriptome divergence index (TDI), while the fruiting process has the highest of both indexes. These findings provide new insights into the regulations of carbon metabolism and gene expressions during the early stages of fungal developmental paths differentiation, and improve our understanding of the evolutionary process of life history and reproductive strategy in fungi.

Effect of depleted uranium on a soil microcosm fungal community and influence of a plant-ectomycorrhizal association

Fungi are one of the most biogeochemically active components of the soil microbiome, becoming particularly important in metal polluted terrestrial environments. There is scant information on the mycobiota of uranium (U) polluted sites and the effect of metallic depleted uranium (DU) stress on fungal communities in soil has not been reported. The present study aimed to establish the effect of DU contamination on a fungal community in soil using a culture-independent approach, fungal ribosomal intergenic spacer analysis (F-RISA). Experimental soil microcosms also included variants with plants (Pinus silvestris) and P. silvestris/Rhizopogon rubescens ectomycorrhizal associations. Soil contamination with DU resulted in the appearance of RISA bands of the ITS fragments of fungal metagenomic DNA that were characteristic of the genus Mortierella (Mortierellomycotina: Mucoromycota) in pine-free microcosms and for ectomycorrhizal fungi of the genus Scleroderma (Basidiomycota) in microcosms with mycorrhizal pines. The precise taxonomic affinity of the ITS fragments from the band appearing for non-mycorrhizal pines combined with DU remained uncertain, the most likely being related to the subphylum Zoopagomycotina. Thus, soil contamination by thermodynamically unstable metallic depleted uranium can cause a significant change in a soil fungal community under experimental conditions. These changes were also strongly affected by the presence of pine seedlings and their mycorrhizal status which impacted on DU biocorrosion and the release of bioavailable uranium species.

A 3-variable PDE model for predicting fungal growth derived from microscopic mechanisms

In this work, we present a new PDE model of the growth of Postia placenta, a species of brown rot fungus. The formulation was derived mainly from the biological mechanisms embedded in our discrete model, validated against experimental data. In order to mimic the growth mechanisms, we propose a new reaction-diffusion formulation, based on three variables: the concentration of tips, the branch density and the total hyphal density. The evolution of tips obeys a reaction-diffusion model, with constant diffusivity, while the evolution of the two other variables results from time integrals. The numerical solution is in excellent agreement with the averaged radial tip/hyphal densities of the mycelial network obtained by the discrete model. Thanks to the efficient exponential Euler method with Krylov subspace approximation, the solution needs only 3.5 s of CPU time to simulate 104-day of mycelium growth, in comparison with 8 hours for the discrete model. The great reduction of the RAM memory and computing time gives the possibility to upscale the simulation. The novelty of the PDE system is that the spatial colonization is formulated as a diffusion mechanism, which is self-standing, contrary to models based on an advection term. The continuous model can also reproduce the radial densities when the growth parameters in the discrete model are varied to adapt to different growth conditions. The correlation constructed between the two models provides us a tool for mutual insights between local biological mechanisms to the global biomass distribution, especially when analyzing experimental data.

Roles of saprotrophic fungi in biodegradation or transformation of organic and inorganic pollutants in co-contaminated sites

For decades, human activities, industrialization, and agriculture have contaminated soils and water with several compounds, including potentially toxic metals and organic persistent xenobiotics. The co-occurrence of those toxicants poses challenging environmental problems, as complicated chemical interactions and synergies can arise and lead to severe and toxic effects on organisms. The use of fungi, alone or with bacteria, for bioremediation purposes is a growing biotechnology with high potential in terms of cost-effectiveness, an environmental-friendly perspective and feasibility, and often representing a sustainable nature-based solution. This paper reviews different ecological, metabolic, and physiological aspects involved in fungal bioremediation of co-contaminated soils and water systems, not only addressing best methods and approaches to assess the simultaneous presence of metals and organic toxic compounds and their consequences on provided ecosystem services but also the interactions between fungi and bacteria, in order to suggest further study directions in this field.

Ecology of Contaminant Biotransformation in the Mycosphere: Role of Transport Processes

Fungi and bacteria often share common microhabitats. Their co-occurrence and coevolution give rise to manifold ecological interactions in the mycosphere, here defined as the microhabitats surrounding and affected by hyphae and mycelia. The extensive structure of mycelia provides ideal "logistic networks" for transport of bacteria and matter in structurally and chemically heterogeneous soil ecosystems. We describe the characteristics of the mycosphere as a unique and highly dynamic bacterial habitat and a hot spot for contaminant biotransformation. In particular, we emphasize the role of the mycosphere for (i) bacterial dispersal and colonization of subsurface interfaces and new habitats, (ii) matter transport processes and contaminant bioaccessibility, and (iii) the functional stability of microbial ecosystems when exposed to environmental fluctuations such as stress or disturbances. Adopting concepts from ecological theory, the chapter disentangles bacterial-fungal impacts on contaminant biotransformation in a systemic approach that interlinks empirical data from microbial ecosystems with simulation data from computational models. This approach provides generic information on key factors, processes, and ecological principles that drive microbial contaminant biotransformation in soil. We highlight that the transport processes create favorable habitat conditions for efficient bacterial contaminant degradation in the mycosphere. In-depth observation, understanding, and prediction of the role of mycosphere transport processes will support the use of bacterial-fungal interactions in nature-based solutions for contaminant biotransformation in natural and man-made ecosystems, respectively.

Build your own soil: exploring microfluidics to create microbial habitat structures

Soil is likely the most complex ecosystem on earth. Despite the global importance and extraordinary diversity of soils, they have been notoriously challenging to study. We show how pioneering microfluidic techniques provide new ways of studying soil microbial ecology by allowing simulation and manipulation of chemical conditions and physical structures at the microscale in soil model habitats.

The Mycelium as a Network

The characteristic growth pattern of fungal mycelia as an interconnected network has a major impact on how cellular events operating on a micron scale affect colony behavior at an ecological scale. Network structure is intimately linked to flows of resources across the network that in turn modify the network architecture itself. This complex interplay shapes the incredibly plastic behavior of fungi and allows them to cope with patchy, ephemeral resources, competition, damage, and predation in a manner completely different from multicellular plants or animals. Here, we try to link network structure with impact on resource movement at different scales of organization to understand the benefits and challenges of organisms that grow as connected networks. This inevitably involves an interdisciplinary approach whereby mathematical modeling helps to provide a bridge between information gleaned by traditional cell and molecular techniques or biophysical approaches at a hyphal level, with observations of colony dynamics and behavior at an ecological level.

The Geomycology of Elemental Cycling and Transformations in the Environment

Geomicrobiology addresses the roles of microorganisms in geological and geochemical processes, and geomycology is a part of this topic focusing on the fungi. Geoactive roles of fungi include organic and inorganic transformations important in nutrient and element cycling, rock and mineral bioweathering, mycogenic biomineral formation, and metal-fungal interactions. Lichens and mycorrhizas are significant geoactive agents. Organic matter decomposition is important for cycling of major biomass-associated elements, e.g., C, H, N, O, P, and S, as well as all other elements found in lower concentrations. Transformations of metals and minerals are central to geomicrobiology, and fungi affect changes in metal speciation, as well as mediate mineral formation or dissolution. Such mechanisms are components of biogeochemical cycles for metals as well as associated elements in biomass, soil, rocks, and minerals, e.g., S, P, and metalloids. Fungi may have the greatest geochemical influence within the terrestrial environment. However, they are also important in the aquatic environment and are significant components of the deep subsurface, extreme environments, and habitats polluted by xenobiotics, metals, and radionuclides. Applications of geomycology include metal and radionuclide bioleaching, biorecovery, detoxification, bioremediation, and the production of biominerals or metal(loid) elements with catalytic or other properties. Adverse effects include biodeterioration of natural and synthetic materials, rock and mineral-based building materials (e.g., concrete), cultural heritage, metals, alloys, and related substances and adverse effects on radionuclide mobility and containment. The ubiquity and importance of fungi in the biosphere underline the importance of geomycology as a conceptual framework encompassing the environmental activities of fungi.

Spiral and Rotor Patterns Produced by Fairy Ring Fungi

A broad class of soil fungi form the annular patterns known as ‘fairy rings’ and provide one of the only means to observe spatio-temporal dynamics of otherwise cryptic fungal growth processes in natural environments. We present observations of novel spiral and rotor patterns produced by fairy ring fungi and explain these behaviors mathematically by first showing that a well known model of fairy ring fungal growth and the Gray-Scott reaction-diffusion model are mathematically equivalent. We then use bifurcation analysis and numerical simulations to identify the conditions under which spiral waves and rotors can arise. We demonstrate that the region of dimensionless parameter space supporting these more complex dynamics is adjacent to that which produces the more familiar fairy rings, and identify experimental manipulations to test the transitions between these spatial modes. These same manipulations could also feasibly induce fungal colonies to transition from rotor/spiral formation to a set of richer, as yet unobserved, spatial patterns.

Automated image-based analysis of spatio-temporal fungal dynamics

Due to their ability to grow in complex environments, fungi play an important role in most ecosystems and have for that reason been the subject of numerous studies. Some of the main obstacles to the study of fungal growth are the heterogeneity of growth environments and the limited scope of laboratory experiments. Given the increasing availability of image capturing techniques, a new approach lies in image analysis. Most previous image analysis studies involve manual labelling of the fungal network, tracking of individual hyphae, or invasive techniques that do not allow for tracking the evolution of the entire fungal network. In response, this work presents a highly versatile tool combining image analysis and graph theory to monitor fungal growth through time and space for different fungal species and image resolutions. In addition, a new experimental set-up is presented that allows for a functional description of fungal growth dynamics and a quantitative mutual comparison of different growth behaviors. The presented method is completely automated and facilitates the extraction of the most studied fungal growth features such as the total length of the mycelium, the area of the mycelium and the fractal dimension. The compactness of the fungal network can also be monitored over time by computing measures such as the number of tips, the node degree and the number of nodes. Finally, the average growth angle and the internodal length can be extracted to study the morphology of the fungi. In summary, the introduced method offers an updated and broader alternative to classical and narrowly focused approaches, thus opening new avenues of investigation in the field of mycology.

Modelling combat strategies in fungal mycelia

Fungal mycelia have a well-established role in nutrient cycling and are widely used as agents in biological control and in the remediation of polluted landscapes. Competition and combat between different fungal communities is common in these contexts and its outcome impacts on local biodiversity and the success of such biotechnological applications. In this investigation a mathematical model representing mycelia as a system of partial differential equations is used to simulate combat between two fungal colonies growing into a nutrient-free domain. The resultant equations are integrated numerically and the model simulates well-established outcomes of combat between fungal communities. The outcome of pairwise combat is shown to depend on numerous factors including the suppression of advancing hyphae in rivals, the degradation of a rival's established biomass and the utilization and redistribution of available nutrient resources. It is demonstrated how non-transitive hierarchies in fungal communities can be established through switching mechanisms, mirroring observations reported in experimental studies, and how specialized defensive structures can emerge through changes in the redistribution of internal resources.

Growth-induced mass flows in fungal networks

Cord-forming fungi form extensive networks that continuously adapt to maintain an efficient transport system. As osmotically driven water uptake is often distal from the tips, and aqueous fluids are incompressible, we propose that growth induces mass flows across the mycelium, whether or not there are intrahyphal concentration gradients. We imaged the temporal evolution of networks formed by Phanerochaete velutina, and at each stage calculated the unique set of currents that account for the observed changes in cord volume, while minimizing the work required to overcome viscous drag. Predicted speeds were in reasonable agreement with experimental data, and the pressure gradients needed to produce these flows are small. Furthermore, cords that were predicted to carry fast-moving or large currents were significantly more likely to increase in size than cords with slow-moving or small currents. The incompressibility of the fluids within fungi means there is a rapid global response to local fluid movements. Hence velocity of fluid flow is a local signal that conveys quasi-global information about the role of a cord within the mycelium. We suggest that fluid incompressibility and the coupling of growth and mass flow are critical physical features that enable the development of efficient, adaptive biological transport networks.

Root growth models: towards a new generation of continuous approaches

Models of root system growth emerged in the early 1970s, and were based on mathematical representations of root length distribution in soil. The last decade has seen the development of more complex architectural models and the use of computer-intensive approaches to study developmental and environmental processes in greater detail. There is a pressing need for predictive technologies that can integrate root system knowledge, scaling from molecular to ensembles of plants. This paper makes the case for more widespread use of simpler models of root systems based on continuous descriptions of their structure. A new theoretical framework is presented that describes the dynamics of root density distributions as a function of individual root developmental parameters such as rates of lateral root initiation, elongation, mortality, and gravitropsm. The simulations resulting from such equations can be performed most efficiently in discretized domains that deform as a result of growth, and that can be used to model the growth of many interacting root systems. The modelling principles described help to bridge the gap between continuum and architectural approaches, and enhance our understanding of the spatial development of root systems. Our simulations suggest that root systems develop in travelling wave patterns of meristems, revealing order in otherwise spatially complex and heterogeneous systems. Such knowledge should assist physiologists and geneticists to appreciate how meristem dynamics contribute to the pattern of growth and functioning of root systems in the field.

Fungal mycelia allow chemotactic dispersal of polycyclic aromatic hydrocarbon-degrading bacteria in water-unsaturated systems

Contaminant biodegradation in soil is frequently limited by hindered physical access of bacteria to the contaminants. In the frame of the development of novel bioremediation approaches based on ecological principles, we tested the hypothesis that fungal networks facilitate the movement of bacteria by providing continuous liquid films in which gradients of chemoattractants can form and chemotactic swimming can take place. Unlike bacteria, filamentous fungi spread with ease in water-unsaturated soil. In a simple laboratory model of a water-unsaturated environment, we studied the movement of polycyclic aromatic hydrocarbon-degrading Pseudomonas putida PpG7 (NAH7) along a mycelium of Pythium ultimum. Some undirected dispersal was observed in the absence of a chemoattractant or when the non-chemotactic derivative strain P. putida G7.C1 (pHG100) was used. The bacterial movement became fourfold more effective and clearly directed when the chemotactic wild type was used and salicylate was present as a chemoattractant. No dispersal of bacteria was found in the absence of the fungus. These findings point at a role of mycelia for the translocation of chemicals and microorganisms. The results suggest that fungi improve the accessibility of contaminants in water-unsaturated environments.

Nonlinear colony extension of Sclerotinia minor and S. sclerotiorum

Fungal colonies initially extend exponentially and reach a constant linear extension rate determined solely by their growth in the peripheral zone. However the radial extension rates of Sclerotinia sclerotiorum and S. minor accelerate over time on PDA. Experiments were conducted to analyze the variable extension rates of the two Sclerotinia species and compare them with those of Verticillium dahliae and Cladosporium sp. In addition, the effects of starter disk size, disk position in the parent colony, the age of the parent colony, the concentration of potato dextrose broth and of incubation temperature also were determined. While the growth of Cladosporium sp. and V. dahliae followed established linear trends, the radial extension of S. sclerotiorum and S. minor colonies continuously accelerated over time until they reached the edge of the (150 mm diam) Petri dish. A polynomial model fitted the radial extension of colonies of Sclerotinia spp. Furthermore the accelerating colony extension rate was partly due to increasing colony radius. The rates of extension from mycelial disks transferred from the parental colony were positively correlated with the radius of the mycelial disks transferred. The rates of extension also were dependent on where the transferred disks were taken from parent colonies and the age and radius of the parent colony. On potato dextrose agar medium the extension rates of colonies of S. sclerotiorum and S. minor also were affected by broth concentration and temperature. With increasing nutrient concentration colony extension rates increased and were highest at 25 C. This study revealed a novel pattern of radial growth for Sclerotinia spp. that diverged from the established growth patterns of fungal colonies. Knowledge of the differences in growth behavior may be exploited in the laboratory studies on fungal competition and hyperparasitism and potentially in disease control strategies.

Modelling mycelial networks in structured environments

The growth habitat of most filamentous fungi is complex and displays a range of nutritional, structural, and temporal heterogeneities. There are inherent difficulties in obtaining and interpreting experimental data from such systems, and hence in this article a cellular automaton model is described to augment experimental investigation. The model, which explicitly includes nutrient uptake, translocation, and anastomosis, is calibrated for Rhizoctonia solani and is used to simulate growth in a range of three-dimensional domains, including those exhibiting soil-like characteristics. Results are compared with experimental data, and it is shown how the structure of the growth domain significantly influences key properties of the model mycelium. Thus, predictions are made of how environmental structure can influence the growth of fungal mycelia.

Growth model for arbuscular mycorrhizal fungi

In order to quantify the contribution of arbuscular mycorrhizal (AM) fungi to plant phosphorus nutrition, the development and extent of the external fungal mycelium and its nutrient uptake capacity are of particular importance. We develop and analyse a model of the growth of AM fungi associated with plant roots, suitable for describing mechanistically the effects of the fungi on solute uptake by plants. The model describes the development and distribution of the fungal mycelium in soil in terms of the creation and death of hyphae, tip-tip and tip-hypha anastomosis, and the nature of the root-fungus interface. It is calibrated and corroborated using published experimental data for hyphal length densities at different distances away from root surfaces. A good agreement between measured and simulated values was found for three fungal species with different morphologies: Scutellospora calospora (Nicol. & Gerd.) Walker & Sanders; Glomus sp.; and Acaulospora laevis Gerdemann & Trappe associated with Trifolium subterraneum L. The model and findings are expected to contribute to the quantification of the role of AM fungi in plant mineral nutrition and the interpretation of different foraging strategies among fungal species.

Modelling in vitro growth of dense root networks

Hairy roots are plants genetically transformed by Agrobacterium rhizogenes, which do not produce shoots and are composed mainly by roots. Hairy roots of Ophiorrhiza mungos Linn. are currently gaining interest of pharmacologists, since a secondary product of their metabolism, camptothecin, is used in chemotherapy. To optimize the production of valuable secondary metabolites it is necessary to understand the metabolism and growth of these roots systems. In this work, a mathematical model for description of apical growth of a dense root network (e.g. hairy roots) is derived. A continuous approach is used to define densities of root tips and root volume. Equations are posed to describe the evolution of these and are coupled to the distribution of nutrient concentration in the medium and inside the network. Following the principles of irreversible thermodynamics, growth velocity is defined as the sum over three different driving forces: nutrient concentration gradients, space gradients and root tip diffusion. A finite volume scheme was used for the simulation and parameters were chosen to fit experimental data from O. mungos Linn. hairy roots. Internal nutrient concentration determines short-term growth. Long-term behavior is limited by the total nutrient amount in the medium. Therefore, mass yield could be increased by guaranteeing a constant supply of nutrients. Increasing the initial mass of inoculation did not result in higher mass yields, since nutrient consumption due to metabolism also rose. Four different growth strategies are compared and their properties discussed. This allowed to understand which strategy might be the best to increase mass production optimally. The model is able to describe very well the temporal evolution of mass increase and nutrient uptake. Our results provide further understanding of growth and density distribution of hairy root network and therefore it is a sound base for future applications to describe, e.g., secondary metabolite production.

A fungal growth model fitted to carbon-limited dynamics of Rhizoctonia solani

Here, a quasi-steady-state approximation was used to simplify a mathematical model for fungal growth in carbon-limiting systems, and this was fitted to growth dynamics of the soil-borne plant pathogen and saprotroph Rhizoctonia solani. The model identified a criterion for invasion into carbon-limited environments with two characteristics driving fungal growth, namely the carbon decomposition rate and a measure of carbon use efficiency. The dynamics of fungal spread through a population of sites with either low (0.0074 mg) or high (0.016 mg) carbon content were well described by the simplified model with faster colonization for the carbon-rich environment. Rhizoctonia solani responded to a lower carbon availability by increasing the carbon use efficiency and the carbon decomposition rate following colonization. The results are discussed in relation to fungal invasion thresholds in terms of carbon nutrition.

Within-host disease ecology in the sea fan Gorgonia ventalina: modeling the spatial immunodynamics of a coral-pathogen interaction

We develop a spatially explicit model for the within-host interactions between a fungal pathogen and the immune response by its coral host. The model is parameterized for the recent epizootic of Aspergillus sydowii in the sea fan Gorgonia ventalina, but its structure is adaptable to many other diseases attacking corals worldwide, fungal infections in other invertebrates and plants, and opportunistic fungal infections in vertebrates. Model processes include pathogen growth and spread through consumption of host tissue, chemotactic attraction of undifferentiated host amoebocytes to infections, and amoebocyte differentiation into various cell types that attack the pathogen. Sensitivity analysis shows that the spread rate of a single localized infection is determined primarily by the pathogen's potential rate of host tissue consumption and by the host's ability to replenish the pool of undifferentiated amoebocytes and sustain a long-term response. The spatial localization of immune responses creates potentially strong indirect interactions between distant lesions, allowing new infections to grow rapidly while host resources are concentrated at older, larger infections. These findings provide possible mechanistic explanations for effects of environmental stressors (e.g., ocean warming, nutrient enrichment) on aspergillosis prevalence and severity and for the observed high spatial and between-host variability in disease impacts.

Geomycology: biogeochemical transformations of rocks, minerals, metals and radionuclides by fungi, bioweathering and bioremediation

The study of the role that fungi have played and are playing in fundamental geological processes can be termed 'geomycology' and this article seeks to emphasize the fundamental importance of fungi in several key areas. These include organic and inorganic transformations and element cycling, rock and mineral transformations, bioweathering, mycogenic mineral formation, fungal-clay interactions, metal-fungal interactions, and the significance of such processes in the environment and their relevance to areas of environmental biotechnology such as bioremediation. Fungi are intimately involved in biogeochemical transformations at local and global scales, and although such transformations occur in both aquatic and terrestrial habitats, it is the latter environment where fungi probably have the greatest influence. Within terrestrial aerobic ecosystems, fungi may exert an especially profound influence on biogeochemical processes, particularly when considering soil, rock and mineral surfaces, and the plant root-soil interface. The geochemical transformations that take place can influence plant productivity and the mobility of toxic elements and substances, and are therefore of considerable socio-economic relevance, including human health. Of special significance are the mutualistic symbioses, lichens and mycorrhizas. Some of the fungal transformations discussed have beneficial applications in environmental biotechnology, e.g. in metal leaching, recovery and detoxification, and xenobiotic and organic pollutant degradation. They may also result in adverse effects when these processes are associated with the degradation of foodstuffs, natural products, and building materials, including wood, stone and concrete. It is clear that a multidisciplinary approach is essential to understand fully all the phenomena encompassed within geomycology, and it is hoped that this review will serve to catalyse further research, as well as stimulate interest in an area of mycology of global significance.

The Development of Fungal Networks in Complex Environments

Fungi are of fundamental importance in terrestrial ecosystems playing important roles in decomposition, nutrient cycling, plant symbiosis and pathogenesis, and have significant potential in several areas of environmental biotechnology such as biocontrol and bioremediation. In all of these contexts, the fungi are growing in environments exhibiting spatio-temporal nutritional and structural heterogeneities. In this work, a discrete mathematical model is derived that allows detailed understanding of how events at the hyphal level are influenced by the nature of various environmental heterogeneities. Mycelial growth and function is simulated in a range of environments including homogeneous conditions, nutritionally-heterogeneous conditions and structurally-heterogeneous environments, the latter emulating porous media such as soils. Our results provide further understanding of the crucial processes involved in fungal growth, nutrient translocation and concomitant functional consequences, e.g. acidification, and have implications for the biotechnological application of fungi.

Biomass recycling and the origin of phenotype in fungal mycelia

Fungi are one of the most important and widespread components of the biosphere, and are essential for the growth of over 90% of all vascular plants. Although they are a separate kingdom of life, we know relatively little about the origins of their ubiquitous existence. This reflects a wider ignorance arising from their status as indeterminate organisms epitomized by extreme phenotypic plasticity that is essential for survival in complex environments. Here we show that the fungal phenotype may have its origins in the defining characteristic of indeterminate organisms, namely their ability to recycle locally immobilized internal resources into a mobilized form capable of being directed to new internal sinks. We show that phenotype can be modelled as an emergent phenomenon resulting from the interplay between simple local processes governing uptake and remobilization of internal resources, and macroscopic processes associated with their transport. Observed complex growth forms are reproduced and the sensitive dependence of phenotype on environmental context may be understood in terms of nonlinearities associated with regulation of the recycling apparatus.

Basidiomycete mycelia in forest soils: dimensions, dynamics and roles in nutrient distribution

Basidiomycete mycelia are ubiquitous in forest soils where they fulfil a range of key ecological functions. Population studies, based largely on basidiome collections, indicate that mycelia of many ectomycorrhizal and saprotrophic basidiomycetes can spread vegetatively for considerable distances through soil, but the extent to which these become physically or physiologically fragmented is unclear. This review considers aspects of the distribution, dynamics and translocatory activities of individual basidiomycete mycelia in forest soil, highlighting current gaps in our understanding and possible ways to address these.

Simulating colonial growth of fungi with the Neighbour-Sensing model of hyphal growth

The Neighbour-Sensing model brings together the basic essentials of hyphal growth kinetics into a vector-based mathematical model in which the growth vector of each virtual hyphal tip is calculated by reference to the surrounding virtual mycelium. The model predicts the growth pattern of many hyphae into three spatial dimensions and has been used to simulate complex fungal fruit body shapes. In this paper we show how the Neighbour-Sensing model can simulate growth in semi-solid substrata like agar or soil, enabling realistic simulation of mycelial colonies of filamentous fungi grown in 'Petri-dish style' experimental conditions. Newly implemented capabilities in the model include: a measurement and logging system within the program that maintains basic statistics about the mycelium it is simulating, this facilitates kinetic experimentation; inclusion of 'substrates' in the data space causing positive or negative tropisms for the growing mycelium; a horizontal plane tropism that provides a way of simulating colonies growing in or on a substratum like agar or soil by imposing a horizontal constraint on the data space the cyberhyphal tips can explore; three categories of hypha--standard hyphae are those that start the simulation, leading hyphae can emerge from the colony peripheral growth zone to take on a leading role, and secondary hyphae are branches that can arise late, far behind the peripheral growth zone, when mature hyphal segments resume branching to in-fill the older parts of the colony. We show how the model can be used to investigate hyphal growth kinetics in silico in experimental scenarios that would be difficult or impracticable in vivo. We also show that the Neighbour-Sensing model can generate sufficiently realistic cord-like structures to encourage the belief that this model is now sufficiently advanced for parameters to be defined that simulate specific in silico cyberfungi. The potential utility of these cyberspecies is that they provide a means to model the morphogenetic effects of a variety of factors, from environmental and nutritional features to mutations, in experimentally realistic situations, offering a valuable addition to the experimental toolkit of all those interested in fungal growth and morphology.

Reaction–diffusion models of development with state-dependent chemical diffusion coefficients

Reaction-diffusion models are widely used to model developmental processes. The great majority of current models invoke constant diffusion coefficients. However, the diffusion of metabolites or signals through tissues is frequently such that this assumption may reasonably be questioned. We consider several different physical mechanisms leading to effective diffusion coefficients in biological tissues which vary with the local conditions, including models in which juxtacrine signaling results in the diffusion of a signal in the absence of material transport. We develop a mathematical formalism for transforming local transport laws into diffusive terms. This procedure is appropriate when the typical length scale over which the concentrations change significantly is much greater than the dimensions of a cell. We review previous developmental models which considered the possibility of state-dependent diffusion coefficients. We also provide a few new motivating examples.

Mathematical modeling of regulatory mechanisms in yeast colony development

In the present study, yeast colony development serves as a model system to study growth of fungal populations with negligible nutrient and signal transport within the mycelium. Mathematical simulations address the question whether colony development is governed by diffusional limitation of nutrients. A hybrid one-dimensional cellular automaton model was developed that describes growth of discrete cells based upon microscopic interaction rules in a continuous field of nutrient and messenger. The model is scaled for the geometry of the experimental setup, cell size, growth- and substrate uptake rates. Therefore, calculated cell density profiles and nutrient distributions can be compared to experimental results and the model assumptions can be verified. In the physiologically relevant parameter range, simulations show an exponentially declining cell density along the median axis of the colonies in case of a diffusion limited growth scenario. These results are in good agreement with cell density profiles obtained in cultivations of the yeast Candida boidinii with glucose as the limiting carbon source but stand in contrast to the constant cell density profile estimated for Yarrowia lipolytica grown under the same conditions. While from the comparison of experimental results and simulations a diffusion limited growth mechanism is proposed for glucose limited C. boidinii colonies, this hypothesis is rejected for the growth of Y. lipolytica. As an alternative, a quorum sensing model was developed that can explain the evolution of constant cell density profiles based on the effect of a not further characterized unstable or volatile messenger.

Translocation of carbon by Rhizoctonia solani in nutritionally-heterogeneous microcosms

Responses of Rhizoctonia solani to spatial heterogeneity in sources of carbon, and associated translocation of carbon (C), were studied in a simple microcosm system comprising two discrete domains of agar gels separated on a glass slide and overlain with a porous membrane. Two arrangements of the gel pairs were used, one containing two equally large resources (representing 'homogeneous' conditions) and one containing a large and a negligible resource (representing 'heterogeneous' conditions). The nutrient sources were a standard mineral salt medium with or without glucose as sole C source. The fungus was inoculated onto one domain and growth responses determined by direct measurement of biomass. Translocation of C was quantified by use of 13C-enriched glucose. This substrate was either added to the agar at the outset, when studying newly developing colonies, or as a pulse into already established colonies. When growing in heterogeneous conditions, the fungus actively translocated C from a glucose-containing domain to sustain growth in the adjacent region lacking such a resource. In homogeneous conditions there was evidence of passive translocation (diffusion), but the fungus preferentially used local resource to maintain growth. Active translocation was only observed in newly growing colonies, whereas passive translocation occurred in both growing and established colonies. When the fungus was pulsed with a 13C-enriched glucose solution after 10 d growth, 2.5 times more 13C was taken up by the fungus grown in heterogeneous than homogeneous conditions, suggesting uptake exceeded local demands. In heterogeneous conditions, the total amount of 13C enriched glucose taken up by the fungus was independent of the location of the enriched glucose in the underlying medium. When the nylon membrane was replaced by Cellophane (an additional C source), degradation of the membrane and an increase in biomass occurred only in the heterogeneous system. The possible implications for these results in soil systems are discussed.