Computational analysis of Ciona intestinalis operons

Operons are clusters of genes that are co-regulated from a common promoter. Operons are typically associated with prokaryotes, although a small number of eukaryotes have been shown to possess them. Among metazoans, operons have been extensively characterized in the nematode Caenorhabditis elegans in which ～15% of the total genes are organized into operons. The most recent genome assembly for the ascidian Ciona intestinalis placed ～20% of the genes (2909 total) into 1310 operons. The majority of these operons are composed of two genes, while the largest are composed of six. Here is reported a computational analysis of the genes that comprise the Ciona operons. Gene ontology (GO) terms were identified for about two-thirds of the operon-encoded genes. Using the extensive collection of public EST libraries, estimates of temporal patterns of gene expression were generated for the operon-encoded genes. Lastly, conservation of operons was analyzed by determining how many operon-encoded genes were present in the ascidian Ciona savignyi and whether these genes were organized in orthologous operons. Over 68% of the operon-encoded genes could be assigned one or more GO terms and 697 of the 1310 operons contained genes in which all genes had at least one GO term. Of these 697 operons, GO terms were shared by all of the genes within 146 individual operons, suggesting that most operons encode genes with unrelated functions. An analysis of operon gene expression from nine different EST libraries indicated that for 587 operons, all of the genes that comprise an individual operon were expressed together in at least one EST library, suggesting that these genes may be co-regulated. About 50% (74/146) of the operons with shared GO terms also showed evidence of gene co-regulation. Comparisons with the C. savignyi genome identified orthologs for 1907 of 2909 operon genes. About 38% (504/1310) of the operons are conserved between the two Ciona species. These results suggest that like C. elegans, operons in Ciona are comprised of a variety of genes that are not necessarily related in function. The genes in only 50% of the operons appear to be co-regulated, suggesting that more complex gene regulatory mechanisms are likely operating.     

Specific association of the CRP-cyclic AMP complex with the control region of the lactose operon of Escherichia coli K 12: a direct fluorimetric study using DNA fragments of different lengths

Specific-site binding of the cAMP . CRP complex to the control region of the lactose operon of E. coli was measured directly. All of the protein molecules did bind specifically, and the binding constant for the major CRP site was not dependent on the length (62, 219 or 301 base pairs) of the DNA fragments used. Comparing the values of the binding constant measured for the major site and for the weaker "operator" CRP site, and referring to the published "consensus sequence" derived from the known CRP sites in a series of operons, we suggest that two sub-sites support additive contribution to the total binding free energy.     

Conservation of ribosomal protein gene ordering in 16 complete genomes

The organization of ribosomal proteins in 16 prokaryotic genomes was studied as an example of comparative genome analyses of gene systems. Hypothetical ribosomal protein-containing operons were constructed. These operons also contained putative genes and other non-ribosomal genes. The correspondences among these genes across different organisms were clarified by sequence homology computations. In this way a cross tabulation of 70 ribosomal proteins genes was constructed. On average, these were organized into 9-14 operons in each genome. There were also 25 non-ribosomal or putative genes in these mainly ribosomal protein operons. Hence the table contains 95 genes in total. It was found that: (i) the conservation of the block of about 20 r-proteins in the L3 and L4 operons across almost the entire eubacteria and ar-chaebacteria is remarkable; (ii) some operons only belong to eubacteria or archaebacte-ria; (iii) although the ribosomal protein operons are highly conserved within domain, there are fine variat     

Comparative analysis of the whole set of rRNA operons between an enterohemorrhagic Escherichia coli O157:H7 Sakai strain and an Escherichia coli K-12 strain MG1655

Two primer sets for direct sequence determination of all seven rRNA operons (rrn) of Escherichia coli have been developed; one is for specific-amplification of each rrn operon and the other is for direct sequencing of the amplified operons. Using these primer sets, we determined the nucleotide sequences of seven rrn operons, including promoter and terminator regions, of an enterohemorrhagic E. coli (EHEC) O157:H7 Sakai strain. To elucidate the intercistronic or intraspecific variation of rrn operons, their sequences were compared with those for the K-12 rrn operons. The rrn genes and the internal transcribed spacer regions showed a higher similarity to each other in each strain than between the corresponding operons of the two strains. However, the degree of intercistronic homogeneity was much higher in the EHEC strain than in K-12. In contrast, promoter and terminator regions in each operons were conserved between the corresponding operons of the two strains, which exceeded intercistronic similarity.     

Conservation of ribosomal protein gene ordering in 16 complete genomes

The organization of ribosomal proteins in 16 prokaryotic genomes was studied as an example of comparative genome analyses of gene systems. Hypothetical ribosomal protein-containing operons were constructed. These operons also contained putative genes and other non-ribosomal genes. The correspondences among these genes across different organisms were clarified by sequence homology computations. In this way a cross tabulation of 70 ribosomal proteins genes was constructed. On average, these were organized into 9-14 operons in each genome. There were also 25 non-ribosomal or putative genes in these mainly ribosomal protein operons. Hence the table contains 95 genes in total. It was found that: (i) the conservation of the block of about 20 r-proteins in the L3 and L4 operons across almost the entire eubacteria and archaebacteria is remarkable; (ii) some operons only belong to eubacteria or archaebacteria; (iii) although the ribosomal protein operons are highly conserved within domain, there are fine variations in some operons across different organisms within each domain, and these variations are informative on the evolutionary relations among the organisms. This method provides a new potential for studying the origin and evolution of old species.     

Functional analysis of promoters involved in quorum sensing-based regulation of bacteriocin production in Lactobacillus

Bacteriocin production in Lactobacillus sake LTH673 involves at least four operons: a regulatory operon (sppIPKR); two operons encoding bacteriocins and their immunity proteins (sppAiA and orfX); and an operon needed for secretion (sppTE). We show here that the response regulator encoded by sppR in L. sake LTH673, as well as the homologous response regulators encoded by plnC and plnD in Lactobacillus plantarum C11, bind to characteristic repeats found in the -80 to -40 regions of spp operons. The promoters controlling bacteriocin operons are strictly regulated, and their activity is increased more than 1000-fold upon activation. Constitutive expression for the regulatory and transport operons is driven, at least in part, by promoters upstream of the -80 to -40 regions. Peak promoter activity of the regulatory and transporter operons precedes that of the two bacteriocin operons. The results reveal how promoters involved in quorum sensing-based regulation of bacteriocin production in Lactobacillus differ in strength, leakiness and timing of their activity.     

Gene regulation at-a-distance in E. coli: new insights

The number of E. coli genes/operons regulated from sites distant from the gene, though limited, steadily increases. The regulation of the ula genes, in charge of L-ascorbate utilization, as well as the negative autoregulation of the non-related lambdaCI and 186CI repressors, for efficient switching of the corresponding phages from lysogeny to lysis, are recent examples. The interaction between the two GalR dimers, separated by 114 bp, undetectable in vitro, has been genetically mapped. lac repressor-operator loops might insulate a gene and its expression from the genomic environment. The genes in charge of nitrogen assimilation sequentially react to ammonia deprivation, via an increasing intracellular NRI concentration. Other sigma54-dependent genes are activated in response to various stimuli.     

Conserved and hybrid meta-cleavage operons from PAH-degrading Burkholderia RP007.

We have compared the sequence and gene order of meta-cleavage pathway operons from alpha- and gamma-subgroups of the Proteobacteria with operons from Burkholderia sp. strain RP007 which belongs to the beta-subgroup of the Proteobacteria. Burkholderia RP007 was isolated for its ability to degrade phenanthrene and contains two meta-cleavage operons. One exhibits a comparable gene order to previously characterised gamma-subgroup Proteobacterial (Pseudomonas) meta operons, whilst the other has distinctive features present in both alpha- and gamma-subgroup Proteobacterial (Sphingomonas and Pseudomonas) meta operons. Gene sequence conservation, highlighted by examining the phylogeny of Proteobacterial catechol 2,3-dioxygenase sequences, reveals that sequences generally cluster in a manner which correlates with the taxonomic grouping of the Proteobacterial subgroup from which they originated.     

Statistical Estimate of the Total Number of Operons Specific for Bacillus subtilis Sporulation

As an alternative to exhaustive mapping, an attempt has been made to obtain a rough estimate of the total number of sporulation operons by statistical treatment. Sixteen sporulation mutants taken at random were characterized biochemically and morphologically. The mutations they contained were mapped to determine whether they fell into any one of 23 known operons. From the proportion that do so (10/16), it is calculated that the most probable number of sporulation operons is 37 (68% confidence limits of 31 and 46). If allowance is made for the fact that two of the operons apparently contain mutagenic "hot spots" and the calculation is amended accordingly, the most probable numbers of operons becomes 42 (limits 33 and 59).     

Comparison of mammalian cell entry operons of mycobacteria: in silico analysis and expression profiling

Mammalian cell entry (mce) operons, implicated in the entry of mycobacteria into host cells, are present in pathogenic and saprophytic species. It is likely that the genes in these operons have functions other than those required for entry into host cells. Using in silico analysis we have identified domains within the mce operons that might justify their occurrence in saprophytic species like Mycobacterium smegmatis. Our analysis identified in addition to the mce domain, the presence of the Ttg2B and Ttg2C domains, typical of proteins involved in transport. We have also analysed and compared the expression profile between mce operons of Mycobacterium tuberculosis, Mycobacterium bovis and M. smegmatis under different growth conditions. In case of M. smegmatis, each operon presented domain truncation for at least one gene. We observe differential expression among the operons in M. smegmatis growing under different culture conditions. Bacilli growing in nutritionally rich medium with aeration, only the mce4 operon was expressed while during stationary phase of a standing culture, all four mce operons were expressed. In M. bovis, in addition to the absence of the mce3 operon, several protein domains encoded by the other operons were truncated. We detected expression of the mce2 operon in the exponential and stationary growth phase, while the mce1 operon was only expressed in the stationary growth phase. Differential expression of mce operons and their redundancy in the genome of the majority members of mycobacteria are discussed in view of our results.     

Sequence heterogeneity of the ten rRNA operons in Clostridium perfringens

We have cloned and sequenced rRNA operons of Clostridium perfringens strain 13 and analyzed the sequence structure in view of the phylogenesis. The organism had ten copies of rRNA operons all of that comprised of 16S, 23S and 5S rDNAs except for one operon. The operons clustered around the origin of replication, ranging within one-third of the whole genome sequence as it is arranged in a circle. Seven operons were transcribed in clockwise direction, and the remaining three were transcribed in counter clockwise direction assuming that the gyrA was transcribed in clockwise direction. Two of the counter clockwise operons contained tRNA(Ile) genes between the 16S and 23S rDNAs, and the other had a tRNA(Ile) genes between the 16S and 23S rDNAs and a tRNA(Asn) gene in the place of the 5S rDNA. Microheterogeneity was found within the rRNA structural genes and spacer regions. The length of each 16S, 23S and 5S rDNA were almost identical among the ten operons, however, the intergenic spacer region of 16S-23S and 23S-5S were variable in the length depending on loci of the rRNA operons on the chromosome. Nucleotide sequences of the helix 19, helix 19a, helix 20 and helix 21 of 23S rDNA were divergent and the diversity appeared to be correlated with the loci of the rRNA operons on the chromosome.