Protein signatures distinctive of alpha proteobacteria and its subgroups and a model for alpha-proteobacterial evolution

S. Andrade R, El-Hani C. On the origin of mitochondria: a multilayer network approach

Backgound

The endosymbiotic theory is widely accepted to explain the origin of mitochondria from a bacterial ancestor. While ample evidence supports the intimate connection of Alphaproteobacteria to the mitochondrial ancestor, pinpointing its closest relative within sampled Alphaproteobacteria is still an open evolutionary debate. Many different phylogenetic methods and approaches have been used to answer this challenging question, further compounded by the heterogeneity of sampled taxa, varying evolutionary rates of mitochondrial proteins, and the inherent biases in each method, all factors that can produce phylogenetic artifacts. By harnessing the simplicity and interpretability of protein similarity networks, herein we re-evaluated the origin of mitochondria within an enhanced multilayer framework, which is an extension and improvement of a previously developed method.

Methods

We used a dataset of eight proteins found in mitochondria (N = 6 organisms) and bacteria (N = 80 organisms). The sequences were aligned and resulting identity matrices were combined to generate an eight-layer multiplex network. Each layer corresponded to a protein network, where nodes represented organisms and edges were placed following mutual sequence identity. The Multi-Newman-Girvan algorithm was applied to evaluate community structure, and bifurcation events linked to network partition allowed to trace patterns of divergence between studied taxa.

Results

In our network-based analysis, we first examined the topology of the 8-layer multiplex when mitochondrial sequences disconnected from the main alphaproteobacterial cluster. The resulting topology lent firm support toward an Alphaproteobacteria-sister placement for mitochondria, reinforcing the hypothesis that mitochondria diverged from the common ancestor of all Alphaproteobacteria. Additionally, we observed that the divergence of Rickettsiales was an early event in the evolutionary history of alphaproteobacterial clades.

Conclusion

By leveraging complex networks methods to the challenging question of circumscribing mitochondrial origin, we suggest that the entire Alphaproteobacteria clade is the closest relative to mitochondria (Alphaproteobacterial-sister hypothesis), echoing recent findings based on different datasets and methodologies.

Ancient Bacterial Class Alphaproteobacteria Cytochrome P450 Monooxygenases Can Be Found in Other Bacterial Species

Cytochrome P450 monooxygenases (CYPs/P450s), heme-thiolate proteins, are well-known players in the generation of chemicals valuable to humans and as a drug target against pathogens. Understanding the evolution of P450s in a bacterial population is gaining momentum. In this study, we report comprehensive analysis of P450s in the ancient group of the bacterial class Alphaproteobacteria. Genome data mining and annotation of P450s in 599 alphaproteobacterial species belonging to 164 genera revealed the presence of P450s in only 241 species belonging to 82 genera that are grouped into 143 P450 families and 214 P450 subfamilies, including 77 new P450 families. Alphaproteobacterial species have the highest average number of P450s compared to Firmicutes species and cyanobacterial species. The lowest percentage of alphaproteobacterial species P450s (2.4%) was found to be part of secondary metabolite biosynthetic gene clusters (BGCs), compared other bacterial species, indicating that during evolution large numbers of P450s became part of BGCs in other bacterial species. Our study identified that some of the P450 families found in alphaproteobacterial species were passed to other bacterial species. This is the first study to report on the identification of CYP125 P450, cholesterol and cholest-4-en-3-one hydroxylase in alphaproteobacterial species (Phenylobacterium zucineum) and to predict cholesterol side-chain oxidation capability (based on homolog proteins) by P. zucineum.

Cell cycle progression in Caulobacter requires a nucleoid-associated protein with high AT sequence recognition

Faithful cell cycle progression in the dimorphic bacterium Caulobacter crescentus requires spatiotemporal regulation of gene expression and cell pole differentiation. We discovered an essential DNA-associated protein, GapR, that is required for Caulobacter growth and asymmetric division. GapR interacts with adenine and thymine (AT)-rich chromosomal loci, associates with the promoter regions of cell cycle-regulated genes, and shares hundreds of recognition sites in common with known master regulators of cell cycle-dependent gene expression. GapR target loci are especially enriched in binding sites for the transcription factors GcrA and CtrA and overlap with nearly all of the binding sites for MucR1, a regulator that controls the establishment of swarmer cell fate. Despite constitutive synthesis, GapR accumulates preferentially in the swarmer compartment of the predivisional cell. Homologs of GapR, which are ubiquitous among the α-proteobacteria and are encoded on multiple bacteriophage genomes, also accumulate in the predivisional cell swarmer compartment when expressed in Caulobacter The Escherichia coli nucleoid-associated protein H-NS, like GapR, selectively associates with AT-rich DNA, yet it does not localize preferentially to the swarmer compartment when expressed exogenously in Caulobacter, suggesting that recognition of AT-rich DNA is not sufficient for the asymmetric accumulation of GapR. Further, GapR does not silence the expression of H-NS target genes when expressed in E. coli, suggesting that GapR and H-NS have distinct functions. We propose that Caulobacter has co-opted a nucleoid-associated protein with high AT recognition to serve as a mediator of cell cycle progression.

Discrete Responses to Limitation for Iron and Manganese in Agrobacterium tumefaciens: Influence on Attachment and Biofilm Formation

Transition metals such as iron and manganese are crucial trace nutrients for the growth of most bacteria, functioning as catalytic cofactors for many essential enzymes. Dedicated uptake and regulatory systems have evolved to ensure their acquisition for growth, while preventing toxicity. Transcriptomic analysis of the iron- and manganese-responsive regulons of Agrobacterium tumefaciens revealed that there are discrete regulatory networks that respond to changes in iron and manganese levels. Complementing earlier studies, the iron-responsive gene network is quite large and includes many aspects of iron-dependent metabolism and the iron-sparing response. In contrast, the manganese-responsive network is restricted to a limited number of genes, many of which can be linked to transport and utilization of the transition metal. Several of the target genes predicted to drive manganese uptake are required for growth under manganese-limited conditions, and an A. tumefaciens mutant with a manganese transport deficiency is attenuated for plant virulence. Iron and manganese limitation independently inhibit biofilm formation by A. tumefaciens, and several candidate genes that could impact biofilm formation were identified in each regulon. The biofilm-inhibitory effects of iron and manganese do not rely on recognized metal-responsive transcriptional regulators, suggesting alternate mechanisms influencing biofilm formation. However, under low-manganese conditions the dcpA operon is upregulated, encoding a system that controls levels of the cyclic di-GMP second messenger. Mutation of this regulatory pathway dampens the effect of manganese limitation.

IMPORTANCE

Responses to changes in transition metal levels, such as those of manganese and iron, are important for normal metabolism and growth in bacteria. Our study used global gene expression profiling to understand the response of the plant pathogen Agrobacterium tumefaciens to changes of transition metal availability. Among the properties that are affected by both iron and manganese levels are those required for normal surface attachment and biofilm formation, but the requirement for each of these transition metals is mechanistically independent from the other.

What are the Evolutionary Origins of Mitochondria? A Complex Network Approach

Mitochondria originated endosymbiotically from an Alphaproteobacteria-like ancestor. However, it is still uncertain which extant group of Alphaproteobacteria is phylogenetically closer to the mitochondrial ancestor. The proposed groups comprise the order Rickettsiales, the family Rhodospirillaceae, and the genus Rickettsia. In this study, we apply a new complex network approach to investigate the evolutionary origins of mitochondria, analyzing protein sequences modules in a critical network obtained through a critical similarity threshold between the studied sequences. The dataset included three ATP synthase subunits (4, 6, and 9) and its alphaproteobacterial homologs (b, a, and c). In all the subunits, the results gave no support to the hypothesis that Rickettsiales are closely related to the mitochondrial ancestor. Our findings support the hypothesis that mitochondria share a common ancestor with a clade containing all Alphaproteobacteria orders, except Rickettsiales.

New insights into 1-aminocyclopropane-1-carboxylate (ACC) deaminase phylogeny, evolution and ecological significance

The main objective of this work is the study of the phylogeny, evolution and ecological importance of the enzyme 1-aminocyclopropane-1-carboxylate (ACC) deaminase, the activity of which represents one of the most important and studied mechanisms used by plant growth–promoting microorganisms. The ACC deaminase gene and its regulatory elements presence in completely sequenced organisms was verified by multiple searches in diverse databases, and based on the data obtained a comprehensive analysis was conducted. Strain habitat, origin and ACC deaminase activity were taken into account when analyzing the results. In order to unveil ACC deaminase origin, evolution and relationships with other closely related pyridoxal phosphate (PLP) dependent enzymes a phylogenetic analysis was also performed. The data obtained show that ACC deaminase is mostly prevalent in some Bacteria, Fungi and members of Stramenopiles. Contrary to previous reports, we show that ACC deaminase genes are predominantly vertically inherited in various bacterial and fungal classes. Still, results suggest a considerable degree of horizontal gene transfer events, including interkingdom transfer events. A model for ACC deaminase origin and evolution is also proposed. This study also confirms the previous reports suggesting that the Lrp-like regulatory protein AcdR is a common mechanism regulating ACC deaminase expression in Proteobacteria, however, we also show that other regulatory mechanisms may be present in some Proteobacteria and other bacterial phyla. In this study we provide a more complete view of the role for ACC deaminase than was previously available. The results show that ACC deaminase may not only be related to plant growth promotion abilities, but may also play multiple roles in microorganism's developmental processes. Hence, exploring the origin and functioning of this enzyme may be the key in a variety of important agricultural and biotechnological applications.

The ζ subunit of the F1FO-ATP synthase of α-proteobacteria controls rotation of the nanomotor with a different structure

The ζ subunit is a novel natural inhibitor of the α-proteobacterial F1FO-ATPase described originally in Paracoccus denitrificans. To characterize the mechanism by which this subunit inhibits the F1FO nanomotor, the ζ subunit of Paracoccus denitrificans (Pd-ζ) was analyzed by the combination of kinetic, biochemical, bioinformatic, proteomic, and structural approaches. The ζ subunit causes full inhibition of the sulfite-activated PdF1-ATPase with an apparent IC50 of 270 nM by a mechanism independent of the ε subunit. The inhibitory region of the ζ subunit resides in the first 14 N-terminal residues of the protein, which protrude from the 4-α-helix bundle structure of the isolated ζ subunit, as resolved by NMR. Cross-linking experiments show that the ζ subunit interacts with rotor (γ) and stator (α, β) subunits of the F1-ATPase, indicating that the ζ subunit hinders rotation of the central stalk. In addition, a putatively regulatory nucleotide-binding site was found in the ζ subunit by isothermal titration calorimetry. Together, the data show that the ζ subunit controls the rotation of F1FO-ATPase by a mechanism reminiscent of, but different from, those described for mitochondrial IF1 and bacterial ε subunits where the 4-α-helix bundle of ζ seems to work as an anchoring domain that orients the N-terminal inhibitory domain to hinder rotation of the central stalk.

Whole-proteome analysis of twelve species of alphaproteobacteria links four pathogens

Thousands of whole-genome and whole-proteome sequences have been made available through advances in sequencing technology, and sequences of millions more organisms will become available in the coming years. This wealth of genetic information will provide numerous opportunities to enhance our understanding of these organisms including a greater understanding of relationships among species. Researchers have used 16S rRNA and other gene sequences to study the evolutionary origins of bacteria, but these strategies do not provide insight into the sharing of genes among bacteria via horizontal transfer. In this work we use an open source software program called pClust to cluster proteins from the complete proteomes of twelve species of Alphaproteobacteria and generate a dendrogram from the resulting orthologous protein clusters. We compare the results with dendrograms constructed using the 16S rRNA gene and multiple sequence alignment of seven housekeeping genes. Analysis of the whole proteomes of these pathogens grouped Rickettsia typhi with three other animal pathogens whereas conventional sequence analysis failed to group these pathogens together. We conclude that whole-proteome analysis can give insight into relationships among species beyond their phylogeny, perhaps reflecting the effects of horizontal gene transfer and potentially providing insight into the functions of shared genes by means of shared phenotypes.

Meta-analysis of general bacterial subclades in whole-genome phylogenies using tree topology profiling

In the last two decades, a large number of whole-genome phylogenies have been inferred to reconstruct the Tree of Life (ToL). Underlying data models range from gene or functionality content in species to phylogenetic gene family trees and multiple sequence alignments of concatenated protein sequences. Diversity in data models together with the use of different tree reconstruction techniques, disruptive biological effects and the steadily increasing number of genomes have led to a huge diversity in published phylogenies. Comparison of those and, moreover, identification of the impact of inference properties (underlying data model, inference technique) on particular reconstructions is almost impossible. In this work, we introduce tree topology profiling as a method to compare already published whole-genome phylogenies. This method requires visual determination of the particular topology in a drawn whole-genome phylogeny for a set of particular bacterial clans. For each clan, neighborhoods to other bacteria are collected into a catalogue of generalized alternative topologies. Particular topology alternatives found for an ordered list of bacterial clans reveal a topology profile that represents the analyzed phylogeny. To simulate the inhomogeneity of published gene content phylogenies we generate a set of seven phylogenies using different inference techniques and the SYSTERS-PhyloMatrix data model. After tree topology profiling on in total 54 selected published and newly inferred phylogenies, we separate artefactual from biologically meaningful phylogenies and associate particular inference results (phylogenies) with inference background (inference techniques as well as data models). Topological relationships of particular bacterial species groups are presented. With this work we introduce tree topology profiling into the scientific field of comparative phylogenomics.

Molecular characterization of NAD+-dependent DNA ligase from Wolbachia endosymbiont of lymphatic filarial parasite Brugia malayi

The lymphatic filarial parasite, Brugia malayi contains Wolbachia endobacteria that are essential for development, viability and fertility of the parasite. Therefore, wolbachial proteins have been currently seen as the potential antifilarial drug targets. NAD(+)-dependent DNA ligase is characterized as a promising drug target in several organisms due to its crucial, indispensable role in DNA replication, recombination and DNA repair. We report here the cloning, expression and purification of NAD(+)-dependent DNA ligase of Wolbachia endosymbiont of B. malayi (wBm-LigA) for its molecular characterization. wBm-LigA has all the domains that are present in nearly all the eubacterial NAD(+)-dependent DNA ligases such as N-terminal adenylation domain, OB fold, helix-hairpin-helix (HhH) and BRCT domain except zinc-binding tetracysteine domain. The purified recombinant protein (683-amino acid) was found to be biochemically active and was present in its native form as revealed by the circular dichroism and fluorescence spectra. The purified recombinant enzyme was able to catalyze intramolecular strand joining on a nicked DNA as well as intermolecular joining of the cohesive ends of BstEII restricted lamda DNA in an in vitro assay. The enzyme was localized in the various life-stages of B. malayi parasites by immunoblotting and high enzyme expression was observed in Wolbachia within B. malayi microfilariae and female adult parasites along the hypodermal chords and in the gravid portion as evident by the confocal microscopy. Ours is the first report on this enzyme of Wolbachia and these findings would assist in validating the antifilarial drug target potential of wBm-LigA in future studies.

A review of Brucella infection in marine mammals, with special emphasis on Brucella pinnipedialis in the hooded seal (Cystophora cristata)

Brucella spp. were isolated from marine mammals for the first time in 1994. Two novel species were later included in the genus; Brucella ceti and Brucella pinnipedialis, with cetaceans and seals as their preferred hosts, respectively. Brucella spp. have since been isolated from a variety of marine mammals. Pathological changes, including lesions of the reproductive organs and associated abortions, have only been registered in cetaceans. The zoonotic potential differs among the marine mammal Brucella strains. Many techniques, both classical typing and molecular microbiology, have been utilised for characterisation of the marine mammal Brucella spp. and the change from the band-based approaches to the sequence-based approaches has greatly increased our knowledge about these strains. Several clusters have been identified within the B. ceti and B. pinnipedialis species, and multiple studies have shown that the hooded seal isolates differ from other pinniped isolates. We describe how different molecular methods have contributed to species identification and differentiation of B. ceti and B. pinnipedialis, with special emphasis on the hooded seal isolates. We further discuss the potential role of B. pinnipedialis for the declining Northwest Atlantic hooded seal population.

Viral distribution and life strategies in the Bach Dang Estuary, Vietnam

Although the structure and dynamics of planktonic viruses in freshwater and seawater environments are relatively well documented, little is known about the occurrence and activity of these viruses in estuaries, especially in the tropics. Viral abundance, life strategies, and morphotype distribution were examined in the Bach Dang Estuary (Vietnam) during the dry season in 2009. The abundance of both viruses and their prokaryotic hosts decreased significantly from upstream to downstream, probably as the result of nutrient dilution and osmotic stress faced by the freshwater communities. The antibiotic mitomycin-C revealed that the fraction of lysogenic cells was substantially higher in the lower seawater part of the estuary (max 27.1%) than in the upper freshwater area where no inducible lysogens were observed. The question of whether there is a massive, continuous induction of marine lysogens caused by the mixing with freshwater is considered. Conversely, the production of lytic viruses declined as salinity increased, indicating a spatial succession of viral life strategies in this tropical estuary. Icosahedral tailless viruses with capsids smaller than 60 nm dominated the viral assemblage throughout the estuary (63.0% to 72.1% of the total viral counts), and their distribution was positively correlated with that of viral lytic production. Interestingly, the gamma-proteobacteria explained a significant portion of the variance in the <60 nm and 60 to 90 nm tailless viruses (92% and 80%, respectively), and in the Myoviridae (73%). Also, 60% of the variance of the tailless larger viruses (>90 nm) was explained by the beta-proteobacteria. Overall, these results support the view that the environment, through selection mechanisms, probably shapes the structure of the prokaryotic community. This might be in turn a source of selection for the virioplankton community via specific affiliation favoring particular morphotypes and life strategies.

Molecular signatures for the main phyla of photosynthetic bacteria and their subgroups

The bacterial groups corresponding to different photosynthetic prokaryotes are presently identified mainly on the basis of their branching in phylogenetic trees. The availability of genome sequences is enabling identification of many molecular signatures that are specific for different groups of photosynthetic bacteria. Our recent work has identified large numbers of signatures consisting of conserved inserts or deletions (indels) in widely distributed proteins, as well as whole proteins that are specific for various sequenced species/strains from Cyanobacteria, Chlorobi, and Proteobacteria phyla. Based upon these signatures, it is now possible to identify/distinguish bacteria from these phyla of photosynthetic bacteria as well as their major subclades in clear molecular terms. The use of these signatures in conjunction with phylogenomic analyses, summarized here, is leading to a holistic picture concerning the branching order and evolutionary relationships among the above groups of photosynthetic bacteria. Although detailed studies in this regard have not yet been carried on Chloroflexi and Heliobacteriaceae, we have identified some conserved indels that are specific for these groups. Some of the conserved indels for the photosynthetic bacteria are present in photosynthesis-related proteins. These include a 4 aa insert in the pyruvate flavodoxin/ferridoxin oxidoreductase that is specific for the genus Chloroflexus, a 2 aa insert in magnesium chelatase that is uniquely shared by all Cyanobacteria except the deepest branching Clade A (Gloebacterales), a 6 aa insert in an A-type flavoprotein that is specific for various marine unicellular Cyanobacteria, a 2 aa insert in heme oxygenase that is specific for various Prochlorococcus strains/isolates, and 1 aa deletion in the protein protochlorophyllide oxidoreductase that is commonly shared by various Prochlorococcus strains except the deepest branching isolates MIT 9303 and MIT 9313. The identified CSIs are located in the structures of these proteins in surface loops indicating that they may be important in mediating protein-protein interactions. The cellular functions of these conserved indels, or most of the signature proteins are presently unknown, but they provide valuable means for discovering novel properties that are unique to different groups of photosynthetic bacteria.

Genista tinctoria microsymbionts from Poland are new members of Bradyrhizobium japonicum bv. genistearum

The phylogeny and taxonomic position of slow-growing Genista tinctoria rhizobia from Poland, Ukraine and England were estimated by comparative 16S rDNA, atpD, and dnaK sequence analyses, PCR-RFLP of 16S rDNA, DNA G+C content, and DNA-DNA hybridization. Each core gene studied placed the G. tinctoria rhizobia in the genus Bradyrhizobium cluster with unequivocal bootstrap support. G. tinctoria symbionts and bradyrhizobial strains shared 96-99% similarity in 16S rDNA sequences. Their similarity for atpD and dnaK sequences was 93-99% and 89-99%, respectively. These data clearly showed that G. tinctoria rhizobia belonged to the genus Bradyrhizobium. 16S rDNA sequence analysis was in good agreement with the results of the PCR-RFLP of the 16S rRNA gene. Although the tested strains formed separate lineages to the reference bradyrhizobia their RFLP 16S rDNA patterns were quite similar. The genomic DNA G+C content of three G. tinctoria rhizobia was in the range from 60.64 to 62.83 mol%. Data for species identification were obtained from DNA-DNA hybridization experiments. G. tinctoria microsymbionts from Poland were classified within Bradyrhizobium japonicum genomospecies based on 56-82% DNA-DNA similarity.

The diversity and evolution of cell cycle regulation in alpha-proteobacteria: a comparative genomic analysis

Background

In the bacterium Caulobacter crescentus, CtrA coordinates DNA replication, cell division, and polar morphogenesis and is considered the cell cycle master regulator. CtrA activity varies during cell cycle progression and is modulated by phosphorylation, proteolysis and transcriptional control. In a phosphorylated state, CtrA binds specific DNA sequences, regulates the expression of genes involved in cell cycle progression and silences the origin of replication. Although the circuitry regulating CtrA is known in molecular detail in Caulobacter, its conservation and functionality in the other alpha-proteobacteria are still poorly understood.

Results

Orthologs of Caulobacter factors involved in the regulation of CtrA were systematically scanned in genomes of alpha-proteobacteria. In particular, orthologous genes of the divL-cckA-chpT-ctrA phosphorelay, the divJ-pleC-divK two-component system, the cpdR-rcdA-clpPX proteolysis system, the methyltransferase ccrM and transcriptional regulators dnaA and gcrA were identified in representative genomes of alpha-proteobacteria. CtrA, DnaA and GcrA binding sites and CcrM putative methylation sites were predicted in promoter regions of all these factors and functions controlled by CtrA in all alphas were predicted.

Conclusions

The regulatory cell cycle architecture was identified in all representative alpha-proteobacteria, revealing a high diversification of circuits but also a conservation of logical features. An evolutionary model was proposed where ancient alphas already possessed all modules found in Caulobacter arranged in a variety of connections. Two schemes appeared to evolve: a complex circuit in Caulobacterales and Rhizobiales and a simpler one found in Rhodobacterales.

Evolutionary placement of Xanthomonadales based on conserved protein signature sequences

Xanthomonadales comprises one of the largest phytopathogenic bacterial groups, and is currently classified within the gamma-proteobacteria. However, the phylogenetic placement of this group is not clearly resolved, and the results of different studies contradict one another. In this work, the evolutionary position of Xanthomonadales was determined by analyzing the presence of shared insertions and deletions (INDELs) in highly conserved proteins. Several distinctive insertions found in most of the members of the gamma-proteobacteria are absent in Xanthomonadales and groups such as Legionelalles, Chromatiales, Methylococcales, Thiotrichales and Cardiobacteriales. These INDELs were most likely introduced after the branching of Xanthomonadales from most of the gamma-proteobacteria and provide evidence for the phylogenetic placement of the early gamma-proteobacteria. Moreover, other proteins contain insertions exclusive to the Xanthomonadales order, confirming that this is a monophyletic group and provide important specific genetic markers. Thus, the data presented clearly support the Xanthomonadales group as an independent subdivision, and constitute one of the deepest branching lineage within the gamma-proteobacteria clade.

Evolutionary roots of iodine and thyroid hormones in cell-cell signaling

In vertebrates, thyroid hormones (THs, thyroxine, and triiodothyronine) are critical cell signaling molecules. THs regulate and coordinate physiology within and between cells, tissues, and whole organisms, in addition to controlling embryonic growth and development, via dose-dependent regulatory effects on essential genes. While invertebrates and plants do not have thyroid glands, many utilize THs for development, while others store iodine as TH derivatives or TH precursor molecules (iodotyrosines)-or produce similar hormones that act in analogous ways. Such common developmental roles for iodotyrosines across kingdoms suggest that a common endocrine signaling mechanism may account for coordinated evolutionary change in all multi-cellular organisms. Here, I expand my earlier hypothesis for the role of THs in vertebrate evolution by proposing a critical evolutionary role for iodine, the essential ingredient in all iodotyrosines and THs. Iodine is known to be crucial for life in many unicellular organisms (including evolutionarily ancient cyanobacteria), in part, because it acts as a powerful antioxidant. I propose that during the last 3-4 billion years, the ease with which various iodine species become volatile, react with simple organic compounds, and catalyze biochemical reactions explains why iodine became an essential constituent of life and the Earth's atmosphere-and a potential marker for the origins of life. From an initial role as membrane antioxidant and biochemical catalyst, spontaneous coupling of iodine with tyrosine appears to have created a versatile, highly reactive and mobile molecule, which over time became integrated into the machinery of energy production, gene function, and DNA replication in mitochondria. Iodotyrosines later coupled together to form THs, the ubiquitous cell-signaling molecules used by all vertebrates. Thus, due to their evolutionary history, THs, and their derivative and precursors molecules not only became essential for communicating within and between cells, tissues and organs, and for coordinating development and whole-body physiology in vertebrates, but they can also be shared between organisms from different kingdoms.

The differential interaction of Brucella and ochrobactrum with innate immunity reveals traits related to the evolution of stealthy pathogens

Background

During evolution, innate immunity has been tuned to recognize pathogen-associated molecular patterns. However, some α-Proteobacteria are stealthy intracellular pathogens not readily detected by this system. Brucella members follow this strategy and are highly virulent, but other Brucellaceae like Ochrobactrum are rhizosphere inhabitants and only opportunistic pathogens. To gain insight into the emergence of the stealthy strategy, we compared these two phylogenetically close but biologically divergent bacteria.

Methodology/Principal Findings

In contrast to Brucella abortus, Ochrobactrum anthropi did not replicate within professional and non-professional phagocytes and, whereas neutrophils had a limited action on B. abortus, they were essential to control O. anthropi infections. O. anthropi triggered proinflammatory responses markedly lower than Salmonella enterica but higher than B. abortus. In macrophages and dendritic cells, the corresponding lipopolysaccharides reproduced these grades of activation, and binding of O. anthropi lipopolysaccharide to the TLR4 co-receptor MD-2 and NF-κB induction laid between those of B. abortus and enteric bacteria lipopolysaccharides. These differences correlate with reported variations in lipopolysaccharide core sugars, sensitivity to bactericidal peptides and outer membrane permeability.

Conclusions/Significance

The results suggest that Brucellaceae ancestors carried molecules not readily recognized by innate immunity, so that non-drastic variations led to the emergence of stealthy intracellular parasites. They also suggest that some critical envelope properties, like selective permeability, are profoundly altered upon modification of pathogen-associated molecular patterns, and that this represents a further adaptation to the host. It is proposed that this adaptive trend is relevant in other intracellular α-Proteobacteria like Bartonella, Rickettsia, Anaplasma, Ehrlichia and Wolbachia.

dnaJ is a useful phylogenetic marker for alphaproteobacteria

In the past, bacterial phylogeny relied almost exclusively on 16S rRNA gene sequence analysis. More recently, multilocus sequence analysis has been used to infer organismal phylogenies. In this study, the dnaJ chaperone gene was investigated as a marker for phylogeny studies in alphaproteobacteria. Preliminary analysis of G+C contents and G+C3s contents (the G+C content of the synonymous third codon position) showed no clear evidence of horizontal transfer of this gene in proteobacteria. dnaJ-based phylogenies were then analysed at three taxonomic levels: the Proteobacteria, the Alphaproteobacteria and the genus Mesorhizobium. Dendrograms based on DnaJ and 16S rRNA gene sequences revealed the same topology described previously for the Proteobacteria. These results indicate that the DnaJ phylogenetic signal is able to reproduce the accepted relationships among the five classes of the Proteobacteria. At a lower taxonomic level, using 20 alphaproteobacteria, the 16S rRNA gene-based phylogeny is distinct from the one based on DnaJ sequence analysis. Although the same clusters are generated, only the topology of the DnaJ tree is consistent with broader phylogenies from recent studies based on concatenated alignments of multiple core genes. For example, the DnaJ tree shows the two clusters within the Rhizobiales as closely related, as expected, while the 16S rRNA gene-based phylogeny shows them as distantly related. In order to evaluate the phylogenetic performance of dnaJ at the genus level, a multilocus analysis based on five housekeeping genes (atpD, gapA, gyrB, recA and rplB) was performed for ten Mesorhizobium species. In contrast to the 16S rRNA gene, the DnaJ sequence analysis generated a tree similar to the multilocus dendrogram. For identification of chickpea mesorhizobium isolates, a dnaJ nucleotide sequence-based tree was used. Despite different topologies, 16S rRNA gene- and dnaJ-based trees led to the same species identification. This study suggests that the dnaJ gene is a good phylogenetic marker, particularly for the class Alphaproteobacteria, since its phylogeny is consistent with phylogenies based on multilocus approaches.

Phylogenomics and signature proteins for the alpha proteobacteria and its main groups

Background

Alpha proteobacteria are one of the largest and most extensively studied groups within bacteria. However, for these bacteria as a whole and for all of its major subgroups (viz. Rhizobiales, Rhodobacterales, Rhodospirillales, Rickettsiales, Sphingomonadales and Caulobacterales), very few or no distinctive molecular or biochemical characteristics are known.

Results

We have carried out comprehensive phylogenomic analyses by means of Blastp and PSI-Blast searches on the open reading frames in the genomes of several α-proteobacteria (viz. Bradyrhizobium japonicum, Brucella suis, Caulobacter crescentus, Gluconobacter oxydans, Mesorhizobium loti, Nitrobacter winogradskyi, Novosphingobium aromaticivorans, Rhodobacter sphaeroides 2.4.1, Silicibacter sp. TM1040, Rhodospirillum rubrum and Wolbachia (Drosophila) endosymbiont). These studies have identified several proteins that are distinctive characteristics of all α-proteobacteria, as well as numerous proteins that are unique repertoires of all of its main orders (viz. Rhizobiales, Rhodobacterales, Rhodospirillales, Rickettsiales, Sphingomonadales and Caulobacterales) and many families (viz. Rickettsiaceae, Anaplasmataceae, Rhodospirillaceae, Acetobacteraceae, Bradyrhiozobiaceae, Brucellaceae and Bartonellaceae). Many other proteins that are present at different phylogenetic depths in α-proteobacteria provide important information regarding their evolution. The evolutionary relationships among α-proteobacteria as deduced from these studies are in excellent agreement with their branching pattern in the phylogenetic trees and character compatibility cliques based on concatenated sequences for many conserved proteins. These studies provide evidence that the major groups within α-proteobacteria have diverged in the following order: (Rickettsiales(Rhodospirillales (Sphingomonadales (Rhodobacterales (Caulobacterales-Parvularculales (Rhizobiales)))))). We also describe two conserved inserts in DNA Gyrase B and RNA polymerase beta subunit that are distinctive characteristics of the Sphingomonadales and Rhodosprilllales species, respectively. The results presented here also provide support for the grouping of Hyphomonadaceae and Parvularcula species with the Caulobacterales and the placement of Stappia aggregata with the Rhizobiaceae group.

Conclusion

The α-proteobacteria-specific proteins and indels described here provide novel and powerful means for the taxonomic, biochemical and molecular biological studies on these bacteria. Their functional studies should prove helpful in identifying novel biochemical and physiological characteristics that are unique to these bacteria.

A robust species tree for the alphaproteobacteria

The branching order and coherence of the alphaproteobacterial orders have not been well established, and not all studies have agreed that mitochondria arose from within the Rickettsiales. A species tree for 72 alphaproteobacteria was produced from a concatenation of alignments for 104 well-behaved protein families. Coherence was upheld for four of the five orders with current standing that were represented here by more than one species. However, the family Hyphomonadaceae was split from the other Rhodobacterales, forming an expanded group with Caulobacterales that also included Parvularcula. The three earliest-branching alphaproteobacterial orders were the Rickettsiales, followed by the Rhodospirillales and then the Sphingomonadales. The principal uncertainty is whether the expanded Caulobacterales group is more closely associated with the Rhodobacterales or the Rhizobiales. The mitochondrial branch was placed within the Rickettsiales as a sister to the combined Anaplasmataceae and Rickettsiaceae, all subtended by the Pelagibacter branch. Pelagibacter genes will serve as useful additions to the bacterial outgroup in future evolutionary studies of mitochondrial genes, including those that have transferred to the eukaryotic nucleus.

Genome analysis of DNA repair genes in the alpha proteobacterium Caulobacter crescentus

Background

The integrity of DNA molecules is fundamental for maintaining life. The DNA repair proteins protect organisms against genetic damage, by removal of DNA lesions or helping to tolerate them. DNA repair genes are best known from the gamma-proteobacterium Escherichia coli, which is the most understood bacterial model. However, genome sequencing raises questions regarding uniformity and ubiquity of these DNA repair genes and pathways, reinforcing the need for identifying genes and proteins, which may respond to DNA damage in other bacteria.

Results

In this study, we employed a bioinformatic approach, to analyse and describe the open reading frames potentially related to DNA repair from the genome of the alpha-proteobacterium Caulobacter crescentus. This was performed by comparison with known DNA repair related genes found in public databases. As expected, although C. crescentus and E. coli bacteria belong to separate phylogenetic groups, many of their DNA repair genes are very similar. However, some important DNA repair genes are absent in the C. crescentus genome and other interesting functionally related gene duplications are present, which do not occur in E. coli. These include DNA ligases, exonuclease III (xthA), endonuclease III (nth), O₆-methylguanine-DNA methyltransferase (ada gene), photolyase-like genes, and uracil-DNA-glycosylases. On the other hand, the genes imuA and imuB, which are involved in DNA damage induced mutagenesis, have recently been described in C. crescentus, but are absent in E. coli. Particularly interesting are the potential atypical phylogeny of one of the photolyase genes in alpha-proteobacteria, indicating an origin by horizontal transfer, and the duplication of the Ada orthologs, which have diverse structural configurations, including one that is still unique for C. crescentus.

Conclusion

The absence and the presence of certain genes are discussed and predictions are made considering the particular aspects of the C. crescentus among other known DNA repair pathways. The observed differences enlarge what is known for DNA repair in the Bacterial world, and provide a useful framework for further experimental studies in this organism.

A new alpha-proteobacterial clade of Bdellovibrio-like predators: implications for the mitochondrial endosymbiotic theory

Bdellovibrio-and-like organisms (BALOs) are peculiar, ubiquitous, small-sized, highly motile Gram-negative bacteria that are obligatory predators of other bacteria. Typically, these predators invade the periplasm of their prey where they grow and replicate. To date, BALOs constitute two highly diverse families affiliated with the delta-proteobacteria class. In this study, Micavibrio spp., a BALO lineage of epibiotic predators, were isolated from soil. These bacteria attach to digest and grow at the expense of other prokaryotes, much like other BALOs. Multiple phylogenetic analyses based on six genes revealed that they formed a deep branch within the alpha-proteobacteria, not affiliated with any of the alpha-proteobacterial orders. The presence of BALOs deep among the alpha-proteobacteria suggests that their peculiar mode of parasitism maybe an ancestral character in this proteobacterial class. The origin of the mitochondrion from an alpha-proteobacterium endosymbiont is strongly supported by molecular phylogenies. Accumulating data suggest that the endosymbiont's host was also a prokaryote. As prokaryotes are unable to phagocytose, the means by which the endosymbiont gained access into its host remains mysterious. We here propose a scenario based on the BALO feeding-mode to hypothesize a mechanism at play at the origin of the mitochondrial endosymbiosis.

Application of the character compatibility approach to generalized molecular sequence data: branching order of the proteobacterial subdivisions

The character compatibility approach, which removes all homoplasic characters and involves finding the largest clique of compatible characters in a dataset, in principle, provides a powerful means for obtaining correct topology in difficult to resolve cases. However, the usefulness of this approach to generalized molecular sequence data for phylogeny determination has not been studied in the past. We have used this approach to determine the topology of 23 proteobacterial species (6 each of alpha-, beta- and gamma-, 3 delta-, and 2 epsilon-proteobacteria) using sequence data for 10 conserved proteins (Hsp60, Hsp70, EF-Tu, EF-G, alanyl-tRNA synthetase, RecA, GyrA, GyrB, RpoB and RpoC). All sites in the sequence alignments of these proteins where only two amino acids were found, with each amino acid present in at least two species, were selected. Mutual compatibility determination on these binary state sites was carried out by two means. In one case, all of these sites were combined into a large dataset (Set A; 957 characters) prior to compatibility analysis. In the second case, compatibility analysis was carried out on characters from individual proteins and all compatible sites were combined into a large dataset (Set B; 398 characters) for further studies. Upon compatibility analyses, the largest cliques that were obtained from Sets A and B consisted of 337 and 323 compatible characters, respectively. In these cliques, all proteobacterial subgroups were clearly distinguished and branching orders of most of the species were also resolved. The epsilon-proteobacteria exhibited the earliest branching, whereas the beta- and gamma-subgroups were found to have emerged last. The relative placement of the alpha- and delta-subgroups, however, was not resolved. The topology of these species was also determined based on 16S rRNA sequences and a concatenated dataset of sequences for all 10 proteins by means of neighbor-joining, maximum likelihood, and maximum parsimony methods. In the protein trees, all proteobacterial groups were reliably resolved and they branched in the following order: (epsilon(delta(alpha(beta,gamma)))). However, in the rRNA trees, the gamma- and beta-subgroups exhibited polyphyletic branching and many internal nodes were not resolved. These results indicate that the character compatibility analysis using generalized molecular sequence data provides a powerful means for evolutionary studies. Based on molecular sequences, it should be possible to obtain very large datasets of compatible characters that should prove very helpful in clarifying difficult to resolve phylogenetic relationships.

Computational reconstruction of iron- and manganese-responsive transcriptional networks in alpha-proteobacteria

We used comparative genomics to investigate the distribution of conserved DNA-binding motifs in the regulatory regions of genes involved in iron and manganese homeostasis in alpha-proteobacteria. Combined with other computational approaches, this allowed us to reconstruct the metal regulatory network in more than three dozen species with available genome sequences. We identified several classes of cis-acting regulatory DNA motifs (Irr-boxes or ICEs, RirA-boxes, Iron-Rhodo-boxes, Fur-alpha-boxes, Mur-box or MRS, MntR-box, and IscR-boxes) in regulatory regions of various genes involved in iron and manganese uptake, Fe-S and heme biosynthesis, iron storage, and usage. Despite the different nature of the iron regulons in selected lineages of alpha-proteobacteria, the overall regulatory network is consistent with, and confirmed by, many experimental observations. This study expands the range of genes involved in iron homeostasis and demonstrates considerable interconnection between iron-responsive regulatory systems. The detailed comparative and phylogenetic analyses of the regulatory systems allowed us to propose a theory about the possible evolution of Fe and Mn regulons in alpha-proteobacteria. The main evolutionary event likely occurred in the common ancestor of the Rhizobiales and Rhodobacterales, where the Fur protein switched to regulating manganese transporters (and hence Fur had become Mur). In these lineages, the role of global iron homeostasis was taken by RirA and Irr, two transcriptional regulators that act by sensing the physiological consequence of the metal availability rather than its concentration per se, and thus provide for more flexible regulation.

The availability of hundreds of complete genomes allows one to use comparative genomics to describe key metabolic processes and regulatory gene networks. Genome context analyses and comparisons of transcription factor binding sites between genomes offer a powerful approach for functional gene annotation. Reconstruction of transcriptional regulatory networks allows for better understanding of cellular processes, which can be substantiated by direct experimentation. Iron homeostasis in bacteria is conferred by the regulation of various iron uptake transporters, iron storage ferritins, and iron-containing enzymes. In high concentrations, iron is poisonous for the cell, so strict control of iron homeostasis is maintained, mostly at the level of transcription by iron-responsive regulators. Despite their general importance, iron regulatory networks in most bacterial species are not well-understood. In this study, Rodionov and colleagues applied comparative genomic approaches to describe the regulatory network formed by genes involved in iron homeostasis in the alpha subclass of proteobacteria, which have extremely versatile lifestyles. These networks are mediated by a set of various DNA motifs (or regulatory signals) that occur in 5′ gene regions and involve at least six different metal-responsive regulators. This study once again shows the power of comparative genomics in the analysis of complex regulatory networks and their evolution.

Signature proteins that are distinctive characteristics of Actinobacteria and their subgroups

The Actinobacteria constitute one of the main phyla of Bacteria. Presently, no morphological and very few molecular characteristics are known which can distinguish species of this highly diverse group. In this work, we have analyzed the genomes of four actinobacteria (viz. Mycobacterium leprae TN, Leifsonia xyli subsp. xyli str. CTCB07, Bifidobacterium longum NCC2705 and Thermobifida fusca YX) to search for proteins that are unique to Actinobacteria. Our analyses have identified 233 actinobacteria-specific proteins, homologues of which are generally not present in any other bacteria. These proteins can be grouped as follows: (i) 29 proteins uniquely present in most sequenced actinobacterial genomes; (ii) 6 proteins present in almost all actinobacteria except Bifidobacterium longum and another 37 proteins absent in B. longum and few other species; (iii) 11 proteins which are mainly present in Corynebacterium, Mycobacterium and Nocardia (CMN) subgroup as well as Streptomyces, T. fusca and Frankia sp., but they are not found in Bifidobacterium and Micrococcineae; (iv) 8 proteins that are specific for T. fusca and Streptomyces species, plus 2 proteins also present in the Frankia species; (v) 13 proteins that are specific for the Corynebacterineae or the CMN group; (vi) 14 proteins only found in Mycobacterium and Nocardia; (vii) 24 proteins unique to different Mycobacterium species; (viii) 8 proteins specific to the Micrococcineae; (ix) 85 proteins which are distributed sporadically in actinobacterial species. Additionally, many examples of lateral gene transfer from Actinobacteria to Magnetospirillum magnetotacticum have also been identified. The identified proteins provide novel molecular means for defining and circumscribing the Actinobacteria phylum and a number of subgroups within it. The distribution of these proteins also provides useful information regarding interrelationships among the actinobacterial subgroups. Most of these proteins are of unknown function and studies aimed at understanding their cellular functions should reveal common biochemical and physiological characteristics unique to either all actinobacteria or particular subgroups of them. The identified proteins also provide potential targets for development of drugs that are specific for actinobacteria.

BLAST screening of chlamydial genomes to identify signature proteins that are unique for the Chlamydiales, Chlamydiaceae, Chlamydophila and Chlamydia groups of species

Background

Chlamydiae species are of much importance from a clinical viewpoint. Their diversity both in terms of their numbers as well as clinical involvement are presently believed to be significantly underestimated. The obligate intracellular nature of chlamydiae has also limited their genetic and biochemical studies. Thus, it is of importance to develop additional means for their identification and characterization.

Results

We have carried out analyses of available chlamydiae genomes to identify sets of unique proteins that are either specific for all Chlamydiales genomes, or different Chlamydiaceae family members, or members of the Chlamydia and Chlamydophila genera, or those unique to Protochlamydia amoebophila, but which are not found in any other bacteria. In total, 59 Chlamydiales-specific proteins, 79 Chlamydiaceae-specific proteins, 20 proteins each that are specific for both Chlamydia and Chlamydophila and 445 ORFs that are Protochlamydia-specific were identified. Additionally, 33 cases of possible gene loss or lateral gene transfer were also detected.

Conclusion

The identified chlamydiae-lineage specific proteins, many of which are highly conserved, provide novel biomarkers that should prove of much value in the diagnosis of these bacteria and in exploration of their prevalence and diversity. These conserved protein sequences (CPSs) also provide novel therapeutic targets for drugs that are specific for these bacteria. Lastly, functional studies on these chlamydiae or chlamydiae subgroup-specific proteins should lead to important insights into lineage-specific adaptations with regards to development, infectivity and pathogenicity.

Signature proteins that are distinctive of alpha proteobacteria

Background

The 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 α-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.

Results

Blast analyses of protein sequences from the above genomes have led to the identification of 61 proteins which are distinctive characteristics of α-proteobacteria and are generally not found in any other bacteria. These α-proteobacterial signature proteins are generally of hypothetical functions and they can be classified as follows: (i) Six proteins (CC2102, CC3292, CC3319, CC1887, CC1725 and CC1365) which are uniquely present in most sequenced α-proteobacterial genomes; (ii) Ten proteins (CC1211, CC1886, CC2245, CC3470, CC0520, CC0365, CC0366, CC1977, CC3010 and CC0100) which are present in all α-proteobacteria except the Rickettsiales; (iii) Five proteins (CC2345, CC3115, CC3401, CC3467 and CC1021) not found in the intracellular bacteria belonging to the order Rickettsiales and the Bartonellaceae family; (iv) Four proteins (CC1652, CC2247, CC3295 and CC1035) that are absent from various Rickettsiales as well as Rhodobacterales; (v) Three proteins (RP104, RP105 and RP106) that are unique to the order Rickettsiales and four proteins (RP766, RP192, RP030 and RP187) which are specific for the Rickettsiaceae family; (vi) Six proteins (BQ00140, BQ00720, BQ03880, BQ12030, BQ07670 and BQ11900) which are specific to the order Rhizobiales; (vii) Four proteins (BQ01660, BQ02450, BQ03770 and BQ13470) which are specific for the order Rhizobiales excluding the family Bradyrhizobiaceae; (viii) Nine proteins (BQ12190, BQ11460, BQ11450, BQ11430, BQ11380, BQ11160, BQ11120, BQ11100 and BQ11030 which are distinctive of the Bartonellaceae family;(ix) Six proteins (CC0189, CC0569, CC0331, CC0349, CC2323 and CC2637) which show sporadic distribution in α-proteobacteria, (x) Four proteins (CC2585, CC0226, CC2790 and RP382) in which lateral gene transfers are indicated to have occurred between α-proteobacteria and a limited number of other bacteria.

Conclusion

The identified proteins provide novel means for defining and identifying the α-proteobacteria and many of its subgroups in clear molecular terms and in understanding the evolution of this group of species. These signature proteins, together with the large number of α-proteobacteria specific indels that have recently been identified , provide evidence that all species from this diverse group share many unifying and distinctive characteristics. Functional studies on these proteins should prove very helpful in the identification of such characteristics.