Diversity and Structure of Soil Fungal Communities across Experimental Everglades Tree Islands

Fungi play prominent roles in ecosystem services (e.g., nutrient cycling, decomposition) and thus have increasingly garnered attention in restoration ecology. However, it is unclear how most management decisions impact fungal communities, making it difficult to protect fungal diversity and utilize fungi to improve restoration success. To understand the effects of restoration decisions and environmental variation on fungal communities, we sequenced soil fungal microbiomes from 96 sites across eight experimental Everglades tree islands approximately 15 years after restoration occurred. We found that early restoration decisions can have enduring consequences for fungal communities. Factors experimentally manipulated in 2003–2007 (e.g., type of island core) had significant legacy effects on fungal community composition. Our results also emphasized the role of water regime in fungal diversity, composition, and function. As the relative water level decreased, so did fungal diversity, with an approximately 25% decline in the driest sites. Further, as the water level decreased, the abundance of the plant pathogen–saprotroph guild increased, suggesting that low water may increase plant-pathogen interactions. Our results indicate that early restoration decisions can have long-term consequences for fungal community composition and function and suggest that a drier future in the Everglades could reduce fungal diversity on imperiled tree islands.


Introduction
Fungi play important roles in many ecosystem functions and services, especially those that involve soil [1,2], where they make up an estimated 55-85% of the microbial biomass [3,4]. These fungi are crucial for the decomposition of organic carbon, cycling of nitrogen and phosphorus, and belowground carbon sequestration [1,[5][6][7]. Soil fungi also indirectly contribute to ecosystem function through their interactions with primary producers. For instance, they affect plant growth and community composition through pathogenic attacks on particular plant taxa, changes to plant-plant competition, and beneficial interactions that ameliorate environmental stress [8][9][10][11][12]. Given the important ecological roles of fungi and their interactions with primary producers, soil fungal communities can be a valuable tool for improving terrestrial restoration but are often overlooked in both restoration planning and the assessment of success [13]. Research has shown that fungal amendment can improve the fertility increasing tree resilience to natural and anthropogenic stress [9], and fungal pathogens may be crucial for restoring a natural tree community composition by regulating the abundance of dominant plant species [37]. While fungal communities in the peat and periphyton of the Everglades wetlands have been investigated in a few cases [38][39][40], the terrestrial fungal communities of Everglades tree islands are unexplored. In fact, we are aware of no studies of the diversity or composition of soil fungal communities on Everglades tree islands. To improve the protection, and even utilization, of fungal diversity in the restoration of tree islands, studies of the factors that structure these fungal communities are needed. Here, we use a set of eight experimental Everglades tree islands, in which abiotic and biotic factors have been manipulated to understand tree island formation, to better identify how microhabitat environmental variation and experimental restoration decisions may influence fungal diversity, composition, and functional guilds.

Study Site and Environmental Data Collection
We studied tree island soil fungal communities at the approximately 32-ha Loxahatchee Impoundment Landscape Assessment (LILA) facility in Loxahatchee National Wildlife Refuge in Palm Beach County, Florida, USA (26.489 • N, 80.219 • W). LILA contains eight approximately 2500 m 2 experimental tree islands that were constructed in 2003 through a collaboration between the South Florida Water Management District and the Army Corps of Engineers. LILA is especially suitable for investigating how the soil microbiome responds to different restoration decisions, because the construction of tree islands allowed landscape-level experimental manipulations of tree island characteristics. For instance, the islands' cores (i.e., bases of the islands) were manipulated to represent two common island types in the Everglades: peat core and limestone core islands [41]. Using locally sourced peat and limestone derived from the habitat adjacent to the current LILA facility, half of the experimental islands were constructed with peat cores and half were constructed with limestone cores ( Figure S1); then, all the islands were covered with a top layer of peat substrate [42]. Similarly, each of the eight experimental islands was split into quadrants in which a mixture of 10 tree species were planted at four different densities (1 m, 1.67 m, 2.33 m, and 3 m spacing) in 2006 and 2007 ( Figures S1 and S2). Given the importance of plant-fungal interactions, tree planting density may be a biotic factor that strongly affects fungal community composition. In addition, LILA has other readily available data on microhabitat variation within and among islands that can help inform our understanding of the soil microbiome. For example, because the water level is continuously monitored within the macrocosms and the topography of each island is well mapped, it is possible to track the hydrology of microhabitats using the DBHydro project from the South Florida Water Management District [42]. Using surface water-level data and the elevation of each site, we calculated the average site-specific 'relative water level' for all plots in 2018 (the year fungi were sampled) following the methods in [43]. Further, understory plant communities that recruited to these experimental tree islands have been monitored since 2009 [42], allowing for consideration of the understory plant community in structuring soil microbiomes. For these understory plant communities, we calculated Shannon diversity, Pielou's evenness, and plant richness for use in our analyses. During the monitoring of each site, understory plant biomass was estimated (by applying allometric equations developed from separate biomass collection plots to estimates of plant cover). To gain insight into the light environment (another abiotic feature of interest), hemispherical canopy photographs were also taken using a digital camera (Nikon Coolpix 995; Nikon, Tokyo, Japan) equipped with a hemispherical lens (Nikon Fisheye Converter FC-E8 0.21×, Tokyo, Japan) at each site, and canopy openness was determined using the software Gap Light Analyzer (GLA), version 2.0, Burnaby, BC, Canada [44].

Soil Collection
In June 2018, we collected soil microbiomes from 96 sites. We colocated our fungal community sampling with 96 understory plant community plots established in 2009. These plots were evenly distributed among islands (12 per island) and among experimental treatments, i.e., island core type treatments (48 plots on peat core versus 48 plots on limestone core island) and tree planting density treatments (24 plots per density). On each island, the three plots were randomly assigned to locations within each tree density quadrant (e.g., Figure S2). For each sampling point, we used a soil corer to aseptically collect approximately 50 mL of soil at 15-25 cm below the soil surface and approximately 0.25 m from the center of the 1 m 2 plot (indicated by permanent PVC site marker). Soil samples were placed in sterile conical tubes and transported on ice to University of Miami's campus where they were kept at −80 • C.

DNA Sequencing and Data Processing
We extracted DNA from 0.25 g of each soil sample using the E.Z.N.A. Soil DNA kit (Omega, D5625-01, Norcross, GA, USA). We removed PCR inhibitors from DNA samples with agarose gel electrophoresis and extracted genomic bands (i.e., >15 kilobases; E.Z.N.A Gel Extraction kit, Omega, D2500-01, Norcross, GA, USA). We performed quality control by endpoint PCR of the fungal ribosomal ITS1 region to confirm that enough PCR inhibitors were removed to obtain DNA of sufficient quality for amplification. The ITS1 region was amplified via PCR using two-step dual-indexing [45], and amplicon libraries were sequenced at the University of Minnesota Genomics Center (UMGC) using the Illumina MiSeq platform (v3, 300 bp paired-end, San Diego, CA, USA). We also sequenced three negative controls in which ultrapure water (G-Biosciences, St. Louis, MO, USA; 786-293) was used in place of soil during extractions to confirm that samples were not contaminated during the DNA extraction process. UMGC demultiplexed reads using bcl2fastq. We denoised reads, joined paired ends, and grouped reads into exact sequence variants (ESVs; i.e., operational taxonomic units at 100% sequence similarity) with the QIIME2 pipeline (v2019.1, [46]). We rarefied sequences to 4000 reads per sample to account for unequal sequencing depths and excluded samples with less than 4000 reads (12 of 96 samples). We classified ESVs into fungal "species" if they shared at least 97% sequence similarity to reference sequences in the UNITE database (version v7_01.12.2017, [47]). For the fungal community at each site, richness was determined based on the number of rarefied ESVs, and Shannon diversity was calculated using sequence reads of ESVs as estimates of abundance. Note that using ESVs, compared to OTUs (Operational Taxonomic Units) with 97% sequence similarity, may result in higher absolute fungal richness, but ESVs and 97% OTUs are often highly correlated, such that comparisons among treatments are robust to sequence binning choices [48].

Data Analyses
To understand which abiotic and biotic factors explain fungal community diversity and richness, we performed model selection with the abiotic and biotic factors in Table 1 as explanatory variables and identified the best model based on AIC using dredge (R package MuMin; [49]). We also included tree island identity as a random effect in the model selection to account for understory plant community variation between islands. We confirmed that model variables were not highly correlated with each other before performing model selection (all correlation coefficient <0.7; Figure 2). To investigate the possible causes of variation in fungal community composition, we performed a distance-based redundancy analysis (dbRDA, Bray-Curtis) with the same abiotic and biotic factors as explanatory variables using capscale (R package vegan; [50]). Indicator taxa analysis was used to identify which of the 500 most common fungal taxa were associated with differences in community composition between treatments for significant categorical variables (e.g., limestone versus peat core). We calculated 'indicator values' of the fungal taxa based on the relative abundance and consistency in each treatment and determined the significance of these values with permutation tests (R package indicspecies [51,52]), which were FDR (False Discovery Rate) corrected for multiple comparisons. We also conducted a Mantel test between fungal and plant community distance matrices using mantel (R package vegan, [50]) to examine whether the variation in fungal community composition can be explained by variation in understory plant community composition. To investigate the distribution of functional guilds of fungi, we used the FUNGuild database (version v1.1; [53]), which assigned each taxa identified at the species level via UNITE to a functional guild (e.g., leaf saprotroph or fungal parasite) and a trophic mode (e.g., symbiotrophic, pathotrophic) ( Figure S3A,B). ESVs that were not identified to the species level were not included in this analysis. We also filtered to include only taxa that had guild and trophic mode assignments with the confidence rankings "Highly Probable" and identified the three most common guilds based on their normalized ESV counts. For our analyses, we focused on functional guild assignments using model selection to determine which biotic and abiotic variables explained variation in the number of different guilds present at a site (guild richness) and the relative abundances of the top three most common functional guilds. As above, our model selection used linear mixed effects models with abiotic and biotic factors as explanatory variables and tree island identity as a random effect. All statistical analyses were performed in R ver 3.6.0.

Sequencing and ESV Taxonomic Assignments Statistics
After quality control filtering, 84 of the 96 samples had sufficient read counts to reach our saturation cutoff (≥4000 reads) required to be included in analyses. Across these samples, there were over 9000 ESVs, with an average of 243 ± 11 ESVs per sample. Using the UNITE database, we identified approximately 35% of all ESVs and an average of 100 ± 5 taxa per sample (mean ± s.e) to at least the order level. Approximately 20% of all ESVs were identified in UNITE to the species level, allowing them to be included in our functional analyses using FUNGuild. Fusarium keratoplasticum was most common across sites, occurring in 71 of the 84 sites, followed by Penicillium sp. (69 of 84 sites) and Retroconis fusiformis (69 of 84 sites). An unidentified fungal taxa had the greatest number of reads on average (1149) followed in abundance by an ESV from the order Tremellales (614) and Inocybe curvipes (597). Demultiplexed sequence data, the ESV table, and the taxonomy table have been submitted to NCBI (BioProject: PRJNA639837).

Fungal Community Diversity and Richness
The relative water level and evenness of the understory plant community were identified by model selection as key environmental characteristics explaining variation in fungal diversity (Table 1). Relative water level and understory plant community evenness as well as canopy openness, understory plant richness, core type, and tree density were also included in the best model for fungal richness (Table 1). Further investigation of the significant explanatory variables in the best fit models revealed a positive relationship with several environmental factors for both diversity and richness. For example, fungal diversity increased with increasing relative water level (F 1,68 = 8.63, p = 0.004; Figure 1A) and with increasing understory plant community evenness (F 1,68 = 8.83, p = 0.004; Figure 1B), suggesting the maintenance of higher water levels and even plant understories on tree islands could improve fungal diversity and therefore the variety of ecosystem services they provide (as supported by the functional guild analysis below). Fungal community richness also increased with increasing canopy openness (F 1,65 = 9.41, p = 0.003; Figure 1C) as well as with increasing understory evenness (F 1,65 = 7.66, p = 0.007; Figure 1D).
Diversity 2020, 12, x FOR PEER REVIEW 6 of 17 supported by the functional guild analysis below). Fungal community richness also increased with increasing canopy openness (F1,65= 9.41, p = 0.003; Figure 1C) as well as with increasing understory evenness (F1,65 = 7.66, p = 0.007; Figure 1D). In all panels, each point represents a site in the experimental tree islands, and lines are based on a linear model relationship between these variables with the shaded area indicating the 95% confidence intervals around the line. In (A), note that relative water level values equal to zero represent sites where the water level was on average at the soil surface. Values greater than zero represent the sites that were on average inundated, and values less than zero represent sites that on average had water below the soil surface.

Fungal Community Composition
In addition to its importance in fungal diversity, relative water level also explained a significant amount of variation in fungal community composition (db-RDA, F1,68 = 2.95, p = 0.001; Figure 2). The core type, one of the manipulative restoration treatments implemented during the construction of the experimental tree islands, affected fungal community composition as well (F1,68 = 1.87, p = 0.002). We found that limestone cores resulted in a significantly different community composition than peat In all panels, each point represents a site in the experimental tree islands, and lines are based on a linear model relationship between these variables with the shaded area indicating the 95% confidence intervals around the line. In (A), note that relative water level values equal to zero represent sites where the water level was on average at the soil surface. Values greater than zero represent the sites that were on average inundated, and values less than zero represent sites that on average had water below the soil surface.

Fungal Community Composition
In addition to its importance in fungal diversity, relative water level also explained a significant amount of variation in fungal community composition (db-RDA, F 1,68 = 2.95, p = 0.001; Figure 2). The core type, one of the manipulative restoration treatments implemented during the construction of the experimental tree islands, affected fungal community composition as well (F 1,68 = 1.87, p = 0.002). We found that limestone cores resulted in a significantly different community composition than peat cores on both of the first two axes of fungal community composition (CAPs 1 and 2; Figure 2), with a 218% higher value on CAP 1 (F 1,76 = 3.83, p = 0.054) and a 185% smaller value on CAP2 (F 1,76 = 61.55, p = 2.17 × 10 −11 ) for limestone compared to peat core islands. To gain insight into the difference between the communities in peat and limestone core islands, we performed an indicator taxa analysis. After correcting for multiple comparisons, we determined that Thelephoraceae, Eurotiomycetes, and Inocybe curvipes were indicative of peat communities, while Agaricomycetes and Archaeorhizomyces were indicative of limestone communities (Table S1). We also found that there was a significant relationship between the inter-site variations in the fungal and understory plant community composition (Mantel test: R = 0.11, p = 0.003).
Diversity 2020, 12, x FOR PEER REVIEW 7 of 17 cores on both of the first two axes of fungal community composition (CAPs 1 and 2; Figure 2), with a 218% higher value on CAP 1 (F1,76 = 3.83, p = 0.054) and a 185% smaller value on CAP2 (F1,76 = 61.55, p = 2.17 × 10 −11 ) for limestone compared to peat core islands. To gain insight into the difference between the communities in peat and limestone core islands, we performed an indicator taxa analysis. After correcting for multiple comparisons, we determined that Thelephoraceae, Eurotiomycetes, and Inocybe curvipes were indicative of peat communities, while Agaricomycetes and Archaeorhizomyces were indicative of limestone communities (Table S1). We also found that there was a significant relationship between the inter-site variations in the fungal and understory plant community composition (Mantel test: R = 0.11, p = 0.003).

Distribution of Fungal Functional Guilds
In total, we found 15 fungal functional guilds in our dataset with each of the eight islands hosting an average of approximately 5 guilds (±0.22) ( Figure S3A). The most common functional guilds were ectomycorrhizal fungi, dung-wood saprotroph, and plant pathogen-wood saprotroph, which Ordination demonstrating the relationship between fungal community composition and environmental variables. Each point represents the fungal community composition at a site on the experimental tree islands. Points are colored by relative water level (in meters) from high water (light colors) to low water (dark colors). The graphs on either side of the ordination indicate the mean and standard error of tree density and island core type along the first two axes of community composition (CAP1 and CAP2). Different lowercase letters denote significant differences. The table (in upper right corner) details correlations between continuous environmental variables.

Distribution of Fungal Functional Guilds
In total, we found 15 fungal functional guilds in our dataset with each of the eight islands hosting an average of approximately 5 guilds (±0.22) ( Figure S3A). The most common functional guilds were ectomycorrhizal fungi, dung-wood saprotroph, and plant pathogen-wood saprotroph, which occurred in 63%, 77%, and 67% of the sites (respectively) and on all islands ( Figure S3A). When examining factors that influence guild richness, we found that the number of functional guilds among fungi increased in sites that had a more even understory plant community (F 1,63 = 10.12, p = 0.0023; Figure 3A) with an approximately 70% increase in guilds in the most even compared to the least even five sites. Core type and relative water level were also included in the best model for the number of functional guilds (core type: F 1,63 = 2.92, p = 0.1378; relative water level: F 1,63 = 0.87, p = 0.3521), with an 18% greater number of guilds in peat compared to limestone core islands and an approximately 8% decrease in the driest compared to wettest sites. When we examined the factors that influence abundances of the three most common guilds from our dataset, the abundance of only one-plant pathogen-wood saprotroph-could be explained by any of our measured variables. As the relative water level decreased, the abundance of the plant pathogen-wood saprotroph guild increased (F 1,68 = 5.47, p = 0.022; Figure 3B). fungi increased in sites that had a more even understory plant community (F1,63 = 10.12, p = 0.0023; Figure 3A) with an approximately 70% increase in guilds in the most even compared to the least even five sites. Core type and relative water level were also included in the best model for the number of functional guilds (core type: F1,63 = 2.92, p = 0.1378; relative water level: F1,63 = 0.87, p = 0.3521), with an 18% greater number of guilds in peat compared to limestone core islands and an approximately 8% decrease in the driest compared to wettest sites. When we examined the factors that influence abundances of the three most common guilds from our dataset, the abundance of only one-plant pathogen-wood saprotroph-could be explained by any of our measured variables. As the relative water level decreased, the abundance of the plant pathogen-wood saprotroph guild increased (F1,68 = 5.47, p = 0.022; Figure 3B).

Discussion
To better incorporate fungal communities into restoration planning in terms of both protecting their diversity and utilizing these hidden players to improve restoration success, a more complete understanding of which environmental factors and management decisions affect soil fungal

Discussion
To better incorporate fungal communities into restoration planning in terms of both protecting their diversity and utilizing these hidden players to improve restoration success, a more complete understanding of which environmental factors and management decisions affect soil fungal communities is required. Our study contributes to this goal by showing that (1) fungal diversity and composition were affected by restoration decision factors that were manipulated on experimental tree islands during their construction, and (2) variation in several important aspects of fungal communities was explained by microhabitat differences in other environmental variables of interest for management, including hydrology and properties of the naturally recruiting plant understory community. Below, we address these results in more detail by first examining possible mechanisms through which abiotic factors, and then biotic factors, can contribute to fungal community composition and diversity. We conclude by briefly discussing directions for future research and the implications for the Everglades.

Abiotic Factors: The Role of Hydrology, Island Core Type, and Light Environment in Fungal Diversity and Composition
Abiotic stressors play important roles in limiting organismal niches [54], and thus abiotic contexts often contribute substantially to the diversity and composition of communities in nature [55]. Previous work in the Everglades has identified abiotic features that affect communities of macro-organisms on tree islands [56,57]. For instance, the tree island plant community structure is largely driven by water level, nutrient availability, and disturbance (e.g., fire, windstorms, and drought) [30]. Our study found that the fungal community composition and/or diversity on tree islands also appear to respond strongly to all three of the abiotic factors investigated-relative water level, light environment, and core type. In particular, hydrology was important for the composition and diversity of the soil fungal community, with greater fungal diversity at sites with higher water levels. We saw a 34% increase in fungal diversity when comparing the five wettest sites to the five driest sites. This shift in fungal diversity and composition is likely due to the sensitivity of some species of fungi to dry conditions [58,59] or the increased dispersal of spores to wetter sites [38,60]. Hydrology is one of the most anthropogenically altered aspects of the Everglades, with much of the historic sheet flow from central to southern Florida now diverted through approximately 1800 miles of levees and canals [61]. This change has generally led to a change of inundation periods for tree islands [62]. Given the substantial decrease in fungal diversity we found at dry sites in our study, it is likely that these major alterations to tree island inundation are also having extensive consequences for the fungal communities affected by water diversions. As a result of the importance of fungi in ecosystem functions [1], a loss of fungal diversity on unnaturally dry tree islands may underpin mechanisms causing tree island loss (e.g., reduced availability of the beneficial symbionts on which trees depend, or a loss of high-quality decomposers required for nutrient recycling). These shifts in diversity and composition also suggest the importance of hydrologic restoration targets that consider fungal ecology, and tree island restoration and creation success may be assisted by a more careful consideration of these important constituents of soil biota.
In addition, our investigation of abiotic factors demonstrated that early restoration decisions can have abiding effects on fungal communities years later. In 2003, the experimental islands at LILA were constructed with two core types-peat and limestone-to represent the natural variation in Everglades tree islands. We found that differences in the choice of core type affected the overall fungal community composition, including a significant difference between limestone core and peat core communities along both axes of variation in community composition. Further research into the indicator taxa identified for limestone versus peat core islands revealed functional differences among these taxa. Indicator taxa of peat communities included ectomycorrhizal fungi [63,64], while those in limestone communities included taxa described as putative wood-decayers and saprotrophs [65][66][67]. The distinctions in fungal communities on the two core types are presumably the result of the different soil environments they provide [68]. For instance, previous studies of the tree islands at LILA have suggested increased water retention on peat islands compared to limestone [69,70]. The presence of limestone may also change the mineral content of soils through inputs of calcium carbonate, which has been associated with significant shifts in soil microbial communities in other systems [71]. In addition, the difference in some of these indicator taxa could also be the result of the presence of the tree Morella cerifera, which is the only surveyed species on the islands that has been reported to associate with ectomycorrhizal fungi (in addition to its association with arbuscular mycorrhizal fungi) [72]. Morella cerifera had a somewhat lower survival rate on peat core islands, indicating they are more stressed on these islands, which may result in an increased reliance on ectomycorrhizal fungi to ameliorate stress [42]. Therefore, we examined whether the presence of M. cerifera was predictive of the estimated abundance (total number of reads) of the ectomycorrhizal indicator taxa. While we did not find support for a clear relationship between M. cerifera presence and these ectomycorrhizal fungi (p > 0.05 in a distribution-free randomization test with the presence of M. cerifera and island identity in the model), additional studies that investigate unmeasured aspects of M. cerifera biology, such as tree size, and the root colonization of this tree may help elucidate if and how it influences ectomycorrhizal abundance. More generally, the effect of core type in our study indicates that initial decisions in tree island construction have cascading effects on the fungal community, even 15 years later.
We also found that as light availability increased (as indicated by canopy openness), fungal richness also increased. While light environments can have consequences for the soil microbiome through effects on soil temperature and moisture [73], light may act indirectly through its effects on the understory plant community [74]. Further investigation of this possibility using an RDA (Redundancy Analysis) with plant understory community composition data as the response and the same explanatory abiotic and biotic factors used in the fungal analyses showed that canopy openness did not explain variation in understory composition (F 1,68 = 1.24, p = 0.251). While not conclusive, this result suggests that the light environment's effects on the soil microbiome is unlikely to be driven by changes in the plant understory composition.

Biotic Factors: The Role of Understory Plant Communities and Tree Density in Fungal Community Diversity and Composition
Our study showed that biotic factors may also have consequences for the Everglades tree islands soil fungi. We found that the characteristics of the understory plant community explained variation in fungal community diversity, richness, and composition. For instance, as understory plant community evenness increased, so did both fungal diversity and richness. This relationship may stem from differences among plant species in their priming of soil microbial communities [75] and/or their contributions to the leaf litter [76]. Plants can actively manipulate the soil microbial community through the release of root exudates [77] and allelochemicals [78], and this priming can be species-specific. Similarly, plants can have species-specific differences in their leaf litter, and the composition of leaf litter can influence fungal communities [79]. Soil priming and litter deposition by a more even understory plant community may promote a more even distribution of soil microhabitats for fungal taxa to use, thereby fostering the greater fungal diversity we found. We also demonstrated that differences in plant community composition among sites are significantly related to differences in fungal community composition. Unlike the effects of manipulated variables such as core type and tree density, we cannot determine the direction of interaction (i.e., plant understory community composition is driving fungal community composition, the reverse, or both are responding to an external pressure). However, it is well known that fungal and plant communities are interlinked [80]. Understory plant communities may change fungal communities through their root exudates [77], leaf litter [81,82], or by acting as hosts for pathogenic or mutualistic fungi [83]. Furthermore, a thriving understory community may provide shade, thus changing the microhabitat conditions of soil moisture. On the other hand, fungal composition may alter plant composition through nutrient availability [84], decomposition [85], water availability [9], pathogen load [86], or mutualistic interactions [87]. Most likely, both the fungal and plant communities are affecting one another via plant-soil feedbacks [88] and are also mutually responding to major abiotic factors such as hydrology. Future work that manipulates these components (i.e., abiotic factors, microbial and plant community characteristics) using factorial mesocosm experiments is needed to fully determine the strength and direction in which these communities are affecting one another, and how much outside habitat characteristics contribute to this relationship. In addition, the initial tree planting density was identified by model selection as part of the best model for fungal community richness. While its effect was not significant within the best model (nor was tree density as commonly implicated as the understory plant community in our analyses), tree density's inclusion in the best model hints at a role for this experimentally manipulated factor in influencing fungal communities. Future studies investigating the effects of tree composition on fungal communities could provide additional insights into consequences for fungi that investigating tree density alone may miss.

Fungal Functional Guilds on Tree Islands
In addition to characterizing fungal diversity and composition, our study provided initial insight into fungal contributions to ecosystem services by characterizing the distribution of fungal functional guilds on tree islands. While not all fungal taxa could be characterized into functional groups by FUNGuild and functions sometimes differ between ecosystems [89], this approach represents a first step in understanding how different restoration decisions and environmental variation may impact fungal roles within Everglades tree islands. For example, we found that the number of guilds was significantly greater in sites with a more even plant understory community, which is in line with the greater fungal taxa diversity at these sites. The arbuscular mycorrhizal fungal guild (fungal mutualists that associate with plant roots and trade water/nutrients acquired through hyphal networks for photosynthetic carbon) was present at all five of the sites with the highest plant understory evenness, while the five sites with the lowest plant evenness shared the 'tri-guild' of dung saprotroph-wood saprotroph-undefined saprotroph in common. More broadly, we found that some guilds were often missing from the sites with the low functional richness. In particular, we noted that 5 of the 9 saprotroph guilds (i.e., dung, leaf, soil, wood, undefined, and undefined-wood saprotroph guilds) were present in none of the 15 most guild-poor sites. While two of these guilds-leaf and undefined-wood saprotroph-were also commonly absent in guild-rich sites, this general loss of saprotroph guilds in guild-poor sites suggests that functional losses may first affect decomposition services.
When we examined the top three most common individual guilds of fungi, only pathogen-wood saprotrophs were affected by any abiotic or biotic factors. Specifically, as the relative water level increased, pathogen-wood saprotrophs' relative abundance decreased, which indicates that changes in water management may have cascading effects on pathogen-wood saprotroph abundance. In addition to the possible direct effects of water on pathogens, increases in tree stress caused by low water levels may make trees more susceptible to pathogens, increasing the abundance of this guild. While some tree island plants prefer lower water levels, flood-tolerant trees are often more stressed when water levels drop (stress indicated by leaf loss; Sah, unpublished observation), which could increase the pathogen load in drier sites. The relationship between water level and pathogen abundance detected in this study is important to consider when managing tree islands, as lower water levels could possibly result in more fungal pathogens and lower overall fungal diversity. Due to the possible dual function of the species in this guild, future studies are needed to determine if plant pathogen-wood saprotrophs taxa are acting more as pathogens or decomposers across tree islands. This additional information on fungal function could help inform management decisions, since plant pathogens could be detrimental to tree island stability, while wood saprotrophs may be important for soil formation [36,90]. Alternatively, members of this dual-function guild may be plant pathogens that decompose the woody vegetation of their host once it declines and/or dies. A more nuanced understanding of the function and life histories of taxa within this guild may be especially important, as manipulations of Everglades hydrology as well as increases in hydroperiod have led to the loss of soil on tree islands [41,91].
While these results provide a crucial first look at the consequences of management decisions and environmental variation for the fungal community function, substantial gaps still remain to be filled. For example, we found that many taxa were not identifiable at the species level, making it difficult to use the FUNGuild database to characterize the functions of these diverse taxa (likely including taxa unique to the imperiled Everglades). Further, using amplicon sequencing alone makes it difficult to gain a full perspective into community function [92]. We suggest that future work investigate function through complementary approaches, such as functional assays, metatranscriptomics, and the direct assessment of fungal structures. For example, growth under different resource conditions and functional gene assays could inform our understanding of the primary mode of nutrition and assess some functional effects (e.g., [93,94]). Metatranscriptomic studies could further help elucidate functional responses by determining which functional genes have expression that is actively being up-and down-regulated within communities experiencing different management decisions [95]. In addition, studies of fungal structures would provide more insight into fungal diversity and how environmental variables influence the distribution of life history stages (i.e., hyphae, spores, fruiting bodies) of taxa.

Conclusions and Future Work
This study demonstrates that soil fungal communities on Everglades tree islands can be driven by abiotic and biotic factors, some of which are determined by management decisions. Our results indicate that early restoration decisions can have long-term consequences for fungal communities and suggest that a drier future in the Everglades could reduce fungal diversity on imperiled tree islands. In addition to hydrology and tree establishment, future restoration projects may want to consider the understory plant community, as our study shows a relationship between the understory community and fungal diversity and function. In our opinion, there are at least two general types of studies that are likely to be profitable going forward. First, additional manipulative experiments at macrocosm and mesocosm levels are required to tease out specific fungal functions and directly test the causality of important management decisions and environmental variables not manipulated here (e.g., understory composition). Second, field studies characterizing fungal communities on natural healthy tree islands and ghost islands would determine if fungal diversity on constructed islands closely reflects natural communities and help identify pathogens or other taxa involved in reduced ecosystem function. Collectively, the results of the work presented here and those of the suggested studies will help inform management that increases the benefit of fungal communities for the restoration and conservation of fungal diversity in this threatened ecosystem.
Supplementary Materials: The following are available online at http://www.mdpi.com/1424-2818/12/9/0324/s1: Figure S1: The layout of LILA includes eight tree islands, half of which were constructed with peat cores and half of which were constructed with limestone cores. Limestone core islands are indicated with hashed lines. Each island was divided into four quadrants in 2006 and randomly assigned to one of four tree planting density treatments. Initial tree planting density is indicated by the color of quadrants on each island; red: 1 m, green:1.67 m, blue:2.33 m, purple: 3 m spacing (i.e., distance between trees), Figure S2: An example map of one tree island within LILA. Each colored point represents one of 8 different tree species, initially planted in four density treatments. Each white plus indicates locations where soil samples were collected. The following tree species were planted in 2006 and are included in the tree species legend above: AG = Annona glabra; AR = Acer rubrum; BS = Bursera simaruba; CI = Chrysobolanus icaco; FA = Ficus aurea; IC = Ilex cassine; MC = Morella cerifera; PP = Persea palustris. In 2007, Eugenia axillaris and Myrsine floridana were added to replace trees that did not survive initial plantings, Figure S3: Distribution of fungal functional guilds (A) and trophic modes (B) across the eight experimental tree islands. (A) Each color represents the relative abundance of one of the 15 guilds to which FUNGuild assigned fungal taxa. (B) Each color represents the relative abundance of one of the four trophic modes that fungal taxa were assigned by FUNGuild. In both graphs, each vertical bar represents an individual site, and sites are grouped based on tree island (represented by numbers 1-8, island numbers 1, 4, 6, and 7 were limestone core islands and 2, 3, 5 and 8 were peat core islands), Table S1: Indicator taxa for limestone (L) and peat (P) core island communities identified by the indicator taxa analysis. A is the probability that a site belongs to a core type, given the taxa has been found at that site. B is the probability of finding that taxa in sites with that core type. The indicator values are between 0 and 1 with greater indicator values demonstrating greater specificity of a taxa to that core type.