Optimal gene partition into operons correlates with gene functional order

Gene arrangement into operons varies between bacterial species. Genes in a given system can be on one operon in some organisms and on several operons in other organisms. Existing theories explain why genes that work together should be on the same operon, since this allows for advantageous lateral gene transfer and accurate stoichiometry. But what causes the frequent separation into multiple operons of co-regulated genes that act together in a pathway? Here we suggest that separation is due to benefits made possible by differential regulation of each operon. We present a simple mathematical model for the optimal distribution of genes into operons based on a balance of the cost of operons and the benefit of regulation that provides 'just-when-needed' temporal order. The analysis predicts that genes are arranged such that genes on the same operon do not skip functional steps in the pathway. This prediction is supported by genomic data from 137 bacterial genomes. Our work suggests that gene arrangement is not only the result of random historical drift, genome re-arrangement and gene transfer, but has elements that are solutions of an evolutionary optimization problem. Thus gene functional order may be inferred by analyzing the operon structure across different genomes.     

Regulation of bacterial respiration: Comparison of microarray and comparative genomics data

A comparison was made of the structures of the Fnr and ArcA modulons and regulons. The data on modulon composition were taken from published microarray assays, whereas regulons were characterized using comparative genomic approaches. The regulatory cascade involving Fnr and ArcA contributes greatly to the extension of the Fnr modulon over the Fnr regulon by adding operons of the ArcA modulon. The Fnr and ArcA regulons were shown to contain 26 and 16 operons, respectively. Ten operons had high-score and highly conserved sites for both Fnr and ArcA and were isolated as a so-called core of regulons.     

Prediction and experimental testing of ferric uptake regulator regulons in vibrios

Iron homeostasis is in many bacteria regulated by the ferric uptake regulator (Fur). Despite the available information on Fur regulons, it is likely that there are Fur-regulated genes and operons that are unique to vibrios, and knowledge into these can potentially provide new insights into vibrio virulence and pathogenesis. We constructed a vibrio-specific alignment matrix based on Fur-binding sites from the literature and used existing software (Patser) to search five published vibrio genomes and the Vibrio salmonicida draft genome for Fur-regulated genes. The consensus Fur-binding site from our matrix is 5'-AATGANAATNATTNTCATT-3'. Fur-binding motifs were found associated with 50-61 single genes and 16-20 operons in each genome. Predictions were tested by monitoring the expression of a subset of genes and operons in V. salmonicida. Six previously undescribed Fur-regulated genes showed increased expression under iron-restrictive conditions. Our work provides a comprehensive list of predicted Fur regulons in six vibrio genomes, which may be used to generate new hypotheses for future experiments.     

Nebulon: a system for the inference of functional relationships of gene products from the rearrangement of predicted operons

Since operons are unstable across Prokaryotes, it has been suggested that perhaps they re-combine in a conservative manner. Thus, genes belonging to a given operon in one genome might re-associate in other genomes revealing functional relationships among gene products. We developed a system to build networks of functional relationships of gene products based on their organization into operons in any available genome. The operon predictions are based on inter-genic distances. Our system can use different kinds of thresholds to accept a functional relationship, either related to the prediction of operons, or to the number of non-redundant genomes that support the associations. We also work by shells, meaning that we decide on the number of linking iterations to allow for the complementation of related gene sets. The method shows high reliability benchmarked against knowledge-bases of functional interactions. We also illustrate the use of Nebulon in finding new members of regulons, and of other functional groups of genes. Operon rearrangements produce thousands of high-quality new interactions per prokaryotic genome, and thousands of confirmations per genome to other predictions, making it another important tool for the inference of functional interactions from genomic context.     

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.     

Three-dimensional genome organization in interphase and its relation to genome function

Higher order chromatin structure, i.e. the three-dimensional (3D) organization of the genome in the interphase nucleus, is an important component in the orchestration of gene expression in the mammalian genome. In this review we describe principles of higher order chromatin structure discussing three organizational parameters, i.e. chromatin folding, chromatin compaction and the nuclear position of the chromatin fibre. We argue that principles of 3D genome organization are probabilistic traits, reflected in a considerable cell-to-cell variation in 3D genome structure. It will be essential to understand how such higher order organizational aspects contribute to genome function to unveil global genome regulation.     

A unified framework for inferring the multi-scale organization of chromatin domains from Hi-C

Chromosomes are giant chain molecules organized into an ensemble of three-dimensional structures characterized with its genomic state and the corresponding biological functions. Despite the strong cell-to-cell heterogeneity, the cell-type specific pattern demonstrated in high-throughput chromosome conformation capture (Hi-C) data hints at a valuable link between structure and function, which makes inference of chromatin domains (CDs) from the pattern of Hi-C a central problem in genome research. Here we present a unified method for analyzing Hi-C data to determine spatial organization of CDs over multiple genomic scales. By applying statistical physics-based clustering analysis to a polymer physics model of the chromosome, our method identifies the CDs that best represent the global pattern of correlation manifested in Hi-C. The multi-scale intra-chromosomal structures compared across different cell types uncover the principles underlying the multi-scale organization of chromatin chain: (i) Sub-TADs, TADs, and meta-TADs constitute a robust hierarchical structure. (ii) The assemblies of compartments and TAD-based domains are governed by different organizational principles. (iii) Sub-TADs are the common building blocks of chromosome architecture. Our physically principled interpretation and analysis of Hi-C not only offer an accurate and quantitative view of multi-scale chromatin organization but also help decipher its connections with genome function.     

The evolutionary dynamics of functional modules and the extraordinary plasticity of regulons: the Escherichia coli perspective

Using profiles of phylogenetic profiles (P-cubic) we compared the evolutionary dynamics of different kinds of functional associations. Ordered from most to least evolutionarily stable, these associations were genes in the same operons, genes whose products participate in the same biochemical pathway, genes coding for physically interacting proteins and genes in the same regulons. Regulons showed the most plastic functional interactions with evolutionary stabilities barely better than those of unrelated genes. Further regulon analyses showed that global regulators contain less evolutionarily stable associations than local regulators. Genes co-repressed by global regulators had a higher evolutionary conservation than genes co-activated by global regulators. However, the reverse was true for genes co-repressed and co-activated by local regulators. Of all the regulon-related associations, the relationship between regulators and their target genes showed the most evolutionary stability. Different negative data sets built to contrast against each of the analysed kinds of modules also differed in evolutionary conservation revealing further underlying genome organization. Applying P-cubic analyses to other genomes might help visualize genome organization, understand the evolutionary importance and plasticity of functional associations and compare the quality of data sets expected to reflect functional interactions, such as those coming from high-throughput experiments.     

Conservation of the mosaic structure of the four internal transcribed spacers and localisation of the rrn operons on the Streptococcus pneumoniae genome

The detection of heterogeneity of the 16S-23S ribosomal intergenic transcribed spacer (ITS) region has become rather common over the past years for identification and typing purposes of bacteria. The ITS not only varies in sequence and length, but also in number of alleles per genome and in their position on the chromosome together with the ribosomal clusters. The ITS characterisation has allowed discrimination of several species within a genus and variation in ITS sequences between the multiple rrn operons present within a genome may be as high or greater than between strains of the same species or subspecies. It is important to understand the variability of ITS sequences in a given genome to gain insights into bacterial physiology and taxonomy. The present study describes the possibility to type Streptococcus pneumoniae by PCR-ribotyping of the spacer region, the determination of the molecular structure of the ITS, and the determination of the number and localisation of rrn operons in this microorganism. Our results show that the genome of S. pneumoniae contains four ribosomal operons, showing the same genomic organisation among strains, each containing a single ITS allele of 270 bp. The ITS sequence presents a mosaic organisation of blocks highly conserved intra- and inter-species within the genus Streptococcus, giving no possibility for variations to arise.