Search

ABSTRACT Targeted amplicon sequencing is widely used in microbial ecology studies. However, sequencing artifacts and amplification biases are of great concern. To identify sources of these artifacts, a systematic analysis was performed using mock communities comprised of 16S rRNA genes from 33 bacterial strains. Our results indicated that while sequencing errors were generally isolated to low-abundance operational taxonomic units, chimeric sequences were a major source of artifacts. Singleton and doubleton sequences were primarily chimeras. Formation of chimeric sequences was significantly correlated with the GC content of the targeted sequences. Low-GC-content mock community members exhibited lower rates of chimeric sequence formation. GC content also had a large impact on sequence recovery. The quantitative capacity was notably limited, with substantial recovery variations and weak correlation between anticipated and observed strain abundances. The mock community strains with higher GC content had higher recovery rates than strains with lower GC content. Amplification bias was also observed due to the differences in primer affinity. A two-step PCR strategy reduced the number of chimeric sequences by half. In addition, comparative analyses based on the mock communities showed that several widely used sequence processing pipelines/methods, including DADA2, Deblur, UCLUST, UNOISE, and UPARSE, had different advantages and disadvantages in artifact removal and rare species detection. These results are important for improving sequencing quality and reliability and developing new algorithms to process targeted amplicon sequences.IMPORTANCEAmplicon sequencing of targeted genes is the predominant approach to estimate the membership and structure of microbial communities. However, accurate reconstruction of community composition is difficult due to sequencing errors, and other methodological biases and effective approaches to overcome these challenges are essential. Using a mock community of 33 phylogenetically diverse strains, this study evaluated the effect of GC content on sequencing results and tested different approaches to improve overall sequencing accuracy while characterizing the pros and cons of popular amplicon sequence data processing approaches. The sequencing results from this study can serve as a benchmarking data set for future algorithmic improvements. Furthermore, the new insights on sequencing error, chimera formation, and GC bias from this study will help enhance the quality of amplicon sequencing studies and support the development of new data analysis approaches.

ABSTRACT Sulfate-reducing bacteria (SRB) are essential functional microbial taxa for degrading organic matter (OM) in anoxic marine environments. However, there are little experimental data regarding how SRB regulates microbial communities. Here, we applied a top-down microbial community management approach by inhibiting SRB to elucidate their contributions to the microbial community during OM degradation. Based on the highly replicated microcosms (n = 20) of five different incubation stages, we found that many microbial community properties were influenced after inhibiting SRB, including the composition, structure, network, and community assembly processes. We also found a strong coexistence pattern between SRB and other abundant phylogenetic lineages via positive frequency-dependent selection. The relative abundances of the families Synergistaceae, Peptostreptococcaceae, Dethiosulfatibacteraceae, Prolixibacteraceae, Marinilabiliaceae, and Marinifilaceae were simultaneously suppressed after inhibiting SRB during OM degradation. A close association between SRB and the order Marinilabiliales among coexisting taxa was most prominent. They contributed to preserved modules during network successions, were keystone nodes mediating the networked community, and contributed to homogeneous ecological selection. The molybdate tolerance test of the isolated strains of Marinilabiliales showed that inhibited SRB (not the inhibitor of SRB itself) triggered a decrease in the relative abundance of Marinilabiliales. We also found that inhibiting SRB resulted in reduced pH, which is unsuitable for the growth of most Marinilabiliales strains, while the addition of pH buffer (HEPES) in SRB-inhibited treatment microcosms restored the pH and the relative abundances of these bacteria. These data supported that SRB could modify niches to affect species coexistence. IMPORTANCE Our model offers insight into the ecological properties of SRB and identifies a previously undocumented dimension of OM degradation. This targeted inhibition approach could provide a novel framework for illustrating how functional microbial taxa associate the composition and structure of the microbial community, molecular ecological network, and community assembly processes. These findings emphasize the importance of SRB during OM degradation. Our results proved the feasibility of the proposed study framework, inhibiting functional taxa at the community level, for illustrating when and to what extent functional taxa can contribute to ecosystem services.

ABSTRACT Land conversion for intensive agriculture produces unfavorable changes to soil ecosystems, causing global concern. Soil bacterial communities mediate essential terrestrial ecosystem processes, making it imperative to understand their responses to agricultural perturbations. Here, we used high-throughput sequencing coupled with a functional gene array to study temporal dynamics of soil bacterial communities over 1 year under different disturbance intensities across a U.S. Southern Plains agroecosystem, including tallgrass prairie, Old World bluestem pasture, no-tillage (NT) canola, and conventional tillage (CT) wheat. Land use had the greatest impact on bacterial taxonomic diversity, whereas sampling time and its interaction with land use were central to functional diversity differences. The main drivers of taxonomic diversity were tillage > sampling time > temperature, while all measured factors explained similar amounts of variations in functional diversity. Temporal differences had the strongest correlation with total nitrogen > rainfall > nitrate. Within land uses, community variations for CT wheat were attributed to nitrogen levels, whereas soil organic matter and soil water content explained community variations for NT canola. In comparison, all measured factors contributed almost equally to variations in grassland bacterial communities. Finally, functional diversity had a stronger relationship with taxonomic diversity for CT wheat compared to phylogenetic diversity in the prairie. These findings reinforce that tillage management has the greatest impact on bacterial community diversity, with sampling time also critical. Hence, our study highlights the importance of the interaction between temporal dynamics and land use in influencing soil microbiomes, providing support for reducing agricultural disturbance to conserve soil biodiversity. IMPORTANCE Agricultural sustainability relies on healthy soils and microbial diversity. Agricultural management alters soil conditions and further influences the temporal dynamics of soil microbial communities essential to ecosystem functions, including organic matter dynamics, nutrient cycling, and plant nutrient availability. Yet, the responses to agricultural management are also dependent on soil type and climatic region, emphasizing the importance of assessing sustainability at local scales. To evaluate the impact of agricultural management practices, we examined bacterial communities across a management disturbance gradient over 1 year in a U.S. Southern Plains agroecosystem and determined that intensive management disturbance and sampling time critically impacted bacterial structural diversity, while their interactive effect influenced functional diversity and other soil health indicators. Overall, this study provides insights into how reducing soil disturbance can positively impact microbial community diversity and soil properties in the U.S. Southern Plains.

ABSTRACT Hexavalent chromium [Cr(VI)] is a common environmental pollutant. However, little is known about the genetic basis of microbial evolution under Cr(VI) stress and the influence of the prior evolution histories on the subsequent evolution under Cr(VI) stress. In this study, Desulfovibrio vulgaris Hildenborough (DvH), a model sulfate-reducing bacterium, was experimentally evolved for 600 generations. By evolving the replicate populations of three genetically diverse DvH clones, including ancestor (AN, without prior experimental evolution history), non-stress-evolved EC3-10, and salt stress-evolved ES9-11, the contributions of adaptation, chance, and pre-existing genetic divergence to the evolution under Cr(VI) stress were able to be dissected. Significantly decreased lag phases under Cr(VI) stress were observed in most evolved populations, while increased Cr(VI) reduction rates were primarily observed in populations evolved from EC3-10 and ES9-11. The pre-existing genetic divergence in the starting clones showed strong influences on the changes in lag phases, growth rates, and Cr(VI) reduction rates. Additionally, the genomic mutation spectra in populations evolved from different starting clones were significantly different. A total of 14 newly mutated genes obtained mutations in at least two evolved populations, suggesting their importance in Cr(VI) adaptation. An in-frame deletion mutation of one of these genes, the chromate transporter gene DVU0426, demonstrated that it played an important role in Cr(VI) tolerance. Overall, our study identified potential key functional genes for Cr(VI) tolerance and demonstrated the important role of pre-existing genetic divergence in evolution under Cr(VI) stress conditions. IMPORTANCE Chromium is one of the most common heavy metal pollutants of soil and groundwater. The potential of Desulfovibrio vulgaris Hildenborough in heavy metal bioremediation such as Cr(VI) reduction was reported previously; however, experimental evidence of key functional genes involved in Cr(VI) resistance are largely unknown. Given the genetic divergence of microbial populations in nature, knowledge on how this divergence affects the microbial adaptation to a new environment such as Cr(VI) stress is very limited. Taking advantage of our previous study, three groups of genetically diverse D. vulgaris Hildenborough populations with or without prior experimental evolution histories were propagated under Cr(VI) stress for 600 generations. Whole-population genome resequencing of the evolved populations revealed the genomic changes underlying the improved Cr(VI) tolerance. The strong influence of the pre-existing genetic divergence in the starting clones on evolution under Cr(VI) stress conditions was demonstrated at both phenotypic and genetic levels.

ABSTRACT While functional gene arrays (FGAs) have greatly expanded our understanding of complex microbial systems, specificity, sensitivity, and quantitation challenges remain. We developed a new generation of FGA, GeoChip 5.0, using the Agilent platform. Two formats were created, a smaller format (GeoChip 5.0S), primarily covering carbon-, nitrogen-, sulfur-, and phosphorus-cycling genes and others providing ecological services, and a larger format (GeoChip 5.0M) containing the functional categories involved in biogeochemical cycling of C, N, S, and P and various metals, stress response, microbial defense, electron transport, plant growth promotion, virulence, gyrB, and fungus-, protozoan-, and virus-specific genes. GeoChip 5.0M contains 161,961 oligonucleotide probes covering >365,000 genes of 1,447 gene families from broad, functionally divergent taxonomic groups, including bacteria (2,721 genera), archaea (101 genera), fungi (297 genera), protists (219 genera), and viruses (167 genera), mainly phages. Computational and experimental evaluation indicated that designed probes were highly specific and could detect as little as 0.05 ng of pure culture DNAs within a background of 1 μg community DNA (equivalent to 0.005% of the population). Additionally, strong quantitative linear relationships were observed between signal intensity and amount of pure genomic (∼99% of probes detected; r > 0.9) or soil (∼97%; r > 0.9) DNAs. Application of the GeoChip to a contaminated groundwater microbial community indicated that environmental contaminants (primarily heavy metals) had significant impacts on the biodiversity of the communities. This is the most comprehensive FGA to date, capable of directly linking microbial genes/populations to ecosystem functions. IMPORTANCE The rapid development of metagenomic technologies, including microarrays, over the past decade has greatly expanded our understanding of complex microbial systems. However, because of the ever-expanding number of novel microbial sequences discovered each year, developing a microarray that is representative of real microbial communities, is specific and sensitive, and provides quantitative information remains a challenge. The newly developed GeoChip 5.0 is the most comprehensive microarray available to date for examining the functional capabilities of microbial communities important to biogeochemistry, ecology, environmental sciences, and human health. The GeoChip 5 is highly specific, sensitive, and quantitative based on both computational and experimental assays. Use of the array on a contaminated groundwater sample provided novel insights on the impacts of environmental contaminants on groundwater microbial communities.

ABSTRACT Contamination from anthropogenic activities has significantly impacted Earth’s biosphere. However, knowledge about how environmental contamination affects the biodiversity of groundwater microbiomes and ecosystem functioning remains very limited. Here, we used a comprehensive functional gene array to analyze groundwater microbiomes from 69 wells at the Oak Ridge Field Research Center (Oak Ridge, TN), representing a wide pH range and uranium, nitrate, and other contaminants. We hypothesized that the functional diversity of groundwater microbiomes would decrease as environmental contamination (e.g., uranium or nitrate) increased or at low or high pH, while some specific populations capable of utilizing or resistant to those contaminants would increase, and thus, such key microbial functional genes and/or populations could be used to predict groundwater contamination and ecosystem functioning. Our results indicated that functional richness/diversity decreased as uranium (but not nitrate) increased in groundwater. In addition, about 5.9% of specific key functional populations targeted by a comprehensive functional gene array (GeoChip 5) increased significantly (P < 0.05) as uranium or nitrate increased, and their changes could be used to successfully predict uranium and nitrate contamination and ecosystem functioning. This study indicates great potential for using microbial functional genes to predict environmental contamination and ecosystem functioning. IMPORTANCE Disentangling the relationships between biodiversity and ecosystem functioning is an important but poorly understood topic in ecology. Predicting ecosystem functioning on the basis of biodiversity is even more difficult, particularly with microbial biomarkers. As an exploratory effort, this study used key microbial functional genes as biomarkers to provide predictive understanding of environmental contamination and ecosystem functioning. The results indicated that the overall functional gene richness/diversity decreased as uranium increased in groundwater, while specific key microbial guilds increased significantly as uranium or nitrate increased. These key microbial functional genes could be used to successfully predict environmental contamination and ecosystem functioning. This study represents a significant advance in using functional gene markers to predict the spatial distribution of environmental contaminants and ecosystem functioning toward predictive microbial ecology, which is an ultimate goal of microbial ecology.