Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Associations between an Invasive Plant (Taeniatherum caput-medusae, Medusahead) and Soil Microbial Communities

  • Elise S. Gornish ,

    egornish@ucdavis.edu

    Affiliation Department of Plant Sciences, University of California Davis, Davis, California, United States of America

    ⨯
  • Noah Fierer,

    Affiliations Cooperative Institute for Research in Environmental Sciences, University of Colorado, Boulder, Colorado, United States of America, Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, Colorado, United States of America

    ⨯
  • Albert Barberán

    Affiliation Cooperative Institute for Research in Environmental Sciences, University of Colorado, Boulder, Colorado, United States of America

    ⨯

Abstract

Understanding plant-microbe relationships can be important for developing management strategies for invasive plants, particularly when these relationships interact with underlying variables, such as habitat type and seedbank density, to mediate control efforts. In a field study located in California, USA, we investigated how soil microbial communities differ across the invasion front of Taeniatherum caput-medusae (medusahead), an annual grass that has rapidly invaded most of the western USA. Plots were installed in habitats where medusahead invasion is typically successful (open grassland) and typically not successful (oak woodland). Medusahead was seeded into plots at a range of densities (from 0–50,000 seeds/m2) to simulate different levels of invasion. We found that bacterial and fungal soil community composition were significantly different between oak woodland and open grassland habitats. Specifically, ectomycorrhizal fungi were more abundant in oak woodlands while arbuscular mycorrhizal fungi and plant pathogens were more abundant in open grasslands. We did not find a direct effect of medusahead density on soil microbial communities across the simulated invasion front two seasons after medusahead were seeded into plots. Our results suggest that future medusahead management initiatives might consider plant-microbe interactions.

Citation: Gornish ES, Fierer N, Barberán A (2016) Associations between an Invasive Plant (Taeniatherum caput-medusae, Medusahead) and Soil Microbial Communities. PLoS ONE 11(9): e0163930. https://doi.org/10.1371/journal.pone.0163930

Editor: Zhili He, University of Oklahoma, UNITED STATES

Received: May 12, 2016; Accepted: September 16, 2016; Published: September 29, 2016

Copyright: © 2016 Gornish et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Soil sample information, phylotype abundance tables, and bacterial and fungal representative sequences are publicly available in FigShare (10.6084/m9.figshare.3113125).

Funding: ESG was supported by the Department of Plant Sciences and the Division of Agriculture and Natural Resources at University of California Davis. Funding to support this work was provided to NF from the National Science Foundation (DEB-0953331). AB was supported by a James S. McDonnell Postdoctoral Fellowship. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Plant communities are typically composed of a combination of native and non-native species. The majority of these non-native species are benign, demonstrating little to no negative effect on neighboring organisms. However, a small fraction of these non-native plants are characterized as invasive because they are able to profoundly modify local plant and animal communities, nutrient cycling, hydrological regimes and fire frequency [1–2]. Not only do these impacts erode biodiversity and devalue ecosystem services, but they can also enhance further invasion by con- and heterospecific exotics (e.g. [3]).

Soil microbial communities might mediate relationships between invasive plant species and their ecosystem impacts [4–6]. Soil microbial communities, which are typically dominated by fungi and bacteria, can be altered by invasive plants directly through growth facilitation or inhibition near the root zone [7], and indirectly through changes in abiotic conditions (e.g. pH or nutrient availability) that occur in tandem with weed establishment [8]. For example, species-specific effects of non-native grasses on soil nutrients have been shown to subsequently modify soil microbial community composition, biomass, and bacterial:fungal ratios [9]. In addition to being vulnerable to impacts from aboveground plant dynamics, soil microbial communities may also play an important role in mediating the success of plant invasions [4]. For example, extant soil biota have been shown to enhance invasion success of some of the world’s most noxious invasive plants, such as exotic knotweeds (Fallopia spp.) [10]. In general, it has been concluded that invasive species may be differentially affected by soil bacterial or fungal pathogens as compared to native plant species [11, 8], but see [12], which could be important for developing reliable control strategies for invasive plants that demonstrate resistance to current management efforts.

The invasive annual winter grass medusahead (Taeniatherum caput-medusae [L.] Nevski) has invaded much of the western USA and has been shown to decrease soil carbon stocks, reduce native plant diversity, and enhance fire frequency [13]. A recent meta-analysis of medusahead control outcomes in annual grassland and intermountain regions identified large variance in the effectiveness of conventional approaches for managing medusahead [14], suggesting that underlying variables, such as habitat type and seedbank density, might mediate control efforts. Despite increasing recognition that bacterial and fungal communities can influence plant invasion dynamics, only two published studies have investigated the direct relationship between medusahead and soil microbial communities [15–16]. These studies have conflicting results, suggesting both that the interaction between medusahead and soil microorganisms might or might not enhance its own invasion. Understanding medusahead effects on the soil microbial community is critical for enhancing predictions of invasion effects and for developing effective management strategies.

We investigated how soil microbial communities differ across the invasion front of medusahead in experimental plots in open grassland and oak woodland habitat in the Sierra Foothill region of California, USA. We attempted to understand (1) if medusahead modifies soil microbial communities across the invasion front (simulated by differences in seed density) within systems; and (2) how soil microbial communities differ between areas where medusahead invasion is successful (open grassland habitat) and not successful (oak woodland habitat), and the factors that could be responsible for these differences. We hypothesized that medusahead would modify the soil microbial communities within each habitat. We expected this for two reasons. First, early work on this species [15], as well as more recent work on other invasive annual grasses with similar invasion dynamics to medusahead, have demonstrated linkages between the soil microbes and invasion success [17]. Second, plant-soil interactions are common in the savannah/oak woodlands of California [18], so we would expect strong effects from the extant soil microorganisms. We also hypothesized that soil microbial communities (in particular symbiotic and pathogenic fungi) would differ between areas where medusahead invasion is typically successful (open grasslands) and typically not successful (oak woodlands). In California, invasion of winter annual grasses can be strongly limited within oak canopies [19], possibly due to the different microbial communities associated with oak trees compared to adjacent open grasslands [20–21]. Through shading, litter input, and hydraulic lift, Mediterranean oak trees can also modify a wide variety of soil edaphic factors, such as pH and organic matter concentrations, which directly influence soil microbial communities (e.g. [22]). At present, we do not know what role, if any, soil microbial communities play in mediating the likelihood of successful medusahead invasions.

Materials and Methods

Study site and experiment

The study site was located on a research reserve in Yuba County, California (39°14´N, 121°18´W), which experiences a Mediterranean climate of hot, dry summers and cool, wet winters. Permission to perform this experiment was granted by the owner of the property, the University of California Division of Agriculture and Natural Resources. The field study did not involve endangered or protected species.

Mean annual precipitation is 75cm and mean annual temperature is 17.8°C. Soils at the site are fine-loamy, mixed, superactive Ultic Haploxeralfs and fine, mixed, superactive Typic Rhodoxeralfs. Soil pH ranges between 5.7 and 6.2. The experiment is located in an annual grassland system that is irregularly interrupted by small patches of winter-deciduous blue oak (Quercus douglasii), and evergreen interior live oak (Q. wislizeni) that provide approximately 40% shade. The area has experienced seasonal low intensity grazing by livestock since the 1960’s.

The experimental set-up is described fully in [19, 23]. Briefly, plots (1m2, separated by 2m) were installed in open grassland habitat and paired oak woodland habitat. The two habitats differ in the identity of dominant herbaceous species and the presence of leaf litter [19]. Average soil temperatures and soil moisture in the top 5cm of soil were also lower in oak woodland plots (18.4°C and 1.6%, respectively) compared to open grassland plots (22.2°C and 6.4%, respectively).

The experimental site was mowed, solarized to enhance seedbank germination and then sprayed with glyphosate herbicide to kill existing and newly germinated plants. In September 2013, fully replicated (n = 4) plots were hand seeded with one of five densities of field-collected medusahead (0, 100, 1000, 10000, and 50000 seeds/m2), mixed in with 500 grams of medusahead thatch. Immediately following the addition of medusahead seed, 6,000 seeds each of neighboring grass species (annual rye, and Blando brome) and 4,000 seeds of a clover mix were added—for a total of 16,000 neighbor seeds—to maintain a realistic competitive environment. Medusahead tiller density the following season reflected the seeding rate. This treatment was expected to simulate differences in medusahead invasion intensity from low to high infestation [24].

A defoliation treatment was applied to half the plots in April 2014. This treatment was intended to simulate a mowing or grazing regime included in a typical management program. Defoliation was applied when 75% of the medusahead tillers were in the ‘boot’ stage within a plot. All standing biomass in treatment plots was clipped using electronic shears positioned approximately 15cm above the soil surface.

To understand how factors associated with the oak woodland habitat might contribute to medusahead invasion, an additional set of replicated 1m2 plots were installed in the open grassland sites. Treatments were deployed in order to simulate environmental factors associated with the oak woodland habitat, and included the presence of shading, the presence of oak litter, and the presence of both shading and litter. Shading was applied via 50% shade cloth suspended over the plots, and litter was applied by collecting 500g of litter from under paired oak canopies and distributing it evenly over treatment plots. Medusahead seeds (at a density of 50000 seeds/m2) were then introduced to these plots.

Soil sampling

In April 2015, we collected surface soil cores (7cm depth), where the majority of fungal and biomass is present [25] in four random locations in each plot. Soil from each plot was mixed together in the field, sieved and placed in a plastic bag. Soil samples were frozen in the field and transported at -20°C to the University of Colorado, Boulder for microbial analyses.

Molecular analyses

Microbial diversity was assessed using high-throughput sequencing methods to describe the composition of taxonomic marker gene sequences. For bacterial analyses, we sequenced the V4 hypervariable region of the 16S rRNA gene, using the 515-F (GTGCCAGCMGCCGCGGTAA) and 806-R (GGACTACHVGGGTWTCTAAT) primer pair [26]. For the fungal analyses, we sequenced the first internal transcribed spacer (ITS1) region of the rRNA operon, using the ITS1-F (CTTGGTCATTTAGAGGAAGTAA) and ITS2 (GCTGCGTTCTTCATCGATGC) primer pair [27]. The primers included Illumina adapters and an error-correcting 12-bp barcode unique to each sample. PCR products were quantified using the PicoGreen dsDNA assay, and pooled together in equimolar concentrations for sequencing on an Illumina MiSeq instrument. All sequencing runs were conducted at the University of Colorado Next Generation Sequencing Facility.

Reads were demultiplexed using a custom Python script (https://github.com/leffj/helper-code-for-uparse), with quality filtering and phylotype clustering conducted using UPARSE [28]. For quality filtering, we used a maxee value of 0.5 (that is, a maximum of 0.5 nucleotides incorrectly assigned in every sequence). Singleton sequences were removed prior to phylotype clustering. Quality filtered sequence reads were then mapped to phylotypes at the 97% similarity threshold. Phylotype taxonomy was assigned using the Ribosomal Database Project (RDP) classifier with a confidence threshold of 0.5 [29] trained on the 16S rRNA Greengenes database [30] or the ITS rRNA UNITE database [31], for bacteria and fungi respectively. Sequences representing any phylotypes unclassified at the domain level or classified as mitochondria, chloroplasts, archaea or protists were removed. Subsequently, we removed potential contaminants (i.e. phylotypes with abundances greater than 1% in the blanks and no-template controls [32]), and we normalized the sequence counts using a cumulative-sum scaling approach [33]. We used FUNGuild to identify fungi functional guilds [34]. Soil sample information, phylotype abundance tables, and bacterial and fungal representative sequences are publicly available in FigShare (10.6084/m9.figshare.3113125).

Statistical analyses

Soil microbial community similarity patterns were represented by non-metric multidimensional scaling (NMDS) using the Bray-Curtis distance metric. We used nested analysis of variance (ANOVA) and permutational multivariate analysis of variance (PERMANOVA) based on 1,000 permutations [35] to assess the explanatory power of the different treatments on soil microbial richness and community similarity patterns, respectively. Differences in the proportion of taxa and fungal functional guilds were tested using non-parametric Wilcoxon tests after false discovery rate (FDR) correction [36]. All multivariate statistical analyses were implemented in the R environment (www.r-project.org) using the vegan package (vegan.r-forge.r-project.org).

Results

The total number of phylotypes across all soil samples was 5604 and 4349, for bacteria and fungi respectively (S1 Fig). The average number of phylotypes per soil sample was 1781 and 687, for bacteria and fungi respectively. Oak woodland soil samples tended to harbor more bacterial and fungal phylotypes than open grassland samples (although the differences were not statistically significant; ANOVA P > 0.05; Fig 1A and 1B respectively; see S2 Fig for Shannon diversity). Both bacterial and fungal soil community composition were significantly different between open grasslands and oak woodlands (PERMANOVA R2 = 0.24, P < 0.001, R2 = 0.28, P < 0.001, respectively; Fig 1C and 1D respectively). Bacteria from the proteobacterial classes alpha, beta, gamma and delta, and Acidobacteria subgroup 6 as well as fungi from the Pezizomycetes, Agaricomycetes and Eurotiomycetes classes were more abundant in oak soils than in grassland soils (Wilcoxon test P < 0.01 after FDR; S3 Fig; see S1 Table for results at the genus level). Acidobacteria from the classes Solibacteres and Acidobacteriia, Spartobacteria, Planctomycetia and Phycisphaerae, and fungi from the classes Dothideomycetes, Sordariomycetes and Glomeromycetes were more abundant in open grassland soils (Wilcoxon test P < 0.01 after FDR; S3 Fig; see S1 Table for results at the genus level). As expected from the contrasted fungal taxonomic composition, oak woodlands were different from open grasslands based on the overall abundance of fungal functional guilds (S4 Fig). Ectomycorrhizal fungi (in particular, species from the Tuber and Tomentella genera) were significantly more abundant in oak woodlands while arbuscular mycorrhizal fungi and plant pathogens (in particular, Drechslera and Fusarium species) were more abundant in open grasslands (Wilcoxon test P < 0.01 after FDR; Fig 2).

thumbnail
Download:
Fig 1.

Observed richness (number of different phylotypes per sample) for bacterial and fungal soil communities (A, B). Community similarity patterns for bacterial and fungal soil communities using non-metric multidimensional scaling (C, D).

https://doi.org/10.1371/journal.pone.0163930.g001

thumbnail
Download:
Fig 2. Significant differences (Wilcoxon test P < 0.01 after false discovery rate correction) in the abundance of fungal functional guilds between oak woodland soil samples (green) and open grassland soil samples (blue).

https://doi.org/10.1371/journal.pone.0163930.g002

We found no significant effects of seed density treatments or clipping treatments on soil microbial richness (ANOVA P > 0.05 for both bacteria and fungi; Fig 1A and 1B) or microbial community composition (PERMANOVA P > 0.05 for both bacteria and fungi; Fig 1C and 1D) within each habitat type. Likewise, for those plots where we simulated the oak environment, we observed no significant effects of the litter or shade treatments on the richness of soil bacterial communities (Wilcoxon test P > 0.05; Fig 3A) or bacterial community composition (PERMANOVA P > 0.05; Fig 3C). For fungi, we detected weak significant effects of litter on richness (Wilcoxon test P < 0.01; Fig 3B) and on fungal community composition (PERMANOVA R2 = 0.09 P < 0.01; Fig 3D).

thumbnail
Download:
Fig 3.

Observed richness (number of different phylotypes per sample) for bacterial and fungal soil communities (A, B). Community similarity patterns for bacterial and fungal soil communities using non-metric multidimensional scaling (C, D).

https://doi.org/10.1371/journal.pone.0163930.g003

Discussion

Soil biota has been implicated in the facilitation of invasive plant dominance [5, 37]. However, not all invasive species support plant-soil microbe feedbacks as a driver of invasion [38–39]. We attempted to identify how relationships between the weedy annual grass medusahead and soil microorganisms might mediate invasion success. Unexpectedly, we did not find evidence for medusahead density effects on the soil microbial communities across a simulated invasion front in either habitat. This supports other studies that have documented instances where soil communities are unresponsive to the presence of invasive weeds in arid grasslands [40] and in other systems [41]. However, several aspects of our experiment could hinder our ability to capture an existing relationship between medusahead seed density and soil microbial communities. First, we assessed the relationship between medusahead seed density and bulk soil microbial communities and not rhizosphere communities. Rhizosphere microbial communities are different from bulk soil communities [42], so the potential effects of medusahead invasion intensity on microbial communities may be observable at smaller spatial scales in the rhizosphere (but see [43]). Second, although this study was conducted across two growing seasons, there may have been insufficient time for soil microbial communities to respond to different seed densities of medusahead [44]. Although microbial communities have been shown to respond to changes in aboveground plant communities in as little as a month [45], these communities could be especially slow to respond to the presence of weeds in environments where soil edaphic factors are slow to change in response to invasion. Moreover, in California grasslands, soil appears to be particularly buffered from aboveground changes [40]. Third, extracellular microbial DNA and DNA from dead cells can persist in soils for years and thus, obscure DNA-based present estimates of soil microbial composition [46]. Finally, as this study only assessed the composition of the microbial communities, we cannot eliminate the possibility that medusahead seed density can influence the activity and function of belowground soil communities.

Using reciprocal soil transplant experiments, [47] reported higher medusahead biomass in introduced soil than in native soil, which suggests that medusahead success is partially due to release from native soil pathogens [48–49]. In addition to escaping from soil pathogens, certain plant invasive species have been shown to accumulate local pathogens [11, 49]. Although previous studies have reported that medusahead is sensitive to antagonistic fungi [50–51], we observed important differences in bacterial and fungal community composition between open grassland sites where medusahead is typically found in high densities and oak woodland sites where medusahead is typically found in low densities. Specifically, we found significantly higher abundances of fungal pathogens in open grasslands compared to oak woodland habitats. Environmental conditions in the grassland habitat are likely more ideal for both soil and foliar fugal pathogens, which can be important drivers of above ground plant dynamics (e.g. [52]). Indeed, these pathogens have been documented in grasslands in other studies (e.g. [53]). These results collectively highlight the potential contribution of microbial mechanisms (e.g., via pathogen accumulation) to medusahead dominance in California grasslands. Because the differences in bacterial and fungal communities exist in the absence of medusahead, it is more likely that favorable grassland soil microbial communities facilitate medusahead establishment instead of resulting from the invasion itself.

In addition to negative interactions, a large number of plant species establish symbiotic associations with soil microorganisms (in particular with mycorrhizal fungi and nitrogen-fixing bacteria, [6, 54]). In this experiment we detected a higher proportion of ectomycorrhizal fungi in soil samples from oak woodland habitats. This result is expected as ectomycorrhizal fungi are important to oak trees for acquiring nutrients and for increasing root absorptive area [18, 55]. Our results also show that oak litter (rather than shade) influence overall soil fungal community composition and richness, but not soil bacterial community composition and richness. Although plant litter inputs can change important environmental conditions for soil bacteria such as pH and base cation content [42], soil fungi are key decomposers of plant necromass and depend more directly on leaf litter than bacteria [56]. We also observed a significantly higher proportion of arbuscular mycorrhizal fungi in open grasslands. Given the generally non-specific interactions with arbuscular mycorrhizal fungi, it has been proposed that several invasive plants make use of these fungi to enhance their success [57–60]; but see [61]. Other invasive plants (for example, the garlic mustard Alliaria petiolata) inhibit mycorrhizal fungi on which natives depend [62]. Collectively, this work suggests that biocontrol and management initiatives should consider the potentially beneficial plant-microbe interactions rather than just focusing on antagonistic relationships.

The context dependency associated with invasion success and weed management efficiency is well documented for both medusahead as well as other weedy species (e.g. [63]). The presence and abundance of bacterial and fungal groups potentially underlie this context dependency in several instances. Environmental changes, such as exacerbated drought conditions, might modify suitability of oak woodland habitat and perhaps enhance invasibility of previously resistant systems. Therefore, given the complex relationships between aboveground and belowground biota [64], understanding the potential mechanisms mediating the association between invasive plant species and soil microorganisms could provide practical information for developing effective management strategies, as well as insight into the ecology of plant-soil food webs and diversity.

Supporting Information

S1 Fig. Sample-based phylotype accumulation curves for bacterial and fungal soil communities.

https://doi.org/10.1371/journal.pone.0163930.s001

(EPS)

S2 Fig. Shannon diversity for bacterial and fungal soil communities (A, B).

Differences between open grassland samples and oak woodland samples were statistically significant (ANOVA P < 0.05).

https://doi.org/10.1371/journal.pone.0163930.s002

(EPS)

S3 Fig. Significant differences (Wilcoxon test P < 0.01 after false discovery rate correction) in the abundance of bacterial and fungal classes between oak woodland soil samples (green) and open grassland soil samples (blue).

https://doi.org/10.1371/journal.pone.0163930.s003

(EPS)

S4 Fig. Abundance of fungal functional guilds between oak woodland soil samples (green) and open grassland soil samples (blue).

Note that the y axis is squared.

https://doi.org/10.1371/journal.pone.0163930.s004

(PDF)

S1 Table. Significant differences (Wilcoxon test P < 0.01 after false discovery rate correction) in the abundance of bacterial and fungal genera between oak woodland soil samples and open grassland soil samples.

https://doi.org/10.1371/journal.pone.0163930.s005

(DOCX)

Acknowledgments

We are thankful to Xavier Rojas for lab mentoring.

Author Contributions

  1. Conceptualization: ESG NF AB.
  2. Data curation: ESG AB.
  3. Formal analysis: AB.
  4. Funding acquisition: ESG NF AB.
  5. Investigation: ESG NF AB.
  6. Methodology: ESG NF AB.
  7. Project administration: ESG NF AB.
  8. Resources: NF AB.
  9. Software: NF AB.
  10. Supervision: ESG NF AB.
  11. Validation: NF AB.
  12. Visualization: AB.
  13. Writing – original draft: ESG AB.
  14. Writing – review & editing: ESG NF AB.

References

  1. 1. Grosholz E. Ecological and evolutionary consequences of coastal invasions. Trends Ecol. Evol. 2002; 17: 22–27.
  2. 2. Gaertner M, Biggs R, Te Beest M, Hui C, Molofsky J, Richardson DM. Invasive plants as drivers of regime shifts: identifying high-priority invaders that alter feedback relationships. Divers. Distrib. 2014; 20: 733–744.
  3. 3. Ehrenfeld JG, Kourtev P, Huang W. Changes in soil functions following invasions of exotic understory plants in deciduous forests. Ecol. Appl. 2001; 11: 1287–1300.
  4. 4. Wolfe BE, Klironomos JN. Breaking new ground: Soil communities and exotic plant invasion. Bioscience. 2005; 55: 477–487.
  5. 5. Reinhart KO, Callaway RM. Soil biota and invasive plants. New Phytol. 2006; 170: 445–457. pmid:16626467
  6. 6. Pringle A, Bever JD, Gardes NM, Oarrebt JL, Rillig MC, Klironomos JN. Mycorrhizal symbioses and plant invasions. Annu. Rev. Ecol. Evol. System. 2009; 40: 699–715.
  7. 7. van der Putten WH, Klironomos JN, Wardle DA. Microbial ecology of biological invasions. ISME J. 2007; 1: 28–37. pmid:18043611
  8. 8. Inderjit , van der Putten WH. Impacts of soil microbial communities on exotic plant invasions. Trends Ecol. Evol. 2010; 25: 512–519. pmid:20638747
  9. 9. Perkins LB, Nowak RS. Native and non-native grasses generate common types of plant-soil feedbacks by altering soil nutrients and microbial communities. Oikos. 2013; 122: 199–208.
  10. 10. Parepa M, Schaffner U, Bossdorf O. Help from under ground: soil biota facilitate knotweed invasion. Ecosphere 2013; 4: 31.
  11. 11. Eppinga MB, Rietkerk M, Dekker SC, De Ruiter PC, van der Putten W. Accumulation of local pathogens: a new hypothesis to explain exotic plant invasions. Oikos. 2006; 114: 168–176.
  12. 12. Bunn RA, Ramsey PW, Lekberg Y. Do native and invasive plants differ in their interactions with arbuscular mycorrhizal fungi? A meta-analysis. J Ecol. 2015; 103: 1547–1556.
  13. 13. Nafus AM, Davies KW. Medusahead ecology and management: California annual grasslands to the intermountain west. Invas. Plant Sci. Manag. 2014; 7: 210–221.
  14. 14. James JJ, Gornish ES, DiTomaso JM, Davy J, Doran M, Becchetti T, et al. Managing medusahead (Taeniatherum caput-medusae) on rangeland: A meta analysis of control effects and assessment of stakeholders needs. Range Ecol. Manag. 2015; 68: 215–223.
  15. 15. Trent JD, Young JA, Blank RR. Potential role of soil microorganisms in medusahead invasion. Proceedings—Ecology and management of annual rangelands. US. For. Serv. 1994; 140–142.
  16. 16. Blank RR, Sforza R, Morgan T. Shrublands under fire: Disturbance and recovery in a changing world. US. For. Serv. RMRS-P. 2008; 52: 73–76.
  17. 17. Belnap J, Phillips SL. Soil biota in an ungrazed grassland: Response to annual grass (Bromus tectorum) invasion. Ecol. Appl. 2001; 11: 1261–1275.
  18. 18. Querejeta JI, Egerton-Warburton LM, Allen MF. Topographic position modulates the mycorrhizal response of oak trees to interannual rainfall variability. Ecology. 2009; 90: 649–662. pmid:19341136
  19. 19. Gornish ES, James JJ, Laca EA. The value of oak woodland habitats as control for Medusahead (Taeniatherum caput-medusae). In: Proceedings of the Seventh California Oak Symposium: Managing Oak Woodlands in a Dynamic World, November 3–6, 2014, Visalia, CA. USDA For. Serv. Gen. Tech. Rept. PSW-GTR-251. 2015; 579 p.
    • 20. Stark JM, Firestone MK. Kinetic characteristics of ammonium-oxidizer communities in a California oak woodland-annual grassland. Soil Biochem. 1996; 28: 1307–1317.
    • 21. Waldrop MP, Firestone MK. Microbial community utilization of recalcitrant and simple carbon compounds: Impact of oak-woodland plant communities. Oecologia. 2004; 138: 275–284. pmid:14614618
    • 22. Dahlgren RA, Singer MJ, Huang X. Oak tree and grazing impacts on soil properties and nutrients in a California oak woodland. Biogeochemistry. 1997; 29: 45–64.
    • 23. Gornish ES, James JJ. Interactions among habitat, management, and demography for an invasive annual grass. Plant Ecol. 2016.
    • 24. Marañón T, Bartolome JW. Seed and seedling populations in two contrasting communities: open grassland and oak (Quercus agrifolia) understory in California. Oecologia. 1989; 10: 147–158.
    • 25. Shahin O, Martin-St Paul N, Rambal S, Joffre R, Richard F. Ectomycorrhizal fungal diversity in Quercus ilex Mediterranean woodlands: Variation among sites and over soil depth profiles in hyphal exploration types, species richness and community composition. Symbiosis. 2013; 61: 1–12.
    • 26. Fierer N, Leff JW, Adams BJ, Nielsen UN, Bates ST, Lauber CL, et al. Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. P Natl. Acad. Sci. USA. 2012; 109: 21390–21395. pmid:23236140
    • 27. McGuire KL, Payne SG, Palmer MI, Gillikin CM, Keefe D, Kim SJ et al. Digging the New York City skyline: soil fungal communities in green roofs and city parks. PLOS ONE. 2013; 8: e58020. pmid:23469260
    • 28. Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods. 2013; 10: 996–998. pmid:23955772
    • 29. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007; 73: 5261–5267. pmid:17586664
    • 30. McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2013; 6: 610–618. pmid:22134646
    • 31. Abarenkov K, Nilsson RH, Larsson KH, Alexander IJ, Eberhardt U, Erland S et al. The UNITE database for molecular identification of fungi—recent updates and future perspectives. New Phytol. 2010; 186: 281–285. pmid:20409185
    • 32. Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF, et al. Reagent and laboratory contamination can critically impact sequence-based microbiome analysis. BMC Biol. 2014; 12: 87. pmid:25387460
    • 33. Paulson JN, Stine OC, Bravo HC, Pop M. Differential abundance analysis for microbial marker-gene surveys. Nat. Methods. 2013; 10: 1200–1202. pmid:24076764
    • 34. Nguyen NH, Song Z, Bates ST, Branco S, Tedersoo L, Menke J, et al. FUNGuild: an open annotation tool for parsing fugal community datasets by ecological guild. Fungal Ecol. 2015;
    • 35. Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001; 26: 32–46.
    • 36. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Statist. Soc. B. 1995; 57: 289–300.
    • 37. Vestergard M, Ronn R, Ekelund F. Above-belowground interactions govern the course and impact of biological invasions. AOB Plants 2015; 7: plv025. pmid:25854693
    • 38. Shannon S, Flory SL, Reynolds H. Competitive context alters plant-soil feedback in an experimental woodland community. Oecologia. 2012; 169: 235–243. pmid:22101382
    • 39. Bansal S, Sheley RL. Annual grass invasion in sagebrush steppe: the relative importance of climate, soil properties and biotic interactions. Comm Ecol 2016; pmid:26920900
    • 40. Carey CJ, Beman JM, Eviner VT, Malmstrom C, Hart SC. Soil microbial community structure is unaltered by plant invasion, vegetation clipping, and nitrogen fertilization in experimental semi-arid grasslands. Front. Microb. 2012; 6: 466. pmid:26042104
    • 41. Kao-Kniffin J, Balser TC. Soil fertility and the impact of exotic invasion on microbial communities in Hawaiian forests. Microb. Ecol. 2008; 56: 55–63. pmid:17912579
    • 42. Prescott CE, Grayston SJ. Tree species influence on microbial communities in litter and soil: current knowledge and research needs. Forest Ecol. Manag. 2013; 309: 19–27.
    • 43. Kourtev PS, Ehrenfeld JG, Haggblom M. Exotic plant species alter the microbial community structure and function in the soil. Ecology. 2002; 83: 3152–3166.
    • 44. Potthoff M, Steenwerth K, Jackson LE, Drenovsky RE, Scow KM, Joergensen RG. Soil microbial community composition as affected by restoration practices in California grassland. Soil Biol. Biochem. 2006; 38: 1851–1860.
    • 45. Taylor MM, Hild AL, Shaw NL, Norton U, Collier TR. Plant recruitment and soil microbial characteristics of rehabilitation seedlings following wildfire in northern Utah. Restor. Ecol. 2014; 22: 598–607.
    • 46. Carini P, Mardsden PJ, Leff JW, Morgan EE, Strickland MS, Fierer N. Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. BioRxiv. 2016; http://dx.doi.org/10.1101/043372
    • 47. Blank RR, Sforza R. Plant-soil relationships of the invasive annual grass Taeniatherum caput-medusae: a reciprocal transplant experiment. Plant Soil 2007; 298: 7–19.
    • 48. Klironomos JN. Feedback with soil biota contributes to plant rarity and invasiveness in communities. Nature. 2002; 417: 67–70. pmid:11986666
    • 49. Mitchell CE, Power AG. Release of invasive plants from fungal and viral pathogens. Nature. 2003; 421: 625–627. pmid:12571594
    • 50. Grey WE, Quimby PC, Mathre DE, Young JA. Potential for biological control of downy brome (Bromus tectorum) and medusahead (Taeniatherum caput-medusae) with crown and root rot fungi. Weed Technol. 1995; 9: 362–265.
    • 51. Siegwart M, Bon MC, Widmer TL, Credpy N, Sforza R. First report of Fusarium arthrosporioides on medusaead (Taeniatherum caput-medusae) and preliminary tests for host-specificity. Plant Pathol. 2003; 52: 416.
    • 52. Allan E, van Ruijven J, Crawley MJ. Foliar fungal pathogens and grassland bibodiversity. Ecology 2010; 91: 2572–2582. pmid:20957952
    • 53. Olivera NL, Prieto L, Bertiller MB, Ferrero MA. Sheep grazing and soil bacterial diversity in shrublands of the Patagonian Monte, Argentina. J Arid Environ 2016; 125: 16–20.
    • 54. Franche C, Lindström K, Elmerich C. Nitrogen-fixing bacteria associated with leguminous and non-leguminous plants. Plant Soil 2009; 321: 35–59.
    • 55. Lehto T, Zwiazek JJ. Ecotomycorrhizas and water relations of trees: a review. Mycorr. 2011; 21: 71–90.
    • 56. Osono T. Ecology of ligninolytic fungi associated with leaf litter decomposition. Ecol. Res. 2007; 22: 955–974.
    • 57. Marler MJ, Zabinski CA, Callaway RM. Mycorrhizae indirectly enhance competitive effects of an invasive forb on a native bunchgrass. Ecology. 1999; 80: 1180–1186.
    • 58. Richardson DM, Allsopp N, D’Antonio CM, Milton SJ, Rejmanek M. Plant invasions—the role of mutualisms. Biol. Rev. 2000; 75: 65–93. pmid:10740893
    • 59. Batten KM, Scow KM, Davies KF, Harrison SP. Two invasive plants alter soil microbial community composition in serpentine grasslands. Biol. Invasions 2006; 8: 217–230.
    • 60. Lekberg Y, Gibbons SM, Rosendahl S, Ramsey PW. Severe plant invasions can increase mycorrhizal fungal abundance and diversity. ISME J. 2013; 7: 1424–1433. pmid:23486251
    • 61. Volgelsand KM, Bever JD. Mycorrhizal densities decline in association with nonnative plants and contribute to plant invasion. Ecology. 2009; 90: 399–407. pmid:19323224
    • 62. Stinson KA, Campbell SA, Powell JR, Wolfe BE, Callaway RM, Thelen GC, et al. Invasive plant suppresses the growth of native tree seedlings by disrupting belowground mutualisms. PLOS ONE 2006; 12: e1001817. pmid:16623597
    • 63. Kyser GB, Doran MP, McDougald NK, Orloff SB, Vargas RN, Wilson RG, et al. Site Characteristics Determine the Success of Prescribed Burning for Medusahead (Taeniatherum caput-medusae) Control. Invas. Plant Sci. Manag. 2008; 1: 376–384.
    • 64. Barberán A, McGuire KL, Wolf JA, Jones AF, Wright SJ, Turner BL, et al. Relating belowground microbial composition to the taxonomic, phylogenetic, and functional trait distributions of trees in a tropical forest. Ecol. Lett. 2015; 18: 1397–1405. pmid:26472095