Expression profiling of hypothetical genes in Desulfovibrio vulgaris leads to improved functional annotation

Gene Sets and Mechanisms of Sulfate-Reducing Bacteria Biofilm Formation and Quorum Sensing With Impact on Corrosion

Sulfate-reducing bacteria (SRB) have a unique ability to respire under anaerobic conditions using sulfate as a terminal electron acceptor, reducing it to hydrogen sulfide. SRB thrives in many natural environments (freshwater sediments and salty marshes), deep subsurface environments (oil wells and hydrothermal vents), and processing facilities in an industrial setting. Owing to their ability to alter the physicochemical properties of underlying metals, SRB can induce fouling, corrosion, and pipeline clogging challenges. Indigenous SRB causes oil souring and associated product loss and, subsequently, the abandonment of impacted oil wells. The sessile cells in biofilms are 1,000 times more resistant to biocides and induce 100-fold greater corrosion than their planktonic counterparts. To effectively combat the challenges posed by SRB, it is essential to understand their molecular mechanisms of biofilm formation and corrosion. Here, we examine the critical genes involved in biofilm formation and microbiologically influenced corrosion and categorize them into various functional categories. The current effort also discusses chemical and biological methods for controlling the SRB biofilms. Finally, we highlight the importance of surface engineering approaches for controlling biofilm formation on underlying metal surfaces.

Genetic Basis of Chromate Adaptation and the Role of the Pre-existing Genetic Divergence during an Experimental Evolution Study with Desulfovibrio vulgaris Populations

Hexavalent chromium [Cr(VI)] is a common environmental pollutant. However, little is known about the genetic basis of microbial evolution under Cr(VI) stress and the influence of the prior evolution histories on the subsequent evolution under Cr(VI) stress. In this study, Desulfovibrio vulgaris Hildenborough (DvH), a model sulfate-reducing bacterium, was experimentally evolved for 600 generations. By evolving the replicate populations of three genetically diverse DvH clones, including ancestor (AN, without prior experimental evolution history), non-stress-evolved EC3-10, and salt stress-evolved ES9-11, the contributions of adaptation, chance, and pre-existing genetic divergence to the evolution under Cr(VI) stress were able to be dissected. Significantly decreased lag phases under Cr(VI) stress were observed in most evolved populations, while increased Cr(VI) reduction rates were primarily observed in populations evolved from EC3-10 and ES9-11. The pre-existing genetic divergence in the starting clones showed strong influences on the changes in lag phases, growth rates, and Cr(VI) reduction rates. Additionally, the genomic mutation spectra in populations evolved from different starting clones were significantly different. A total of 14 newly mutated genes obtained mutations in at least two evolved populations, suggesting their importance in Cr(VI) adaptation. An in-frame deletion mutation of one of these genes, the chromate transporter gene DVU0426, demonstrated that it played an important role in Cr(VI) tolerance. Overall, our study identified potential key functional genes for Cr(VI) tolerance and demonstrated the important role of pre-existing genetic divergence in evolution under Cr(VI) stress conditions.

IMPORTANCE Chromium is one of the most common heavy metal pollutants of soil and groundwater. The potential of Desulfovibrio vulgaris Hildenborough in heavy metal bioremediation such as Cr(VI) reduction was reported previously; however, experimental evidence of key functional genes involved in Cr(VI) resistance are largely unknown. Given the genetic divergence of microbial populations in nature, knowledge on how this divergence affects the microbial adaptation to a new environment such as Cr(VI) stress is very limited. Taking advantage of our previous study, three groups of genetically diverse D. vulgaris Hildenborough populations with or without prior experimental evolution histories were propagated under Cr(VI) stress for 600 generations. Whole-population genome resequencing of the evolved populations revealed the genomic changes underlying the improved Cr(VI) tolerance. The strong influence of the pre-existing genetic divergence in the starting clones on evolution under Cr(VI) stress conditions was demonstrated at both phenotypic and genetic levels.

Translatomics combined with transcriptomics and proteomics reveals novel functional, recently evolved orphan genes in Escherichia coli O157:H7 (EHEC)

Background

Genomes of E. coli, including that of the human pathogen Escherichia coli O157:H7 (EHEC) EDL933, still harbor undetected protein-coding genes which, apparently, have escaped annotation due to their small size and non-essential function. To find such genes, global gene expression of EHEC EDL933 was examined, using strand-specific RNAseq (transcriptome), ribosomal footprinting (translatome) and mass spectrometry (proteome).

Results

Using the above methods, 72 short, non-annotated protein-coding genes were detected. All of these showed signals in the ribosomal footprinting assay indicating mRNA translation. Seven were verified by mass spectrometry. Fifty-seven genes are annotated in other enterobacteriaceae, mainly as hypothetical genes; the remaining 15 genes constitute novel discoveries. In addition, protein structure and function were predicted computationally and compared between EHEC-encoded proteins and 100-times randomly shuffled proteins. Based on this comparison, 61 of the 72 novel proteins exhibit predicted structural and functional features similar to those of annotated proteins. Many of the novel genes show differential transcription when grown under eleven diverse growth conditions suggesting environmental regulation. Three genes were found to confer a phenotype in previous studies, e.g., decreased cattle colonization.

Conclusions

These findings demonstrate that ribosomal footprinting can be used to detect novel protein coding genes, contributing to the growing body of evidence that hypothetical genes are not annotation artifacts and opening an additional way to study their functionality. All 72 genes are taxonomically restricted and, therefore, appear to have evolved relatively recently de novo.

Electronic supplementary material

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

A Post-Genomic View of the Ecophysiology, Catabolism and Biotechnological Relevance of Sulphate-Reducing Prokaryotes

Dissimilatory sulphate reduction is the unifying and defining trait of sulphate-reducing prokaryotes (SRP). In their predominant habitats, sulphate-rich marine sediments, SRP have long been recognized to be major players in the carbon and sulphur cycles. Other, more recently appreciated, ecophysiological roles include activity in the deep biosphere, symbiotic relations, syntrophic associations, human microbiome/health and long-distance electron transfer. SRP include a high diversity of organisms, with large nutritional versatility and broad metabolic capacities, including anaerobic degradation of aromatic compounds and hydrocarbons. Elucidation of novel catabolic capacities as well as progress in the understanding of metabolic and regulatory networks, energy metabolism, evolutionary processes and adaptation to changing environmental conditions has greatly benefited from genomics, functional OMICS approaches and advances in genetic accessibility and biochemical studies. Important biotechnological roles of SRP range from (i) wastewater and off gas treatment, (ii) bioremediation of metals and hydrocarbons and (iii) bioelectrochemistry, to undesired impacts such as (iv) souring in oil reservoirs and other environments, and (v) corrosion of iron and concrete. Here we review recent advances in our understanding of SRPs focusing mainly on works published after 2000. The wealth of publications in this period, covering many diverse areas, is a testimony to the large environmental, biogeochemical and technological relevance of these organisms and how much the field has progressed in these years, although many important questions and applications remain to be explored.

Independence of nitrate and nitrite inhibition of Desulfovibrio vulgaris Hildenborough and use of nitrite as a substrate for growth

Sulfate-reducing microbes, such as Desulfovibrio vulgaris Hildenborough, cause “souring” of petroleum reservoirs through produced sulfide and precipitate heavy metals, either as sulfides or by alteration of the metal reduction state. Thus, inhibitors of these microbes, including nitrate and nitrite ions, are studied in order to limit their impact. Nitrite is a potent inhibitor of sulfate reducers, and it has been suggested that nitrate does not inhibit these microbes directly but by reduction to nitrite, which serves as the ultimate inhibitor. Here we provide evidence that nitrate inhibition of D. vulgaris can be independent of nitrite production. We also show that D. vulgaris can use nitrite as a nitrogen source or terminal electron acceptor for growth. Moreover, we report that use of nitrite as a terminal electron acceptor requires nitrite reductase (nrfA) as a D. vulgaris nrfA mutant cannot respire nitrite but remains capable of utilizing nitrite as a nitrogen source. These results illuminate previously uncharacterized metabolic abilities of D. vulgaris that may allow niche expansion in low-sulfate environments. Understanding these abilities may lead to better control of sulfate-reducing bacteria in industrial settings and more accurate prediction of their interactions in the environment.

Metallomics of two microorganisms relevant to heavy metal bioremediation reveal fundamental differences in metal assimilation and utilization

Although as many as half of all proteins are thought to require a metal cofactor, the metalloproteomes of microorganisms remain relatively unexplored. Microorganisms from different environments are likely to vary greatly in the metals that they assimilate, not just among the metals with well-characterized roles but also those lacking any known function. Herein we investigated the metal utilization of two microorganisms that were isolated from very similar environments and are of interest because of potential roles in the immobilization of heavy metals, such as uranium and chromium. The metals assimilated and their concentrations in the cytoplasm of Desulfovibrio vulgaris strain Hildenborough (DvH) and Enterobacter cloacae strain Hanford (EcH) varied dramatically, with a larger number of metals present in Enterobacter. For example, a total of 9 and 19 metals were assimilated into their cytoplasmic fractions, respectively, and DvH did not assimilate significant amounts of zinc or copper whereas EcH assimilated both. However, bioinformatic analysis of their genome sequences revealed a comparable number of predicted metalloproteins, 813 in DvH and 953 in EcH. These allowed some rationalization of the types of metal assimilated in some cases (Fe, Cu, Mo, W, V) but not in others (Zn, Nd, Ce, Pr, Dy, Hf and Th). It was also shown that U binds an unknown soluble protein in EcH but this incorporation was the result of extracellular U binding to cytoplasmic components after cell lysis.

Comparison of strand-specific transcriptomes of enterohemorrhagic Escherichia coli O157:H7 EDL933 (EHEC) under eleven different environmental conditions including radish sprouts and cattle feces

Background

Multiple infection sources for enterohemorrhagic Escherichia coli O157:H7 (EHEC) are known, including animal products, fruit and vegetables. The ecology of this pathogen outside its human host is largely unknown and one third of its annotated genes are still hypothetical. To identify genetic determinants expressed under a variety of environmental factors, we applied strand-specific RNA-sequencing, comparing the SOLiD and Illumina systems.

Results

Transcriptomes of EHEC were sequenced under 11 different biotic and abiotic conditions: LB medium at pH4, pH7, pH9, or at 15°C; LB with nitrite or trimethoprim-sulfamethoxazole; LB-agar surface, M9 minimal medium, spinach leaf juice, surface of living radish sprouts, and cattle feces. Of 5379 annotated genes in strain EDL933 (genome and plasmid), a surprising minority of only 144 had null sequencing reads under all conditions. We therefore developed a statistical method to distinguish weakly transcribed genes from background transcription. We find that 96% of all genes and 91.5% of the hypothetical genes exhibit a significant transcriptional signal under at least one condition. Comparing SOLiD and Illumina systems, we find a high correlation between both approaches for fold-changes of the induced or repressed genes. The pathogenicity island LEE showed highest transcriptional activity in LB medium, minimal medium, and after treatment with antibiotics. Unique sets of genes, including many hypothetical genes, are highly up-regulated on radish sprouts, cattle feces, or in the presence of antibiotics. Furthermore, we observed induction of the shiga-toxin carrying phages by antibiotics and confirmed active biofilm related genes on radish sprouts, in cattle feces, and on agar plates.

Conclusions

Since only a minority of genes (2.7%) were not active under any condition tested (null reads), we suggest that the assumption of significant genome over-annotations is wrong. Environmental transcriptomics uncovered hitherto unknown gene functions and unique regulatory patterns in EHEC. For instance, the environmental function of azoR had been elusive, but this gene is highly active on radish sprouts. Thus, NGS-transcriptomics is an appropriate technique to propose new roles of hypothetical genes and to guide future research.

Electronic supplementary material

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

Genetic basis for nitrate resistance in Desulfovibrio strains

Nitrate is an inhibitor of sulfate-reducing bacteria (SRB). In petroleum production sites, amendments of nitrate and nitrite are used to prevent SRB production of sulfide that causes souring of oil wells. A better understanding of nitrate stress responses in the model SRB, Desulfovibrio vulgaris Hildenborough and Desulfovibrio alaskensis G20, will strengthen predictions of environmental outcomes of nitrate application. Nitrate inhibition of SRB has historically been considered to result from the generation of small amounts of nitrite, to which SRB are quite sensitive. Here we explored the possibility that nitrate might inhibit SRB by a mechanism other than through nitrite inhibition. We found that nitrate-stressed D. vulgaris cultures grown in lactate-sulfate conditions eventually grew in the presence of high concentrations of nitrate, and their resistance continued through several subcultures. Nitrate consumption was not detected over the course of the experiment, suggesting adaptation to nitrate. With high-throughput genetic approaches employing TnLE-seq for D. vulgaris and a pooled mutant library of D. alaskensis, we determined the fitness of many transposon mutants of both organisms in nitrate stress conditions. We found that several mutants, including homologs present in both strains, had a greatly increased ability to grow in the presence of nitrate but not nitrite. The mutated genes conferring nitrate resistance included the gene encoding the putative Rex transcriptional regulator (DVU0916/Dde_2702), as well as a cluster of genes (DVU0251-DVU0245/Dde_0597-Dde_0605) that is poorly annotated. Follow-up studies with individual D. vulgaris transposon and deletion mutants confirmed high-throughput results. We conclude that, in D. vulgaris and D. alaskensis, nitrate resistance in wild-type cultures is likely conferred by spontaneous mutations. Furthermore, the mechanisms that confer nitrate resistance may be different from those that confer nitrite resistance.

High-throughput isolation and characterization of untagged membrane protein complexes: outer membrane complexes of Desulfovibrio vulgaris

Cell membranes represent the “front line” of cellular defense and the interface between a cell and its environment. To determine the range of proteins and protein complexes that are present in the cell membranes of a target organism, we have utilized a “tagless” process for the system-wide isolation and identification of native membrane protein complexes. As an initial subject for study, we have chosen the Gram-negative sulfate-reducing bacterium Desulfovibrio vulgaris. With this tagless methodology, we have identified about two-thirds of the outer membrane- associated proteins anticipated. Approximately three-fourths of these appear to form homomeric complexes. Statistical and machine-learning methods used to analyze data compiled over multiple experiments revealed networks of additional protein–protein interactions providing insight into heteromeric contacts made between proteins across this region of the cell. Taken together, these results establish a D. vulgaris outer membrane protein data set that will be essential for the detection and characterization of environment-driven changes in the outer membrane proteome and in the modeling of stress response pathways. The workflow utilized here should be effective for the global characterization of membrane protein complexes in a wide range of organisms.

Systematic characterization of hypothetical proteins in Synechocystis sp. PCC 6803 reveals proteins functionally relevant to stress responses

We described here a global detection and functional inference of hypothetical proteins involved in stress response in Synechocystis sp. PCC 6803. In the study, we first applied an iTRAQ-LC-MS/MS based quantitative proteomics to the Synechocystis cells grown under five stress conditions. The analysis detected a total of 807 hypothetical proteins with high confidence. Among them, 480 were differentially regulated. We then applied a Weighted Gene Co-expression Network Analysis approach to construct transcriptional networks for Synechocystis under nutrient limitation and osmotic stress conditions using transcriptome datasets. The analysis showed that 305 and 467 coding genes of hypothetical proteins were functionally relevant to nutrient limitation and osmotic stress, respectively. A comparison of responsive hypothetical proteins to all stress conditions allowed identification of 22 hypothetical proteins commonly responsive to all stresses, suggesting they may be part of the core stress responses in Synechocystis. Finally, functional inference of these core stress responsive proteins using both sequence similarity and non-similarity approaches was conducted. The study provided new insights into the stress response networks in Synechocystis, and also demonstrated that a combination of experimental "OMICS" and bioinformatics methodologies could improve functional annotation for hypothetical proteins.

Transcriptomic and proteomic analyses of Desulfovibrio vulgaris biofilms: carbon and energy flow contribute to the distinct biofilm growth state

Background

Desulfovibrio vulgaris Hildenborough is a sulfate-reducing bacterium (SRB) that is intensively studied in the context of metal corrosion and heavy-metal bioremediation, and SRB populations are commonly observed in pipe and subsurface environments as surface-associated populations. In order to elucidate physiological changes associated with biofilm growth at both the transcript and protein level, transcriptomic and proteomic analyses were done on mature biofilm cells and compared to both batch and reactor planktonic populations. The biofilms were cultivated with lactate and sulfate in a continuously fed biofilm reactor, and compared to both batch and reactor planktonic populations.

Results

The functional genomic analysis demonstrated that biofilm cells were different compared to planktonic cells, and the majority of altered abundances for genes and proteins were annotated as hypothetical (unknown function), energy conservation, amino acid metabolism, and signal transduction. Genes and proteins that showed similar trends in detected levels were particularly involved in energy conservation such as increases in an annotated ech hydrogenase, formate dehydrogenase, pyruvate:ferredoxin oxidoreductase, and rnf oxidoreductase, and the biofilm cells had elevated formate dehydrogenase activity. Several other hydrogenases and formate dehydrogenases also showed an increased protein level, while decreased transcript and protein levels were observed for putative coo hydrogenase as well as a lactate permease and hyp hydrogenases for biofilm cells. Genes annotated for amino acid synthesis and nitrogen utilization were also predominant changers within the biofilm state. Ribosomal transcripts and proteins were notably decreased within the biofilm cells compared to exponential-phase cells but were not as low as levels observed in planktonic, stationary-phase cells. Several putative, extracellular proteins (DVU1012, 1545) were also detected in the extracellular fraction from biofilm cells.

Conclusions

Even though both the planktonic and biofilm cells were oxidizing lactate and reducing sulfate, the biofilm cells were physiologically distinct compared to planktonic growth states due to altered abundances of genes/proteins involved in carbon/energy flow and extracellular structures. In addition, average expression values for multiple rRNA transcripts and respiratory activity measurements indicated that biofilm cells were metabolically more similar to exponential-phase cells although biofilm cells are structured differently. The characterization of physiological advantages and constraints of the biofilm growth state for sulfate-reducing bacteria will provide insight into bioremediation applications as well as microbially-induced metal corrosion.

Towards a rigorous network of protein-protein interactions of the model sulfate reducer Desulfovibrio vulgaris Hildenborough

Protein-protein interactions offer an insight into cellular processes beyond what may be obtained by the quantitative functional genomics tools of proteomics and transcriptomics. The aforementioned tools have been extensively applied to study Escherichia coli and other aerobes and more recently to study the stress response behavior of Desulfovibrio vulgaris Hildenborough, a model obligate anaerobe and sulfate reducer and the subject of this study. Here we carried out affinity purification followed by mass spectrometry to reconstruct an interaction network among 12 chromosomally encoded bait and 90 prey proteins based on 134 bait-prey interactions identified to be of high confidence. Protein-protein interaction data are often plagued by the lack of adequate controls and replication analyses necessary to assess confidence in the results, including identification of potential false positives. We addressed these issues through the use of biological replication, exponentially modified protein abundance indices, results from an experimental negative control, and a statistical test to assign confidence to each putative interacting pair applicable to small interaction data studies. We discuss the biological significance of metabolic features of D. vulgaris revealed by these protein-protein interaction data and the observed protein modifications. These include the distinct role of the putative carbon monoxide-induced hydrogenase, unique electron transfer routes associated with different oxidoreductases, and the possible role of methylation in regulating sulfate reduction.

Shewanella oneidensis MR-1 sensory box protein involved in aerobic and anoxic growth

Although little is known of potential function for conserved signaling proteins, it is hypothesized that such proteins play important roles to coordinate cellular responses to environmental stimuli. In order to elucidate the function of a putative sensory box protein (PAS domains) in Shewanella oneidensis MR-1, the physiological role of SO3389 was characterized. The predicted open reading frame (ORF) encodes a putative sensory box protein that has PAS, GGDEF, and EAL domains, and an in-frame deletion mutant was constructed (ΔSO3389) with approximately 95% of the ORF deleted. Under aerated conditions, wild-type and mutant cultures had similar growth rates, but the mutant culture had a lower growth rate under static, aerobic conditions. Oxygen consumption rates were lower for mutant cultures (1.5-fold), and wild-type cultures also maintained lower dissolved oxygen concentrations under aerated growth conditions. When transferred to anoxic conditions, the mutant did not grow with fumarate, iron(III), or dimethyl sulfoxide (DMSO) as electron acceptors. Biochemical assays demonstrated the expression of different c-type cytochromes as well as decreased fumarate reductase activity in the mutant transferred to anoxic growth conditions. Transcriptomic studies showed the inability of the mutant to up-express and down-express genes, including c-type cytochromes (e.g., SO4047/SO4048, SO3285/SO3286), reductases (e.g., SO0768, SO1427), and potential regulators (e.g., SO1329). The complemented strain was able to grow when transferred from aerobic to anoxic growth conditions with the tested electron acceptors. The modeled structure for the SO3389 PAS domains was highly similar to the crystal structures of FAD-binding PAS domains that are known O2/redox sensors. Based on physiological, genomic, and bioinformatic results, we suggest that the sensory box protein, SO3389, is an O2/redox sensor that is involved in optimization of aerobic growth and transitions to anoxia in S. oneidensis MR-1.

How sulphate-reducing microorganisms cope with stress: lessons from systems biology

Sulphate-reducing microorganisms (SRMs) are a phylogenetically diverse group of anaerobes encompassing distinct physiologies with a broad ecological distribution. As SRMs have important roles in the biogeochemical cycling of carbon, nitrogen, sulphur and various metals, an understanding of how these organisms respond to environmental stresses is of fundamental and practical importance. In this Review, we highlight recent applications of systems biology tools in studying the stress responses of SRMs, particularly Desulfovibrio spp., at the cell, population, community and ecosystem levels. The syntrophic lifestyle of SRMs is also discussed, with a focus on system-level analyses of adaptive mechanisms. Such information is important for understanding the microbiology of the global sulphur cycle and for developing biotechnological applications of SRMs for environmental remediation, energy production, biocorrosion control, wastewater treatment and mineral recovery.

The anaerobe-specific orange protein complex of Desulfovibrio vulgaris hildenborough is encoded by two divergent operons coregulated by σ54 and a cognate transcriptional regulator

Analysis of sequenced bacterial genomes revealed that the genomes encode more than 30% hypothetical and conserved hypothetical proteins of unknown function. Among proteins of unknown function that are conserved in anaerobes, some might be determinants of the anaerobic way of life. This study focuses on two divergent clusters specifically found in anaerobic microorganisms and mainly composed of genes encoding conserved hypothetical proteins. We show that the two gene clusters DVU2103-DVU2104-DVU2105 (orp2) and DVU2107-DVU2108-DVU2109 (orp1) form two divergent operons transcribed by the σ(54)-RNA polymerase. We further demonstrate that the σ(54)-dependent transcriptional regulator DVU2106, located between orp1 and orp2, collaborates with σ(54)-RNA polymerase to orchestrate the simultaneous expression of the divergent orp operons. DVU2106, whose structural gene is transcribed by the σ(70)-RNA polymerase, negatively retrocontrols its own expression. By using an endogenous pulldown strategy, we identify a physiological complex composed of DVU2103, DVU2104, DVU2105, DVU2108, and DVU2109. Interestingly, inactivation of DVU2106, which is required for orp operon transcription, induces morphological defects that are likely linked to the absence of the ORP complex. A putative role of the ORP proteins in positioning the septum during cell division is discussed.

Integrating multiple 'omics' analysis for microbial biology: application and methodologies

Recent advances in various 'omics' technologies enable quantitative monitoring of the abundance of various biological molecules in a high-throughput manner, and thus allow determination of their variation between different biological states on a genomic scale. Several popular 'omics' platforms that have been used in microbial systems biology include transcriptomics, which measures mRNA transcript levels; proteomics, which quantifies protein abundance; metabolomics, which determines abundance of small cellular metabolites; interactomics, which resolves the whole set of molecular interactions in cells; and fluxomics, which establishes dynamic changes of molecules within a cell over time. However, no single 'omics' analysis can fully unravel the complexities of fundamental microbial biology. Therefore, integration of multiple layers of information, the multi-'omics' approach, is required to acquire a precise picture of living micro-organisms. In spite of this being a challenging task, some attempts have been made recently to integrate heterogeneous 'omics' datasets in various microbial systems and the results have demonstrated that the multi-'omics' approach is a powerful tool for understanding the functional principles and dynamics of total cellular systems. This article reviews some basic concepts of various experimental 'omics' approaches, recent application of the integrated 'omics' for exploring metabolic and regulatory mechanisms in microbes, and advances in computational and statistical methodologies associated with integrated 'omics' analyses. Online databases and bioinformatic infrastructure available for integrated 'omics' analyses are also briefly discussed.