# Unique Evolution ofSymbiobacterium thermophilum Suggested from Gene Content and Orthologous Protein Sequence Comparisons

Authors: Kenro Oshima; Kenji Ueda; Teruhiko Beppu; Hiromi Nishida
Journal: International Journal of Evolutionary Biology
 (2011)

---

## Abstract

Comparisons of gene content and orthologous protein sequence constitute a major strategy in whole-genome comparison studies. It is expected that horizontal gene transfer between phylogenetically distant organisms and lineage-specific gene loss have greater influence on gene content-based phylogenetic analysis than orthologous protein sequence-based phylogenetic analysis. To determine the evolution of the syntrophic bacteriumSymbiobacterium thermophilum, we analyzed phylogenetic relationships among Clostridia on the basis of gene content and orthologous protein sequence comparisons. These comparisons revealed that these 2 phylogenetic relationships are topologically different. Our results suggest that each Clostridia has a species-specific gene content because frequent genetic exchanges or gene losses have occurred during evolution. Specifically, the phylogenetic positions of syntrophic Clostridia were different between these 2 phylogenetic analyses, suggesting that large diversity in the living environments may cause the observed species-specific gene content. S. thermophilum occupied the most distant position from the other syntrophic Clostridia in the gene content-based phylogenetic tree. We identified 32 genes (14 under relaxed selection and 18 under functional constraint) evolving under Symbiobacterium-specific selection on the basis of synonymous-to-nonsynonymous substitution ratios. Five of the 14 genes under relaxed selection are related to transcription. In contrast, none of the 18 genes under functional constraint is related to transcription.


---

## Body

## 1. Introduction

Symbiobacterium thermophilum is a phylogenetically unique bacterium that effectively grows only in coculture with a cognate Geobacillus sp. [1]. 16S rDNA-based phylogenetic analysis has shown that it is actually a Gram-positive bacterium [2]. Although S. thermophilum phylogenetically belongs to Clostridia (low GC-content bacterial group), the genome of S. thermophilum has a high GC content (68.7%) [3]. Furthermore, 2 recent independent analyses concluded that Symbiobacterium affiliates with Clostridia (a class of Firmicutes): Ding et al. [4] carried out genome-context network analysis of 195 fully sequenced representative species, including S. thermophilum, and we analyzed the concatenated alignment of ribosomal protein sequences [5].In a previous phylogenetic analysis that was based on ribosomal protein sequence comparisons [5], S. thermophilum was closely related to 6 recently sequenced Clostridia that have distinct properties, that is,Carboxydothermus hydrogenoformans,Desulfitobacterium hafniense,Moorella thermoacetica, Pelotomaculum thermopropionicum,Desulfotomaculum reducens, and Syntrophomonas wolfei.  Symbiobacterium is dependent on the multiple functions of Geobacillus, including the supply of CO2 [1]. C. hydrogenoformans [6] grows by utilizing CO as a sole carbon source and water as an electron acceptor, which produces CO2 and hydrogen as waste products. D. hafniense [7] carries out anaerobic dechlorination of tetrachloroethene (PCE). M. thermoacetica [8] is an acetogenic bacterium that has been widely used to study the Wood-Ljungdahl pathway of CO and CO2 fixation (reductive acetyl-CoA pathway). P. thermopropionicum [9] is a member of a complex anaerobic microbial consortium where it catalyzes the intermediate bottleneck step by digesting volatile fatty acids (VFAs) and alcohols produced by upstream fermenting bacteria and it supplies acetate, hydrogen, and CO2 to downstream methanogenic archaea. D. reducens is an anaerobic sulfate-reducing bacterium [10]. S. wolfei is a fatty-acid-degrading hydrogen/formate-producing anaerobic bacterium [11].Comparisons of gene content and orthologous protein sequence constitute the major strategy in the whole-genome comparison study [12]. Clostridia have the large amount of bacteria. The phylogenetic position of Symbiobacterium remains uncertain in Clostridia. In this study, we reconstructed phylogenetic trees of Clostridia on the basis of the 2 different methods and compared them.


## 2. Methods

### 2.1. Phylogenetic Analysis on the Basis of Gene Content Comparisons

We used the following 51 bacteria (50 Clostridia and 1Bacillus belonging to Firmicutes) in this analysis: Alkaliphilus metalliredigens, Alkaliphilus oremlandii, Ammonifex degensii, Anaerocellum thermophilum, Anaerococcus prevotii, Bacillus subtilis, Caldicellulosiruptor saccharolyticus, CandidatusDesulforudis audaxviator, Carboxydothermus hydrogenoformans, Clostridium acetobutylicum, Clostridium beijerinckii, Clostridium botulinumA ATCC 19397, C. botulinumA ATCC 3502, C. botulinumA Hall, C. botulinumA2, C. botulinumA3 Loch Maree, C. botulinumB Eklund 17B, C. botulinumB1 Okra, C. botulinumBa4, C. botulinumE3, C. botulinumF Langeland, Clostridium cellulolyticum, Clostridium difficile630, C. difficileCD196, Clostridium kluyveriDSM 555, C. kluyveriNBRC 12016, Clostridium novyi, Clostridium perfringensATCC 13124, C. perfringensSM101, C. perfringens13, Clostridium phytofermentans, Clostridium tetani E88, Clostridium thermocellum, Coprothermobacter proteolyticus, Desulfitobacterium hafnienseDCB-2, D. hafnienseY51, Desulfotomaculum acetoxidans, Desulfotomaculum reducens, Eubacterium eligens, Eubacterium rectale, Finegoldia magna, Halothermothrix orenii, Heliobacterium modesticaldum, Moorella thermoacetica, Natranaerobius thermophilus, Pelotomaculum thermopropionicum, Symbiobacterium thermophilum, Syntrophomonas wolfei, Thermoanaerobacter pseudethanolicus, Thermoanaerobactersp. X514, and Thermoanaerobacter tengcongensis. Ortholog cluster analysis among the above 51 bacteria was performed using the MBGD [13] (Microbial Genome Database for Comparative Analysis; http://mbgd.nibb.ac.jp/). The analysis (minimum cluster size, 2) provided a gene presence/absence data matrix (10,636 genes × 51 organisms), which served as the basis for a distance matrix between all pairs of the 51 organisms. The distance was calculated from the different ratios between the presence/absence patterns of the 10,636 genes. On the basis of distance matrix, a neighbor-joining tree was reconstructed using MEGA software version 4 [14]. The bootstrap was performed with 1000 replicates.


### 2.2. Phylogenetic Analysis on the Basis of 112 Orthologous Protein Sequence Comparisons

We used the following 55 bacteria (54 Clostridia and 1Bacillus) in this analysis: Acidaminococcus fermentans, A. metalliredigens, A. degensii, A. thermophilum, A. prevotii, B. subtilis, C. saccharolyticus, CandidatusD. audaxviator, C. hydrogenoformans, Clostridiales genomosp. BVAB3 UPII9-5, C. acetobutylicum, C. beijerinckii, C. botulinumA ATCC 19397, C. botulinumA ATCC 3502, C. botulinumA Hall, C. botulinumA2 Kyoto, C. botulinumA3 Loch Maree, C. botulinumB Eklund 17B, C. botulinumB1 Okra, C. botulinumBa4 657, C. botulinumE3 Alaska E43, C. botulinumF Langeland, C. cellulolyticum, C. difficile630, C. difficileCD196, C. difficileR20291, C. kluyveriDSM 555, C. kluyveriNBRC 12016, C. novyi, C. perfringensATCC 13124, C. perfringensSM101, C. perfringens13, C. phytofermentans, C. tetani, C. thermocellum, C. proteolyticus, D. hafnienseDCB-2, D. hafnienseY51, D. acetoxidans, D. reducens, E. eligens, E. rectale, F. magna, H. orenii, H. modesticaldum, M. thermoacetica, N. thermophilus, P. thermopropionicum, S. thermophilum, S. wolfei, Thermoanaerobacter italicus, T. pseudethanolicus, T.sp.X514, T. tengcongensis, and Veillonella parvula. From the above 55 bacteria, 112 proteins were extracted as orthologous proteins by using a previously described method [15]. Thus, we constructed 112 multiple alignments using Clustal W [16]. Then, a concatenated multiple alignment of the 112 multiple alignments was generated. The complete multiple alignment had 52,204 amino acid sites, including 19,818 gap/insertion sites. Hence, phylogenetic analyses were performed on the basis of 32,386 amino acid sites without the gap/insertion sites. The neighbor-joining tree was reconstructed using MEGA software version 4 [14]. The bootstrap was performed with 1000 replicates. The rate variation among sites was considered to have a gamma-distributed rate (α=1). The other default parameters (e.g., Poisson distance) were not changed.


### 2.3. Extraction of Genes Evolving under Symbiobacterium–Specific Selection among Syntrophic Clostridia

AmongBacillus subtilis, Carboxydothermus hydrogenoformans,Desulfitobacterium hafniense,Moorella thermoacetica, Pelotomaculum thermopropionicum,Desulfotomaculum reducens, Symbiobacterium thermophilum, and Syntrophomonas wolfei, 472 genes were extracted as orthologous genes by the previously described method [15]. Synonymous substitution occurs more frequently than nonsynonymous substitution in protein-coding sequences because of relaxed functional constraints (nonsynonymous-to-synonymous ratio ω<1) [17], whereas they occur equally in noncoding regions and pseudogenes (ω=1). We calculated the likelihood of both the codon substitution model allowing for one ω (model R1) and the S. thermophilum branch-specific model allowing for 2 ratios (ω0 and ω1; model R2), using PAML version 3.14 [18]. In model R2, the branches of the gene tree were partitioned into the Symbiobacterium branch (ω1) and other related branches (ω0). Likelihood ratio test statistics were calculated as twice the difference between the 2 log-likelihoods (2Δln) and compared with a χ2 distribution with degrees of freedom equal to the difference in the number of parameters between the 2 models [19]. According to this method, the genes evolving under the Symbiobacterium-specific selection among Bacillus and 7 Clostridia were extracted.


## 2.1. Phylogenetic Analysis on the Basis of Gene Content Comparisons

We used the following 51 bacteria (50 Clostridia and 1Bacillus belonging to Firmicutes) in this analysis: Alkaliphilus metalliredigens, Alkaliphilus oremlandii, Ammonifex degensii, Anaerocellum thermophilum, Anaerococcus prevotii, Bacillus subtilis, Caldicellulosiruptor saccharolyticus, CandidatusDesulforudis audaxviator, Carboxydothermus hydrogenoformans, Clostridium acetobutylicum, Clostridium beijerinckii, Clostridium botulinumA ATCC 19397, C. botulinumA ATCC 3502, C. botulinumA Hall, C. botulinumA2, C. botulinumA3 Loch Maree, C. botulinumB Eklund 17B, C. botulinumB1 Okra, C. botulinumBa4, C. botulinumE3, C. botulinumF Langeland, Clostridium cellulolyticum, Clostridium difficile630, C. difficileCD196, Clostridium kluyveriDSM 555, C. kluyveriNBRC 12016, Clostridium novyi, Clostridium perfringensATCC 13124, C. perfringensSM101, C. perfringens13, Clostridium phytofermentans, Clostridium tetani E88, Clostridium thermocellum, Coprothermobacter proteolyticus, Desulfitobacterium hafnienseDCB-2, D. hafnienseY51, Desulfotomaculum acetoxidans, Desulfotomaculum reducens, Eubacterium eligens, Eubacterium rectale, Finegoldia magna, Halothermothrix orenii, Heliobacterium modesticaldum, Moorella thermoacetica, Natranaerobius thermophilus, Pelotomaculum thermopropionicum, Symbiobacterium thermophilum, Syntrophomonas wolfei, Thermoanaerobacter pseudethanolicus, Thermoanaerobactersp. X514, and Thermoanaerobacter tengcongensis. Ortholog cluster analysis among the above 51 bacteria was performed using the MBGD [13] (Microbial Genome Database for Comparative Analysis; http://mbgd.nibb.ac.jp/). The analysis (minimum cluster size, 2) provided a gene presence/absence data matrix (10,636 genes × 51 organisms), which served as the basis for a distance matrix between all pairs of the 51 organisms. The distance was calculated from the different ratios between the presence/absence patterns of the 10,636 genes. On the basis of distance matrix, a neighbor-joining tree was reconstructed using MEGA software version 4 [14]. The bootstrap was performed with 1000 replicates.


## 2.2. Phylogenetic Analysis on the Basis of 112 Orthologous Protein Sequence Comparisons

We used the following 55 bacteria (54 Clostridia and 1Bacillus) in this analysis: Acidaminococcus fermentans, A. metalliredigens, A. degensii, A. thermophilum, A. prevotii, B. subtilis, C. saccharolyticus, CandidatusD. audaxviator, C. hydrogenoformans, Clostridiales genomosp. BVAB3 UPII9-5, C. acetobutylicum, C. beijerinckii, C. botulinumA ATCC 19397, C. botulinumA ATCC 3502, C. botulinumA Hall, C. botulinumA2 Kyoto, C. botulinumA3 Loch Maree, C. botulinumB Eklund 17B, C. botulinumB1 Okra, C. botulinumBa4 657, C. botulinumE3 Alaska E43, C. botulinumF Langeland, C. cellulolyticum, C. difficile630, C. difficileCD196, C. difficileR20291, C. kluyveriDSM 555, C. kluyveriNBRC 12016, C. novyi, C. perfringensATCC 13124, C. perfringensSM101, C. perfringens13, C. phytofermentans, C. tetani, C. thermocellum, C. proteolyticus, D. hafnienseDCB-2, D. hafnienseY51, D. acetoxidans, D. reducens, E. eligens, E. rectale, F. magna, H. orenii, H. modesticaldum, M. thermoacetica, N. thermophilus, P. thermopropionicum, S. thermophilum, S. wolfei, Thermoanaerobacter italicus, T. pseudethanolicus, T.sp.X514, T. tengcongensis, and Veillonella parvula. From the above 55 bacteria, 112 proteins were extracted as orthologous proteins by using a previously described method [15]. Thus, we constructed 112 multiple alignments using Clustal W [16]. Then, a concatenated multiple alignment of the 112 multiple alignments was generated. The complete multiple alignment had 52,204 amino acid sites, including 19,818 gap/insertion sites. Hence, phylogenetic analyses were performed on the basis of 32,386 amino acid sites without the gap/insertion sites. The neighbor-joining tree was reconstructed using MEGA software version 4 [14]. The bootstrap was performed with 1000 replicates. The rate variation among sites was considered to have a gamma-distributed rate (α=1). The other default parameters (e.g., Poisson distance) were not changed.


## 2.3. Extraction of Genes Evolving under Symbiobacterium–Specific Selection among Syntrophic Clostridia

AmongBacillus subtilis, Carboxydothermus hydrogenoformans,Desulfitobacterium hafniense,Moorella thermoacetica, Pelotomaculum thermopropionicum,Desulfotomaculum reducens, Symbiobacterium thermophilum, and Syntrophomonas wolfei, 472 genes were extracted as orthologous genes by the previously described method [15]. Synonymous substitution occurs more frequently than nonsynonymous substitution in protein-coding sequences because of relaxed functional constraints (nonsynonymous-to-synonymous ratio ω<1) [17], whereas they occur equally in noncoding regions and pseudogenes (ω=1). We calculated the likelihood of both the codon substitution model allowing for one ω (model R1) and the S. thermophilum branch-specific model allowing for 2 ratios (ω0 and ω1; model R2), using PAML version 3.14 [18]. In model R2, the branches of the gene tree were partitioned into the Symbiobacterium branch (ω1) and other related branches (ω0). Likelihood ratio test statistics were calculated as twice the difference between the 2 log-likelihoods (2Δln) and compared with a χ2 distribution with degrees of freedom equal to the difference in the number of parameters between the 2 models [19]. According to this method, the genes evolving under the Symbiobacterium-specific selection among Bacillus and 7 Clostridia were extracted.


## 3. Results and Discussion

Phylogenetic relationships among Clostridia on the basis of gene content comparison (Figure1) were topologically different from those generated on the basis of orthologous protein sequence comparison (Figure 2). For example, in the gene content-based phylogenetic tree, Alkaliphilus, Clostridium (except for C. cellulolyticum and C. thermocellum), Desulfitobacterium, and Eubacterium formed a monophyletic lineage with 85% bootstrap support (Figure 1). In contrast, in the 112 orthologous protein sequence-based phylogenetic relationships, Alkaliphilus, Anaerococcus,Clostridium (except for C. cellulolyticum and C. thermocellum), Eubacterium, and Finegoldia formed a monophyletic lineage with 98% bootstrap support (Figure 2). Thus, the phylogenetic positions of Anaerococcus, Desulfitobacterium, and Finegoldiawere different between these 2 trees. In addition, Coprothermobacter proteolyticus was positioned differently in the 2 trees. Moreover, the very long branch in the orthologous protein-based tree suggests that C. proteolyticus has a substitution pattern that is different from other related Clostridia.Figure 1
Phylogenetic relationships on the basis of gene content comparisons among 50 Clostridia andBacillus subtilis. The ortholog cluster analysis (minimum cluster size, 2) among the 51 bacteria was performed using the MBGD [13]. This analysis produced the gene presence/absence data matrix (10,636 genes × 51 organisms), which was used to generate the distance matrix between all pairs of the 51 bacteria. On the basis of the distance matrix, a neighbor-joining tree was reconstructed using MEGA software version 4 [14]. The bootstrap was performed with 1000 replicates. The bar indicates a 200-gene difference.Figure 2
Phylogenetic relationships on the basis of 112 orthologous protein sequence comparisons among 54 Clostridia andB. subtilis. The 112 proteins were extracted as orthologous proteins from the 55 bacteria by a previously described method [15]. We constructed the 112 multiple alignments by using Clustal W [16]. Then, a concatenated multiple alignment of the 112 multiple alignments was generated. The complete multiple alignment had 52,204 amino acid sites, including 19,818 gap/insertion sites. Hence, phylogenetic analyses were performed on the basis of 32,386 amino acid sites without the gap/insertion sites. The neighbor-joining tree was reconstructed using MEGA software version 4 [14]. The bootstrap was performed with 1000 replicates. The rate variation among sites was assumed to have a gamma distributed rate (α=1). No other default parameters were changed. The bar indicates a 10% difference.We expected horizontal gene transfer between phylogenetically distant organisms and lineage-specific gene loss to have greater influence on the gene content-based phylogenetic analysis than the orthologous protein-based analysis [12, 20]. Bacteria make their gene content suitable for the living environment by changing it through gene acquisition and loss.The phylogenetic positions of 2D.  hafniense strains are located near those of Alkaliphilus, Clostridium (except for C. cellulolyticum and C. thermocellum), and Eubacterium in the gene content-based phylogenetic tree (Figure 1). However, those phylogenetic positions were located in the phylogenetic lineage of syntrophic Clostridia in the orthologous protein-based tree (Figure 2). The gene content-based phylogenetic tree (Figure 1) indicates that Symbiobacterium branched off at the earliest stage of Clostridia species diversification. In contrast, Natranaerobius branched off at the earliest species diversification stage in the orthologous protein sequence-based phylogenetic tree (Figure 2).AlthoughS. thermophilum occupied the most basal position in the gene content-based Clostridia lineage (Figure 1), it was located in the syntrophic Clostridia lineage on the basis of orthologous protein sequence comparisons (Figure 2). Syntrophic bacteria evolved to acquire different sets of genes despite their close phylogenetic relationship. Thus, although Symbiobacterium clusters with syntrophic Clostridia, its gene content is very different. S. thermophilum has the most distant position from the other syntrophic Clostridia in the phylogenetic tree on the basis of gene content comparisons.Although the physiological reason for the high CO2 requirement of S. thermophilum is not yet known, we assumed that it is related to the carbonic anhydrase deficiency (the ubiquitous enzyme catalyzing interconversion between CO2 and bicarbonate; EC 4.2.1.1), as deficiency of this enzyme results in the need for high CO2 levels in several model microorganisms [1]. S. thermophilum lost this enzyme in the course of evolution [5]. In this previous analysis, we inferred that C. hydrogenoformans and M. thermoacetica have also lost the gene for carbonic anhydrase; however, we recently noticed that C. hydrogenoformans had 2 potential carbonic anhydrase coding genes with structures different from the other syntrophic Clostridia carbonic anhydrases. Therefore, only Moorella has lost the carbonic anhydrase gene, in addition to Symbiobacterium. However, according to our results, these two bacteria are not closely related to each other (Figures 1 and 2), suggesting that the gene loss in these 2 species occurred independently during evolution.Our results imply that each syntrophic Clostridial organism, especiallySymbiobacterium, would have genes that evolved in an organism-specific manner. We expect that characterization of such genes will provide useful information with regard to the evolutionary history and physiological features specific to the corresponding organism [21, 22]. We identified 32 genes evolving under Symbiobacterium-specific selection (Table 1). The analysis revealed that the likelihood of model R2 was significantly higher (P<.05) than that of model R1 in the 32 genes. Of these, 14 genes showed ω1/ω0>1 and 18 showed ω1/ω0<1.Table 1
Genes evolving underSymbiobacterium-specific selection.
Geneω1ω1/ω02Δlnω1/ω0>1hfq (RNA chaperone, STH1746)0.534724.30465.7413spoIIAB (anti-sigma F factor, STH1813)0.39675.97448.7835flgG (flagellar hook protein, STH2995)0.377475.48006.4323ilvC (ketol-acid reductoisomerase, STH2688)0.22403.467510.7272rplL (50S ribosomal protein L7/L12, STH3086)0.21838.364013.4750argD (N-acetylornithine aminotransferase, STH2881)0.20842.02924.1224rplK (50S ribosomal protein L11, STH3090)0.18699.34504.2192ylmE (alanine racemase domain-containing protein, STH1227)0.152624.222215.4681proJ (gamma-glutamyl kinase, STH2540)0.149726.26324.4715sigA (RNA polymerase sigma factor, STH0588)0.13153.767917.9996rpoC (RNA polymerase subunit beta', STH3084)0.08382.14877.2876glmS (glucosamine-fructose-6-phosphate aminotransferase, STH1279)0.01562.785713.0700aroE (3-phosphoshikimate 1 carboxyvinyltransferase, STH1419)0.01252.84094.4748cheY (two-component response regulator involved in modulation of flagellar, STH1540)0.00442.93336.7786ω1/ω0<1flgD (flagellar hook assembly protein, STH2996)0.01230.07154.4609fliS (flagellar protein FliS, STH2976)0.00730.08854.0842yloM (ribosomal RNA small subunit methyltransferase B, STH1349)0.00450.044112.0081ftsH (cell division protease, STH3198)0.00400.065511.9908spoVFB (dipicolinate synthase subunit B, STH1546)0.00390.05916.5852rplW (50S ribosomal protein L23, STH3073)0.00390.14293.9835trmD (tRNA methyltransferase, STH1470)0.00380.05745.7865argC (N-acetyl-gamma-glutamyl-phosphate reductase, STH2892)0.00380.07214.1368rpsC (30S ribosomal protein S3, STH3069)0.00370.15044.1064prfA (peptide chain release factor RF-1, STH0073)0.00350.07506.3618ligA (NAD-dependent DNA ligase, STH2825)0.00340.06544.5717spo0J (ParB-like nuclease domain-containing protein, STH3332)0.00340.039710.1363ftsE (cell division ATP-binding protein, STH0139)0.00270.04075.6285metG (methionyl-tRNA synthetase, STH3252)0.00270.04704.1885rplC (50S ribosomal protein L3, STH3075)0.00230.09204.3304rplB (50S ribosomal protein L2, STH3072)0.00200.06173.9836rpsH (30S ribosomal protein S8, STH3061)0.00180.10474.5829infA (translation initiation factor IF-1, STH3052)0.00040.02344.7842Among the 32 genes evolving underSymbiobacterium-specific selection, the RNA chaperone Hfq-coding gene has the highest ω1 value (0.5347) (Table 1). Hfq facilitates pairing interactions between small regulatory RNAs and their mRNA targets, which has a variety of functions in bacteria [23]. Among 73 conserved amino acid sites of Hfq (Figure 3), S. thermophilum has more specific sites (7 sites) than the outgroup Bacillus (4 sites), indicating that the Hfq gene is one of the genes evolving under Symbiobacterium-specific selection.Figure 3
Alignment of amino acid sequences of Hfq.Sth,  Symbiobacterium thermophilum; Bsu,  Bacillus subtilis; Chy,  Carboxydothermus hydrogenoformans; Dre,  Desulfotomaculum reducens; Dha, Desulfitobacterium hafniense; Mth, Moorella thermoacetica; Pth, Pelotomaculum thermopropionicum; Swo,  Syntrophomonas wolfei. Red and blue sites indicate Symbiobacterium- and Bacillus-specific sites, respectively. The dots represent identical residues of S. thermophilum amino acid.Two genes related to transcription,sigA (RNA polymerase sigma factor coding gene) and rpoC (RNA polymerase subunit beta’ coding gene) have evolved under relaxed selection (Table 1). These results could be related to the high GC content of Symbiobacterium genes. Thus, we hypothesized that the GC bias of the promoter sequence induced Symbiobacterium-specific SigA, a DNA-binding protein, which led to the structural change of RNA polymerase complex (including RpoC). We discussed the relationships between the GC content and phylogeny of the Symbiobacterium genes [24].In addition,spoIIAB and cheY are also related to transcription. Thus, 5 of the 14 genes under more relaxed selection than other Clostridia are related to transcription. However, none of the 18 genes under functional constraint is related to transcription. Those results suggest that, under relaxed selection, the transcription system may be related to S. thermophilum-specific gene content. In fact, Symbiobacterium lost the transcriptional regulator genes arsR, GntR, and Lrp compared to other syntrophic Clostridia (See in the Supplementary Material available online at doi: 10.4061/2011/376831 Table S1.).It is noteworthy that some functionally related genes exhibited opposite nucleotide substitution patterns inS. thermophilum (Table 1). For example, argD (N-acetylornithine aminotransferase coding gene) has evolved under relaxed selection whereas argC (N-acetyl-gamma-glutamyl-phosphate reductase coding gene) has evolved under functional constraint. Another example is the genes encoding flagella-associated proteins; flgG (flagellar hook protein coding gene) has evolved under relaxed selection, whereas flgD (flagellar hook assembly protein coding gene) and fliS (flagellar protein coding gene) have evolved under functional constraint. flgG exhibited the highest ω1/ω0 value (75.48) (Table 1). Flagella mediate interactions between P. thermopropionicum and methanogenic archaea [25]. Similar specialized functions in syntrophic association could have been a limiting factor for the evolution of the above 2 flagellum genes in Symbiobacterium.In conclusion, our results suggest thatS. thermophilum has evolved in a unique manner compared to other syntrophic Clostridia from the perspective of gene content. Codon substitution analysis also suggests several unique genes that evolved in a Symbiobacterium-specific manner. Although speculative, the gene loss or relaxed evolution of several transcriptional regulator genes implies that environmental response might be involved in Symbiobacterium-specific evolution.



---

*Source: 376831-2010-12-22.xml*