Construction of a phylogenetic tree of photosynthetic prokaryotes based on average similarities of whole genome sequences

Comparative genomic analysis of Chryseobacterium species: deep insights into plant-growth-promoting and halotolerant capacities

Members of the genus Chryseobacterium have attracted great interest as beneficial bacteria that can promote plant growth and biocontrol. Given the recent risks of climate change, it is important to develop tolerance strategies for efficient applications of plant-beneficial bacteria in saline environments. However, the genetic determinants of plant-growth-promoting and halotolerance effects in Chryseobacterium have not yet been investigated at the genomic level. Here, a comparative genomic analysis was conducted with seven Chryseobacterium species. Phylogenetic and phylogenomic analyses revealed niche-specific evolutionary distances between soil and freshwater Chryseobacterium species, consistent with differences in genomic statistics, indicating that the freshwater bacteria have smaller genome sizes and fewer genes than the soil bacteria. Phosphorus- and zinc-cycling genes (required for nutrient acquisition in plants) were universally present in all species, whereas nitrification and sulphite reduction genes (required for nitrogen- and sulphur-cycling, respectively) were distributed only in soil bacteria. A pan-genome containing 6842 gene clusters was constructed, which reflected the general features of the core, accessory and unique genomes. Halotolerant species with an accessory genome shared a Kdp potassium transporter and biosynthetic pathways for branched-chain amino acids and the carotenoid lycopene, which are associated with countermeasures against salt stress. Protein-protein interaction network analysis was used to define the genetic determinants of Chryseobacterium salivictor NBC122 that reduce salt damage in bacteria and plants. Sixteen hub genes comprised the aromatic compound degradation and Por secretion systems, which are required to cope with complex stresses associated with saline environments. Horizontal gene transfer and CRISPR-Cas analyses indicated that C. salivictor NBC122 underwent more evolutionary events when interacting with different environments. These findings provide deep insights into genomic adaptation to dynamic interactions between plant-growth-promoting Chryseobacterium and salt stress.

Phylogeny analysis of whole protein-coding genes in metagenomic data detected an environmental gradient for the microbiota

Environmental factors affect the growth of microorganisms and therefore alter the composition of microbiota. Correlative analysis of the relationship between metagenomic composition and the environmental gradient can help elucidate key environmental factors and establishment principles for microbial communities. However, a reasonable method to quantitatively compare whole metagenomic data and identify the primary environmental factors for the establishment of microbiota has not been reported so far. In this study, we developed a method to compare whole proteomes deduced from metagenomic shotgun sequencing data, and quantitatively display their phylogenetic relationships as metagenomic trees. We called this method Metagenomic Phylogeny by Average Sequence Similarity (MPASS). We also compared one of the metagenomic trees with dendrograms of environmental factors using a comparison tool for phylogenetic trees. The MPASS method correctly constructed metagenomic trees of simulated metagenomes and soil and water samples. The topology of the metagenomic tree of samples from the Kirishima hot springs area in Japan was highly similarity to that of the dendrograms based on previously reported environmental factors for this area. The topology of the metagenomic tree also reflected the dynamics of microbiota at the taxonomic and functional levels. Our results strongly suggest that MPASS can successfully classify metagenomic shotgun sequencing data based on the similarity of whole protein-coding sequences, and will be useful for the identification of principal environmental factors for the establishment of microbial communities. Custom Perl script for the MPASS pipeline is available at https://github.com/s0sat/MPASS.

Heterologous complementation systems verify the mosaic distribution of three distinct protoporphyrinogen IX oxidase in the cyanobacterial phylum

The pathways for synthesizing tetrapyrroles, including heme and chlorophyll, are well-conserved among organisms, despite the divergence of several enzymes in these pathways. Protoporphyrinogen IX oxidase (PPOX), which catalyzes the last common step of the heme and chlorophyll biosynthesis pathways, is encoded by three phylogenetically-unrelated genes, hemY, hemG and hemJ. All three types of homologues are present in the cyanobacterial phylum, showing a mosaic phylogenetic distribution. Moreover, a few cyanobacteria appear to contain two types of PPOX homologues. Among the three types of cyanobacterial PPOX homologues, only a hemJ homologue has been experimentally verified for its functionality. An objective of this study is to provide experimental evidence for the functionality of the cyanobacterial PPOX homologues by using two heterologous complementation systems. First, we introduced hemY and hemJ homologues from Gloeobacter violaceus PCC7421, hemY homologue from Trichodesmium erythraeum, and hemG homologue from Prochlorococcus marinus MIT9515 into a ΔhemG strain of E. coli. hemY homologues from G. violaceus and T. erythraeum, and the hemG homologue of P. marinus complimented the E. coli strain. Subsequently, we attempted to replace the endogenous hemJ gene of the cyanobacterium Synechocystis sp. PCC6803 with the four PPOX homologues mentioned above. Except for hemG from P. marinus, the other PPOX homologues substituted the function of hemJ in Synechocystis. These results show that all four homologues encode functional PPOX. The transformation of Synechocystis with G. violaceus hemY homologue rendered the cells sensitive to an inhibitor of the HemY-type PPOX, acifluorfen, indicating that the hemY homologue is sensitive to this inhibitor, while the wild-type G. violaceus was tolerant to it, most likely due to the presence of HemJ protein. These results provide an additional level of evidence that G. violaceus contains two types of functional PPOX.

Photosynthetic Systems Suggest an Evolutionary Pathway to Diderms

Bacteria are divided primarily into monoderms (with one cell membrane, and usually Gram-positive, due to a thick peptidoglycan layer) and diderms (with two cell membranes, and mostly Gram-negative, due to a thin peptidoglycan layer sandwiched between the two membranes). Photosynthetic species are spread among the taxonomic groups, some having type I reaction centers (RCI in monoderm phylum Firmicutes; and diderm phyla Acidobacteria and Chlorobi), others with type II reaction centers (RCII in monoderm phylum Chloroflexi; and diderm taxa Gemmatimonadetes, and alpha-, beta-, and gamma-Proteobacteria), and some containing both (RCI and RCII, only in diderm phylum Cyanobacteria). In most bacterial phylograms, photosystem types and diderm taxa are polyphyletic. A more parsimonious arrangement, which is supported by photosystem evolution, as well as additional sets of molecular characters, suggests that endosymbiotic events resulted in the formation of the diderms. In the model presented, monoderms readily form a monophyletic group, while diderms are produced by at least two endosymbiotic events, followed by additional evolutionary changes.

Exploring the Diversity Within the Genus Francisella - An Integrated Pan-Genome and Genome-Mining Approach

Pan-genome analysis is a powerful method to explore genomic heterogeneity and diversity of bacterial species. Here we present a pan-genome analysis of the genus Francisella, comprising a dataset of 63 genomes and encompassing clinical as well as environmental isolates from distinct geographic locations. To determine the evolutionary relationship within the genus, we performed phylogenetic whole-genome studies utilizing the average nucleotide identity, average amino acid identity, core genes and non-recombinant loci markers. Based on the analyses, the phylogenetic trees obtained identified two distinct clades, A and B and a diverse cluster designated C. The sizes of the pan-, core-, cloud-, and shell-genomes of Francisella were estimated and compared to those of two other facultative intracellular pathogens, Legionella and Piscirickettsia. Francisella had the smallest core-genome, 692 genes, compared to 886 and 1,732 genes for Legionella and Piscirickettsia respectively, while the pan-genome of Legionella was more than twice the size of that of the other two genera. Also, the composition of the Francisella Type VI secretion system (T6SS) was analyzed. Distinct differences in the gene content of the T6SS were identified. In silico approaches performed to identify putative substrates of these systems revealed potential effectors targeting the cell wall, inner membrane, cellular nucleic acids as well as proteins, thus constituting attractive targets for site-directed mutagenesis. The comparative analysis performed here provides a comprehensive basis for the assessment of the phylogenomic relationship of members of the genus Francisella and for the identification of putative T6SS virulence traits.

In vivo Proinflammatory Cytokine Production by CD-1 Mice in Response to Equipotential Doses of Rhodobacter capsulatus PG and Salmonella enterica Lipopolysaccharides

The capacities of relatively nontoxic lipopolysaccharide (LPS) from Rhodobacter capsulatus PG and highly potent LPS from Salmonella enterica serovar Typhimurium to evoke proinflammatory cytokine production have been compared in vivo. Intravenous administration of S. enterica LPS at a relatively low dose (1 mg/kg body weight) led to upregulation of TNF-α, IL-6, and IFN-γ production by non-sensitized CD-1 mice. LPS from R. capsulatus PG used at a four-times higher dose than that from S. enterica elicited production of almost the same amount of systemic TNF-α; therefore, the doses of 4 mg/kg LPS from R. capsulatus PG and 1 mg/kg LPS from S. enterica were considered to be approximately equipotential doses with respect to the LPS-dependent TNF-α production by CD-1 mice. Rhodobacter capsulatus PG LPS was a weaker inducer of the production of TNF-α, IL-6, and IFN-γ, as compared to the equipotential dose of S. enterica LPS. Administration of R. capsulatus PG LPS before S. enterica LPS decreased production of IFN-γ, but not of TNF-α and IL-6, induced by S. enterica LPS. Rhodobacter capsulatus PG LPS also suppressed IFN-γ production induced by S. enterica LPS when R. capsulatus PG LPS had been injected as little as 10 min after S. enterica LPS, but to a much lesser extent. Rhodobacter capsulatus PG LPS did not affect TNF-α and IL-6 production induced by the equipotential dose of S. enterica LPS. In order to draw conclusion on the endotoxic activity of particular LPSs, species-specific structure or arrangement of the animal or human immune systems should be considered.

Comparative analyses of whole-genome protein sequences from multiple organisms

Phylogenies based on entire genomes are a powerful tool for reconstructing the Tree of Life. Several methods have been proposed, most of which employ an alignment-free strategy. Average sequence similarity methods are different than most other whole-genome methods, because they are based on local alignments. However, previous average similarity methods fail to reconstruct a correct phylogeny when compared against other whole-genome trees. In this study, we developed a novel average sequence similarity method. Our method correctly reconstructs the phylogenetic tree of in silico evolved E. coli proteomes. We applied the method to reconstruct a whole-proteome phylogeny of 1,087 species from all three domains of life, Bacteria, Archaea, and Eucarya. Our tree was automatically reconstructed without any human decisions, such as the selection of organisms. The tree exhibits a concentric circle-like structure, indicating that all the organisms have similar total branch lengths from their common ancestor. Branching patterns of the members of each phylum of Bacteria and Archaea are largely consistent with previous reports. The topologies are largely consistent with those reconstructed by other methods. These results strongly suggest that this approach has sufficient taxonomic resolution and reliability to infer phylogeny, from phylum to strain, of a wide range of organisms.

Nutrient Acquisition and the Metabolic Potential of Photoferrotrophic Chlorobi

Anoxygenic photosynthesis evolved prior to oxygenic photosynthesis and harnessed energy from sunlight to support biomass production on the early Earth. Models that consider the availability of electron donors predict that anoxygenic photosynthesis using Fe(II), known as photoferrotrophy, would have supported most global primary production before the proliferation of oxygenic phototrophs at approximately 2.3 billion years ago. These photoferrotrophs have also been implicated in the deposition of banded iron formations, the world's largest sedimentary iron ore deposits that formed mostly in late Archean and early Proterozoic Eons. In this work we present new data and analyses that illuminate the metabolic capacity of photoferrotrophy in the phylum Chlorobi. Our laboratory growth experiments and biochemical analyses demonstrate that photoferrotrophic Chlorobi are capable of assimilatory sulfate reduction and nitrogen fixation under sulfate and nitrogen limiting conditions, respectively. Furthermore, the evolutionary histories of key enzymes in both sulfur (CysH and CysD) and nitrogen fixation (NifDKH) pathways are convoluted; protein phylogenies, however, suggest that early Chlorobi could have had the capacity to assimilate sulfur and fix nitrogen. We argue, then, that the capacity for photoferrotrophic Chlorobi to acquire these key nutrients enabled them to support primary production and underpin global biogeochemical cycles in the Precambrian.

Diversity in photosynthetic electron transport under [CO2]-limitation: the cyanobacterium Synechococcus sp. PCC 7002 and green alga Chlamydomonas reinhardtii drive an O2-dependent alternative electron flow and non-photochemical quenching of chlorophyll fluorescence during CO2-limited photosynthesis

Some cyanobacteria, but not all, experience an induction of alternative electron flow (AEF) during CO₂-limited photosynthesis. For example, Synechocystis sp. PCC 6803 (S. 6803) exhibits AEF, but Synechococcus elongatus sp. PCC 7942 does not. This difference is due to the presence of flavodiiron 2 and 4 proteins (FLV2/4) in S. 6803, which catalyze electron donation to O₂. In this study, we observed a low-[CO₂] induced AEF in the marine cyanobacterium Synechococcus sp. PCC 7002 that lacks FLV2/4. The AEF shows high affinity for O₂, compared with AEF mediated by FLV2/4 in S. 6803, and can proceed under extreme low [O₂] (about a few µM O₂). Further, the transition from CO₂-saturated to CO₂-limited photosynthesis leads a preferential excitation of PSI to PSII and increased non-photochemical quenching of chlorophyll fluorescence. We found that the model green alga Chlamydomonas reinhardtii also has an O₂-dependent AEF showing the same affinity for O₂ as that in S. 7002. These data represent the diverse molecular mechanisms to drive AEF in cyanobacteria and green algae. In this paper, we further discuss the diversity, the evolution, and the physiological function of strategy to CO₂-limitation in cyanobacterial and green algal photosynthesis.

Origin of Bacteriochlorophyll a and the Early Diversification of Photosynthesis

Photosynthesis originated in the domain Bacteria billions of years ago; however, the identity of the last common ancestor to all phototrophic bacteria remains undetermined and speculative. Here I present the evolution of BchF or 3-vinyl-bacteriochlorophyll hydratase, an enzyme exclusively found in bacteria capable of synthetizing bacteriochlorophyll a. I show that BchF exists in two forms originating from an early divergence, one found in the phylum Chlorobi, including its paralogue BchV, and a second form that was ancestral to the enzyme found in the remaining anoxygenic phototrophic bacteria. The phylogeny of BchF is consistent with bacteriochlorophyll a evolving in an ancestral phototrophic bacterium that lived before the radiation event that gave rise to the phylum Chloroflexi, Chlorobi, Acidobacteria, Proteobacteria, and Gemmatimonadetes, but only after the divergence of Type I and Type II reaction centers. Consequently, it is suggested that the lack of phototrophy in many groups of extant bacteria is a derived trait.

A fresh look at the evolution and diversification of photochemical reaction centers

In this review, I reexamine the origin and diversification of photochemical reaction centers based on the known phylogenetic relations of the core subunits, and with the aid of sequence and structural alignments. I show, for example, that the protein folds at the C-terminus of the D1 and D2 subunits of Photosystem II, which are essential for the coordination of the water-oxidizing complex, were already in place in the most ancestral Type II reaction center subunit. I then evaluate the evolution of reaction centers in the context of the rise and expansion of the different groups of bacteria based on recent large-scale phylogenetic analyses. I find that the Heliobacteriaceae family of Firmicutes appears to be the earliest branching of the known groups of phototrophic bacteria; however, the origin of photochemical reaction centers and chlorophyll synthesis cannot be placed in this group. Moreover, it becomes evident that the Acidobacteria and the Proteobacteria shared a more recent common phototrophic ancestor, and this is also likely for the Chloroflexi and the Cyanobacteria. Finally, I argue that the discrepancies among the phylogenies of the reaction center proteins, chlorophyll synthesis enzymes, and the species tree of bacteria are best explained if both types of photochemical reaction centers evolved before the diversification of the known phyla of phototrophic bacteria. The primordial phototrophic ancestor must have had both Type I and Type II reaction centers.