Genome-enabled studies of anaerobic, nitrate-dependent iron oxidation in the chemolithoautotrophic bacterium Thiobacillus denitrificans.

Thiobacillus denitrificans is a chemolithoautotrophic bacterium capable of anaerobic, nitrate-dependent U(IV) and Fe(II) oxidation, both of which can strongly influence the long-term efficacy of in situ reductive immobilization of uranium in contaminated aquifers. We previously identified two c-type cytochromes involved in nitrate-dependent U(IV) oxidation in T. denitrificans and hypothesized that c-type cytochromes would also catalyze Fe(II) oxidation, as they have been found to play this role in anaerobic phototrophic Fe(II)-oxidizing bacteria. Here we report on efforts to identify genes associated with nitrate-dependent Fe(II) oxidation, namely (a) whole-genome transcriptional studies [using FeCO3, Fe(2) (+), and U(IV) oxides as electron donors under denitrifying conditions], (b) Fe(II) oxidation assays performed with knockout mutants targeting primarily highly expressed or upregulated c-type cytochromes, and (c) random transposon-mutagenesis studies with screening for Fe(II) oxidation. Assays of mutants for 26 target genes, most of which were c-type cytochromes, indicated that none of the mutants tested were significantly defective in nitrate-dependent Fe(II) oxidation. The non-defective mutants included the c 1-cytochrome subunit of the cytochrome bc 1 complex (complex III), which has relevance to a previously proposed role for this complex in nitrate-dependent Fe(II) oxidation and to current concepts of reverse electron transfer. A transposon mutant with a disrupted gene associated with NADH:ubiquinone oxidoreductase (complex I) was ~35% defective relative to the wild-type strain; this strain was similarly defective in nitrate reduction with thiosulfate as the electron donor. Overall, our results indicate that nitrate-dependent Fe(II) oxidation in T. denitrificans is not catalyzed by the same c-type cytochromes involved in U(IV) oxidation, nor have other c-type cytochromes yet been implicated in the process.

A Method to Detect Differential Gene expression in Cross-Species Hybridization Experiments at Gene and Probe Level.

MOTIVATION: Whole genome microarrays are increasingly becoming the method of choice to study responses in model organisms to disease, stressors or other stimuli. However, whole genome sequences are available for only some model organisms, and there are still many species whose genome sequences are not yet available. Cross-species studies, where arrays developed for one species are used to study gene expression in a closely related species, have been used to address this gap, with some promising results. Current analytical methods have included filtration of some probes or genes that showed low hybridization activities. But consensus filtration schemes are still not available. RESULTS: A novel masking procedure is proposed based on currently available target species sequences to filter out probes and study a cross-species data set using this masking procedure and gene-set analysis. Gene-set analysis evaluates the association of some priori defined gene groups with a phenotype of interest. Two methods, Gene Set Enrichment Analysis (GSEA) and Test of Test Statistics (ToTS) were investigated. The results showed that masking procedure together with ToTS method worked well in our data set. The results from an alternative way to study cross-species hybridization experiments without masking are also presented. We hypothesize that the multi-probes structure of Affymetrix microarrays makes it possible to aggregate the effects of both well-hybridized and poorly-hybridized probes to study a group of genes. The principles of gene-set analysis were applied to the probe-level data instead of gene-level data. The results showed that ToTS can give valuable information and thus can be used as a powerful technique for analyzing cross-species hybridization experiments. AVAILABILITY: Software in the form of R code is available at http://anson.ucdavis.edu/~ychen/cross-species.html. SUPPLEMENTARY DATA: Supplementary data are available at http://anson.ucdavis.edu/~ychen/cross-species.html.

SpaK/SpaR two-component system characterized by a structure-driven domain-fusion method and in vitro phosphorylation studies.

Here we introduce a quantitative structure-driven computational domain-fusion method, which we used to predict the structures of proteins believed to be involved in regulation of the subtilin pathway in Bacillus subtilis, and used to predict a protein-protein complex formed by interaction between the proteins. Homology modeling of SpaK and SpaR yielded preliminary structural models based on a best template for SpaK comprising a dimer of a histidine kinase, and for SpaR a response regulator protein. Our LGA code was used to identify multi-domain proteins with structure homology to both modeled structures, yielding a set of domain-fusion templates then used to model a hypothetical SpaK/SpaR complex. The models were used to identify putative functional residues and residues at the protein-protein interface, and bioinformatics was used to compare functionally and structurally relevant residues in corresponding positions among proteins with structural homology to the templates. Models of the complex were evaluated in light of known properties of the functional residues within two-component systems involving His-Asp phosphorelays. Based on this analysis, a phosphotransferase complexed with a beryllofluoride was selected as the optimal template for modeling a SpaK/SpaR complex conformation. In vitro phosphorylation studies performed using wild type and site-directed SpaK mutant proteins validated the predictions derived from application of the structure-driven domain-fusion method: SpaK was phosphorylated in the presence of (32)P-ATP and the phosphate moiety was subsequently transferred to SpaR, supporting the hypothesis that SpaK and SpaR function as sensor and response regulator, respectively, in a two-component signal transduction system, and furthermore suggesting that the structure-driven domain-fusion approach correctly predicted a physical interaction between SpaK and SpaR. Our domain-fusion algorithm leverages quantitative structure information and provides a tool for generation of hypotheses regarding protein function, which can then be tested using empirical methods.

Comparative transcriptome analysis of Methylibium petroleiphilum PM1 exposed to the fuel oxygenates methyl tert-butyl ether and ethanol.

High-density whole-genome cDNA microarrays were used to investigate substrate-dependent gene expression of Methylibium petroleiphilum PM1, one of the best-characterized aerobic methyl tert-butyl ether (MTBE)-degrading bacteria. Differential gene expression profiling was conducted with PM1 grown on MTBE and ethanol as sole carbon sources. Based on microarray high scores and protein similarity analysis, an MTBE regulon located on the megaplasmid was identified for further investigation. Putative functions for enzymes encoded in this regulon are described with relevance to the predicted MTBE degradation pathway. A new unique dioxygenase enzyme system that carries out the hydroxylation of tert-butyl alcohol to 2-methyl-2-hydroxy-1-propanol in M. petroleiphilum PM1 was discovered. Hypotheses regarding the acquisition and evolution of MTBE genes as well as the involvement of IS elements in these complex processes were formulated. The pathways for toluene, phenol, and alkane oxidation via toluene monooxygenase, phenol hydroxylase, and propane monooxygenase, respectively, were upregulated in MTBE-grown cells compared to ethanol-grown cells. Four out of nine putative cyclohexanone monooxygenases were also upregulated in MTBE-grown cells. The expression data allowed prediction of several hitherto-unknown enzymes of the upper MTBE degradation pathway in M. petroleiphilum PM1 and aided our understanding of the regulation of metabolic processes that may occur in response to pollutant mixtures and perturbations in the environment.

Whole-genome analysis of the methyl tert-butyl ether-degrading beta-proteobacterium Methylibium petroleiphilum PM1.

Methylibium petroleiphilum PM1 is a methylotroph distinguished by its ability to completely metabolize the fuel oxygenate methyl tert-butyl ether (MTBE). Strain PM1 also degrades aromatic (benzene, toluene, and xylene) and straight-chain (C(5) to C(12)) hydrocarbons present in petroleum products. Whole-genome analysis of PM1 revealed an approximately 4-Mb circular chromosome and an approximately 600-kb megaplasmid, containing 3,831 and 646 genes, respectively. Aromatic hydrocarbon and alkane degradation, metal resistance, and methylotrophy are encoded on the chromosome. The megaplasmid contains an unusual t-RNA island, numerous insertion sequences, and large repeated elements, including a 40-kb region also present on the chromosome and a 29-kb tandem repeat encoding phosphonate transport and cobalamin biosynthesis. The megaplasmid also codes for alkane degradation and was shown to play an essential role in MTBE degradation through plasmid-curing experiments. Discrepancies between the insertion sequence element distribution patterns, the distributions of best BLASTP hits among major phylogenetic groups, and the G+C contents of the chromosome (69.2%) and plasmid (66%), together with comparative genome hybridization experiments, suggest that the plasmid was recently acquired and apparently carries the genetic information responsible for PM1's ability to degrade MTBE. Comparative genomic hybridization analysis with two PM1-like MTBE-degrading environmental isolates (approximately 99% identical 16S rRNA gene sequences) showed that the plasmid was highly conserved (ca. 99% identical), whereas the chromosomes were too diverse to conduct resequencing analysis. PM1's genome sequence provides a foundation for investigating MTBE biodegradation and exploring the genetic regulation of multiple biodegradation pathways in M. petroleiphilum and other MTBE-degrading beta-proteobacteria.

Whole-genome transcriptional analysis of chemolithoautotrophic thiosulfate oxidation by Thiobacillus denitrificans under aerobic versus denitrifying conditions.

Thiobacillus denitrificans is one of the few known obligate chemolithoautotrophic bacteria capable of energetically coupling thiosulfate oxidation to denitrification as well as aerobic respiration. As very little is known about the differential expression of genes associated with key chemolithoautotrophic functions (such as sulfur compound oxidation and CO2 fixation) under aerobic versus denitrifying conditions, we conducted whole-genome, cDNA microarray studies to explore this topic systematically. The microarrays identified 277 genes (approximately 10% of the genome) as differentially expressed using RMA (robust multiarray average) statistical analysis and a twofold cutoff. Genes upregulated (ca. 6- to 150-fold) under aerobic conditions included a cluster of genes associated with iron acquisition (e.g., siderophore-related genes), a cluster of cytochrome cbb3 oxidase genes, cbbL and cbbS (encoding the large and small subunits of form I ribulose 1,5-bisphosphate carboxylase/oxygenase, or RubisCO), and multiple molecular chaperone genes. Genes upregulated (ca. 4- to 95-fold) under denitrifying conditions included nar, nir, and nor genes (associated, respectively, with nitrate reductase, nitrite reductase, and nitric oxide reductase, which catalyze successive steps of denitrification), cbbM (encoding form II RubisCO), and genes involved with sulfur compound oxidation (including two physically separated but highly similar copies of sulfide:quinone oxidoreductase and of dsrC, associated with dissimilatory sulfite reductase). Among genes associated with denitrification, relative expression levels (i.e., degree of upregulation with nitrate) tended to decrease in the order nar > nir > nor > nos. Reverse transcription-quantitative PCR analysis was used to validate these trends.

The genome sequence of the obligately chemolithoautotrophic, facultatively anaerobic bacterium Thiobacillus denitrificans.

The complete genome sequence of Thiobacillus denitrificans ATCC 25259 is the first to become available for an obligately chemolithoautotrophic, sulfur-compound-oxidizing, beta-proteobacterium. Analysis of the 2,909,809-bp genome will facilitate our molecular and biochemical understanding of the unusual metabolic repertoire of this bacterium, including its ability to couple denitrification to sulfur-compound oxidation, to catalyze anaerobic, nitrate-dependent oxidation of Fe(II) and U(IV), and to oxidize mineral electron donors. Notable genomic features include (i) genes encoding c-type cytochromes totaling 1 to 2 percent of the genome, which is a proportion greater than for almost all bacterial and archaeal species sequenced to date, (ii) genes encoding two [NiFe]hydrogenases, which is particularly significant because no information on hydrogenases has previously been reported for T. denitrificans and hydrogen oxidation appears to be critical for anaerobic U(IV) oxidation by this species, (iii) a diverse complement of more than 50 genes associated with sulfur-compound oxidation (including sox genes, dsr genes, and genes associated with the AMP-dependent oxidation of sulfite to sulfate), some of which occur in multiple (up to eight) copies, (iv) a relatively large number of genes associated with inorganic ion transport and heavy metal resistance, and (v) a paucity of genes encoding organic-compound transporters, commensurate with obligate chemolithoautotrophy. Ultimately, the genome sequence of T. denitrificans will enable elucidation of the mechanisms of aerobic and anaerobic sulfur-compound oxidation by beta-proteobacteria and will help reveal the molecular basis of this organism's role in major biogeochemical cycles (i.e., those involving sulfur, nitrogen, and carbon) and groundwater restoration.