A comparative study of the ribosomal RNA operons of Streptomyces coelicolor A3(2) and sequence analysis of rrnA

Genome mining of biosynthetic and chemotherapeutic gene clusters in Streptomyces bacteria

Streptomyces bacteria are known for their prolific production of secondary metabolites, many of which have been widely used in human medicine, agriculture and animal health. To guide the effective prioritization of specific biosynthetic gene clusters (BGCs) for drug development and targeting the most prolific producer strains, knowledge about phylogenetic relationships of Streptomyces species, genome-wide diversity and distribution patterns of BGCs is critical. We used genomic and phylogenetic methods to elucidate the diversity of major classes of BGCs in 1,110 publicly available Streptomyces genomes. Genome mining of Streptomyces reveals high diversity of BGCs and variable distribution patterns in the Streptomyces phylogeny, even among very closely related strains. The most common BGCs are non-ribosomal peptide synthetases, type 1 polyketide synthases, terpenes, and lantipeptides. We also found that numerous Streptomyces species harbor BGCs known to encode antitumor compounds. We observed that strains that are considered the same species can vary tremendously in the BGCs they carry, suggesting that strain-level genome sequencing can uncover high levels of BGC diversity and potentially useful derivatives of any one compound. These findings suggest that a strain-level strategy for exploring secondary metabolites for clinical use provides an alternative or complementary approach to discovering novel pharmaceutical compounds from microbes.

Inter- and intracellular colonization of Arabidopsis roots by endophytic actinobacteria and the impact of plant hormones on their antimicrobial activity

Many actinobacteria live in close association with eukaryotes such as fungi, insects, animals and plants. Plant-associated actinobacteria display (endo)symbiotic, saprophytic or pathogenic life styles, and can make up a substantial part of the endophytic community. Here, we characterised endophytic actinobacteria isolated from root tissue of Arabidopsis thaliana (Arabidopsis) plants grown in soil from a natural ecosystem. Many of these actinobacteria belong to the family of Streptomycetaceae with Streptomyces olivochromogenes and Streptomyces clavifer as well represented species. When seeds of Arabidopsis were inoculated with spores of Streptomyces strain coa1, which shows high similarity to S. olivochromogenes, roots were colonised intercellularly and, unexpectedly, also intracellularly. Subsequent exposure of endophytic isolates to plant hormones typically found in root and shoot tissues of Arabidopsis led to altered antibiotic production against Escherichia coli and Bacillus subtilis. Taken together, our work reveals remarkable colonization patterns of endophytic streptomycetes with specific traits that may allow a competitive advantage inside root tissue.

Genome-wide analysis of in vivo binding of the master regulator DasR in Streptomyces coelicolor identifies novel non-canonical targets

Streptomycetes produce a wealth of natural products, including over half of all known antibiotics. It was previously demonstrated that N-acetylglucosamine and secondary metabolism are closely entwined in streptomycetes. Here we show that DNA recognition by the N-acetylglucosamine-responsive regulator DasR is growth-phase dependent, and that DasR can bind to sites in the S. coelicolor genome that have no obvious resemblance to previously identified DasR-responsive elements. Thus, the regulon of DasR extends well beyond what was previously predicted and includes a large number of genes with functions far removed from N-acetylglucosamine metabolism, such as genes for small RNAs and DNA transposases. Conversely, the DasR regulon during vegetative growth largely correlates to the presence of canonical DasR-responsive elements. The changes in DasR binding in vivo following N-acetylglucosamine induction were studied in detail and a possible molecular mechanism by which the influence of DasR is extended is discussed. Discussion of DasR binding was further informed by a parallel transcriptome analysis of the respective cultures. Evidence is provided that DasR binds directly to the promoters of all genes encoding pathway-specific regulators of antibiotic production in S. coelicolor, thereby providing an exquisitely simple link between nutritional control and secondary metabolism.

An engineered strong promoter for streptomycetes

Well-characterized promoters are essential tools for metabolic engineering and synthetic biology. In Streptomyces coelicolor, the native kasOp is a temporally expressed promoter strictly controlled by two regulators, ScbR and ScbR2. In this work, first, kasOp was engineered to remove a common binding site of ScbR and ScbR2 upstream of its core region, thus generating a stronger promoter, kasOp3. Second, another ScbR binding site internal to the kasOp3 core promoter region was abolished by random mutation and screening of the mutant library to obtain the strongest promoter, kasOp* (where the asterisk is used to distinguish the engineered promoter from the native promoter). The activities of kasOp* were compared with those of two known strong promoters, ermEp* and SF14p, in three Streptomyces species. kasOp* showed the highest activity at the transcription and protein levels in all three hosts. Furthermore, relative to ermEp* and SF14p, kasOp* was shown to confer the highest actinorhodin production level when used to drive the expression of actII-ORF4 in S. coelicolor. Therefore, kasOp* is a simple and well-defined strong promoter useful for gene overexpression in streptomycetes.

Detection of potential suberinase-encoding genes in Streptomyces scabiei strains and other actinobacteria

Streptomyces scabiei causes common scab, an economically important disease of potato tubers. Some authors have previously suggested that S. scabiei penetration into host plant tissue is facilitated by secretion of esterase enzymes degrading suberin, a lipidic biopolymer of the potato periderm. In the present study, S. scabiei EF-35 showed high esterase activity in suberin-containing media. This strain also exhibited esterase activity in the presence of other biopolymers, such as lignin, cutin, or xylan, but at a much lower level. In an attempt to identify the esterases involved in suberin degradation, translated open reading frames of S. scabiei 87-22 were examined for the presence of protein sequences corresponding to extracellular esterases of S. scabiei FL1 and of the fungus Coprinopsis cinerea VTT D-041011, which have previously been shown to be produced in the presence of suberin. Two putative extracellular suberinase genes, estA and sub1, were identified. The presence of these genes in several actinobacteria was investigated by Southern blot hybridization, and both genes were found in most common-scab-inducing strains. Moreover, reverse transcription - polymerase chain reaction performed with S. scabiei EF-35 showed that estA was expressed in the presence of various biopolymers, including suberin, whereas the sub1 gene appeared to be specifically expressed in the presence of suberin and cutin.

Can whole genome analysis refine the taxonomy of the genus Rhodococcus?

The current systematics of the genus Rhodococcus is unclear, partly because many members were originally included before the application of a polyphasic taxonomic approach, central to which is the acquisition of 16S rRNA sequence data. This has resulted in the reclassification and description of many new species. Hence, the literature is replete with new species names that have not been brought together in an organized and easily interpreted form. This taxonomic confusion has been compounded by assigning many xenobiotic degrading isolates with phylogenetic positions but without formal taxonomic descriptions. In order to provide a framework for a taxonomic approach based on multiple genetic loci, a survey was undertaken of the known genome characteristics of members of the genus Rhodococcus including: (i) genetics of cell envelope biosynthesis; (ii) virulence genes; (iii) gene clusters involved in metabolic degradation and industrially relevant pathways; (iv) genetic analysis tools; (v) rapid identification of bacteria including rhodococci with specific gene RFLPs; (vi) genomic organization of rrn operons. Genes encoding virulence factors have been characterized for Rhodococcus equi and Rhodococcus fascians. Based on peptide signature comparisons deduced from gene sequences for cytochrome P-450, mono- and dioxygenases, alkane degradation, nitrile metabolism, proteasomes and desulfurization, phylogenetic relationships can be deduced for Rhodococcus erythropolis, Rhodococcus globerulus, Rhodococcus ruber and a number of undesignated Rhodococcus spp. that may distinguish the genus Rhodococcus into two further genera. The linear genome topologies that exist in some Rhodococcus species may alter a previously proposed model for the analysis of genomic fingerprinting techniques used in bacterial systematics.

First complete nucleotide sequence and heterologous gene organization of the two rRNA operons in the phytoplasma genome

Phytoplasmas are cell-wallless Gram-positive low G + C bacteria belonging to the Mollicutes that inhabit the cytoplasm of plants and insects. Although phytoplasmas possess two ribosomal RNA (rrn) operons, only one has been fully sequenced. Here, we determined the complete nucleotide sequence of both rrn operons (designated rrnA and rrnB) of onion yellows (OY) phytoplasma. Both operons have rRNA genes organized as 5'-16S-23S-5S-3' with very highly conserved sequences; the 16S, 23S, and 5S rRNA genes are 99.9, 99.8, and 99.1% identical between the two operons. However, the organization of tRNA genes in the upstream region from 16S rRNA gene and in the downstream region from 5S rRNA gene differs markedly. Several promoter candidates were detected upstream from both operons, which suggests that both operons are functional. Interestingly, both have a tRNA(Ile) gene in the 16S-23S spacer region, while the reported rrnB operon of loofah witches' broom phytoplasma does not, indicating heterogenous gene organization of rrnB within phytoplasmas. The phytoplasma tRNA gene organization is similar to that of acholeplasmas, a closely related mollicute, and different from that of mycoplasmas, another mollicute. Moreover, the organization suggests that the rrn operons were derived from that of a related nonmollicute bacterium, Bacillus subtilis. This data should shed light on the evolutionary relationships and phylogeny of the mollicutes.

Evidence that a single EF-Ts suffices for the recycling of multiple and divergent EF-Tu species in Streptomyces coelicolor A3(2) and Streptomyces ramocissimus

The tsf genes from Streptomyces coelicolor A3(2) and Streptomyces ramocissimus, encoding the guanine-nucleotide exchange factor EF-Ts, were cloned and sequenced. Streptomycetes have multiple and highly divergent EF-Tu species, with EF-Tu1 and EF-Tu3 showing only about 65% amino acid sequence identity, and yet these can apparently interact with a single EF-Ts species. tsf lies in an operon with rpsB, which encodes ribosomal protein S2. The amino acid sequence of S2 from S. coelicolor differs from most other bacterial S2 homologues in having a C-terminal extension of 70 aa residues with a highly repetitive organization, the function of which is unknown. Transcription analysis of the rpsB-tsf operon of S. coelicolor by promoter probing, nuclease S1 mapping and Northern blotting revealed that the genes give rise to a bicistronic transcript from a single promoter upstream of rpsB. An attenuator was identified in the rpsB-tsf intergenic region; it results in an approximately 2:1 ratio of rpsB vs tsf transcripts. Although tuf1, encoding the major EF-Tu, is located in the rpsL ribosomal protein operon, an additional promoter in the fus-tuf1 intergenic region leads to a significant excess of EF-Tu over ribosomes. Most amino acid residues known from the Escherichia coli crystal structure of the EF-Tu-EF-Ts complex to be directly involved in interaction between the two elongation factors are conserved between E. coli and Streptomyces. However, whenever interaction residues in the EF-Tu moiety show divergence among Streptomyces EF-Tu1, EF-Tu2 and EF-Tu3, the single Streptomyces EF-Ts exhibits compensatory substitutions of the corresponding residues. These apparently enable productive interaction to occur with all three EF-Tus.

Molecular cloning and comparative sequence analyses of rRNA operons in Streptomyces nodosus ATCC 14899

The genome of Streptomyces nodosus contains six ribosomal RNA (rRNA) operons. Four of the rRNA operons; rrnB, rrnD, rrnE and rrnF were cloned. We have completely sequenced all four operons, including a region 750 base pairs (bp) upstream of the 16S rRNA gene. The three rRNA genes present in each operon were closely linked in the order 16S-23S-5S. A sequence comparison of the four operons showed more than 99% sequence similarity between the corresponding 16S and 23S rRNA genes, and more than 97% similarity between 5S rRNA genes. The sequence differences observed between 23S rRNA genes appeared to be localized in two specific regions. Substantial sequence differences were found in the region upstream of the 16S rRNA gene as well as in the internal transcribed spacers. No tRNA gene was found in the 16S-23S spacer regions.

Roles of multiple promoters in transcription of ribosomal DNA: effects of growth conditions on precursor rRNA synthesis in mycobacteria

The roles of multiple promoters in the synthesis of rRNA under different conditions of growth were investigated, using two mycobacterial species as model organisms. When Mycobacterium smegmatis was grown under optimal conditions, its two rRNA operons contributed equally, with two promoters, one from each operon, being responsible for most transcripts. In stationary-phase growth or balanced growth under carbon starvation conditions, one operon (rrnAf) dominated and its three promoters contributed more equally to the generation of transcripts. Mycobacterium tuberculosis has a single operon with two promoters, one of which generated 80% of transcripts, at all stages of growth. We infer that each promoter functions independently according to its intrinsic strength when cells are growing slowly so that one operon with three promoters is roughly equivalent to three operons with one promoter; at high growth rates, occlusion effects reduce the efficiency of multiple promoters to that of a single promoter.

Ribosomal RNA synthesis in Streptomyces lividans under heat shock conditions

Clones containing rRNA genes were isolated from a gene library of Streptomyces lividans when RNA produced under heat shock conditions was used as a probe. Two of the clones carried entire rRNA operons rrnA and rrnF, respectively, the expression of both operons being under the control of four different promoters. At least two of the promoters were fully functional when the temperature increased from 30 to 45 degrees C, ensuring transcription of the rRNA genes under the heat shock. A third clone carried a partial rRNA operon in which expression was controlled by a main promoter that was functional at both 30 and 45 degrees C.

The actinomycete Thermobispora bispora contains two distinct types of transcriptionally active 16S rRNA genes

Here we present the first description of the presence of two distinct types of 16S rRNA genes in the genome of a (eu)bacterium, Thermobispora bispora. We cloned and determined the nucleotide sequences of all four rRNA operons of T. bispora. Sequence comparisons revealed that the genome of T. bispora contains two distinct types of 16S rRNA genes, each type consisting of two identical or nearly identical copies, and three identical copies of the 23S RNA gene. The nucleotide sequences of the two types of 16S rRNA genes differ at 98 nucleotide positions (6.4% of total nucleotides) together with six regions of deletion-insertions. None of the base substitutions or insertion-deletions corresponds to any of the approximately 600 evolutionarily invariable or rarely variable nucleotides, indicating that both genes are functional. Both types of 16S rRNA genes are transcribed and processed as determined by Northern (RNA) hybridization and reverse transcriptase-mediated PCR.

Organization and nucleotide sequence analysis of the ribosomal gene set (rrnB) from Streptomyces lividans

Streptomyces lividans (Sl) contains six ribosomal RNA (rRNA) gene sets, rrnA-F (Suzuki, Y., Ono, Y., Nagata, A. and Yamada, T. (1988) Molecular cloning and characterization of an rRNA operon in Streptomyces lividans TK21. J. Bacteriol. 170, 1631-1636). We have cloned the rrnB gene cluster. Physical mapping revealed that rrnB gene set is located on a 290 kb Asel fragment in the 11 to 12 o'clock region of the S. coelicolor A3(2) chromosome. The complete nucleotide (nt) sequence of Sl 23S rRNA has been determined. The structural gene of the Sl 23S rRNA codes for the 3108 nt RNA chain. The G+C content of the 23S rRNA is 57.3 mol%. The length of the spacer region between the 23S and 5S genes is 99 bp. Analysis of the sequences between the 16S and 23S genes and downstream of the 5S rRNA gene failed to identify any tRNA-like sequences. A secondary structure model of Sl 23 rRNA is proposed, based on the earlier published model of Gutell and Fox (Nucleic Acids Res. 16 (1988) 175-269).

Sequence heterogeneities of genes encoding 16S rRNAs in Paenibacillus polymyxa detected by temperature gradient gel electrophoresis

Sequence heterogeneities in 16S rRNA genes from individual strains of Paenibacillus polymyxa were detected by sequence-dependent separation of PCR products by temperature gradient gel electrophoresis (TGGE). A fragment of the 16S rRNA genes, comprising variable regions V6 to V8, was used as a target sequence for amplifications. PCR products from P. polymyxa (type strain) emerged as a well-defined pattern of bands in the gradient gel. Six plasmids with different inserts, individually demonstrating the migration characteristics of single bands of the pattern, were obtained by cloning the PCR products. Their sequences were analyzed as a representative sample of the total heterogeneity. An amount of 10 variant nucleotide positions in the fragment of 347 bp was observed, with all substitutions conserving the relevant secondary structures of the V6 and V8 regions in the RNA molecules. Hybridizations with specifically designed probes demonstrated different chromosomal locations of the respective rRNA genes. Amplifications of reverse-transcribed rRNA from ribosome preparations, as well as whole-cell hybridizations, revealed a predominant representation of particular sequences in ribosomes of exponentially growing laboratory cultures. Different strains of P. polymyxa showed not only remarkably differing patterns of PCR products in TGGE analysis but also discriminative whole-cell labeling with the designed oligonucleotide probes, indicating the different representation of individual sequences in active ribosomes. Our results demonstrate the usefulness of TGGE for the structural analysis of heterogeneous rRNA genes together with their expression, stress problems of the generation of meaningful data for 16S rRNA sequences and probe designs, and might have consequences for evolutionary concepts.

Sequence heterogeneity in the two 16S rRNA genes of Phormium yellow leaf phytoplasma

Phormium yellow leaf (PYL) phytoplasma causes a lethal disease of the monocotyledon, New Zealand flax (Phormium tenax). The 16S rRNA genes of PYL phytoplasma were amplified from infected flax by PCR and cloned, and the nucleotide sequences were determined. DNA sequencing and Southern hybridization analysis of genomic DNA indicated the presence of two copies of the 16S rRNA gene. The two 16S rRNA genes exhibited sequence heterogeneity in 4 nucleotide positions and could be distinguished by the restriction enzymes BpmI and BsrI. This is the first record in which sequence heterogeneity in the 16S rRNA genes of a phytoplasma has been determined by sequence analysis. A phylogenetic tree based on 16S rRNA gene sequences showed that PYL phytoplasma is most closely related to the stolbur and German grapevine yellows phytoplasmas, which form the stolbur subgroup of the aster yellows group. This phylogenetic position of PYL phytoplasma was supported by 16S/23S spacer region sequence data.

Mapping of genes involved in macromolecular synthesis on the chromosome of Streptomyces coelicolor A3(2)

The genes for the beta, beta', and seven sigma factor subunits of RNA polymerase, for elongation factors EF-Tu1 and EF-Tu3, and for six rRNA operons were mapped on the combined genetic and physical map of the Streptomyces coelicolor chromosome. Like the previously mapped tRNA genes, the RNA polymerase and rRNA genes map to scattered positions. The lack of rRNA operons in the immediate vicinity of the origin of replication (oriC) and the absence of tRNA genes in any of the rRNA operons are novel features of the Streptomyces chromosome.

Identification of cis-acting regulatory regions upstream of the rRNA operons of Rhodobacter sphaeroides

The promoter region(s) for the rRNA operons of Rhodobacter sphaeroides was identified. By utilizing both in vivo and in vitro techniques, the transcriptional start sites of all three operons were identified. Upstream of the transcriptional start, -10 and -35 promoter regions that bear little similarity to typical Escherichia coli promoter sequences were identified. In addition to the promoter sequences, probable Fis protein-binding sites were identified upstream of all three rRNA operons. Transcriptional fusions of the promoter regions from rrnA and rrnB were constructed by utilizing the reporter molecule encoded by xylE and analyzed under various growth conditions, in both a wild-type background and an rrnBC mutant background. Production of the xylE gene product (catechol 2,3-dioxygenase) was always greatest under photosynthetic growth conditions. However, the upstream region of rrnB, when fused with xylE, produced significantly more catechol 2,3-dioxygenase than did analogous regions of rrnA, suggesting that the promoters of the rrn operons differ in strength. These results were further confirmed by the study of mutant strains altered for the rrn operons either singly or in combination. Little or no expression of the xylE gene was manifest in E. coli when directed by rDNA sequences derived from R. sphaeroides.

Gene organization and primary structure of a ribosomal RNA gene cluster from Streptomyces griseus subsp. griseus

The Streptomyces griseus subsp. griseus KCTC 9080 genome contains six rRNA-encoding gene (rDNA) clusters. One rDNA cluster (rrnE), contained in an 8.7-kb BamHI fragment, was cloned and sequenced. The rDNA were arranged in the order 16S-23S-5S, and separated by small intergenic spacers. No tRNA-encoding sequences were found in the spacer regions. The lengths of the mature 16S, 23S and 5S rRNAs were 1528, 3120 and 120 nucleotides (nt), respectively. The G + C content of the gene cluster was lower than that of the chromosomal DNA. In general, the primary and secondary structures of the three rRNAs showed good agreement with those from other Streptomyces species. However, in comparison with Escherichia coli, two noticeable changes (mismatches and deletions) and two large insertions were identified in the 16S and 23S rRNAs, respectively. On the other hand, regions showing considerable heterogeneity, even within the genus Streptomyces, were found in both rRNAs. Putative primers and processing signals showing high sequence similarity to those from other Streptomyces species were located in the region upstream from the 5' end of the mature 16S rRNA. A potential hairpin loop structure reminiscent of a Rho-independent terminator was located just downstream from the 5S rRNA. A considerable degree of sequence conservation and variation within rDNA gene clusters was revealed in this study, both at the infra- and suprageneric levels.