Liyou Wu

Spatial and resource factors influencing high microbial diversity in soil

ABSTRACT To begin defining the key determinants that drive microbial community structure in soil, we examined 29 soil samples from four geographically distinct locations taken from the surface, vadose zone, and saturated subsurface using a small-subunit rRNA-based cloning approach. While microbial communities in low-carbon, saturated, subsurface soils showed dominance, microbial communities in low-carbon surface soils showed remarkably uniform distributions, and all species were equally abundant. Two diversity indices, the reciprocal of Simpson’s index (1/D) and the log series index, effectively distinguished between the dominant and uniform diversity patterns. For example, the uniform profiles characteristic of the surface communities had diversity index values that were 2 to 3 orders of magnitude greater than those for the high-dominance, saturated, subsurface communities. In a site richer in organic carbon, microbial communities consistently exhibited the uniform distribution pattern regardless of soil water content and depth. The uniform distribution implies that competition does not shape the structure of these microbial communities. Theoretical studies based on mathematical modeling suggested that spatial isolation could limit competition in surface soils, thereby supporting the high diversity and a uniform community structure. Carbon resource heterogeneity may explain the uniform diversity patterns observed in the high-carbon samples even in the saturated zone. Very high levels of chromium contamination (e.g., >20%) in the high-organic-matter soils did not greatly reduce the diversity. Understanding mechanisms that may control community structure, such as spatial isolation, has important implications for preservation of biodiversity, management of microbial communities for bioremediation, biocontrol of root diseases, and improved soil fertility.

GeoChip: a comprehensive microarray for investigating biogeochemical, ecological and environmental processes

Owing to their vast diversity and as-yet uncultivated status, detection, characterization and quantification of microorganisms in natural settings are very challenging, and linking microbial diversity to ecosystem processes and functions is even more difficult. Microarray-based genomic technology for detecting functional genes and processes has a great promise of overcoming such obstacles. Here, a novel comprehensive microarray, termed GeoChip, has been developed, containing 24 243 oligonucleotide (50 mer) probes and covering >10 000 genes in >150 functional groups involved in nitrogen, carbon, sulfur and phosphorus cycling, metal reduction and resistance, and organic contaminant degradation. The developed GeoChip was successfully used for tracking the dynamics of metal-reducing bacteria and associated communities for an in situ bioremediation study. This is the first comprehensive microarray currently available for studying biogeochemical processes and functional activities of microbial communities important to human health, agriculture, energy, global climate change, ecosystem management, and environmental cleanup and restoration. It is particularly useful for providing direct linkages of microbial genes/populations to ecosystem processes and functions.

Development and Evaluation of Functional Gene Arrays for Detection of Selected Genes in the Environment

ABSTRACT To determine the potential of DNA array technology for assessing functional gene diversity and distribution, a prototype microarray was constructed with genes involved in nitrogen cycling: nitrite reductase (nirS and nirK) genes, ammonia mono-oxygenase (amoA) genes, and methane mono-oxygenase (pmoA) genes from pure cultures and those cloned from marine sediments. In experiments using glass slide microarrays, genes possessing less than 80 to 85% sequence identity were differentiated under hybridization conditions of high stringency (65°C). The detection limit fornirS genes was approximately 1 ng of pure genomic DNA and 25 ng of soil community DNA using our optimized protocol. A linear quantitative relationship (r2 = 0.89 to 0.94) was observed between signal intensity and target DNA concentration over a range of 1 to 100 ng for genomic DNA (or genomic DNA equivalent) from both pure cultures and mixed communities. However, the quantitative capacity of microarrays for measuring the relative abundance of targeted genes in complex environmental samples is less clear due to divergent target sequences. Sequence divergence and probe length affected hybridization signal intensity within a certain range of sequence identity and size, respectively. This prototype functional gene array did reveal differences in the apparent distribution ofnir and amoA and pmoA gene families in sediment and soil samples. Our results indicate that glass-based microarray hybridization has potential as a tool for revealing functional gene composition in natural microbial communities; however, more work is needed to improve sensitivity and quantitation and to understand the associated issue of specificity.

Nitrite Reductase Genes (nirK and nirS) as Functional Markers To Investigate Diversity of Denitrifying Bacteria in Pacific Northwest Marine Sediment Communities

ABSTRACT Genetic heterogeneity of denitrifying bacteria in sediment samples from Puget Sound and two sites on the Washington continental margin was studied by PCR approaches amplifying nirK andnirS genes. These structurally different but functionally equivalent single-copy genes coding for nitrite reductases, a key enzyme of the denitrification process, were used as a molecular marker for denitrifying bacteria. nirS sequences could be amplified from samples of both sampling sites, whereas nirKsequences were detected only in samples from the Washington margin. To assess the underlying nir gene structure, PCR products of both genes were cloned and screened by restriction fragment length polymorphism (RFLP). Rarefraction analysis revealed a high level of diversity especially for nirS clones from Puget Sound and a slightly lower level of diversity for nirK andnirS clones from the Washington margin. One group dominated within nirK clones, but no dominance and only a few redundant clones were seen between sediment samples fornirS clones in both habitats. Hybridization and sequencing confirmed that all but one of the 228 putative nirS clones were nirS with levels of nucleotide identities as low as 45.3%. Phylogenetic analysis grouped nirS clones into three distinct subclusters within the nirS gene tree which corresponded to the two habitats from which they were obtained. These sequences had little relationship to any strain with knownnirS sequences or to isolates (mostly close relatives ofPseudomonas stutzeri) from the Washington margin sediment samples. nirK clones were more closely related to each other than were the nirS clones, with 78.6% and higher nucleotide identities; clones showing only weak hybridization signals were not related to known nirK sequences. AllnirK clones were also grouped into a distinct cluster which could not be placed with any strain with known nirKsequences. These findings show a very high diversity of nirsequences within small samples and that these novel nirclusters, some very divergent from known sequences, are not known in cultivated denitrifiers.

Evaluation of PCR-Generated Chimeras, Mutations, and Heteroduplexes with 16S rRNA Gene-Based Cloning

ABSTRACT To evaluate PCR-generated artifacts (i.e., chimeras, mutations, and heteroduplexes) with the 16S ribosomal DNA (rDNA)-based cloning approach, a model community of four species was constructed from alpha, beta, and gamma subdivisions of the division Proteobacteriaas well as gram-positive bacterium, all of which could be distinguished by HhaI restriction digestion patterns. The overall PCR artifacts were significantly different among the three TaqDNA polymerases examined: 20% for Z-Taq, with the highest processitivity; 15% for LA-Taq, with the highest fidelity and intermediate processitivity; and 7% for the conventionally used DNA polymerase, AmpliTaq. In contrast to the theoretical prediction, the frequency of chimeras for both Z-Taq(8.7%) and LA-Taq (6.2%) was higher than that for AmpliTaq (2.5%). The frequencies of chimeras and of heteroduplexes for Z-Taq were almost three times higher than those of AmpliTaq. The total PCR artifacts increased as PCR cycles and template concentrations increased and decreased as elongation time increased. Generally the frequency of chimeras was lower than that of mutations but higher than that of heteroduplexes. The total PCR artifacts as well as the frequency of heteroduplexes increased as the species diversity increased. PCR artifacts were significantly reduced by using AmpliTaq and fewer PCR cycles (fewer than 20 cycles), and the heteroduplexes could be effectively removed from PCR products prior to cloning by polyacrylamide gel purification or T7 endonuclease I digestion. Based upon these results, an optimal approach is proposed to minimize PCR artifacts in 16S rDNA-based microbial community studies.

Simultaneous Recovery of RNA and DNA from Soils and Sediments

ABSTRACT Recovery of mRNA from environmental samples for measurement of in situ metabolic activities is a significant challenge. A robust, simple, rapid, and effective method was developed for simultaneous recovery of both RNA and DNA from soils of diverse composition by adapting our previous grinding-based cell lysis method (Zhou et al., Appl. Environ. Microbiol. 62:316–322, 1996) for DNA extraction. One of the key differences is that the samples are ground in a denaturing solution at a temperature below 0°C to inactivate nuclease activity. Two different methods were evaluated for separating RNA from DNA. Among the methods examined for RNA purification, anion exchange resin gave the best results in terms of RNA integrity, yield, and purity. With the optimized protocol, intact RNA and high-molecular-weight DNA were simultaneously recovered from 19 soil and stream sediment samples of diverse composition. The RNA yield from these samples ranged from 1.4 to 56 μg g of soil−1 dry weight), whereas the DNA yield ranged from 23 to 435 μg g−1. In addition, studies with the same soil sample showed that the DNA yield was, on average, 40% higher than that in our previous procedure and 68% higher than that in a commercial bead milling method. For the majority of the samples, the DNA and RNA recovered were of sufficient purity for nuclease digestion, microarray hybridization, and PCR or reverse transcription-PCR amplification.

Transcriptome dynamics of Deinococcus radiodurans recovering from ionizing radiation

Deinococcus radiodurans R1 (DEIRA) is a bacterium best known for its extreme resistance to the lethal effects of ionizing radiation, but the molecular mechanisms underlying this phenotype remain poorly understood. To define the repertoire of DEIRA genes responding to acute irradiation (15 kGy), transcriptome dynamics were examined in cells representing early, middle, and late phases of recovery by using DNA microarrays covering ≈94% of its predicted genes. At least at one time point during DEIRA recovery, 832 genes (28% of the genome) were induced and 451 genes (15%) were repressed 2-fold or more. The expression patterns of the majority of the induced genes resemble the previously characterized expression profile of recA after irradiation. DEIRA recA, which is central to genomic restoration after irradiation, is substantially up-regulated on DNA damage (early phase) and down-regulated before the onset of exponential growth (late phase). Many other genes were expressed later in recovery, displaying a growth-related pattern of induction. Genes induced in the early phase of recovery included those involved in DNA replication, repair, and recombination, cell wall metabolism, cellular transport, and many encoding uncharacterized proteins. Collectively, the microarray data suggest that DEIRA cells efficiently coordinate their recovery by a complex network, within which both DNA repair and metabolic functions play critical roles. Components of this network include a predicted distinct ATP-dependent DNA ligase and metabolic pathway switching that could prevent additional genomic damage elicited by metabolism-induced free radicals.

Detection of genes involved in biodegradation and biotransformation in microbial communities by using 50-mer oligonucleotide microarrays

ABSTRACT To effectively monitor biodegrading populations, a comprehensive 50-mer-based oligonucleotide microarray was developed based on most of the 2,402 known genes and pathways involved in biodegradation and metal resistance. This array contained 1,662 unique and group-specific probes with <85% similarity to their nontarget sequences. Based on artificial probes, our results showed that under hybridization conditions of 50°C and 50% formamide, the 50-mer microarray hybridization can differentiate sequences having <88% similarity. Specificity tests with representative pure cultures indicated that the designed probes on the arrays appeared to be specific to their corresponding target genes. The detection limit was ∼5 to 10 ng of genomic DNA in the absence of background DNA and 50 to 100 ng of pure-culture genomic DNA in the presence of background DNA or 1.3 × 107 cells in the presence of background RNA. Strong linear relationships between the signal intensity and the target DNA and RNA were observed (r2 = 0.95 to 0.99). Application of this type of microarray to analyze naphthalene-amended enrichment and soil microcosms demonstrated that microflora changed differently depending on the incubation conditions. While the naphthalene-degrading genes from Rhodococcus-type microorganisms were dominant in naphthalene-degrading enrichments, the genes involved in naphthalene (and polyaromatic hydrocarbon and nitrotoluene) degradation from gram-negative microorganisms, such as Ralstonia, Comamonas, and Burkholderia, were most abundant in the soil microcosms. In contrast to general conceptions, naphthalene-degrading genes from Pseudomonas were not detected, although Pseudomonas is widely known as a model microorganism for studying naphthalene degradation. The real-time PCR analysis with four representative genes showed that the microarray-based quantification was very consistent with real-time PCR (r2 = 0.74). In addition, application of the arrays to both polyaromatic-hydrocarbon- and benzene-toluene-ethylbenzene-xylene-contaminated and uncontaminated soils indicated that the developed microarrays appeared to be useful for profiling differences in microbial community structures. Our results indicate that this technology has potential as a specific, sensitive, and quantitative tool in revealing a comprehensive picture of the compositions of biodegradation genes and the microbial community in contaminated environments, although more work is needed to improve detection sensitivity.

Microarray-based analysis of subnanogram quantities of microbial community DNAs by using whole-community genome amplification.

ABSTRACT Microarray technology provides the opportunity to identify thousands of microbial genes or populations simultaneously, but low microbial biomass often prevents application of this technology to many natural microbial communities. We developed a whole-community genome amplification-assisted microarray detection approach based on multiple displacement amplification. The representativeness of amplification was evaluated using several types of microarrays and quantitative indexes. Representative detection of individual genes or genomes was obtained with 1 to 100 ng DNA from individual or mixed genomes, in equal or unequal abundance, and with 1 to 500 ng community DNAs from groundwater. Lower concentrations of DNA (as low as 10 fg) could be detected, but the lower template concentrations affected the representativeness of amplification. Robust quantitative detection was also observed by significant linear relationships between signal intensities and initial DNA concentrations ranging from (i) 0.04 to 125 ng (r2 = 0.65 to 0.99) for DNA from pure cultures as detected by whole-genome open reading frame arrays, (ii) 0.1 to 1,000 ng (r2 = 0.91) for genomic DNA using community genome arrays, and (iii) 0.01 to 250 ng (r2 = 0.96 to 0.98) for community DNAs from ethanol-amended groundwater using 50-mer functional gene arrays. This method allowed us to investigate the oligotrophic microbial communities in groundwater contaminated with uranium and other metals. The results indicated that microorganisms containing genes involved in contaminant degradation and immobilization are present in these communities, that their spatial distribution is heterogeneous, and that microbial diversity is greatly reduced in the highly contaminated environment.

Reproducibility and quantitation of amplicon sequencing-based detection

To determine the reproducibility and quantitation of the amplicon sequencing-based detection approach for analyzing microbial community structure, a total of 24 microbial communities from a long-term global change experimental site were examined. Genomic DNA obtained from each community was used to amplify 16S rRNA genes with two or three barcode tags as technical replicates in the presence of a small quantity (0.1% wt/wt) of genomic DNA from Shewanella oneidensis MR-1 as the control. The technical reproducibility of the amplicon sequencing-based detection approach is quite low, with an average operational taxonomic unit (OTU) overlap of 17.2%±2.3% between two technical replicates, and 8.2%±2.3% among three technical replicates, which is most likely due to problems associated with random sampling processes. Such variations in technical replicates could have substantial effects on estimating β-diversity but less on α-diversity. A high variation was also observed in the control across different samples (for example, 66.7-fold for the forward primer), suggesting that the amplicon sequencing-based detection approach could not be quantitative. In addition, various strategies were examined to improve the comparability of amplicon sequencing data, such as increasing biological replicates, and removing singleton sequences and less-representative OTUs across biological replicates. Finally, as expected, various statistical analyses with preprocessed experimental data revealed clear differences in the composition and structure of microbial communities between warming and non-warming, or between clipping and non-clipping. Taken together, these results suggest that amplicon sequencing-based detection is useful in analyzing microbial community structure even though it is not reproducible and quantitative. However, great caution should be taken in experimental design and data interpretation when the amplicon sequencing-based detection approach is used for quantitative analysis of the β-diversity of microbial communities.