Applications of functional gene microarrays for profiling microbial communities

Unlocking secrets of microbial ecotoxicology: recent achievements and future challenges

Environmental pollution is one of the main challenges faced by humanity. By their ubiquity and vast range of metabolic capabilities, microorganisms are affected by pollution with consequences on their host organisms and on the functioning of their environment. They also play key roles in the fate of pollutants through the degradation, transformation, and transfer of organic or inorganic compounds. Thus, they are crucial for the development of nature-based solutions to reduce pollution and of bio-based solutions for environmental risk assessment of chemicals. At the intersection between microbial ecology, toxicology, and biogeochemistry, microbial ecotoxicology is a fast-expanding research area aiming to decipher the interactions between pollutants and microorganisms. This perspective paper gives an overview of the main research challenges identified by the Ecotoxicomic network within the emerging One Health framework and in the light of ongoing interest in biological approaches to environmental remediation and of the current state of the art in microbial ecology. We highlight prevailing knowledge gaps and pitfalls in exploring complex interactions among microorganisms and their environment in the context of chemical pollution and pinpoint areas of research where future efforts are needed.

Oyster culture changed the phosphorus speciation in sediments through biodeposition

Phosphorus speciation in the sediments is regulated by a series of physicochemical and microbial processes, and directly affects water phosphorus pool. However, the influence of culture activities and microbial metabolism on the sedimentary phosphorus speciation is poorly studied. In this study, we compared the abundance of distinguishable phosphorus phases and other physicochemical properties of sediments from oyster-farming areas and reference areas. The Geochip 5.0 technique was introduced to reveal the microbiological mechanisms of phosphorus metabolic alteration. The results showed that oyster culture enhanced the bioavailability of phosphorus in sediments. The free organic phosphorus was reduced significantly, whereas the free inorganic phosphorus and iron-bound phosphorus greatly increased in the oyster culture area (P ≤ 0.05). Moreover, the results of Geochip showed that the oyster culture reshaped the microbial network structure in sediments, with typical phosphate-solubilizing and phosphorus-accumulating microbes being enriched by 17.76% and 10.60%. The abundance of functional genes related to the main phosphorus cycle pathways were also significantly increased (P ≤ 0.05) in the culture area compared to the reference area. This work suggested that oyster culture can greatly improve the microbial phosphorus metabolism and provided insights into the environmental recovery and reconstruction from marine aquaculture activities.

Oyster Biodeposition Alleviates Sediment Nutrient Overload: A Case Study at Shenzhen Bay, China

Oysters are ecological engineers, and previous studies have examined their role as competent facilitators of ecological restoration. However, the decisive role of oysters in the aquatic environment is still debatable because oyster biodeposition (OBD) may also increase the nutrients enriched in sediments. In order to better interpret this problem, we sampled sediment cores from representative oyster culture areas and uncultured areas in Shenzhen Bay. The results have shown that the TOC (total organic carbon) and TN (total nitrogen) decreased significantly (p < 0.05) at the surface sediment layer (0-20-cm deep) and the sediment layer (20-40-cm deep) of the oyster site compared with the reference site. The decreased TOC and TN were also observed at 60- to 100-cm sediment depth in the oyster site. This indicated that the OBD significantly impacted the concentration of TOC and TN in the sediment. To confirm the alleviative role of OBD, we conducted stable isotope (δ¹³C and δ¹⁵N) analyses, which further demonstrated the presence of heavier and less lighter forms of organic carbon and nitrogen sediment. The surface sediment layer (0-20 cm) at the oyster site showed 8% more δ¹³C‰ compared with the control site (p < 0.05), reflecting the reduction in the TOC. In order to reveal the potential microbial mechanisms involved in OBD, we performed a functional analysis using the Geochip5 advanced microarray technology. Regarding carbon metabolism, we observed that genes (encoding pullulanase, glucoamylase, exoglucanase, cellobiase, and xylanase) involved in the degradation of relatively labile C-based molecules (e.g., starch, cellulose, and hemicellulose) were highly represented in an experimental area (p < 0.05). In addition, microbes in the experimental area exhibited a greater capacity for degrading recalcitrant C (e.g., lignin), which involves glyoxal oxidase (glx), manganese peroxidase (mnp), and phenol oxidase. Among the genes controlling nitrogen metabolism, the genes involved in denitrification, assimilation, ammonification, and nitrification were differentially expressed compared with the control area. These results indicated that microbial metabolic roles might have enhanced the C/N-flux speed and reduced the overall nutrient status. We concluded that OBD alleviates sediment nutrient overload under oyster farming from a microbial ecological perspective in a rapidly urbanized coastal area.

Polystyrene microplastics alter the intestinal microbiota function and the hepatic metabolism status in marine medaka (Oryzias melastigma)

To assess the potential effects of microplastics (MPs) on gut microbiome, a simple investigation of gut microbial structure is not sufficient, and the function and association of gut microbial structure with host health should also be taken into account. Here, the effects of two particle sizes (2 and 200 μm) of polystyrene MPs (PS-MPs) on the gut microbiota of medaka were evaluated following oral administration at 0.3 and 3.0 μg/mg for 28 days. No change in body length and gut histopathology damage were observed. However, the exposure to PS-MPs significantly decreased fish body weight and disrupted the liver anti-oxidative status. The PS-MPs caused a shift in the gut microbial structure of medaka accompanied by changes in community function, including significant environmental stress, increased carbon degradation/fixation activities, and partially modified nitrogen/phosphorus/sulfur metabolic abilities. Furthermore, the PS-MPs exposure disturbed the glycolipid/tyrosine/energy metabolism and the endocrine balance. A potential correlation between the gut microecology and host response to PS-MPs exposure was also observed. These results indicated that the PS-MPs may contribute to gut-liver axis disruption, which could be the underlying toxicological mechanisms of PS-MPs exposure. This work has improved our knowledge about the relationship between gut microbiota dysbiosis and host metabolic disorders following MPs exposure.

Microbial functional genes commonly respond to elevated carbon dioxide

Atmospheric CO₂ concentration is increasing, largely due to anthropogenic activities. Previous studies of individual free-air CO₂ enrichment (FACE) experimental sites have shown significant impacts of elevated CO₂ (eCO₂) on soil microbial communities; however, no common microbial response patterns have yet emerged, challenging our ability to predict ecosystem functioning and sustainability in the future eCO₂ environment. Here we analyzed 66 soil microbial communities from five FACE sites, and showed common microbial response patterns to eCO₂, especially for key functional genes involved in carbon and nitrogen fixation (e.g., pcc/acc for carbon fixation, nifH for nitrogen fixation), carbon decomposition (e.g., amyA and pulA for labile carbon decomposition, mnp and lcc for recalcitrant carbon decomposition), and greenhouse gas emissions (e.g., mcrA for methane production, norB for nitrous oxide production) across five FACE sites. Also, the relative abundance of those key genes was generally increased and directionally associated with increased biomass, soil carbon decomposition, and soil moisture. In addition, a further literature survey of more disparate FACE experimental sites indicated increased biomass, soil carbon decay, nitrogen fixation, methane and nitrous oxide emissions, plant and soil carbon and nitrogen under eCO₂. A conceptual framework was developed to link commonly responsive functional genes with ecosystem processes, such as pcc/acc vs. soil carbon storage, amyA/pulA/mnp/lcc vs. soil carbon decomposition, and nifH vs. nitrogen availability, suggesting that such common responses of microbial functional genes may have the potential to predict ecosystem functioning and sustainability in the future eCO₂ environment.

Advances in monitoring soil microbial community dynamic and function

Microorganisms are vital to the overall ecosystem functioning, stability, and sustainability. Soil fertility and health depend on chemical composition and also on the qualitative and quantitative nature of microorganisms inhabiting it. Historically, denaturing gradient gel electrophoresis (DGGE) and temperature gradient gel electrophoresis (TGGE), single-strand conformation polymorphism, DNA amplification fingerprinting, amplified ribosomal DNA restriction analysis, terminal restriction fragment length polymorphism, length heterogeneity PCR, and ribosomal intergenic spacer analysis were used to assess soil microbial community structure (SMCS), abundance, and diversity. However, these methods had significant shortcomings and limitations for application in land reclamation monitoring. SMCS has been primarily determined by phospholipid fatty acid (PLFA) analysis. This method provides a direct measure of viable biomass in addition to a biochemical profile of the microbial community. PLFA has limitations such as overlap in the composition of microorganisms and the specificity of PLFAs signature. In recent years, high-throughput next-generation sequencing has dramatically increased the resolution and detectable spectrum of diverse microbial phylotypes from environmental samples and it plays a significant role in microbial ecology studies. Next-generation sequencings using 454, Illumina, SOLiD, and Ion Torrent platforms are rapid and flexible. The two methods, PLFA and next-generation sequencing, are useful in detecting changes in microbial community diversity and structure in different ecosystems. Single-molecule real-time (SMRT) and nanopore sequencing technologies represent third-generation sequencing (TGS) platforms that have been developed to address the shortcomings of second-generation sequencing (SGS). Enzymatic and soil respiration analyses are performed to further determine soil quality and microbial activities. Other valuable methods that are being recently applied to microbial function and structures include NanoSIM, GeoChip, and DNA stable staple isotope probing (DNA-SIP) technologies. They are powerful metagenomics tool for analyzing microbial communities, including their structure, metabolic potential, diversity, and their impact on ecosystem functions. This review is a critical analysis of current methods used in monitoring soil microbial community dynamic and functions.

Functional Gene Array-Based Ultrasensitive and Quantitative Detection of Microbial Populations in Complex Communities

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.

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.

Microbial functional gene patterns related to soil greenhouse gas emissions in oil contaminated areas

Linking microbial community structure to physiology and ecological processes is a critical focus of microbial ecology. To understand the microbial functional gene patterns related to soil greenhouse gas [carbon dioxide (CO₂), methane (CH₄) and nitrous oxide (N₂O)] emissions under oil contamination, we used functional gene array (GeoChip 5.0) analysis and network methods to investigate the feedback responses of soil microbial functional gene patterns and identify keystone genes in Shengli Oilfield, China. The microbial functional gene number, relative abundance and diversity involved in carbon degradation and nitrogen cycling decreased consistently with the reduced CO₂ and N₂O flux in oil contaminated soils, whereas the gene number and relative abundance of methane-production related genes increased with contamination. Functional molecular ecological networks were built based on random matrix theory, where network structures and properties showed significantly variation between oil contaminated and uncontaminated soils (P<0.05). Network nodes, connectivity and complexity all reduced under oil contamination. The sensitive and the highest connective genes in the network were identified as keystone genes, based on Mann-Whitney U tests and network analysis. Our findings improved the understanding of the microbe-mediated mechanisms affecting soil greenhouse gas emissions.

Microbial Functional Gene Diversity Predicts Groundwater Contamination and Ecosystem Functioning

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.

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.

Elucidation of the methanogenic potential from coalbed microbial communities amended with volatile fatty acids

The potential for modern coalfield methanogenesis was assessed using formation water from the Illinois Basin, Powder River Basin and Cook Inlet gas field as inocula for nutrient-replete incubations amended with C1-C5 fatty acids as presumed intermediates formed during anaerobic coal biodegradation. Instead of the expected rapid mineralization of these substrates, methanogenesis was inordinately slow (∼1 μmol day-1), following long lag periods (>100 days), and methane yields typically did not reach stoichiometrically expected levels. However, a gene microarray confirmed the potential for a wide variety of microbiological functions, including methanogenesis, at all sites. The Cook Inlet incubations produced methane at a relatively rapid rate when amended with butyrate (r = 0.98; p = 0.001) or valerate (r = 0.84; p = 0.04), a result that significantly correlated with the number of positive mcr gene sequence probes from the functional gene microarray and was consistent with the in situ detection of C4-C5 alkanoic acids. This finding highlighted the role of syntrophy for the biodegradation of the softer lignite and subbituminous coal in this formation, but methanogenesis from the harder subbituminous and bituminous coals in the other fields was less apparent. We conclude that coal methanogenesis is probably not limited by the inherent lack of metabolic potential, the presence of alternate electron acceptors or the lack of available nutrients, but more likely restricted by the inherent recalcitrance of the coal itself.

The shifts of sediment microbial community phylogenetic and functional structures during chromium (VI) reduction

The Lanzhou reach of the Yellow River, located at the upstream of Lanzhou, has been contaminated by heavy metals and polycyclic aromatic hydrocarbons over a long-time. We hypothesized that indigenous microbial communities would remediate those contaminants and some unique populations could play an important role in this process. In this study, we investigated the sediment microbial community structure and function from the Lanzhou reach. Sediment samples were collected from two nearby sites (site A and site B) in the Lanzhou reach along the Yellow River. Sediment geochemical property data showed that site A sediment samples contained significantly (p < 0.05) higher heavy metals than site B, such as chromium (Cr), manganese (Mn), and copper (Cu). Both site A and B samples were incubated with or without hexavalent chromium (Cr (VI)) for 30 days in the laboratory, and Cr (VI) reduction was only observed in site A sediment samples. After incubation, MiSeq sequencing of 16S rRNA gene amplicons revealed that the phylogenetic composition and structure of microbial communities changed in both samples, and especially Proteobacteria, as the most abundant phylum increased from 45.1 % to 68.2 % in site A, and 50.1 % to 71.3 % in site B, respectively. Some unique OTUs and populations affiliated with Geobacter, Clostridium, Desulfosporosinus and Desulfosporosinus might be involved in Cr (VI) reduction in site A. Furthermore, GeoChip 4.0 (a comprehensive functional gene array) data showed that genes involved in carbon and nitrogen cycling and metal resistance significantly (p < 0.05) increased in site A sediment samples. All the results indicated that indigenous sediment microbial communities might be able to remediate contaminants like Cr (VI), and this information provides possible strategies for future bioremediation of the Lanzhou reach.

HuMiChip2 for strain level identification and functional profiling of human microbiomes

With the massive data generated by the Human Microbiome Project, how to transform such data into useful information and knowledge remains challenging. Here, with currently available sequencing information (reference genomes and metagenomes), we have developed a comprehensive microarray, HuMiChip2, for strain-level identification and functional characterization of human microbiomes. HuMiChip2 was composed of 29,467 strain-specific probes targeting 2063 microbial strains/species and 133,924 sequence- and group-specific probes targeting 157 key functional gene families involved in various metabolic pathways and host-microbiome interaction processes. Computational evaluation of strain-specific probes suggested that they were not only specific to mock communities of sequenced microorganisms and metagenomes from different human body sites but also to non-sequenced microbial strains. Experimental evaluation of strain-specific probes using single strains/species and mock communities suggested a high specificity of these probes with their corresponding targets. Application of HuMiChip2 to human gut microbiome samples showed the patient microbiomes of alcoholic liver cirrhosis significantly (p < 0.05) shifted their functional structure from the healthy individuals, and the relative abundance of 21 gene families significantly (p < 0.1) differed between the liver cirrhosis patients and healthy individuals. At the strain level, five Bacteroides strains were significantly (p < 0.1) and more frequently detected in liver cirrhosis patients. These results suggest that the developed HuMiChip2 is a useful microbial ecological microarray for both strain-level identification and functional profiling of human microbiomes.

Sensitivity and correlation of hypervariable regions in 16S rRNA genes in phylogenetic analysis

Background

Prokaryotic 16S ribosomal RNA (rRNA) sequences are widely used in environmental microbiology and molecular evolution as reliable markers for the taxonomic classification and phylogenetic analysis of microbes. Restricted by current sequencing techniques, the massive sequencing of 16S rRNA gene amplicons encompassing the full length of genes is not yet feasible. Thus, the selection of the most efficient hypervariable regions for phylogenetic analysis and taxonomic classification is still debated. In the present study, several bioinformatics tools were integrated to build an in silico pipeline to evaluate the phylogenetic sensitivity of the hypervariable regions compared with the corresponding full-length sequences.

Results

The correlation of seven sub-regions was inferred from the geodesic distance, a parameter that is applied to quantitatively compare the topology of different phylogenetic trees constructed using the sequences from different sub-regions. The relationship between different sub-regions based on the geodesic distance indicated that V4-V6 were the most reliable regions for representing the full-length 16S rRNA sequences in the phylogenetic analysis of most bacterial phyla, while V2 and V8 were the least reliable regions.

Conclusions

Our results suggest that V4-V6 might be optimal sub-regions for the design of universal primers with superior phylogenetic resolution for bacterial phyla. A potential relationship between function and the evolution of 16S rRNA is also discussed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0992-y) contains supplementary material, which is available to authorized users.

Long-Term Oil Contamination Alters the Molecular Ecological Networks of Soil Microbial Functional Genes

With knowledge on microbial composition and diversity, investigation of within-community interactions is a further step to elucidate microbial ecological functions, such as the biodegradation of hazardous contaminants. In this work, microbial functional molecular ecological networks were studied in both contaminated and uncontaminated soils to determine the possible influences of oil contamination on microbial interactions and potential functions. Soil samples were obtained from an oil-exploring site located in South China, and the microbial functional genes were analyzed with GeoChip, a high-throughput functional microarray. By building random networks based on null model, we demonstrated that overall network structures and properties were significantly different between contaminated and uncontaminated soils (P < 0.001). Network connectivity, module numbers, and modularity were all reduced with contamination. Moreover, the topological roles of the genes (module hub and connectors) were altered with oil contamination. Subnetworks of genes involved in alkane and polycyclic aromatic hydrocarbon degradation were also constructed. Negative co-occurrence patterns prevailed among functional genes, thereby indicating probable competition relationships. The potential "keystone" genes, defined as either "hubs" or genes with highest connectivities in the network, were further identified. The network constructed in this study predicted the potential effects of anthropogenic contamination on microbial community co-occurrence interactions.

Elevated carbon dioxide accelerates the spatial turnover of soil microbial communities

Although elevated CO2 (eCO2 ) significantly affects the α-diversity, composition, function, interaction and dynamics of soil microbial communities at the local scale, little is known about eCO2 impacts on the geographic distribution of micro-organisms regionally or globally. Here, we examined the β-diversity of 110 soil microbial communities across six free air CO2 enrichment (FACE) experimental sites using a high-throughput functional gene array. The β-diversity of soil microbial communities was significantly (P < 0.05) correlated with geographic distance under both CO2 conditions, but declined significantly (P < 0.05) faster at eCO2 with a slope of -0.0250 than at ambient CO2 (aCO2 ) with a slope of -0.0231 although it varied within each individual site, indicating that the spatial turnover rate of soil microbial communities was accelerated under eCO2 at a larger geographic scale (e.g. regionally). Both distance and soil properties significantly (P < 0.05) contributed to the observed microbial β-diversity. This study provides new hypotheses for further understanding their assembly mechanisms that may be especially important as global CO2 continues to increase.

The functional gene composition and metabolic potential of coral-associated microbial communities

The phylogenetic diversity of coral-associated microbes has been extensively examined, but some contention remains regarding whether coral-associated microbial communities are species-specific or site-specific. It is suggested that corals may associate with microbes in terms of function, although little is known about the differences in coral-associated microbial functional gene composition and metabolic potential among coral species. Here, 16S rRNA Illumina sequencing and functional gene array (GeoChip 5.0) were used to assess coral-associated microbial communities. Our results indicate that both host species and environmental variables significantly correlate with shifts in the microbial community structure and functional potential. Functional genes related to key biogeochemical cycles including carbon, nitrogen, sulfur and phosphorus cycling, metal homeostasis, organic remediation, antibiotic resistance and secondary metabolism were shown to significantly vary between and among the four study corals (Galaxea astreata, Porites lutea, Porites andrewsi and Pavona decussata). Genes specific for anammox were also detected for the first time in the coral holobiont and positively correlated with ammonium. This study reveals that variability in the functional potential of coral-associated microbial communities is largely driven by changes in environmental factors and further demonstrates the importance of linking environmental parameters with genomic data in complex environmental systems.

Planting increases the abundance and structure complexity of soil core functional genes relevant to carbon and nitrogen cycling

Plants have an important impact on soil microbial communities and their functions. However, how plants determine the microbial composition and network interactions is still poorly understood. During a four-year field experiment, we investigated the functional gene composition of three types of soils (Phaeozem, Cambisols and Acrisol) under maize planting and bare fallow regimes located in cold temperate, warm temperate and subtropical regions, respectively. The core genes were identified using high-throughput functional gene microarray (GeoChip 3.0), and functional molecular ecological networks (fMENs) were subsequently developed with the random matrix theory (RMT)-based conceptual framework. Our results demonstrated that planting significantly (P < 0.05) increased the gene alpha-diversity in terms of richness and Shannon - Simpson's indexes for all three types of soils and 83.5% of microbial alpha-diversity can be explained by the plant factor. Moreover, planting had significant impacts on the microbial community structure and the network interactions of the microbial communities. The calculated network complexity was higher under maize planting than under bare fallow regimes. The increase of the functional genes led to an increase in both soil respiration and nitrification potential with maize planting, indicating that changes in the soil microbial communities and network interactions influenced ecological functioning.

The performance of a combined nitritation-anammox reactor treating anaerobic digestion supernatant under various C/N ratios

A combined nitritation-anammox reactor was developed to treat the digestion supernatant under various C/N ratios. Due to the difficulties for heterotroph to utilize the refractory organics, the reactor presented relatively stable performance with increasing supernatant addition. Nevertheless, the adverse effects of supernatant would accumulate during the long-term operation and thus weakened the activity and shock resistance of microbes, which further led to the gradual decrease of reactor performance after 92 days' operation. Under this circumstance, supernatant with volatile fatty acids (VFAs) residuals was further introduced into the reactor to investigate the performance of combined nitritation-anammox process with VFA addition. With the appearance of VFAs, the nitrogen removal performance gradually restored and the reactor finally achieved stable and efficient performance with C/N ratio of 0.35. The VFA residuals within 150 mg/L in the supernatant served as the extra electron donors and stimulated the heterotrophic denitrification process, which was vital for the enhancement of reactor. The nitrogen removal rate and total nitrogen removal efficiency reached 0.49 kg N/(m3·day) and 88.8% after 140 days' operation, respectively. The combined nitritation-anammox reactor was proved suitable to treat digestion supernatant.

Impacts of the Three Gorges Dam on microbial structure and potential function

The Three Gorges Dam has significantly altered ecological and environmental conditions within the reservoir region, but how these changes affect bacterioplankton structure and function is unknown. Here, three widely accepted metagenomic tools were employed to study the impact of damming on the bacterioplankton community in the Xiangxi River. Our results indicated that bacterioplankton communities were both taxonomically and functionally different between backwater and riverine sites, which represent communities with and without direct dam effects, respectively. There were many more nitrogen cycling Betaproteobacteria (e.g., Limnohabitans), and a higher abundance of functional genes and KEGG orthology (KO) groups involved in nitrogen cycling in the riverine sites, suggesting a higher level of bacterial activity involved in generating more nitrogenous nutrients for the growth of phytoplankton. Additionally, the KO categories involved in carbon and sulfur metabolism, as well as most of the detected functional genes also showed clear backwater and riverine patterns. As expected, these diversity patterns all significantly correlated with environmental characteristics, confirming that the bacterioplankton communities in the Xiangxi River were really affected by environmental changes from the Three Gorges Dam. This study provides a first comparative metagenomic insight for evaluating the impacts of the large dam on microbial function.

GeoChip-based analysis of microbial community of a combined nitritation-anammox reactor treating anaerobic digestion supernatant

A combined nitritation-anammox reactor was established to treat anaerobic digestion supernatant. The reactor achieved a nitrogen loading rate of 0.5 kg N/(m(3)·d) and total nitrogen removal efficiency of 85% after 140 days' operation. To examine the microbial community responsible for the process, GeoChip 4.0, a high-throughput, microarray-based metagenomic tool, was adopted to measure microbial functional potential under different percentages of digestion supernatant. Intriguingly, our results showed that microbial community composition in a stably functioning bioreactor were significantly different under varying environmental conditions. Functional gene diversities decreased with increasing percentages of digestion supernatant. Genes involved in organic remediation and metal resistance were highly abundant, revealing new metabolic potentials in addition to nitrogen and carbon removal. Compared to the significant decrease of genes involved in denitrification and nitrification caused by inhibition of the digestion supernatant, relative abundances of genes for anammox remained relatively stable. This could be partially attributed to the protection of biofilm, which was vital for the stable performance of nitrogen removal. In addition, nitrogen compounds, C/N ratio and the operation parameters (pH and temperature) were the key variables shaping the microbial community, contributing to a total of 76.64% of the variance of the reactor.

Functional potential of soil microbial communities in the maize rhizosphere

Microbial communities in the rhizosphere make significant contributions to crop health and nutrient cycling. However, their ability to perform important biogeochemical processes remains uncharacterized. Here, we identified important functional genes that characterize the rhizosphere microbial community to understand metabolic capabilities in the maize rhizosphere using the GeoChip-based functional gene array method. Significant differences in functional gene structure were apparent between rhizosphere and bulk soil microbial communities. Approximately half of the detected gene families were significantly (p<0.05) increased in the rhizosphere. Based on the detected gyrB genes, Gammaproteobacteria, Betaproteobacteria, Firmicutes, Bacteroidetes and Cyanobacteria were most enriched in the rhizosphere compared to those in the bulk soil. The rhizosphere niche also supported greater functional diversity in catabolic pathways. The maize rhizosphere had significantly enriched genes involved in carbon fixation and degradation (especially for hemicelluloses, aromatics and lignin), nitrogen fixation, ammonification, denitrification, polyphosphate biosynthesis and degradation, sulfur reduction and oxidation. This research demonstrates that the maize rhizosphere is a hotspot of genes, mostly originating from dominant soil microbial groups such as Proteobacteria, providing functional capacity for the transformation of labile and recalcitrant organic C, N, P and S compounds.

Functional gene arrays-based analysis of fecal microbiomes in patients with liver cirrhosis

Background

Human gut microbiota plays an important role in the pathogenesis of cirrhosis complications. Although the phylogenetic diversity of intestinal microbiota in patients with liver cirrhosis has been examined in several studies, little is known about their functional composition and structure.

Results

To characterize the functional gene diversity of the gut microbiome in cirrhotic patients, we recruited a total of 42 individuals, 12 alcoholic cirrhosis patients, 18 hepatitis B virus (HBV)-related cirrhosis patients, and 12 normal controls. We determined the functional structure of these samples using a specific functional gene array, which is a combination of GeoChip for monitoring biogeochemical processes and HuMiChip specifically designed for analyzing human microbiomes. Our experimental data showed that the microbial community functional composition and structure were dramatically distinctive in the alcoholic cirrhosis. Various microbial functional genes involved in organic remediation, stress response, antibiotic resistance, metal resistance, and virulence were highly enriched in the alcoholic cirrhosis group compared to the control group and HBV-related cirrhosis group. Cirrhosis may have distinct influences on metabolic potential of fecal microbial communities. The abundance of functional genes relevant to nutrient metabolism, including amino acid metabolism, lipid metabolism, nucleotide metabolism, and isoprenoid biosynthesis, were significantly decreased in both alcoholic cirrhosis group and HBV-related cirrhosis group. Significant correlations were observed between functional gene abundances and Child-Pugh scores, such as those encoding aspartate-ammonia ligase, transaldolase, adenylosuccinate synthetase and IMP dehydrogenase.

Conclusions

Functional gene array was utilized to study the gut microbiome in alcoholic and HBV-related cirrhosis patients and controls in this study. Our array data indicated that the functional composition of fecal microbiomes was heavily influenced by cirrhosis, especially by alcoholic cirrhosis. This study provides new insights into the functional potentials and activity of gut microbiota in cirrhotic patients with different etiologies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-753) contains supplementary material, which is available to authorized users.

Microbial community structure and function of nitrobenzene reduction biocathode in response to carbon source switchover

The stress of poised cathode potential condition and carbon source switchover for functional biocathode microbial community influences is poorly understood. Using high-throughput functional gene array (GeoChip v4.2) and Illumina 16S rRNA gene MiSeq sequencing, we investigated the phylogenetic and functional microbial community of the initial inoculum and biocathode for bioelectrochemical reduction of nitrobenzene to less toxic aniline in response to carbon source switchover (from organic glucose to inorganic bicarbonate). Selective transformation of nitrobenzene to aniline maintained in the bicarbonate fed biocathode although nitrobenzene reduction rate and aniline formation rate were significantly decreased compared to those of the glucose-fed biocathode. When the electrical circuit of the glucose-fed biocathode was disconnected, both rates of nitrobenzene reduction and of aniline formation were markedly decreased, confirming the essential role of an applied electric field for the enhancement of nitrobenzene reduction. The stress of poised cathode potential condition led to clear succession of microbial communities from the initial inoculum to biocathode and the carbon source switchover obviously changed the microbial community structure of biocathode. Most of the dominant genera were capable of reducing nitroaromatics to the corresponding aromatic amines regardless of the performance mode. Heterotrophic Enterococcus was dominant in the glucose-fed biocathode while autotrophic Paracoccus and Variovorax were dominant in the bicarbonate-fed biocathode. Relatively higher intensity of diverse multi-heme cytochrome c (putatively involved in electrons transfer) and carbon fixation genes was observed in the biocarbonate-fed biocathode, likely met the requirement of the energy conservation and maintained the nitrobenzene selective reduction capability after carbon source switchover. Extracellular pilin, which are important for biofilm formation and potential conductivity, had a higher gene abundance in the glucose-fed biocathode might explain the enhancement of electro-catalysis activity for nitrobenzene reduction with glucose supply. Dominant nitroaromatics-reducing or electrochemically active bacteria and diverse functional genes related to electrons transfer and nitroaromatics reduction were associated with nitrobenzene reduction efficiency of biocathode communities in response to carbon source switchover.

Phylogenetic and functional gene structure shifts of the oral microbiomes in periodontitis patients

Determining the composition and function of subgingival dental plaque is crucial to understanding human periodontal health and disease, but it is challenging because of the complexity of the interactions between human microbiomes and human body. Here, we examined the phylogenetic and functional gene differences between periodontal and healthy individuals using MiSeq sequencing of 16S rRNA gene amplicons and a specific functional gene array (a combination of GeoChip 4.0 for biogeochemical processes and HuMiChip 1.0 for human microbiomes). Our analyses indicated that the phylogenetic and functional gene structure of the oral microbiomes were distinctly different between periodontal and healthy groups. Also, 16S rRNA gene sequencing analysis indicated that 39 genera were significantly different between healthy and periodontitis groups, and Fusobacterium, Porphyromonas, Treponema, Filifactor, Eubacterium, Tannerella, Hallella, Parvimonas, Peptostreptococcus and Catonella showed higher relative abundances in the periodontitis group. In addition, functional gene array data showed that a lower gene number but higher signal intensity of major genes existed in periodontitis, and a variety of genes involved in virulence factors, amino acid metabolism and glycosaminoglycan and pyrimidine degradation were enriched in periodontitis, suggesting their potential importance in periodontal pathogenesis. However, the genes involved in amino acid synthesis and pyrimidine synthesis exhibited a significantly lower relative abundance compared with healthy group. Overall, this study provides new insights into our understanding of phylogenetic and functional gene structure of subgingival microbial communities of periodontal patients and their importance in pathogenesis of periodontitis.

Elevated nitrate enriches microbial functional genes for potential bioremediation of complexly contaminated sediments

Nitrate is an important nutrient and electron acceptor for microorganisms, having a key role in nitrogen (N) cycling and electron transfer in anoxic sediments. High-nitrate inputs into sediments could have a significant effect on N cycling and its associated microbial processes. However, few studies have been focused on the effect of nitrate addition on the functional diversity, composition, structure and dynamics of sediment microbial communities in contaminated aquatic ecosystems with persistent organic pollutants (POPs). Here we analyzed sediment microbial communities from a field-scale in situ bioremediation site, a creek in Pearl River Delta containing a variety of contaminants including polybrominated diphenyl ethers (PBDEs) and polycyclic aromatic hydrocarbons (PAHs), before and after nitrate injection using a comprehensive functional gene array (GeoChip 4.0). Our results showed that the sediment microbial community functional composition and structure were markedly altered, and that functional genes involved in N-, carbon (C)-, sulfur (S)-and phosphorus (P)- cycling processes were highly enriched after nitrate injection, especially those microorganisms with diverse metabolic capabilities, leading to potential in situ bioremediation of the contaminated sediment, such as PBDE and PAH reduction/degradation. This study provides new insights into our understanding of sediment microbial community responses to nitrate addition, suggesting that indigenous microorganisms could be successfully stimulated for in situ bioremediation of POPs in contaminated sediments with nitrate addition.

GeoChip 4: a functional gene-array-based high-throughput environmental technology for microbial community analysis

Micro-organisms play critical roles in many important biogeochemical processes in the Earth's biosphere. However, understanding and characterizing the functional capacity of microbial communities are still difficult due to the extremely diverse and often uncultivable nature of most micro-organisms. In this study, we developed a new functional gene array, GeoChip 4, for analysing the functional diversity, composition, structure, metabolic potential/activity and dynamics of microbial communities. GeoChip 4 contained approximately 82 000 probes covering 141 995 coding sequences from 410 functional gene families related to microbial carbon (C), nitrogen (N), sulphur (S), and phosphorus (P) cycling, energy metabolism, antibiotic resistance, metal resistance/reduction, organic remediation, stress responses, bacteriophage and virulence. A total of 173 archaeal, 4138 bacterial, 404 eukaryotic and 252 viral strains were targeted, providing the ability to analyse targeted functional gene families of micro-organisms included in all four domains. Experimental assessment using different amounts of DNA suggested that as little as 500 ng environmental DNA was required for good hybridization, and the signal intensities detected were well correlated with the DNA amount used. GeoChip 4 was then applied to study the effect of long-term warming on soil microbial communities at a Central Oklahoma site, with results indicating that microbial communities respond to long-term warming by enriching carbon degradation, nutrient cycling (nitrogen and phosphorous) and stress response gene families. To the best of our knowledge, GeoChip 4 is the most comprehensive functional gene array for microbial community analysis.

Development of HuMiChip for functional profiling of human microbiomes

Understanding the diversity, composition, structure, function, and dynamics of human microbiomes in individual human hosts is crucial to reveal human-microbial interactions, especially for patients with microbially mediated disorders, but challenging due to the high diversity of the human microbiome. Here we have developed a functional gene-based microarray for profiling human microbiomes (HuMiChip) with 36,802 probes targeting 50,007 protein coding sequences for 139 key functional gene families. Computational evaluation suggested all probes included are highly specific to their target sequences. HuMiChip was used to analyze human oral and gut microbiomes, showing significantly different functional gene profiles between oral and gut microbiome. Obvious shifts of microbial functional structure and composition were observed for both patients with dental caries and periodontitis from moderate to advanced stages, suggesting a progressive change of microbial communities in response to the diseases. Consistent gene family profiles were observed by both HuMiChip and next generation sequencing technologies. Additionally, HuMiChip was able to detect gene families at as low as 0.001% relative abundance. The results indicate that the developed HuMiChip is a useful and effective tool for functional profiling of human microbiomes.

Distinct responses of soil microbial communities to elevated CO2 and O3 in a soybean agro-ecosystem

The concentrations of atmospheric carbon dioxide (CO2) and tropospheric ozone (O3) have been rising due to human activities. However, little is known about how such increases influence soil microbial communities. We hypothesized that elevated CO2 (eCO2) and elevated O3 (eO3) would significantly affect the functional composition, structure and metabolic potential of soil microbial communities, and that various functional groups would respond to such atmospheric changes differentially. To test these hypotheses, we analyzed 96 soil samples from a soybean free-air CO2 enrichment (SoyFACE) experimental site using a comprehensive functional gene microarray (GeoChip 3.0). The results showed the overall functional composition and structure of soil microbial communities shifted under eCO2, eO3 or eCO2+eO3. Key functional genes involved in carbon fixation and degradation, nitrogen fixation, denitrification and methane metabolism were stimulated under eCO2, whereas those involved in N fixation, denitrification and N mineralization were suppressed under eO3, resulting in the fact that the abundance of some eO3-supressed genes was promoted to ambient, or eCO2-induced levels by the interaction of eCO2+eO3. Such effects appeared distinct for each treatment and significantly correlated with soil properties and soybean yield. Overall, our analysis suggests possible mechanisms of microbial responses to global atmospheric change factors through the stimulation of C and N cycling by eCO2, the inhibition of N functional processes by eO3 and the interaction by eCO2 and eO3. This study provides new insights into our understanding of microbial functional processes in response to global atmospheric change in soybean agro-ecosystems.

Saliva microbiota carry caries-specific functional gene signatures

Human saliva microbiota is phylogenetically divergent among host individuals yet their roles in health and disease are poorly appreciated. We employed a microbial functional gene microarray, HuMiChip 1.0, to reconstruct the global functional profiles of human saliva microbiota from ten healthy and ten caries-active adults. Saliva microbiota in the pilot population featured a vast diversity of functional genes. No significant distinction in gene number or diversity indices was observed between healthy and caries-active microbiota. However, co-presence network analysis of functional genes revealed that caries-active microbiota was more divergent in non-core genes than healthy microbiota, despite both groups exhibited a similar degree of conservation at their respective core genes. Furthermore, functional gene structure of saliva microbiota could potentially distinguish caries-active patients from healthy hosts. Microbial functions such as Diaminopimelate epimerase, Prephenate dehydrogenase, Pyruvate-formate lyase and N-acetylmuramoyl-L-alanine amidase were significantly linked to caries. Therefore, saliva microbiota carried disease-associated functional signatures, which could be potentially exploited for caries diagnosis.

Microbial diversity in the era of omic technologies

Human life and activity depends on microorganisms, as they are responsible for providing basic elements of life. Although microbes have such a key role in sustaining basic functions for all living organisms, very little is known about their biology since only a small fraction (average 1%) can be cultured under laboratory conditions. This is even more evident when considering that >88% of all bacterial isolates belong to four bacterial phyla, the Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes. Advanced technologies, developed in the last years, promise to revolutionise the way that we characterize, identify, and study microbial communities. In this review, we present the most advanced tools that microbial ecologists can use for the study of microbial communities. Innovative microbial ecological DNA microarrays such as PhyloChip and GeoChip that have been developed for investigating the composition and function of microbial communities are presented, along with an overview of the next generation sequencing technologies. Finally, the Single Cell Genomics approach, which can be used for obtaining genomes from uncultured phyla, is outlined. This tool enables the amplification and sequencing of DNA from single cells obtained directly from environmental samples and is promising to revolutionise microbiology.

StressChip as a high-throughput tool for assessing microbial community responses to environmental stresses

Microbial community responses to environmental stresses are critical for microbial growth, survival, and adaptation. To fill major gaps in our ability to discern the influence of environmental changes on microbial communities from engineered and natural environments, a functional gene-based microarray, termed StressChip, has been developed. First, 46 functional genes involved in microbial responses to environmental stresses such as changes to temperature, osmolarity, oxidative status, nutrient limitation, or general stress response were selected and curated. A total of 22,855 probes were designed, covering 79,628 coding sequences from 985 bacterial, 76 archaeal, and 59 eukaryotic species/strains. Probe specificity was computationally verified. Second, the usefulness of functional genes as indicators of stress response was examined by surveying their distribution in metagenome data sets. The abundance of individual stress response genes is consistent with expected distributions based on respective habitats. Third, the StressChip was used to analyze marine microbial communities from the Deepwater Horizon oil spill. That functional stress response genes were detected in higher abundance (p < 0.05) in oil plume compared to nonplume samples indicated shifts in community composition and structure, consistent with previous results. In summary, StressChip provides a new tool for accessing microbial community functional structure and responses to environmental changes.

Microbial electricity generation enhances decabromodiphenyl ether (BDE-209) degradation

Due to environmental persistence and biotoxicity of polybrominated diphenyl ethers (PBDEs), it is urgent to develop potential technologies to remediate PBDEs. Introducing electrodes for microbial electricity generation to stimulate the anaerobic degradation of organic pollutants is highly promising for bioremediation. However, it is still not clear whether the degradation of PBDEs could be promoted by this strategy. In this study, we hypothesized that the degradation of PBDEs (e.g., BDE-209) would be enhanced under microbial electricity generation condition. The functional compositions and structures of microbial communities in closed-circuit microbial fuel cell (c-MFC) and open-circuit microbial fuel cell (o-MFC) systems for BDE-209 degradation were detected by a comprehensive functional gene array, GeoChip 4.0, and linked with PBDE degradations. The results indicated that distinctly different microbial community structures were formed between c-MFCs and o-MFCs, and that lower concentrations of BDE-209 and the resulting lower brominated PBDE products were detected in c-MFCs after 70-day performance. The diversity and abundance of a variety of functional genes in c-MFCs were significantly higher than those in o-MFCs. Most genes involved in chlorinated solvent reductive dechlorination, hydroxylation, methoxylation and aromatic hydrocarbon degradation were highly enriched in c-MFCs and significantly positively correlated with the removal of PBDEs. Various other microbial functional genes for carbon, nitrogen, phosphorus and sulfur cycling, as well as energy transformation process, were also significantly increased in c-MFCs. Together, these results suggest that PBDE degradation could be enhanced by introducing the electrodes for microbial electricity generation and by specifically stimulating microbial functional genes.

Strain/species-specific probe design for microbial identification microarrays

Specific identification of microorganisms in the environment is important but challenging, especially at the species/strain level. Here, we have developed a novel k-mer-based approach to select strain/species-specific probes for microbial identification with diagnostic microarrays. Application of this approach to human microbiome genomes showed that multiple (≥10 probes per strain) strain-specific 50-mer oligonucleotide probes could be designed for 2,012 of 3,421 bacterial strains of the human microbiome, and species-specific probes could be designed for most of the other strains. The method can also be used to select strain/species-specific probes for sequenced genomes in any environments, such as soil and water.

Elevated CO₂ influences microbial carbon and nitrogen cycling

Background

Elevated atmospheric CO₂ (eCO₂) has been shown to have significant effects on terrestrial ecosystems. However, little is known about its influence on the structure, composition, and functional potential of soil microbial communities, especially carbon (C) and nitrogen (N) cycling. A high-throughput functional gene array (GeoChip 3.0) was used to examine the composition, structure, and metabolic potential of soil microbial communities from a grassland field experiment after ten-year field exposure to ambient and elevated CO₂ concentrations.

Results

Distinct microbial communities were established under eCO₂. The abundance of three key C fixation genes encoding ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco), carbon monoxide dehydrogenase (CODH) and propionyl-CoA/acetyl-CoA carboxylase (PCC/ACC), significantly increased under eCO₂, and so did some C degrading genes involved in starch, cellulose, and hemicellulose. Also, nifH and nirS involved in N cycling were significantly stimulated. In addition, based on variation partitioning analysis (VPA), the soil microbial community structure was largely shaped by direct and indirect eCO₂-driven factors.

Conclusions

These findings suggest that the soil microbial community structure and their ecosystem functioning for C and N cycling were altered dramatically at eCO₂. This study provides new insights into our understanding of the feedback response of soil microbial communities to elevated CO₂ and global change.

Functional ecology of an Antarctic Dry Valley

The McMurdo Dry Valleys are the largest ice-free region in Antarctica and are critically at risk from climate change. The terrestrial landscape is dominated by oligotrophic mineral soils and extensive exposed rocky surfaces where biota are largely restricted to microbial communities, although their ability to perform the majority of geobiological processes has remained largely uncharacterized. Here, we identified functional traits that drive microbial survival and community assembly, using a metagenomic approach with GeoChip-based functional gene arrays to establish metabolic capabilities in communities inhabiting soil and rock surface niches in McKelvey Valley. Major pathways in primary metabolism were identified, indicating significant plasticity in autotrophic, heterotrophic, and diazotrophic strategies supporting microbial communities. This represents a major advance beyond biodiversity surveys in that we have now identified how putative functional ecology drives microbial community assembly. Significant differences were apparent between open soil, hypolithic, chasmoendolithic, and cryptoendolithic communities. A suite of previously unappreciated Antarctic microbial stress response pathways, thermal, osmotic, and nutrient limitation responses were identified and related to environmental stressors, offering tangible clues to the mechanisms behind the enduring success of microorganisms in this seemingly inhospitable terrain. Rocky substrates exposed to larger fluctuations in environmental stress supported greater functional diversity in stress-response pathways than soils. Soils comprised a unique reservoir of genes involved in transformation of organic hydrocarbons and lignin-like degradative pathways. This has major implications for the evolutionary origin of the organisms, turnover of recalcitrant substrates in Antarctic soils, and predicting future responses to anthropogenic pollution.

Shifts of functional gene representation in wheat rhizosphere microbial communities under elevated ozone

Although the influence of ozone (O(3)) on plants has been well studied in agroecosystems, little is known about the effect of elevated O(3) (eO(3)) on soil microbial functional communities. Here, we used a comprehensive functional gene array (GeoChip 3.0) to investigate the functional composition, and structure of rhizosphere microbial communities of Yannong 19 (O(3)-sensitive) and Yangmai 16 (O(3)-relatively sensitive) wheat (Triticum aestivum L.) cultivars under eO(3). Compared with ambient O(3) (aO(3)), eO(3) led to an increase in soil pH and total carbon (C) percentages in grain and straw of wheat plants, and reduced grain weight and soil dissolved organic carbon (DOC). Based on GeoChip hybridization signal intensities, although the overall functional structure of rhizosphere microbial communities did not significantly change by eO(3) or cultivars, the results showed that the abundance of specific functional genes involved in C fixation and degradation, nitrogen (N) fixation, and sulfite reduction did significantly (P<0.05) alter in response to eO(3) and/or wheat cultivars. Also, Yannong 19 appeared to harbor microbial functional communities in the rhizosphere more sensitive in response to eO(3) than Yangmai 16. Additionally, canonical correspondence analysis suggested that the functional structure of microbial community involved in C cycling was largely shaped by soil and plant properties including pH, DOC, microbial biomass C, C/N ratio and grain weight. This study provides new insight into our understanding of the influence of eO(3) and wheat cultivars on soil microbial communities.

Biased diversity metrics revealed by bacterial 16S pyrotags derived from different primer sets

In recent years, PCR-based pyrosequencing of 16S rRNA genes has continuously increased our understanding of complex microbial communities in various environments of the Earth. However, there is always concern on the potential biases of diversity determination using different 16S rRNA gene primer sets and covered regions. Here, we first report how bacterial 16S rRNA gene pyrotags derived from a series of different primer sets resulted in the biased diversity metrics. In total, 14 types of pyrotags were obtained from two-end pyrosequencing of 7 amplicon pools generated by 7 primer sets paired by 1 of 4 forward primers (V1F, V3F, V5F, and V7F) and 1 of 4 reverse primers (V2R, V4R, V6R, and V9R), respectively. The results revealed that: i) the activated sludge exhibited a large bacterial diversity that represented a broad range of bacterial populations and served as a good sample in this methodology research; ii) diversity metrics highly depended on the selected primer sets and covered regions; iii) paired pyrotags obtained from two-end pyrosequencing of each short amplicon displayed different diversity metrics; iv) relative abundance analysis indicated the sequencing depth affected the determination of rare bacteria but not abundant bacteria; v) the primer set of V1F and V2R significantly underestimated the diversity of activated sludge; and vi) the primer set of V3F and V4R was highly recommended for future studies due to its advantages over other primer sets. All of these findings highlight the significance of this methodology research and offer a valuable reference for peer researchers working on microbial diversity determination.

GeoChip-based analysis of the functional gene diversity and metabolic potential of soil microbial communities of mangroves

Mangroves are unique and highly productive ecosystems and harbor very special microbial communities. Although the phylogenetic diversity of sediment microbial communities of mangrove habitats has been examined extensively, little is known regarding their functional gene diversity and metabolic potential. In this study, a high-throughput functional gene array (GeoChip 4.0) was used to analyze the functional diversity, composition, structure, and metabolic potential of microbial communities in mangrove habitats from mangrove national nature reserves in China. GeoChip data indicated that these microbial communities were functionally diverse as measured by the number of genes detected, unique genes, and various diversity indices. Almost all key functional gene categories targeted by GeoChip 4.0 were detected in the mangrove microbial communities, including carbon (C) fixation, C degradation, methane generation, nitrogen (N) fixation, nitrification, denitrification, ammonification, N reduction, sulfur (S) metabolism, metal resistance, antibiotic resistance, and organic contaminant degradation. Detrended correspondence analysis (DCA) of all detected genes showed that Spartina alterniflora (HH), an invasive species, did not harbor significantly different microbial communities from Aegiceras corniculatum (THY), a native species, but did differ from other species, Kenaelia candel (QQ), Aricennia marina (BGR), and mangrove-free mud flat (GT). Canonical correspondence analysis (CCA) results indicated the microbial community structure was largely shaped by surrounding environmental variables, such as total nitrogen (TN), total carbon (TC), pH, C/N ratio, and especially salinity. This study presents a comprehensive survey of functional gene diversity of soil microbial communities from different mangrove habitats/species and provides new insights into our understanding of the functional potential of microbial communities in mangrove ecosystems.