Periodic Transcriptional Organization of the E.coli Genome

Transcriptional regulation shapes the organization of genes on bacterial chromosomes

Transcription factors (TFs) are the key elements responsible for controlling the expression of genes in bacterial genomes and when visualized on a genomic scale form a dense network of transcriptional interactions among themselves and with other protein coding genes. Although the structure of transcriptional regulatory networks (TRNs) is well understood, it is not clear what constrains govern them. Here, we explore this question using the TRNs of model prokaryotes and provide a link between the transcriptional hierarchy of regulons and their genome organization. We show that, to drive the kinetics and concentration gradients, TFs belonging to big and small regulons, depending on the number of genes they regulate, organize themselves differently on the genome with respect to their targets. We then propose a conceptual model that can explain how the hierarchical structure of TRNs might be ultimately governed by the dynamic biophysical requirements for targeting DNA-binding sites by TFs. Our results suggest that the main parameters defining the position of a TF in the network hierarchy are the number and chromosomal distances of the genes they regulate and their protein concentration gradients. These observations give insights into how the hierarchical structure of transcriptional networks can be encoded on the chromosome to drive the kinetics and concentration gradients of TFs depending on the number of genes they regulate and could be a common theme valid for other prokaryotes, proposing the role of transcriptional regulation in shaping the organization of genes on a chromosome.

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.

Differential replication dynamics for large and small Vibrio chromosomes affect gene dosage, expression and location

Background

Replication of bacterial chromosomes increases copy numbers of genes located near origins of replication relative to genes located near termini. Such differential gene dosage depends on replication rate, doubling time and chromosome size. Although little explored, differential gene dosage may influence both gene expression and location. For vibrios, a diverse family of fast growing gammaproteobacteria, gene dosage may be particularly important as they harbor two chromosomes of different size.

Results

Here we examined replication dynamics and gene dosage effects for the separate chromosomes of three Vibrio species. We also investigated locations for specific gene types within the genome. The results showed consistently larger gene dosage differences for the large chromosome which also initiated replication long before the small. Accordingly, large chromosome gene expression levels were generally higher and showed an influence from gene dosage. This was reflected by a higher abundance of growth essential and growth contributing genes of which many locate near the origin of replication. In contrast, small chromosome gene expression levels were low and appeared independent of gene dosage. Also, species specific genes are highly abundant and an over-representation of genes involved in transcription could explain its gene dosage independent expression.

Conclusion

Here we establish a link between replication dynamics and differential gene dosage on one hand and gene expression levels and the location of specific gene types on the other. For vibrios, this relationship appears connected to a polarisation of genetic content between its chromosomes, which may both contribute to and be enhanced by an improved adaptive capacity.

Dissecting the logical types of network control in gene expression profiles

Background

In the bacterium Escherichia coli the transcriptional regulation of gene expression involves both dedicated regulators binding specific DNA sites with high affinity and also global regulators – abundant DNA architectural proteins of the bacterial nucleoid binding multiple sites with a wide range of affinities and thus modulating the superhelical density of DNA. The first form of transcriptional regulation is predominantly pairwise and specific, representing digitial control, while the second form is (in strength and distribution) continuous, representing analog control.

Results

Here we look at the properties of effective networks derived from significant gene expression changes under variation of the two forms of control and find that upon limitations of one type of control (caused e.g. by mutation of a global DNA architectural factor) the other type can compensate for compromised regulation. Mutations of global regulators significantly enhance the digital control, whereas in the presence of global DNA architectural proteins regulation is mostly of the analog type, coupling spatially neighboring genomic loci. Taken together our data suggest that two logically distinct – digital and analog – types of control are balancing each other.

Conclusion

By revealing two distinct logical types of control, our approach provides basic insights into both the organizational principles of transcriptional regulation and the mechanisms buffering genetic flexibility. We anticipate that the general concept of distinguishing logical types of control will apply to many complex biological networks.

How gene order is influenced by the biophysics of transcription regulation

What are the forces that shape the structure of prokaryotic genomes: the order of genes, their proximity, and their orientation? Coregulation and coordinated horizontal gene transfer are believed to promote the proximity of functionally related genes and the formation of operons. However, forces that influence the structure of the genome beyond the level of a single operon remain unknown. Here, we show that the biophysical mechanism by which regulatory proteins search for their sites on DNA can impose constraints on genome structure. Using simulations, we demonstrate that rapid and reliable gene regulation requires that the transcription factor (TF) gene be close to the site on DNA the TF has to bind, thus promoting the colocalization of TF genes and their targets on the genome. We use parameters that have been measured in recent experiments to estimate the relevant length and times scales of this process and demonstrate that the search for a cognate site may be prohibitively slow if a TF has a low copy number and is not colocalized. We also analyze TFs and their sites in a number of bacterial genomes, confirm that they are colocalized significantly more often than expected, and show that this observation cannot be attributed to the pressure for coregulation or formation of selfish gene clusters, thus supporting the role of the biophysical constraint in shaping the structure of prokaryotic genomes. Our results demonstrate how spatial organization can influence timing and noise in gene expression.

The correlation between architecture and mRNA abundance in the genetic regulatory network of Escherichia coli

Background

Two aspects of genetic regulatory networks are the static architecture that describes the overall connectivity between the genes and the dynamics that describes the sequence of genes active at any one time as deduced from mRNA abundances. The nature of the relationship between these two aspects of these networks is a fundamental question. To address it, we have used the static architecture of the connectivity of the regulatory proteins of Escherichia coli to analyse their relationship to the abundance of the mRNAs encoding these proteins. In this we build on previous work which uses Boolean network models, but impose biological constraints that cannot be deduced from the mRNA abundances alone.

Results

For a cell population of E. coli, we find that there is a strong and statistically significant linear dependence between the abundance of mRNA encoding a regulatory protein and the number of genes regulated by this protein. We use this result, together with the ratio of regulatory repressors to promoters, to simulate numerically a genetic regulatory network of a single cell. The resulting model exhibits similar correlations to that of E. coli.

Conclusion

This analysis clarifies the relationship between the static architecture of a regulatory network and the consequences for the dynamics of its pattern of mRNA abundances. It also provides the constraints on the architecture required to construct a model network to simulate mRNA production.

Chromosomal periodicity of evolutionarily conserved gene pairs

Chromosomes are compacted hundreds of times to fit in the cell, packaged into dynamic folds whose structures are largely unknown. Here, we examine patterns in gene locations to infer large-scale features of bacterial chromosomes. Specifically, we analyzed >100 genomes and identified thousands of gene pairs that display two types of evolutionary correlations: a tendency to co-occur and a tendency to be located close together in many genomes. We then analyzed the detailed distribution of these pairs in Escherichia coli and found that genes in a pair tend to be separated by integral multiples of 117 kb along the genome and to be positioned in a 117-kb grid of genomic locations. In addition, the most pair-dense locations coincide with regions of intense transcriptional activity and the positions of top transcribed and conserved genes. These patterns suggest that the E. coli chromosome may be organized into a 117-kb helix-like topology that localizes a subset of the most essential and highly transcribed genes along a specific face of this structure. Our approach indicates an evolutionarily maintained preference in the spacing of genes along the chromosome and offers a general comparative genomics framework for studying chromosome structure, broadly applicable to other organisms.

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.

Locational distribution of gene functional classes in Arabidopsis thaliana

Background

We are interested in understanding the locational distribution of genes and their functions in genomes, as this distribution has both functional and evolutionary significance. Gene locational distribution is known to be affected by various evolutionary processes, with tandem duplication thought to be the main process producing clustering of homologous sequences. Recent research has found clustering of protein structural families in the human genome, even when genes identified as tandem duplicates have been removed from the data. However, this previous research was hindered as they were unable to analyse small sample sizes. This is a challenge for bioinformatics as more specific functional classes have fewer examples and conventional statistical analyses of these small data sets often produces unsatisfactory results.

Results

We have developed a novel bioinformatics method based on Monte Carlo methods and Greenwood's spacing statistic for the computational analysis of the distribution of individual functional classes of genes (from GO). We used this to make the first comprehensive statistical analysis of the relationship between gene functional class and location on a genome. Analysis of the distribution of all genes except tandem duplicates on the five chromosomes of A. thaliana reveals that the distribution on chromosomes I, II, IV and V is clustered at P = 0.001. Many functional classes are clustered, with the degree of clustering within an individual class generally consistent across all five chromosomes. A novel and surprising result was that the locational distribution of some functional classes were significantly more evenly spaced than would be expected by chance.

Conclusion

Analysis of the A. thaliana genome reveals evidence of unexplained order in the locational distribution of genes. The same general analysis method can be applied to any genome, and indeed any sequential data involving classes.

In vivo DNase I sensitivity of the Streptomyces coelicolor chromosome correlates with gene expression: implications for bacterial chromosome structure

For a bacterium, Streptomyces coelicolor A3(2) contains a relatively large genome (8.7 Mb) with a complex and adaptive pattern of gene regulation. We discovered a correlation between the physical structure of the S.coelicolor genome and the transcriptional activity of the genes therein. Twelve genes were surveyed throughout 72 h of growth for both in vivo sensitivity to DNase I digestion and levels of transcription. DNase I-sensitivity correlated positively with transcript levels, implying that it was predictive of gene expression, and indicating increased accessibility of transcribed DNA. The genome was fractionated based on the sensitivity to DNase I digestion, with the low molecular weight (frequently cut) fraction highly enriched for actively transcribed sequences when compared to the infrequently cut fraction, which was representative of the entire genome. This approach will allow comparison of nucleoid proteins, and any modifications thereof, associated with transcriptionally active and inactive regions of the bacterial genome.

Guidelines: From artificial evolution to computational evolution: a research agenda

Computational scientists have developed algorithms inspired by natural evolution for at least 50 years. These algorithms solve optimization and design problems by building solutions that are 'more fit' relative to desired properties. However, the basic assumptions of this approach are outdated. We propose a research programme to develop a new field: computational evolution. This approach will produce algorithms that are based on current understanding of molecular and evolutionary biology and could solve previously unimaginable or intractable computational and biological problems.

Bacterial chromatin

Recent studies have revealed that the bacterial nucleoid is a dynamic entity that alters its overall structure in response to changes in both growth rate and growth phase. These structural changes are correlated with, and might be driven by, changes in the distribution and utilization of DNA supercoiling. In turn, these parameters in addition to the delimitation of topological domains are dependent both on the relative proportions of the abundant nucleoid-associated proteins and on transcriptional activity. The domain structure itself is dynamic.

Regional organization of gene expression in Streptomyces coelicolor

Based on the chromosomal locations of genes inferred from sequence analysis to be essential for the viability of Streptomyces coelicolor, Bentley et al. [Bentley, S.D., et al. 2002. Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2), Nature 417, 141-147.] have suggested that a 4.9 Mb central region of the linear S. coelicolor chromosome encodes 'core' functions expressed during vegetative growth of this species, while 1.5 Mb and 2.3 Mb chromosomal DNA segments lateral to this core encode auxiliary functions proposed to be required under other growth conditions. To examine this hypothesis and experimentally identify genes expressed during vegetative growth of S. coelicolor cultures, we used DNA microarrays to measure globally the abundance of S. coelicolor transcripts in cells growing in liquid medium. We found that, overall, genes corresponding to the 4.9 Mb core region of the S. coelicolor M145 chromosome were more highly expressed under non-limiting growth conditions than genes in the 1.5 Mb left and 2.3 Mb right chromosome arms, supporting the notion of the core versus auxiliary organization of genes on the chromosome. To examine how this chromosomal distribution of transcripts changes under other growth conditions, we also measured gene expression changes during stationary phase and several stress conditions. During stationary phase, the composition of S. coelicolor transcripts appears to shift from large quantities of growth-related transcripts encoded in the core region to those of less characterized genes, which may be essential for differentiation and other physiological responses, encoded throughout the chromosome. After temperature and osmotic upshifts, we found that S. coelicolor transiently induces a set of several hundred genes located throughout the chromosome, which may function in response mechanisms common to the two stress conditions.

Decoding the nucleoid organisation of Bacillus subtilis and Escherichia coli through gene expression data

Background

Although the organisation of the bacterial chromosome is an area of active research, little is known yet on that subject. The difficulty lies in the fact that the system is dynamic and difficult to observe directly. The advent of massive hybridisation techniques opens the way to further studies of the chromosomal structure because the genes that are co-expressed, as identified by microarray experiments, probably share some spatial relationship. The use of several independent sets of gene expression data should make it possible to obtain an exhaustive view of the genes co-expression and thus a more accurate image of the structure of the chromosome.

Results

For both Bacillus subtilis and Escherichia coli the co-expression of genes varies as a function of the distance between the genes along the chromosome. The long-range correlations are surprising: the changes in the level of expression of any gene are correlated (positively or negatively) to the changes in the expression level of other genes located at well-defined long-range distances. This property is true for all the genes, regardless of their localisation on the chromosome.

We also found short-range correlations, which suggest that the location of these co-expressed genes corresponds to DNA turns on the nucleoid surface (14–16 genes).

Conclusion

The long-range correlations do not correspond to the domains so far identified in the nucleoid. We explain our results by a model of the nucleoid solenoid structure based on two types of spirals (short and long). The long spirals are uncoiled expressed DNA while the short ones correspond to coiled unexpressed DNA.

L’épigénétique comme aspect de la postgénomique

A cell transmits to its progeny the activity level of many of its genes, not just their sequence. Just like the sequence may vary through a mutation, the gene activity level may change through an "epimutation" (an epigenetic modification) which is heritable and does not entail any concomitant genetic alteration. An epimutation can have important phenotypic consequences, that eventually survive to the loss of the environmental conditions that triggered it. For instance, epimutations are responsible for the divergence between a neuron and an epithelial cell that both come from the same egg and contain the same genome complement. This phenotypic difference is much larger than the one between the neurons from two animal species with dissimilar genotypes, thereby underlining the importance of epimutations. Tradition opposes the genetic and epigenetic visions, the latter being often adequated to the DNA methylation phenomenon. However, epimutations display a rich spectrum of modes that can all fit in a unique reference system based on correlated chemical, spatial and temporal scales. This reference system allows the integration of purely genetic mutations at one of its ends, thus paving the way to a new, gradual vision that encompasses the genome and the epigenome. At the other end can be found two types of epimutations that are both wide-ranging in space and rapid in producing phenotypic alterations. Firstly, long-range rearrangements of the three-dimensional structure of the chromosome may influence gene expression in an heritable fashion. Such rearrangements seem to result from the collective dynamics of DNA-related activities, particularly transcription. Lastly, heritable regulatory states, e.g. a differentiated state that results from tipping a regulatory "toggle switch", involve components that are distributed throughout the nucleus or the cytoplasm, and possibly all the way to cell confines.

DNA supercoiling — a global transcriptional regulator for enterobacterial growth?

A fundamental principle of exponential bacterial growth is that no more ribosomes are produced than are necessary to support the balance between nutrient availability and protein synthesis. Although this conclusion was first expressed more than 40 years ago, a full understanding of the molecular mechanisms involved remains elusive and the issue is still controversial. There is currently agreement that, although many different systems are undoubtedly involved in fine-tuning this balance, an important control, and in our opinion perhaps the main control, is regulation of the rate of transcription initiation of the stable (ribosomal and transfer) RNA transcriptons. In this review, we argue that regulation of DNA supercoiling provides a coherent explanation for the main modes of transcriptional control - stringent control, growth-rate control and growth-phase control - during the normal growth of Escherichia coli.