Detection of rare variants among nuclei populating the arbuscular mycorrhizal fungal model species Rhizophagus irregularis DAOM197198.

Identifying genuine polymorphic variants is a significant challenge in sequence data analysis, although detecting low-frequency variants in sequence data is essential for estimating demographic parameters and investigating genetic processes, such as selection, within populations. Arbuscular mycorrhizal (AM) fungi are multinucleate organisms, in which individual nuclei collectively operate as a population, and the extent of genetic variation across nuclei has long been an area of scientific interest. In this study, we investigated the patterns of polymorphism discovery and the alternate allele frequency distribution by comparing polymorphism discovery in 2 distinct genomic sequence datasets of the AM fungus model species, Rhizophagus irregularis strain DAOM197198. The 2 datasets used in this study are publicly available and were generated either from pooled spores and hyphae or amplified single nuclei from a single spore. We also estimated the intraorganismal variation within the DAOM197198 strain. Our results showed that the 2 datasets exhibited different frequency patterns for discovered variants. The whole-organism dataset showed a distribution spanning low-, intermediate-, and high-frequency variants, whereas the single-nucleus dataset predominantly featured low-frequency variants with smaller proportions in intermediate and high frequencies. Furthermore, single nucleotide polymorphism density estimates within both the whole organism and individual nuclei confirmed the low intraorganismal variation of the DAOM197198 strain and that most variants are rare. Our study highlights the methodological challenges associated with detecting low-frequency variants in AM fungal whole-genome sequence data and demonstrates that alternate alleles can be reliably identified in single nuclei of AM fungi.

Stochastic nuclear organization and host-dependent allele contribution in Rhizophagus irregularis.

BACKGROUND: Arbuscular mycorrhizal (AM) fungi are arguably the most important symbionts of plants, offering a range of benefits to their hosts. However, the provisioning of these benefits does not appear to be uniform among AM fungal individuals, with genetic variation between fungal symbionts having a substantial impact on plant performance. Interestingly, genetic variation has also been reported within fungal individuals, which contain millions of haploid nuclei sharing a common cytoplasm. In the model AM fungus, Rhizophagus irregularis, several isolates have been reported to be dikaryotes, containing two genetically distinct types of nuclei recognized based on their mating-type (MAT) locus identity. However, their extremely coenocytic nature and lack of a known single nucleus stage has raised questions on the origin, distribution and dynamics of this genetic variation. RESULTS: Here we performed DNA and RNA sequencing at the mycelial individual, single spore and single nucleus levels to gain insight into the dynamic genetic make-up of the dikaryote-like R. irregularis C3 isolate and the effect of different host plants on its genetic variation. Our analyses reveal that parallel spore and root culture batches can have widely variable ratios of two main genotypes in C3. Additionally, numerous polymorphisms were found with frequencies that deviated significantly from the general genotype ratio, indicating a diverse population of slightly different nucleotypes. Changing host plants did not show consistent host effects on nucleotype ratio's after multiple rounds of subculturing. Instead, we found a major effect of host plant-identity on allele-specific expression in C3. CONCLUSION: Our analyses indicate a highly dynamic/variable genetic organization in different isolates of R. irregularis. Seemingly random fluctuations in nucleotype ratio's upon spore formation, recombination events, high variability of non-tandemly repeated rDNA sequences and host-dependent allele expression all add levels of variation that may contribute to the evolutionary success of these widespread symbionts.

Forest Fire Influence on Tomicus piniperda-Associated Fungal Communities and Phloem Nutrient Availability of Colonized Pinus sylvestris.

Forest fire is known to positively affect bark beetle populations by providing fire-damaged trees with impaired defenses for infestation. Tomicus piniperda, the common pine shoot beetle, breeds and lays eggs under the bark of stressed pine trees and is considered a serious forest pest within its native range. Wood-colonizing fungi have been hypothesized to improve substrate quality and detoxify tree defensive chemistry to indirectly facilitate tree colonization by beetles. While some bark beetle species form symbiotic associations with fungi and actively vector their partners when colonizing new trees, T. piniperda does not have mycangia or body hairs for specific vectoring of fungi. To explore the T. piniperda-associated fungal community for signs of specific association, we used ITS metabarcoding to separately characterize fungal communities associated with surface and gut of male and female beetles. We also characterized the temporal changes in fungal community and nutrient status of pine phloem with and without beetle galleries. Sampling was performed 2 years after a natural forest fire and included both burnt and unburnt sites. In our study system, we find that forest fire significantly impacts the fungal community composition associated with T. piniperda and that fire may also indirectly change nutrient availability in phloem to beetle galleries. We conclude that T. piniperda can vector fungi to newly colonized trees but the absence of positive effects on substrate quality and minimal effects of sex indicate that vectoring of associated fungal communities is not a strategy associated with the T. piniperda life cycle.

Whole genome analyses based on single, field collected spores of the arbuscular mycorrhizal fungus Funneliformis geosporum.

Arbuscular mycorrhizal (AM) fungi are ubiquitous mutualistic symbionts of most terrestrial plants and many complete their lifecycles underground. Whole genome analysis of AM fungi has long been restricted to species and strains that can be maintained under controlled conditions that facilitate collection of biological samples. There is some evidence suggesting that AM fungi can adapt to culture resulting in phenotypic and possibly also genotypic changes in the fungi. In this study, we used field isolated spores of AM fungi and identified them as Funneliformis geosporum based on morphology and phylogenetic analyses. We separately assembled the genomes of two representative spores using DNA sequences of 19 and 22 individually amplified nuclei. The genomes were compared with previously published data from other members of Glomeraceae including two strains of F. mosseae. No significant differences were observed among the species in terms of gene content, while the single nucleotide polymorphism density was higher in the strains of F. geosporum than in the strains of F. mosseae. In this study, we demonstrate that it is possible to sequence and assemble genomes from AM fungal spores sampled in the field, which opens up the possibility to include uncultured AM fungi in phylogenomic and comparative genomic analysis and to study genomic variation in natural populations of these important plant symbionts.

Effects of operational taxonomic unit inference methods on soil microeukaryote community analysis using long-read metabarcoding.

Long amplicon metabarcoding has opened the door for phylogenetic analysis of the largely unknown communities of microeukaryotes in soil. Here, we amplified and sequenced the ITS and LSU regions of the rDNA operon (around 1500 bp) from grassland soils using PacBio SMRT sequencing. We tested how three different methods for generation of operational taxonomic units (OTUs) effected estimated richness and identified taxa, and how well large-scale ecological patterns associated with shifting environmental conditions were recovered in data from the three methods. The field site at Kungsängen Nature Reserve has drawn frequent visitors since Linnaeus's time, and its species rich vegetation includes the largest population of Fritillaria meleagris in Sweden. To test the effect of different OTU generation methods, we sampled soils across an abrupt moisture transition that divides the meadow community into a Carex acuta dominated plant community with low species richness in the wetter part, which is visually distinct from the mesic-dry part that has a species rich grass-dominated plant community including a high frequency of F. meleagris. We used the moisture and plant community transition as a framework to investigate how detected belowground microeukaryotic community composition was influenced by OTU generation methods. Soil communities in both moisture regimes were dominated by protists, a large fraction of which were taxonomically assigned to Ciliophora (Alveolata) while 30%-40% of all reads were assigned to kingdom Fungi. Ecological patterns were consistently recovered irrespective of OTU generation method used. However, different methods strongly affect richness estimates and the taxonomic and phylogenetic resolution of the characterized community with implications for how well members of the microeukaryotic communities can be recognized in the data.

Soil fungal communities of ectomycorrhizal dominated woodlands across West Africa.

Forests and woodlands in the West African Guineo-Sudanian transition zone contain many tree species that form symbiotic interactions with ectomycorrhizal (ECM) fungi. These fungi facilitate plant growth by increasing nutrient and water uptake and include many fruiting body-forming fungi, including some edible mushrooms. Despite their importance for ecosystem functioning and anthropogenic use, diversity and distribution of ECM fungi is severely under-documented in West Africa. We conducted a broad regional sampling across five West African countries using soil eDNA to characterize the ECM as well as the total soil fungal community in gallery forests and savanna woodlands dominated by ECM host tree species. We subsequently sequenced the entire ITS region and much of the LSU region to infer a phylogeny for all detected soil fungal species. Utilizing a long read sequencing approach allows for higher taxonomic resolution by using the full ITS region, while the highly conserved LSU gene allows for a more accurate higher-level assignment of species hypotheses, including species without ITS-based taxonomy assignments. We detect no overall difference in species richness between gallery forests and woodlands. However, additional gallery forest plots and more samples per plot would have been needed to firmly conclude this pattern. Based on both abundance and richness, species from the families Russulaceae and Inocybaceae dominate the ECM fungal soil communities across both vegetation types. The community structure of both total soil fungi and ECM fungi was significantly influenced by vegetation types and showed strong correlation within plots. However, we found no significant difference in fungal community structure between samples collected adjacent to different host tree species within each plot. We conclude that within plots, the fungal community is structured more by the overall ECM host plant community than by the species of the individual host tree that each sample was collected from.

Long- and short-read metabarcoding technologies reveal similar spatiotemporal structures in fungal communities.

Fungi form diverse communities and play essential roles in many terrestrial ecosystems, yet there are methodological challenges in taxonomic and phylogenetic placement of fungi from environmental sequences. To address such challenges, we investigated spatiotemporal structure of a fungal community using soil metabarcoding with four different sequencing strategies: short-amplicon sequencing of the ITS2 region (300-400 bp) with Illumina MiSeq, Ion Torrent Ion S5 and PacBio RS II, all from the same PCR library, as well as long-amplicon sequencing of the full ITS and partial LSU regions (1200-1600 bp) with PacBio RS II. Resulting community structure and diversity depended more on statistical method than sequencing technology. The use of long-amplicon sequencing enables construction of a phylogenetic tree from metabarcoding reads, which facilitates taxonomic identification of sequences. However, long reads present issues for denoising algorithms in diverse communities. We present a solution that splits the reads into shorter homologous regions prior to denoising, and then reconstructs the full denoised reads. In the choice between short and long amplicons, we suggest a hybrid approach using short amplicons for sampling breadth and depth, and long amplicons to characterize the local species pool for improved identification and phylogenetic analyses.

Tremella macrobasidiata and Tremella variae have abundant and widespread yeast stages in Lecanora lichens.

Dimorphism is a widespread feature of tremellalean fungi in general, but a little-studied aspect of the biology of lichen-associated Tremella. We show that Tremella macrobasidiata and Tremella variae have an abundant and widespread yeast stage in their life cycles that occurs in Lecanora lichens. Their sexual filamentous stage is restricted to a specific lichen: T. macrobasidiata only forms basidiomata on Lecanora chlarotera hymenia and T. variae only on Lecanora varia thalli. However, the yeast stage of T. macrobasidiata is less specific and can occur in L. varia lichens, whilst all life stages of T. variae may be specific to L. varia. Contrary to the hyphal stages, the yeasts are distributed across the thalli and hymenia of Lecanora lichens, and not limited to specimens with basidiomata. Tremella macrobasidiata was present in all studied L. chlarotera, and in 59% of L. varia specimens. Only in 8% of the L. varia thalli could none of the two Tremella species be detected. Our results indicate that lichen-associated Tremella may be much more abundant and widespread than previously assumed leading to skewed estimations about their distribution ranges and lichen specificity, and raise new questions about their biology, ecology and function in the symbiosis.

In-depth Phylogenomic Analysis of Arbuscular Mycorrhizal Fungi Based on a Comprehensive Set of de novo Genome Assemblies.

Morphological characters and nuclear ribosomal DNA (rDNA) phylogenies have so far been the basis of the current classifications of arbuscular mycorrhizal (AM) fungi. Improved understanding of the evolutionary history of AM fungi requires extensive ortholog sampling and analyses of genome and transcriptome data from a wide range of taxa. To circumvent the need for axenic culturing of AM fungi we gathered and combined genomic data from single nuclei to generate de novo genome assemblies covering seven families of AM fungi. We successfully sequenced the genomes of 15 AM fungal species for which genome data was not previously available. Comparative analysis of the previously published Rhizophagus irregularis DAOM197198 assembly confirm that our novel workflow generates genome assemblies suitable for phylogenomic analysis. Predicted genes of our assemblies, together with published protein sequences of AM fungi and their sister clades, were used for phylogenomic analyses. We evaluated the phylogenetic placement of Glomeromycota in relation to its sister phyla (Mucoromycota and Mortierellomycota), and found no support to reject a polytomy. Finally, we explored the phylogenetic relationships within Glomeromycota. Our results support family level classification from previous phylogenetic studies, and the polyphyly of the order Glomerales with Claroideoglomeraceae as the sister group to Glomeraceae and Diversisporales.

Variation in hyphal production rather than turnover regulates standing fungal biomass in temperate hardwood forests.

Soil fungi link above- and belowground carbon (C) fluxes through their interactions with plants and contribute to C and nutrient dynamics through the production, turnover, and activity of fungal hyphae. Despite their importance to ecosystem processes, estimates of hyphal production and turnover rates are relatively uncommon, especially in temperate hardwood forests. We sequentially harvested hyphal ingrowth bags to quantify the rates of Dikarya (Ascomycota and Basidiomycota) hyphal production and turnover in three hardwood forests in the Midwestern United States, where plots differed in their abundance of arbuscular (AM)- vs. ectomycorrhizal (ECM)-associated trees. Hyphal production rates increased linearly with the percentage of ECM trees and annual production rates were 66% higher in ECM- than AM-dominated plots. Hyphal turnover rates did not differ across the mycorrhizal gradient (plots varying in their abundance of AM vs. ECM trees), suggesting that the greater fungal biomass in ECM-dominated plots relates to greater fungal production rather than slower fungal turnover. Differences in hyphal production across the gradient aligned with distinctly different fungal communities and activities. As ECM trees increased in dominance, fungi inside ingrowth bags produced more extracellular enzymes involved in degrading nitrogen (N)-bearing relative to C-bearing compounds, suggesting greater fungal (and possibly plant) N demand in ECM-dominated soils. Collectively, our results demonstrate that shifts in temperate tree species composition that result in changes in the dominant type of mycorrhizal association may have strong impacts on Dikarya hyphal production, fungal community composition and extracellular enzyme activity, with important consequences for soil C and N cycling.

Naming the untouchable - environmental sequences and niche partitioning as taxonomical evidence in fungi.

Due to their submerged and cryptic lifestyle, the vast majority of fungal species are difficult to observe and describe morphologically, and many remain known to science only from sequences detected in environmental samples. The lack of practices to delimit and name most fungal species is a staggering limitation to communication and interpretation of ecology and evolution in kingdom Fungi. Here, we use environmental sequence data as taxonomical evidence and combine phylogenetic and ecological data to generate and test species hypotheses in the class Archaeorhizomycetes (Taphrinomycotina, Ascomycota). Based on environmental amplicon sequencing from a well-studied Swedish pine forest podzol soil, we generate 68 distinct species hypotheses of Archaeorhizomycetes, of which two correspond to the only described species in the class. Nine of the species hypotheses represent 78% of the sequenced Archaeorhizomycetes community, and are supported by long read data that form the backbone for delimiting species hypothesis based on phylogenetic branch lengths.Soil fungal communities are shaped by environmental filtering and competitive exclusion so that closely related species are less likely to co-occur in a niche if adaptive traits are evolutionarily conserved. In soil profiles, distinct vertical horizons represent a testable niche dimension, and we found significantly differential distribution across samples for a well-supported pair of sister species hypotheses. Based on the combination of phylogenetic and ecological evidence, we identify two novel species for which we provide molecular diagnostics and propose names. While environmental sequences cannot be automatically translated to species, they can be used to generate phylogenetically distinct species hypotheses that can be further tested using sequences as ecological evidence. We conclude that in the case of abundantly and frequently observed species, environmental sequences can support species recognition in the absences of physical specimens, while rare taxa remain uncaptured at our sampling and sequencing intensity.

Building de novo reference genome assemblies of complex eukaryotic microorganisms from single nuclei.

The advent of novel sequencing techniques has unraveled a tremendous diversity on Earth. Genomic data allow us to understand ecology and function of organisms that we would not otherwise know existed. However, major methodological challenges remain, in particular for multicellular organisms with large genomes. Arbuscular mycorrhizal (AM) fungi are important plant symbionts with cryptic and complex multicellular life cycles, thus representing a suitable model system for method development. Here, we report a novel method for large scale, unbiased nuclear sorting, sequencing, and de novo assembling of AM fungal genomes. After comparative analyses of three assembly workflows we discuss how sequence data from single nuclei can best be used for different downstream analyses such as phylogenomics and comparative genomics of single nuclei. Based on analysis of completeness, we conclude that comprehensive de novo genome assemblies can be produced from six to seven nuclei. The method is highly applicable for a broad range of taxa, and will greatly improve our ability to study multicellular eukaryotes with complex life cycles.

Distribution patterns of fungal taxa and inferred functional traits reflect the non-uniform vertical stratification of soil microhabitats in a coastal pine forest.

In boreal systems, soil profiles typically consist of distinct stratified horizons, with organic layers at the surface overlying deeper mineral horizons providing microhabitat variation along a depth gradient, and vertical stratification of fungal communities along such soil profiles is commonly observed. We studied fungal community structure in a coastal pine forest along a gradient of decreasing influence from the coast. In this system, the vertical stratification pattern of soil microhabitats (defined here as organic, mineral with roots and mineral without roots: O, MR and MN, respectively) is non-uniform; organic horizons are sometimes buried under drifting sand dunes. Our results show that soil microhabitats are distinct with respect to physiochemical characteristics, community composition and OTU richness. While community composition was partly related to depth and distance from the coastal forest edge, microhabitat appeared to have the strongest influence. A closer inspection of the OTUs with the highest relative sequence abundance within each microhabitat revealed that microhabitats support functionally distinct fungal communities with respect to trophic mode and growth morphology. These results suggest that in coastal pine forests, variation in soil microhabitats contributes to the high fungal diversity found belowground and may play an important role in optimizing nutrient cycling.

Considerations and consequences of allowing DNA sequence data as types of fungal taxa.

Nomenclatural type definitions are one of the most important concepts in biological nomenclature. Being physical objects that can be re-studied by other researchers, types permanently link taxonomy (an artificial agreement to classify biological diversity) with nomenclature (an artificial agreement to name biological diversity). Two proposals to amend the International Code of Nomenclature for algae, fungi, and plants (ICN), allowing DNA sequences alone (of any region and extent) to serve as types of taxon names for voucherless fungi (mainly putative taxa from environmental DNA sequences), have been submitted to be voted on at the 11th International Mycological Congress (Puerto Rico, July 2018). We consider various genetic processes affecting the distribution of alleles among taxa and find that alleles may not consistently and uniquely represent the species within which they are contained. Should the proposals be accepted, the meaning of nomenclatural types would change in a fundamental way from physical objects as sources of data to the data themselves. Such changes are conducive to irreproducible science, the potential typification on artefactual data, and massive creation of names with low information content, ultimately causing nomenclatural instability and unnecessary work for future researchers that would stall future explorations of fungal diversity. We conclude that the acceptance of DNA sequences alone as types of names of taxa, under the terms used in the current proposals, is unnecessary and would not solve the problem of naming putative taxa known only from DNA sequences in a scientifically defensible way. As an alternative, we highlight the use of formulas for naming putative taxa (candidate taxa) that do not require any modification of the ICN.

Specificity in Arabidopsis thaliana recruitment of root fungal communities from soil and rhizosphere.

Biotic and abiotic conditions in soil pose major constraints on growth and reproductive success of plants. Fungi are important agents in plant soil interactions but the belowground mycobiota associated with plants remains poorly understood. We grew one genotype each from Sweden and Italy of the widely-studied plant model Arabidopsis thaliana. Plants were grown under controlled conditions in organic topsoil local to the Swedish genotype, and harvested after ten weeks. Total DNA was extracted from three belowground compartments: endosphere (sonicated roots), rhizosphere and bulk soil, and fungal communities were characterized from each by amplification and sequencing of the fungal barcode region ITS2. Fungal species diversity was found to decrease from bulk soil to rhizosphere to endosphere. A significant effect of plant genotype on fungal community composition was detected only in the endosphere compartment. Despite A. thaliana being a non-mycorrhizal plant, it hosts a number of known mycorrhiza fungi in its endosphere compartment, which is also colonized by endophytic, pathogenic and saprotrophic fungi. Species in the Archaeorhizomycetes were most abundant in rhizosphere samples suggesting an adaptation to environments with high nutrient turnover for some of these species. We conclude that A. thaliana endosphere fungal communities represent a selected subset of fungi recruited from soil and that plant genotype has small but significant quantitative and qualitative effects on these communities.

Dominant mycorrhizal association of trees alters carbon and nutrient cycling by selecting for microbial groups with distinct enzyme function.

While it is well established that plants associating with arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) fungi cycle carbon (C) and nutrients in distinct ways, we have a limited understanding of whether varying abundance of ECM and AM plants in a stand can provide integrative proxies for key biogeochemical processes. We explored linkages between the relative abundance of AM and ECM trees and microbial functioning in three hardwood forests in southern Indiana, USA. Across each site's 'mycorrhizal gradient', we measured fungal biomass, fungal : bacterial (F : B) ratios, extracellular enzyme activities, soil carbon : nitrogen ratio, and soil pH over a growing season. We show that the percentage of AM or ECM trees in a plot promotes microbial communities that both reflect and determine the C to nutrient balance in soil. Soils dominated by ECM trees had higher F : B ratios and more standing fungal biomass than AM stands. Enzyme stoichiometry in ECM soils shifted to higher investment in extracellular enzymes needed for nitrogen and phosphorus acquisition than in C-acquisition enzymes, relative to AM soils. Our results suggest that knowledge of mycorrhizal dominance at the stand or landscape scale may provide a unifying framework for linking plant and microbial community dynamics, and predicting their effects on ecological function.

Phosphorus cycling in deciduous forest soil differs between stands dominated by ecto- and arbuscular mycorrhizal trees.

Although much is known about how trees and their associated microbes influence nitrogen cycling in temperate forest soils, less is known about biotic controls over phosphorus (P) cycling. Given that mycorrhizal fungi are instrumental for P acquisition and that the two dominant associations - arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) fungi - possess different strategies for acquiring P, we hypothesized that P cycling would differ in stands dominated by trees associated with AM vs ECM fungi. We quantified soil solution P, microbial biomass P, and sequentially extracted inorganic and organic P pools from May to November in plots dominated by trees forming either AM or ECM associations in south-central Indiana, USA. Overall, fungal communities in AM and ECM plots were functionally different and soils exhibited fundamental differences in P cycling. Organic forms of P were more available in ECM plots than in AM plots. Yet inorganic P decreased and organic P accumulated over the growing season in both ECM and AM plots, resulting in increasingly P-limited microbial biomass. Collectively, our results suggest that P cycling in hardwood forests is strongly influenced by biotic processes in soil and that these are driven by plant-associated fungal communities.

FUNGAL BIOGEOGRAPHY. Comment on "Global diversity and geography of soil fungi".

Tedersoo et al. (Research Article, 28 November 2014, p. 1078) present a compelling study regarding patterns of biodiversity of fungi, carried out at a scale unprecedented to date for fungal biogeographical studies. The study demonstrates strong global biogeographic patterns in richness and community composition of soil fungi. What concerns us with the study is what we do not see. Unfortunately, this study underestimates the fungal diversity of one key group of soil fungi due to reliance on a single primer with known flaws.

Archaeorhizomyces borealis sp. nov. and a sequence-based classification of related soil fungal species.

The class Archaeorhizomycetes (Taphrinomycotina, Ascomycota) was introduced to accommodate an ancient lineage of soil-inhabiting fungi found in association with plant roots. Based on environmental sequencing data Archaeorhizomycetes may comprise a significant proportion of the total fungal community in soils. Yet the only species described and cultivated in this class is Archaeorhizomyces finlayi. In this paper, we describe a second species from a pure culture, Archaeorhizomyces borealis NS99-600(T) (=CBS138755(ExT)) based on morphological, physiological, and multi-locus molecular characterization. Archaeorhizomyces borealis was isolated from a root tip of a Pinus sylvestris seedling grown in a forest nursery in Lithuania. Analysis of Archaeorhizomycete species from environmental samples shows that it has a Eurasian distribution and is the most commonly observed species. Archaeorhizomyces borealis shows slow growth in culture and forms yellowish creamy colonies, characteristics that distinguish A. borealis from its closest relative A. finlayi. Here we also propose a sequence-based taxonomic classification of Archaeorhizomycetes and predict that approximately 500 species in this class remain to be isolated and described.

Effect of rainfall-induced soil geochemistry dynamics on grassland soil microbial communities.

In Mediterranean-type grassland ecosystems, the timing of rainfall events controls biogeochemical cycles, as well as the phenology and productivity of plants and animals. Here, we investigate the effect of short-term (days) soil environmental conditions on microbial community structure and composition during a natural wetting and drying cycle. Soil samples were collected from a meadow in Northern California at four time points after the first two rainfall events of the rainy season. We used 16S rRNA microarrays (PhyloChip) to track changes in bacterial and archaeal community composition. Microbial communities at time points 1 and 3 were significantly different than communities at time points 2 and 4. Based on ordination analysis, the available carbon, soil moisture, and temperature explained most of the variation in community structure. For the first time, a complementary and more comprehensive approach using linear regression and generalized logical networks were used to identify linear and nonlinear associations among environmental variables and with the relative abundance of subfamilies. Changes in soil moisture and available carbon were correlated with the relative abundance of many phyla. Only the phylum Actinobacteria showed a lineage-specific relationship to soil moisture but not to carbon or nitrogen. The results indicate that the use of a high taxonomic rank in correlations with nutritional indicators might obscure divergent subfamily-level responses to environmental parameters. An important implication of this research is that there is short-term variation in microbial community composition driven in part by rainfall fluctuation that may not be evident in long-term studies with coarser time resolution.