Analysis of gene expression in operons of Streptomyces coelicolor

Autoinducer-2 promotes adherence of Aeromonas veronii through facilitating the expression of MSHA type IV pili genes mediated by c-di-GMP

IMPORTANCE

Aeromonas veronii can adhere to host cells through different adherence factors including outer-membrane proteins (OMPs), lipopolysaccharide (LPS), and pili, but its adherence mechanisms are still unclear. Here, we evaluated the effect of autoinducer-2 (AI-2) on adherence of A. veronii and its regulation mechanism. After determination of the promotion effect of AI-2 on adherence, we investigated which adherence factor was regulated by AI-2, and the results show that AI-2 only limits the formation of pili. Among the four distinct pili systems, only the mannose-sensitive hemagglutinin (MSHA) type IV pili genes were significantly downregulated after deficiency of AI-2. MshE, an ATPase belonged to MSHA type IV pilin, was confirmed as c-di-GMP receptor, that can bind with c-di-GMP which is positively regulated by AI-2, and the increase of c-di-GMP can promote the expression of MSHA type IV pili genes and adherence of A. veronii. Therefore, this study confirms that c-di-GMP positively regulated by AI-2 binds with MshE, then increases the expression of MSHA pili genes, finally promoting adherence of A. veronii, suggesting a multilevel positive regulatory adhesion mechanism that is responsible for A. veronii adherence.

RNA helicase-regulated processing of the Synechocystis rimO-crhR operon results in differential cistron expression and accumulation of two sRNAs

The arrangement of functionally-related genes in operons is a fundamental element of how genetic information is organized in prokaryotes. This organization ensures coordinated gene expression by co-transcription. Often, however, alternative genetic responses to specific stress conditions demand the discoordination of operon expression. During cold temperature stress, accumulation of the gene encoding the sole Asp-Glu-Ala-Asp (DEAD)-box RNA helicase in Synechocystis sp. PCC 6803, crhR (slr0083), increases 15-fold. Here, we show that crhR is expressed from a dicistronic operon with the methylthiotransferase rimO/miaB (slr0082) gene, followed by rapid processing of the operon transcript into two monocistronic mRNAs. This cleavage event is required for and results in destabilization of the rimO transcript. Results from secondary structure modeling and analysis of RNase E cleavage of the rimO-crhR transcript in vitro suggested that CrhR plays a role in enhancing the rate of the processing in an auto-regulatory manner. Moreover, two putative small RNAs are generated from additional processing, degradation, or both of the rimO transcript. These results suggest a role for the bacterial RNA helicase CrhR in RNase E-dependent mRNA processing in Synechocystis and expand the known range of organisms possessing small RNAs derived from processing of mRNA transcripts.

In vivo functional analysis of a class A β-lactamase-related protein essential for clavulanic acid biosynthesis in Streptomyces clavuligerus

In Streptomyces clavuligerus, the gene cluster involved in the biosynthesis of the clinically used β-lactamase inhibitor clavulanic acid contains a gene (orf12 or cpe) encoding a protein with a C-terminal class A β-lactamase-like domain. The cpe gene is essential for clavulanic acid production, and the recent crystal structure of its product (Cpe) was shown to also contain an N-terminal isomerase/cyclase-like domain, but the function of the protein remains unknown. In the current study, we show that Cpe is a cytoplasmic protein and that both its N- and C-terminal domains are required for in vivo clavulanic acid production in S. clavuligerus. Our results along with those from previous studies allude towards a biosynthetic role for Cpe during the later stages of clavulanic acid production in S. clavuligerus. Amino acids from Cpe essential for biosynthesis were also identified, including one (Lys₈₉) from the recently described N-terminal isomerase-like domain of unknown function. Homologues of Cpe from other clavulanic acid-producing Streptomyces spp. were shown to be functionally equivalent to the S. clavuligerus protein, whereas those from non-producers containing clavulanic acid-like gene clusters were not. The suggested in vivo involvement of an isomerase-like domain recruited by an ancestral β-lactamase related protein, supports a previous hypothesis that Cpe could be involved in a step requiring the opening and modification of the clavulanic acid core during its biosynthesis from 5S precursors.

Kinetic modelling and meta-analysis of the B. subtilis SigA regulatory network during spore germination and outgrowth

This study describes the meta-analysis and kinetic modelling of gene expression control by sigma factor SigA of Bacillus subtilis during germination and outgrowth based on microarray data from 14 time points. The analysis computationally models the direct interaction among SigA, SigA-controlled sigma factor genes (sigM, sigH, sigD, sigX), and their target genes. Of the >800 known genes in the SigA regulon, as extracted from databases, 311 genes were analysed, and 190 were confirmed by the kinetic model as being controlled by SigA. For the remaining genes, alternative regulators satisfying kinetic constraints were suggested. The kinetic analysis suggested another 214 genes as potential SigA targets. The modelling was able to (i) create a particular SigA-controlled gene expression network that is active under the conditions for which the expression time series was obtained, and where SigA is the dominant regulator, (ii) suggest new potential SigA target genes, and (iii) find other possible regulators of a given gene or suggest a new mechanism of its control by identifying a matching profile of unknown regulator(s). Selected predicted regulatory interactions were experimentally tested, thus validating the model.

Elucidation of the regulatory role of the fructose operon reveals a novel target for enhancing the NADPH supply in Corynebacterium glutamicum

The performance of Corynebacterium glutamicum cell factories producing compounds which rely heavily on NADPH has been reported to depend on the sugar being metabolized. While some aspects of this phenomenon have been elucidated, there are still many unresolved questions as to how sugar metabolism is linked to redox and to the general metabolism. We here provide new insights into the regulation of the metabolism of this important platform organism by systematically characterizing mutants carrying various lesions in the fructose operon. Initially, we found that a strain where the dedicated fructose uptake system had been inactivated (KO-ptsF) was hampered in growth on sucrose minimal medium, and suppressor mutants appeared readily. Comparative genomic analysis in conjunction with enzymatic assays revealed that suppression was linked to inactivation of the pfkB gene, encoding a fructose-1-phosphate kinase. Detailed characterization of KO-ptsF, KO-pfkB and double knock-out (DKO) derivatives revealed a strong role for sugar-phosphates, especially fructose-1-phosphate (F1P), in governing sugar as well as redox metabolism due to effects on transcriptional regulation of key genes. These findings allowed us to propose a simple model explaining the correlation between sugar phosphate concentration, gene expression and ultimately the observed phenotype. To guide us in our analysis and help us identify bottlenecks in metabolism we debugged an existing genome-scale model onto which we overlaid the transcriptome data. Based on the results obtained we managed to enhance the NADPH supply and transform the wild-type strain into delivering the highest yield of lysine ever obtained on sucrose and fructose, thus providing a good example of how regulatory mechanisms can be harnessed for bioproduction.

Moraxella catarrhalis AcrAB-OprM efflux pump contributes to antimicrobial resistance and is enhanced during cold shock response

Moraxella catarrhalis is a common pathogen of the human respiratory tract. Multidrug efflux pumps play a major role in antibiotic resistance and virulence in many Gram-negative organisms. In the present study, the role of the AcrAB-OprM efflux pump in antibiotic resistance was investigated by constructing mutants that lack the acrA, acrB, and oprM genes in M. catarrhalis strain O35E. We observed a moderate (1.5-fold) decrease in the MICs of amoxicillin and cefotaxime and a marked (4.7-fold) decrease in the MICs of clarithromycin for acrA, acrB, and oprM mutants in comparison with the wild-type O35E strain. Exposure of the M. catarrhalis strains O35E and 300 to amoxicillin triggered an increased transcription of all AcrAB-OprM pump genes, and exposure of strains O35E, 300, and 415 to clarithromycin enhanced the expression of acrA and oprM mRNA. Inactivation of the AcrAB-OprM efflux pump genes demonstrated a decreased ability to invade epithelial cells compared to the parental strain, suggesting that acrA, acrB, and oprM are required for efficient invasion of human pharyngeal epithelial cells. Cold shock increases the expression of AcrAB-OprM efflux pump genes in all three M. catarrhalis strains tested. Increased expression of AcrAB-OprM pump genes after cold shock leads to a lower accumulation of Hoechst 33342 (H33342), a substrate of AcrAB-OprM efflux pumps, indicating that cold shock results in increased efflux activity. In conclusion, the AcrAB-OprM efflux pump appears to play a role in the antibiotic resistance and virulence of M. catarrhalis and is involved in the cold shock response.

Changes in activity of metabolic and regulatory pathways during germination of S. coelicolor

BACKGROUND

Bacterial spore germination is a developmental process during which all required metabolic pathways are restored to transfer cells from their dormant state into vegetative growth. Streptomyces are soil dwelling filamentous bacteria with complex life cycle, studied mostly for they ability to synthesize secondary metabolites including antibiotics.

RESULTS

Here, we present a systematic approach that analyzes gene expression data obtained from 13 time points taken over 5.5 h of Streptomyces germination. Genes whose expression was significantly enhanced/diminished during the time-course were identified, and classified to metabolic and regulatory pathways. The classification into metabolic pathways revealed timing of the activation of specific pathways during the course of germination. The analysis also identified remarkable changes in the expression of specific sigma factors over the course of germination. Based on our knowledge of the targets of these factors, we speculate on their possible roles during germination. Among the factors whose expression was enhanced during the initial part of germination, SigE is though to manage cell wall reconstruction, SigR controls protein re-aggregation, and others (SigH, SigB, SigI, SigJ) control osmotic and oxidative stress responses.

CONCLUSIONS

From the results, we conclude that most of the metabolic pathway mRNAs required for the initial phases of germination were synthesized during the sporulation process and stably conserved in the spore. After rehydration in growth medium, the stored mRNAs are being degraded and resynthesized during first hour. From the analysis of sigma factors we conclude that conditions favoring germination evoke stress-like cell responses.

A comparison of key aspects of gene regulation in Streptomyces coelicolor and Escherichia coli using nucleotide-resolution transcription maps produced in parallel by global and differential RNA sequencing

Streptomyces coelicolor is a model for studying bacteria renowned as the foremost source of natural products used clinically. Post-genomic studies have revealed complex patterns of gene expression and links to growth, morphological development and individual genes. However, the underlying regulation remains largely obscure, but undoubtedly involves steps after transcription initiation. Here we identify sites involved in RNA processing and degradation as well as transcription within a nucleotide-resolution map of the transcriptional landscape. This was achieved by combining RNA-sequencing approaches suited to the analysis of GC-rich organisms. Escherichia coli was analysed in parallel to validate the methodology and allow comparison. Previously, sites of RNA processing and degradation had not been mapped on a transcriptome-wide scale for E. coli. Through examples, we show the value of our approach and data sets. This includes the identification of new layers of transcriptional complexity associated with several key regulators of secondary metabolism and morphological development in S. coelicolor and the identification of host-encoded leaderless mRNA and rRNA processing associated with the generation of specialized ribosomes in E. coli. New regulatory small RNAs were identified for both organisms. Overall the results illustrate the diversity in mechanisms used by different bacterial groups to facilitate and regulate gene expression.

Metabolic engineering of Corynebacterium glutamicum for L-arginine production

L-arginine is an important amino acid for diverse industrial and health product applications. Here we report the development of metabolically engineered Corynebacterium glutamicum ATCC 21831 for the production of L-arginine. Random mutagenesis is first performed to increase the tolerance of C. glutamicum to L-arginine analogues, followed by systems metabolic engineering for further strain improvement, involving removal of regulatory repressors of arginine operon, optimization of NADPH level, disruption of L-glutamate exporter to increase L-arginine precursor and flux optimization of rate-limiting L-arginine biosynthetic reactions. Fed-batch fermentation of the final strain in 5 l and large-scale 1,500 l bioreactors allows production of 92.5 and 81.2 g l(-1) of L-arginine with the yields of 0.40 and 0.35 g L-arginine per gram carbon source (glucose plus sucrose), respectively. The systems metabolic engineering strategy described here will be useful for engineering Corynebacteria strains for the industrial production of L-arginine and related products.

Escherichia coli achieves faster growth by increasing catalytic and translation rates of proteins

Regulation levels of the gene expression cascade controlling protein levels and metabolic fluxes for cells to achieve faster growth have not been elaborated in acceptable detail. Furthermore, there is need for specific growth rate (μ) dependent absolute quantitative transcriptome and proteome data to understand the molecular relationships for enabling cells to modify μ. We address these questions, for the first time, by presenting regulatory strategies for more efficient metabolism of Escherichia coli at higher μ by statistical covariance analysis of genome-wide intracellular mRNA and protein concentrations coupled to metabolic flux analysis in the steady state range of μ = 0.11-0.49 h(-1). Our analyses show dominating post-transcriptional control of protein abundances and post-translational control of flux rates. On average, E. coli achieved five-times faster growth through 3.7-fold increase of apparent catalytic rates of enzymes (kapp) and 2.5-fold increased translation rates, demonstrating the relevance of post-translational regulation for increasing flux throughput. Interestingly, pathways carrying the highest flux showed both high protein abundance and kapp values. Furthermore, co-regulation analysis of enzymatic capacities revealed tightly coupled regulatory dependencies of protein synthesis and RNA precursor synthesis, substrate utilization, biosynthetic and energy generation pathways carrying the highest flux. We also observed metabolic pathway and COG specific protein and metabolic flux control levels, protein expression costs and genome-wide principles for translation efficiency and transcription unit polarity. This work contributes to the much needed quantitative understanding of coordinated gene expression regulation and metabolic flux control. Our findings will also advance modeling and metabolic engineering of industrial strains.

Interplay of gene expression noise and ultrasensitive dynamics affects bacterial operon organization

Bacterial chromosomes are organized into polycistronic cotranscribed operons, but the evolutionary pressures maintaining them are unclear. We hypothesized that operons alter gene expression noise characteristics, resulting in selection for or against maintaining operons depending on network architecture. Mathematical models for 6 functional classes of network modules showed that three classes exhibited decreased noise and 3 exhibited increased noise with same-operon cotranscription of interacting proteins. Noise reduction was often associated with a decreased chance of reaching an ultrasensitive threshold. Stochastic simulations of the lac operon demonstrated that the predicted effects of transcriptional coupling hold for a complex network module. We employed bioinformatic analysis to find overrepresentation of noise-minimizing operon organization compared with randomized controls. Among constitutively expressed physically interacting protein pairs, higher coupling frequencies appeared at lower expression levels, where noise effects are expected to be dominant. Our results thereby suggest an important role for gene expression noise, in many cases interacting with an ultrasensitive switch, in maintaining or selecting for operons in bacterial chromosomes.

In some species, most notably bacteria, chromosomal genes are arranged into clusters called operons. In operons, the process of transcription is physically coupled: a single pass of the RNA polymerase enzyme reading that region of the chromosome simultaneously produces messenger RNA encoding multiple proteins. So far, we do not have a satisfying explanation for what evolutionary forces have maintained operons on bacterial chromosomes. We hypothesized that different types of interactions between operon-coded proteins affect how strongly operons are selected for between two genes. The proposed mechanism for this effect is that operons correlate gene expression noise, changing how it manifests in the post-translational network depending on the type of protein interaction. Mathematical models demonstrate that operons reduce noise for some types of interactions but not others. We found that operon-dependent noise reduction has an underlying dependence on surprisingly high sensitivity of the network to the ratio of proteins from each gene. Databases of genetic information show that E. coli has operons more frequently than random if operons reduce noise for the type of interaction various gene pairs have, but not otherwise. Our study thus provides an example of how the architecture of post-translational networks affects bacterial evolution.

Binding Motifs in Bacterial Gene Promoters Modulate Transcriptional Effects of Global Regulators CRP and ArcA

Bacterial gene regulation involves transcription factors (TF) that bind to DNA recognition sequences in operon promoters. These recognition sequences, many of which are palindromic, are known as regulatory elements or transcription factor binding sites (TFBS). Some TFs are global regulators that can modulate the expression of hundreds of genes. In this study we examine global regulator half-sites, where a half-site, which we shall call a binding motif (BM), is one half of a palindromic TFBS. We explore the hypothesis that the number of BMs plays an important role in transcriptional regulation, examining empirical data from transcriptional profiling of the CRP and ArcA regulons. We compare the power of BM counts and of full TFBS characteristics to predict induced transcriptional activity. We find that CRP BM counts have a nonlinear effect on CRP-dependent transcriptional activity and predict this activity better than full TFBS quality or location.

Multi-tier regulation of the streptomycete morphogenetic peptide SapB

Streptomyces coelicolor is a morphologically complex bacterium requiring the secretion of surface-active proteins to progress through its life cycle. SapB represents an important class of these biosurfactants, as illustrated by its ability to restore aerial hyphae formation when applied exogenously to developmental mutants. However, such aerial hyphae fail to sporulate, exemplifying the need to co-ordinate the timing of SapB production with other developmental events. SapB has an unusual lantibiotic structure. Its structural gene, ramS, is only 38 nucleotides downstream of the gene encoding its putative modification enzyme, RamC. Transient, co-ordinated expression of the operon was thought to be controlled by the response regulator RamR. However, we show that ramS is transcribed throughout the cell cycle with a dual expression profile dissimilar to the tightly controlled ramC expression. Surprisingly, post-translational modification relies on prior membrane localization of the precursor peptide, RamS, as demonstrated by the absence of RamS modification in S. coelicolor hyphae treated with the Bacillus subtilis lipoprotein surfactin. Our results demonstrate that interspecies interaction can also be mediated by interference of post-translational events. Further, temporal and spatial regulation of irreversible post-translational modification of a surface-active morphogenetic peptide suggests a new model for the control of key developmental events.

Bacterial transcriptomics: what is beyond the RNA horiz-ome?

Over the past 3 years, bacterial transcriptomics has undergone a massive revolution. Increased sequencing capacity and novel tools have made it possible to explore the bacterial transcriptome to an unprecedented depth, which has revealed that the transcriptome is more complex and dynamic than expected. Alternative transcripts within operons challenge the classic operon definition, and many small RNAs involved in the regulation of transcription, translation and pathogenesis have been discovered. Furthermore, mRNAs may localize to specific areas in the cell, and the spatial organization and dynamics of the chromosome have been shown to be important for transcription. Epigenetic modifications of DNA also affect transcription, and RNA processing affects translation. Therefore, transcription in bacteria resembles that in eukaryotes in terms of complexity more closely than was previously thought. Here we will discuss the contribution of 'omics' approaches to these discoveries as well as the possible impact that they are expected to have in the future.

Expression of a polycistronic messenger RNA involved in antibiotic production in an rnc mutant of Streptomyces coelicolor

RNase III is a double strand specific endoribonuclease that is involved in the regulation of gene expression in bacteria. In Streptomyces coelicolor, an RNase III (rnc) null mutant manifests decreased ability to synthesize antibiotics, suggesting that RNase III globally regulates antibiotic production in that species. As RNase III is involved in the processing of ribosomal RNAs in S. coelicolor and other bacteria, an alternative explanation for the effects of the rnc mutation on antibiotic production would involve the formation of defective ribosomes in the absence of RNase III. Those ribosomes might be unable to translate the long polycistronic messenger RNAs known to be produced by operons containing genes for antibiotic production. To examine this possibility, we have constructed a reporter plasmid whose insert encodes an operon derived from the actinorhodin cluster of S. coelicolor. We show that an rnc null mutant of S. coelicolor is capable of translating the polycistronic message transcribed from the operon. We show further that RNA species with the mobilities expected for mature 16S and 23S ribosomal RNAs are produced in the rnc mutant even though the mutant contains higher levels of the 30S rRNA precursor than the wild-type strain.

Properties of CsnR, the transcriptional repressor of the chitosanase gene, csnA, of Streptomyces lividans

A palindromic sequence is present in the intergenic region preceding the chitosanase gene csnA (SSPG_06922) of Streptomyces lividans TK24. This sequence was also found in front of putative chitosanase genes in several other actinomycete genomes and upstream genes encoding putative transcriptional regulators of the ROK family, including csnR (SSPG_04872) in S. lividans. The latter was examined as a possible transcriptional regulator (CsnR) of chitosanase gene expression. In vitro, purified CsnR bound strongly to the palindromic sequences of the csnA and csnR genes (equilibrium dissociation constant [K(D)] = 0.032 and 0.040 nM, respectively). Binding was impaired in the presence of chitosan oligosaccharides and d-glucosamine, and chitosan dimer was found to be the best effector, as determined by an equilibrium competition experiment and 50% inhibitory concentration (IC(50)) determination, while glucose, N-acetyl-glucosamine, and galactosamine had no effect. In vivo, comparison of the S. lividans wild type and ΔCsnR strains using β-lactamase reporter genes showed that CsnR represses the expression of csnA and of its own gene, which was confirmed by quantitative PCR (qPCR). CsnR is localized at the beginning of a gene cluster, possibly an operon, the organization of which is conserved through many actinomycete genomes. The CsnR-mediated chitosanase regulation mechanism seems to be widespread among actinomycetes.

Enhanced expression of recX in Mycobacterium tuberculosis owing to a promoter internal to recA

RecX is a small protein that interacts with, and modulates the activity of, RecA protein. In mycobacteria the recX gene is located immediately downstream of the recA gene, and the coding regions overlap. It has previously been shown that these two genes are co-transcribed in Mycobacterium smegmatis. In this study we examine the expression of recX in Mycobacterium tuberculosis. In addition to being co-transcribed with recA from the DNA-damage inducible recA promoters, we identify a constitutive recX promoter located within the recA coding sequence that is strong enough to make a significant contribution to the expression level of recX in the absence of DNA damage. Intriguingly, this promoter is inactivated in M. smegmatis by a critical base change in the -10 promoter motif, which probably accounts for the lower level of expression of recX relative to recA that we observed in that species. It is possible that this difference in relative expression influences RecA functions including the response to DNA damage of LexA-regulated genes.

WebGeSTer DB--a transcription terminator database

We present WebGeSTer DB, the largest database of intrinsic transcription terminators (http://pallab.serc.iisc.ernet.in/gester). The database comprises of a million terminators identified in 1060 bacterial genome sequences and 798 plasmids. Users can obtain both graphic and tabular results on putative terminators based on default or user-defined parameters. The results are arranged in different tiers to facilitate retrieval, as per the specific requirements. An interactive map has been incorporated to visualize the distribution of terminators across the whole genome. Analysis of the results, both at the whole-genome level and with respect to terminators downstream of specific genes, offers insight into the prevalence of canonical and non-canonical terminators across different phyla. The data in the database reinforce the paradigm that intrinsic termination is a conserved and efficient regulatory mechanism in bacteria. Our database is freely accessible.

Regulation of a novel gene cluster involved in secondary metabolite production in Streptomyces coelicolor

Antibiotic biosynthesis in the streptomycetes is a complex and highly regulated process. Here, we provide evidence for the contribution of a novel genetic locus to antibiotic production in Streptomyces coelicolor. The overexpression of a gene cluster comprising four protein-encoding genes (abeABCD) and an antisense RNA-encoding gene (α-abeA) stimulated the production of the blue-pigmented metabolite actinorhodin on solid medium. Actinorhodin production also was enhanced by the overexpression of an adjacent gene (abeR) encoding a predicted Streptomyces antibiotic regulatory protein (SARP), while the deletion of this gene impaired actinorhodin production. We found the abe genes to be differentially regulated and controlled at multiple levels. Upstream of abeA was a promoter that directed the transcription of abeABCD at a low but constitutive level. The expression of abeBCD was, however, significantly upregulated at a time that coincided with the initiation of aerial development and the onset of secondary metabolism; this expression was activated by the binding of AbeR to four heptameric repeats upstream of a promoter within abeA. Expressed divergently to the abeBCD promoter was α-abeA, whose expression mirrored that of abeBCD but did not require activation by AbeR. Instead, α-abeA transcript levels were subject to negative control by the double-strand-specific RNase, RNase III.

Use of an inducible promoter for antibiotic production in a heterologous host

The biosynthetic gene cluster of the aminocoumarin antibiotic novobiocin comprises 20 coding sequences. Sixteen of them code for essential enzymes for novobiocin formation, transcribed in the form of a single 18-kb polycistronic mRNA. In the present study, we replaced the genuine promoter of this operon by the tetracycline-inducible promoter tcp830 and at the same time deleting the two pathway-specific positive regulator genes of novobiocin biosynthesis. The heterologous producer Streptomyces coelicolor M512 harboring the modified gene cluster produced, upon addition of 2 mg L(-1) of the inducer compound anhydrotetracyline, 3.4-fold more novobiocin than strains carrying the unmodified cluster. A second tcp830 promoter was inserted in the middle of the 18-kb operon in order to ensure adequate transcription of the rearmost genes. However, this did not lead to a further increase of novobiocin formation, showing that a single tcp830 promoter was sufficient to achieve high transcription of all 16 genes of the operon. A high induction of novobiocin formation was achieved within a wide range of anhydrotetracyline concentrations (0.25-2.0 mg L(-1)). Growth of the strains was not affected by these concentrations. The inducer compound could be added either at the time of inoculation or at any other time up to mid-growth phase, always achieving a similar final antibiotic production. Therefore, the tcp830 promoter presents a robust, easy-to-use system for the inducible expression of biosynthetic gene clusters in heterologous hosts, independent from the genuine regulatory network.

Transcriptional regulation of the novobiocin biosynthetic gene cluster

The aminocoumarin antibiotic novobiocin is a gyrase inhibitor formed by a Streptomyces strain. The biosynthetic gene cluster of novobiocin spans 23.4 kb and contains 20 coding sequences, among them the two regulatory genes novE and novG. We investigated the location of transcriptional promoters within this cluster by insertion of transcriptional terminator cassettes and RT-PCR analysis of the resulting mutants. The cluster was found to contain eight DNA regions with promoter activity. The regulatory protein NovG binds to a previously identified binding site within the promoter region located upstream of novH, but apparently not to any of the other seven promoters. Quantitative real-time PCR was used to compare the number of transcripts in a strain carrying an intact novobiocin cluster with strains carrying mutated clusters. Both in-frame deletion of the regulatory gene novG and insertion of a terminator cassette into the biosynthetic gene novH led to a strong reduction of the number of transcripts of the genes located between novH and novW. This suggested that these 16 biosynthetic genes form a single operon. Three internal promoters are located within this operon but appear to be of minor importance, if any, under our experimental conditions. Transcription of novG was found to depend on the presence of NovE, suggesting that the two regulatory genes, novE and novG, act in a cascade-like mechanism. The resistance gene gyrB(R), encoding an aminocoumarin-resistant gyrase B subunit, may initially be co-transcribed with the genes from novH to novW. However, when the gyrase inhibitor novobiocin accumulates in the cultures, gyrB(R) is transcribed from its own promoter. Previous work has suggested that this promoter is controlled by the superhelical density of chromosomal DNA.

Inferring the transcriptional network of Bacillus subtilis

The adaptation of bacteria to the vigorous environmental changes they undergo is crucial to their survival. They achieve this adaptation partly via intricate regulation of the transcription of their genes. In this study, we infer the transcriptional network of the Gram-positive model organism, Bacillus subtilis. We use a data integration workflow, exploiting both motif and expression data, towards the generation of condition-dependent transcriptional modules. In building the motif data, we rely on both known and predicted information. Known motifs were derived from DBTBS, while predicted motifs were generated by a de novo motif detection method that utilizes comparative genomics. The expression data consists of a compendium of microarrays across different platforms. Our results indicate that a considerable part of the B. subtilis network is yet undiscovered; we could predict 417 new regulatory interactions for known regulators and 453 interactions for yet uncharacterized regulators. The regulators in our network showed a preference for regulating modules in certain environmental conditions. Also, substantial condition-dependent intra-operonic regulation seems to take place. Global regulators seem to require functional flexibility to attain their roles by acting as both activators and repressors.

A universal operon predictor for prokaryotic genomes

Identification of operons at the genome scale of prokaryotic organisms represents a key step in deciphering of their transcriptional regulation machinery, biological pathways, and networks. While numerous computational methods have been shown to be effective in predicting operons for well-studied organisms such as Escherichia coli K12 and Bacillus subtilis 168, these methods generally do not generalize well to genomes other than the ones used to train the methods, or closely related genomes because they rely on organism-specific information. Several methods have been explored to address this problem through utilizing only genomic structural information conserved across multiple organisms, but they all suffer from the issue of low prediction sensitivity. In this paper, we report a novel operon prediction method that is applicable to any prokaryotic genome with high prediction accuracy. The key idea of the method is to predict operons through identification of conserved gene clusters across multiple genomes and through deriving a key parameter relevant to the distribution of intergenic distances in genomes. We have implemented this method using a graph-theoretic approach, to calculate a set of maximum gene clusters in the target genome that are conserved across multiple reference genomes. Our computational results have shown that this method has higher prediction sensitivity as well as specificity than most of the published methods. We have carried out a preliminary study on operons unique to archaea and bacteria, respectively, and derived a number of interesting new insights about operons between these two kingdoms. The software and predicted operons of 365 prokaryotic genomes are available at http://csbl.bmb.uga.edu/~dongsheng/UNIPOP.

Occurrence, divergence and evolution of intrinsic terminators across eubacteria

In Escherichia coli, the canonical intrinsic terminator of transcription includes a palindrome followed by a U-trail on the transcript. The apparent underrepresentation of such terminators in eubacterial genomes led us to develop a rapid and accurate algorithm, GeSTer, to predict putative intrinsic terminators. Now, we have analyzed 378 genome sequences with an improved version of GeSTer. Our results indicate that the canonical E. coli type terminators are not overwhelmingly abundant in eubacteria. The atypical structures, having stem-loop structures but lacking 'U' trail, occur downstream of genes in all the analyzed genomes but different phyla show conserved preference for different types of terminators. This propensity correlates with genomic GC content and presence of the factor, Rho. 60-70% of identified terminators in all the genomes show "optimized" stem-length and DeltaG. These results provide evidence that eubacteria extensively rely on the mechanism of intrinsic termination, with a considerable divergence in their structure, positioning and prevalence. The software and detailed results for individual genomes are freely available on request.

Development and application of versatile high density microarrays for genome-wide analysis of Streptomyces coelicolor: characterization of the HspR regulon

BACKGROUND

DNA microarrays are a key resource for global analysis of genome content, gene expression and the distribution of transcription factor binding sites. We describe the development and application of versatile high density ink-jet in situ-synthesized DNA arrays for the G+C rich bacterium Streptomyces coelicolor. High G+C content DNA probes often perform poorly on arrays, yielding either weak hybridization or non-specific signals. Thus, more than one million 60-mer oligonucleotide probes were experimentally tested for sensitivity and specificity to enable selection of optimal probe sets for the genome microarrays. The heat-shock HspR regulatory system of S. coelicolor, a well-characterized repressor with a small number of known targets, was exploited to test and validate the arrays for use in global chromatin immunoprecipitation-on-chip (ChIP-chip) and gene expression analysis.

RESULTS

In addition to confirming dnaK, clpB and lon as in vivo targets of HspR, it was revealed, using a novel ChIP-chip data clustering method, that HspR also apparently interacts with ribosomal RNA (rrnD operon) and specific transfer RNA genes (the tRNAGln/tRNAGlu cluster). It is suggested that enhanced synthesis of Glu-tRNAGlu may reflect increased demand for tetrapyrrole biosynthesis following heat-shock. Moreover, it was found that heat-shock-induced genes are significantly enriched for Gln/Glu codons relative to the whole genome, a finding that would be consistent with HspR-mediated control of the tRNA species.

CONCLUSIONS

This study suggests that HspR fulfils a broader, unprecedented role in adaptation to stresses than previously recognized -- influencing expression of key components of the translational apparatus in addition to molecular chaperone and protease-encoding genes. It is envisaged that these experimentally optimized arrays will provide a key resource for systems level studies of Streptomyces biology.

Predicted transcription factor binding sites as predictors of operons in Escherichia coli and Streptomyces coelicolor

Background

As a polycistronic transcriptional unit of one or more adjacent genes, operons play a key role in regulation and function in prokaryotic biology, and a better understanding of how they are constituted and controlled is needed. Recent efforts have attempted to predict operonic status in sequenced genomes using a variety of techniques and data sources. To date, non-homology based operon prediction strategies have mainly used predicted promoters and terminators present at the extremities of transcriptional unit as predictors, with reasonable success. However, transcription factor binding sites (TFBSs), typically found upstream of the first gene in an operon, have not yet been evaluated.

Results

Here we apply a method originally developed for the prediction of TFBSs in Escherichia coli that minimises the need for prior knowledge and tests its ability to predict operons in E. coli and the 'more complex', pharmaceutically important, Streptomyces coelicolor. We demonstrate that through building genome specific TFBS position-specific-weight-matrices (PSWMs) it is possible to predict operons in E. coli and S. coelicolor with 83% and 93% accuracy respectively, using only TFBS as delimiters of operons. Additionally, the 'palindromicity' of TFBS footprint data of E. coli is characterised.

Conclusion

TFBS are proposed as novel independent features for use in prokaryotic operon prediction (whether alone or as part of a set of features) given their efficacy as operon predictors in E. coli and S. coelicolor. We also show that TFBS footprint data in E. coli generally contains inverted repeats with significantly (p < 0.05) greater palindromicity than random sequences. Consequently, the palindromicity of putative TFBSs predicted can also enhance operon predictions.

Transcriptome dynamics-based operon prediction and verification in Streptomyces coelicolor

Streptomyces spp. produce a variety of valuable secondary metabolites, which are regulated in a spatio-temporal manner by a complex network of inter-connected gene products. Using a compilation of genome-scale temporal transcriptome data for the model organism, Streptomyces coelicolor, under different environmental and genetic perturbations, we have developed a supervised machine-learning method for operon prediction in this microorganism. We demonstrate that, using features dependent on transcriptome dynamics and genome sequence, a support vector machines (SVM)-based classification algorithm can accurately classify >90% of gene pairs in a set of known operons. Based on model predictions for the entire genome, we verified the co-transcription of more than 250 gene pairs by RT-PCR. These results vastly increase the database of known operons in S. coelicolor and provide valuable information for exploring gene function and regulation to harness the potential of this differentiating microorganism for synthesis of natural products.

New pleiotropic effects of eliminating a rare tRNA from Streptomyces coelicolor, revealed by combined proteomic and transcriptomic analysis of liquid cultures

BACKGROUND

In Streptomyces coelicolor, bldA encodes the only tRNA for a rare leucine codon, UUA. This tRNA is unnecessary for growth, but is required for some aspects of secondary metabolism and morphological development. We describe a transcriptomic and proteomic analysis of the effects of deleting bldA on cellular processes during submerged culture: conditions relevant to the industrial production of antibiotics.

RESULTS

At the end of rapid growth, a co-ordinated transient up-regulation of about 100 genes, including many for ribosomal proteins, was seen in the parent strain but not the DeltabldA mutant. Increased basal levels of the signal molecule ppGpp in the mutant strain may be responsible for this difference. Transcripts or proteins from a further 147 genes classified as bldA-influenced were mostly expressed late in culture in the wild-type, though others were significantly transcribed during exponential growth. Some were involved in the biosynthesis of seven secondary metabolites; and some have probable roles in reorganising metabolism after rapid growth. Many of the 147 genes were "function unknown", and may represent unknown aspects of Streptomyces biology. Only two of the 147 genes contain a TTA codon, but some effects of bldA could be traced to TTA codons in regulatory genes or polycistronic operons. Several proteins were affected post-translationally by the bldA deletion. There was a statistically significant but weak positive global correlation between transcript and corresponding protein levels. Different technical limitations of the two approaches were a major cause of discrepancies in the results obtained with them.

CONCLUSION

Although deletion of bldA has very conspicuous effects on the gross phenotype, the bldA molecular phenotype revealed by the "dualomic" approach has shown that only about 2% of the genome is affected; but this includes many previously unknown effects at a variety of different levels, including post-translational changes in proteins and global cellular physiology.

A predicted operon map for Mycobacterium tuberculosis

The prediction of operons in Mycobacterium tuberculosis (MTB) is a first step toward understanding the regulatory network of this pathogen. Here we apply a statistical model using logistic regression to predict operons in MTB. As predictors, our model incorporates intergenic distance and the correlation of gene expression calculated for adjacent gene pairs from over 474 microarray experiments with MTB RNA. We validate our findings with known examples from the literature and experimentation. From this model, we rank each potential operon pair by the strength of evidence for cotranscription, choose a classification threshold with a true positive rate of over 90% at a false positive rate of 9.1%, and use it to construct an operon map for the MTB genome.