Small and mighty: adaptation of superphylum Patescibacteria to groundwater environment drives their genome simplicity.

BACKGROUND: The newly defined superphylum Patescibacteria such as Parcubacteria (OD1) and Microgenomates (OP11) has been found to be prevalent in groundwater, sediment, lake, and other aquifer environments. Recently increasing attention has been paid to this diverse superphylum including > 20 candidate phyla (a large part of the candidate phylum radiation, CPR) because it refreshed our view of the tree of life. However, adaptive traits contributing to its prevalence are still not well known. RESULTS: Here, we investigated the genomic features and metabolic pathways of Patescibacteria in groundwater through genome-resolved metagenomics analysis of > 600 Gbp sequence data. We observed that, while the members of Patescibacteria have reduced genomes (~ 1 Mbp) exclusively, functions essential to growth and reproduction such as genetic information processing were retained. Surprisingly, they have sharply reduced redundant and nonessential functions, including specific metabolic activities and stress response systems. The Patescibacteria have ultra-small cells and simplified membrane structures, including flagellar assembly, transporters, and two-component systems. Despite the lack of CRISPR viral defense, the bacteria may evade predation through deletion of common membrane phage receptors and other alternative strategies, which may explain the low representation of prophage proteins in their genomes and lack of CRISPR. By establishing the linkages between bacterial features and the groundwater environmental conditions, our results provide important insights into the functions and evolution of this CPR group. CONCLUSIONS: We found that Patescibacteria has streamlined many functions while acquiring advantages such as avoiding phage invasion, to adapt to the groundwater environment. The unique features of small genome size, ultra-small cell size, and lacking CRISPR of this large lineage are bringing new understandings on life of Bacteria. Our results provide important insights into the mechanisms for adaptation of the superphylum in the groundwater environments, and demonstrate a case where less is more, and small is mighty.

Large Circular Plasmids from Groundwater Plasmidomes Span Multiple Incompatibility Groups and Are Enriched in Multimetal Resistance Genes.

Naturally occurring plasmids constitute a major category of mobile genetic elements responsible for harboring and transferring genes important in survival and fitness. A targeted evaluation of plasmidomes can reveal unique adaptations required by microbial communities. We developed a model system to optimize plasmid DNA isolation procedures targeted to groundwater samples which are typically characterized by low cell density (and likely variations in the plasmid size and copy numbers). The optimized method resulted in successful identification of several hundred circular plasmids, including some large plasmids (11 plasmids more than 50 kb in size, with the largest being 1.7 Mb in size). Several interesting observations were made from the analysis of plasmid DNA isolated in this study. The plasmid pool (plasmidome) was more conserved than the corresponding microbiome distribution (16S rRNA based). The circular plasmids were diverse as represented by the presence of seven plasmid incompatibility groups. The genes carried on these groundwater plasmids were highly enriched in metal resistance. Results from this study confirmed that traits such as metal, antibiotic, and phage resistance along with toxin-antitoxin systems are encoded on abundant circular plasmids, all of which could confer novel and advantageous traits to their hosts. This study confirms the ecological role of the plasmidome in maintaining the latent capacity of a microbiome, enabling rapid adaptation to environmental stresses.IMPORTANCE Plasmidomes have been typically studied in environments abundant in bacteria, and this is the first study to explore plasmids from an environment characterized by low cell density. We specifically target groundwater, a significant source of water for human/agriculture use. We used samples from a well-studied site and identified hundreds of circular plasmids, including one of the largest sizes reported in plasmidome studies. The striking similarity of the plasmid-borne ORFs in terms of taxonomical and functional classifications across several samples suggests a conserved plasmid pool, in contrast to the observed variability in the 16S rRNA-based microbiome distribution. Additionally, the stress response to environmental factors has stronger conservation via plasmid-borne genes as marked by abundance of metal resistance genes. Last, identification of novel and diverse plasmids enriches the existing plasmid database(s) and serves as a paradigm to increase the repertoire of biological parts that are available for modifying novel environmental strains.

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.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.

Peculiar citric acid cycle of hydrothermal vent chemolithoautotroph Hydrogenovibrio crunogenus, and insights into carbon metabolism by obligate autotrophs.

The genome sequence of the obligate chemolithoautotroph Hydrogenovibrio crunogenus paradoxically predicts a complete oxidative citric acid cycle (CAC). This prediction was tested by multiple approaches including whole cell carbon assimilation to verify obligate autotrophy, phylogenetic analysis of CAC enzyme sequences and enzyme assays. Hydrogenovibrio crunogenus did not assimilate any of the organic compounds provided (acetate, succinate, glucose, yeast extract, tryptone). Enzyme activities confirmed that its CAC is mostly uncoupled from the NADH pool. 2-Oxoglutarate:ferredoxin oxidoreductase activity is absent, though pyruvate:ferredoxin oxidoreductase is present, indicating that sequence-based predictions of substrate for this oxidoreductase were incorrect, and that H. crunogenus may have an incomplete CAC. Though the H. crunogenus CAC genes encode uncommon enzymes, the taxonomic distribution of their top matches suggests that they were not horizontally acquired. Comparison of H. crunogenus CAC genes to those present in other 'Proteobacteria' reveals that H. crunogenus and other obligate autotrophs lack the functional redundancy for the steps of the CAC typical for facultative autotrophs and heterotrophs, providing another possible mechanism for obligate autotrophy.

Construction of Viable Soil Defined Media Using Quantitative Metabolomics Analysis of Soil Metabolites.

Exometabolomics enables analysis of metabolite utilization of low molecular weight organic substances by soil bacteria. Environmentally-based defined media are needed to examine ecologically relevant patterns of substrate utilization. Here, we describe an approach for the construction of defined media using untargeted characterization of water soluble soil microbial metabolites from a saprolite soil collected from the Oak Ridge Field Research Center (ORFRC). To broadly characterize metabolites, both liquid chromatography mass spectrometry (LC/MS) and gas chromatography mass spectrometry (GC/MS) were used. With this approach, 96 metabolites were identified, including amino acids, amino acid derivatives, sugars, sugar alcohols, mono- and di-carboxylic acids, nucleobases, and nucleosides. From this pool of metabolites, 25 were quantified. Molecular weight cut-off filtration determined the fraction of carbon accounted for by the quantified metabolites and revealed that these soil metabolites have an uneven quantitative distribution (e.g., trehalose accounted for 9.9% of the <1 kDa fraction). This quantitative information was used to formulate two soil defined media (SDM), one containing 23 metabolites (SDM1) and one containing 46 (SDM2). To evaluate the viability of the SDM, we examined the growth of 30 phylogenetically diverse soil bacterial isolates from the ORFRC field site. The simpler SDM1 supported the growth of 13 isolates while the more complex SDM2 supported 15 isolates. To investigate SDM1 substrate preferences, one isolate, Pseudomonas corrugata strain FW300-N2E2 was selected for a time-series exometabolomics analysis. Interestingly, it was found that this organism preferred lower-abundance substrates such as guanine, glycine, proline and arginine and glucose and did not utilize the more abundant substrates maltose, mannitol, trehalose and uridine. These results demonstrate the viability and utility of using exometabolomics to construct a tractable environmentally relevant media. We anticipate that this approach can be expanded to other environments to enhance isolation and characterization of diverse microbial communities.

Rapid detection of microbial cell abundance in aquatic systems.

The detection and quantification of naturally occurring microbial cellular densities is an essential component of environmental systems monitoring. While there are a number of commonly utilized approaches for monitoring microbial abundance, capacitance-based biosensors represent a promising approach because of their low-cost and label-free detection of microbial cells, but are not as well characterized as more traditional methods. Here, we investigate the applicability of enhanced alternating current electrokinetics (ACEK) capacitive sensing as a new application for rapidly detecting and quantifying microbial cellular densities in cultured and environmentally sourced aquatic samples. ACEK capacitive sensor performance was evaluated using two distinct and dynamic systems - the Great Australian Bight and groundwater from the Oak Ridge Reservation in Oak Ridge, TN. Results demonstrate that ACEK capacitance-based sensing can accurately determine microbial cell counts throughout cellular concentrations typically encountered in naturally occurring microbial communities (10(3)-10(6) cells/mL). A linear relationship was observed between cellular density and capacitance change correlations, allowing a simple linear curve fitting equation to be used for determining microbial abundances in unknown samples. This work provides a foundation for understanding the limits of capacitance-based sensing in natural environmental samples and supports future efforts focusing on evaluating the robustness ACEK capacitance-based within aquatic environments.

Molybdenum Availability Is Key to Nitrate Removal in Contaminated Groundwater Environments.

The concentrations of molybdenum (Mo) and 25 other metals were measured in groundwater samples from 80 wells on the Oak Ridge Reservation (ORR) (Oak Ridge, TN), many of which are contaminated with nitrate, as well as uranium and various other metals. The concentrations of nitrate and uranium were in the ranges of 0.1 µM to 230 mM and <0.2 nM to 580 µM, respectively. Almost all metals examined had significantly greater median concentrations in a subset of wells that were highly contaminated with uranium (≥126 nM). They included cadmium, manganese, and cobalt, which were 1,300- to 2,700-fold higher. A notable exception, however, was Mo, which had a lower median concentration in the uranium-contaminated wells. This is significant, because Mo is essential in the dissimilatory nitrate reduction branch of the global nitrogen cycle. It is required at the catalytic site of nitrate reductase, the enzyme that reduces nitrate to nitrite. Moreover, more than 85% of the groundwater samples contained less than 10 nM Mo, whereas concentrations of 10 to 100 nM Mo were required for efficient growth by nitrate reduction for two Pseudomonas strains isolated from ORR wells and by a model denitrifier, Pseudomonas stutzeri RCH2. Higher concentrations of Mo tended to inhibit the growth of these strains due to the accumulation of toxic concentrations of nitrite, and this effect was exacerbated at high nitrate concentrations. The relevance of these results to a Mo-based nitrate removal strategy and the potential community-driving role that Mo plays in contaminated environments are discussed.

Natural bacterial communities serve as quantitative geochemical biosensors.

UNLABELLED: Biological sensors can be engineered to measure a wide range of environmental conditions. Here we show that statistical analysis of DNA from natural microbial communities can be used to accurately identify environmental contaminants, including uranium and nitrate at a nuclear waste site. In addition to contamination, sequence data from the 16S rRNA gene alone can quantitatively predict a rich catalogue of 26 geochemical features collected from 93 wells with highly differing geochemistry characteristics. We extend this approach to identify sites contaminated with hydrocarbons from the Deepwater Horizon oil spill, finding that altered bacterial communities encode a memory of prior contamination, even after the contaminants themselves have been fully degraded. We show that the bacterial strains that are most useful for detecting oil and uranium are known to interact with these substrates, indicating that this statistical approach uncovers ecologically meaningful interactions consistent with previous experimental observations. Future efforts should focus on evaluating the geographical generalizability of these associations. Taken as a whole, these results indicate that ubiquitous, natural bacterial communities can be used as in situ environmental sensors that respond to and capture perturbations caused by human impacts. These in situ biosensors rely on environmental selection rather than directed engineering, and so this approach could be rapidly deployed and scaled as sequencing technology continues to become faster, simpler, and less expensive. IMPORTANCE: Here we show that DNA from natural bacterial communities can be used as a quantitative biosensor to accurately distinguish unpolluted sites from those contaminated with uranium, nitrate, or oil. These results indicate that bacterial communities can be used as environmental sensors that respond to and capture perturbations caused by human impacts.

Advances in monitoring environmental microbes.

Culture-independent approaches, such as next-generation sequencing and microarray-based tools, provide insight into the identity and functional diversity of microbial communities. Although these approaches are potentially powerful tools in understanding microbial structure and function, there are a number of limitations that may bias conclusions. In order to mitigate these biases, an understanding of potential biases within each stage of the experimental process is necessary. This review focuses on the biases associated with sample collection, nucleic acid extraction, processing, sequencing analyses, and Chip technologies used in microbial ecology studies.

In-silico identification of phenotype-biased functional modules.

BACKGROUND: Phenotypes exhibited by microorganisms can be useful for several purposes, e.g., ethanol as an alternate fuel. Sometimes, the target phenotype maybe required in combination with other phenotypes, in order to be useful, for e.g., an industrial process may require that the organism survive in an anaerobic, alcohol rich environment and be able to feed on both hexose and pentose sugars to produce ethanol. This combination of traits may not be available in any existing organism or if they do exist, the mechanisms involved in the phenotype-expression may not be efficient enough to be useful. Thus, it may be required to genetically modify microorganisms. However, before any genetic modification can take place, it is important to identify the underlying cellular subsystems responsible for the expression of the target phenotype. RESULTS: In this paper, we develop a method to identify statistically significant and phenotypically-biased functional modules. The method can compare the organismal network information from hundreds of phenotype expressing and phenotype non-expressing organisms to identify cellular subsystems that are more prone to occur in phenotype-expressing organisms than in phenotype non-expressing organisms. We have provided literature evidence that the phenotype-biased modules identified for phenotypes such as hydrogen production (dark and light fermentation), respiration, gram-positive, gram-negative and motility, are indeed phenotype-related. CONCLUSION: Thus we have proposed a methodology to identify phenotype-biased cellular subsystems. We have shown the effectiveness of our methodology by applying it to several target phenotypes. The code and all supplemental files can be downloaded from (http://freescience.org/cs/phenotype-biased-biclusters/).

NIBBS-search for fast and accurate prediction of phenotype-biased metabolic systems.

Understanding of genotype-phenotype associations is important not only for furthering our knowledge on internal cellular processes, but also essential for providing the foundation necessary for genetic engineering of microorganisms for industrial use (e.g., production of bioenergy or biofuels). However, genotype-phenotype associations alone do not provide enough information to alter an organism's genome to either suppress or exhibit a phenotype. It is important to look at the phenotype-related genes in the context of the genome-scale network to understand how the genes interact with other genes in the organism. Identification of metabolic subsystems involved in the expression of the phenotype is one way of placing the phenotype-related genes in the context of the entire network. A metabolic system refers to a metabolic network subgraph; nodes are compounds and edges labels are the enzymes that catalyze the reaction. The metabolic subsystem could be part of a single metabolic pathway or span parts of multiple pathways. Arguably, comparative genome-scale metabolic network analysis is a promising strategy to identify these phenotype-related metabolic subsystems. Network Instance-Based Biased Subgraph Search (NIBBS) is a graph-theoretic method for genome-scale metabolic network comparative analysis that can identify metabolic systems that are statistically biased toward phenotype-expressing organismal networks. We set up experiments with target phenotypes like hydrogen production, TCA expression, and acid-tolerance. We show via extensive literature search that some of the resulting metabolic subsystems are indeed phenotype-related and formulate hypotheses for other systems in terms of their role in phenotype expression. NIBBS is also orders of magnitude faster than MULE, one of the most efficient maximal frequent subgraph mining algorithms that could be adjusted for this problem. Also, the set of phenotype-biased metabolic systems output by NIBBS comes very close to the set of phenotype-biased subgraphs output by an exact maximally-biased subgraph enumeration algorithm ( MBS-Enum ). The code (NIBBS and the module to visualize the identified subsystems) is available at http://freescience.org/cs/NIBBS.

SPICE: discovery of phenotype-determining component interplays.

BACKGROUND: A latent behavior of a biological cell is complex. Deriving the underlying simplicity, or the fundamental rules governing this behavior has been the Holy Grail of systems biology. Data-driven prediction of the system components and their component interplays that are responsible for the target system's phenotype is a key and challenging step in this endeavor. RESULTS: The proposed approach, which we call System Phenotype-related Interplaying Components Enumerator (SPICE), iteratively enumerates statistically significant system components that are hypothesized (1) to play an important role in defining the specificity of the target system's phenotype(s); (2) to exhibit a functionally coherent behavior, namely, act in a coordinated manner to perform the phenotype-specific function; and (3) to improve the predictive skill of the system's phenotype(s) when used collectively in the ensemble of predictive models. SPICE can be applied to both instance-based data and network-based data. When validated, SPICE effectively identified system components related to three target phenotypes: biohydrogen production, motility, and cancer. Manual results curation agreed with the known phenotype-related system components reported in literature. Additionally, using the identified system components as discriminatory features improved the prediction accuracy by 10% on the phenotype-classification task when compared to a number of state-of-the-art methods applied to eight benchmark microarray data sets. CONCLUSION: We formulate a problem--enumeration of phenotype-determining system component interplays--and propose an effective methodology (SPICE) to address this problem. SPICE improved identification of cancer-related groups of genes from various microarray data sets and detected groups of genes associated with microbial biohydrogen production and motility, many of which were reported in literature. SPICE also improved the predictive skill of the system's phenotype determination compared to individual classifiers and/or other ensemble methods, such as bagging, boosting, random forest, nearest shrunken centroid, and random forest variable selection method.

Efficient α, ß-motif finder for identification of phenotype-related functional modules.

BACKGROUND: Microbial communities in their natural environments exhibit phenotypes that can directly cause particular diseases, convert biomass or wastewater to energy, or degrade various environmental contaminants. Understanding how these communities realize specific phenotypic traits (e.g., carbon fixation, hydrogen production) is critical for addressing health, bioremediation, or bioenergy problems. RESULTS: In this paper, we describe a graph-theoretical method for in silico prediction of the cellular subsystems that are related to the expression of a target phenotype. The proposed (α, ß)-motif finder approach allows for identification of these phenotype-related subsystems that, in addition to metabolic subsystems, could include their regulators, sensors, transporters, and even uncharacterized proteins. By comparing dozens of genome-scale networks of functionally associated proteins, our method efficiently identifies those statistically significant functional modules that are in at least α networks of phenotype-expressing organisms but appear in no more than ß networks of organisms that do not exhibit the target phenotype. It has been shown via various experiments that the enumerated modules are indeed related to phenotype-expression when tested with different target phenotypes like hydrogen production, motility, aerobic respiration, and acid-tolerance. CONCLUSION: Thus, we have proposed a methodology that can identify potential statistically significant phenotype-related functional modules. The functional module is modeled as an (α, ß)-clique, where α and ß are two criteria introduced in this work. We also propose a novel network model, called the two-typed, divided network. The new network model and the criteria make the problem tractable even while very large networks are being compared. The code can be downloaded from http://www.freescience.org/cs/ABClique/

DENSE: efficient and prior knowledge-driven discovery of phenotype-associated protein functional modules.

BACKGROUND: Identifying cellular subsystems that are involved in the expression of a target phenotype has been a very active research area for the past several years. In this paper, cellular subsystem refers to a group of genes (or proteins) that interact and carry out a common function in the cell. Most studies identify genes associated with a phenotype on the basis of some statistical bias, others have extended these statistical methods to analyze functional modules and biological pathways for phenotype-relatedness. However, a biologist might often have a specific question in mind while performing such analysis and most of the resulting subsystems obtained by the existing methods might be largely irrelevant to the question in hand. Arguably, it would be valuable to incorporate biologist's knowledge about the phenotype into the algorithm. This way, it is anticipated that the resulting subsytems would not only be related to the target phenotype but also contain information that the biologist is likely to be interested in. RESULTS: In this paper we introduce a fast and theoretically guranteed method called DENSE (Dense and ENriched Subgraph Enumeration) that can take in as input a biologist's prior knowledge as a set of query proteins and identify all the dense functional modules in a biological network that contain some part of the query vertices. The density (in terms of the number of network egdes) and the enrichment (the number of query proteins in the resulting functional module) can be manipulated via two parameters γ and µ, respectively. CONCLUSION: This algorithm has been applied to the protein functional association network of Clostridium acetobutylicum ATCC 824, a hydrogen producing, acid-tolerant organism. The algorithm was able to verify relationships known to exist in literature and also some previously unknown relationships including those with regulatory and signaling functions. Additionally, we were also able to hypothesize that some uncharacterized proteins are likely associated with the target phenotype. The DENSE code can be downloaded from http://www.freescience.org/cs/DENSE/