Gene Organization: Selection, Selfishness, and Serendipity

Cross-evaluation of E. coli's operon structures via a whole-cell model suggests alternative cellular benefits for low- versus high-expressing operons

Many bacteria use operons to coregulate genes, but it remains unclear how operons benefit bacteria. We integrated E. coli's 788 polycistronic operons and 1,231 transcription units into an existing whole-cell model and found inconsistencies between the proposed operon structures and the RNA-seq read counts that the model was parameterized from. We resolved these inconsistencies through iterative, model-guided corrections to both datasets, including the correction of RNA-seq counts of short genes that were misreported as zero by existing alignment algorithms. The resulting model suggested two main modes by which operons benefit bacteria. For 86% of low-expression operons, adding operons increased the co-expression probabilities of their constituent proteins, whereas for 92% of high-expression operons, adding operons resulted in more stable expression ratios between the proteins. These simulations underscored the need for further experimental work on how operons reduce noise and synchronize both the expression timing and the quantity of constituent genes. A record of this paper's transparent peer review process is included in the supplemental information.

Position preference of essential genes in prokaryotic operons

Essential genes, which form the basis of life activities, are crucial for the survival of organisms. Essential genes tend to be located in operons, but how they are distributed in operons is still unclear for most prokaryotes. In order to clarify the general rule of position preference of essential genes in operons, an index of the average position of genes in an operon was proposed, and the distributions of essential and non-essential genes in operons in 51 bacterial genomes and two archaeal genomes were analyzed based on this new index. Consequently, essential genes were found to preferentially occupy the front positions of the operons, which tend to be expressed at higher levels.

Transcription Inhibitors with XRE DNA-Binding and Cupin Signal-Sensing Domains Drive Metabolic Diversification in Pseudomonas

Bacteria of the Pseudomonas genus, including the major human pathogen Pseudomonas aeruginosa, are known for their complex regulatory networks and high number of transcription factors, which contribute to their impressive adaptive ability. However, even in the most studied species, most of the regulators are still uncharacterized.

Transcription factors (TFs) are instrumental in the bacterial response to new environmental conditions. They can act as direct signal sensors and subsequently induce changes in gene expression leading to physiological adaptation. Here, by combining transcriptome sequencing (RNA-seq) and cistrome determination (DAP-seq), we studied a family of eight TFs in Pseudomonas aeruginosa. This family, encompassing TFs with XRE-like DNA-binding and cupin signal-sensing domains, includes the metabolic regulators ErfA, PsdR, and PauR and five so-far-unstudied TFs. The genome-wide delineation of their regulons identified 39 regulatory interactions with genes mostly involved in metabolism. We found that the XRE-cupin TFs are inhibitors of their neighboring genes, forming local, functional units encoding proteins with functions in condition-specific metabolic pathways. Growth phenotypes of isogenic mutants highlighted new roles for PauR and PA0535 in polyamines and arginine metabolism. The phylogenetic analysis of this family of regulators across the bacterial kingdom revealed a wide diversity of such metabolic regulatory modules and identified species with potentially higher metabolic versatility. Numerous genes encoding uncharacterized XRE-cupin TFs were found near metabolism-related genes, illustrating the need of further systematic characterization of transcriptional regulatory networks in order to better understand the mechanisms of bacterial adaptation to new environments.

IMPORTANCE Bacteria of the Pseudomonas genus, including the major human pathogen Pseudomonas aeruginosa, are known for their complex regulatory networks and high number of transcription factors, which contribute to their impressive adaptive ability. However, even in the most studied species, most of the regulators are still uncharacterized. With the recent advances in high-throughput sequencing methods, it is now possible to fill this knowledge gap and help the understanding of how bacteria adapt and thrive in new environments. By leveraging these methods, we provide an example of a comprehensive analysis of an entire family of transcription factors and bring new insights into metabolic and regulatory adaptation in the Pseudomonas genus.

The Divided Bacterial Genome: Structure, Function, and Evolution

Approximately 10% of bacterial genomes are split between two or more large DNA fragments, a genome architecture referred to as a multipartite genome. This multipartite organization is found in many important organisms, including plant symbionts, such as the nitrogen-fixing rhizobia, and plant, animal, and human pathogens, including the genera Brucella, Vibrio, and Burkholderia. The availability of many complete bacterial genome sequences means that we can now examine on a broad scale the characteristics of the different types of DNA molecules in a genome. Recent work has begun to shed light on the unique properties of each class of replicon, the unique functional role of chromosomal and nonchromosomal DNA molecules, and how the exploitation of novel niches may have driven the evolution of the multipartite genome. The aims of this review are to (i) outline the literature regarding bacterial genomes that are divided into multiple fragments, (ii) provide a meta-analysis of completed bacterial genomes from 1,708 species as a way of reviewing the abundant information present in these genome sequences, and (iii) provide an encompassing model to explain the evolution and function of the multipartite genome structure. This review covers, among other topics, salient genome terminology; mechanisms of multipartite genome formation; the phylogenetic distribution of multipartite genomes; how each part of a genome differs with respect to genomic signatures, genetic variability, and gene functional annotation; how each DNA molecule may interact; as well as the costs and benefits of this genome structure.

Footprints of optimal protein assembly strategies in the operonic structure of prokaryotes

In this work, we investigate optimality principles behind synthesis strategies for protein complexes using a dynamic optimization approach. We show that the cellular capacity of protein synthesis has a strong influence on optimal synthesis strategies reaching from a simultaneous to a sequential synthesis of the subunits of a protein complex. Sequential synthesis is preferred if protein synthesis is strongly limited, whereas a simultaneous synthesis is optimal in situations with a high protein synthesis capacity. We confirm the predictions of our optimization approach through the analysis of the operonic organization of protein complexes in several hundred prokaryotes. Thereby, we are able to show that cellular protein synthesis capacity is a driving force in the dissolution of operons comprising the subunits of a protein complex. Thus, we also provide a tested hypothesis explaining why the subunits of many prokaryotic protein complexes are distributed across several operons despite the presumably less precise co-regulation.

Natural selection for operons depends on genome size

In prokaryotes, genome size is associated with metabolic versatility, regulatory complexity, effective population size, and horizontal transfer rates. We therefore analyzed the covariation of genome size and operon conservation to assess the evolutionary models of operon formation and maintenance. In agreement with previous results, intraoperonic pairs of essential and of highly expressed genes are more conserved. Interestingly, intraoperonic pairs of genes are also more conserved when they encode proteins at similar cell concentrations, suggesting a role of cotranscription in diminishing the cost of waste and shortfall in gene expression. Larger genomes have fewer and smaller operons that are also less conserved. Importantly, lower conservation in larger genomes was observed for all classes of operons in terms of gene expression, essentiality, and balanced protein concentration. We reached very similar conclusions in independent analyses of three major bacterial clades (α- and β-Proteobacteria and Firmicutes). Operon conservation is inversely correlated to the abundance of transcription factors in the genome when controlled for genome size. This suggests a negative association between the complexity of genetic networks and operon conservation. These results show that genome size and/or its proxies are key determinants of the intensity of natural selection for operon organization. Our data fit better the evolutionary models based on the advantage of coregulation than those based on genetic linkage or stochastic gene expression. We suggest that larger genomes with highly complex genetic networks and many transcription factors endure weaker selection for operons than smaller genomes with fewer alternative tools for genetic regulation.

Conserved patterns in bacterial genomes: a conundrum physically tailored by evolutionary tinkering

The proper functioning of bacteria is encoded in their genome at multiple levels or scales, each of which is constrained by specific physical forces. At the smallest spatial scales, interatomic forces dictate the folding and function of proteins and nucleic acids. On longer length scales, stochastic forces emerging from the thermal jiggling of proteins and RNAs impose strong constraints on the organization of genes along chromosomes, more particularly in the context of the building of nucleoprotein complexes and the operational mode of regulatory agents. At the cellular level, transcription, replication and cell division activities generate forces that act on both the internal structure and cellular location of chromosomes. The overall result is a complex multi-scale organization of genomes that reflects the evolutionary tinkering of bacteria. The goal of this review is to highlight avenues for deciphering this complexity by focusing on patterns that are conserved among evolutionarily distant bacteria. To this end, I discuss three different organizational scales: the protein structures, the chromosomal organization of genes and the global structure of chromosomes.

Identification and functional characterization of bidirectional gene pairs and their intergenic regions in maize

Background

Bidirectional gene pairs exist as a specific form of gene organization in microorganisms and mammals as well as in model plant species, such as Arabidopsis and rice. Little is known about bidirectional gene pairs in maize, which has a large genome and is one of the most important grain crops.

Results

We conducted a genome-wide search in maize using genome sequencing results from the inbred line B73. In total, 1696 bidirectional transcript pairs were identified using a modified search model. We functionally characterized the promoter activity of the intergenic regions of most of the bidirectional transcript pairs that were expressed in embryos using a maize embryo transient expression system. A comparative study of bidirectional gene pairs performed for three monocot (Zea mays, Sorghum bicolor and Oryza sativa) and two dicot (Arabidopsis thaliana and Glycine max) plant genomes showed that bidirectional gene pairs were abundant in the five plant species. Orthologous bidirectional gene pairs were clearly distinguishable between the monocot and dicot species although the total numbers of orthologous bidirectional genes were similar. Analysis of the gene pairs using the Blast2GO software suite showed that the molecular functions (MF), cellular components (CC), and biological processes (BP) associated with the bidirectional transcripts were similar among the five plant species.

Conclusions

The evolutionary analysis of the function and structure of orthologous bidirectional gene pairs in various plant species revealed a potential pathway of their origin, which may be required for the evolution of a new species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-338) contains supplementary material, which is available to authorized users.

The precarious prokaryotic chromosome

Evolutionary selection for optimal genome preservation, replication, and expression should yield similar chromosome organizations in any type of cells. And yet, the chromosome organization is surprisingly different between eukaryotes and prokaryotes. The nuclear versus cytoplasmic accommodation of genetic material accounts for the distinct eukaryotic and prokaryotic modes of genome evolution, but it falls short of explaining the differences in the chromosome organization. I propose that the two distinct ways to organize chromosomes are driven by the differences between the global-consecutive chromosome cycle of eukaryotes and the local-concurrent chromosome cycle of prokaryotes. Specifically, progressive chromosome segregation in prokaryotes demands a single duplicon per chromosome, while other "precarious" features of the prokaryotic chromosomes can be viewed as compensations for this severe restriction.

A global analysis of adaptive evolution of operons in cyanobacteria

Operons are an important feature of prokaryotic genomes. Evolution of operons is hypothesized to be adaptive and has contributed significantly towards coordinated optimization of functions. Two conflicting theories, based on (i) in situ formation to achieve co-regulation and (ii) horizontal gene transfer of functionally linked gene clusters, are generally considered to explain why and how operons have evolved. Furthermore, effects of operon evolution on genomic traits such as intergenic spacing, operon size and co-regulation are relatively less explored. Based on the conservation level in a set of diverse prokaryotes, we categorize the operonic gene pair associations and in turn the operons as ancient and recently formed. This allowed us to perform a detailed analysis of operonic structure in cyanobacteria, a morphologically and physiologically diverse group of photoautotrophs. Clustering based on operon conservation showed significant similarity with the 16S rRNA-based phylogeny, which groups the cyanobacterial strains into three clades. Clade C, dominated by strains that are believed to have undergone genome reduction, shows a larger fraction of operonic genes that are tightly packed in larger sized operons. Ancient operons are in general larger, more tightly packed, better optimized for co-regulation and part of key cellular processes. A sub-clade within Clade B, which includes Synechocystis sp. PCC 6803, shows a reverse trend in intergenic spacing. Our results suggest that while in situ formation and vertical descent may be a dominant mechanism of operon evolution in cyanobacteria, optimization of intergenic spacing and co-regulation are part of an ongoing process in the life-cycle of operons.

An effort to make sense of antisense transcription in bacteria

Analysis of bacterial transcriptomes have shown the existence of a genome-wide process of overlapping transcription due to the presence of antisense RNAs, as well as mRNAs that overlapped in their entire length or in some portion of the 5′- and 3′-UTR regions. The biological advantages of such overlapping transcription are unclear but may play important regulatory roles at the level of transcription, RNA stability and translation. In a recent report, the human pathogen Staphylococcus aureus is observed to generate genome-wide overlapping transcription in the same bacterial cells leading to a collection of short RNA fragments generated by the endoribonuclease III, RNase III. This processing appears most prominently in Gram-positive bacteria. The implications of both the use of pervasive overlapping transcription and the processing of these double stranded templates into short RNAs are explored and the consequences discussed.

Genomics of the proteorhodopsin-containing marine flavobacterium Dokdonia sp. strain MED134

Proteorhodopsin phototrophy is expected to have considerable impact on the ecology and biogeochemical roles of marine bacteria. However, the genetic features contributing to the success of proteorhodopsin-containing bacteria remain largely unknown. We investigated the genome of Dokdonia sp. strain MED134 (Bacteroidetes) for features potentially explaining its ability to grow better in light than darkness. MED134 has a relatively high number of peptidases, suggesting that amino acids are the main carbon and nitrogen sources. In addition, MED134 shares with other environmental genomes a reduction in gene copies at the expense of important ones, like membrane transporters, which might be compensated by the presence of the proteorhodopsin gene. The genome analyses suggest Dokdonia sp. MED134 is able to respond to light at least partly due to the presence of a strong flavobacterial consensus promoter sequence for the proteorhodopsin gene. Moreover, Dokdonia sp. MED134 has a complete set of anaplerotic enzymes likely to play a role in the adaptation of the carbon anabolism to the different sources of energy it can use, including light or various organic matter compounds. In addition to promoting growth, proteorhodopsin phototrophy could provide energy for the degradation of complex or recalcitrant organic matter, survival during periods of low nutrients, or uptake of amino acids and peptides at low concentrations. Our analysis suggests that the ability to harness light potentially makes MED134 less dependent on the amount and quality of organic matter or other nutrients. The genomic features reported here may well be among the keys to a successful photoheterotrophic lifestyle.

Metazoan operons accelerate recovery from growth-arrested states

Existing theories explain why operons are advantageous in prokaryotes, but their occurrence in metazoans is an enigma. Nematode operon genes, typically consisting of growth genes, are significantly upregulated during recovery from growth-arrested states. This expression pattern is anticorrelated to nonoperon genes, consistent with a competition for transcriptional resources. We find that transcriptional resources are initially limiting during recovery and that recovering animals are highly sensitive to any additional decrease in transcriptional resources. We provide evidence that operons become advantageous because, by clustering growth genes into operons, fewer promoters compete for the limited transcriptional machinery, effectively increasing the concentration of transcriptional resources and accelerating recovery. Mathematical modeling reveals how a moderate increase in transcriptional resources can substantially enhance transcription rate and recovery. This design principle occurs in different nematodes and the chordate C. intestinalis. As transition from arrest to rapid growth is shared by many metazoans, operons could have evolved to facilitate these processes.

The operonic location of auto-transcriptional repressors is highly conserved in bacteria

Bacterial genes are commonly encoded in clusters, known as operons, which share transcriptional regulatory control and often encode functionally related proteins that take part in certain biological pathways. Operons that are coregulated are known to colocalize in the genome, suggesting that their spatial organization is under selection for efficient expression regulation. However, the internal order of genes within operons is believed to be poorly conserved, and hence expression requirements are claimed to be too weak to oppose gene rearrangements. In light of these opposing views, we set out to investigate whether the internal location of the regulatory genes within operons is under selection. Our analysis shows that transcription factors (TFs) are preferentially encoded as either first or last in their operons, in the two diverged model bacteria Escherichia coli and Bacillus subtilis. In a higher resolution, we find that TFs that repress transcription of the operon in which they are encoded (autorepressors), contribute most of this signal by specific preference of the first operon position. We show that this trend is strikingly conserved throughout highly diverged bacterial phyla. Moreover, these autorepressors regulate operons that carry out highly diverse biological functions. We propose a model according to which autorepressors are selected to be located first in their operons in order to optimize transcription regulation. Specifically, the first operon position helps autorepressors to minimize leaky transcription of the operon structural genes, thus minimizing energy waste. Our analysis provides statistically robust evidence for a paradigm of bacterial autorepressor preferential operonic location. Corroborated with our suggested model, an additional layer of operon expression control that is common throughout the bacterial domain is revealed.

Allostery and intrinsic disorder mediate transcription regulation by conditional cooperativity

Regulation of the phd/doc toxin-antitoxin operon involves the toxin Doc as co- or derepressor depending on the ratio between Phd and Doc, a phenomenon known as conditional cooperativity. The mechanism underlying this observed behavior is not understood. Here we show that monomeric Doc engages two Phd dimers on two unrelated binding sites. The binding of Doc to the intrinsically disordered C-terminal domain of Phd structures its N-terminal DNA-binding domain, illustrating allosteric coupling between highly disordered and highly unstable domains. This allosteric effect also couples Doc neutralization to the conditional regulation of transcription. In this way, higher levels of Doc tighten repression up to a point where the accumulation of toxin triggers the production of Phd to counteract its action. Our experiments provide the basis for understanding the mechanism of conditional cooperative regulation of transcription typical of toxin-antitoxin modules. This model may be applicable for the regulation of other biological systems.

Origins and diversification of a complex signal transduction system in prokaryotes

The molecular machinery that controls chemotaxis in bacteria is substantially more complex than any other signal transduction system in prokaryotes, and its origins and variability among living species are unknown. We found that this multiprotein "chemotaxis system" is present in most prokaryotic species and evolved from simpler two-component regulatory systems that control prokaryotic transcription. We discovered, through genomic analysis, signaling systems intermediate between two-component systems and chemotaxis systems. Evolutionary genomics established central and auxiliary components of the chemotaxis system. While tracing its evolutionary history, we also developed a classification scheme that revealed more than a dozen distinct classes of chemotaxis systems, enabling future predictive modeling of chemotactic behavior in unstudied species.

Comparative analysis of sequence periodicity among prokaryotic genomes points to differences in nucleoid structure and a relationship to gene expression

Regular spacing of short runs of A or T nucleotides in DNA sequences with a period close to the helical period of the DNA double helix has been associated with intrinsic DNA bending and nucleosome positioning in eukaryotes. Analogous periodic signals were also observed in prokaryotic genomes. While the exact role of this periodicity in prokaryotes is not known, it has been proposed to facilitate the DNA packaging in the prokaryotic nucleoid and/or to promote negative or positive supercoiling. We developed a methodology for assessments of intragenomic heterogeneity of these periodic patterns and applied it in analysis of 1,025 prokaryotic chromosomes. This technique allows more detailed analysis of sequence periodicity than previous methods where sequence periodicity was assessed in an integral form across the whole chromosome. We found that most genomes have the periodic signal confined to several chromosomal segments while most of the chromosome lacks a strong sequence periodicity. Moreover, there are significant differences among different prokaryotes in both the intensity and persistency of sequence periodicity related to DNA curvature. We proffer that the prokaryotic nucleoid consists of relatively rigid sections stabilized by short intrinsically bent DNA segments and characterized by locally strong periodic patterns alternating with regions featuring a weak periodic signal, which presumably permits higher structural flexibility. This model applies to most bacteria and archaea. In genomes with an exceptionally persistent periodic signal, highly expressed genes tend to concentrate in aperiodic sections, suggesting that structural heterogeneity of the nucleoid is related to local differences in transcriptional activity.

The Drosophila gene CheB42a is a novel modifier of Deg/ENaC channel function

Degenerin/epithelial Na⁺ channels (DEG/ENaC) represent a diverse family of voltage-insensitive cation channels whose functions include Na⁺ transport across epithelia, mechanosensation, nociception, salt sensing, modification of neurotransmission, and detecting the neurotransmitter FMRFamide. We previously showed that the Drosophila melanogaster Deg/ENaC gene lounge lizard (llz) is co-transcribed in an operon-like locus with another gene of unknown function, CheB42a. Because operons often encode proteins in the same biochemical or physiological pathway, we hypothesized that CHEB42A and LLZ might function together. Consistent with this hypothesis, we found both genes expressed in cells previously implicated in sensory functions during male courtship. Furthermore, when coexpressed, LLZ coprecipitated with CHEB42A, suggesting that the two proteins form a complex. Although LLZ expressed either alone or with CHEB42A did not generate ion channel currents, CHEB42A increased current amplitude of another DEG/ENaC protein whose ligand (protons) is known, acid-sensing ion channel 1a (ASIC1a). We also found that CHEB42A was cleaved to generate a secreted protein, suggesting that CHEB42A may play an important role in the extracellular space. These data suggest that CHEB42A is a modulatory subunit for sensory-related Deg/ENaC signaling. These results are consistent with operon-like transcription of CheB42a and llz and explain the similar contributions of these genes to courtship behavior.

Conditions for the evolution of gene clusters in bacterial genomes

Genes encoding proteins in a common pathway are often found near each other along bacterial chromosomes. Several explanations have been proposed to account for the evolution of these structures. For instance, natural selection may directly favour gene clusters through a variety of mechanisms, such as increased efficiency of coregulation. An alternative and controversial hypothesis is the selfish operon model, which asserts that clustered arrangements of genes are more easily transferred to other species, thus improving the prospects for survival of the cluster. According to another hypothesis (the persistence model), genes that are in close proximity are less likely to be disrupted by deletions. Here we develop computational models to study the conditions under which gene clusters can evolve and persist. First, we examine the selfish operon model by re-implementing the simulation and running it under a wide range of conditions. Second, we introduce and study a Moran process in which there is natural selection for gene clustering and rearrangement occurs by genome inversion events. Finally, we develop and study a model that includes selection and inversion, which tracks the occurrence and fixation of rearrangements. Surprisingly, gene clusters fail to evolve under a wide range of conditions. Factors that promote the evolution of gene clusters include a low number of genes in the pathway, a high population size, and in the case of the selfish operon model, a high horizontal transfer rate. The computational analysis here has shown that the evolution of gene clusters can occur under both direct and indirect selection as long as certain conditions hold. Under these conditions the selfish operon model is still viable as an explanation for the evolution of gene clusters.

Genes involved in a common pathway or function are frequently found near each other on bacterial chromosomes. A number of hypotheses have been previously presented to explain this observation. A particularly influential theory is the selfish operon model, which posits that horizontal transfer could promote gene clustering by favouring transfer of arrangements of genes that are close together. Subsequent theoretical development and analysis of genomic data have contributed to the debate about the plausibility of this model. Here, by re-examining the evolutionary dynamics of gene clusters, we provide and discuss conditions under which gene clusters can evolve. We find that first, some form of bias for clustering is required for clusters to evolve. This bias can be in the form of bias in horizontal transfer towards genes that are close together, or direct natural selection for gene proximity. Our computational work does not present a theoretical obstacle to the selfish operon model as a possible explanation for the evolution of gene clusters.

Relationship between operon preference and functional properties of persistent genes in bacterial genomes

BACKGROUND

Genes in bacteria may be organised into operons, leading to strict co-expression of the genes that participate in the same operon. However, comparisons between different bacterial genomes have shown that much of the operon structure is dynamic on an evolutionary time scale. This indicates that there are opposing effects influencing the tendency for operon formation, and these effects may be reflected in properties like evolutionary rate, complex formation, metabolic pathways and gene fusion.

RESULTS

We have used multi-species protein-protein comparisons to generate a high-quality set of genes that are persistent in bacterial genomes (i.e. they have close to universal distribution). We have analysed these genes with respect to operon participation and important functional properties, including evolutionary rate and protein-protein interactions.

CONCLUSIONS

Genes for ribosomal proteins show a very slow rate of evolution. This is consistent with a strong tendency for the genes to participate in operons and for their proteins to be involved in essential and well defined complexes. Persistent genes for non-ribosomal proteins can be separated into two classes according to tendency to participate in operons. Those with a strong tendency for operon participation make proteins with fewer interaction partners that seem to participate in relatively static complexes and possibly linear pathways. Genes with a weak tendency for operon participation tend to produce proteins with more interaction partners, but possibly in more dynamic complexes and convergent pathways. Genes that are not regulated through operons are therefore more evolutionary constrained than the corresponding operon-associated genes and will on average evolve more slowly.

Conjugative plasmids: vessels of the communal gene pool

Comparative whole-genome analyses have demonstrated that horizontal gene transfer (HGT) provides a significant contribution to prokaryotic genome innovation. The evolution of specific prokaryotes is therefore tightly linked to the environment in which they live and the communal pool of genes available within that environment. Here we use the term supergenome to describe the set of all genes that a prokaryotic 'individual' can draw on within a particular environmental setting. Conjugative plasmids can be considered particularly successful entities within the communal pool, which have enabled HGT over large taxonomic distances. These plasmids are collections of discrete regions of genes that function as 'backbone modules' to undertake different aspects of overall plasmid maintenance and propagation. Conjugative plasmids often carry suites of 'accessory elements' that contribute adaptive traits to the hosts and, potentially, other resident prokaryotes within specific environmental niches. Insight into the evolution of plasmid modules therefore contributes to our knowledge of gene dissemination and evolution within prokaryotic communities. This communal pool provides the prokaryotes with an important mechanistic framework for obtaining adaptability and functional diversity that alleviates the need for large genomes of specialized 'private genes'.

Role of translational coupling in robustness of bacterial chemotaxis pathway

Evolutionary selection for robustness of signaling output in the face of stochastic variations in protein expression may explain the organization of bacterial chemotaxis genes.

Chemotaxis allows bacteria to colonize their environment more efficiently and to find optimal growth conditions, and is consequently under strong evolutionary selection. Theoretical and experimental analyses of bacterial chemotaxis suggested that the pathway has been evolutionarily optimized to produce robust output under conditions of such physiological perturbations as stochastic intercellular variations in protein levels while at the same time minimizing complexity and cost of protein expression. Pathway topology in Escherichia coli apparently evolved to produce an invariant output under concerted variations in protein levels, consistent with experimentally observed transcriptional coupling of chemotaxis genes. Here, we show that the pathway robustness is further enhanced through the pairwise translational coupling of adjacent genes. Computer simulations predicted that the robustness of the pathway against the uncorrelated variations in protein levels can be enhanced by a selective pairwise coupling of individual chemotaxis genes on one mRNA, with the order of genes in E. coli ranking among the best in terms of noise compensation. Translational coupling between chemotaxis genes was experimentally confirmed, and coupled expression of these genes was shown to improve chemotaxis. Bioinformatics analysis further revealed that E. coli gene order corresponds to consensus in sequenced bacterial genomes, confirming evolutionary selection for noise reduction. Since polycistronic gene organization is common in bacteria, translational coupling between adjacent genes may provide a general mechanism to enhance robustness of their signaling and metabolic networks. Moreover, coupling between expression of neighboring genes is also present in eukaryotes, and similar principles of noise reduction might thus apply to all cellular networks.

All cellular networks are subject to fluctuations in the levels of their components. Robustness of the network output in the face of stochastic gene expression, or gene expression noise, is therefore essential to ensure proper function. Selection for robustness might thus have shaped much of the cellular evolution. We have used Escherichia coli chemotaxis, one of the most thoroughly studied model systems for signal transduction, to analyze the role of gene organization in robustness. Our mathematical modeling predicted that coupling the expression of chemotaxis proteins with opposing functions should buffer the output of the signaling pathway against stochastic variations in protein production. Consistent with this model, protein coexpression was indeed observed to improve chemotaxis and to be under selection during chemotaxis-driven spreading of a cell population. We show that tight coexpression is ensured by both transcriptional and translational gene coupling. We conclude that evolutionary selection for pathway robustness in the presence of gene expression noise can explain, not only the polycistronic organization of chemotaxis genes, but also the gene order within chemotaxis operons. Selection on the gene order was further confirmed by the observation of a strong bias towards specific pairwise occurrences of chemotaxis genes in sequenced prokaryotic genomes.

Darwinian evolution in the light of genomics

Comparative genomics and systems biology offer unprecedented opportunities for testing central tenets of evolutionary biology formulated by Darwin in the Origin of Species in 1859 and expanded in the Modern Synthesis 100 years later. Evolutionary-genomic studies show that natural selection is only one of the forces that shape genome evolution and is not quantitatively dominant, whereas non-adaptive processes are much more prominent than previously suspected. Major contributions of horizontal gene transfer and diverse selfish genetic elements to genome evolution undermine the Tree of Life concept. An adequate depiction of evolution requires the more complex concept of a network or ‘forest’ of life. There is no consistent tendency of evolution towards increased genomic complexity, and when complexity increases, this appears to be a non-adaptive consequence of evolution under weak purifying selection rather than an adaptation. Several universals of genome evolution were discovered including the invariant distributions of evolutionary rates among orthologous genes from diverse genomes and of paralogous gene family sizes, and the negative correlation between gene expression level and sequence evolution rate. Simple, non-adaptive models of evolution explain some of these universals, suggesting that a new synthesis of evolutionary biology might become feasible in a not so remote future.

The organization of the bacterial genome

Many bacterial cellular processes interact intimately with the chromosome. Such interplay is the major driving force of genome structure or organization. Interactions take place at different scales-local for gene expression, global for replication-and lead to the differentiation of the chromosome into organizational units such as operons, replichores, or macrodomains. These processes are intermingled in the cell and create complex higher-level organizational features that are adaptive because they favor the interplay between the processes. The surprising result of selection for genome organization is that gene repertoires change much more quickly than chromosomal structure. Comparative genomics and experimental genomic manipulations are untangling the different cellular and evolutionary mechanisms causing such resilience to change. Since organization results from cellular processes, a better understanding of chromosome organization will help unravel the underlying cellular processes and their diversity.

Evolution of genome architecture

Charles Darwin believed that all traits of organisms have been honed to near perfection by natural selection. The empirical basis underlying Darwin's conclusions consisted of numerous observations made by him and other naturalists on the exquisite adaptations of animals and plants to their natural habitats and on the impressive results of artificial selection. Darwin fully appreciated the importance of heredity but was unaware of the nature and, in fact, the very existence of genomes. A century and a half after the publication of the "Origin", we have the opportunity to draw conclusions from the comparisons of hundreds of genome sequences from all walks of life. These comparisons suggest that the dominant mode of genome evolution is quite different from that of the phenotypic evolution. The genomes of vertebrates, those purported paragons of biological perfection, turned out to be veritable junkyards of selfish genetic elements where only a small fraction of the genetic material is dedicated to encoding biologically relevant information. In sharp contrast, genomes of microbes and viruses are incomparably more compact, with most of the genetic material assigned to distinct biological functions. However, even in these genomes, the specific genome organization (gene order) is poorly conserved. The results of comparative genomics lead to the conclusion that the genome architecture is not a straightforward result of continuous adaptation but rather is determined by the balance between the selection pressure, that is itself dependent on the effective population size and mutation rate, the level of recombination, and the activity of selfish elements. Although genes and, in many cases, multigene regions of genomes possess elaborate architectures that ensure regulation of expression, these arrangements are evolutionarily volatile and typically change substantially even on short evolutionary scales when gene sequences diverge minimally. Thus, the observed genome architectures are, mostly, products of neutral processes or epiphenomena of more general selective processes, such as selection for genome streamlining in successful lineages with large populations. Selection for specific gene arrangements (elements of genome architecture) seems only to modulate the results of these processes.

Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

The first bacterial genome was sequenced in 1995, and the first archaeal genome in 1996. Soon after these breakthroughs, an exponential rate of genome sequencing was established, with a doubling time of approximately 20 months for bacteria and approximately 34 months for archaea. Comparative analysis of the hundreds of sequenced bacterial and dozens of archaeal genomes leads to several generalizations on the principles of genome organization and evolution. A crucial finding that enables functional characterization of the sequenced genomes and evolutionary reconstruction is that the majority of archaeal and bacterial genes have conserved orthologs in other, often, distant organisms. However, comparative genomics also shows that horizontal gene transfer (HGT) is a dominant force of prokaryotic evolution, along with the loss of genetic material resulting in genome contraction. A crucial component of the prokaryotic world is the mobilome, the enormous collection of viruses, plasmids and other selfish elements, which are in constant exchange with more stable chromosomes and serve as HGT vehicles. Thus, the prokaryotic genome space is a tightly connected, although compartmentalized, network, a novel notion that undermines the ‘Tree of Life’ model of evolution and requires a new conceptual framework and tools for the study of prokaryotic evolution.

GadX/GadW-dependent regulation of the Escherichia coli acid fitness island: transcriptional control at the gadY-gadW divergent promoters and identification of four novel 42 bp GadX/GadW-specific binding sites

Escherichia coli has the remarkable ability to resist severe acid stress for several hours. With the notable exception of the gadBC operon, the most important genes involved in acid resistance are present within the acid fitness island (AFI), a 15 kb H-NS-repressed and RpoS-controlled genome region. The AraC/XylS-like transcriptional regulators GadX and GadW are also encoded within this region. In this article, we show that gadW transcription occurs from two native promoters, which are affected by the transcription of the divergently transcribed and GadX-dependent gadY small RNA, and from the gadX promoter. The gadXW dicistronic transcript is subjected to post-transcriptional processing in which GadY is involved. In contrast, gadW transcription negatively affects gadY transcription. By aligning the GadX/GadW binding site on the gadY promoter with the GadX/GadW binding sites previously identified in the gadA and gadBC 5' regulatory regions, we generated a 42 bp GadX/GadW consensus sequence. DNase I footprinting analyses confirmed that a 42 bp GadX/GadW binding site, which matched the consensus sequence 5'-WANDNCTDWTWKTRAYATWAWMATG KCTGATNTTTWYNTYAK-3', is also present in the regulatory region of the slp-yhiF, hdeAB and gadE-mtdEF operons, all of which belong to the AFI. The presence of five GadX/GadW-specific binding sites in the AFI suggests that GadX and GadW may act as H-NS counter-silencers.

Multi-gene gateway clone design for expression of multiple heterologous genes in living cells: eukaryotic clones containing two and three ORF multi-gene cassettes expressed from a single promoter

Two types of eukaryotic operon-type Expression clones were constructed using the Multisite Gateway system employing six types of att signals. These clones harbored a DNA cassette containing two heterologous ORFs (cDNAs) or three heterologous ORFs in tandem downstream of a single promoter. The most promoter-proximal ORF was translated via a Kozak signal and the downstream one or two ORF(s) were translated as directed by internal ribosome entry site(s) (IRES). These clones were observed to produce two or three different proteins at levels that depended on the activities of the translational initiation signals used. With the intention of modulating the expression level of the first ORF, the translational initiation signals including a Kozak sequence and 11 different IRESs were investigated for their efficiency using a single ORF. The translational activity of these signals varied within a 10-fold magnitude. Using these results, expression at pre-described relative levels was achieved from the optional IRES of the respective ORFs in the cassette. Controllable expression at desired levels of two different ORFs directed by optional IRESs on a bicistronic construct, transcribed from a single promoter, was demonstrated.

Comparative genome analysis of Lactobacillus reuteri and Lactobacillus fermentum reveal a genomic island for reuterin and cobalamin production

Lactobacillus reuteri is a heterofermentative lactic acid bacterium that naturally inhabits the gut of humans and other animals. The probiotic effects of L. reuteri have been proposed to be largely associated with the production of the broad-spectrum antimicrobial compound reuterin during anaerobic metabolism of glycerol. We determined the complete genome sequences of the reuterin-producing L. reuteri JCM 1112(T) and its closely related species Lactobacillus fermentum IFO 3956. Both are in the same phylogenetic group within the genus Lactobacillus. Comparative genome analysis revealed that L. reuteri JCM 1112(T) has a unique cluster of 58 genes for the biosynthesis of reuterin and cobalamin (vitamin B(12)). The 58-gene cluster has a lower GC content and is apparently inserted into the conserved region, suggesting that the cluster represents a genomic island acquired from an anomalous source. Two-dimensional nuclear magnetic resonance (2D-NMR) with (13)C(3)-glycerol demonstrated that L. reuteri JCM 1112(T) could convert glycerol to reuterin in vivo, substantiating the potential of L. reuteri JCM 1112(T) to produce reuterin in the intestine. Given that glycerol is shown to be naturally present in feces, the acquired ability to produce reuterin and cobalamin is an adaptive evolutionary response that likely contributes to the probiotic properties of L. reuteri.

Spatial effects on the speed and reliability of protein-DNA search

Strong experimental and theoretical evidence shows that transcription factors (TFs) and other specific DNA-binding proteins find their sites using a two-mode search: alternating between three-dimensional (3D) diffusion through the cell and one-dimensional (1D) sliding along the DNA. We show that, due to the 1D component of the search process, the search time of a TF can depend on the initial position of the TF. We formalize this effect by discriminating between two types of searches: global and local. Using analytical calculations and simulations, we estimate how close a TF and binding site need to be to make a local search likely. We then use our model to interpret the wide range of experimental measurements of this parameter. We also show that local and global searches differ significantly in average search time and the variability of search time. These results lead to a number of biological implications, including suggestions of how prokaryotes achieve rapid gene regulation and the relationship between the search mechanism and noise in gene expression. Lastly, we propose a number of experiments to verify the existence and quantify the extent of spatial effects on the TF search process in prokaryotes.

Complete genome sequence of the myxobacterium Sorangium cellulosum

The genus Sorangium synthesizes approximately half of the secondary metabolites isolated from myxobacteria, including the anti-cancer metabolite epothilone. We report the complete genome sequence of the model Sorangium strain S. cellulosum So ce56, which produces several natural products and has morphological and physiological properties typical of the genus. The circular genome, comprising 13,033,779 base pairs, is the largest bacterial genome sequenced to date. No global synteny with the genome of Myxococcus xanthus is apparent, revealing an unanticipated level of divergence between these myxobacteria. A large percentage of the genome is devoted to regulation, particularly post-translational phosphorylation, which probably supports the strain's complex, social lifestyle. This regulatory network includes the highest number of eukaryotic protein kinase-like kinases discovered in any organism. Seventeen secondary metabolite loci are encoded in the genome, as well as many enzymes with potential utility in industry.

Chromatin remodelling is a major source of coexpression of linked genes in yeast

In diverse organisms, neighbouring genes in the genome tend to be positively coexpressed more than expected by chance. When the similarity of transcription regulation is controlled for, adjacent genes have much higher coexpression rates than unlinked genes, supporting a role for chromatin modelling. Consequently, many incidences of low-to-moderate level coexpression of linked genes might well be spurious rather than an indication of functional coordination. These results have implications for gene therapy and for understanding gene order evolution, suggesting that chromosomal proximity alone is adequate to achieve some level of coexpression.

Cell- and gene-specific regulation of primary target genes by the androgen receptor

The androgen receptor (AR) mediates the physiologic and pathophysiologic effects of androgens including sexual differentiation, prostate development, and cancer progression by binding to genomic androgen response elements (AREs), which influence transcription of AR target genes. The composition and context of AREs differ between genes, thus enabling AR to confer multiple regulatory functions within a single nucleus. We used expression profiling of an immortalized human prostate epithelial cell line to identify 205 androgen-responsive genes (ARGs), most of them novel. In addition, we performed chromatin immunoprecipitation to identify 524 AR binding regions and validated in reporter assays the ARE activities of several such regions. Interestingly, 67% of our AREs resided within approximately 50 kb of the transcription start sites of 84% of our ARGs. Indeed, most ARGs were associated with two or more AREs, and ARGs were sometimes themselves linked in gene clusters containing up to 13 AREs and 12 ARGs. AREs appeared typically to be composite elements, containing AR binding sequences adjacent to binding motifs for other transcriptional regulators. Functionally, ARGs were commonly involved in prostate cell proliferation, communication, differentiation, and possibly cancer progression. Our results provide new insights into cell- and gene-specific mechanisms of transcriptional regulation of androgen-responsive gene networks.

Conservation of the links between gene transcription and chromosomal organization in the highly reduced genome of Buchnera aphidicola

BACKGROUND

Genomic studies on bacteria have clearly shown the existence of chromosomal organization as regards, for example, to gene localization, order and orientation. Moreover, transcriptomic analyses have demonstrated that, in free-living bacteria, gene transcription levels and chromosomal organization are mutually influenced. We have explored the possible conservation of relationships between mRNA abundances and chromosomal organization in the highly reduced genome of Buchnera aphidicola, the primary endosymbiont of the aphids, and a close relative to Escherichia coli.

RESULTS

Using an oligonucleotide-based microarray, we normalized the transcriptomic data by genomic DNA signals in order to have access to inter-gene comparison data. Our analysis showed that mRNA abundances, gene organization (operon) and gene essentiality are correlated in Buchnera (i.e., the most expressed genes are essential genes organized in operons) whereas no link between mRNA abundances and gene strand bias was found. The effect of Buchnera genome evolution on gene expression levels has also been analysed in order to assess the constraints imposed by the obligate symbiosis with aphids, underlining the importance of some gene sets for the survival of the two partners. Finally, our results show the existence of spatial periodic transcriptional patterns in the genome of Buchnera.

CONCLUSION

Despite an important reduction in its genome size and an apparent decay of its capacity for regulating transcription, this work reveals a significant correlation between mRNA abundances and chromosomal organization of the aphid-symbiont Buchnera.

Distinctive topologies of partner-switching signaling networks correlate with their physiological roles

Regulatory networks controlling bacterial gene expression often evolve from common origins and share homologous proteins and similar network motifs. However, when functioning in different physiological contexts, these motifs may be re-arranged with different topologies that significantly affect network performance. Here we analyze two related signaling networks in the bacterium Bacillus subtilis in order to assess the consequences of their different topologies, with the aim of formulating design principles applicable to other systems. These two networks control the activities of the general stress response factor sigma(B) and the first sporulation-specific factor sigma(F). Both networks have at their core a "partner-switching" mechanism, in which an anti-sigma factor forms alternate complexes either with the sigma factor, holding it inactive, or with an anti-anti-sigma factor, thereby freeing sigma. However, clear differences in network structure are apparent: the anti-sigma factor for sigma(F) forms a long-lived, "dead-end" complex with its anti-anti-sigma factor and ADP, whereas the genes encoding sigma(B) and its network partners lie in a sigma(B)-controlled operon, resulting in positive and negative feedback loops. We constructed mathematical models of both networks and examined which features were critical for the performance of each design. The sigma(F) model predicts that the self-enhancing formation of the dead-end complex transforms the network into a largely irreversible hysteretic switch; the simulations reported here also demonstrate that hysteresis and slow turn off kinetics are the only two system properties associated with this complex formation. By contrast, the sigma(B) model predicts that the positive and negative feedback loops produce graded, reversible behavior with high regulatory capacity and fast response time. Our models demonstrate how alterations in network design result in different system properties that correlate with regulatory demands. These design principles agree with the known or suspected roles of similar networks in diverse bacteria.

The SsrA-SmpB ribosome rescue system is important for growth of Bacillus subtilis at low and high temperatures

Bacillus subtilis has multiple stress response systems whose integrated action promotes growth and survival under unfavorable conditions. Here we address the function and transcriptional organization of a five-gene cluster containing ssrA, previously known to be important for growth at high temperature because of the role of its tmRNA product in rescuing stalled ribosomes. Reverse transcription-PCR experiments detected a single message for the secG-yvaK-rnr-smpB-ssrA cluster, suggesting that it constitutes an operon. However, rapid amplification of cDNA ends-PCR and lacZ fusion experiments indicated that operon transcription is complex, with at least five promoters controlling different segments of the cluster. One sigma(A)-like promoter preceded secG (P(1)), and internal sigma(A)-like promoters were found in both the rnr-smpB (P(2)) and smpB-ssrA intervals (P(3) and P(HS)). Another internal promoter lay in the secG-yvaK intercistronic region, and this activity (P(B)) was dependent on the general stress factor sigma(B). Null mutations in the four genes downstream from P(B) were tested for their effects on growth. Loss of yvaK (carboxylesterase E) or rnr (RNase R) caused no obvious phenotype. By contrast, smpB was required for growth at high temperature (52 degrees C), as anticipated if its product (a small ribosomal binding protein) is essential for tmRNA (ssrA) function. Notably, smpB and ssrA were also required for growth at low temperature (16 degrees C), a phenotype not previously associated with tmRNA activity. These results extend the known high-temperature role of ssrA and indicate that the ribosome rescue system is important at both extremes of the B. subtilis temperature range.

Population genetics and linkage analysis of loci within the FCT region of Streptococcus pyogenes

The FCT regions of Streptococcus pyogenes strains encode a variety of cell wall-anchored surface proteins that bind the extracellular matrix of the human host and/or give rise to pilus-like appendages. Strong linkage is evident between transcription-regulatory loci positioned within the FCT and emm regions and the emm pattern genotype marker for preferred infection of the throat or skin. These findings provide a basis for the hypothesis that FCT region gene products contribute to tissue-specific infection. In an initial series of steps to address this possibility, the FCT regions of 13 strains underwent comparative sequence analysis, the gene content of the FCT region was characterized for 113 strains via PCR, and genetic linkage was assessed. A history of extensive recombination within FCT regions was evident. The emm pattern D-defined skin specialist strains were highly homogenous in their FCT region gene contents, whereas the emm pattern A-C-defined throat specialist strains displayed a greater variety of forms. Most pattern A-C strains harbored prtF1 (75%) but lacked cpa (75%). In contrast, the majority of emm pattern D strains had cpa (92%) but lacked prtF1 (79%). Models based on FCT and emm region genotypes revealed the most parsimonious pathways of evolution. Using niche-determining candidate genes to infer phylogeny, emm pattern E strains--the so-called generalists, which lack a strong tissue site preference--occupied a transition zone separating most throat specialists from skin specialists. Overall, population genetic analysis supports the possibility that the FCT region gives rise to surface proteins that are largely necessary, but not always sufficient, to confer tissue site preference for infection.

Chlamydial type III secretion system is encoded on ten operons preceded by sigma 70-like promoter elements

Many gram-negative bacterial pathogens employ type III secretion systems for infectious processes. Chlamydiae are obligate intracellular bacteria that encode a conserved type III secretion system that is likely requisite for growth. Typically, genes encoding type III secretion systems are located in a single locus; however, for chlamydiae these genes are scattered throughout the genome. Little is known regarding the gene regulatory mechanisms for this essential virulence determinant. To facilitate identification of cis-acting transcriptional regulatory elements, the operon structure was determined. This analysis revealed 10 operons that contained 37 genes associated with the type III secretion system. Linkage within these operons suggests a role in type III secretion for each of these genes, including 13 genes encoding proteins with unknown function. The transcriptional start site for each operon was determined. In conjunction with promoter activity assays, this analysis revealed that the type III secretion system operons encode sigma(70)-like promoter elements. Transcriptional initiation by a sigma factor responsible for constitutive gene expression indicates that undefined activators or repressors regulate developmental stage-specific expression of chlamydial type III secretion system genes.

Whole genome sequencing of a novel temperate bacteriophage ofP. aeruginosa: evidence of tRNA gene mediating integration of the phage genome into the host bacterial chromosome

Whole genome sequencing of a novel Pseudomonas aeruginosa temperate bacteriophage PaP3 has been completed. The genome contains 45 503 bp with GC content of 52.1%, without more than 100 bp sequence hitting homologue in all sequenced phage genomes. A total of 256 open reading frames (ORFs) are found in the genome, and 71 ORFs are predicated as coding sequence (CDS). All 71 CDS are divided into the two opposite direction groups, and both groups meet at the bidirectional terminator site locating the near middle of the genome. The genome is dsDNA with 5'-protruded cohesive ends and cohesive sequence is 'GCCGGCCCCTTTCCGCGTTA' (20 mer). There are four tRNA genes (tRNA(Asn), tRNA(Asp), tRNA(Tyr) and tRNA(Pro)) clustering at the 5'-terminal of the genome. Analysis of integration site of PaP3 in the host bacterial genome confirmed that the core sequence of (GGTCGTAGGTTCGAATCCTAC-21mer) locates at tRNA(Pro) gene within the attP region and at tRNA(Lys) gene in the attB region. The results indicated that 3'-end of tRNA(Pro) gene of the PaP3 genome is involved in the integration reaction and 5'-end of tRNA(Lys) gene of host bacteria genome is hot spot of the integration.

An integrative method for accurate comparative genome mapping

We present MAGIC, an integrative and accurate method for comparative genome mapping. Our method consists of two phases: preprocessing for identifying "maximal similar segments," and mapping for clustering and classifying these segments. MAGIC's main novelty lies in its biologically intuitive clustering approach, which aims towards both calculating reorder-free segments and identifying orthologous segments. In the process, MAGIC efficiently handles ambiguities resulting from duplications that occurred before the speciation of the considered organisms from their most recent common ancestor. We demonstrate both MAGIC's robustness and scalability: the former is asserted with respect to its initial input and with respect to its parameters' values. The latter is asserted by applying MAGIC to distantly related organisms and to large genomes. We compare MAGIC to other comparative mapping methods and provide detailed analysis of the differences between them. Our improvements allow a comprehensive study of the diversity of genetic repertoires resulting from large-scale mutations, such as indels and duplications, including explicitly transposable and phagic elements. The strength of our method is demonstrated by detailed statistics computed for each type of these large-scale mutations. MAGIC enabled us to conduct a comprehensive analysis of the different forces shaping prokaryotic genomes from different clades, and to quantify the importance of novel gene content introduced by horizontal gene transfer relative to gene duplication in bacterial genome evolution. We use these results to investigate the breakpoint distribution in several prokaryotic genomes.

Is optimal gene order impossible?

Recent evidence suggests that yeast genes encoding proteins that are present in the same protein complex tend to be linked and to be co-expressed. More generally, we found that genes that are close to each other in the protein interaction network tend to be linked more often than expected and are often co-expressed. Unexpectedly, we found that linked genes in network proximity have unusually high recombination rates. Because high recombination rates are associated with high rates of genome re-organization, our findings might explain why the clustering of genes in proximity in the network is such a weak effect: there could be a co-evolutionary cycle of physical linkage for co-expression, upwards modification of the recombination rate and concomitant break-up of a cluster. Under such a model an "optimal" gene order is never stable.

Identification of the syr-syp box in the promoter regions of genes dedicated to syringomycin and syringopeptin production by Pseudomonas syringae pv. syringae B301D

The phytotoxins syringopeptin and syringomycin are synthesized by nonribosomal peptide synthetases which are encoded by the syringomycin (syr) and syringopeptin (syp) genomic island of Pseudomonas syringae pv. syringae. Previous studies demonstrated that expression of the syr-syp genes was controlled by the salA-syrF regulatory pathway, which in turn was induced by plant signal molecules. In this study, the 132-kb syr-syp genomic island was found to be organized into five polycistronic operons along with eight individual genes based on reverse transcriptional PCR and bioinformatic analysis. The transcriptional start sites of the salA gene and operons III and IV were located 63, 75, and 104 bp upstream of the start codons of salA, syrP, and syrB1, respectively, using primer extension analysis. The predicted -10/-35 promoter region of operon IV was confirmed based on deletion and site-directed mutagenesis analyses of the syrB1::uidA reporter with beta-glucuronidase assays. A 20-bp conserved sequence (TGtCccgN(6)cggGaCA, termed the syr-syp box) with dyad symmetry around the -35 region was identified via computer analysis for the syr-syp genes/operons responsible for biosynthesis and secretion of syringomycin and syringopeptin. Expression of the syrB1::uidA fusion was decreased 59% when 6 bp was deleted from the 5' end of the syr-syp box in the promoter region of operon IV. These results demonstrate that the conserved promoter sequences of the syr-syp genes contribute to the coregulation of syringomycin and syringopeptin production.

Clusters of adjacent and similarly expressed genes across normal human tissues complicate comparative transcriptomic discovery

Transcriptomic techniques are valuable tools with which to validate genetic and biological hypotheses and are now widely available for research. However, with the exception of tumor biology, comparative genomics analyses have been difficult to use as discovery engines to describe biologically relevant expression changes. We propose that physical proximity of human genes correlates with similar mRNA expression, so that increased expression might include a disease-relevant gene and many other genes in the adjacent region. To increase the efficiency of combining susceptibility gene mapping and interpretation of transcriptomics, we developed a method to identify clusters of adjacent and similarly expressed genes. Gene expression profiles for 28,945 genes across 101 normal human tissues were obtained from the Gene Logic BioExpress system. The expression similarity for genes in sliding-windows was measured using average pair-wise Pearson correlation coefficients. We identified 187 clusters (p < 10e-4) of co-regulated genes, including 2648 genes, or 9.1% of all genes considered and termed these "clusters of adjacent and similarly expressed genes" (CASEGs). Genes in 15 (8.2%) of these clusters demonstrate a significant co-expression enrichment (p < 10e-10). This study demonstrates the coordinate expression of neighboring genes and provides a comprehensive view of expression-based compartmentalization of the human genome, which can be overlaid on genetic susceptibility gene maps.