Insights into the Ecological Diversification of the Hymenochaetales based on Comparative Genomics and Phylogenomics With an Emphasis on Coltricia.

To elucidate the genomic traits of ecological diversification in the Hymenochaetales, we sequenced 15 new genomes, with attention to ectomycorrhizal (EcM) Coltricia species. Together with published data, 32 genomes, including 31 Hymenochaetales and one outgroup, were comparatively analyzed in total. Compared with those of parasitic and saprophytic members, EcM species have significantly reduced number of plant cell wall degrading enzyme genes, and expanded transposable elements, genome sizes, small secreted proteins, and secreted proteases. EcM species still retain some of secreted carbohydrate-active enzymes (CAZymes) and have lost the key secreted CAZymes to degrade lignin and cellulose, while possess a strong capacity to degrade a microbial cell wall containing chitin and peptidoglycan. There were no significant differences in secreted CAZymes between fungi growing on gymnosperms and angiosperms, suggesting that the secreted CAZymes in the Hymenochaetales evolved before differentiation of host trees into gymnosperms and angiosperms. Nevertheless, parasitic and saprophytic species of the Hymenochaetales are very similar in many genome features, which reflect their close phylogenetic relationships both being white rot fungi. Phylogenomic and molecular clock analyses showed that the EcM genus Coltricia formed a clade located at the base of the Hymenochaetaceae and divergence time later than saprophytic species. And Coltricia remains one to two genes of AA2 family. These indicate that the ancestors of Coltricia appear to have originated from saprophytic ancestor with the ability to cause a white rot. This study provides new genomic data for EcM species and insights into the ecological diversification within the Hymenochaetales based on comparative genomics and phylogenomics analyses.

Revising the taxonomic placement of Laetiporus persicinus within the Laetiporaceae.

The fungus currently known as Laetiporus persicinus is a recognizable brown-rot decayer that is widespread on oak hosts in the southeastern United States. This species was first described as Polyporus persicinus in 1872 based on collections by Henry W. Ravenel from South Carolina. In this study, we elucidate the phylogenetic relationships of Laetiporus persicinus based on maximum likelihood and Bayesian inference analyses of a four-locus data set (18S, 28S, rpb2, and tef1) from taxa within the Fomitopsidaceae and Laetiporaceae. The internal transcribed spacer (ITS) region was analyzed separately because it was not possible to align this locus across a diverse data set that included taxa from multiple families. Our analysis and previous studies indicate that Laetiporus persicinus does not belong to Laetiporus sensu stricto, and we found a strongly supported relationship between Laetiporus persicinus and the African species Kusaghiporia usambarensis, despite the fact that the 28S phylogeny resolved a different (but unsupported) topology. Here, we propose Kusaghiporia persicinus, comb. nov., based on a combination of morphological and molecular data. Laetiporus persicinus shares many morphological features with K. usambarensis that are missing in other Laetiporus species, including centrally stipitate basidiomata, a brown to pinkish pileus surface, and a pore layer that bruises when touched. However, K. usambarensis and L. persicinus differ in basidiospore size and shape as well as their geographic distributions. We provide a revised taxonomic treatment for this common wood-decay fungus.

Disease-related population declines in bats demonstrate non-exchangeability in generalist predators.

The extent to which persisting species may fill the functional role of extirpated or declining species has profound implications for the structure of biological communities and ecosystem functioning. In North America, arthropodivorous bats are threatened on a continent-wide scale by the spread of white-nose syndrome (WNS), a disease caused by the fungus Pseudogymnoascus destructans. We tested whether bat species that display lower mortality from this disease can partially fill the functional role of other bat species experiencing population declines. Specifically, we performed high-throughput amplicon sequencing of guano from two generalist predators: the little brown bat (Myotis lucifugus) and big brown bat (Eptesicus fuscus). We then compared changes in prey consumption before versus after population declines related to WNS. Dietary niches contracted for both species after large and abrupt declines in little brown bats and smaller declines in big brown bats, but interspecific dietary overlap did not change. Furthermore, the incidence and taxonomic richness of agricultural pest taxa detected in diet samples decreased following bat population declines. Our results suggest that persisting generalist predators do not necessarily expand their dietary niches following population declines in other predators, providing further evidence that the functional roles of different generalist predators are ecologically distinct.

The sexual spore pigment asperthecin is required for normal ascospore production and protection from UV light in Aspergillus nidulans.

Many fungi develop both asexual and sexual spores that serve as propagules for dissemination and/or recombination of genetic traits. Asexual spores are often heavily pigmented and this pigmentation provides protection from UV light. However, little is known about any purpose pigmentation that may serve for sexual spores. The model Ascomycete Aspergillus nidulans produces both green pigmented asexual spores (conidia) and red pigmented sexual spores (ascospores). Here we find that the previously characterized red pigment, asperthecin, is the A. nidulans ascospore pigment. The asperthecin biosynthetic gene cluster is composed of three genes: aptA, aptB, and aptC, where deletion of either aptA (encoding a polyketide synthase) or aptB (encoding a thioesterase) yields small, mishappen hyaline ascospores; while deletion of aptC (encoding a monooxygenase) yields morphologically normal but purple ascospores. ∆aptA and ∆aptB but not ∆aptC or wild type ascospores are extremely sensitive to UV light. We find that two historical ascospore color mutants, clA6 and clB1, possess mutations in aptA and aptB sequences, respectively.

Coarse Woody Debris Decomposition Assessment Tool: Model validation and application.

Coarse woody debris (CWD) is a significant component of the forest biomass pool; hence a model is warranted to predict CWD decomposition and its role in forest carbon (C) and nutrient cycling under varying management and climatic conditions. A process-based model, CWDDAT (Coarse Woody Debris Decomposition Assessment Tool) was calibrated and validated using data from the FACE (Free Air Carbon Dioxide Enrichment) Wood Decomposition Experiment utilizing pine (Pinus taeda), aspen (Populous tremuloides) and birch (Betula papyrifera) on nine Experimental Forests (EF) covering a range of climate, hydrology, and soil conditions across the continental USA. The model predictions were evaluated against measured FACE log mass loss over 6 years. Four widely applied metrics of model performance demonstrated that the CWDDAT model can accurately predict CWD decomposition. The R2 (squared Pearson's correlation coefficient) between the simulation and measurement was 0.80 for the model calibration and 0.82 for the model validation (P<0.01). The predicted mean mass loss from all logs was 5.4% lower than the measured mass loss and 1.4% lower than the calculated loss. The model was also used to assess the decomposition of mixed pine-hardwood CWD produced by Hurricane Hugo in 1989 on the Santee Experimental Forest in South Carolina, USA. The simulation reflected rapid CWD decomposition of the forest in this subtropical setting. The predicted dissolved organic carbon (DOC) derived from the CWD decomposition and incorporated into the mineral soil averaged 1.01 g C m-2 y-1 over the 30 years. The main agents for CWD mass loss were fungi (72.0%) and termites (24.5%), the remainder was attributed to a mix of other wood decomposers. These findings demonstrate the applicability of CWDDAT for large-scale assessments of CWD dynamics, and fine-scale considerations regarding the fate of CWD carbon.

Coarse woody debris decomposition assessment tool: Model development and sensitivity analysis.

Coarse woody debris (CWD) is an important component in forests, hosting a variety of organisms that have critical roles in nutrient cycling and carbon (C) storage. We developed a process-based model using literature, field observations, and expert knowledge to assess woody debris decomposition in forests and the movement of wood C into the soil and atmosphere. The sensitivity analysis was conducted against the primary ecological drivers (wood properties and ambient conditions) used as model inputs. The analysis used eighty-nine climate datasets from North America, from tropical (14.2° N) to boreal (65.0° N) zones, with large ranges in annual mean temperature (26.5°C in tropical to -11.8°C in boreal), annual precipitation (6,143 to 181 mm), annual snowfall (0 to 612 kg m-2), and altitude (3 to 2,824 m above mean see level). The sensitivity analysis showed that CWD decomposition was strongly affected by climate, geographical location and altitude, which together regulate the activity of both microbial and invertebrate wood-decomposers. CWD decomposition rate increased with increments in temperature and precipitation, but decreased with increases in latitude and altitude. CWD decomposition was also sensitive to wood size, density, position (standing vs downed), and tree species. The sensitivity analysis showed that fungi are the most important decomposers of woody debris, accounting for over 50% mass loss in nearly all climatic zones in North America. The model includes invertebrate decomposers, focusing mostly on termites, which can have an important role in CWD decomposition in tropical and some subtropical regions. The role of termites in woody debris decomposition varied widely, between 0 and 40%, from temperate areas to tropical regions. Woody debris decomposition rates simulated for eighty-nine locations in North America were within the published range of woody debris decomposition rates for regions in northern hemisphere from 1.6° N to 68.3° N and in Australia.

Predator preferences shape the diets of arthropodivorous bats more than quantitative local prey abundance.

Although most predators are generalists, the majority of studies on the association between prey availability and prey consumption have focused on specialist predators. To investigate the role of highly generalist predators in a complex food web, we measured the relationships between prey consumption and prey availability in two common arthropodivorous bats. Specifically, we used high-throughput amplicon sequencing coupled with a known mock community to characterize seasonal changes in little brown and big brown bat diets. We then linked spatiotemporal variation in prey consumption with quantitative prey availability estimated from intensive prey community sampling. We found that although quantitative prey availability fluctuated substantially over space and time, the most commonly consumed prey items were consistently detected in bat diets independently of their respective abundance. Positive relationships between prey abundance and probability of consumption were found only among prey groups that were less frequently detected in bat diets. While the probability of prey consumption was largely unrelated to abundance, the community structure of prey detected in bat diets was influenced by the local or regional abundance of prey. Observed patterns suggest that while little brown and big brown bats maintain preferences for particular prey independently of quantitative prey availability, total dietary composition may reflect some degree of opportunistic foraging. Overall, our findings suggest that generalist predators can display strong prey preferences that persist despite quantitative changes in prey availability.

Major histocompatibility complex variation is similar in little brown bats before and after white-nose syndrome outbreak.

White-nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans (Pd), has driven alarming declines in North American hibernating bats, such as little brown bat (Myotis lucifugus). During hibernation, infected little brown bats are able to initiate anti-Pd immune responses, indicating pathogen-mediated selection on the major histocompatibility complex (MHC) genes. However, such immune responses may not be protective as they interrupt torpor, elevate energy costs, and potentially lead to higher mortality rates. To assess whether WNS drives selection on MHC genes, we compared the MHC DRB gene in little brown bats pre- (Wisconsin) and post- (Michigan, New York, Vermont, and Pennsylvania) WNS (detection spanning 2014-2015). We genotyped 131 individuals and found 45 nucleotide alleles (27 amino acid alleles) indicating a maximum of 3 loci (1-5 alleles per individual). We observed high allelic admixture and a lack of genetic differentiation both among sampling sites and between pre- and post-WNS populations, indicating no signal of selection on MHC genes. However, post-WNS populations exhibited decreased allelic richness, reflecting effects from bottleneck and drift following rapid population declines. We propose that mechanisms other than adaptive immunity are more likely driving current persistence of little brown bats in affected regions.

A trait-based understanding of wood decomposition by fungi.

As the primary decomposers of organic material in terrestrial ecosystems, fungi are critical agents of the global carbon cycle. Yet our ability to link fungal community composition to ecosystem functioning is constrained by a limited understanding of the factors accounting for different wood decomposition rates among fungi. Here we examine which traits best explain fungal decomposition ability by combining detailed trait-based assays on 34 saprotrophic fungi from across North America in the laboratory with a 5-y field study comprising 1,582 fungi isolated from 74 decomposing logs. Fungal growth rate (hyphal extension rate) was the strongest single predictor of fungal-mediated wood decomposition rate under laboratory conditions, and accounted for up to 27% of the in situ variation in decomposition in the field. At the individual level, decomposition rate was negatively correlated with moisture niche width (an indicator of drought stress tolerance) and with the production of nutrient-mineralizing extracellular enzymes. Together, these results suggest that decomposition rates strongly align with a dominance-tolerance life-history trade-off that was previously identified in these isolates, forming a spectrum from slow-growing, stress-tolerant fungi that are poor decomposers to fast-growing, highly competitive fungi with fast decomposition rates. Our study illustrates how an understanding of fungal trait variation could improve our predictive ability of the early and midstages of wood decay, to which our findings are most applicable. By mapping our results onto the biogeographic distribution of the dominance-tolerance trade-off across North America, we approximate broad-scale patterns in intrinsic fungal-mediated wood decomposition rates.

Fungal functional ecology: bringing a trait-based approach to plant-associated fungi.

Fungi play many essential roles in ecosystems. They facilitate plant access to nutrients and water, serve as decay agents that cycle carbon and nutrients through the soil, water and atmosphere, and are major regulators of macro-organismal populations. Although technological advances are improving the detection and identification of fungi, there still exist key gaps in our ecological knowledge of this kingdom, especially related to function. Trait-based approaches have been instrumental in strengthening our understanding of plant functional ecology and, as such, provide excellent models for deepening our understanding of fungal functional ecology in ways that complement insights gained from traditional and -omics-based techniques. In this review, we synthesize current knowledge of fungal functional ecology, taxonomy and systematics and introduce a novel database of fungal functional traits (FunFun ). FunFun is built to interface with other databases to explore and predict how fungal functional diversity varies by taxonomy, guild, and other evolutionary or ecological grouping variables. To highlight how a quantitative trait-based approach can provide new insights, we describe multiple targeted examples and end by suggesting next steps in the rapidly growing field of fungal functional ecology.

Relationships among wood-boring beetles, fungi, and the decomposition of forest biomass.

A prevailing paradigm in forest ecology is that wood-boring beetles facilitate wood decay and carbon cycling, but empirical tests have yielded mixed results. We experimentally determined the effects of wood borers on fungal community assembly and wood decay within pine trunks in the southeastern United States. Pine trunks were made either beetle-accessible or inaccessible. Fungal communities were compared using culturing and high-throughput amplicon sequencing (HTAS) of DNA and RNA. Prior to beetle infestation, living pines had diverse fungal endophyte communities. Endophytes were displaced by beetle-associated fungi in beetle-accessible trees, whereas some endophytes persisted as saprotrophs in beetle-excluded trees. Beetles increased fungal diversity several fold. Over forty taxa of Ascomycota were significantly associated with beetles, but beetles were not consistently associated with any known wood-decaying fungi. Instead, increasing ambrosia beetle infestations caused reduced decay, consistent with previous in vitro experiments that showed beetle-associated fungi reduce decay rates by competing with decay fungi. No effect of bark-inhabiting beetles on decay was detected. Platypodines carried significantly more fungal taxa than scolytines. Molecular results were validated by synthetic and biological mock communities and were consistent across methodologies. RNA sequencing confirmed that beetle-associated fungi were biologically active in the wood. Metabarcode sequencing of the LSU/28S marker recovered important fungal symbionts that were missed by ITS2, though community-level effects were similar between markers. In contrast to the current paradigm, our results indicate ambrosia beetles introduce diverse fungal communities that do not extensively decay wood, but instead reduce decay rates by competing with wood decay fungi.

An improved method for utilizing high-throughput amplicon sequencing to determine the diets of insectivorous animals.

DNA analysis of predator faeces using high-throughput amplicon sequencing (HTS) enhances our understanding of predator-prey interactions. However, conclusions drawn from this technique are constrained by biases that occur in multiple steps of the HTS workflow. To better characterize insectivorous animal diets, we used DNA from a diverse set of arthropods to assess PCR biases of commonly used and novel primer pairs for the mitochondrial gene, cytochrome oxidase C subunit 1 (COI). We compared diversity recovered from HTS of bat guano samples using a commonly used primer pair "ZBJ" to results using the novel primer pair "ANML." To parameterize our bioinformatics pipeline, we created an arthropod mock community consisting of single-copy (cloned) COI sequences. To examine biases associated with both PCR and HTS, mock community members were combined in equimolar amounts both pre- and post-PCR. We validated our system using guano from bats fed known diets and using composite samples of morphologically identified insects collected in pitfall traps. In PCR tests, the ANML primer pair amplified 58 of 59 arthropod taxa (98%), whereas ZBJ amplified 24-40 of 59 taxa (41%-68%). Furthermore, in an HTS comparison of field-collected samples, the ANML primers detected nearly fourfold more arthropod taxa than the ZBJ primers. The additional arthropods detected include medically and economically relevant insect groups such as mosquitoes. Results revealed biases at both the PCR and sequencing levels, demonstrating the pitfalls associated with using HTS read numbers as proxies for abundance. The use of an arthropod mock community allowed for improved bioinformatics pipeline parameterization.

Characterization of PdCP1, a serine carboxypeptidase from Pseudogymnoascus destructans, the causal agent of White-nose Syndrome.

Pseudogymnoascus destructans is a pathogenic fungus responsible for White-nose Syndrome (WNS), a disease afflicting multiple species of North American bats. Pseudogymnoascus destructans infects susceptible bats during hibernation, invading dermal tissue and causing extensive tissue damage. In contrast, other Pseudogymnoascus species are non-pathogenic and cross-species comparisons may therefore reveal factors that contribute to virulence. In this study, we compared the secretome of P. destructans with that from several closely related Pseudogymnoascus species. A diverse set of hydrolytic enzymes were identified, including a putative serine peptidase, PdCP1, that was unique to the P. destructans secretome. A recombinant form of PdCP1 was purified and substrate preference determined using a multiplexed-substrate profiling method based on enzymatic degradation of a synthetic peptide library and analysis by mass spectrometry. Most peptide substrates were sequentially truncated from the carboxyl-terminus revealing that this enzyme is a bona fide carboxypeptidase. Peptides with arginine located close to the carboxyl-terminus were rapidly cleaved, and a fluorescent substrate containing arginine was therefore used to characterize PdCP1 activity and to screen a selection of peptidase inhibitors. Antipain and leupeptin were found to be the most potent inhibitors of PdCP1 activity.

Non-biological synthetic spike-in controls and the AMPtk software pipeline improve mycobiome data.

High-throughput amplicon sequencing (HTAS) of conserved DNA regions is a powerful technique to characterize microbial communities. Recently, spike-in mock communities have been used to measure accuracy of sequencing platforms and data analysis pipelines. To assess the ability of sequencing platforms and data processing pipelines using fungal internal transcribed spacer (ITS) amplicons, we created two ITS spike-in control mock communities composed of cloned DNA in plasmids: a biological mock community, consisting of ITS sequences from fungal taxa, and a synthetic mock community (SynMock), consisting of non-biological ITS-like sequences. Using these spike-in controls we show that: (1) a non-biological synthetic control (e.g., SynMock) is the best solution for parameterizing bioinformatics pipelines, (2) pre-clustering steps for variable length amplicons are critically important, (3) a major source of bias is attributed to the initial polymerase chain reaction (PCR) and thus HTAS read abundances are typically not representative of starting values. We developed AMPtk, a versatile software solution equipped to deal with variable length amplicons and quality filter HTAS data based on spike-in controls. While we describe herein a non-biological SynMock community for ITS sequences, the concept and AMPtk software can be widely applied to any HTAS dataset to improve data quality.

Draft Genome Sequence of Burkholderia cepacia ATCC 17759, a Polyhydroxybutyrate-Co-Valerate Copolymer-Producing Bacterium.

Burkholderia cepacia ATCC 17759, isolated from forest soils in Trinidad, accumulates large amounts of polyhydroxyalkanoate copolymers when grown on xylose, mannose, arabinose, other carbohydrates, and organic acid cosubstrates. This 8.72-Mb draft genome sequence of B. cepacia ATCC 17759 will provide better insight into this organism's utility in lignocellulose bioconversion.

Detecting Symbioses in Complex Communities: the Fungal Symbionts of Bark and Ambrosia Beetles Within Asian Pines.

Separating symbioses from incidental associations is a major obstacle in symbiosis research. In this survey of fungi associated with Asian bark and ambrosia beetles, we used quantitative culture and DNA barcode identification to characterize fungal communities associated with co-infesting beetle species in pines (Pinus) of China and Vietnam. To quantitatively discern likely symbioses from coincidental associations, we used multivariate analysis and multilevel pattern analysis (a type of indicator species analysis). Nearly half of the variation in fungal community composition in beetle galleries and on beetle bodies was explained by beetle species. We inferred a spectrum of ecological strategies among beetle-associated fungi: from generalist multispecies associates to highly specialized single-host symbionts that were consistently dominant within the mycangia of their hosts. Statistically significant fungal associates of ambrosia beetles were typically only found with one beetle species. In contrast, bark beetle-associated fungi were often associated with multiple beetle species. Ambrosia beetles and their galleries were frequently colonized by low-prevalence ambrosia fungi, suggesting that facultative ambrosial associations are commonplace, and ecological mechanisms such as specialization and competition may be important in these dynamic associations. The approach used here could effectively delimit symbiotic interactions in any system where symbioses are obscured by frequent incidental associations. It has multiple advantages including (1) powerful statistical tests for non-random associations among potential symbionts, (2) simultaneous evaluation of multiple co-occurring host and symbiont associations, and (3) identifying symbionts that are significantly associated with multiple host species.

Extreme sensitivity to ultraviolet light in the fungal pathogen causing white-nose syndrome of bats.

Bat white-nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans, has decimated North American hibernating bats since its emergence in 2006. Here, we utilize comparative genomics to examine the evolutionary history of this pathogen in comparison to six closely related nonpathogenic species. P. destructans displays a large reduction in carbohydrate-utilizing enzymes (CAZymes) and in the predicted secretome (~50%), and an increase in lineage-specific genes. The pathogen has lost a key enzyme, UVE1, in the alternate excision repair (AER) pathway, which is known to contribute to repair of DNA lesions induced by ultraviolet (UV) light. Consistent with a nonfunctional AER pathway, P. destructans is extremely sensitive to UV light, as well as the DNA alkylating agent methyl methanesulfonate (MMS). The differential susceptibility of P. destructans to UV light in comparison to other hibernacula-inhabiting fungi represents a potential "Achilles' heel" of P. destructans that might be exploited for treatment of bats with WNS.

Phylogenetics of a Fungal Invasion: Origins and Widespread Dispersal of White-Nose Syndrome.

Globalization has facilitated the worldwide movement and introduction of pathogens, but epizoological reconstructions of these invasions are often hindered by limited sampling and insufficient genetic resolution among isolates. Pseudogymnoascus destructans, a fungal pathogen causing the epizootic of white-nose syndrome in North American bats, has exhibited few genetic polymorphisms in previous studies, presenting challenges for both epizoological tracking of the spread of this fungus and for determining its evolutionary history. We used single nucleotide polymorphisms (SNPs) from whole-genome sequencing and microsatellites to construct high-resolution phylogenies of P. destructans Shallow genetic diversity and the lack of geographic structuring among North American isolates support a recent introduction followed by expansion via clonal reproduction across the epizootic zone. Moreover, the genetic relationships of isolates within North America suggest widespread mixing and long-distance movement of the fungus. Genetic diversity among isolates of P. destructans from Europe was substantially higher than in those from North America. However, genetic distance between the North American isolates and any given European isolate was similar to the distance between the individual European isolates. In contrast, the isolates we examined from Asia were highly divergent from both European and North American isolates. Although the definitive source for introduction of the North American population has not been conclusively identified, our data support the origin of the North American invasion by P. destructans from Europe rather than Asia.IMPORTANCE This phylogenetic study of the bat white-nose syndrome agent, P. destructans, uses genomics to elucidate evolutionary relationships among populations of the fungal pathogen to understand the epizoology of this biological invasion. We analyze hypervariable and abundant genetic characters (microsatellites and genomic SNPs, respectively) to reveal previously uncharacterized diversity among populations of the pathogen from North America and Eurasia. We present new evidence supporting recent introduction of the fungus to North America from a diverse Eurasian population, with limited increase in genetic variation in North America since that introduction.

Amplicon-Based Sequencing of Soil Fungi from Wood Preservative Test Sites.

Soil samples were collected from field sites in two AWPA (American Wood Protection Association) wood decay hazard zones in North America. Two field plots at each site were exposed to differing preservative chemistries via in-ground installations of treated wood stakes for approximately 50 years. The purpose of this study is to characterize soil fungal species and to determine if long term exposure to various wood preservatives impacts soil fungal community composition. Soil fungal communities were compared using amplicon-based DNA sequencing of the internal transcribed spacer 1 (ITS1) region of the rDNA array. Data show that soil fungal community composition differs significantly between the two sites and that long-term exposure to different preservative chemistries is correlated with different species composition of soil fungi. However, chemical analyses using ICP-OES found levels of select residual preservative actives (copper, chromium and arsenic) to be similar to naturally occurring levels in unexposed areas. A list of indicator species was compiled for each treatment-site combination; functional guild analyses indicate that long-term exposure to wood preservatives may have both detrimental and stimulatory effects on soil fungal species composition. Fungi with demonstrated capacity to degrade industrial pollutants were found to be highly correlated with areas that experienced long-term exposure to preservative testing.