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

Bacterial nucleoid is a riddle wrapped in a mystery inside an enigma

Bacterial chromosome, the nucleoid, is traditionally modeled as a rosette of DNA mega-loops, organized around proteinaceous central scaffold by nucleoid-associated proteins (NAPs), and mixed with the cytoplasm by transcription and translation. Electron microscopy of fixed cells confirms dispersal of the cloud-like nucleoid within the ribosome-filled cytoplasm. Here, I discuss evidence that the nucleoid in live cells forms DNA phase separate from riboprotein phase, the "riboid." I argue that the nucleoid-riboid interphase, where DNA interacts with NAPs, transcribing RNA polymerases, nascent transcripts, and ssRNA chaperones, forms the transcription zone. An active part of phase separation, transcription zone enforces segregation of the centrally positioned information phase (the nucleoid) from the surrounding action phase (the riboid), where translation happens, protein accumulates, and metabolism occurs. I speculate that HU NAP mostly tiles up the nucleoid periphery-facilitating DNA mobility but also supporting transcription in the interphase. Besides extruding plectonemically supercoiled DNA mega-loops, condensins could compact them into solenoids of uniform rings, while HU could support rigidity and rotation of these DNA rings. The two-phase cytoplasm arrangement allows the bacterial cell to organize the central dogma activities, where (from the cell center to its periphery) DNA replicates and segregates, DNA is transcribed, nascent mRNA is handed over to ribosomes, mRNA is translated into proteins, and finally, the used mRNA is recycled into nucleotides at the inner membrane. The resulting information-action conveyor, with one activity naturally leading to the next one, explains the efficiency of prokaryotic cell design-even though its main intracellular transportation mode is free diffusion.

Examination of Large Chromosomal Inversions in the Genome of Erwinia amylovora Strains Reveals Worldwide Distribution and North America-Specific Types

Erwinia amylovora is a relatively homogeneous species with low genetic diversity at the nucleotide level. However, phenotypic differences and genomic structural variations among E. amylovora strains have been documented. In this study, we identified 10 large chromosomal inversion (LCI) types in the Spiraeoideae-infecting (SI) E. amylovora strains by combining whole genome sequencing and PCR-based molecular markers. It was found that LCIs were mainly caused by homologous recombination events among seven rRNA operons (rrns) in SI E. amylovora strains. Although ribotyping results identified inter- and intra-variations in the internal transcribed spacer (ITS1 and ITS2) regions among rrns, LCIs tend to occur between rrns transcribed in the opposite directions and with the same tRNA content (tRNA-Glu or tRNA-Ile/Ala) in ITS1. Based on the LCI types, physical/estimated replichore imbalance (PRI/ERI) was examined and calculated. Among the 117 SI strains evaluated, the LCI types of Ea1189, CFBP1430, and Ea273 were the most common, with ERI values at 1.31, 7.87, and 4.47°, respectively. These three LCI types had worldwide distribution, whereas the remaining seven LCI types were restricted to North America (or certain regions of the United States). Our results indicated ongoing chromosomal recombination events in the SI E. amylovora population and showed that LCI events are mostly symmetrical, keeping the ERI less than 15°. These findings provide initial evidence about the prevalence of certain LCI types in E. amylovora strains, how LCI occurs, and its potential evolutionary advantage and history, which might help track the movement of the pathogen.

Spatiotemporal Coupling of DNA Supercoiling and Genomic Sequence Organization-A Timing Chain for the Bacterial Growth Cycle?

In this article we describe the bacterial growth cycle as a closed, self-reproducing, or autopoietic circuit, reestablishing the physiological state of stationary cells initially inoculated in the growth medium. In batch culture, this process of self-reproduction is associated with the gradual decline in available metabolic energy and corresponding change in the physiological state of the population as a function of "travelled distance" along the autopoietic path. We argue that this directional alteration of cell physiology is both reflected in and supported by sequential gene expression along the chromosomal OriC-Ter axis. We propose that during the E. coli growth cycle, the spatiotemporal order of gene expression is established by coupling the temporal gradient of supercoiling energy to the spatial gradient of DNA thermodynamic stability along the chromosomal OriC-Ter axis.

Identification of the key functional genes in salt-stress tolerance of Cyanobacterium Phormidium tenue using in silico analysis

The development of artificial biocrust using cyanobacterium Phormidium tenue has been suggested as an effective strategy to prevent soil degradation. Here, a combination of in silico approaches with growth rate, photosynthetic pigment, morphology, and transcript analysis was used to identify specific genes and their protein products in response to 500 mM NaCl in P. tenue. The results show that 500 mM NaCl induces the expression of genes encoding glycerol-3-phosphate dehydrogenase (glpD) as a Flavoprotein, ribosomal protein S12 methylthiotransferase (rimO), and a hypothetical protein (sll0939). The constructed co-expression network revealed a group of abiotic stress-responsive genes. Using the Basic Local Alignment Search Tool (BLAST), the homologous proteins of rimO, glpD, and sll0939 were identified in the P. tenue genome. Encoded proteins of glpD, rimO, and DUF1622 genes, respectively, contain (DAO and DAO C), (UPF0004, Radical SAM and TRAM 2), and (DUF1622) domains. The predicted ligand included 22B and MG for DUF1622, FS5 for rimO, and FAD for glpD protein. There was no direct disruption in ligand-binding sites of these proteins by Na+, Cl-, or NaCl. The growth rate, photosynthetic pigment, and morphology of P. tenue were investigated, and the result showed an acceptable tolerance rate of this microorganism under salt stress. The quantitative real-time polymerase chain reaction (qRT-PCR) results revealed the up-regulation of glpD, rimO, and DUF1622 genes under salt stress. This is the first report on computational and experimental analyses of the glpD, rimO, and DUF1622 genes in P. tenue under salt stress to the best of our knowledge.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03050-w.

Universal and idiosyncratic characteristic lengths in bacterial genomes

In condensed matter physics, simplified descriptions are obtained by coarse-graining the features of a system at a certain characteristic length, defined as the typical length beyond which some properties are no longer correlated. From a physics standpoint, in vitro DNA has thus a characteristic length of 300 base pairs (bp), the Kuhn length of the molecule beyond which correlations in its orientations are typically lost. From a biology standpoint, in vivo DNA has a characteristic length of 1000 bp, the typical length of genes. Since bacteria live in very different physico-chemical conditions and since their genomes lack translational invariance, whether larger, universal characteristic lengths exist is a non-trivial question. Here, we examine this problem by leveraging the large number of fully sequenced genomes available in public databases. By analyzing GC content correlations and the evolutionary conservation of gene contexts (synteny) in hundreds of bacterial chromosomes, we conclude that a fundamental characteristic length around 10-20 kb can be defined. This characteristic length reflects elementary structures involved in the coordination of gene expression, which are present all along the genome of nearly all bacteria. Technically, reaching this conclusion required us to implement methods that are insensitive to the presence of large idiosyncratic genomic features, which may co-exist along these fundamental universal structures.

Chromosomal organization of transcription: in a nutshell

Early studies of transcriptional regulation focused on individual gene promoters defined specific transcription factors as central agents of genetic control. However, recent genome-wide data propelled a different view by linking spatially organized gene expression patterns to chromosomal dynamics. Therefore, the major problem in contemporary molecular genetics concerned with transcriptional gene regulation is to establish a unifying model that reconciles these two views. This problem, situated at the interface of polymer physics and network theory, requires development of an integrative methodology. In this review, we discuss recent achievements in classical model organism E. coli and provide some novel insights gained from studies of a bacterial plant pathogen, D. dadantii. We consider DNA topology and the basal transcription machinery as key actors of regulation, in which activation of functionally relevant genes is coupled to and coordinated with the establishment of extended chromosomal domains of coherent transcription. We argue that the spatial organization of genome plays a fundamental role in its own regulation.

Insights into the Mechanisms of Basal Coordination of Transcription Using a Genome-Reduced Bacterium

Coordination of transcription in bacteria occurs at supra-operonic scales, but the extent, specificity, and mechanisms of such regulation are poorly understood. Here, we tackle this problem by profiling the transcriptome of the model organism Mycoplasma pneumoniae across 115 growth conditions. We identify three qualitatively different levels of co-expression corresponding to distinct relative orientations and intergenic properties of adjacent genes. We reveal that the degree of co-expression between co-directional adjacent operons, and more generally between genes, is tightly related to their capacity to be transcribed en bloc into the same mRNA. We further show that this genome-wide pervasive transcription of adjacent genes and operons is specifically repressed by DNA regions preferentially bound by RNA polymerases, by intrinsic terminators, and by large intergenic distances. Taken together, our findings suggest that the basal coordination of transcription is mediated by the physical entities and mechanical properties of the transcription process itself, and that operon-like behaviors may strongly vary from condition to condition.

Conserved Units of Co-Expression in Bacterial Genomes: An Evolutionary Insight into Transcriptional Regulation

Genome-wide measurements of transcriptional activity in bacteria indicate that the transcription of successive genes is strongly correlated beyond the scale of operons. Here, we analyze hundreds of bacterial genomes to identify supra-operonic segments of genes that are proximal in a large number of genomes. We show that these synteny segments correspond to genomic units of strong transcriptional co-expression. Structurally, the segments contain operons with specific relative orientations (co-directional or divergent) and nucleoid-associated proteins are found to bind at their boundaries. Functionally, operons inside a same segment are highly co-expressed even in the apparent absence of regulatory factors at their promoter regions. Remote operons along DNA can also be co-expressed if their corresponding segments share a transcriptional or sigma factor, without requiring these factors to bind directly to the promoters of the operons. As evidence that these results apply across the bacterial kingdom, we demonstrate them both in the Gram-negative bacterium Escherichia coli and in the Gram-positive bacterium Bacillus subtilis. The underlying process that we propose involves only RNA-polymerases and DNA: it implies that the transcription of an operon mechanically enhances the transcription of adjacent operons. In support of a primary role of this regulation by facilitated co-transcription, we show that the transcription en bloc of successive operons as a result of transcriptional read-through is strongly and specifically enhanced in synteny segments. Finally, our analysis indicates that facilitated co-transcription may be evolutionary primitive and may apply beyond bacteria.

GREAT: a web portal for Genome Regulatory Architecture Tools

GREAT (Genome REgulatory Architecture Tools) is a novel web portal for tools designed to generate user-friendly and biologically useful analysis of genome architecture and regulation. The online tools of GREAT are freely accessible and compatible with essentially any operating system which runs a modern browser. GREAT is based on the analysis of genome layout -defined as the respective positioning of co-functional genes- and its relation with chromosome architecture and gene expression. GREAT tools allow users to systematically detect regular patterns along co-functional genomic features in an automatic way consisting of three individual steps and respective interactive visualizations. In addition to the complete analysis of regularities, GREAT tools enable the use of periodicity and position information for improving the prediction of transcription factor binding sites using a multi-view machine learning approach. The outcome of this integrative approach features a multivariate analysis of the interplay between the location of a gene and its regulatory sequence. GREAT results are plotted in web interactive graphs and are available for download either as individual plots, self-contained interactive pages or as machine readable tables for downstream analysis. The GREAT portal can be reached at the following URL https://absynth.issb.genopole.fr/GREAT and each individual GREAT tool is available for downloading.

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.

Spatial organization of transcription in bacterial cells

Prokaryotic transcription has been extensively studied over the past half a century. However, there often exists a gap between the structural, mechanistic description of transcription obtained from in vitro biochemical studies, and the cellular, phenomenological observations from in vivo genetic studies. It is now accepted that a living bacterial cell is a complex entity; the heterogeneous cellular environment is drastically different from the homogenous, well-mixed situation in vitro. Where molecules are inside a cell may be important for their function; hence, the spatial organization of different molecular components may provide a new means of transcription regulation in vivo, possibly bridging this gap. In this review, we survey current evidence for the spatial organization of four major components of transcription [genes, transcription factors, RNA polymerase (RNAP) and RNAs] and critically analyze their biological significance.

A model for Escherichia coli chromosome packaging supports transcription factor-induced DNA domain formation

What physical mechanism leads to organization of a highly condensed and confined circular chromosome? Computational modeling shows that confinement-induced organization is able to overcome the chromosome's propensity to mix by the formation of topological domains. The experimentally observed high precision of separate subcellular positioning of loci (located on different chromosomal domains) in Escherichia coli naturally emerges as a result of entropic demixing of such chromosomal loops. We propose one possible mechanism for organizing these domains: regulatory control defined by the underlying E. coli gene regulatory network requires the colocalization of transcription factor genes and target genes. Investigating this assumption, we find the DNA chain to self-organize into several topologically distinguishable domains where the interplay between the entropic repulsion of chromosomal loops and their compression due to the confining geometry induces an effective nucleoid filament-type of structure. Thus, we propose that the physical structure of the chromosome is a direct result of regulatory interactions. To reproduce the observed precise ordering of the chromosome, we estimate that the domain sizes are distributed between 10 and 700 kb, in agreement with the size of topological domains identified in the context of DNA supercoiling.

Genomic arrangement of regulons in bacterial genomes

Regulons, as groups of transcriptionally co-regulated operons, are the basic units of cellular response systems in bacterial cells. While the concept has been long and widely used in bacterial studies since it was first proposed in 1964, very little is known about how its component operons are arranged in a bacterial genome. We present a computational study to elucidate of the organizational principles of regulons in a bacterial genome, based on the experimentally validated regulons of E. coli and B. subtilis. Our results indicate that (1) genomic locations of transcriptional factors (TFs) are under stronger evolutionary constraints than those of the operons they regulate so changing a TF's genomic location will have larger impact to the bacterium than changing the genomic position of any of its target operons; (2) operons of regulons are generally not uniformly distributed in the genome but tend to form a few closely located clusters, which generally consist of genes working in the same metabolic pathways; and (3) the global arrangement of the component operons of all the regulons in a genome tends to minimize a simple scoring function, indicating that the global arrangement of regulons follows simple organizational principles.

Phenotypic landscape of a bacterial cell

The explosion of sequence information in bacteria makes developing high-throughput, cost-effective approaches to matching genes with phenotypes imperative. Using E. coli as proof of principle, we show that combining large-scale chemical genomics with quantitative fitness measurements provides a high-quality data set rich in discovery. Probing growth profiles of a mutant library in hundreds of conditions in parallel yielded > 10,000 phenotypes that allowed us to study gene essentiality, discover leads for gene function and drug action, and understand higher-order organization of the bacterial chromosome. We highlight new information derived from the study, including insights into a gene involved in multiple antibiotic resistance and the synergy between a broadly used combinatory antibiotic therapy, trimethoprim and sulfonamides. This data set, publicly available at http://ecoliwiki.net/tools/chemgen/, is a valuable resource for both the microbiological and bioinformatic communities, as it provides high-confidence associations between hundreds of annotated and uncharacterized genes as well as inferences about the mode of action of several poorly understood drugs.

Chromosomal rearrangements formed by rrn recombination do not improve replichore balance in host-specific Salmonella enterica serovars

Background

Most of the ∼2,600 serovars of Salmonella enterica have a broad host range as well as a conserved gene order. In contrast, some Salmonella serovars are host-specific and frequently exhibit large chromosomal rearrangements from recombination between rrn operons. One hypothesis explaining these rearrangements suggests that replichore imbalance introduced from horizontal transfer of pathogenicity islands and prophages drives chromosomal rearrangements in an attempt to improve balance.

Methodology/Principal Findings

This hypothesis was directly tested by comparing the naturally-occurring chromosomal arrangement types to the theoretically possible arrangement types, and estimating their replichore balance using a calculator. In addition to previously characterized strains belonging to host-specific serovars, the arrangement types of 22 serovar Gallinarum strains was also determined. Only 48 out of 1,440 possible arrangement types were identified in 212 host-specific strains. While the replichores of most naturally-occurring arrangement types were well-balanced, most theoretical arrangement types had imbalanced replichores. Furthermore, the most common types of rearrangements did not change replichore balance.

Conclusions/Significance

The results did not support the hypothesis that replichore imbalance causes these rearrangements, and suggest that the rearrangements could be explained by aspects of a host-specific lifestyle.

Replication-biased genome organisation in the crenarchaeon Sulfolobus

BACKGROUND

Species of the crenarchaeon Sulfolobus harbour three replication origins in their single circular chromosome that are synchronously initiated during replication.

RESULTS

We demonstrate that global gene expression in two Sulfolobus species is highly biased, such that early replicating genome regions are more highly expressed at all three origins. The bias by far exceeds what would be anticipated by gene dosage effects alone. In addition, early replicating regions are denser in archaeal core genes (enriched in essential functions), display lower intergenic distances, and are devoid of mobile genetic elements.

CONCLUSION

The strong replication-biased structuring of the Sulfolobus chromosome implies that the multiple replication origins serve purposes other than simply shortening the time required for replication. The higher-level chromosomal organisation could be of importance for minimizing the impact of DNA damage, and may also be linked to transcriptional regulation.

Genomic arrangement of bacterial operons is constrained by biological pathways encoded in the genome

It is generally known that bacterial genes working in the same biological pathways tend to group into operons, possibly to facilitate cotranscription and to provide stoichiometry. However, very little is understood about what may determine the global arrangement of bacterial genes in a genome beyond the operon level. Here we present evidence that the global arrangement of operons in a bacterial genome is largely influenced by the tendency that a bacterium keeps its operons encoding the same biological pathway in nearby genomic locations, and by the tendency to keep operons involved in multiple pathways in locations close to the other members of their participating pathways. We also observed that the activation frequencies of pathways also influence the genomic locations of their encoding operons, tending to have operons of the more frequently activated pathways more tightly clustered together. We have quantitatively assessed the influences on the global genomic arrangement of operons by different factors. We found that the current arrangements of operons in most of the bacterial genomes we studied tend to minimize the overall distance between consecutive operons of a same pathway across all pathways encoded in the genome.

Systemic analysis of the symbiotic function of Buchnera aphidicola, the primary endosymbiont of the pea aphid Acyrthosiphon pisum

Buchnera aphidicola is the primary obligate intracellular symbiont of most aphid species. B. aphidicola and aphids have been evolving in parallel since their association started, about 150 Myr ago. Both partners have lost their autonomy, and aphid diversification has been confined to smaller ecological niches by this co-evolution. B. aphidicola has undergone major genomic and biochemical changes as a result of adapting to intracellular life. Several genomes of B. aphidicola from different aphid species have been sequenced in the last decade, making it possible to carry out analyses and comparative studies using system-level in silico methods. This review attempts to provide a systemic description of the symbiotic function of aphid endosymbionts, particularly of B. aphidicola from the pea aphid Acyrthosiphon pisum, by analyzing their structural genomic properties, as well as their genetic and metabolic networks.

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.

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.

Chromosome-specific spatial periodicities in gene expression revealed by spectral analysis

Recent years have seen an unprecedented surge of research activity in studies of gene expression. This extensive work, however, has been almost uniformly focused on genome-wide gene expression and has largely ignored the fundamental fact that every gene has a specific chromosome location. We propose a novel method of spectral analysis for detecting hidden periodicities in gene expression signals ordered along the length of each chromosome. Using this method, we have discovered that each chromosome in rodents and humans has a unique periodic pattern of gene expression. The uncovered spatial periodicities in gene expression are tissue-specific in the sense that the largest differences in humans were observed between two normal tissues (brain and mammary gland) as well as between their tumor counterparts (glioma and breast cancer). The smallest differences resulted from the comparison of tumors (glioma and breast cancer) with their normal counterparts. All such effects do not extend to all chromosomes but are limited to only some of them. The estimated periods and amplitudes are identical for the genes located on the positive and negative DNA strands. While precise molecular mechanisms of chromosome-specific periodicities in gene expression have yet to be unraveled, their universal presence in different tissues adds another dimension to the current understanding of the genome organization.

Analyzing stochastic transcription to elucidate the nucleoid's organization

BACKGROUND

The processes of gene transcription, translation, as well as the reactions taking place between gene products, are subject to stochastic fluctuations. These stochastic events are being increasingly examined as it emerges that they can be crucial in the cell's survival. In a previous study we had examined the transcription patterns of two bacterial species (Escherichia coli and Bacillus subtilis) to elucidate the nucleoid's organization. The basic idea is that genes that share transcription patterns, must share some sort of spatial relationship, even if they are not close to each other on the chromosome. We had found that picking any gene at random, its transcription will be correlated with genes at well-defined short - as well as long-range distances, leaving the explanation of the latter an open question. In this paper we study the transcription correlations when the only transcription taking place is stochastic, in other words, no active or "deterministic" transcription takes place. To this purpose we use transcription data of Sinorhizobium meliloti.

RESULTS

Even when only stochastic transcription takes place, the co-expression of genes varies as a function of the distance between genes: we observe again the short-range as well as the regular, long-range correlation patterns.

CONCLUSION

We explain these latter with a model based on the physical constraints acting on the DNA, forcing it into a conformation of groups of a few successive large and transcribed loops, which are evenly spaced along the chromosome and separated by small, non-transcribed loops. We discuss the question about the link between shared transcription patterns and physiological relationship and come to the conclusion that when genes are distantly placed along the chromosome, the transcription correlation does not imply a physiological relationship.

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.

DNA-protein interactions and bacterial chromosome architecture

Bacteria, like eukaryotic organisms, must compact the DNA molecule comprising their genome and form a functional chromosome. Yet, bacteria do it differently. A number of factors contribute to genome compaction and organization in bacteria, including entropic effects, supercoiling and DNA-protein interactions. A gamut of new experimental techniques have allowed new advances in the investigation of these factors, and spurred much interest in the dynamic response of the chromosome to environmental cues, segregation, and architecture, during both exponential and stationary phases. We review these recent developments with emphasis on the multifaceted roles that DNA-protein interactions play.