Benefits of symbiotic ectomycorrhizal fungi to plant water relations depend on plant genotype in pinyon pine

Rhizosphere microbes, such as root-associated fungi, can improve plant access to soil resources, affecting plant health, productivity, and stress tolerance. While mycorrhizal associations are ubiquitous, plant-microbe interactions can be species specific. Here we show that the specificity of the effects of microbial symbionts on plant function can go beyond species level: colonization of roots by ectomycorrhizal fungi (EMF) of the genus Geopora has opposite effects on water uptake, and stomatal control of desiccation in drought tolerant and intolerant genotypes of pinyon pine (Pinus edulis Engelm.). These results demonstrate, for the first time, that microorganisms can have significant and opposite effects on important plant functional traits like stomatal control of desiccation that are associated with differential mortality and growth in nature. They also highlight that appropriate pairing of plant genotypes and microbial associates will be important for mitigating climate change impacts on vegetation.

Seasonal variation and potential roles of dark septate fungi in an arid grassland

High temperatures and extended drought in temperate and tropical arid ecosystems promote the colonization of diverse microenvironments by dark septate fungi (DSF). These fungi contribute to soil nutrient cycling, soil stabilization, and plant survival, but the roles of individual DSF species, their distributions, and their community diversity are poorly understood. The objective of this study was to evaluate the distribution, seasonal variation, and potential roles of DSF on plant growth. We collected biocrust (lichen-, moss-, and cyanobacterium-dominated biocrusts) soils at different depths and rhizosphere soils from two grasses, Bromus tectorum and Pleuraphis jamesii, in an arid grassland near Moab, Utah, USA. Seasonal variation of DSF was evaluated using culture-based approaches and compared with fungal community profiles from next-generation sequencing (NGS). Culturing showed that DSF were 30% more abundant in biocrusts compared with the focal rhizospheres. The abundance of DSF varied seasonally in belowground samples (rhizosphere and below-biocrust), with a significant increase during the summer months. Pleosporales was the dominant order (35%) in both biocrust and rhizosphere soils out of 817 isolated fungi. Dominant DSF genera in culture included Alternaria, Preussia, Cladosporium, Phoma, and an unknown Pleosporales. Similar results were observed in biocrust and rhizosphere soils NGS. Further, seed germination experiments using dominant taxa were conducted to determine their potential roles on germination and seedling growth using maize as a model plant. Cladosporium and unknown Pleosporales isolates showed plant growth-promoting ability. The variation in abundance of DSF, their differential occurrence in different microenvironments, and their ability to grow in a xerotolerant medium reflect adaptations to summer environmental conditions and to changes in the abundance of organic matter, as well as a potential increase in plant investment in these fungi when heat and drought stresses are more severe.

Resistance, Resilience, and Recovery of Dryland Soil Bacterial Communities Across Multiple Disturbances

Dryland ecosystems are sensitive to perturbations and generally slow to recover post disturbance. The microorganisms residing in dryland soils are especially important as they contribute to soil structure and nutrient cycling. Disturbance can have particularly strong effects on dryland soil structure and function, yet the natural resistance and recovery of the microbial components of dryland soils has not been well documented. In this study, the recovery of surface soil bacterial communities from multiple physical and environmental disturbances is assessed. Samples were collected from three field sites in the vicinity of Moab, UT, United States, 6 to 7 years after physical and climate disturbance manipulations had been terminated, allowing for the assessment of community recovery. Additionally, samples were collected in a transect that included three habitat patches: the canopy zone soils under the dominant shrubs, the interspace soils that are colonized by biological soil crusts, and edge soils at the plot borders. Field site and habitat patch were significant factors structuring the bacterial communities, illustrating that sites and habitats harbored unique soil microbiomes. Across the different sites and disturbance treatments, there was evidence of significant bacterial community recovery, as bacterial biomass and diversity were not significantly different than control plots. There was, however, a small number of 16S rRNA gene amplicon sequence variants that distinguished particular treatments, suggesting that legacy effects of the disturbances still remained. Taken together, these data suggest that dryland bacterial communities may possess a previously unappreciated potential to recover within years of the original disturbance.

Ectomycorrhizal and Dark Septate Fungal Associations of Pinyon Pine Are Differentially Affected by Experimental Drought and Warming

Changing climates can cause shifts in temperature and precipitation, resulting in warming and drought in some regions. Although each of these factors has been shown to detrimentally affect forest ecosystems worldwide, information on the impacts of the combined effects of warming and drought is lacking. Forest trees rely on mutualistic root-associated fungi that contribute significantly to plant health and protection against climate stresses. We used a six-year, ecosystem-scale temperature and precipitation manipulation experiment targeted to simulate the climate in 2100 in the Southwestern United States to quantify the effects of drought, warming and combined drought and warming on the root colonization (abundance), species composition and diversity of ectomycorrhizal fungi (EMF), and dark septate fungal endophytes in a widespread woodland tree, pinyon pine (Pinus edulis E.). Our results show that pinyon shoot growth after 6 years of these treatments was reduced more by drought than warming. The combined drought and warming treatment reduced the abundance and diversity of EMF more than either treatment alone. Individual ectomycorrhizal fungal taxa, including the drought tolerant Cenococcum geophilum, were present in all treatments but the combined drought and warming treatment. The combined drought and warming treatment also reduced the abundance of dark septate endophytes (DSE), but did not affect their diversity or species composition. The current year shoot growth of the trees correlated positively with ectomycorrhizal fungal diversity, highlighting the importance of diversity in mutualistic relationships to plant growth. Our results suggest that EMF may be more important than DSE to aboveground growth in P. edulis, but also more susceptible to the negative effects of combined climate stressors.

Keratinophilic fungi: Specialized fungal communities in a desert ecosystem identified using cultured-based and Illumina sequencing approaches

Soil fungi in desert ecosystems are adapted to harsh environmental conditions such as high soil surface temperatures and limited organic matter and water. Given limited carbon inputs from plant material, heterotrophic fungi likely use unconventional sources of carbon in these systems. A baiting method was used to culture keratinophilic fungi from biocrust and rhizosphere soils in an arid grassland in Utah, USA. Fungi were baited using llama and sheep wool, horsehair, and snakeskin on two media, and pure cultures were identified using ITS and LSU rRNA sequences. One hundred-eighteen fungal colonies were grown, representing a total of 32 Operational Taxonomic Units (OTUs) at 97 % similarity. Cultures were dominated by the phylum Ascomycota (88 %) followed by Mucoromycota (8.6 %) and Basidiomycota (3.4 %). The orders Pleosporales, Eurotiales, Hypocreales, and Sordariales were commonly isolated, with the dominant taxa Alternaria (27 %), Aspergillus (22 %), Fusarium (11 %), and Chaetomium (8%). Thirty percent of the fungi isolated have the capacity to degrade keratin in vitro using a keratin azure assay, with Penicillium showing the highest degradation followed by Geomyces, Alternaria, and Fusarium. Although keratin degraders can be infectious, dermatophytes associated with skin infections were not isolated in culture or detected in Illumina sequencing. Illumina sequencing was used to determine general patterns in seasonal variation and habitat preference of keratinophiles. Alternaria was the most abundant genus with >70 % of the sequences. The combination of Illumina data with culture-dependent approaches facilitated the characterization of a specialized community and confirmed the low abundance of dermatophytes in this arid site.

Genomes and secretomes of Ascomycota fungi reveal diverse functions in plant biomass decomposition and pathogenesis

BACKGROUND

The dominant fungi in arid grasslands and shrublands are members of the Ascomycota phylum. Ascomycota fungi are important drivers in carbon and nitrogen cycling in arid ecosystems. These fungi play roles in soil stability, plant biomass decomposition, and endophytic interactions with plants. They may also form symbiotic associations with biocrust components or be latent saprotrophs or pathogens that live on plant tissues. However, their functional potential in arid soils, where organic matter, nutrients and water are very low or only periodically available, is poorly characterized.

RESULTS

Five Ascomycota fungi were isolated from different soil crust microhabitats and rhizosphere soils around the native bunchgrass Pleuraphis jamesii in an arid grassland near Moab, UT, USA. Putative genera were Coniochaeta, isolated from lichen biocrust, Embellisia from cyanobacteria biocrust, Chaetomium from below lichen biocrust, Phoma from a moss microhabitat, and Aspergillus from the soil. The fungi were grown in replicate cultures on different carbon sources (chitin, native bunchgrass or pine wood) relevant to plant biomass and soil carbon sources. Secretomes produced by the fungi on each substrate were characterized. Results demonstrate that these fungi likely interact with primary producers (biocrust or plants) by secreting a wide range of proteins that facilitate symbiotic associations. Each of the fungal isolates secreted enzymes that degrade plant biomass, small secreted effector proteins, and proteins involved in either beneficial plant interactions or virulence. Aspergillus and Phoma expressed more plant biomass degrading enzymes when grown in grass- and pine-containing cultures than in chitin. Coniochaeta and Embellisia expressed similar numbers of these enzymes under all conditions, while Chaetomium secreted more of these enzymes in grass-containing cultures.

CONCLUSIONS

This study of Ascomycota genomes and secretomes provides important insights about the lifestyles and the roles that Ascomycota fungi likely play in arid grassland, ecosystems. However, the exact nature of those interactions, whether any or all of the isolates are true endophytes, latent saprotrophs or opportunistic phytopathogens, will be the topic of future studies.

Presence and distribution of insect-associated and entomopathogenic fungi in a temperate pine forest soil: An integrated approach

For decades entomopathogenic fungi have garnered interest as possible alternatives to chemical pesticides. However, their ecology outside of agroecosystems demands further study. We assessed the diversity and abundance of entomopathogenic and insect-associated fungi at a loblolly pine forest in North Carolina, USA using culture-dependent and next-generation sequencing libraries. Fungi were isolated using Galleriamellonella larvae, as well as from soil dilutions plated on a selective medium. Isolates were identified using Sanger sequencing of the ITS and LSU rRNA gene regions, and represented 36 OTUs including Metarhizium, Lecanicillium, and Paecilomyces. Additionally, we assessed the chitinolytic potential of isolates and found widespread, variable ability to degrade chitin within and between genera. Phylogenetic analyses resolved several isolates to genus, with some forming clades with other insect-associated taxa, as well as with fungi associated with plant tissues. Saprophytes were widely distributed in soil, while entomopathogens were less abundant and present primarily in the top two cm of the soil. The similarity between culture-dependent and next-generation sequencing results demonstrates that both methods can be used concurrently in this system to study the ecology of entomopathogenic fungi.

Chronic Physical Disturbance Substantially Alters the Response of Biological Soil Crusts to a Wetting Pulse, as Characterized by Metatranscriptomic Sequencing

Biological soil crusts (biocrusts) are microbial communities that are a feature of arid surface soils worldwide. In drylands where precipitation is pulsed and ephemeral, the ability of biocrust microbiota to rapidly initiate metabolic activity is critical to their survival. Community gene expression was compared after a short duration (1 h) wetting pulse in both intact and soils disturbed by chronic foot trampling. Across the metatranscriptomes the majority of transcripts were cyanobacterial in origin, suggesting that cyanobacteria accounted for the bulk of the transcriptionally active cells. Chronic trampling substantially altered the functional profile of the metatranscriptomes, specifically resulting in a significant decrease in transcripts for nitrogen fixation. Soil depth (biocrust and below crust) was a relatively small factor in differentiating the metatranscriptomes, suggesting that the metabolically active bacteria were similar between shallow soil horizons. The dry samples were consistently enriched for hydrogenase genes, indicating that molecular hydrogen may serve as an energy source for the desiccated soil communities. The water pulse was associated with a restructuring of the metatranscriptome, particularly for the biocrusts. Biocrusts increased transcripts for photosynthesis and carbon fixation, suggesting a rapid resuscitation upon wetting. In contrast, the trampled surface soils showed a much smaller response to wetting, indicating that trampling altered the metabolic response of the community. Finally, several biogeochemical cycling genes in carbon and nitrogen cycling were assessed for their change in abundance due to wetting in the biocrusts. Different transcripts encoding the same gene product did not show a consensus response, with some more abundant in dry or wet biocrusts, highlighting the challenges in relating transcript abundance to biogeochemical cycling rates. These observations demonstrate that metatranscriptome sequencing was able to distinguish alterations in the function of arid soil microbial communities at two varying temporal scales, a long-term ecosystems disturbance through foot trampling, and a short term wetting pulse. Thus, community metatranscriptomes have the potential to inform studies on the response and resilience of biocrusts to various environmental perturbations.

Short-Term Transcriptional Response of Microbial Communities to Nitrogen Fertilization in a Pine Forest Soil

Numerous studies have examined the long-term effect of experimental nitrogen (N) deposition in terrestrial ecosystems; however, N-specific mechanistic markers are difficult to disentangle from responses to other environmental changes. The strongest picture of N-responsive mechanistic markers is likely to arise from measurements over a short (hours to days) time scale immediately after inorganic N deposition. Therefore, we assessed the short-term (3-day) transcriptional response of microbial communities in two soil strata from a pine forest to a high dose of N fertilization (ca. 1 mg/g of soil material) in laboratory microcosms. We hypothesized that N fertilization would repress the expression of fungal and bacterial genes linked to N mining from plant litter. However, despite N suppression of microbial respiration, the most pronounced differences in functional gene expression were between strata rather than in response to the N addition. Overall, ∼4% of metabolic genes changed in expression with N addition, while three times as many (∼12%) were significantly different across the different soil strata in the microcosms. In particular, we found little evidence of N changing expression levels of metabolic genes associated with complex carbohydrate degradation (CAZymes) or inorganic N utilization. This suggests that direct N repression of microbial functional gene expression is not the principle mechanism for reduced soil respiration immediately after N deposition. Instead, changes in expression with N addition occurred primarily in general cell maintenance areas, for example, in ribosome-related transcripts. Transcriptional changes in functional gene abundance in response to N addition observed in longer-term field studies likely result from changes in microbial composition.IMPORTANCE Ecosystems are receiving increased nitrogen (N) from anthropogenic sources, including fertilizers and emissions from factories and automobiles. High levels of N change ecosystem functioning. For example, high inorganic N decreases the microbial decomposition of plant litter, potentially reducing nutrient recycling for plant growth. Understanding how N regulates microbial decomposition can improve the prediction of ecosystem functioning over extended time scales. We found little support for the conventional view that high N supply represses the expression of genes involved in decomposition or alters the expression of bacterial genes for inorganic N cycling. Instead, our study of pine forest soil 3 days after N addition showed changes in microbial gene expression related to cell maintenance and stress response. This highlights the challenge of establishing predictive links between microbial gene expression levels and measures of ecosystem function.

Shared features of cryptic plasmids from environmental and pathogenic Francisella species

The Francisella genus includes several recognized species, additional potential species, and other representatives that inhabit a range of incredibly diverse ecological niches, but are not closely related to the named species. Francisella species have been obtained from a wide variety of clinical and environmental sources; documented species include highly virulent human and animal pathogens, fish pathogens, opportunistic human pathogens, tick endosymbionts, and free-living isolates inhabiting brackish water. While more than 120 Francisella genomes have been sequenced to date, only a few contain plasmids, and most of these appear to be cryptic, with unknown benefit to the host cell. We have identified several putative cryptic plasmids in the sequenced genomes of three Francisella novicida and F. novicida-like strains (TX07-6608, AZ06-7470, DPG_3A-IS) and two new Francisella species (F. frigiditurris CA97-1460 and F. opportunistica MA06-7296). These plasmids were compared to each other and to previously identified plasmids from other Francisella species. Some of the plasmids encoded functions potentially involved in replication, conjugal transfer and partitioning, environmental survival (transcriptional regulation, signaling, metabolism), and hypothetical proteins with no assignable functions. Genomic and phylogenetic comparisons of these new plasmids to the other known Francisella plasmids revealed some similarities that add to our understanding of the evolutionary relationships among the diverse Francisella species.

Bifiguratus adelaidae, gen. et sp. nov., a new member of Mucoromycotina in endophytic and soil-dwelling habitats

Illumina amplicon sequencing of soil in a temperate pine forest in the southeastern United States detected an abundant, nitrogen (N)-responsive fungal genotype of unknown phylogenetic affiliation. Two isolates with ribosomal sequences consistent with that genotype were subsequently obtained. Examination of records in GenBank revealed that a genetically similar fungus had been isolated previously as an endophyte of moss in a pine forest in the southwestern United States. The three isolates were characterized using morphological, genomic, and multilocus molecular data (18S, internal transcribed spacer [ITS], and 28S rRNA sequences). Phylogenetic and maximum likelihood phylogenomic reconstructions revealed that the taxon represents a novel lineage in Mucoromycotina, only preceded by Calcarisporiella, the earliest diverging lineage in the subphylum. Sequences for the novel taxon are frequently detected in environmental sequencing studies, and it is currently part of UNITE's dynamic list of most wanted fungi. The fungus is dimorphic, grows best at room temperature, and is associated with a wide variety of bacteria. Here, a new monotypic genus, Bifiguratus, is proposed, typified by Bifiguratus adelaidae.

Ribosomal RNA gene detection and targeted culture of novel nitrogen-responsive fungal taxa from temperate pine forest soil

Soil fungal communities are responsible for carbon and nitrogen (N) cycling. The high complexity of the soil fungal community and the high proportion of taxonomically unidentifiable sequences confound ecological interpretations in field studies because physiological information is lacking for many organisms known only by their rRNA sequences. This situation forces experimental comparisons to be made at broader taxonomic racks where functions become difficult to infer. The objective of this study was to determine OTU (operational taxonomic units) level responses of the soil fungal community to N enrichment in a temperate pine forest experiment and to use the sequencing data to guide culture efforts of novel N-responsive fungal taxa. Replicate samples from four soil horizons (up to 10 cm depth) were obtained from ambient, enriched CO₂ and N-fertilization plots. Through a fungal large subunit rRNA gene (LSU) sequencing survey, we identified two novel fungal clades that were abundant in our soil sampling (representing up to 27% of the sequences in some samples) and responsive to changes in soil N. The two N-responsive taxa with no predicted taxonomic association were targeted for isolation and culturing from specific soil samples where their sequences were abundant. Representatives of both OTUs were successfully cultured using a filtration approach. One taxon (OTU6) was most closely related to Saccharomycotina; the second taxon (OTU69) was most closely related to Mucoromycotina. Both taxa likely represent novel species. This study shows how observation of specific OTUs level responses to altered N status in a large rRNA gene field survey provided the impetus to design targeted culture approaches for isolation of novel N-responsive fungal taxa.

Assembly of Active Bacterial and Fungal Communities Along a Natural Environmental Gradient

Dormancy is thought to promote biodiversity within microbial communities, but how assembly of the active community responds to changes in environmental conditions is unclear. To measure the active and dormant communities of bacteria and fungi colonizing decomposing litter in maple forests, we targeted ribosomal genes and transcripts across a natural environmental gradient. Within bacterial and fungal communities, the active and dormant communities were phylogenetically distinct, but patterns of phylogenetic clustering varied. For bacteria, active communities were significantly more clustered than dormant communities, while the reverse was found for fungi. The proportion of operational taxonomic units (OTUs) classified as active and the degree of phylogenetic clustering of the active bacterial communities declined with increasing pH and decreasing C/N. No significant correlations were found for the fungal community. The opposing pattern of phylogenetic clustering in dormant and active communities and the differential response of active communities to environmental gradients suggest that dormancy differentially structures bacterial and fungal communities.

A new fungal large subunit ribosomal RNA primer for high-throughput sequencing surveys

The inclusion of phylogenetic metrics in community ecology has provided insights into important ecological processes, particularly when combined with high-throughput sequencing methods; however, these approaches have not been widely used in studies of fungal communities relative to other microbial groups. Two obstacles have been considered: (1) the internal transcribed spacer (ITS) region has limited utility for constructing phylogenies and (2) most PCR primers that target the large subunit (LSU) ribosomal unit generate amplicons that exceed current limits of high-throughput sequencing platforms. We designed and tested a PCR primer (LR22R) to target approximately 300-400 bp region of the D2 hypervariable region of the fungal LSU for use with the Illumina MiSeq platform. Both in silico and empirical analyses showed that the LR22R-LR3 pair captured a broad range of fungal taxonomic groups with a small fraction of non-fungal groups. Phylogenetic placement of publically available LSU D2 sequences showed broad agreement with taxonomic classification. Comparisons of the LSU D2 and the ITS2 ribosomal regions from environmental samples and known communities showed similar discriminatory abilities of the two primer sets. Together, these findings show that the LR22R-LR3 primer pair has utility for phylogenetic analyses of fungal communities using high-throughput sequencing methods.

Soil bacterial and fungal community responses to nitrogen addition across soil depth and microhabitat in an arid shrubland

Arid shrublands are stressful environments, typified by alkaline soils low in organic matter, with biologically-limiting extremes in water availability, temperature, and UV radiation. The widely-spaced plants and interspace biological soil crusts in these regions provide soil nutrients in a localized fashion, creating a mosaic pattern of plant- or crust-associated microhabitats with distinct nutrient composition. With sporadic and limited rainfall, nutrients are primarily retained in the shallow surface soil, patterning biological activity. We examined soil bacterial and fungal community responses to simulated nitrogen (N) deposition in an arid Larrea tridentata-Ambrosia dumosa field experiment in southern Nevada, USA, using high-throughput sequencing of ribosomal RNA genes. To examine potential interactions among the N application, microhabitat and soil depth, we sampled soils associated with shrub canopies and interspace biological crusts at two soil depths (0-0.5 or 0-10 cm) across the N-amendment gradient (0, 7, and 15 kg ha(-1) yr(-1)). We hypothesized that localized compositional differences in soil microbiota would constrain the impacts of N addition to a microhabitat distribution that would reflect highly localized geochemical conditions and microbial community composition. The richness and community composition of both bacterial and fungal communities differed significantly by microhabitat and with soil depth in each microhabitat. Only bacterial communities exhibited significant responses to the N addition. Community composition correlated with microhabitat and depth differences in soil geochemical features. Given the distinct roles of soil bacteria and fungi in major nutrient cycles, the resilience of fungi and sensitivity of bacteria to N amendments suggests that increased N input predicted for many arid ecosystems could shift nutrient cycling toward pathways driven primarily by fungal communities.

Soil microbial responses to nitrogen addition in arid ecosystems

The N cycle of arid ecosystems is influenced by low soil organic matter, high soil pH, and extremes in water potential and temperature that lead to open canopies and development of biological soil crusts (biocrusts). We investigated the effects of N amendment on soil microbial dynamics in a Larrea tridentata-Ambrosia dumosa shrubland site in southern Nevada USA. Sites were fertilized with a NO₃-NH₄ mix at 0, 7, and 15 kg N ha^-1 y^-1 from March 2012 to March 2013. In March 2013, biocrust (0–0.5 cm) and bulk soils (0–10 cm) were collected beneath Ambrosia canopies and in the interspaces between plants. Biomass responses were assessed as bacterial and fungal SSU rRNA gene copy number and chlorophyll a concentration. Metabolic responses were measured by five ecoenzyme activities and rates of N transformation. By most measures, nutrient availability, microbial biomass, and process rates were greater in soils beneath the shrub canopy compared to the interspace between plants, and greater in the surface biocrust horizon compared to the deeper 10 cm soil profile. Most measures responded positively to experimental N addition. Effect sizes were generally greater for bulk soil than biocrust. Results were incorporated into a meta-analysis of arid ecosystem responses to N amendment that included data from 14 other studies. Effect sizes were calculated for biomass and metabolic responses. Regressions of effect sizes, calculated for biomass, and metabolic responses, showed similar trends in relation to N application rate and N load (rate × duration). The critical points separating positive from negative treatment effects were 88 kg ha^-1 y^-1 and 159 kg ha^-1, respectively, for biomass, and 70 kg ha^-1 y^-1 and 114 kg ha^-1, respectively, for metabolism. These critical values are comparable to those for microbial biomass, decomposition rates and respiration reported in broader meta-analyses of N amendment effects in mesic ecosystems. However, large effect sizes at low N addition rates indicate that arid ecosystems are sensitive to modest increments in anthropogenic N deposition.

Forest floor community metatranscriptomes identify fungal and bacterial responses to N deposition in two maple forests

Anthropogenic N deposition alters patterns of C and N cycling in temperate forests, where forest floor litter decomposition is a key process mediated by a diverse community of bacteria and fungi. To track forest floor decomposer activity we generated metatranscriptomes that simultaneously surveyed the actively expressed bacterial and eukaryote genes in the forest floor, to compare the impact of N deposition on the decomposers in two natural maple forests in Michigan, USA, where replicate field plots had been amended with N for 16 years. Site and N amendment responses were compared using about 74,000 carbohydrate active enzyme transcript sequences (CAZymes) in each metatranscriptome. Parallel ribosomal RNA (rRNA) surveys of bacterial and fungal biomass and taxonomic composition showed no significant differences in either biomass or OTU richness between the two sites or in response to N. Site and N amendment were not significant variables defining bacterial taxonomic composition, but they were significant for fungal community composition, explaining 17 and 14% of the variability, respectively. The relative abundance of expressed bacterial and fungal CAZymes changed significantly with N amendment in one of the forests, and N-response trends were also identified in the second forest. Although the two ambient forests were similar in community biomass, taxonomic structure and active CAZyme profile, the shifts in active CAZyme profiles in response to N-amendment differed between the sites. One site responded with an over-expression of bacterial CAZymes, and the other site responded with an over-expression of both fungal and different bacterial CAZymes. Both sites showed reduced representation of fungal lignocellulose degrading enzymes in N-amendment plots. The metatranscriptome approach provided a holistic assessment of eukaryote and bacterial gene expression and is applicable to other systems where eukaryotes and bacteria interact.

Elevated atmospheric CO2 stimulates soil fungal diversity through increased fine root production in a semiarid shrubland ecosystem

Soil fungal communities are likely to be central in mediating microbial feedbacks to climate change through their effects on soil carbon (C) storage, nutrient cycling, and plant health. Plants often produce increased fine root biomass in response to elevated atmospheric carbon dioxide (CO2 ), but the responses of soil microbial communities are variable and uncertain, particularly in terms of species diversity. In this study, we describe the responses of the soil fungal community to free air CO2 enrichment (FACE) in a semiarid chaparral shrubland in Southern California (dominated by Adenomstoma fasciculatum) using large subunit rRNA gene sequencing. Community composition varied greatly over the landscape and responses to FACE were subtle, involving a few specific groups. Increased frequency of Sordariomycetes and Leotiomycetes, the latter including the Helotiales, a group that includes many dark septate endophytes known to associate positively with roots, was observed in the FACE plots. Fungal diversity, both in terms of richness and evenness, increased consistently in the FACE treatment, and was relatively high compared to other studies that used similar methods. Increases in diversity were observed across multiple phylogenetic levels, from genus to class, and were distributed broadly across fungal lineages. Diversity was also higher in samples collected close to (5 cm) plants compared to samples in canopy gaps (30 cm away from plants). Fungal biomass correlated well with soil organic matter (SOM) content, but patterns of diversity were correlated with fine root production rather than SOM. We conclude that the fungal community in this ecosystem is tightly linked to plant fine root production, and that future changes in the fungal community in response to elevated CO2 and other climatic changes will be primarily driven by changes in plant belowground allocation. Potential feedbacks mediated by soil fungi, such as soil C sequestration, nutrient cycling, and pathogenesis, are discussed.

Comparative analyses of a putative Francisella conjugative element

A large circular plasmid detected in Francisella novicida-like strain PA10-7858, designated pFNPA10, was sequenced completely and analyzed. This 41,013-bp plasmid showed no homology to any of the previously sequenced Francisella plasmids and was 8-10 times larger in size than them. A total of 57 ORFs were identified within pFNPA10 and at least 9 of them encoded putative proteins with homology to different conjugal transfer proteins. The presence of iteron-like direct repeats and an ORF encoding a putative replication protein within pFNPA10 suggested that it replicated by the theta mode. Phylogenetic analyses indicated that pFNPA10 had no near neighbors in the databases and that it may have originated within an environmental Francisella lineage. Based on its features, pFNPA10 appears to be a novel extra-chromosomal genetic element within the genus Francisella. The suitability of pFNPA10 as a vector for transformation of species of Francisella by conjugation remains to be explored.

Ribosomal Database Project: data and tools for high throughput rRNA analysis

Ribosomal Database Project (RDP; http://rdp.cme.msu.edu/) provides the research community with aligned and annotated rRNA gene sequence data, along with tools to allow researchers to analyze their own rRNA gene sequences in the RDP framework. RDP data and tools are utilized in fields as diverse as human health, microbial ecology, environmental microbiology, nucleic acid chemistry, taxonomy and phylogenetics. In addition to aligned and annotated collections of bacterial and archaeal small subunit rRNA genes, RDP now includes a collection of fungal large subunit rRNA genes. RDP tools, including Classifier and Aligner, have been updated to work with this new fungal collection. The use of high-throughput sequencing to characterize environmental microbial populations has exploded in the past several years, and as sequence technologies have improved, the sizes of environmental datasets have increased. With release 11, RDP is providing an expanded set of tools to facilitate analysis of high-throughput data, including both single-stranded and paired-end reads. In addition, most tools are now available as open source packages for download and local use by researchers with high-volume needs or who would like to develop custom analysis pipelines.

High bacterial diversity of biological soil crusts in water tracks over permafrost in the high arctic polar desert

In this study we report the bacterial diversity of biological soil crusts (biocrusts) inhabiting polar desert soils at the northern land limit of the Arctic polar region (83° 05 N). Employing pyrosequencing of bacterial 16S rRNA genes this study demonstrated that these biocrusts harbor diverse bacterial communities, often as diverse as temperate latitude communities. The effect of wetting pulses on the composition of communities was also determined by collecting samples from soils outside and inside of permafrost water tracks, hill slope flow paths that drain permafrost-affected soils. The intermittent flow regime in the water tracks was correlated with altered relative abundance of phylum level taxonomic bins in the bacterial communities, but the alterations varied between individual sampling sites. Bacteria related to the Cyanobacteria and Acidobacteria demonstrated shifts in relative abundance based on their location either inside or outside of the water tracks. Among cyanobacterial sequences, the proportion of sequences belonging to the family Oscillatoriales consistently increased in relative abundance in the samples from inside the water tracks compared to those outside. Acidobacteria showed responses to wetting pulses in the water tracks, increasing in abundance at one site and decreasing at the other two sites. Subdivision 4 acidobacterial sequences tended to follow the trends in the total Acidobacteria relative abundance, suggesting these organisms were largely responsible for the changes observed in the Acidobacteria. Taken together, these data suggest that the bacterial communities of these high latitude polar biocrusts are diverse but do not show a consensus response to intermittent flow in water tracks over high Arctic permafrost.

Changes in Fungal Community Composition in Response to Elevated Atmospheric CO2 and Nitrogen Fertilization Varies with Soil Horizon

Increasing levels of atmospheric carbon dioxide (CO₂) and rates of nitrogen (N)-deposition to forest ecosystems are predicted to alter the structure and function of soil fungal communities, but the spatially heterogeneous distribution of soil fungi has hampered investigations aimed at understanding such impacts. We hypothesized that soil physical and chemical properties and fungal community composition would be differentially impacted by elevated atmospheric CO₂ (eCO₂) and N-fertilization in spatially separated field samples, in the forest floor, 0–2, 2–5, and 5–10 cm depth intervals in a loblolly pine Free-Air Carbon Dioxide Enrichment (FACE) experiment. In all soils, quantitative PCR-based estimates of fungal biomass were highest in the forest floor. Fungal richness, based on pyrosequencing of the fungal ribosomal large subunit gene, increased in response to N-fertilization in 0–2 cm and forest floor intervals. Composition shifted in forest floor, 0–2 and 2–5 cm intervals in response to N-fertilization, but the shift was most distinct in the 0–2 cm interval, in which the largest number of statistically significant changes in soil chemical parameters (i.e., phosphorus, organic matter, calcium, pH) was also observed. In the 0–2 cm interval, increased recovery of sequences from the Thelephoraceae, Tricholomataceae, Hypocreaceae, Clavicipitaceae, and Herpotrichiellaceae families and decreased recovery of sequences from the Amanitaceae correlated with N-fertilization. In this same depth interval, Amanitaceae, Tricholomataceae, and Herpotriciellaceae sequences were recovered less frequently from soils exposed to eCO₂ relative to ambient conditions. These results demonstrated that vertical stratification should be taken into consideration in future efforts to elucidate environmental impacts on fungal communities and their feedbacks on ecosystem processes.

Response of biological soil crust diazotrophs to season, altered summer precipitation, and year-round increased temperature in an arid grassland of the colorado plateau, USA

Biological soil crusts (biocrusts), which supply significant amounts of fixed nitrogen into terrestrial ecosystems worldwide (∼33 Tg y(-1)), are likely to respond to changes in temperature and precipitation associated with climate change. Using nifH gene-based surveys, we explored variation in the diazotrophic community of biocrusts of the Colorado Plateau, USA in response to season (autumn vs. spring), as well as field manipulations that increased the frequency of small volume precipitation events and year-round soil temperature. Abundance of nifH genes in biocrusts ranged from 3 × 10(6) to 1 × 10(8) g(-1) soil, and nifH from heterocystous cyanobacteria closely related to Scytonema hyalinum, Spirirestis rafaelensis, and Nostoc commune comprised >98% of the total. Although there was no apparent seasonal effect on total nifH gene abundance in the biocrusts, T-RFLP analysis revealed a strong seasonal pattern in nifH composition. SpirirestisnifH abundance was estimated to oscillate 1 to >2 orders of magnitude between autumn (low) and spring (high). A year-round increase of soil temperature (2-3°C) had little effect on the diazotroph community structure over 2 years. Altered summer precipitation had little impact on diazotroph community structure over the first 1.5 years of the study, when natural background patterns across years and seasons superseded any treatment effects. However, after the second summer of treatments, nifH abundance was 2.6-fold lower in biocrusts receiving altered precipitation. Heterocystous cyanobacteria were apparently more resilient to altered precipitation than other cyanobacteria. The results demonstrate that diazotrophic community composition of biocrusts in this semi-arid grassland undergoes strong seasonal shifts and that the abundance of its dominant members decreased in response to more frequent, small volume precipitation events.

Rapid phylogenetic and functional classification of short genomic fragments with signature peptides

BACKGROUND

Classification is difficult for shotgun metagenomics data from environments such as soils, where the diversity of sequences is high and where reference sequences from close relatives may not exist. Approaches based on sequence-similarity scores must deal with the confounding effects that inheritance and functional pressures exert on the relation between scores and phylogenetic distance, while approaches based on sequence alignment and tree-building are typically limited to a small fraction of gene families. We describe an approach based on finding one or more exact matches between a read and a precomputed set of peptide 10-mers.

RESULTS

At even the largest phylogenetic distances, thousands of 10-mer peptide exact matches can be found between pairs of bacterial genomes. Genes that share one or more peptide 10-mers typically have high reciprocal BLAST scores. Among a set of 403 representative bacterial genomes, some 20 million 10-mer peptides were found to be shared. We assign each of these peptides as a signature of a particular node in a phylogenetic reference tree based on the RNA polymerase genes. We classify the phylogeny of a genomic fragment (e.g., read) at the most specific node on the reference tree that is consistent with the phylogeny of observed signature peptides it contains. Using both synthetic data from four newly-sequenced soil-bacterium genomes and ten real soil metagenomics data sets, we demonstrate a sensitivity and specificity comparable to that of the MEGAN metagenomics analysis package using BLASTX against the NR database. Phylogenetic and functional similarity metrics applied to real metagenomics data indicates a signal-to-noise ratio of approximately 400 for distinguishing among environments. Our method assigns ~6.6 Gbp/hr on a single CPU, compared with 25 kbp/hr for methods based on BLASTX against the NR database.

CONCLUSIONS

Classification by exact matching against a precomputed list of signature peptides provides comparable results to existing techniques for reads longer than about 300 bp and does not degrade severely with shorter reads. Orders of magnitude faster than existing methods, the approach is suitable now for inclusion in analysis pipelines and appears to be extensible in several different directions.

Genetic diversity within the genus Francisella as revealed by comparative analyses of the genomes of two North American isolates from environmental sources

Background

Francisella tularensis is an intracellular pathogen that causes tularemia in humans and the public health importance of this bacterium has been well documented in recent history. Francisella philomiragia, a distant relative of F. tularensis, is thought to constitute an environmental lineage along with Francisella novicida. Nevertheless, both F. philomiragia and F. novicida have been associated with human disease, primarily in immune-compromised individuals. To understand the genetic relationships and evolutionary contexts among different lineages within the genus Francisella, the genome of Francisella spp. strain TX07-7308 was sequenced and compared to the genomes of F. philomiragia strains ATCC 25017 and 25015, F. novicida strain U112, and F. tularensis strain Schu S4.

Results

The size of strain ATCC 25017 chromosome was 2,045,775 bp and contained 1,983 protein-coding genes. The size of strain TX07-7308 chromosome was 2,035,931 bp and contained 1,980 protein-coding genes. Pairwise BLAST comparisons indicated that strains TX07-7308 and ATCC 25017 contained 1,700 protein coding genes in common. NUCmer analyses revealed that the chromosomes of strains TX07-7308 and ATCC 25017 were mostly collinear except for a few gaps, translocations, and/or inversions. Using the genome sequence data and comparative analyses with other members of the genus Francisella (e.g., F. novicida strain U112 and F. tularensis strain Schu S4), several strain-specific genes were identified. Strains TX07-7308 and ATCC 25017 contained an operon with six open reading frames encoding proteins related to enzymes involved in thiamine biosynthesis that was absent in F. novicida strain U112 and F. tularensis strain Schu S4. Strain ATCC 25017 contained an operon putatively involved in lactose metabolism that was absent in strain TX07-7308, F. novicida strain U112, and F. tularensis strain Schu S4. In contrast, strain TX07-7308 contained an operon putatively involved in glucuronate metabolism that was absent in the genomes of strain ATCC 25017, F. novicida strain U112, and F. tularensis strain Schu S4. The polymorphic nature of polysaccharide biosynthesis/modification gene clusters among different Francisella strains was also evident from genome analyses.

Conclusions

From genome comparisons, it appeared that genes encoding novel functions have contributed to the metabolic enrichment of the environmental lineages within the genus Francisella. The inability to acquire new genes coupled with the loss of ancestral traits and the consequent reductive evolution may be a cause for, as well as an effect of, niche selection of F. tularensis. Sequencing and comparison of the genomes of more isolates are required to obtain further insights into the ecology and evolution of different species within the genus Francisella.

A robust PCR primer design platform applied to the detection of Acidobacteria Group 1 in soil

Environmental biosurveillance and microbial ecology studies use PCR-based assays to detect and quantify microbial taxa and gene sequences within a complex background of microorganisms. However, the fragmentary nature and growing quantity of DNA-sequence data make group-specific assay design challenging. We solved this problem by developing a software platform that enables PCR-assay design at an unprecedented scale. As a demonstration, we developed quantitative PCR assays for a globally widespread, ecologically important bacterial group in soil, Acidobacteria Group 1. A total of 33,684 Acidobacteria 16S rRNA gene sequences were used for assay design. Following 1 week of computation on a 376-core cluster, 83 assays were obtained. We validated the specificity of the top three assays, collectively predicted to detect 42% of the Acidobacteria Group 1 sequences, by PCR amplification and sequencing of DNA from soil. Based on previous analyses of 16S rRNA gene sequencing, Acidobacteria Group 1 species were expected to decrease in response to elevated atmospheric CO(2). Quantitative PCR results, using the Acidobacteria Group 1-specific PCR assays, confirmed the expected decrease and provided higher statistical confidence than the 16S rRNA gene-sequencing data. These results demonstrate a powerful capacity to address previously intractable assay design challenges.

Targeted and shotgun metagenomic approaches provide different descriptions of dryland soil microbial communities in a manipulated field study

The extent to which different sequence-based approaches describe environmental microbial communities in comparative studies is an important consideration when deriving inferences from ecological studies. The ability of a targeted metagenomic approach [small subunit (SSU) rRNA pyrosequencing] and shotgun metagenome approaches were compared to identify distinguishing features in dryland soil microbial communities from two different habitats: biological soil crusts (biocrusts) and creosote bush root zones. A parallel comparison was conducted to determine the ability of each approach to detect community differences potentially arising from a more subtle experimental treatment, long-term elevated atmospheric carbon dioxide. As expected, the biocrust datasets were clearly differentiated from root zone datasets using either of the sequencing approaches. However, the composition described by each approach was significantly different. The magnitude of comparative differences due to habitat or elevated CO2 treatment was larger with pyrosequenced SSU datasets or SSU reads recruited from shotgun metagenomes, than from SEED-classified shotgun metagenome reads. Finally, based on prior knowledge of the biocrust communities, the SSU-based datasets more accurately identified the dominant biocrust cyanobacteria sequences compared to the shotgun metagenome datasets.

Common bacterial responses in six ecosystems exposed to 10 years of elevated atmospheric carbon dioxide

Six terrestrial ecosystems in the USA were exposed to elevated atmospheric CO(2) in single or multifactorial experiments for more than a decade to assess potential impacts. We retrospectively assessed soil bacterial community responses in all six-field experiments and found ecosystem-specific and common patterns of soil bacterial community response to elevated CO(2) . Soil bacterial composition differed greatly across the six ecosystems. No common effect of elevated atmospheric CO(2) on bacterial biomass, richness and community composition across all of the ecosystems was identified, although significant responses were detected in individual ecosystems. The most striking common trend across the sites was a decrease of up to 3.5-fold in the relative abundance of Acidobacteria Group 1 bacteria in soils exposed to elevated CO(2) or other climate factors. The Acidobacteria Group 1 response observed in exploratory 16S rRNA gene clone library surveys was validated in one ecosystem by 100-fold deeper sequencing and semi-quantitative PCR assays. Collectively, the 16S rRNA gene sequencing approach revealed influences of elevated CO(2) on multiple ecosystems. Although few common trends across the ecosystems were detected in the small surveys, the trends may be harbingers of more substantive changes in less abundant, more sensitive taxa that can only be detected by deeper surveys. Representative bacterial 16S rRNA gene clone sequences were deposited in GenBank with Accession No. JQ366086–JQ387568.

Comparative assessment of fungal cellobiohydrolase I richness and composition in cDNA generated using oligo(dT) primers or random hexamers

Understanding soil fungal distribution and activities, particularly at the level of gene expression, is important in unveiling mechanisms regulating their activities in situ. Recent identification of fungal genes involved in carbon cycling has provided the foundation for developing reverse-transcriptase PCR assays to monitor spatiotemporal gene expression patterns in soils and other complex microbial systems. The polyadenylated 3' ends of eukaryotic mRNA transcripts enables the use of oligo(dT) primers for cDNA synthesis, but this can result in the overrepresentation of the 3' end of transcripts in cDNA pools. In an effort to increase the uniformity of transcripts represented in cDNA pools, random hexamers have been used. The use of both priming methods is abundant in the literature, but we do not know how these methods perform relative to each other. We performed comparative richness and compositional analyses of the fungal glycosyl hydrolase family 7 cellobiohydrolase I gene cbhI amplified from soil cDNAs that had been generated using either oligo(dT) primers or random hexamers. Our results demonstrate that similar cbhI richness and composition were recovered using both approaches. Richness estimates and compositional profiles of cbhI sequence libraries generated from random hexamer-primed cDNA were more variable than from libraries generated from oligo(dT) primed cDNA. However, our overall results indicate that, on average, comparable richness and composition were recovered from soil cDNAs when either priming method was used.

Response and resilience of soil biocrust bacterial communities to chronic physical disturbance in arid shrublands

The impact of 10 years of annual foot trampling on soil biocrusts was examined in replicated field experiments at three cold desert sites of the Colorado Plateau, USA. Trampling detrimentally impacted lichens and mosses, and the keystone cyanobacterium, Microcoleus vaginatus, resulting in increased soil erosion and reduced C and N concentrations in surface soils. Trampled biocrusts contained approximately half as much extractable DNA and 20-52% less chlorophyll a when compared with intact biocrusts at each site. Two of the three sites also showed a decline in scytonemin-containing, diazotrophic cyanobacteria in trampled biocrusts. 16S rRNA gene sequence and terminal restriction fragment length polymorphism (T-RFLP) analyses of soil bacteria from untrampled and trampled biocrusts demonstrated a reduced proportion (23-65% reduction) of M. vaginatus and other Cyanobacteria in trampled plots. In parallel, other soil bacterial species that are natural residents of biocrusts, specifically members of the Actinobacteria, Chloroflexi and Bacteroidetes, became more readily detected in trampled than in untrampled biocrusts. Replicate 16S rRNA T-RFLP profiles from trampled biocrusts at all three sites contained significantly more fragments (n = 17) than those of untrampled biocrusts (n≤6) and exhibited much higher variability among field replicates, indicating transition to an unstable disturbed state. Despite the dramatic negative impacts of trampling on biocrust physical structure and composition, M. vaginatus could still be detected in surface soils after 10 years of annual trampling, suggesting the potential for biocrust re-formation over time. Physical damage of biocrusts, in concert with changing temperature and precipitation patterns, has potential to alter performance of dryland ecosystems for decades.

Biological consequences of ancient gene acquisition and duplication in the large genome of Candidatus Solibacter usitatus Ellin6076

Members of the bacterial phylum Acidobacteria are widespread in soils and sediments worldwide, and are abundant in many soils. Acidobacteria are challenging to culture in vitro, and many basic features of their biology and functional roles in the soil have not been determined. Candidatus Solibacter usitatus strain Ellin6076 has a 9.9 Mb genome that is approximately 2–5 times as large as the other sequenced Acidobacteria genomes. Bacterial genome sizes typically range from 0.5 to 10 Mb and are influenced by gene duplication, horizontal gene transfer, gene loss and other evolutionary processes. Our comparative genome analyses indicate that the Ellin6076 large genome has arisen by horizontal gene transfer via ancient bacteriophage and/or plasmid-mediated transduction, and widespread small-scale gene duplications, resulting in an increased number of paralogs. Low amino acid sequence identities among functional group members, and lack of conserved gene order and orientation in regions containing similar groups of paralogs, suggest that most of the paralogs are not the result of recent duplication events. The genome sizes of additional cultured Acidobacteria strains were estimated using pulsed-field gel electrophoresis to determine the prevalence of the large genome trait within the phylum. Members of subdivision 3 had larger genomes than those of subdivision 1, but none were as large as the Ellin6076 genome. The large genome of Ellin6076 may not be typical of the phylum, and encodes traits that could provide a selective metabolic, defensive and regulatory advantage in the soil environment.

Reverse transcription-PCR methods significantly impact richness and composition measures of expressed fungal cellobiohydrolase I genes in soil and litter

The importance of soil fungi in complex carbon degradation and the recent identification of genes involved in this process have sparked considerable interest in examining fungal gene expression in situ. Expression of target eukaryotic genes is commonly examined using reverse transcription (RT)-PCR, during which single-stranded (ss) complementary DNA (cDNA) is synthesized from an oligo (dT) primer and the gene of interest is subsequently amplified by PCR using gene specific primers. Another method that is being increasingly employed in environmental gene expression studies is SMART PCR, which generates and amplifies double-stranded (ds) complementary DNA (cDNA) from sscDNA using PCR, prior to gene-specific PCR. We performed a replicated comparison of these two methods using RNA extracted from forest soil and litter to determine if the two approaches yielded comparable results. Richness, composition and reproducibility of gene expression profiles of the fungal glycosyl hydrolase family 7 (GH7) cellobiohydrolase I gene (cbhI) were examined when amplified from sscDNA or dscDNA synthesized using SMART PCR. In the dscDNA libraries from soil or litter samples, richness was significantly reduced and the composition was altered relative to sscDNA libraries. Library composition was significantly more reproducible among replicate sscDNA libraries than among parallel dscDNA libraries from litter. In sum, the reduced richness and altered composition produced in the dscDNA libraries could substantially influence ecological interpretations of the data. Defining the factors underpinning the methodological biases will potentially aid in optimizing the design of gene expression studies in soils and other complex environmental samples.

A comparison of taxon co-occurrence patterns for macro- and microorganisms

We examine co-occurrence patterns of microorganisms to evaluate community assembly "rules". We use methods previously applied to macroorganisms, both to evaluate their applicability to microorganisms and to allow comparison of co-occurrence patterns observed in microorganisms to those found in macroorganisms. We use a null model analysis of 124 incidence matrices from microbial communities, including bacteria, archaea, fungi, and algae, and we compare these results to previously published findings from a meta-analysis of almost 100 macroorganism data sets. We show that assemblages of microorganisms demonstrate nonrandom patterns of co-occurrence that are broadly similar to those found in assemblages of macroorganisms. These results suggest that some taxon co-occurrence patterns may be general characteristics of communities of organisms from all domains of life. We also find that co-occurrence in microbial communities does not vary among taxonomic groups or habitat types. However, we find that the degree of co-occurrence does vary among studies that use different methods to survey microbial communities. Finally, we discuss the potential effects of the undersampling of microbial communities on our results, as well as processes that may contribute to nonrandom patterns of co-occurrence in both macrobial and microbial communities such as competition, habitat filtering, historical effects, and neutral processes.

Three distinct clades of cultured heterocystous cyanobacteria constitute the dominant N2-fixing members of biological soil crusts of the Colorado Plateau, USA

The identity of the numerically dominant N(2)-fixing bacteria in biological soil crusts of the Colorado Plateau region and two outlying areas was determined using multiple approaches, to link the environmental diversity of nifH gene sequences to cultured bacterial isolates from the regions. Of the nifH sequence-types detected in soil crusts of the Colorado Plateau, 89% (421/473) were most closely related to nifH signature sequences from cyanobacteria of the order Nostocales. N(2)-fixing cyanobacterial strains were cultured from crusts and their morphotypes, 16S rRNA gene and nifH gene sequences were characterized. The numerically dominant diazotrophs in the Colorado Plateau crusts fell within three clades of heterocystous cyanobacteria. Two clades are well-represented by phylogenetically and morphologically coherent strains, corresponding to the descriptions of Nostoc commune and Scytonema hyalinum, which are widely recognized as important N(2)-fixing components of soil crusts. A third, previously-overlooked clade was represented by a phylogenetically coherent but morphologically diverse group of strains that encompass the morphogenera Tolypothrix and Spirirestis. Many of the strains in each of these groups contained at least two nifH copies that represent different clusters in the nifH environmental survey.

Acidobacteria phylum sequences in uranium-contaminated subsurface sediments greatly expand the known diversity within the phylum

The abundance and composition of bacteria of the phylum Acidobacteria were surveyed in subsurface sediments from uranium-contaminated sites using amplification of 16S rRNA genes followed by clone/sequence analysis. Analysis of sequences from this study and public databases produced a revised and greatly expanded phylogeny of the Acidobacteria phylum consisting of 26 subgroups.

Environmental survey for four pathogenic bacteria and closely related species using phylogenetic and functional genes

Bacterial species with high DNA sequence similarity to pathogens could affect the specificity of assays designed to detect biological threat agents in environmental samples. The natural presence of four pathogenic bacteria, Bacillus anthracis, Clostridium perfringens, Francisella tularensis, and Yersinia pestis and their closely related species, was determined for a large collection of soil and aerosol samples. Polymerase chain reaction (PCR) and gene sequencing were used using group-specific 16S rRNA primers to identify pathogens and related species, and pathogen-specific virulence genes. Close relatives of B. anthracis (B. cereus group species) were detected in 37% of the soils and 25% of the aerosol samples. The B. anthracis protective antigen (pag) gene or a close homolog was detected in 16 of these samples. For the other three pathogen groups, the frequency of detection was much lower, and none of the samples were positive with both the phylogenetic and virulence gene primer sets.

Composition of bacillus species in aerosols from 11 U.S. cities

A PCR-based heteroduplex assay was used to determine the presence and composition of Bacillus species in 11,059 Environmental Protection Agency PM2.5 aerosol samples from 11 U.S. cities. The assay differentiated three groups: Type A containing Bacillus anthracis and very closely related, often pathogenic, Bacillus cereus and Bacillus thuringiensis strains; Type B containing other B. cereus and B. thuringiensis strains; and a third group of more-distantly related Bacillus species. Eight of the 11 cities were positive for Bacillus species in 50% or more of the samples, and the percent of aerosol samples that contained the HD Type A group ranged from 3% to 32%. Cities from the eastern half of the United States generally contained a higher frequency and broader diversity of Bacillus species than the western half of the United States. Positive samples were detected throughout the year. These results have implications for pathogen detection in environmental samples, understanding the natural evolution of new pathogenic strains, and incidence of infection caused by strains of the B. cereus subgroup.

Current and emerging technologies for the study of bacteria in the outdoor air

Pathogenic bacteria dispersed into the outdoor air from natural events or human activity are an important concern affecting public health, agriculture, ecological health and international security. Recent studies in atmospheric microbiology have contributed to our ability to detect bacterial pathogens in air samples, increased our knowledge of the spatial and temporal dispersal of aerosolized bacteria and furthered our understanding of the natural bacterial diversity and composition of air. Significant questions remain about bacterial dispersal and diversity in outdoor air that will benefit from continued technological advances. These include determining spatial and temporal patterns of aerosol dispersal, identifying the links between detection, concentration and viability of bacterial pathogens and disease incidence, and documenting the composition and regional movement of bacterial species in air.

Microbial biogeography: putting microorganisms on the map

We review the biogeography of microorganisms in light of the biogeography of macroorganisms. A large body of research supports the idea that free-living microbial taxa exhibit biogeographic patterns. Current evidence confirms that, as proposed by the Baas-Becking hypothesis, 'the environment selects' and is, in part, responsible for spatial variation in microbial diversity. However, recent studies also dispute the idea that 'everything is everywhere'. We also consider how the processes that generate and maintain biogeographic patterns in macroorganisms could operate in the microbial world.

Detection of diverse new Francisella-like bacteria in environmental samples

Following detection of putative Francisella species in aerosol samples from Houston, Texas, we surveyed soil and water samples from the area for the agent of tularemia, Francisella tularensis, and related species. The initial survey used 16S rRNA gene primers to detect Francisella species and related organisms by PCR amplification of DNA extracts from environmental samples. This analysis indicated that sequences related to Francisella were present in one water and seven soil samples. This is the first report of the detection of Francisella-related species in soil samples by DNA-based methods. Cloning and sequencing of PCR products indicated the presence of a wide variety of Francisella-related species. Sequences from two soil samples were 99.9% similar to previously reported sequences from F. tularensis isolates and may represent new subspecies. Additional analyses with primer sets developed for detection and differentiation of F. tularensis subspecies support the finding of very close relatives to known F. tularensis strains in some samples. While the pathogenicity of these organisms is unknown, they have the potential to be detected in F. tularensis-specific assays. Similarly, a potential new subspecies of Francisella philomiragia was identified. The majority of sequences obtained, while more similar to those of Francisella than to any other genus, were phylogenetically distinct from known species and formed several new clades potentially representing new species or genera. The results of this study revise our understanding of the diversity and distribution of Francisella and have implications for tularemia epidemiology and our ability to detect bioterrorist activities.