
==== Front
mSystems
mSystems
msystems
mSystems
2379-5077
American Society for Microbiology 1752 N St., N.W., Washington, DC

39189771
msystems00919-24
10.1128/msystems.00919-24
msystems.00919-24
Research Article
applied-and-industrial-microbiologyApplied and Industrial MicrobiologyFungal elemental profiling unleashed through rapid laser-induced breakdown spectroscopy (LIBS)
https://orcid.org/0000-0002-3207-1466
Rush Tomás A. 1 Conceptualization Data curation Formal analysis Funding acquisition Investigation Methodology Project administration Resources Software Supervision Validation Visualization Writing – original draft Writing – review and editing rushta@ornl.gov

Wymore Ann M. 1 Data curation Methodology Writing – review and editing
Rodríguez, Jr Miguel 1 Data curation Methodology Writing – original draft Writing – review and editing
Jawdy Sara 1 Methodology Resources Writing – review and editing
Vilgalys Rytas J. 2 Methodology Resources Writing – review and editing
Martin Madhavi Z. 1 Conceptualization Data curation Funding acquisition Investigation Methodology Visualization Writing – original draft Writing – review and editing
https://orcid.org/0000-0002-2091-9415
Andrews Hunter B. 3 Conceptualization Data curation Formal analysis Funding acquisition Investigation Methodology Project administration Resources Software Supervision Validation Visualization Writing – original draft Writing – review and editing andrewshb@ornl.gov

1 Biosciences Division, Oak Ridge National Laboratory , Oak Ridge, Tennessee, USA
2 Biology Department, Duke University , Durham, North Carolina, USA
3 Radioisotope Science and Technology Division, Oak Ridge National Laboratory , Oak Ridge, Tennessee, USA
Editor Gilbert Jack A. University of California San Diego , La Jolla, California, USA

Address correspondence to Tomás A. Rush, rushta@ornl.gov
Address correspondence to Hunter B. Andrews, andrewshb@ornl.gov
The authors declare no conflict of interest.

9 2024
27 8 2024
27 8 2024
9 9 e00919-2428 7 2024
30 7 2024
Copyright © 2024 Rush et al.
2024
Rush et al.
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

ABSTRACT

Elemental profiling of fungal species as a phenotyping tool is an understudied topic and is typically performed to examine plant tissue or non-biological materials. Traditional analytical techniques such as inductively coupled plasma–optical emission spectroscopy (ICP-OES) and inductively coupled plasma–mass spectrometry (ICP-MS) have been used to identify elemental profiles of fungi; however, these techniques can be cumbersome due to the difficulty of preparing samples. Additionally, the instruments used for these techniques can be expensive to procure and operate. Laser-induced breakdown spectroscopy (LIBS) is an alternative elemental analytical technique—one that is sensitive across the periodic table, easy to use on various sample types, and is cost-effective in both procurement and operation. LIBS has not been used on axenic filamentous fungal isolates grown in substrate media. In this work, as a proof of concept, we used LIBS on two genetically distinct fungal species grown on a nutrient-rich and nutrient-poor substrate media to determine whether robust elemental profiles can be detected and whether differences between the fungal isolates can be identified. Our results demonstrate a distinct correlation between fungal species and their elemental profile, regardless of the substrate media, as the same strains shared a similar uptake of carbon, zinc, phosphorus, manganese, and magnesium, which could play a vital role in their survival and propagation. Independently, each fungal species exhibited a unique elemental profile. This work demonstrates a unique and valuable approach to rapidly phenotype fungi through optical spectroscopy, and this approach can be critical in understanding these fungi's behavior and interactions with the environment.

IMPORTANCE

Historically, ionomics, the elemental profiling of an organism or materials, has been used to understand the elemental composition in waste materials to identify and recycle heavy metals or rare earth elements, identify the soil composition in space exploration on the moon or Mars, or understand human disorders or disease. To our knowledge, ionomic profiling of microbes, particularly fungi, has not been investigated to answer applied and fundamental biological questions. The reason is that current ionomic analytical techniques can be laborious in sample preparation, fail to measure all potential elements accurately, are cost-prohibitive, or provide inconsistent results across replications. In our previous efforts, we explored whether laser-induced breakdown spectroscopy (LIBS) could be used in determining the elemental profiles of poplar tissue, which was successful. In this proof-of-concept endeavor, we undertook a transdisciplinary effort between applied and fundamental mycology and elemental analytical techniques to address the biological question of how LIBS can used for fungi grown axenically in a nutrient-rich and nutrient-poor environment.

KEYWORDS

fungal behavior
elemental profiling
ionomics
endophytes
poplar
LIBS
U.S. Department of Energy (DOE) DE-AC05-00OR22725 Rush Tomás A. cover-dateSeptember 2024
==== Body
pmcINTRODUCTION

Nutrient acquisition, delivery, and availability dictate microbes' phenotype, lifestyle, and survival. Among the best microbes with which to study the influence of substrate availability is fungi because they are heterotrophic, meaning they cannot produce their food from environmental nutrients; thus, these organisms derive their nutritional requirements from other complex organic plants or animals (1). Tapping into this biological process by manipulating the available substrate has led to discoveries in identifying how fungi interact with a host, environment, or biological products such as alternative sources of protein and natural products (2, 3). Nutrient acquisition and the diverse roles it plays in controlling fungal behavior and producing biological products, as well as the elements obtained and how they differ across fungal species, are underexplored topics in the literature. However, in the past few years, interest has grown in using different analytical techniques for elemental profiling of fungal species.

Previous studies investigating fungal elemental profiling (also known as fungal ionomics) employed traditional analytical techniques such as inductively coupled plasma–optical emission spectroscopy [ICP-OES, also known as atomic emission spectroscopy (AES)] and inductively coupled plasma–mass spectrometry (ICP-MS) (Table 1). Unless laser ablation sampling is used, both of these analytical techniques require complete digestion of the samples, which makes sample preparation laborious (4). However, recently, a few studies have used single cell-inductively coupled plasma-time-of-flight mass spectrometry (SC-ICP-TOF-MS) to analyze individual yeast cell ionomic profiles that do not require digestion (5, 6). Unfortunately, ICP-MS, ICP-TOF-MS, and coupled laser ablation introduction systems can be cost-prohibitive. Even in the case of ICP-OES, the operational costs (e.g., power, argon gas) can be expensive.

TABLE 1 Previous analytical tools used to determine fungal ionomic profilesa

Method	Fungi examined	Number of elements detected	Reference	
HPLC with ICP-MS	Aspergillus fumigatus, Candida albicans, Cryptococcus neoformans, and Saccharomyces cerevisiae	10	7	
ICP-OES	Saccharomyces cerevisiae	12	8	
ICP-OES	Agaricus bisporus	6	9	
ICP-OES	Saccharomyces cerevisiae	13	10	
ICP-MS	Cordyceps kyushuensis	20	11	
ICP-MS	Hemileccinum impolitum, Butyriboletus roseoflavus, and Boletus umbriniporus	19	12	
ICP-MS	Saccharomyces cerevisiae	17	13	
ICP-MS	Schizophyllum commune	4	14	
ICP-MS	16 mushrooms	Various	15	
ICP-MS and ICP-OES	Perenniporia fraxinea	11	16	
ICP-OES	Suillus variegatus	18	17	
ICP-OES	Pleurotus spp.	31	18	
ICP-OES	30 edible mushrooms	15	19	
LIBS	Pleurotus opuntiae	8	20	
LIBS/SEM-EDS	Candida strains	4	21	
LIBS/Raman/MS	Unknown	Various	22	
SC-ICP-MS	Saccharomyces cerevisiae	Various	5	
SC-ICP-TOF-MS	Saccharomyces cerevisiae	6	6	
a Note: HPLC = high-performance liquid chromatography.

We posit that laser-induced breakdown spectroscopy (LIBS) can be a robust optical method for elemental profiling because it requires little-to-no sample preparation, can detect both light and heavy elements, has low procurement and operational costs, and can rapidly produce large data sets with proper replications (23, 24). Recently, a study examined how Pleurotus opuntiae absorbs and removes lead contamination using LIBS (20). However, to our knowledge, no studies have investigated the sole use of LIBS in an untargeted fashion to determine fungal elemental profiles. To address this knowledge gap, we used LIBS to identify the fungal ionomic profiles of two genetically different fungal species grown on nutrient-rich and nutrient-poor media. These measurements were performed in a manner that can be scaled for high-throughput analysis in future studies.

Most fungal ionomics or elemental studies have focused on how edible fungi uptake trace elements, heavy metals, or exhibit tolerance to certain elements that could be harmful to humans if consumed (9, 11, 12, 14–20, 25). Other studies focused on the genetic model yeast organism, Saccharomyces cerevisiae, to correlate genes to the results of elemental profiling (13) in order to test single-cell analytic techniques (5, 6); to establish protocols to examine fungal ionomics (8); to test other detection methods such as hydride generation atomic absorption spectroscopy (HGAAS), which did not produce reliable results (26); or to test strain variation for the uptake and tolerance of various elements (10). Rarely have studies focused on how fungal species or strains' elemental profiles change when grown on different substrate media to answer biological questions about their growth and development. One of the few investigations that has focused on differences in the fungal ionomics profiles between strains was that of Eide et al. (10). They examined 4,385 S. cerevisiae mutant strains to understand differences in their ability to uptake nutrients while grown on different substrate media. This study used ICP-OES to detect 13 elements, and 4.83% of strains had an altered yeast ionome profiles when grown on a nutrient-rich medium (i.e., yeast peptone dextrose). In a separate study, LIBS was coupled with SEM-EDS in a clinical study to discriminate between Candida strains by distinguishing only carbon, hydrogen, nitrogen, and oxygen emission spectra (21). Lastly, LIBS was paired with Raman spectroscopy and mass spectrometry to understand elements obtained by fungi and how they degrade specialized coatings (i.e., epoxy primer-based coating used in marine and aviation applications) (22). Taken together, the literature's lack of attention on the sole use of LIBS to determine the elemental profile of fungi and the significant variation in nutrient uptake observed within one fungal species, as shown in the yeast experiments, emphasizes the importance of comparative study between two genetically distinct species.

Previously, we used LIBS to analyze nutrient elements from fresh samples of poplar roots and soil and to investigate the plant nutrient transport seen under different environmental conditions (27–29); therefore, we predict that the same procedure can be leveraged profitably with fungal biomass to provide a robust ionomic profile. We hypothesize that LIBS can be used as a technique to (i) identify ionomic profile of axenic fungal cultures and (ii) understand how fungi ionomics profiles differ when grown in nutrient-poor or nutrient-rich media. Despite the challenges of using LIBS for simultaneous quantification of multiple elements and low-level (<100 ppm) quantification, it can still be a valuable tool for high-throughput screening. Because our goal is to determine whether LIBS can be used to distinguish axenic fungal cultures, we examined two genetically distinct fungi, Hyaloscypha finlandica and Mucor hiemalis, grown in a nutrient-rich medium called potato dextrose broth (PDB) and a nutrient-poor medium called glucose minimal medium (GMM) broth after 21 days of growth at 25°C. Both fungi were selected because (i) they are commonly isolated as endophytes from Populus trichocarpa x deltoides genotype 52-225 (30), (ii) we can easily grow axenic cultures in liquid broth and solid agar media, (iii) they have different physiological structures and growth rates, and (iv) they would presumably be the first representatives of fungal ionomics profiles. Our study will provide baseline data set results from axenic isolates, specifically those not derived from yeasts or mushrooms collected from a market or nature, which can be colonized by other microbes or infiltrated with debris or soil particles. H. finlandica is an ascomycete and dark septate endophyte that has a mutualistic lifestyle with poplar trees. M. hiemalis is a ubiquitous mucoromycete with a saprotrophic and pathogenic lifestyle in soil and humans. After 21 days of growth at 25°C, we identified 15 elements through LIBS.

MATERIALS AND METHODS

Isolates used and preparation of fungal biomass for LIBS analysis

H. finlandica (syn: Cadophora finlandica) strain PMI 746 and M. hiemalis strain PMI 3043 were used for this experiment. They were grown in sterile 125 mL Erlenmeyer screw cap flasks filled with 50 mL of PDB (BD Difco, New Jersey) or GMM without 1% thiamine (31) at 25°C in the dark at 250 rpm using an Eppendorf New Brunswick Excella E25 incubator shaker. There were five biological replications per condition. The fungi were grown for 21 days because H. finlandica grows slower compared to M. hiemalis in potato dextrose agar (PDA) and GMM agar (Fig. 1), and the goal was to compare the ionomics profiles equally. After 21 days, fungal biomass was collected in a 50 mL sterile Falcon tube and centrifuged at 4,000 rpm for 10 minutes, and then the supernatant was removed. Then, 25 mL of sterile Millipore water was added to the fungal biomass, vortexed for 2 minutes, and centrifuged again, and the supernatant was then removed. This process was carried out twice. After washing the fungal biomass, samples were placed in liquid nitrogen and stored at −80°C. Samples were lyophilized until dried, weighed, and autoclaved on a dry cycle at 121°C, 15 psi for 30 minutes before being removed from the containers, sectioned, and mounted on double-sided carbon tape (Nisshin EM. CO., LTD) affixed to glass slides. Samples were autoclaved to prevent viable fungal spores or hyphal from contaminating the LIBS system. Carbon tape with mounted samples was re-covered with tape backing and compressed for 48 hours prior to analysis.

Fig 1 Growth and development of M. hiemalis and H. finlandica on potato dextrose agar (PDA) or broth (PDB) and glucose minimal medium (GMM) agar or broth. (a–d) Macroscopic photos display different physiologies and color appearances of the same fungi across different substrate media. Scale bar is 1 cm. Hyphal development on the edges (e–h), and scale bar is 1 mm. Closer image of hyphal development (i–l), and scale bar is 200 µm. (m) Growth rate of the fungi on different growth media, (n) development of hyphal tips in H. finlandica, (o) sporulation in M. hiemalis, (p) H. finlandica biomass, and (q) M. hiemalis biomass growth in each substrate that was used for the ionomic profiling. M. hiemalis did not produce distinguishable hyphal tips to be counted, and H. finlandica did not produce spores in our experiment. For growth rate, one-way ANOVA was significant (P-value = 0.0001), and Tukey's post-hoc analysis was used, where different letters indicated significant differences. For hyphal tips and sporulation counts, a Student's t-test was used, where (*) indicates P-value < 0.05; (**) indicates P-value < 0.01; and (ns) indicates not significant. There were five biological replications per condition. The error bars in all figures indicate the standard error of the mean.

Growth of Mucor hiemalis and Hyaloscypha finlandica on PDA and GMM media through various stages and microscopic views. It includes growth area over days, number of hyphal tips, spore count, and dry biomass weight on different media.

Fungal physiology analysis

Both fungi were grown on PDA (BD Difco, New Jersey; cat. no. 247940) for 30 days at 25°C in the dark, prior to the experiment. Afterwards, a sterile core borer with a diameter of 0.5 cm2 was used to cut fungal plugs. The plugs were transferred to PDA or GMM with agar (Sigma-Aldrich; cat. no. 05040) to take physiology measurements. The growth rate was determined by measuring the diameter of the fungi 1-day post-inoculation (dpi), 2 dpi, 3 dpi, 7 dpi, and 14 dpi at 25°C in the dark. The plug's 0.5 cm2 diameter was excluded from the final growth area. Macroscopic and microscopic photos were taken with a Zeiss SteREO Discovery v12. There were three biological replications per fungus and substrate media. Sporangiospores were collected from M. hiemalis grown on PDA and GMM. Next, 5 mL of sterile Millipore water was placed in each plate, and a sterile L-shaped spreader was used to disrupt the sporangiospores. Then, the sporangiospores were collected in a sterile 15 mL Falcon tube and vortexed for 3 minutes at maximum speed. Ten microliters of spore suspension was collected and placed in the Countess Cell Counting Chamber Slide and analyzed by the Countess II FL Automated Cell Counter. We had two technical replications per biological replication. Spore counting was not performed for H. finlandica because spores were not produced in culture. Hyphal tips were counted and analyzed for H. finlandica using a ZEISS SteREO Discovery v12. Every visible hyphal tip within a frame of the image was counted. There were three technical replications per biological replication. Hyphal tips were not analyzed for M. hiemalis because no distinguishable tips were observed.

LIBS analysis

Measurements were performed using a LIBS-8 module from Applied Photonics; this instrument utilizes a 1,064 nm Nd:YAG laser passed through beam expansion optics prior to focusing onto the sample to provide a tighter spot size of approximately 100 µm. A laser pulse energy of 100 mJ was used to provide strong signal-to-noise levels. Prior to testing samples, metal coupons were measured to ensure proper LIBS performance. Each sample was tested in a 4 × 4 grid with one shot per spot, providing 16 spectra per sample. The plasma light was collected using eight collimating optics located above the sample, off-axis from the impending laser pulse. The light was transported to an eight-channel spectrometer with a range of 178–1,022 nm. This spectrometer bank was operated with a delay time of 3 µs and an integration time of 50 µs.

The results of every four shots were averaged to provide four spectra per sample, with five technical replicates, leading to 20 spectra per biological specimen. The spectra were inspected to identify all elemental emission peaks, and the background-corrected peak intensities were tabulated for each spectrum. The elemental emissions were compared to each other to investigate the Pearson correlation coefficients, providing insight into the unique fungi sample profiles. A positive correlation coefficient indicates two variables having a positive linear correlation; conversely, a negative correlation coefficient indicates a negative linear correlation. The linear correlation between two variables is categorized as small, medium, or large if the magnitude (absolute value) of the correlation coefficient value is 0.1–0.3, 0.3–0.5, or 0.5–1.0, respectively. A magnitude value of <0.1 is categorized as having no linear correlation. Principal component analysis (PCA) was applied to the collected spectra, and the principal component (PC) loadings were compared to the elemental profiles to identify the behavior captured by each PC. Lastly, the PCA transformed data were modeled using a linear classification model to distinguish fungi samples using only LIBS elemental profiles.

Statistical analysis

Fungal physiology statistical analyses were performed using the GraphPad Prism software version 10.2.0 (392) (GraphPad, San Diego, California). For fungal growth rates, one-way analysis of variance (ANOVA) was performed per time point and condition. A Tukey's multiple comparison test was performed per time point and condition. Similar letters indicate no significant differences, whereas different letters indicate significant differences (P-value < 0.05). For the number of hyphal or spores counted, a paired t-test was performed. For the dry biomass comparison, an unpaired t-test was performed within a species and between the two-substrate media. For both t-test analyses, asterisks indicate significant differences (P-value < 0.05). The error bars in all figures indicate the standard error of the mean.

RESULTS AND DISCUSSION

First, we determined differences in growth rates, biomass weight, sporulation, number of hyphal branching tips, and appearances in the nutrient-rich and nutrient-poor media between the fungal species (Fig. 1).

Next, we compared the averaged and normalized LIBS spectra for the different sample groups (Fig. 2a). At the larger scale, the LIBS profiles appear nearly identical, but unique spectral fingerprints begin to emerge upon closer investigation (see Fig. 2b and c). The spectra were parsed from 190 nm to 900 nm to identify all element signatures within the sample set. The elements identified include carbon, zinc, phosphorous, manganese, magnesium, silicon, iron, calcium, aluminum, sodium, hydrogen, lithium, potassium, oxygen, and nitrogen.

Fig 2 The averaged and normalized LIBS spectra for each fungi species in the different growth media are shown with (a) an offset to compare the labeled major emission peaks. Multiple peaks may be seen for a single element due to the branching electronic transitions as the plasma cools. A notional diagram of a LIBS measurement and a closer look at the ultraviolet carbon, zinc, and phosphorus emissions are shown in (b). A similar zoomed plot of the 370–400 nm region (c) shows the iron, magnesium, calcium, and aluminum elemental emissions and the carbon–nitrogen (CN) molecular bands. The average relative emission intensities are compared for each element in (d). Data are based on the average of five biological replications per condition.

Normalized intensity LIBS spectra for elemental analysis of Mucor hiemalis and Hyaloscypha finlandica grown on PDB and GMM media. It also depicts relative intensities for various elements for the various samples.

The relative emission intensities of all these elements were compared to one another to investigate trends in the elemental profiles of the fungi species and to discern any effects of the growth media. Several trends were identified when comparing the fungi species. H. finlandica showed consistently more intense emissions from zinc, phosphorus, manganese, magnesium, iron, hydrogen, and oxygen, regardless of the growth media. The mutually elevated hydrogen and oxygen levels likely indicate that H. finlandica retains moisture to a greater extent than M. hiemalis, regardless of the growth media. M. hiemalis exhibited more intense lithium and calcium emissions. M. hiemalis showed lower aluminum levels than H. finlandica grown in the same media. Higher calcium and aluminum levels were detected in fungi grown in PDB, while specimens grown in GMM had more intense carbon emissions.

The emission signal intensities for each sample type were compared to one another through Pearson correlation analysis to investigate how elemental signatures related to one another. The Pearson correlation coefficients for M. hiemalis (Fig. 3a) and H. finlandica (Fig. 3b) are visualized in Fig. 3. For both M. hiemalis and H. finlandica, the correlation coefficients indicate strong positive correlations with carbon, zinc, phosphorus, manganese, and magnesium emissions. The positive correlations between these elements in both species may indicate their vital role in fungi propagation and survival. For simplicity, this group of elements (i.e., carbon, zinc, phosphorus, manganese, and magnesium) is referred to as essential elements for the remainder of the correlation analysis. Similar to the essential elements, hydrogen, nitrogen, and oxygen exhibit strong positive correlations regardless of the species, but this result is likely due to contributions from both sample moisture and the atmospheric contributions from these elements when LIBS is not performed under an inert cover gas. When the correlation coefficients of each fungi species are compared, a few noticeable differences are seen. Firstly, H. finlandica exhibits strong positive correlations among sodium, hydrogen, and the essential element group. This indicates a reliance of H. finlandica on sodium that M. hiemalis does not exhibit. Similar behavior was observed for potassium in H. finlandica, but generally, a medium positive correlation exists between the essential elements and potassium. M. hiemalis shows a robust negative correlation between potassium and essential elements. It shows stronger positive correlations among silicon, iron, and the vital elements; although H. finlandica shows a positive correlation between silicon and calcium, M. hiemalis does not. Both fungi species show negative correlations between lithium and the essential elements, although M. hiemalis shows a stronger negative relationship. This likely indicates that lithium harms fungal growth. The last significant difference between the two fungi species is M. hiemalis' strong positive correlation among aluminum with zinc, magnesium, and iron.

Fig 3 Pearson correlation coefficients between elemental emissions of (a) M. hiemalis and (b) H. finlandica. Each cell represents the Pearson correlation coefficient between the row and column labels. Positive correlations are shown in red, and negative correlations are shown in blue.

Correlation heatmaps for elemental relationships in Mucor hiemalis and Hyaloscypha finlandica. Values indicate the strength and direction of correlations between different elements for each fungus.

In addition to correlation analysis on the emission intensities, PCA was applied to the entire LIBS spectrum to further investigate the key drivers for the sample set's variance. PCA score plots are shown in Fig. 4, with the markers colored by fungi species and growth media. Five PCs were used to explain a total of 99.11% of the signal variance in the data set. The score plots show that the samples form clusters around their fungi species and growth media—meaning PCA loadings can be used to decipher which signals are heavily impacted by these factors. PC1 clearly captures the variance related to the growth media used (Fig. 4a), where a negative PC1 score indicates PDB and a positive score points to GMM. Together, PC1 and PC3 begin to cluster the various fungi species, particularly M. hiemalis and H. finlandica grown in PDB, as seen in Fig. 4b. The score plot shown in Fig. 4c reveals that PC4 is a strong indicator of the fungal species.

Fig 4 PCA score plots. Markers are colored by the sample class, and markers indicate the growth media; the PCA explained that variance ratios were determined to be as follows: PC1, 87.2%; PC2, 8.11%; PC3, 2.06%; PC4, 1.17%; and PC5, 0.57%.

Three PCA score plots depict clusters and separations based on PC1, PC2, and PC3, highlighting the differences in elemental composition and growth conditions between the fungi on different media.

To further investigate the individual PCs, the PCA scores for each sample can be added to the data matrix and correlated directly with the emission intensities. The Pearson correlation coefficient values relating the PCA scores to the emission intensities are tabulated in Fig. 5.

Fig 5 Pearson correlation coefficients between elemental emissions and PCs. Each cell represents the Pearson correlation coefficient between the row and column labels. Positive correlations are shown in red, and negative correlations are shown in blue.

Correlation heatmap depicts relationships between various elements and principal components. Correlation values indicate the strength and direction of each relationship.

PC1 was observed to correlate strongly with sodium, hydrogen, oxygen, nitrogen, and essential elements. It shows a strong negative correlation with calcium and medium negative correlations with lithium and the CN molecular emission. Because a negative PC1 score can be used to discern when the samples are grown in PDB, we would expect more intense calcium, lithium, and CN emissions in the PDB samples. This behavior is verified in the spectra shown in Fig. 2d. PC4, which largely separates fungal species, shows strong positive correlations with manganese, iron, and aluminum and medium positive correlations with magnesium and calcium. Accordingly, M. hiemalis grown in GMM has the strongest manganese, iron, and aluminum emissions (see Fig. 2d), and M. hiemalis consistently has greater magnesium and calcium emissions than H. finlandica grown in the same media.

The PCA analysis allows for a more exhaustive interpretation of the elemental profiles of the samples tested. Both growth media and fungi species can be distinguished visually on score plots. Therefore, the final step was to test the ability to classify fungi samples based on LIBS elemental spectra. The sample spectra (n = 80) were randomly split into thirds: one-third was used to train the PCA classifier, and two-thirds were used to test the classifier. PCA classifier models were built to simultaneously predict growth media, fungi species, or both. The models predicted the test group classes with 100% accuracy for both the growth media and fungi species classifiers. The combined classifier predicted the test group classes with 98% accuracy; only one M. hiemalis grown in GMM was falsely predicted to be H. finlandica grown in PDB. If the average number of shots is increased from four to eight, then the classifier accuracy rises to 100%.

Linking the phenotype to ionomic profiles would indicate that the elemental fold-increase in aluminum, carbon, hydrogen, iron, manganese, magnesium, nitrogen, oxygen, phosphorus, and zinc accumulated in the M. hiemalis biomass when grown in GMM may have influenced its faster growth rate, transparent appearance, and reduced biomass and sporulation in comparison to PDA/PDB. Conversely, the fold-increase in calcium, lithium, and potassium in M. hiemalis biomass when grown in PDA or PDB could have been the cause for its slower growth rate, opaque yellow appearance, and increased biomass and sporulation when compared to GMM. For H. finlandica, differences in nutrient-rich or nutrient-poor media did not produce a significant difference between growth rate and biomass. However, when H. finlandica grew in PDB, a fold-increase of calcium and lithium could have impacted the increased hyphal tips, its uniform appearance, and the gray opaque appearance observed when grown on PDA. Interestingly, more increase-fold of elements like aluminum, carbon, hydrogen, iron, manganese, magnesium, nitrogen, oxygen, phosphorus, potassium, sodium, and zinc could have negatively influenced the growth rate and the number of hyphal tips and its sporadic formation. Moreover, these elements could have contributed to the yellow-green opaque appearance of the fungus.

Conclusion

This study demonstrates that LIBS can be used successfully to generate rapid, high-quality elemental fingerprints for fungal biomass with minimal sample preparation. To our knowledge, this is the first study to examine endophytic fungi grown axenically in liquid broth media of various substrates and that were not yeasts or mushrooms collected from nature or purchased from a market. LIBS can detect fungal elemental composition to the parts per million (ppm) level with minimal sample preparation. On average, between 0.2 and 0.6 mg of fungal biomass was needed to obtain the data sets. We concluded that elemental profiles differ based on the nutrient-rich or nutrient-poor environments or based on the fungal taxonomic assignment for calcium, hydrogen, iron, sodium, and silicon. Despite the nutrient-rich or nutrient-poor environments and taxonomically different fungal species, both fungi will uptake carbon, zinc, phosphorus, manganese, and magnesium, indicating that these elements may be vital for their survival and propagation. The results of this screening exercise could be used to plan future studies targeted at the impacts of specific nutrients.

Although we believe that LIBS will be a tool of the future for high-throughput ionomic profiling of fungi, there are inherent limitations. The main limit is that when used as a high-throughput screening tool as in this study, LIBS is only qualitative. For full quantification, matrix matched standards would be needed, which would likely lead to fungal samples being measured using a standard addition (24). This would require more fungal biomass, limit throughput, and result in only a few elements quantified at a time. Fortunately, qualitative elemental profiles still provide powerful insight into fungal behavior in different environments and additional analytical techniques could be leveraged for further investigation after LIBS screening as required.

This successful demonstration that LIBS can provide biological context about the mechanisms by which a microbe may acquire elements within its biomass sets the stage for future investigations using LIBS as a high-throughput phenotypic tool to provide insights into the elements required for fungal survival and propagation. We are currently undertaking experiments to determine whether ionomic profiling is conserved between species within the same taxonomic class; whether there is a variation among profiles within isolates of the same species; whether fungi that are classified as mutualists, saprotrophs, or pathogens cluster together when growing on the same substrate media; and whether those classifications merge when they are grown on nutrient-poor media.

ACKNOWLEDGMENTS

The authors would like to acknowledge Christopher Orosco for assistance with graphics and John Batson III for editing the manuscript.

This work was supported by the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the U.S. Department of Energy under contract DE-AC-05-000R22725. This research was also funded by the Genomic System Sciences Program, U.S. Department of Energy, Office of Science, Biological and Environmental Research. Fungal isolates were obtained from Plant-Microbe Interface Scientific Focus Area microbial collection (ORNL; http://pmi.ornl.gov).

Notice: This manuscript has been authored by UT-Battelle, L.L.C., under contract DE-AC05-00OR22725 with the US Department of Energy (DOE). The US government retains and the publisher, by accepting the article for publication, acknowledges that the US government retains a non-exclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript or allow others to do so, for US government purposes. DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan).

DATA AVAILABILITY

Raw data for fungal physiology (Fig. 1) are found in Table S1 and compared elemental compositions (Fig. 2) are found in Table S2. Fungal strains are available upon request.

SUPPLEMENTAL MATERIAL

The following material is available online at https://doi.org/10.1128/msystems.00919-24.

10.1128/msystems.00919-24.SuF1 Captions msystems.00919-24-s0001.docx

Captions for supplemental tables.

10.1128/msystems.00919-24.SuF2 Table S1 msystems.00919-24-s0002.xlsx

Raw data for fungal physiology.

10.1128/msystems.00919-24.SuF3 Table S2 msystems.00919-24-s0003.xlsx

Raw data for compared elemental emissions.

ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.
==== Refs
REFERENCES

1 Naranjo-Ortiz MA, Gabaldón T. 2019. Fungal evolution: major ecological adaptations and evolutionary transitions. Biol Rev Camb Philos Soc 94 :1443–1476. doi:10.1111/brv.12510 31021528
2 Keller NP. 2019. Fungal secondary metabolism: regulation, function and drug discovery. Nat Rev Microbiol 17 :167–180. doi:10.1038/s41579-018-0121-1 30531948
3 Rush TA, Shrestha HK, Gopalakrishnan Meena M, Spangler MK, Ellis JC, Labbé JL, Abraham PE. 2021. Bioprospecting Trichoderma: a systematic roadmap to screen genomes and natural products for biocontrol applications. Front Fungal Biol 2 :716511. doi:10.3389/ffunb.2021.716511 37744103
4 Hansen TH, de Bang TC, Laursen KH, Pedas P, Husted S, Schjoerring JK. 2013. Multielement plant tissue analysis using ICP spectrometry. Methods Mol Biol 953 :121–141. doi:10.1007/978-1-62703-152-3_8 23073880
5 Pereira JSF, Álvarez-Fernández García R, Corte-Rodríguez M, Manteca A, Bettmer J, LeBlanc KL, Mester Z, Montes-Bayón M. 2023. Towards single cell ICP-MS normalized quantitative experiments using certified selenized yeast. Talanta 252 :123786. doi:10.1016/j.talanta.2022.123786 36027616
6 Qin W, Stärk HJ, Müller S, Reemtsma T. 2022. Exploring the extent of phosphorus and heavy metal uptake by single cells of Saccharomyces cerevisiae and their effects on intrinsic elements by SC-ICP-TOF-MS. Front Microbiol 13 :870931. doi:10.3389/fmicb.2022.870931 35547146
7 Wehmeier S, Morrison E, Plato A, Raab A, Feldmann J, Bedekovic T, Wilson D, Brand AC. 2020. Multi trace element profiling in pathogenic and non-pathogenic fungi. Fungal Biol 124 :516–524. doi:10.1016/j.funbio.2020.03.001 32389315
8 Thines L, Iserentant A, Morsomme P. 2020. Determination of the cellular ion concentration in Saccharomyces cerevisiae using ICP-AES. Bio Protoc 10 :e3727. doi:10.21769/BioProtoc.3727
9 Ozcan MM, Dursun N, Al Juhaimi FY. 2013. Heavy metals intake by cultured mushrooms growing in model system. Environ Monit Assess 185 :8393–8397. doi:10.1007/s10661-013-3181-8 23591676
10 Eide DJ, Clark S, Nair TM, Gehl M, Gribskov M, Guerinot ML, Harper JF. 2005. Characterization of the yeast ionome: a genome-wide analysis of nutrient mineral and trace element homeostasis in Saccharomyces cerevisiae. Genome Biol 6 :R77. doi:10.1186/gb-2005-6-9-r77 16168084
11 Zhang G, Zhao Y, Liu F, Ling J, Lin J, Zhang C. 2013. Determination of essential and toxic elements in Cordyceps kyushuensis Kawam by inductively coupled plasma mass spectrometry. J Pharm Biomed Anal 72 :172–176. doi:10.1016/j.jpba.2012.08.007 23102516
12 Zhang J, Barałkiewicz D, Hanć A, Falandysz J, Wang Y. 2020. Contents and health risk assessment of elements in three edible ectomycorrhizal fungi (Boletaceae) from polymetallic soils in Yunnan province, SW China. Biol Trace Elem Res 195 :250–259. doi:10.1007/s12011-019-01843-y 31363991
13 Yu D, Danku JMC, Baxter I, Kim S, Vatamaniuk OK, Vitek O, Ouzzani M, Salt DE. 2012. High-resolution genome-wide scan of genes, gene-networks and cellular systems impacting the yeast ionome. BMC Genomics 13 :623. doi:10.1186/1471-2164-13-623 23151179
14 Traxler L, Shrestha J, Richter M, Krause K, Schäfer T, Kothe E. 2022. Metal adaptation and transport in hyphae of the wood-rot fungus Schizophyllum commune. J Hazard Mater 425 :127978. doi:10.1016/j.jhazmat.2021.127978 34896706
15 Tel-Çayan G, Ullah Z, Öztürk M, Yabanlı M, Aydın F, Duru ME. 2018. Heavy metals, trace and major elements in 16 wild mushroom species determined by ICP-MS. AT Spectrosc 39 :29–37. doi:10.46770/AS.2018.01.004
16 Sturini M, Girometta C, Maraschi F, Savino E, Profumo A. 2017. A preliminary investigation on metal bioaccumulation by Perenniporia fraxinea. Bull Environ Contam Toxicol 98 :508–512. doi:10.1007/s00128-017-2038-1 28204838
17 Szubstarska J, Jarzyńska G, Falandysz J. 2012. Trace elements in variegated bolete (Suillus variegatus) fungi. Chem Pap 66 :1026–1031. doi:10.2478/s11696-012-0216-5
18 Siwulski M, Mleczek M, Rzymski P, Budka A, Jasińska A, Niedzielski P, Kalač P, Gąsecka M, Budzyńska S, Mikołajczak P. 2017. Screening the multi-element content of pleurotus mushroom species using inductively coupled plasma optical emission spectrometer (ICP-OES). Food Anal Methods 10 :487–496. doi:10.1007/s12161-016-0608-1
19 Rasalanavho M, Moodley R, Jonnalagadda SB. 2019. Elemental distribution including toxic elements in edible and inedible wild growing mushrooms from South Africa. Environ Sci Pollut Res 26 :7913–7925. doi:10.1007/s11356-019-04223-0
20 Yadav P, Mishra V, Kumar T, Rai AK, Gaur A, Singh MP. 2023. An approach to evaluate Pb tolerance and its removal mechanisms by Pleurotus opuntiae. J Fungi (Basel) 9 :405. doi:10.3390/jof9040405 37108860
21 Manzoor S, Ugena L, Tornero-Lopéz J, Martín H, Molina M, Camacho JJ, Cáceres JO. 2016. Laser induced breakdown spectroscopy for the discrimination of Candida strains. Talanta 155 :101–106. doi:10.1016/j.talanta.2016.04.030 27216662
22 Huang X, Liu L, Lu Y, Dong H, Mao A, Li P, Cui B, Silvain J-F, Lu Y. 2023. Chemical element variation in fungi-induced coating degradation using laser-induced breakdown spectroscopy combined with Raman spectroscopy, mass spectrometry, and multivariate analyses. J Anal At Spectrom 38 :1668–1675. doi:10.1039/D3JA00001J
23 Hahn DW, Omenetto N. 2012. Laser-induced breakdown spectroscopy (LIBS), part II: review of instrumental and methodological approaches to material analysis and applications to different fields. Appl Spectrosc 66 :347–419. doi:10.1366/11-06574 22449322
24 Andrews HB, Wymore AM, Wetter EE, Herndon EM, Li H, Martin SA, Griffiths NA, Yang X, Muchero W, Weston DJ, Martin MZ. 2023. Rapid screening of wood and leaf tissues: investigating silicon-based phytoliths in Populus trichocarpa for carbon storage applications using laser-induced breakdown spectroscopy and scanning electron microscopy–energy dispersive X-ray spectroscopy. J Anal At Spectrom 38 :2353–2364. doi:10.1039/D3JA00186E
25 Řanda Z, Kučera J. 2004. Trace elements in higher fungi (mushrooms) determined by activation analysis. J Radioanal Nucl Chem 259 :99–107. doi:10.1023/B:JRNC.0000015813.27926.32
26 Renard NE, Tompkins TA. 2002. Evaluation of methods for total selenium determination in yeast. Biol Trace Elem Res 88 :185–191. doi:10.1385/BTER:88:2:185 12296428
27 Martin M, Brice D, Martin S, André N, Labbé N. 2021. Inorganic characterization of switchgrass biomass using laser-induced breakdown spectroscopy. Spectrochimica Acta Part B: Atomic Spectroscopy 186 :106323. doi:10.1016/j.sab.2021.106323
28 Martin MZ, Glasgow DC, Tschaplinski TJ, Tuskan GA, Gunter LE, Engle NL, Wymore AM, Weston DJ. 2017. Correlating laser-induced breakdown spectroscopy with neutron activation analysis to determine the elemental concentration in the ionome of the Populus trichocarpa leaf. Spectrochim Acta Part B: At Spectrosc 138 :46–53. doi:10.1016/j.sab.2017.10.008
29 Andrews HB, Martin MZ, Wymore AM, Kalluri UC. 2024. Rapid in situ nutrient element distribution in plants and soils using laser-induced breakdown spectroscopy(LIBS). Plant Soil 495 :3–12. doi:10.1007/s11104-023-05988-7
30 Bonito G, Hameed K, Ventura R, Krishnan J, Schadt CW, Vilgalys R. 2016. Isolating a functionally relevant guild of fungi from the root microbiome of Populus. Fungal Ecol 22 :35–42. doi:10.1016/j.funeco.2016.04.007
31 Kaminskyj SGW. 2001. Fundamentals of growth, storage, genetics and microscopy of Aspergillus nidulans. Fungal Genet Rep 48 :25–31. doi:10.4148/1941-4765.1175
