Predicting regulons and their cis-regulatory motifs by comparative genomics

We have combined and compared three techniques for predicting functional interactions based on comparative genomics (methods based on conserved operons, protein fusions and correlated evolution) and optimized these methods to predict coregulated sets of genes in 24 complete genomes, including Saccharomyces cerevisiae, Caernorhabditis elegans and 22 prokaryotes. The method based on conserved operons was the most useful for this purpose. Upstream regions of the genes comprising these predicted regulons were then used to search for regulatory motifs in 22 prokaryotic genomes using the motif-discovery program AlignACE. Many significant upstream motifs, including five known Escherichia coli regulatory motifs, were identified in this manner. The presence of a significant regulatory motif was used to refine the members of the predicted regulons to generate a final set of predicted regulons that share significant regulatory elements.     

Exploring environmental intra-species diversity through non-redundant pangenome assemblies

At the genome level, microorganisms are highly adaptable both in terms of allele and gene composition. Such heritable traits emerge in response to different environmental niches and can have a profound influence on microbial community dynamics. As a consequence, any individual genome or population will contain merely a fraction of the total genetic diversity of any operationally defined “species”, whose ecological potential can thus be only fully understood by studying all of their genomes and the genes therein. This concept, known as the pangenome, is valuable for studying microbial ecology and evolution, as it partitions genomes into core (present in all the genomes from a species, and responsible for housekeeping and species-level niche adaptation among others) and accessory regions (present only in some, and responsible for intra-species differentiation). Here we present SuperPang, an algorithm producing pangenome assemblies from a set of input genomes of varying quality, including metagenome-assembled genomes (MAGs). SuperPang runs in linear time and its results are complete, non-redundant, preserve gene ordering and contain both coding and non-coding regions. Our approach provides a modular view of the pangenome, identifying operons and genomic islands, and allowing to track their prevalence in different populations. We illustrate this by analysing intra-species diversity in Polynucleobacter, a bacterial genus ubiquitous in freshwater ecosystems, characterized by their streamlined genomes and their ecological versatility. We show how SuperPang facilitates the simultaneous analysis of allelic and gene content variation under different environmental pressures, allowing us to study the drivers of microbial diversification at unprecedented resolution.     

Diverse bacteria isolated from root nodules of <Emphasis Type="Italic">Trifolium</Emphasis>, <Emphasis Type="Italic">Crotalaria</Emphasis> and <Emphasis Type="Italic">Mimosa</Emphasis> grown in the subtropical regions of China

We investigated the regulation of the two of the three groE operons (cpn.1 and cpn.2) of the root-nodulating bacterium R. leguminosarum strain A34. Both are heat inducible, and both have a CIRCE sequence in their upstream regions, suggesting regulation by an HrcA repressor. Mutagenesis of the CIRCE sequence upstream of cpn.1 led to an increase in the levels of cpn.1 mRNA, and knock-out of the hrcA gene increased the level of Cpn60.1 protein (the GroEL homologue encoded by the cpn.1 operon). Inactivation of the hrcA gene also caused increased expression of a 29 kDa protein that was identified as RhiA, a component of a quorum-sensing system. However, neither loss of the upstream CIRCE sequence, nor loss of HrcA function, had any effect on expression from the cpn.2 promoter. Further analysis of the cpn.2 upstream region suggested regulation could be mediated by an RpoH system, and this was confirmed by deleting the rpoH gene from the chromosome, which led to a decreased level of Cpn60.2 expression. Inactivation of RpoH led to a reduction in growth rate which could be partly compensated for by inactivation of HrcA, indicating an overlap in the in vivo function of the proteins regulated by these two systems. Accession numbers: DQ173160 (hrcA operon); DQ173161 (rpoH gene).     

The Evolution of Histidine Biosynthesis in Archaea: Insights into the <Emphasis Type="Italic">his</Emphasis> Genes Structure and Organization in LUCA

The available sequences of genes encoding the enzymes associated with histidine biosynthesis suggest that this is an ancient metabolic pathway that was assembled prior to the diversification of Bacteria, Archaea, and Eucarya. Paralogous duplication, gene elongation, and fusion events of several different his genes have played a major role in shaping this biosynthetic route. We have analyzed the structure and organization of histidine biosynthetic genes from 55 complete archaeal genomes and combined it with phylogenetic inference in order to investigate the mechanisms responsible for the assembly of the his pathway and the origin of his operons. We show that a wide variety of different organizations of his genes exists in Archaea and that some his genes or entire his (sub-)operons have been likely transferred horizontally between Archaea and Bacteria. However, we show that, in most Archaea, his genes are monofunctional (except for hisD) and scattered throughout the genome, suggesting that his operons might have been assembled multiple times during evolution and that in some cases they are the result of recent evolutionary events. An evolutionary model for the structure and organization of his genes in LUCA is proposed.