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.     

Germline Expression Influences Operon Organization in the Caenorhabditis elegans Genome

Operons are found across multiple kingdoms and phyla, from prokaryotes to chordates. In the nematode Caenorhabditis elegans, the genome contains >1000 operons that compose ~15% of the protein-coding genes. However, determination of the force(s) promoting the origin and maintenance of operons in C. elegans has proved elusive. Compared to bacterial operons, genes within a C. elegans operon often show poor coexpression and only sometimes encode proteins with related functions. Using analysis of microarray and large-scale in situ hybridization data, we demonstrate that almost all operon-encoded genes are expressed in germline tissue. However, genes expressed during spermatogenesis are excluded from operons. Operons group together along chromosomes in local clusters that also contain monocistronic germline-expressed genes. Additionally, germline expression of genes in operons is largely independent of the molecular function of the encoded proteins. These analyses demonstrate that mechanisms governing germline gene expression influence operon origination and/or maintenance. Thus, gene expression in a specific tissue can have profound effects on the evolution of genome organization.     

New insights into the molecular mechanisms of evolution: stress increases genetic diversity

The mechanisms of stress-induced mutagenesis in prokaryotes and realization of reserved (preaccumulated) genetic variation in eukaryotes are considered. In prokaryotes, replication becomes error-prone in stress because of the induction of the SOS response and the inactivation of the mismatch repair system; stress also increases the transposition rate and the efficiency of interspecific gene transfer. In eukaryotes, chaperone HSP90, which restores the native folding of mutant proteins (e.g., signal transduction and morphogenetic proteins) in normal conditions, fails to do so in stress, which leads to abrupt expression of multiple mutations earlier reserved in the corresponding genes. The role of these mechanisms in the evolution of prokaryotes and eukaryotes is discussed.     

New Insights into the Molecular Mechanisms of Evolution: Stress Increases Genetic Diversity

The mechanisms of stress-induced mutagenesis in prokaryotes and realization of reserved (preaccumulated) genetic variation in eukaryotes are considered. In prokaryotes, replication becomes error-prone in stress because of the induction of the SOS response and the inactivation of the mismatch repair system; stress also increases the transposition rate and the efficiency of interspecific gene transfer. In eukaryotes, chaperone HSP90, which restores the native folding of mutant proteins (e.g., signal transduction and morphogenetic proteins) in normal conditions, fails to do so in stress, which leads to abrupt expression of multiple mutations earlier reserved in the corresponding genes. The role of these mechanisms in the evolution of prokaryotes and eukaryotes is discussed.     

Sterol regulatory element-binding proteins induce an entire pathway of cholesterol synthesis

To evaluate the effects of sterol regulatory element-binding proteins (SREBPs) on the expression of the individual enzymes in the cholesterol synthetic pathway, we examined expression of these genes in the livers from wild-type and transgenic mice overexpressing nuclear SREBP-1a or -2. As estimated by a Northern blot analysis, overexpression of nuclear SREBP-1a or -2 caused marked increases in mRNA levels of the whole battery of cholesterogenic genes. This SREBP activation covers not only rate-limiting enzymes such as HMG CoA synthase and reductase that have been well established as SREBP targets, but also all the enzyme genes in the cholesterol synthetic pathway tested here. The activated genes include mevalonate kinase, mevalonate pyrophosphate decarboxylase, isopentenyl phosphate isomerase, geranylgeranyl pyrophosphate synthase, farnesyl pyrophosphate synthase, squalene synthase, squalene epoxidase, lanosterol synthase, lanosterol demethylase, and 7-dehydro-cholesterol reductase. These results demonstrate that SREBPs activate every step of cholesterol synthetic pathway, contributing to an efficient cholesterol synthesis.     

Frequent,independent transfers of a catabolic gene from bacteria to contrasted filamentous eukaryotes

Even genetically distant prokaryotes can exchange genes between them, and these horizontal gene transfer events play a central role in adaptation and evolution. While this was long thought to be restricted to prokaryotes, certain eukaryotes have acquired genes of bacterial origin. However, gene acquisitions in eukaryotes are thought to be much less important in magnitude than in prokaryotes. Here, we describe the complex evolutionary history of a bacterial catabolic gene that has been transferred repeatedly from different bacterial phyla to stramenopiles and fungi. Indeed, phylogenomic analysis pointed to multiple acquisitions of the gene in these filamentous eukaryotes—as many as 15 different events for 65 microeukaryotes. Furthermore, once transferred, this gene acquired introns and was found expressed in mRNA databases for most recipients. Our results show that effective inter-domain transfers and subsequent adaptation of a prokaryotic gene in eukaryotic cells can happen at an unprecedented magnitude.     

Regulation of respiration in enterobacteria: comparison of microarray and comparative genomic data

Microarrays are widely used for gene expression profiling. In the case of prokaryotes such arrays usually provide data about composition of modulons, groups of genes whose expression is influenced by a single regulatory system or external stimulus. Unlike modulons, regulons include only genes directly controlled by regulatory systems. Here we compared the structures of the Fnr and ArcA modulons and regulons. The data about modulon composition were taken from published microarray assays, whereas regulons were characterized using comparative genomic approaches. The Fnr and ArcA regulons were shown to contain 26 and 16 operons, respectively. Ten operons had high-score and highly conserved site for both Fnr and ArcA. These genes are the "core of regulons". Remarkably, all "core genes" encode enzymes involved in aerobic respiration and central metabolism. The Fnr-ArcA regulatory cascade plays an important role in expansion of the Fnr modulon.     

The life-cycle of operons

Operons are a major feature of all prokaryotic genomes, but how and why operon structures vary is not well understood. To elucidate the life-cycle of operons, we compared gene order between Escherichia coli K12 and its relatives and identified the recently formed and destroyed operons in E. coli. This allowed us to determine how operons form, how they become closely spaced, and how they die. Our findings suggest that operon evolution may be driven by selection on gene expression patterns. First, both operon creation and operon destruction lead to large changes in gene expression patterns. For example, the removal of lysA and ruvA from ancestral operons that contained essential genes allowed their expression to respond to lysine levels and DNA damage, respectively. Second, some operons have undergone accelerated evolution, with multiple new genes being added during a brief period. Third, although genes within operons are usually closely spaced because of a neutral bias toward deletion and because of selection against large overlaps, genes in highly expressed operons tend to be widely spaced because of regulatory fine-tuning by intervening sequences. Although operon evolution may be adaptive, it need not be optimal: new operons often comprise functionally unrelated genes that were already in proximity before the operon formed.     

Ancient Origin of the Tryptophan Operon and the Dynamics of Evolutionary Change

The seven conserved enzymatic domains required for tryptophan (Trp) biosynthesis are encoded in seven genetic regions that are organized differently (whole-pathway operons, multiple partial-pathway operons, and dispersed genes) in prokaryotes. A comparative bioinformatics evaluation of the conservation and organization of the genes of Trp biosynthesis in prokaryotic operons should serve as an excellent model for assessing the feasibility of predicting the evolutionary histories of genes and operons associated with other biochemical pathways. These comparisons should provide a better understanding of possible explanations for differences in operon organization in different organisms at a genomics level. These analyses may also permit identification of some of the prevailing forces that dictated specific gene rearrangements during the course of evolution. Operons concerned with Trp biosynthesis in prokaryotes have been in a dynamic state of flux. Analysis of closely related organisms among the Bacteria at various phylogenetic nodes reveals many examples of operon scission, gene dispersal, gene fusion, gene scrambling, and gene loss from which the direction of evolutionary events can be deduced. Two milestone evolutionary events have been mapped to the 16S rRNA tree of Bacteria, one splitting the operon in two, and the other rejoining it by gene fusion. The Archaea, though less resolved due to a lesser genome representation, appear to exhibit more gene scrambling than the Bacteria. The trp operon appears to have been an ancient innovation; it was already present in the common ancestor of Bacteria and Archaea. Although the operon has been subjected, even in recent times, to dynamic changes in gene rearrangement, the ancestral gene order can be deduced with confidence. The evolutionary history of the genes of the pathway is discernible in rough outline as a vertical line of descent, with events of lateral gene transfer or paralogy enriching the analysis as interesting features that can be distinguished. As additional genomes are thoroughly analyzed, an increasingly refined resolution of the sequential evolutionary steps is clearly possible. These comparisons suggest that present-day trp operons that possess finely tuned regulatory features are under strong positive selection and are able to resist the disruptive evolutionary events that may be experienced by simpler, poorly regulated operons.     

Ancient origin of the tryptophan operon and the dynamics of evolutionary change.

The seven conserved enzymatic domains required for tryptophan (Trp) biosynthesis are encoded in seven genetic regions that are organized differently (whole-pathway operons, multiple partial-pathway operons, and dispersed genes) in prokaryotes. A comparative bioinformatics evaluation of the conservation and organization of the genes of Trp biosynthesis in prokaryotic operons should serve as an excellent model for assessing the feasibility of predicting the evolutionary histories of genes and operons associated with other biochemical pathways. These comparisons should provide a better understanding of possible explanations for differences in operon organization in different organisms at a genomics level. These analyses may also permit identification of some of the prevailing forces that dictated specific gene rearrangements during the course of evolution. Operons concerned with Trp biosynthesis in prokaryotes have been in a dynamic state of flux. Analysis of closely related organisms among the Bacteria at various phylogenetic nodes reveals many examples of operon scission, gene dispersal, gene fusion, gene scrambling, and gene loss from which the direction of evolutionary events can be deduced. Two milestone evolutionary events have been mapped to the 16S rRNA tree of Bacteria, one splitting the operon in two, and the other rejoining it by gene fusion. The Archaea, though less resolved due to a lesser genome representation, appear to exhibit more gene scrambling than the Bacteria. The trp operon appears to have been an ancient innovation; it was already present in the common ancestor of Bacteria and Archaea. Although the operon has been subjected, even in recent times, to dynamic changes in gene rearrangement, the ancestral gene order can be deduced with confidence. The evolutionary history of the genes of the pathway is discernible in rough outline as a vertical line of descent, with events of lateral gene transfer or paralogy enriching the analysis as interesting features that can be distinguished. As additional genomes are thoroughly analyzed, an increasingly refined resolution of the sequential evolutionary steps is clearly possible. These comparisons suggest that present-day trp operons that possess finely tuned regulatory features are under strong positive selection and are able to resist the disruptive evolutionary events that may be experienced by simpler, poorly regulated operons.     

Experimental demonstration of operon formation catalyzed by insertion sequence

Operons are a hallmark of the genomic and regulatory architecture of prokaryotes. However, the mechanism by which two genes placed far apart gradually come close and form operons remains to be elucidated. Here, we propose a new model of the origin of operons: Mobile genetic elements called insertion sequences can facilitate the formation of operons by consecutive insertion–deletion–excision reactions. This mechanism barely leaves traces of insertion sequences and thus difficult to detect in nature. In this study, as a proof-of-concept, we reproducibly demonstrated operon formation in the laboratory. The insertion sequence IS3 and the insertion sequence excision enhancer are genes found in a broad range of bacterial species. We introduced these genes into insertion sequence-less Escherichia coli and found that, supporting our hypothesis, the activity of the two genes altered the expression of genes surrounding IS3, closed a 2.7 kb gap between a pair of genes, and formed new operons. This study shows how insertion sequences can facilitate the rapid formation of operons through locally increasing the structural mutation rates and highlights how coevolution with mobile elements may shape the organization of prokaryotic genomes and gene regulation.