Exploring the enzymatic repertoires of Bacteria and Archaea and their associations with metabolic maps

The evolution, survival, and adaptation of microbes are consequences of gene duplication, acquisition, and divergence in response to environmental challenges. In this context, enzymes play a central role in the evolution of organisms, because they are fundamental in cell metabolism. Here, we analyzed the enzymatic repertoire in 6,467 microbial genomes, including their abundances, and their associations with metabolic maps. We found that the enzymes follow a power-law distribution, in relation to the genome sizes. Therefore, we evaluated the total proportion enzymatic classes in relation to the genomes, identifying a descending-order proportion: transferases (EC:2.-), hydrolases (EC:3.-), oxidoreductases (EC:1.-), ligases (EC:6.-), lyases (EC:4.-), isomerases (EC:5.-), and translocases (EC:7-.). In addition, we identified a preferential use of enzymatic classes in metabolism pathways for xenobiotics, cofactors and vitamins, carbohydrates, amino acids, glycans, and energy. Therefore, this analysis provides clues about the functional constraints associated with the enzymatic repertoire of functions in Bacteria and Archaea.

Homology-based reconstruction of regulatory networks for bacterial and archaeal genomes

Gene regulation is a key process for all microorganisms, as it allows them to adapt to different environmental stimuli. However, despite the relevance of gene expression control, for only a handful of organisms is there related information about genome regulation. In this work, we inferred the gene regulatory networks (GRNs) of bacterial and archaeal genomes by comparisons with six organisms with well-known regulatory interactions. The references we used are: Escherichia coli K-12 MG1655, Bacillus subtilis 168, Mycobacterium tuberculosis, Pseudomonas aeruginosa PAO1, Salmonella enterica subsp. enterica serovar typhimurium LT2, and Staphylococcus aureus N315. To this end, the inferences were achieved in two steps. First, the six model organisms were contrasted in an all-vs-all comparison of known interactions based on Transcription Factor (TF)-Target Gene (TG) orthology relationships and Transcription Unit (TU) assignments. In the second step, we used a guilt-by-association approach to infer the GRNs for 12,230 bacterial and 649 archaeal genomes based on TF-TG orthology relationships of the six bacterial models determined in the first step. Finally, we discuss examples to show the most relevant results obtained from these inferences. A web server with all the predicted GRNs is available at https://regulatorynetworks.unam.mx/ or http://132.247.46.6/.

Comparative Genomic Analysis Reveals Potential Pathogenicity and Slow-Growth Characteristics of Genus Brevundimonas and Description of Brevundimonas pishanensis sp. nov

The genus Brevundimonas consists of Gram-negative bacteria widely distributed in environment and can cause human infections. However, the genomic characteristics and pathogenicity of Brevundimonas remain poorly studied. Here, the whole-genome features of 24 Brevundimonas type strains were described. Brevundimonas spp. had relatively small genomes (3.13 ± 0.29 Mb) within the family Caulobacteraceae but high G+C contents (67.01 ± 2.19 mol%). Two-dimensional hierarchical clustering divided those genomes into 5 major clades, in which clades II and V contained nine and five species, respectively. Interestingly, phylogenetic analysis showed a one-to-one match between core and accessory genomes, which suggested coevolution of species within the genus Brevundimonas. The unique genes were annotated to biological functions like catalytic activity, signaling and cellular processes, multisubstance metabolism, etc. The majority of Brevundimonas spp. harbored virulence-associated genes icl, tufA, kdsA, htpB, and acpXL, which encoded isocitrate lyase, elongation factor, 2-dehydro-3-deoxyphosphooctonate aldolase, heat shock protein, and acyl carrier protein, respectively. In addition, genomic islands (GIs) and phages/prophages were identified within the Brevundimonas genus. Importantly, a novel Brevundimonas species was identified from the feces of a patient (suffering from diarrhea) by the analyses of biochemical characteristics, phylogenetic tree of 16S rRNA gene, multilocus sequence analysis (MLSA) sequences, and genomic data. The name Brevundimonas pishanensis sp. nov. was proposed, with type strain CHPC 1.3453 (= GDMCC 1.2503^T= KCTC 82824^T). Brevundimonas spp. also showed obvious slow growth compared with that of Escherichia coli. Our study reveals insights into genomic characteristics and potential virulence-associated genes of Brevundimonas spp., and provides a basis for further intensive study of the pathogenicity of Brevundimonas.

IMPORTANCEBrevundimonas spp., a group of bacteria from the family Caulobacteraceae, is associated with nosocomial infections, deserve widespread attention. Our study elucidated genes potentially associated with the pathogenicity of the Brevundimonas genus. We also described some new characteristics of Brevundimonas spp., such as small chromosome size, high G+C content, and slow-growth phenotypes, which made the Brevundimonas genus a good model organism for in-depth studies of growth rate traits. Apart from the comparative analysis of the genomic features of the Brevundimonas genus, we also reported a novel Brevundimonas species, Brevundimonas pishanensis, from the feces of a patient with diarrhea. Our study promotes the understanding of the pathogenicity characteristics of Brevundimonas species bacteria.

Complementary Tendencies in the Use of Regulatory Elements (Transcription Factors, Sigma Factors, and Riboswitches) in Bacteria and Archaea

In prokaryotes, the key players in transcription initiation are sigma factors and transcription factors that bind to DNA to modulate the process, while premature transcription termination at the 5' end of the genes is regulated by attenuation and, in particular, by attenuation associated with riboswitches. In this study, we describe the distribution of these regulators across phylogenetic groups of bacteria and archaea and find that their abundance not only depends on the genome size, as previously described, but also varies according to the phylogeny of the organism. Furthermore, we observed a tendency for organisms to compensate for the low frequencies of a particular type of regulatory element (i.e., transcription factors) with a high frequency of other types of regulatory elements (i.e., sigma factors). This study provides a comprehensive description of the more abundant COG, KEGG, and Rfam families of transcriptional regulators present in prokaryotic genomes.IMPORTANCE In this study, we analyzed the relationship between the relative frequencies of the primary regulatory elements in bacteria and archaea, namely, transcription factors, sigma factors, and riboswitches. In bacteria, we reveal a compensatory behavior for transcription factors and sigma factors, meaning that in phylogenetic groups in which the relative number of transcription factors was low, we found a tendency for the number of sigma factors to be high and vice versa. For most of the phylogenetic groups analyzed here, except for Firmicutes and Tenericutes, a clear relationship with other mechanisms was not detected for transcriptional riboswitches, suggesting that their low frequency in most genomes does not constitute a significant impact on the global variety of transcriptional regulatory elements in prokaryotic organisms.

A progesterone biosensor derived from microbial screening

Bacteria are an enormous and largely untapped reservoir of biosensing proteins. We describe an approach to identify and isolate bacterial allosteric transcription factors (aTFs) that recognize a target analyte and to develop these TFs into biosensor devices. Our approach utilizes a combination of genomic screens and functional assays to identify and isolate biosensing TFs, and a quantum-dot Förster Resonance Energy Transfer (FRET) strategy for transducing analyte recognition into real-time quantitative measurements. We use this approach to identify a progesterone-sensing bacterial aTF and to develop this TF into an optical sensor for progesterone. The sensor detects progesterone in artificial urine with sufficient sensitivity and specificity for clinical use, while being compatible with an inexpensive and portable electronic reader for point-of-care applications. Our results provide proof-of-concept for a paradigm of microbially-derived biosensors adaptable to inexpensive, real-time sensor devices.

The Temporal Expression of Global Regulator Protein CsrA Is Dually Regulated by ClpP During the Biphasic Life Cycle of Legionella pneumophila

Legionella pneumophila, an environmental bacterium that parasitizes protozoa, is the causative pathogen of Legionnaires' disease. L. pneumophila adopts a distinct biphasic life cycle that allows it to adapt to environmental conditions for survival, replication, and transmission. This cycle consists of a non-virulent replicative phase (RP) and a virulent transmissive phase (TP). Timely and fine-tuned expression of growth and virulence factors in a life cycle-dependent manner is crucial. Herein, we report evidence that CsrA, a key regulator of the switch between the RP and the TP, is dually regulated in a ClpP-dependent manner during the biphasic life cycle of L. pneumophila. First, we show that the protein level of CsrA is temporal during the life cycle and is degraded by ClpP during the TP. The ectopic expression of CsrA in a ΔclpP mutant, but not in the wild type, inhibits both the initiation of the RP in vitro and the invasiveness to Acanthamoeba castellanii, indicating that the ClpP-mediated proteolytic pathway regulates the CsrA protein level. We further show that the temporally expressed IHFB is the transcriptional inhibitor of csrA and is degraded via a ClpP-dependent manner during the RP. During the RP, the level of CsrA is increased by promoting the degradation of IHFB and reducing the degradation of the accumulated CsrA via a ClpP-dependent manner. During the TP, the level of CsrA is decreased by inhibiting the degradation of IHFB and promoting the degradation of the accumulated CsrA via a ClpP-dependent manner as well. In conclusion, our results show that the growth-stage-specific expression level of CsrA is dually regulated by ClpP-dependent proteolysis at both the transcription and protein levels during the biphasic life cycle of L. pneumophila.

Diversity, versatility and complexity of bacterial gene regulation mechanisms: opportunities and drawbacks for applications in synthetic biology

Gene expression occurs in two essential steps: transcription and translation. In bacteria, the two processes are tightly coupled in time and space, and highly regulated. Tight regulation of gene expression is crucial. It limits wasteful consumption of resources and energy, prevents accumulation of potentially growth inhibiting reaction intermediates, and sustains the fitness and potential virulence of the organism in a fluctuating, competitive and frequently stressful environment. Since the onset of studies on regulation of enzyme synthesis, numerous distinct regulatory mechanisms modulating transcription and/or translation have been discovered. Mostly, various regulatory mechanisms operating at different levels in the flow of genetic information are used in combination to control and modulate the expression of a single gene or operon. Here, we provide an extensive overview of the very diverse and versatile bacterial gene regulatory mechanisms with major emphasis on their combined occurrence, intricate intertwinement and versatility. Furthermore, we discuss the potential of well-characterized basal expression and regulatory elements in synthetic biology applications, where they may ensure orthogonal, predictable and tunable expression of (heterologous) target genes and pathways, aiming at a minimal burden for the host.

Abundance, diversity and domain architecture variability in prokaryotic DNA-binding transcription factors

Gene regulation at the transcriptional level is a central process in all organisms, and DNA-binding transcription factors, known as TFs, play a fundamental role. This class of proteins usually binds at specific DNA sequences, activating or repressing gene expression. In general, TFs are composed of two domains: the DNA-binding domain (DBD) and an extra domain, which in this work we have named “companion domain” (CD). This latter could be involved in one or more functions such as ligand binding, protein-protein interactions or even with enzymatic activity. In contrast to DBDs, which have been widely characterized both experimentally and bioinformatically, information on the abundance, distribution, variability and possible role of the CDs is scarce. Here, we investigated these issues associated with the domain architectures of TFs in prokaryotic genomes. To this end, 19 families of TFs in 761 non-redundant bacterial and archaeal genomes were evaluated. In this regard we found four main groups based on the abundance and distribution in the analyzed genomes: i) LysR and TetR/AcrR; ii) AraC/XylS, SinR, and others; iii) Lrp, Fis, ArsR, and others; and iv) a group that included only two families, ArgR and BirA. Based on a classification of the organisms according to the life-styles, a major abundance of regulatory families in free-living organisms, in contrast with pathogenic, extremophilic or intracellular organisms, was identified. Finally, the protein architecture diversity associated to the 19 families considering a weight score for domain promiscuity evidenced which regulatory families were characterized by either a large diversity of CDs, here named as “promiscuous” families given the elevated number of variable domains found in those TFs, or a low diversity of CDs. Altogether this information helped us to understand the diversity and distribution of the 19 Prokaryotes TF families. Moreover, initial steps were taken to comprehend the variability of the extra domain in those TFs, which eventually might assist in evolutionary and functional studies.

The functional landscape bound to the transcription factors of Escherichia coli K-12

Motivated by the experimental evidences accumulated in the last ten years and based on information deposited in RegulonDB, literature look up, and sequence analysis, we analyze the repertoire of 304 DNA-binding Transcription factors (TFs) in Escherichia coli K-12. These regulators were grouped in 78 evolutionary families and are regulating almost half of the total genes in this bacterium. In structural terms, 60% of TFs are composed by two-domains, 30% are monodomain, and 10% three- and four-structural domains. As previously noticed, the most abundant DNA-binding domain corresponds to the winged helix-turn-helix, with few alternative DNA-binding structures, resembling the hypothesis of successful protein structures with the emergence of new ones at low scales. In summary, we identified and described the characteristics associated to the DNA-binding TF in E. coli K-12. We also identified twelve functional modules based on a co-regulated gene matrix. Finally, diverse regulons were predicted based on direct associations between the TFs and potential regulated genes. This analysis should increase our knowledge about the gene regulation in the bacterium E. coli K-12, and provide more additional clues for comprehensive modelling of transcriptional regulatory networks in other bacteria.

Increments and duplication events of enzymes and transcription factors influence metabolic and regulatory diversity in prokaryotes

In this work, the content of enzymes and DNA-binding transcription factors (TFs) in 794 non-redundant prokaryotic genomes was evaluated. The identification of enzymes was based on annotations deposited in the KEGG database as well as in databases of functional domains (COG and PFAM) and structural domains (Superfamily). For identifications of the TFs, hidden Markov profiles were constructed based on well-known transcriptional regulatory families. From these analyses, we obtained diverse and interesting results, such as the negative rate of incremental changes in the number of detected enzymes with respect to the genome size. On the contrary, for TFs the rate incremented as the complexity of genome increased. This inverse related performance shapes the diversity of metabolic and regulatory networks and impacts the availability of enzymes and TFs. Furthermore, the intersection of the derivatives between enzymes and TFs was identified at 9,659 genes, after this point, the regulatory complexity grows faster than metabolic complexity. In addition, TFs have a low number of duplications, in contrast to the apparent high number of duplications associated with enzymes. Despite the greater number of duplicated enzymes versus TFs, the increment by which duplicates appear is higher in TFs. A lower proportion of enzymes among archaeal genomes (22%) than in the bacterial ones (27%) was also found. This low proportion might be compensated by the interconnection between the metabolic pathways in Archaea. A similar proportion was also found for the archaeal TFs, for which the formation of regulatory complexes has been proposed. Finally, an enrichment of multifunctional enzymes in Bacteria, as a mechanism of ecological adaptation, was detected.

The repertoire of DNA-binding transcription factors in prokaryotes: functional and evolutionary lessons

The capabilities of organisms to contend with environmental changes depend on their genes and their ability to regulate their expression. DNA-binding transcription factors (TFs) play a central role in this process, because they regulate gene expression positively and/or negatively, depending on the operator context and ligand-binding status. In this review, we summarise recent findings regarding the function and evolution of TFs in prokaryotes. We consider the abundance of TFs in bacteria and archaea, the role of DNA-binding domains and their partner domains, and the effects of duplication events in the evolution of regulatory networks. Finally, a comprehensive picture for how regulatory networks have evolved in prokaryotes is provided.

Joint scaling laws in functional and evolutionary categories in prokaryotic genomes

We propose and study a class-expansion/innovation/loss model of genome evolution taking into account biological roles of genes and their constituent domains. In our model, numbers of genes in different functional categories are coupled to each other. For example, an increase in the number of metabolic enzymes in a genome is usually accompanied by addition of new transcription factors regulating these enzymes. Such coupling can be thought of as a proportional ‘recipe’ for genome composition of the type ‘a spoonful of sugar for each egg yolk’. The model jointly reproduces two known empirical laws: the distribution of family sizes and the non-linear scaling of the number of genes in certain functional categories (e.g. transcription factors) with genome size. In addition, it allows us to derive a novel relation between the exponents characterizing these two scaling laws, establishing a direct quantitative connection between evolutionary and functional categories. It predicts that functional categories that grow faster-than-linearly with genome size to be characterized by flatter-than-average family size distributions. This relation is confirmed by our bioinformatics analysis of prokaryotic genomes. This proves that the joint quantitative trends of functional and evolutionary classes can be understood in terms of evolutionary growth with proportional recipes.

The Regulatory Network of Pseudomonas aeruginosa

Background

Pseudomonas aeruginosa is an important bacterial model due to its metabolic and pathogenic abilities, which allow it to interact and colonize a wide range of hosts, including plants and animals. In this work we compile and analyze the structure and organization of an experimentally supported regulatory network in this bacterium.

Results

The regulatory network consists of 690 genes and 1020 regulatory interactions between their products (12% of total genes: 54% sigma and 16% of transcription factors). This complex interplay makes the third largest regulatory network of those reported in bacteria. The entire network is enriched for activating interactions and, peculiarly, self-activation seems to occur more prominent for transcription factors (TFs), which contrasts with other biological networks where self-repression is dominant. The network contains a giant component of 650 genes organized into 11 hierarchies, encompassing important biological processes, such as, biofilms formation, production of exopolysaccharide alginate and several virulence factors, and of the so-called quorum sensing regulons.

Conclusions

The study of gene regulation in P. aeruginosa is biased towards pathogenesis and virulence processes, all of which are interconnected. The network shows power-law distribution -input degree -, and we identified the top ten global regulators, six two-element cycles, the longest paths have ten steps, six biological modules and the main motifs containing three and four elements. We think this work can provide insights for the design of further studies to cover the many gaps in knowledge of this important bacterial model, and for the design of systems strategies to combat this bacterium.

Diversity and distribution of transcription factors: their partner domains play an important role in regulatory plasticity in bacteria

The ability of bacteria to deal with diverse environmental changes depends on their repertoire of genes and their ability to regulate their expression. In this process, DNA-binding transcription factors (TFs) have a fundamental role because they affect gene expression positively and/or negatively depending on operator context and ligand-binding status. Here, we show an exhaustive analysis of winged helix-turn-helix domains (wHTHs), a class of DNA-binding TFs. These proteins were identified in high proportions and widely distributed in bacteria, representing around half of the total TFs identified so far. In addition, we evaluated the repertoire of wHTHs in terms of their partner domains (PaDos), identifying a similar trend, as with TFs, i.e. they are abundant and widely distributed in bacteria. Based on the PaDos, we defined three main groups of families: (i) monolithic, those families with little PaDo diversity, such as LysR; (ii) promiscuous, those families with a high PaDo diversity; and (iii) monodomain, with families of small sizes, such as MarR. These findings suggest that PaDos have a very important role in the diversification of regulatory responses in bacteria, probably contributing to their regulatory complexity. Thus, the TFs discriminate over longer regions on the DNA through their diverse DNA-binding domains. On the other hand, the PaDos would allow a great flexibility for transcriptional regulation due to their ability to sense diverse stimuli through a variety of ligand-binding compounds.

An overview of prokaryotic transcription factors : a summary of function and occurrence in bacterial genomes

Transcriptional initiation is arguably the most important control point for gene expression. It is regulated by a combination of factors, including DNA sequence and its three-dimensional topology, proteins and small molecules. In this chapter, we focus on the trans-acting factors of bacterial regulation. Initiation begins with the recruitment of the RNA polymerase holoenzyme to a specific locus upstream of the gene known as its promoter. The sigma factor, which is a component of the holoenzyme, provides the most fundamental mechanisms for orchestrating broad changes in gene expression state. It is responsible for promoter recognition as well as recruiting the holoenzyme to the promoter. Distinct sigma factors compete with for binding to a common pool of RNA polymerases, thus achieving condition-dependent differential expression. Another important class of bacterial regulators is transcription factors, which activate or repress transcription of target genes typically in response to an environmental or cellular trigger. These factors may be global or local depending on the number of genes and range of cellular functions that they target. The activities of both global and local transcription factors may be regulated either at a post-transcriptional level via signal-sensing protein domains or at the level of their own expression. In addition to modulating polymerase recruitment to promoters, several global factors are considered as "nucleoid-associated proteins" that impose structural constraints on the chromosome by altering the conformation of the bound DNA, thus influencing other processes involving DNA such as replication and recombination. This chapter concludes with a discussion of how regulatory interactions between transcription factors and their target genes can be represented as a network.

Structural and functional map of a bacterial nucleoid

Mapping global protein binding in the E. coli genome reveals extended domains of high protein occupancy.

Genome-wide mapping of transcription factor-DNA interactions in bacterial chromosomes in vivo has begun to reveal global zones occupied by these factors that serve two purposes: compacting the bacterial DNA and influencing global programs of gene transcription.