Haichun Gao

Global Transcriptome Analysis of the Heat Shock Response of Shewanella oneidensis

Shewanella oneidensis is an important model organism for bioremediation studies because of its diverse respiratory capabilities. However, the genetic basis and regulatory mechanisms underlying the ability of S. oneidensis to survive and adapt to various environmentally relevant stresses is poorly understood. To define this organisms molecular response to elevated growth temperatures, temporal gene expression profiles were examined in cells subjected to heat stress by using whole-genome DNA microarrays for S. oneidensis. Approximately 15% (n = 711) of the total predicted S. oneidensis genes (n = 4,648) represented on the microarray were significantly up- or downregulated (P < 0.05) over a 25-min period after shift to the heat shock temperature. As expected, the majority of the genes that showed homology to known chaperones and heat shock proteins in other organisms were highly induced. In addition, a number of predicted genes, including those encoding enzymes in glycolysis and the pentose cycle, serine proteases, transcriptional regulators (MerR, LysR, and TetR families), histidine kinases, and hypothetical proteins were induced. Genes encoding membrane proteins were differentially expressed, suggesting that cells possibly alter their membrane composition or structure in response to variations in growth temperature. A substantial number of the genes encoding ribosomal proteins displayed downregulated coexpression patterns in response to heat stress, as did genes encoding prophage and flagellar proteins. Finally, a putative regulatory site with high conservation to the Escherichia coli sigma32-binding consensus sequence was identified upstream of a number of heat-inducible genes.

Constructing gene co-expression networks and predicting functions of unknown genes by random matrix theory

BackgroundLarge-scale sequencing of entire genomes has ushered in a new age in biology. One of the next grand challenges is to dissect the cellular networks consisting of many individual functional modules. Defining co-expression networks without ambiguity based on genome-wide microarray data is difficult and current methods are not robust and consistent with different data sets. This is particularly problematic for little understood organisms since not much existing biological knowledge can be exploited for determining the threshold to differentiate true correlation from random noise. Random matrix theory (RMT), which has been widely and successfully used in physics, is a powerful approach to distinguish system-specific, non-random properties embedded in complex systems from random noise. Here, we have hypothesized that the universal predictions of RMT are also applicable to biological systems and the correlation threshold can be determined by characterizing the correlation matrix of microarray profiles using random matrix theory.ResultsApplication of random matrix theory to microarray data of S. oneidensis, E. coli, yeast, A. thaliana, Drosophila, mouse and human indicates that there is a sharp transition of nearest neighbour spacing distribution (NNSD) of correlation matrix after gradually removing certain elements insider the matrix. Testing on an in silico modular model has demonstrated that this transition can be used to determine the correlation threshold for revealing modular co-expression networks. The co-expression network derived from yeast cell cycling microarray data is supported by gene annotation. The topological properties of the resulting co-expression network agree well with the general properties of biological networks. Computational evaluations have showed that RMT approach is sensitive and robust. Furthermore, evaluation on sampled expression data of an in silico modular gene system has showed that under-sampled expressions do not affect the recovery of gene co-expression network. Moreover, the cellular roles of 215 functionally unknown genes from yeast, E. coli and S. oneidensis are predicted by the gene co-expression networks using guilt-by-association principle, many of which are supported by existing information or our experimental verification, further demonstrating the reliability of this approach for gene function prediction.ConclusionOur rigorous analysis of gene expression microarray profiles using RMT has showed that the transition of NNSD of correlation matrix of microarray profile provides a profound theoretical criterion to determine the correlation threshold for identifying gene co-expression networks.

Microarray-Based Analysis of Microbial Community RNAs by Whole-Community RNA Amplification

ABSTRACT A new approach, termed whole-community RNA amplification (WCRA), was developed to provide sufficient amounts of mRNAs from environmental samples for microarray analysis. This method employs fusion primers (six to nine random nucleotides with an attached T7 promoter) for the first-strand synthesis. The shortest primer (T7N6S) gave the best results in terms of the yield and representativeness of amplification. About 1,200- to 1,800-fold amplification was obtained with amounts of the RNA templates ranging from 10 to 100 ng, and very representative detection was obtained with 50 to 100 ng total RNA. Evaluation with a Shewanella oneidensis Δfur strain revealed that the amplification method which we developed could preserve the original abundance relationships of mRNAs. In addition, to determine whether representative detection of RNAs can be achieved with mixed community samples, amplification biases were evaluated with a mixture containing equal quantities of RNAs (100 ng each) from four bacterial species, and representative amplification was also obtained. Finally, the method which we developed was applied to the active microbial populations in a denitrifying fluidized bed reactor used for denitrification of contaminated groundwater and ethanol-stimulated groundwater samples for uranium reduction. The genes expressed were consistent with the expected functions of the bioreactor and groundwater system, suggesting that this approach is useful for analyzing the functional activities of microbial communities. This is one of the first demonstrations that microarray-based technology can be used to successfully detect the activities of microbial communities from real environmental samples in a high-throughput fashion.

Global Transcriptome Analysis of the Cold Shock Response of Shewanella oneidensis MR-1 and Mutational Analysis of Its Classical Cold Shock Proteins

This study presents a global transcriptional analysis of the cold shock response of Shewanella oneidensis MR-1 after a temperature downshift from 30 degrees C to 8 or 15 degrees C based on time series microarray experiments. More than 700 genes were found to be significantly affected (P < or = 0.05) upon cold shock challenge, especially at 8 degrees C. The temporal gene expression patterns of the classical cold shock genes varied, and only some of them, most notably so1648 and so2787, were differentially regulated in response to a temperature downshift. The global response of S. oneidensis to cold shock was also characterized by the up-regulation of genes encoding membrane proteins, DNA metabolism and translation apparatus components, metabolic proteins, regulatory proteins, and hypothetical proteins. Most of the metabolic proteins affected are involved in catalytic processes that generate NADH or NADPH. Mutational analyses confirmed that the small cold shock proteins, So1648 and So2787, are involved in the cold shock response of S. oneidensis. The analyses also indicated that So1648 may function only at very low temperatures.

Metal Reduction and Iron Biomineralization by a Psychrotolerant Fe(III)-Reducing Bacterium, Shewanella sp. Strain PV-4

ABSTRACT A marine psychrotolerant, dissimilatory Fe(III)-reducing bacterium, Shewanella sp. strain PV-4, from the microbial mat at a hydrothermal vent of Loihi Seamount in the Pacific Ocean has been further characterized, with emphases on metal reduction and iron biomineralization. The strain is able to reduce metals such as Fe(III), Co(III), Cr(VI), Mn(IV), and U(VI) as electron acceptors while using lactate, formate, pyruvate, or hydrogen as an electron donor. Growth during iron reduction occurred over the pH range of 7.0 to 8.9, a sodium chloride range of 0.05 to 5%, and a temperature range of 0 to 37°C, with an optimum growth temperature of 18°C. Unlike mesophilic dissimilatory Fe(III)-reducing bacteria, which produce mostly superparamagnetic magnetite (<35 nm), this psychrotolerant bacterium produces well-formed single-domain magnetite (>35 nm) at temperatures from 18 to 37°C. The genome size of this strain is about 4.5 Mb. Strain PV-4 is sensitive to a variety of commonly used antibiotics except ampicillin and can acquire exogenous DNA (plasmid pCM157) through conjugation.

Reduction of nitrate in Shewanella oneidensis depends on atypical NAP and NRF systems with NapB as a preferred electron transport protein from CymA to NapA

In the genome of Shewanella oneidensis, a napDAGHB gene cluster encoding periplasmic nitrate reductase (NapA) and accessory proteins and an nrfA gene encoding periplasmic nitrite reductase (NrfA) have been identified. These two systems seem to be atypical because the genome lacks genes encoding cytoplasmic membrane electron transport proteins, NapC for NAP and NrfBCD/NrfH for NRF, respectively. Here, we present evidence that reduction of nitrate to ammonium in S. oneidensis is carried out by these atypical systems in a two-step manner. Transcriptional and mutational analyses suggest that CymA, a cytoplasmic membrane electron transport protein, is likely to be the functional replacement of both NapC and NrfH in S. oneidensis. Surprisingly, a strain devoid of napB encoding the small subunit of nitrate reductase exhibited the maximum cell density sooner than the wild type. Further characterization of this strain showed that nitrite was not detected as a free intermediate in its culture and NapB provides a fitness gain for S. oneidensis to compete for nitrate in the environments. On the basis results from mutational analyses of napA, napB, nrfA and napBnrfA in-frame deletion mutants, we propose that NapB is able to favor nitrate reduction by routing electrons to NapA exclusively.

Probing regulon of ArcA in Shewanella oneidensis MR-1 by integrated genomic analyses

BackgroundThe Arc two-component system is a global regulator controlling many genes involved in aerobic/anaerobic respiration and fermentative metabolism in Escherichia coli. Shewanella oneidensis MR-1 contains a gene encoding a putative ArcA homolog with ~81% amino acid sequence identity to the E. coli ArcA protein but not a full-length arcB gene.ResultsTo understand the role of ArcA in S. oneidensis, an arcA deletion strain was constructed and subjected to both physiological characterization and microarray analysis. Compared to the wild-type MR-1, the mutant exhibited impaired aerobic growth and a defect in utilizing DMSO in the absence of O2. Microarray analyses on cells grown aerobically and anaerobically on fumarate revealed that expression of 1009 genes was significantly affected (p < 0.05) by the mutation. In contrast to E. coli ArcA, the protein appears to be dispensable in regulation of the TCA cycle in S. oneidensis. To further determine genes regulated by the Arc system, an ArcA recognition weight matrix from DNA-binding data and bioinformatics analysis was generated and used to produce an ArcA sequence affinity map. By combining both techniques, we identified an ArcA regulon of at least 50 operons, of which only 6 were found to be directly controlled by ArcA in E. coli.ConclusionThese results indicate that the Arc system in S. oneidensis differs from that in E. coli substantially in terms of its physiological function and regulon while their binding motif are strikingly similar.

Impacts of Shewanella oneidensis c‐type cytochromes on aerobic and anaerobic respiration

Shewanella are renowned for their ability to utilize a wide range of electron acceptors (EA) for respiration, which has been partially accredited to the presence of a large number of the c‐type cytochromes. To investigate the involvement of c‐type cytochrome proteins in aerobic and anaerobic respiration of Shewanella oneidensis Mr ‐1, 36 in‐frame deletion mutants, among possible 41 predicted, c‐type cytochrome genes were obtained. The potential involvement of each individual c‐type cytochrome in the reduction of a variety of EAs was assessed individually as well as in competition experiments. While results on the well‐studied c‐type cytochromes CymA(SO4591) and MtrC(SO1778) were consistent with previous findings, collective observations were very interesting: the responses of S. oneidensis Mr ‐1 to low and highly toxic metals appeared to be significantly different; CcoO, CcoP and PetC, proteins involved in aerobic respiration in various organisms, played critical roles in both aerobic and anaerobic respiration with highly toxic metals as EA. In addition, these studies also suggested that an uncharacterized c‐type cytochrome (SO4047) may be important to both aerobiosis and anaerobiosis.

Pellicle formation in Shewanella oneidensis

BackgroundAlthough solid surface-associated biofilm development of S. oneidensis has been extensively studied in recent years, pellicles formed at the air-liquid interface are largely overlooked. The goal of this work was to understand basic requirements and mechanism of pellicle formation in S. oneidensis.ResultsWe demonstrated that pellicle formation can be completed when oxygen and certain cations were present. Ca(II), Mn(II), Cu(II), and Zn(II) were essential for the process evidenced by fully rescuing pellicle formation of S. oneidensis from the EDTA treatment while Mg (II), Fe(II), and Fe(III) were much less effective. Proteins rather than DNA were crucial in pellicle formation and the major exopolysaccharides may be rich in mannose. Mutational analysis revealed that flagella were not required for pellicle formation but flagellum-less mutants delayed pellicle development substantially, likely due to reduced growth in static media. The analysis also demonstrated that AggA type I secretion system was essential in formation of pellicles but not of solid surface-associated biofilms in S. oneidensis.ConclusionThis systematic characterization of pellicle formation shed lights on our understanding of biofilm formation in S. oneidensis and indicated that the pellicle may serve as a good research model for studying bacterial communities.

A high-throughput percentage-of-binding strategy to measure binding energies in DNA–protein interactions: application to genome-scale site discovery

Quantifying the binding energy in DNA–protein interactions is of critical importance to understand transcriptional regulation. Based on a simple computational model, this study describes a high-throughput percentage-of-binding strategy to measure the binding energy in DNA–protein interactions between the Shewanella oneidensis ArcA two-component transcription factor protein and a systematic set of mutants in an ArcA-P (phosphorylated ArcA) binding site. The binding energies corresponding to each of the 4 nt at each position in the 15-bp binding site were used to construct a position-specific energy matrix (PEM) that allowed a reliable prediction of ArcA-P binding sites not only in Shewanella but also in related bacterial genomes.