Tracking Fungal Community Responses to Maize Plants by DNA- and RNA-Based Pyrosequencing

We assessed soil fungal diversity and community structure at two sampling times (t1 = 47 days and t2 = 104 days of plant age) in pots associated with four maize cultivars, including two genetically modified (GM) cultivars by high-throughput pyrosequencing of the 18S rRNA gene using DNA and RNA templates. We detected no significant differences in soil fungal diversity and community structure associated with different plant cultivars. However, DNA-based analyses yielded lower fungal OTU richness as compared to RNA-based analyses. Clear differences in fungal community structure were also observed in relation to sampling time and the nucleic acid pool targeted (DNA versus RNA). The most abundant soil fungi, as recovered by DNA-based methods, did not necessary represent the most “active” fungi (as recovered via RNA). Interestingly, RNA-derived community compositions at t1 were highly similar to DNA-derived communities at t2, based on presence/absence measures of OTUs. We recovered large proportions of fungal sequences belonging to arbuscular mycorrhizal fungi and Basidiomycota, especially at the RNA level, suggesting that these important and potentially beneficial fungi are not affected by the plant cultivars nor by GM traits (Bt toxin production). Our results suggest that even though DNA- and RNA-derived soil fungal communities can be very different at a given time, RNA composition may have a predictive power of fungal community development through time.     

Diversity and Spatial Structure of Belowground Plant–Fungal Symbiosis in a Mixed Subtropical Forest of Ectomycorrhizal and Arbuscular Mycorrhizal Plants

Plant–mycorrhizal fungal interactions are ubiquitous in forest ecosystems. While ectomycorrhizal plants and their fungi generally dominate temperate forests, arbuscular mycorrhizal symbiosis is common in the tropics. In subtropical regions, however, ectomycorrhizal and arbuscular mycorrhizal plants co-occur at comparable abundances in single forests, presumably generating complex community structures of root-associated fungi. To reveal root-associated fungal community structure in a mixed forest of ectomycorrhizal and arbuscular mycorrhizal plants, we conducted a massively-parallel pyrosequencing analysis, targeting fungi in the roots of 36 plant species that co-occur in a subtropical forest. In total, 580 fungal operational taxonomic units were detected, of which 132 and 58 were probably ectomycorrhizal and arbuscular mycorrhizal, respectively. As expected, the composition of fungal symbionts differed between fagaceous (ectomycorrhizal) and non-fagaceous (possibly arbuscular mycorrhizal) plants. However, non-fagaceous plants were associated with not only arbuscular mycorrhizal fungi but also several clades of ectomycorrhizal (e.g., Russula) and root-endophytic ascomycete fungi. Many of the ectomycorrhizal and root-endophytic fungi were detected from both fagaceous and non-fagaceous plants in the community. Interestingly, ectomycorrhizal and arbuscular mycorrhizal fungi were concurrently detected from tiny root fragments of non-fagaceous plants. The plant–fungal associations in the forest were spatially structured, and non-fagaceous plant roots hosted ectomycorrhizal fungi more often in the proximity of ectomycorrhizal plant roots. Overall, this study suggests that belowground plant–fungal symbiosis in subtropical forests is complex in that it includes “non-typical” plant–fungal combinations (e.g., ectomycorrhizal fungi on possibly arbuscular mycorrhizal plants) that do not fall within the conventional classification of mycorrhizal symbioses, and in that associations with multiple functional (or phylogenetic) groups of fungi are ubiquitous among plants. Moreover, ectomycorrhizal fungal symbionts of fagaceous plants may “invade” the roots of neighboring non-fagaceous plants, potentially influencing the interactions between non-fagaceous plants and their arbuscular-mycorrhizal fungal symbionts at a fine spatial scale.     

Testing Potential Effects of Maize Expressing the Bacillus thuringiensis Cry1Ab Endotoxin (Bt Maize) on Mycorrhizal Fungal Communities via DNA- and RNA-Based Pyrosequencing and Molecular Fingerprinting

The cultivation of genetically modified (GM) crops has increased significantly over the last decades. However, concerns have been raised that some GM traits may negatively affect beneficial soil biota, such as arbuscular mycorrhizal fungi (AMF), potentially leading to alterations in soil functioning. Here, we test two maize varieties expressing the Bacillus thuringiensis Cry1Ab endotoxin (Bt maize) for their effects on soil AM fungal communities. We target both fungal DNA and RNA, which is new for AM fungi, and we use two strategies as an inclusive and robust way of detecting community differences: (i) 454 pyrosequencing using general fungal rRNA gene-directed primers and (ii) terminal restriction fragment length polymorphism (T-RFLP) profiling using AM fungus-specific markers. Potential GM-induced effects were compared to the normal natural variation of AM fungal communities across 15 different agricultural fields. AM fungi were found to be abundant in the experiment, accounting for 8% and 21% of total recovered DNA- and RNA-derived fungal sequences, respectively, after 104 days of plant growth. RNA- and DNA-based sequence analyses yielded most of the same AM fungal lineages. Our research yielded three major conclusions. First, no consistent differences were detected between AM fungal communities associated with GM plants and non-GM plants. Second, temporal variation in AMF community composition (between two measured time points) was bigger than GM trait-induced variation. Third, natural variation of AMF communities across 15 agricultural fields in The Netherlands, as well as within-field temporal variation, was much higher than GM-induced variation. In conclusion, we found no indication that Bt maize cultivation poses a risk for AMF.     

Strong coupling of plant and fungal community structure across western Amazonian rainforests

The Amazon basin harbors a diverse ecological community that has a critical role in the maintenance of the biosphere. Although plant and animal communities have received much attention, basic information is lacking for fungal or prokaryotic communities. This is despite the fact that recent ecological studies have suggested a prominent role for interactions with soil fungi in structuring the diversity and abundance of tropical rainforest trees. In this study, we characterize soil fungal communities across three major tropical forest types in the western Amazon basin (terra firme, seasonally flooded and white sand) using 454 pyrosequencing. Using these data, we examine the relationship between fungal diversity and tree species richness, and between fungal community composition and tree species composition, soil environment and spatial proximity. We find that the fungal community in these ecosystems is diverse, with high degrees of spatial variability related to forest type. We also find strong correlations between α- and β-diversity of soil fungi and trees. Both fungal and plant community β-diversity were also correlated with differences in environmental conditions. The correlation between plant and fungal richness was stronger in fungal lineages known for biotrophic strategies (for example, pathogens, mycorrhizas) compared with a lineage known primarily for saprotrophy (yeasts), suggesting that this coupling is, at least in part, due to direct plant–fungal interactions. These data provide a much-needed look at an understudied dimension of the biota in an important ecosystem and supports the hypothesis that fungal communities are involved in the regulation of tropical tree diversity.     

Dynamics of Bacterial and Fungal Communities on Decaying Salt Marsh Grass

Both bacteria and fungi play critical roles in decomposition processes in many natural environments, yet only rarely have they been studied as an integrated microbial community. Here we describe the bacterial and fungal assemblages associated with two decomposition stages of Spartina alterniflora detritus in a productive southeastern U.S. salt marsh. 16S rRNA genes and 18S-to-28S internal transcribed spacer (ITS) regions were used to target the bacterial and ascomycete fungal communities, respectively, based on DNA sequence analysis of isolates and environmental clones and by using community fingerprinting based on terminal restriction fragment length polymorphism (T-RFLP) analysis. Seven major bacterial taxa (six affiliated with the α-Proteobacteria and one with the Cytophagales) and four major fungal taxa were identified over five sample dates spanning 13 months. Fungal terminal restriction fragments (T-RFs) were informative at the species level; however, bacterial T-RFs frequently comprised a number of related genera. Amplicon abundances indicated that the salt marsh saprophyte communities have little-to-moderate variability spatially or with decomposition stage, but considerable variability temporally. However, the temporal variability could not be readily explained by either successional shifts or simple relationships with environmental factors. Significant correlations in abundance (both positive and negative) were found among dominant fungal and bacterial taxa that possibly indicate ecological interactions between decomposer organisms. Most associations involved one of four microbial taxa: two groups of bacteria affiliated with the α-Proteobacteria and two ascomycete fungi (Phaeosphaeria spartinicola and environmental isolate “4clt”).     

Interpreting Ecological Diversity Indices Applied to Terminal Restriction Fragment Length Polymorphism Data: Insights from Simulated Microbial Communities

Ecological diversity indices are frequently applied to molecular profiling methods, such as terminal restriction fragment length polymorphism (T-RFLP), in order to compare diversity among microbial communities. We performed simulations to determine whether diversity indices calculated from T-RFLP profiles could reflect the true diversity of the underlying communities despite potential analytical artifacts. These include multiple taxa generating the same terminal restriction fragment (TRF) and rare TRFs being excluded by a relative abundance (fluorescence) threshold. True community diversity was simulated using the lognormal species abundance distribution. Simulated T-RFLP profiles were generated by assigning each species a TRF size based on an empirical or modeled TRF size distribution. With a typical threshold (1%), the only consistently useful relationship was between Smith and Wilson evenness applied to T-RFLP data (TRF-E_var) and true Shannon diversity (H′), with correlations between 0.71 and 0.81. TRF-H′ and true H′ were well correlated in the simulations using the lowest number of species, but this correlation declined substantially in simulations using greater numbers of species, to the point where TRF-H′ cannot be considered a useful statistic. The relationships between TRF diversity indices and true indices were sensitive to the relative abundance threshold, with greatly improved correlations observed using a 0.1% threshold, which was investigated for comparative purposes but is not possible to consistently achieve with current technology. In general, the use of diversity indices on T-RFLP data provides inaccurate estimates of true diversity in microbial communities (with the possible exception of TRF-E_var). We suggest that, where significant differences in T-RFLP diversity indices were found in previous work, these should be reinterpreted as a reflection of differences in community composition rather than a true difference in community diversity.     

Watershed-Scale Fungal Community Characterization along a pH Gradient in a Subsurface Environment Cocontaminated with Uranium and Nitrate

The objective of this study was to characterize fungal communities in a subsurface environment cocontaminated with uranium and nitrate at the watershed scale and to determine the potential contribution of fungi to contaminant transformation (nitrate attenuation). The abundance, distribution, and diversity of fungi in subsurface groundwater samples were determined using quantitative and semiquantitative molecular techniques, including quantitative PCR of eukaryotic small-subunit rRNA genes and pyrosequencing of fungal internal transcribed spacer (ITS) regions. Potential bacterial and fungal denitrification was assessed in sediment-groundwater slurries amended with antimicrobial compounds and in fungal pure cultures isolated from the subsurface. Our results demonstrate that subsurface fungal communities are dominated by members of the phylum Ascomycota, and a pronounced shift in fungal community composition occurs across the groundwater pH gradient at the field site, with lower diversity observed under acidic (pH <4.5) conditions. Fungal isolates recovered from subsurface sediments, including cultures of the genus Coniochaeta, which were detected in abundance in pyrosequence libraries of site groundwater samples, were shown to reduce nitrate to nitrous oxide. Denitrifying fungal isolates recovered from the site were classified and found to be distributed broadly within the phylum Ascomycota and within a single genus of the Basidiomycota. Potential denitrification rate assays with sediment-groundwater slurries showed the potential for subsurface fungi to reduce nitrate to nitrous oxide under in situ acidic pH conditions.     

The annual invasive plant Impatiens glandulifera reduces hyphal biomass of soil fungi in deciduous forests

Soil fungi play a crucial role in ecosystem functioning and there is increasing evidence that exotic plants invading forests can affect soil fungal communities. We examined potential effects of the invasive plant Impatiens glandulifera on hyphal biomass of ectomycorrhizal fungi, their genetic diversity and the diversity of other soil fungi in deciduous forests in Switzerland. We compared invaded patches with patches where I. glandulifera had been removed, by establishing pairs of 3-m long transect lines at the edge of seven areas of either type. Along the transects we assessed the length of ectomycorrhizal fungal hyphae using the ‘ingrowth mesh bag method’, and used terminal restriction fragment length polymorphism (T-RFLP) analysis to examine fungal genetic diversity. The invasive plant reduced fungal hyphal biomass by 30–80%: the reduction was largest in the centre of the patch. I. glandulifera did not alter fungal richness, but affected the composition of fungal communities. This is probably the result of a decrease of mycorrhizal fungi, coupled with an increase of saprotrophic fungi. Our findings demonstrate the adverse impacts of an annual invasive plant species on both fungal hyphal biomass and the composition of soil fungal communities. This may negatively affect forest nutrient and carbon cycling, soil stability and the functionality of the fungal community, with major consequences for forest ecosystem functioning.     

A distinctive fungal community inhabiting cryoconite holes on glaciers in Svalbard

Cryoconite holes on glacier surfaces are ice-cold hot spots of microbial diversity and activity but still little is known about their fungal inhabitants. We provide the first report of distinctive fungal communities in cryoconite debris from three valley glaciers at Kongsfjorden, Svalbard. Multivariate analysis of terminal-restriction fragment length polymorphism (T-RFLP) profiles of rRNA ITS amplicons revealed that quite distinct fungal communities were found in cryoconite holes compared with soils from adjacent moraine and tundra sites, and that communities on glaciers with contrasting ice-surface hydrology also differed. Most of the fungi cultured from cryoconite sediment were basidiomycetous yeasts or filamentous Ascomycota (Helotiales/Pleosporales). The latter included aeroaquatic fungi, such as Articulospora and Varicosporium, implying a role for these important freshwater decomposers in the carbon dynamics of cryoconite holes. Matching of the dominant peaks from T-RFLP analysis to predicted peaks of cultured isolates confirmed the abundance of these aeroaquatic fungi but also revealed that most of the dominant T-RFLP peaks did not match any cultured isolates. Considering the prevalence and endangerment of glacial environments worldwide, these findings would suggest that their potential as reservoirs of fungal diversity should not be overlooked.     

Co-existing grass species have distinctive arbuscular mycorrhizal communities

Arbuscular mycorrhizal (AM) fungi are biotrophic symbionts colonizing the majority of land plants, and are of major importance in plant nutrient supply. Their diversity is suggested to be an important determinant of plant community structure, but the influence of host-plant and environmental factors on AM fungal community in plant roots is poorly documented. Using the terminal restriction fragment length polymorphism (T-RFLP) strategy, the diversity of AM fungi was assessed in 89 roots of three grass species (Agrostis capillaris, Festuca rubra, Poa pratensis) that co-occurred in the same plots of a field experiment. The impact of different soil amendments (nitrogen, lime, nitrogen and lime) and insecticide application on AM fungal community was also studied. The level of diversity found in AM fungal communities using the T-RFLP strategy was consistent with previous studies based on clone libraries. Our results clearly confirm that an AM fungal host-plant preference exists, even between different grass species. AM communities colonizing A. capillaris were statistically different from the others (P < 0.05). Although grass species evenness changed in amended soils, AM fungal community composition in roots of a given grass species remained stable. Conversely, in plots where insecticide was applied, we found higher AM fungal diversity and, in F. rubra roots, a statistically different AM fungal community.     

Endophytic Fungal Communities Associated with Vascular Plants in the High Arctic Zone Are Highly Diverse and Host-Plant Specific

This study assessed the diversity and distribution of endophytic fungal communities associated with the leaves and stems of four vascular plant species in the High Arctic using 454 pyrosequencing with fungal-specific primers targeting the ITS region. Endophytic fungal communities showed high diversity. The 76,691 sequences obtained belonged to 250 operational taxonomic units (OTUs). Of these OTUs, 190 belonged to Ascomycota, 50 to Basidiomycota, 1 to Chytridiomycota, and 9 to unknown fungi. The dominant orders were Helotiales, Pleosporales, Capnodiales, and Tremellales, whereas the common known fungal genera were Cryptococcus, Rhizosphaera, Mycopappus, Melampsora, Tetracladium, Phaeosphaeria, Mrakia, Venturia, and Leptosphaeria. Both the climate and host-related factors might shape the fungal communities associated with the four Arctic plant species in this region. These results suggested the presence of an interesting endophytic fungal community and could improve our understanding of fungal evolution and ecology in the Arctic terrestrial ecosystems.     

Topographical Mapping of the Rainbow Trout (Oncorhynchus mykiss) Microbiome Reveals a Diverse Bacterial Community with Antifungal Properties in the Skin

The mucosal surfaces of wild and farmed aquatic vertebrates face the threat of many aquatic pathogens, including fungi. These surfaces are colonized by diverse symbiotic bacterial communities that may contribute to fight infection. Whereas the gut microbiome of teleosts has been extensively studied using pyrosequencing, this tool has rarely been employed to study the compositions of the bacterial communities present on other teleost mucosal surfaces. Here we provide a topographical map of the mucosal microbiome of an aquatic vertebrate, the rainbow trout (Oncorhynchus mykiss). Using 16S rRNA pyrosequencing, we revealed novel bacterial diversity at each of the five body sites sampled and showed that body site is a strong predictor of community composition. The skin exhibited the highest diversity, followed by the olfactory organ, gills, and gut. Flectobacillus was highly represented within skin and gill communities. Principal coordinate analysis and plots revealed clustering of external sites apart from internal sites. A highly diverse community was present within the epithelium, as demonstrated by confocal microscopy and pyrosequencing. Using in vitro assays, we demonstrated that two Arthrobacter sp. skin isolates, a Psychrobacter sp. strain, and a combined skin aerobic bacterial sample inhibit the growth of Saprolegnia australis and Mucor hiemalis, two important aquatic fungal pathogens. These results underscore the importance of symbiotic bacterial communities of fish and their potential role for the control of aquatic fungal diseases.     

Fungal diversity in galls of baldcypress trees

The baldcypress midge (Taxodiomyia cupressi and Taxodiomyia cupressiananassa) forms a gall that originates from leaf tissue. Female insects may inoculate galls with fungi during oviposition, or endophytes from the leaf tissue may grow into the gall interior. We investigated fungal diversity inside of baldcypress galls, comparing the gall communities to leaves and comparing fungal communities in galls that had successful emergence versus no emergence of midges or parasitoids. Galls of midges that successfully emerged were associated with diverse gall fungal communities, some of which were the same as the fungi found in surrounding leaves. Galls with no insect emergence were characterized by relatively low fungal diversity.     

Airborne and Grain Dust Fungal Community Compositions Are Shaped Regionally by Plant Genotypes and Farming Practices

Chronic exposure to airborne fungi has been associated with different respiratory symptoms and pathologies in occupational populations, such as grain workers. However, the homogeneity in the fungal species composition of these bioaerosols on a large geographical scale and the different drivers that shape these fungal communities remain unclear. In this study, the diversity of fungi in grain dust and in the aerosols released during harvesting was determined across 96 sites at a geographical scale of 560 km² along an elevation gradient of 500 m by tag-encoded 454 pyrosequencing of the internal transcribed spacer (ITS) sequences. Associations between the structure of fungal communities in the grain dust and different abiotic (farming system, soil characteristics, and geographic and climatic parameters) and biotic (wheat cultivar and previous crop culture) factors were explored. These analyses revealed a strong relationship between the airborne and grain dust fungal communities and showed the presence of allergenic and mycotoxigenic species in most samples, which highlights the potential contribution of these fungal species to work-related respiratory symptoms of grain workers. The farming system was the major driver of the alpha and beta phylogenetic diversity values of fungal communities. In addition, elevation and soil CaCO₃ concentrations shaped the alpha diversity, whereas wheat cultivar, cropping history, and the number of freezing days per year shaped the taxonomic beta diversity of these communities.     

Redox Fluctuation Structures Microbial Communities in a Wet Tropical Soil

Frequent high-amplitude redox fluctuation may be a strong selective force on the phylogenetic and physiological composition of soil bacterial communities and may promote metabolic plasticity or redox tolerance mechanisms. To determine effects of fluctuating oxygen regimens, we incubated tropical soils under four treatments: aerobic, anaerobic, 12-h oxic/anoxic fluctuation, and 4-day oxic/anoxic fluctuation. Changes in soil bacterial community structure and diversity were monitored with terminal restriction fragment length polymorphism (T-RFLP) fingerprints. These profiles were correlated with gross N cycling rates, and a Web-based phylogenetic assignment tool was used to infer putative community composition from multiple fragment patterns. T-RFLP ordinations indicated that bacterial communities from 4-day oxic/anoxic incubations were most similar to field communities, whereas those incubated under consistently aerobic or anaerobic regimens developed distinctly different molecular profiles. Terminal fragments found in field soils persisted either in 4-day fluctuation/aerobic conditions or in anaerobic/12-h treatments but rarely in both. Only 3 of 179 total fragments were ubiquitous in all soils. Soil bacterial communities inferred from in silico phylogenetic assignment appeared to be dominated by Actinobacteria (especially Micrococcus and Streptomycetes), “Bacilli,” “Clostridia,” and Burkholderia and lost significant diversity under consistently or frequently anoxic incubations. Community patterns correlated well with redox-sensitive processes such as nitrification, dissimilatory nitrate reduction to ammonium (DNRA), and denitrification but did not predict patterns of more general functions such as N mineralization and consumption. The results suggest that this soil's indigenous bacteria are highly adapted to fluctuating redox regimens and generally possess physiological tolerance mechanisms which allow them to withstand unfavorable redox periods.     

Analyzing salt-marsh fungal diversity: comparing ARISA fingerprinting with clone sequencing and pyrosequencing

In a previous study from our laboratory we used automated ribosomal intergenic spacer analysis (ARISA) to assess salt-marsh fungal diversity (Torzilli et al. 2006). The results demonstrated that different salt-marsh plants harbor distinct fungal communities, thereby supporting the hypothesis that substratum type is an important factor in determining fungal community composition. However, ARISA of several pure cultures of salt-marsh fungi indicated that an operational taxonomic unit (OUT) in an ARISA community profile may represent more than one taxon. To assess the extent to which such ambiguity might have affected the interpretation of our ARISA fingerprinting, we have now fingerprinted and sequenced clones derived from the same fungal DNA used for our ARISA community profiles. Results from this confirmed that an ARISA OTU may represent multiple taxa and that a given taxon may be represented by more than one OTU. Nonetheless, sequencing still confirmed the importance of substratum in determining community composition, and indicated that despite ambiguities associated with OTU's, ARISA may be used to provide a quick snapshot of diversity which can be further refined using sequencing methods. In addition, we compared the fungal diversity from short-form Spartina alterniflora as revealed by clone sequencing with that obtained from pyrosequencing, which avoids the cloning biases of traditional sequencing, and provide greatly expanded depth of coverage. Pyrosequencing significantly enhanced the characterization of fungal diversity compared to traditional clone sequencing.     

Phylogenetic diversity, host-specificity and community profiling of sponge-associated bacteria in the northern Gulf of Mexico

Background

Marine sponges can associate with abundant and diverse consortia of microbial symbionts. However, associated bacteria remain unexamined for the majority of host sponges and few studies use phylogenetic metrics to quantify symbiont community diversity. DNA fingerprinting techniques, such as terminal restriction fragment length polymorphisms (T-RFLP), might provide rapid profiling of these communities, but have not been explicitly compared to traditional methods.

Methodology/Principal Findings

We investigated the bacterial communities associated with the marine sponges Hymeniacidon heliophila and Haliclona tubifera, a sympatric tunicate, Didemnum sp., and ambient seawater from the northern Gulf of Mexico by combining replicated clone libraries with T-RFLP analyses of 16S rRNA gene sequences. Clone libraries revealed that bacterial communities associated with the two sponges exhibited lower species richness and lower species diversity than seawater and tunicate assemblages, with differences in species composition among all four source groups. T-RFLP profiles clustered microbial communities by source; individual T-RFs were matched to the majority (80.6%) of clone library sequences, indicating that T-RFLP analysis can be used to rapidly profile these communities. Phylogenetic metrics of community diversity indicated that the two sponge-associated bacterial communities include dominant and host-specific bacterial lineages that are distinct from bacteria recovered from seawater, tunicates, and unrelated sponge hosts. In addition, a large proportion of the symbionts associated with H. heliophila were shared with distant, conspecific host populations in the southwestern Atlantic (Brazil).

Conclusions/Significance

The low diversity and species-specific nature of bacterial communities associated with H. heliophila and H. tubifera represent a distinctly different pattern from other, reportedly universal, sponge-associated bacterial communities. Our replicated sampling strategy, which included samples that reflect the ambient environment, allowed us to differentiate resident symbionts from potentially transient or prey bacteria. Pairing replicated clone library construction with rapid community profiling via T-RFLP analyses will greatly facilitate future studies of sponge-microbe symbioses.