Genomic organization of evolutionarily correlated genes in bacteria: limits and strategies

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.

The supragenic organization of glycoside hydrolase encoding genes reveals distinct strategies for carbohydrate utilization in bacteria

Glycoside hydrolases (GHs) are carbohydrate-active enzymes essential for many environmental (e.g., carbon cycling) and biotechnological (e.g., biofuels) processes. The complete processing of carbohydrates by bacteria requires many enzymes acting synergistically. Here, I investigated the clustered or scattered distribution of 406,337 GH-genes and their association with transporter genes identified in 15,640 completely sequenced bacterial genomes. Different bacterial lineages displayed conserved levels of clustered or scattered GH-genes, but overall, the GH-genes clustering was generally higher than in randomized genomes. In lineages with highly clustered GH-genes (e.g., Bacteroides, Paenibacillus), clustered genes shared the same orientation. These codirectional gene clusters potentially facilitate the genes' co-expression by allowing transcriptional read-through and, at least in some cases, forming operons. In several taxa, the GH-genes clustered with distinct types of transporter genes. The type of transporter genes and the distribution of the so-called GH:TR-genes clusters were conserved in selected lineages. Globally, the phylogenetically conserved clustering of the GH-genes with transporter genes highlights the central function of carbohydrate processing across bacterial lineages. In addition, in bacteria with the most identified GH-genes, the genomic adaptations for carbohydrate processing also mirrored the broad environmental origin of the sequenced strains (e.g., soil and mammal gut) suggesting that a combination of evolutionary history and the environment selects for the specific supragenic organization of the GH-genes supporting the carbohydrate processing in bacterial genomes.

Using genome-wide expression compendia to study microorganisms

A gene expression compendium is a heterogeneous collection of gene expression experiments assembled from data collected for diverse purposes. The widely varied experimental conditions and genetic backgrounds across samples creates a tremendous opportunity for gaining a systems level understanding of the transcriptional responses that influence phenotypes. Variety in experimental design is particularly important for studying microbes, where the transcriptional responses integrate many signals and demonstrate plasticity across strains including response to what nutrients are available and what microbes are present. Advances in high-throughput measurement technology have made it feasible to construct compendia for many microbes. In this review we discuss how these compendia are constructed and analyzed to reveal transcriptional patterns.

Selection, periodicity and potential function for Highly Iterative Palindrome-1 (HIP1) in cyanobacterial genomes

Highly Iterated Palindrome 1 (HIP1, GCGATCGC) is hyper-abundant in most cyanobacterial genomes. In some cyanobacteria, average HIP1 abundance exceeds one motif per gene. Such high abundance suggests a significant role in cyanobacterial biology. However, 20 years of study have not revealed whether HIP1 has a function, much less what that function might be. We show that HIP1 is 15- to 300-fold over-represented in genomes analyzed. More importantly, HIP1 sites are conserved both within and between open reading frames, suggesting that their overabundance is maintained by selection rather than by continual replenishment by neutral processes, such as biased DNA repair. This evidence for selection suggests a functional role for HIP1. No evidence was found to support a functional role as a peptide or RNA motif or a role in the regulation of gene expression. Rather, we demonstrate that the distribution of HIP1 along cyanobacterial chromosomes is significantly periodic, with periods ranging from 10 to 90 kb, consistent in scale with periodicities reported for co-regulated, co-expressed and evolutionarily correlated genes. The periodicity we observe is also comparable in scale to chromosomal interaction domains previously described in other bacteria. In this context, our findings imply HIP1 functions associated with chromosome and nucleoid structure.

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.

Broad-scale redistribution of mRNA abundance and transcriptional machinery in response to growth rate in Salmonella enterica serovar Typhimurium

We have investigated the connection between the four-dimensional architecture of the bacterial nucleoid and the organism's global gene expression programme. By localizing the transcription machinery and the transcriptional outputs across the genome of the model bacterium Salmonella enterica serovar Typhimurium at different stages of the growth cycle, a surprising disconnection between gene dosage and transcriptional output was revealed. During exponential growth, gene output occurred chiefly in the Ori (origin), Ter (terminus) and NSL (non-structured left) domains, whereas the Left macrodomain remained transcriptionally quiescent at all stages of growth. The apparently high transcriptional output in Ter was correlated with an enhanced stability of the RNA expressed there during exponential growth, suggesting that longer mRNA half-lives compensate for low gene dosage. During exponential growth, RNA polymerase (RNAP) was detected everywhere, whereas in stationary phase cells, RNAP was concentrated in the Ter macrodomain. The alternative sigma factors RpoE, RpoH and RpoN were not required to drive transcription in these growth conditions, consistent with their observed binding to regions away from RNAP and regions of active transcription. Specifically, these alternative sigma factors were found in the Ter macrodomain during exponential growth, whereas they were localized at the Ori macrodomain in stationary phase.

Bacterial pathogen gene regulation: a DNA-structure-centred view of a protein-dominated domain

The mechanisms used by bacterial pathogens to regulate the expression of their genes, especially their virulence genes, have been the subject of intense investigation for several decades. Whole genome sequencing projects, together with more targeted studies, have identified hundreds of DNA-binding proteins that contribute to the patterns of gene expression observed during infection as well as providing important insights into the nature of the gene products whose expression is being controlled by these proteins. Themes that have emerged include the importance of horizontal gene transfer to the evolution of pathogens, the need to impose regulatory discipline upon these imported genes and the important roles played by factors normally associated with the organization of genome architecture as regulatory principles in the control of virulence gene expression. Among these architectural elements is the structure of DNA itself, its variable nature at a topological rather than just at a base-sequence level and its ability to play an active (as well as a passive) part in the gene regulation process.

From multiple pathogenicity islands to a unique organized pathogenicity archipelago

Pathogenicity islands are sets of successive genes in a genome that determine the virulence of a bacterium. In a growing number of studies, bacterial virulence appears to be determined by multiple islands scattered along the genome. This is the case in a family of seven plant pathogens and a human pathogen that, under KdgR regulation, massively secrete enzymes such as pectinases that degrade plant cell wall. Here we show that their multiple pathogenicity islands form together a coherently organized, single “archipelago” at the genome scale. Furthermore, in half of the species, most genes encoding secreted pectinases are expressed from the same DNA strand (transcriptional co-orientation). This genome architecture favors DNA conformations that are conducive to genes spatial co-localization, sometimes complemented by co-orientation. As proteins tend to be synthetized close to their encoding genes in bacteria, we propose that this architecture would favor the efficient funneling of pectinases at convergent points within the cell. The underlying functional hypothesis is that this convergent funneling of the full blend of pectinases constitutes a crucial strategy for successful degradation of the plant cell wall. Altogether, our work provides a new approach to describe and predict, at the genome scale, the full virulence complement.

Analysis tools for the interplay between genome layout and regulation

BACKGROUND

Genome layout and gene regulation appear to be interdependent. Understanding this interdependence is key to exploring the dynamic nature of chromosome conformation and to engineering functional genomes. Evidence for non-random genome layout, defined as the relative positioning of either co-functional or co-regulated genes, stems from two main approaches. Firstly, the analysis of contiguous genome segments across species, has highlighted the conservation of gene arrangement (synteny) along chromosomal regions. Secondly, the study of long-range interactions along a chromosome has emphasised regularities in the positioning of microbial genes that are co-regulated, co-expressed or evolutionarily correlated. While one-dimensional pattern analysis is a mature field, it is often powerless on biological datasets which tend to be incomplete, and partly incorrect. Moreover, there is a lack of comprehensive, user-friendly tools to systematically analyse, visualise, integrate and exploit regularities along genomes.

RESULTS

Here we present the Genome REgulatory and Architecture Tools SCAN (GREAT:SCAN) software for the systematic study of the interplay between genome layout and gene expression regulation.

GREAT

SCAN is a collection of related and interconnected applications currently able to perform systematic analyses of genome regularities as well as to improve transcription factor binding sites (TFBS) and gene regulatory network predictions based on gene positional information.

CONCLUSIONS

We demonstrate the capabilities of these tools by studying on one hand the regular patterns of genome layout in the major regulons of the bacterium Escherichia coli. On the other hand, we demonstrate the capabilities to improve TFBS prediction in microbes. Finally, we highlight, by visualisation of multivariate techniques, the interplay between position and sequence information for effective transcription regulation.

Comparative Genomics of Interreplichore Translocations in Bacteria: A Measure of Chromosome Topology?

Genomes evolve not only in base sequence but also in terms of their architecture, defined by gene organization and chromosome topology. Whereas genome sequence data inform us about the changes in base sequences for a large variety of organisms, the study of chromosome topology is restricted to a few model organisms studied using microscopy and chromosome conformation capture techniques. Here, we exploit whole genome sequence data to study the link between gene organization and chromosome topology in bacteria. Using comparative genomics across ∼250 pairs of closely related bacteria we show that: (a) many organisms show a high degree of interreplichore translocations throughout the chromosome and not limited to the inversion-prone terminus (ter) or the origin of replication (oriC); (b) translocation maps may reflect chromosome topologies; and (c) symmetric interreplichore translocations do not disrupt the distance of a gene from oriC or affect gene expression states or strand biases in gene densities. In summary, we suggest that translocation maps might be a first line in defining a gross chromosome topology given a pair of closely related genome sequences.

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.

Coevolution of the Organization and Structure of Prokaryotic Genomes

The cytoplasm of prokaryotes contains many molecular machines interacting directly with the chromosome. These vital interactions depend on the chromosome structure, as a molecule, and on the genome organization, as a unit of genetic information. Strong selection for the organization of the genetic elements implicated in these interactions drives replicon ploidy, gene distribution, operon conservation, and the formation of replication-associated traits. The genomes of prokaryotes are also very plastic with high rates of horizontal gene transfer and gene loss. The evolutionary conflicts between plasticity and organization lead to the formation of regions with high genetic diversity whose impact on chromosome structure is poorly understood. Prokaryotic genomes are remarkable documents of natural history because they carry the imprint of all of these selective and mutational forces. Their study allows a better understanding of molecular mechanisms, their impact on microbial evolution, and how they can be tinkered in synthetic biology.

From structure to function of bacterial chromosomes: Evolutionary perspectives and ideas for new experiments

The link between chromosome structure and function is a challenging open question because chromosomes in vivo are highly dynamic and arduous to manipulate. Here, we examine several promising approaches to tackle this question specifically in bacteria, by integrating knowledge from different sources. Toward this end, we first provide a brief overview of experimental tools that have provided insights into the description of the bacterial chromosome, including genetic, biochemical and fluorescence microscopy techniques. We then explore the possibility of using comparative genomics to isolate functionally important features of chromosome organization, exploiting the fact that features shared between phylogenetically distant bacterial species reflect functional significance. Finally, we discuss possible future perspectives from the field of experimental evolution. Specifically, we propose novel experiments in which bacteria could be screened and selected on the basis of the structural properties of their chromosomes.

Genomic location of the major ribosomal protein gene locus determines Vibrio cholerae global growth and infectivity

The effects on cell physiology of gene order within the bacterial chromosome are poorly understood. In silico approaches have shown that genes involved in transcription and translation processes, in particular ribosomal protein (RP) genes, localize near the replication origin (oriC) in fast-growing bacteria suggesting that such a positional bias is an evolutionarily conserved growth-optimization strategy. Such genomic localization could either provide a higher dosage of these genes during fast growth or facilitate the assembly of ribosomes and transcription foci by keeping physically close the many components of these macromolecular machines. To explore this, we used novel recombineering tools to create a set of Vibrio cholerae strains in which S10-spec-α (S10), a locus bearing half of the ribosomal protein genes, was systematically relocated to alternative genomic positions. We show that the relative distance of S10 to the origin of replication tightly correlated with a reduction of S10 dosage, mRNA abundance and growth rate within these otherwise isogenic strains. Furthermore, this was accompanied by a significant reduction in the host-invasion capacity in Drosophila melanogaster. Both phenotypes were rescued in strains bearing two S10 copies highly distal to oriC, demonstrating that replication-dependent gene dosage reduction is the main mechanism behind these alterations. Hence, S10 positioning connects genome structure to cell physiology in Vibrio cholerae. Our results show experimentally for the first time that genomic positioning of genes involved in the flux of genetic information conditions global growth control and hence bacterial physiology and potentially its evolution.

Increasing evidence indicates that nucleoid spatiotemporal organization is crucial for bacterial physiology since these microorganism lack compartmentalized nucleus. However, it is still unclear how gene order within the chromosome can influence cell physiology. Here, by systematically relocating ribosomal protein genes to different genomic positions in Vibrio cholerae, we revealed drastic differences in growth rate and infectivity of this isogenic strain set. We show that genomic positioning of ribosomal protein genes is crucial for physiology by providing replication-dependent higher dosage. Therefore it might play a key role in genome evolution of bacterial species. This work will contribute to discover genomic rules governing cell physiology which will be essential in the context of the creation of new artificial life forms.

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.

Integration of syntactic and semantic properties of the DNA code reveals chromosomes as thermodynamic machines converting energy into information

Understanding genetic regulation is a problem of fundamental importance. Recent studies have made it increasingly evident that, whereas the cellular genetic regulation system embodies multiple disparate elements engaged in numerous interactions, the central issue is the genuine function of the DNA molecule as information carrier. Compelling evidence suggests that the DNA, in addition to the digital information of the linear genetic code (the semantics), encodes equally important continuous, or analog, information that specifies the structural dynamics and configuration (the syntax) of the polymer. These two DNA information types are intrinsically coupled in the primary sequence organisation, and this coupling is directly relevant to regulation of the genetic function. In this review, we emphasise the critical need of holistic integration of the DNA information as a prerequisite for understanding the organisational complexity of the genetic regulation system.

Non-equilibrium polar localization of proteins in bacterial cells

Many proteins are observed to localize to the poles within bacterial cells. Some bacteria show unipolar localization, yet under different conditions bipolar patterns can emerge. One mechanism for spontaneous polar localization has been shown to involve the combination of protein aggregation and nucleoid occlusion. Whether the different observed patterns represent global energy minima for the cellular system remains to be determined. In this paper we show that for a model consisting only of protein aggregation along with an excluded volume effect due to the DNA polymer, that unipolar patterns are the global energy ground state regardless of protein concentration and DNA density. We extend the model to allow for proteins to be added to the cellular volume at a constant rate and show that bipolar (or multi-foci) patterns emerge as the result of the system being kinetically trapped in a local energy minimum. Lastly we also consider the situation of a growing cell that starts with a pre-existing aggregate at one of the poles and determine conditions under which either unipolar or bipolar patterns can exist at the point when it is ready to divide. This work sheds new interpretations on recently published experimental data and suggests experiments to test whether such a mechanism can drive patterning in bacteria.

Plasmids as scribbling pads for operon formation and propagation

Many bacterial genes are in operons and the process whereby operons are formed is therefore fundamental. To help elucidate this process, we propose in the Scribbling Pad hypothesis that bacteria have been constantly using plasmids for genetic experimentation and, in particular, for the construction of operons. This hypothesis simultaneously solves the problems of the creation of operons and the way operons are propagated. We cite results in the literature to support the hypothesis and make experimental predictions to test it.

PreCisIon: PREdiction of CIS-regulatory elements improved by gene's positION

Conventional approaches to predict transcriptional regulatory interactions usually rely on the definition of a shared motif sequence on the target genes of a transcription factor (TF). These efforts have been frustrated by the limited availability and accuracy of TF binding site motifs, usually represented as position-specific scoring matrices, which may match large numbers of sites and produce an unreliable list of target genes. To improve the prediction of binding sites, we propose to additionally use the unrelated knowledge of the genome layout. Indeed, it has been shown that co-regulated genes tend to be either neighbors or periodically spaced along the whole chromosome. This study demonstrates that respective gene positioning carries significant information. This novel type of information is combined with traditional sequence information by a machine learning algorithm called PreCisIon. To optimize this combination, PreCisIon builds a strong gene target classifier by adaptively combining weak classifiers based on either local binding sequence or global gene position. This strategy generically paves the way to the optimized incorporation of any future advances in gene target prediction based on local sequence, genome layout or on novel criteria. With the current state of the art, PreCisIon consistently improves methods based on sequence information only. This is shown by implementing a cross-validation analysis of the 20 major TFs from two phylogenetically remote model organisms. For Bacillus subtilis and Escherichia coli, respectively, PreCisIon achieves on average an area under the receiver operating characteristic curve of 70 and 60%, a sensitivity of 80 and 70% and a specificity of 60 and 56%. The newly predicted gene targets are demonstrated to be functionally consistent with previously known targets, as assessed by analysis of Gene Ontology enrichment or of the relevant literature and databases.

Gene direction in living organisms

Gene direction, which is important for function, has not been subjected to statistical testing for randomness and for the degree of evolutionary changes. We analyzed 747 sequenced species and 2,061 genomes/chromosomes and detected clear differences in gene direction between kingdoms. All the archaeans, bacteria, and protozoa analyzed have genes characterized mainly by same-direction neighbors (i.e., in head-to-foot or foot-to-head order), with up to 391 genes in tandem in protozoan Leishmania infantum. Fungi and photosynthetic protists have genes characterized by opposite-direction neighbors, except chromosome VII of Ashbya gossypii, a progenitor fungus. The gene direction analysis suggests that the same-direction dominance originated from the last common ancestor of these living organisms, then was strengthened in protozoa, but weakened or lost in fungi, photosynthetic protists and some plants/animals, giving chromosomes/genomes with gene opposite-direction dominance (i.e., towards the random use of both DNA strands).