Carbonic anhydrase (CA) (E.C. 4.2.1.1) is a ubiquitous enzyme catalysing interconversion between CO(2) and bicarbonate. The irregular distribution of the phylogenetically distinct classes of CA in procaryotic genome suggests its complex evolutionary history in procaryotes. Genetic evidence regarding the dispensability of CA under high-CO(2) air in some model organisms indicates that CA-deficient microorganisms can persist in the natural environment by choosing high-CO(2) niches. In this study, we studied the distribution of CA in the genome of Proteobacteria. While a large majority of the genome-sequenced Proteobacteria retained a CA gene(s), intracellular bacterial genera such as Buchnera and Rickettsia contained CA-defective strains. Comparison between CA-retaining and CA- deficient genomes showed the absence of whole coding sequence in some strains and the presence of frameshifted coding sequence in other strains. The evidence suggests that CA is inactivated and lost in some proteobacteria during the course of evolution based on its dispensability.

The α-proteobacterium Wolbachia pipientis infects more than 65% of insect species worldwide and manipulates the host reproductive machinery to enable its own survival. It can live in mutualistic relationships with hosts that cause human disease, including mosquitoes that carry the Dengue virus. Like many other bacteria, Wolbachia contains disulfide bond forming (Dsb) proteins that introduce disulfide bonds into secreted effector proteins. The genome of the Wolbachia strain wMel encodes two DsbA-like proteins sharing just 21% sequence identity to each other, α-DsbA1 and α-DsbA2, and an integral membrane protein, α-DsbB. α-DsbA1 and α-DsbA2 both have a Cys-X-X-Cys active site that, by analogy with Escherichia coli DsbA, would need to be oxidized to the disulfide form to serve as a disulfide bond donor toward substrate proteins. Here we show that the integral membrane protein α-DsbB oxidizes α-DsbA1, but not α-DsbA2. The interaction between α-DsbA1 and α-DsbB is very specific, involving four essential cysteines located in the two periplasmic loops of α-DsbB. In the electron flow cascade, oxidation of α-DsbA1 by α-DsbB is initiated by an oxidizing quinone cofactor that interacts with the cysteine pair in the first periplasmic loop. Oxidizing power is transferred to the second cysteine pair, which directly interacts with α-DsbA1. This reaction is inhibited by a non-catalytic disulfide present in α-DsbA1, conserved in other α-proteobacterial DsbAs but not in γ-proteobacterial DsbAs. This is the first characterization of the integral membrane protein α-DsbB from Wolbachia and reveals that the non-catalytic cysteines of α-DsbA1 regulate the redox relay system in cooperation with α-DsbB.

Homologous recombination is a key pathway found in nearly all bacterial taxa. The recombination complex not only allows bacteria to repair DNA double-strand breaks but also promotes adaption through the exchange of DNA between cells. In Proteobacteria, this process is mediated by the RecBCD complex, which relies on the recognition of a DNA motif named Chi to initiate recombination. The Chi motif has been characterized in Escherichia coli and analogous sequences have been found in several other species from diverse families, suggesting that this mode of action is widespread across bacteria. However, the sequences of Chi-like motifs are known for only five bacterial species: E. coli, Haemophilus influenzae, Bacillus subtilis, Lactococcus lactis, and Staphylococcus aureus. In this study, we detected putative Chi motifs in a large dataset of Proteobacteria and identified four additional motifs sharing high sequence similarity and similar properties to the Chi motif of E. coli in 85 species of Proteobacteria. Most Chi motifs were detected in Enterobacteriaceae and this motif appears well conserved in this family. However, we did not detect Chi motifs for the majority of Proteobacteria, suggesting that different motifs are used in these species. Altogether these results substantially expand our knowledge on the evolution of Chi motifs and on the recombination process in bacteria.

In E. coli, 10 to 15% of growing bacteria produce dimeric chromosomes during DNA replication. These dimers are resolved by XerC and XerD, two tyrosine recombinases that target the 28-nucleotide motif (dif) associated with the chromosome's replication terminus. In streptococci and lactococci, an alternative system is composed of a unique, Xer-like recombinase (XerS) genetically linked to a dif-like motif (dif(SL)) located at the replication terminus. Preliminary observations have suggested that the dif/Xer system is commonly found in bacteria with circular chromosomes but that assumption has not been confirmed in an exhaustive analysis. The aim of the present study was to extensively characterize the dif/Xer system in the proteobacteria, since this taxon accounts for the majority of genomes sequenced to date. To that end, we analyzed 234 chromosomes from 156 proteobacterial species and showed that most species (87.8%) harbor XerC and XerD-like recombinases and a dif-related sequence which (i) is located in non-coding sequences, (ii) is close to the replication terminus (as defined by the cumulative GC skew) (iii) has a palindromic structure, (iv) is encoded by a low G+C content and (v) contains a highly conserved XerD binding site. However, not all proteobacteria display this dif/XerCD system. Indeed, a sub-group of pathogenic epsilon-proteobacteria (including Helicobacter sp and Campylobacter sp) harbors a different recombination system, composed of a single recombinase (XerH) which is phylogenetically distinct from the other Xer recombinases and a motif (dif(H)) sharing homologies with dif(SL). Furthermore, no homologs to dif or Xer recombinases could be detected in small endosymbiont genomes or in certain bacteria with larger chromosomes like the Legionellales. This raises the question of the presence of other chromosomal deconcatenation systems in these species. Our study highlights the complexity of dif/Xer recombinase systems in proteobacteria and paves the way for systematic detection of these components in prokaryotes.

We have performed a computational comparative analysis of six small non-coding RNA (sRNA) families in α-proteobacteria. Members of these families were first identified in the intergenic regions of the nitrogen-fixing endosymbiont S. meliloti by a combined bioinformatics screen followed by experimental verification. Consensus secondary structures inferred from covariance models for each sRNA family evidenced in some cases conserved motifs putatively relevant to the function of trans-encoded base-pairing sRNAs i.e., Hfq-binding signatures and exposed anti Shine-Dalgarno sequences. Two particular family models, namely αr15 and αr35, shared own sub-structural modules with the Rfam model suhB (RF00519) and the uncharacterized sRNA family αr35b, respectively. A third sRNA family, termed αr45, has homology to the cis-acting regulatory element speF (RF00518). However, new experimental data further confirmed that the S. meliloti αr45 representative is an Hfq-binding sRNA processed from or expressed independently of speF, thus refining the Rfam speF model annotation. All the six families have members in phylogenetically related plant-interacting bacteria and animal pathogens of the order of the Rhizobiales, some occurring with high levels of paralogy in individual genomes. In silico and experimental evidences predict differential regulation of paralogous sRNAs in S. meliloti 1021. The distribution patterns of these sRNA families suggest major contributions of vertical inheritance and extensive ancestral duplication events to the evolution of sRNAs in plant-interacting bacteria.

Microbes are a major source of antibiotics, pharmaceuticals, and other bioactive compounds. The production of many specialized microbial metabolites is encoded in biosynthetic gene clusters (BGCs). A challenge associated with natural product discovery is that many BGCs are not expressed under laboratory growth conditions. Here we report a genome-mining approach to discover BGCs with luxR-type quorum sensing (QS) genes, which code for regulatory proteins that control gene expression. Our results show that BGCs linked to genes coding for LuxR-like proteins are widespread in Proteobacteria. In addition, we show that associations between luxR homolog genes and BGCs have evolved independently many times, with functionally diverse gene clusters. Overall, these clusters may provide a source of new natural products for which there is some understanding about how to elicit production. IMPORTANCE Bacteria biosynthesize specialized metabolites with a variety of ecological functions, including defense against other microbes. Genes that code for specialized metabolite biosynthetic enzymes are frequently clustered together. These BGCs are often regulated by a transcription factor encoded within the cluster itself. These pathway-specific regulators respond to a signal or indirectly through other means of environmental sensing. Many specialized metabolites are not produced under laboratory growth conditions, and one reason for this issue is that laboratory growth media lack environmental cues necessary for BGC expression. Here, we report a bioinformatics study that reveals that BGCs are frequently linked to genes coding for LuxR family QS-responsive transcription factors in the phylum Proteobacteria. The products of these luxR homolog-associated gene clusters may serve as a practical source of bioactive metabolites.

Carboxysomes are protein-based organelles essential for carbon fixation in cyanobacteria and proteobacteria. Previously, we showed that the cyanobacterial nucleoid is used to equally space out β-carboxysomes across cell lengths by a two-component system (McdAB) in the model cyanobacterium Synechococcus elongatus PCC 7942. More recently, we found that McdAB systems are widespread among β-cyanobacteria, which possess β-carboxysomes, but are absent in α-cyanobacteria, which possess structurally and phyletically distinct α-carboxysomes. Cyanobacterial α-carboxysomes are thought to have arisen in proteobacteria and then horizontally transferred into cyanobacteria, which suggests that α-carboxysomes in proteobacteria may also lack the McdAB system. Here, using the model chemoautotrophic proteobacterium Halothiobacillus neapolitanus, we show that a McdAB system distinct from that of β-cyanobacteria operates to position α-carboxysomes across cell lengths. We further show that this system is widespread among α-carboxysome-containing proteobacteria and that cyanobacteria likely inherited an α-carboxysome operon from a proteobacterium lacking the mcdAB locus. These results demonstrate that McdAB is a cross-phylum two-component system necessary for positioning both α- and β-carboxysomes. The findings have further implications for understanding the positioning of other protein-based bacterial organelles involved in diverse metabolic processes. PLAIN LANGUAGE SUMMARY: Cyanobacteria are well known to fix atmospheric CO2 into sugars using the enzyme Rubisco. Less appreciated are the carbon-fixing abilities of proteobacteria with diverse metabolisms. Bacterial Rubisco is housed within organelles called carboxysomes that increase enzymatic efficiency. Here we show that proteobacterial carboxysomes are distributed in the cell by two proteins, McdA and McdB. McdA on the nucleoid interacts with McdB on carboxysomes to equidistantly space carboxysomes from one another, ensuring metabolic homeostasis and a proper inheritance of carboxysomes following cell division. This study illuminates how widespread carboxysome positioning systems are among diverse bacteria. Carboxysomes significantly contribute to global carbon fixation; therefore, understanding the spatial organization mechanism shared across the bacterial world is of great interest.

The alpha (alpha) proteobacteria, a very large and diverse group, are presently characterized solely on the basis of 16S rRNA trees, with no known molecular characteristic that is unique to this group. The genomes of three alpha-proteobacteria, Rickettsia prowazekii (RP), Caulobacter crescentus (CC) and Bartonella quintana (BQ), were analyzed in order to search for proteins that are unique to this group.

Symbiotic ß-proteobacteria not only occur in root nodules of legumes but are also found in leaves of certain Rubiaceae. The discovery of bacteria in plants formerly not implicated in endosymbiosis suggests a wider occurrence of plant-microbe interactions. Several ß-proteobacteria of the genus Burkholderia are detected in close association with tropical plants. This interaction has occurred three times independently, which suggest a recent and open plant-bacteria association. The presence or absence of Burkholderia endophytes is consistent on genus level and therefore implies a predictive value for the discovery of bacteria. Only a single Burkholderia species is found in association with a given plant species. However, the endophyte species are promiscuous and can be found in association with several plant species. Most of the endophytes are part of the plant-associated beneficial and environmental group, but others are closely related to B. glathei. This soil bacteria, together with related nodulating and non-nodulating endophytes, is therefore transferred to a newly defined and larger PBE group within the genus Burkholderia.

A comprehensive evolutionary analysis of bacterial genomes implies to identify the hallmark of vertical and non-vertical signals and to discriminate them from the presence of mere phylogenetic noise. In this report we have addressed the impact of factors like the universal distribution of the genes, their essentiality or their functional role in the cell on the inference of vertical signal through phylogenomic methods.

Evolution of new complex biological behaviour tends to arise by novel combinations of existing building blocks. The functional and evolutionary building blocks of the proteome are protein domains, the function of a protein being dependent on its constituent domains. We clustered completely-sequenced proteomes of prokaryotes on the basis of their protein domain content, as defined by Pfam (release 16.0). This revealed that, although there was a correlation between phylogeny and domain content, other factors also have an influence. This observation motivated an investigation of the relationship between an organism's lifestyle and the complement of domains and domain architectures found within its proteome.

Genome editing tools, through the disruption of an organism's native genetic material or the introduction of non-native DNA, facilitate functional investigations to link genotypes to phenotypes. Transposons have been instrumental genetic tools in microbiology, enabling genome-wide, randomized disruption of genes and insertions of new genetic elements. Due to this randomness, identifying and isolating particular transposon mutants (i.e., those with modifications at a genetic locus of interest) can be laborious, often requiring one to sift through hundreds or thousands of mutants. Programmable, site-specific targeting of transposons became possible with recently described CRISPR-associated transposase (CASTs) systems, allowing the streamlined recovery of desired mutants in a single step. Like other CRISPR-derived systems, CASTs can be programmed by guide-RNA that is transcribed from short DNA sequence(s). Here, we describe a CAST system and demonstrate its function in bacteria from three classes of Proteobacteria. A dual plasmid strategy is demonstrated: (i) CAST genes are expressed from a broad-host-range replicative plasmid and (ii) guide-RNA and transposon are encoded on a high-copy, suicidal pUC plasmid. Using our CAST system, single-gene disruptions were performed with on-target efficiencies approaching 100% in Beta- and Gammaproteobacteria (Burkholderia thailandensis and Pseudomonas putida, respectively). We also report a peak efficiency of 45% in the Alphaproteobacterium Agrobacterium fabrum. In B. thailandensis, we performed simultaneous co-integration of transposons at two different target sites, demonstrating CAST's utility in multilocus strategies. The CAST system is also capable of high-efficiency large transposon insertion totaling over 11 kbp in all three bacteria tested. Lastly, the dual plasmid system allowed for iterative transposon mutagenesis in all three bacteria without loss of efficiency. Given these iterative capabilities and large payload capacity, this system will be helpful for genome engineering experiments across several fields of research.

Replicative Helicase DnaB interacts with DnaA, DnaC, DnaG, and DNA polymerase III to commence replication, increase the movement rate of the replication fork, and to assemble part of the primosome. The formation of the replication fork is limited by the ability to load DnaB to the DNA, thus DnaB has shown to be vital to a large extent. In the absence of DnaB, the replication fork is not maintained and in a state of inactivity the replication fork degrades and collapses. To further understand importance of this enzyme from an evolutionary perspective, a genomic analysis DnaB protein sequences, chosen from five Proteobacteria subclasses was performed. Our analysis indicates that, DnaB replicative helicases of Alphaproteobacteria and Epsilonproteobacteria have diverged at an earlier stage from Betaproteobacteria, Deltaproteobacteria and Gammaproteobacteria as well as from one another. Our results were further supported, when we reanalyzed and reconstructed the phylogenetic tree after the inclusion of sequences from Actinobacteria and Firmicute phylum. In addition, Betaproteobacteria, Deltaproteobacteria, and Gammaproteobacteria appear to share a closer common ancestor than from the other two subclasses. The Dot-plot analysis indicated that, the region between amino acid residues 320 to 400 was strongly conserved among all five subclasses.

Comparative genomics approaches are broadly used for analysis of transcriptional regulation in bacterial genomes. In this work, we identified binding sites and reconstructed regulons for 33 orthologous groups of transcription factors (TFs) in 196 reference genomes from 21 taxonomic groups of Proteobacteria. Overall, we predict over 10 600 TF binding sites and identified more than 15 600 target genes for 1896 TFs constituting the studied orthologous groups of regulators. These include a set of orthologues for 21 metabolism-associated TFs from Escherichia coli and/or Shewanella that are conserved in five or more taxonomic groups and several additional TFs that represent non-orthologous substitutions of the metabolic regulators in some lineages of Proteobacteria. By comparing gene contents of the reconstructed regulons, we identified the core, taxonomy-specific and genome-specific TF regulon members and classified them by their metabolic functions. Detailed analysis of ArgR, TyrR, TrpR, HutC, HypR and other amino-acid-specific regulons demonstrated remarkable differences in regulatory strategies used by various lineages of Proteobacteria. The obtained genomic collection of in silico reconstructed TF regulons contains a large number of new regulatory interactions that await future experimental validation. The collection provides a framework for future evolutionary studies of transcriptional regulatory networks in Bacteria. It can be also used for functional annotation of putative metabolic transporters and enzymes that are abundant in the reconstructed regulons.

Bronchiectasis is a debilitating disease with chronic airway infection. Proteobacteria, the dominant phylum, can be detected with high-throughput sequencing.

Most antibiotics are derived from the soil, but their catabolism there, which is necessary to close the antibiotic carbon cycle, remains uncharacterized. We report the first draft genome sequences of soil Proteobacteria identified for subsisting solely on β-lactams as their carbon sources. The genomes encode multiple β-lactamases, although their antibiotic catabolic pathways remain enigmatic.

Perfect short inverted repeats (IRs) are known to be enriched in a variety of bacterial and eukaryotic genomes. Currently, it is unclear whether perfect IRs are conserved over evolutionary time scales. In this study, we aimed to characterize the prevalence and evolutionary conservation of IRs across 20 proteobacterial strains. We first identified IRs in Escherichia coli K-12 substr MG1655 and showed that they are overabundant. We next aimed to test whether this overabundance is reflected in the conservation of IRs over evolutionary time scales. To this end, for each perfect IR identified in E. coli MG1655, we collected orthologous sequences from related proteobacterial genomes. We next quantified the evolutionary conservation of these IRs, that is, the presence of the exact same IR across orthologous regions. We observed high conservation of perfect IRs: out of the 234 examined orthologous regions, 145 were more conserved than expected, which is statistically significant even after correcting for multiple testing. Our results together with previous experimental findings support a model in which imperfect IRs are corrected to perfect IRs in a preferential manner via a template switching mechanism.

The effectiveness of elimination of slightly deleterious mutations depends mainly on drift and recombination frequency. Here we analyze the influence of these two factors on the strength of the purifying selection in mitochondrial and proteobacterial orthologous genes taking into account the differences in the organism lifestyles.

The ability of plant-associated microbes to produce gibberellin A (GA) phytohormones was first described for the fungal rice pathogen Gibberella fujikuroi in the 1930s. Recently the capacity to produce GAs was shown for several bacteria, including symbiotic alpha-proteobacteria (α-rhizobia) and gamma-proteobacteria phytopathogens. All necessary enzymes for GA production are encoded by a conserved operon, which appears to have undergone horizontal transfer between and within these two phylogenetic classes of bacteria. Here the operon was shown to be present and functional in a third class, the beta-proteobacteria, where it is found in several symbionts (β-rhizobia). Conservation of function was examined by biochemical characterization of the enzymes encoded by the operon from Paraburkholderia mimosarum LMG 23256T. Despite the in-frame gene fusion between the short-chain alcohol dehydrogenase/reductase and ferredoxin, the encoded enzymes exhibited the expected activity. Intriguingly, together these can only produce GA9, the immediate precursor to the bioactive GA4, as the cytochrome P450 (CYP115) that catalyzes the final hydroxylation reaction is missing, similar to most α-rhizobia. However, phylogenetic analysis indicates that the operon from β-rhizobia is more closely related to examples from gamma-proteobacteria, which almost invariably have CYP115 and, hence, can produce bioactive GA4. This indicates not only that β-rhizobia acquired the operon by horizontal gene transfer from gamma-proteobacteria, rather than α-rhizobia, but also that they independently lost CYP115 in parallel to the α-rhizobia, further hinting at the possibility of detrimental effects for the production of bioactive GA4 by these symbionts.

Here, we report the draft genome sequence of Pseudomonas sp. strain DrBHI1. The total assembly length is 5,649,751 bp in 146 contigs. This strain was isolated from zebrafish (Danio rerio) feces.