Comparison of the protein-coding genomes of three deep-sea, sulfur-oxidising bacteria: “Candidatus Ruthia magnifica”, “Candidatus Vesicomyosocius okutanii” and Thiomicrospira crunogena
BMC Research Notes volume 10, Article number: 296 (2017)
“ Candidatus Ruthia magnifica”, “Candidatus Vesicomyosocius okutanii” and Thiomicrospira crunogena are all sulfur-oxidising bacteria found in deep-sea vent environments. Recent research suggests that the two symbiotic organisms, “Candidatus R. magnifica” and “Candidatus V. okutanii”, may share common ancestry with the autonomously living species T. crunogena. We used comparative genomics to examine the genome-wide protein-coding content of all three species to explore their similarities. In particular, we used the OrthoMCL algorithm to sort proteins into groups of putative orthologs on the basis of sequence similarity.
The OrthoMCL inflation parameter was tuned using biological criteria. Using the tuned value, OrthoMCL delimited 1070 protein groups. 63.5% of these groups contained one protein from each species. Two groups contained duplicate protein copies from all three species. 123 groups were unique to T. crunogena and ten groups included multiple copies of T. crunogena proteins but only single copies from the other species. “Candidatus R. magnifica” had one unique group, and had multiple copies in one group where the other species had a single copy. There were no groups unique to “Candidatus V. okutanii”, and no groups in which there were multiple “Candidatus V. okutanii” proteins but only single proteins from the other species. Results align with previous suggestions that all three species share a common ancestor. However this is not definitive evidence to make taxonomic conclusions and the possibility of horizontal gene transfer was not investigated. Methodologically, the tuning of the OrthoMCL inflation parameter using biological criteria provides further methods to refine the OrthoMCL procedure.
“Candidatus Ruthia magnifica” strain Cm, “Candidatus Vesicomyosocius okutanii” HA and Thiomicrospira crunogena XCL-2 are all sulfur-oxidising bacteria found in deep-sea vent environments.
“Candidatus R. magnifica” and “Candidatus V. okutanii” live symbiotically in the gill epithelial cells of giant clam species: “Candidatus R. magnifica” in Calyptogena magnifica  and “Candidatus V. okutanii” in Calyptogena okutanii . It is predicted that they are predominantly transmitted vertically via their host’s eggs [3, 4]. These hosts have reduced or vestigial digestive tracts and are therefore dependent on their symbionts for their nutritional requirements. As both giant clam species reside in deep-sea vent environments their symbionts are able to utilise the sulfur, produced by the vents, to provide their hosts with carbon and other nutrients [1, 2]. The symbionts’ dependence on the host varies, for example “Candidatus R. magnifica” encodes pathways to synthesise 20 amino acids , whereas “Candidatus V. okutanii” encodes pathways for 18 amino acids . It has been hypothesised that missing essential genes in the symbiont may help maintain a stable symbiont population in a host cell .
Recent sequence-based reconstructions of phylogenetic trees suggest that “Candidatus R. magnifica” and “Candidatus V. okutanii” form a clade with each other, and a broader clade with T. crunogena [5, 6]. T. crunogena lives independently, though in the same deep-sea vent environments.
Preliminary to detailed studies on ancestry and adaptation among these three taxa, we can predict paralogs and orthologs across their genomes. Paralogs are genes arising by a duplication event within a species, and orthologs are genes in different taxa whose common ancestor is a gene present in the most recent common ancestral taxon [7, 8]. Although these definitions are explicitly phylogenetic, requiring a gene tree and a species tree, prediction of orthologous groups is often performed on the basis of sequence similarity alone. We investigated the evolution of the protein-coding gene content across all three species using OrthoMCL  and BLAST , to create protein groups based on sequence similarity, and UniProt , to assign functions to these groups.
Compared to an earlier comparative genomics study including the three species , our methodology allows more detailed investigation of variation in gene copy number. In contrast to purely reciprocal-best methods which predict only 1:1 orthologous relationships across taxa, OrthoMCL groups putative paralogs into orthologous groups of two or more sequences, imposing no upper limit on group size and no requirement that each group be present in each species.
The 4273π variant of the Raspbian Linux operating system  was used on a Raspberry Pi computer (Version 1, Model B, Revision 2.0). Genome-wide protein sets for “Candidatus Ruthia magnifica” strain Cm, “Candidatus Vesicomyosocius okutanii” HA and Thiomicrospira crunogena XCL-2 were downloaded, in FASTA format, from the Ensembl genomes database (http://ensemblgenomes.org) [14, 15]. OrthoMCL software (http://orthomcl.org)  and MCL  were used to delimit protein groups based on sequence similarity.
The OrthoMCL procedure was followed as outlined in the OrthoMCL user guide, with the exception of using the substitution matrix BLOSUM45 for the ‘all-versus-all’ NCBI BLAST  and omitting BLAST’s -z parameter; and for our final analysis, the inflation value (I) was set to 1.4 when running MCL.
As the inflation value decreases, sequences are included in fewer, larger groups, reducing the tightness and granularity of the delimited groups. To determine the optimal value of the inflation parameter, first a range of values were tested (I = 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2, 4, 6, 10 and 12). Annotated functions of the first (largest) three protein groups were examined using UniProt (http://uniprot.org) . This revealed that only values 1.2–1.9 gave rise to groups that are both functionally cohesive and inclusive. Within this range, group 1, for example, contains diguanylate phosphodiesterases/cyclases (Table 1) but when I is increased to 2 these are split into different groups. Furthermore, when I = 1.1 proteins with different functions (transcriptional regulators—winged helix family) are also included within this group.
Results using I = 1.2, 1.3 and 1.4 gave rise to a group which was not present using other values. It had multiple copies of “Candidatus R. magnifica” proteins but only single copies of “Candidatus V. okutanii” and T. crunogena proteins. UniProt suggested that these proteins had the same function (histidinol-phosphate aminotransferases) but increasing the I value split them up into different groups. Hence on biological grounds, I > 1.4 was rejected. I = 1.4 was used for the final analysis presented here, as it gave the strongest restraints on group formation while still maintaining this aminotransferase group (Table 1, Group 27).
Once groups were delimited, a Perl script, modified from , was used to count the number of times each species was represented in each protein group . To verify the reliability of the script, our final OrthoMCL groups file was used as input for an independently-written protein-counting script (Kevin Kiesworo, unpublished) and the same counts were obtained. Additionally, the OrthoMCL groups file from a similar study on different taxa (Hannah Currant, unpublished) was used as input for our script, and the same counts were obtained as in that study.
Functions of both the largest and most interesting groups were then inferred by searching for protein accessions in UniProt: group members were searched until a common function was found between at least four proteins, or for smaller groups all members were searched (Table 1).
OrthoMCL predicted 1070 protein groups based on sequence similarity . 63.5% of these contained a single protein from each of the three species. Two groups had duplicate protein copies in all three species (Table 1). T. crunogena had 123 unique groups, and multiple copies in ten groups that only had single copies in the other two species. “Candidatus R. magnifica” had one unique group, and had multiple copies in one group where the other species had a single copy. “Candidatus V. okutanii” had no unique groups. Nor did it have multiple copies in any groups that only had single protein copies from the other species. There were no groups that contained multiple protein copies from two species and one copy in the third (Table 2).
Similarities between all three species
679 of the 1070 protein groups delimited by OrthoMCL (63.5%) contained a single protein copy in each of the three species. This high degree of similarity could be a consequence of common ancestry, horizontal transfer in a shared habitat, or most likely a mixture of both.
All three species inhabit deep-sea vent environments which are highly variable and have constantly fluctuating factors such as sulfur and carbon concentrations . In order to survive, the organisms must possess methods that allow them to deal with such fluctuations. One shared process, for example, is their ability to oxidise the sulfur supplied by deep-sea vents to fix carbon for use in cellular functions.
Two groups were predicted to contain duplicate protein copies in all three species (Table 1, Groups 9 and 10). These are consistent with duplication in a common ancestor, with subsequent speciations, although our current work does not distinguish this from other possibilities such as horizontal transfer.
All three species have a duplicate copy of elongation factor Tu (Table 1, Group 9). These paralogs are found in all proteobacteria and it has been hypothesised, therefore, that this duplication event preceded the divergence of this phylum . It has been shown that the tuf genes that encode these proteins undergo gene conversion  which inhibits any divergence, and therefore sub- or neo-functionalisation, of the two genes . The persistence of the duplicate may therefore indicate high levels of expression. Detection of this known group is promising in regards to the reliability of our methods.
Each species also has a duplicate in the group of nitrogen regulatory proteins (P-II; Table 1, Group 10). One of these copies, in “Candidatus V. okutanii”, is the product of the glnK gene. This gene seems to be commonly duplicated in some sub-divisions of proteobacteria  and its evolution seems to be associated with that of the amtB ammonium transporter gene, to which it is physically and functionally linked . Interestingly, these ammonium transporters make up Group 11 (Table 1) and although there are four copies in T. crunogena, the other species have no duplicates. This would be consistent with genome reduction due to a symbiotic lifestyle , although our current work cannot distinguish this with certainty.
Unique to T. crunogena
A large number of protein groups (123) are found only in T. crunogena (Table 2). Also, in 10 groups, T. crunogena has multiple protein copies where the other species only have one copy each (Table 1). As the only independent-living species studied, T. crunogena may require a larger number of genes and proteins for survival. The other, symbiotic, organisms can rely on their hosts to provide some essential functions and, therefore, loss of some genes could prove to be energetically favourable . There is also a lower total protein count for these species (976 and 937 protein sequences, compared to 2196 in T. crunogena).
Unique to “Candidatus R. magnifica”
“Candidatus R. magnifica” has one unique group that consists of glycosyl transferases (Table 1, Group 746).
There was also one group delimited that had multiple “Candidatus R. magnifica” proteins and only single proteins from the other species (Table 1, Group 27). This is a group of histidinol-phosphate aminotransferases.
Unique to “Candidatus V. okutanii”
“Candidatus V. okutanii” has no unique groups or paralogs.
The number of unique protein groups found in T. crunogena may highlight its independent lifestyle that is very different from the other, symbiont, species. On the other hand, all three species shared many groups in common that could be indicative of a shared common ancestor, as was previously hypothesised [5, 6]. However, sequence-based orthology prediction is not sufficient to resolve taxonomy .
Methodologically, our work extends comparative genomics on the low-cost Raspberry Pi computer in two main ways. Firstly, three species were used, as opposed to two species in earlier studies [17, 25]. With faster, more recent versions of the hardware, such as the Raspberry Pi 3, even larger numbers of species would be possible.
Secondly, in our current study the OrthoMCL inflation parameter has been tuned using the biological criterion of functional coherence of the first (largest) three protein groups. This contrasts with algorithmic criteria used by, for example,  and , and may be generalizable to other methods for delimiting groups that also use MCL , for example Orthofinder . There may be no universally optimal way to set the inflation parameter. However, biological criteria will always be valuable, whether used alone or to verify an algorithmic approach. Methodologically, combining biological criteria to guide the choice of inflation parameter with other refinements in family prediction (e.g. ) may be a promising future direction.
The method used only utilises sequence-based orthology prediction to produce protein groups, without phylogeny reconstruction, and so it is not sufficient to resolve taxonomy . In accordance with the International Code of Nomenclature of Bacteria , other information such as metabolic and reproductive features must be known before formal taxonomy can be assigned.
We are also unable to rule out the possibility that the similarities in the protein coding content of these three genomes were due to horizontal gene transfer. It is thought that these events are less common as symbionts of vesicomyid clams (such as Calyptogena magnifica and Calyptogena okutanii) are found in their host oocytes—suggesting that vertical transmission is predominant [3, 4], presumably reducing opportunities for horizontal gene transfer. However, there is also evidence that lateral transmission, and therefore horizontal gene transfer events, can occur . Analysis of horizontal transfer among these species and their relatives, including investigation of the detail of horizontal transfers , would be a promising future direction.
- “Candidatus R. magnifica”:
“Candidatus Ruthia magnifica”
- “Candidatus V. okutanii”:
“Candidatus Vesicomyosocius okutanii”
- T. crunogena :
Newton ILG, Woyke T, Auchtung TA, Dilly GF, Dutton RJ, Fisher MC, Fontanez KM, Lau E, Stewart FJ, Richardson PM, Barry KW, Saunders E, Detter JC, Wu D, Eisen JA, Cavanaugh CM. The Calyptogena magnifica chemoautotrophic symbiont genome. Science. 2007;315:998–1000.
Kuwahara H, Yoshida T, Takaki Y, Shimamura S, Nishi S, Harada M, Matsuyama K, Takishita K, Kawato M, Uematsu K, Fujiwara Y, Sato T, Kato C, Kitagawa M, Kato I, Maruyama T. Reduced genome of the thioautotrophic intracellular symbiont in a deep-sea clam, Calyptogena okutanii. Curr Biol. 2007;17:881–6.
Endow K, Ohta S. Occurrence of bacteria in the primary oocytes of vesicomyid clam Calyptogena soyoae. Mar Ecol Prog Ser. 1990;64:309–11.
Cary SC, Giovannoni SJ. Transovarial inheritance of endosymbiotic bacteria in clams inhabiting deep-sea hydrothermal vents and cold seep. Proc Natl Acad Sci USA. 1993;90:5695–9.
Williams KP, Gillespie JJ, Sobral BWS, Nordberg EK, Snyder EE, Shallom JM, Dickerman AW. Phylogeny of gammaproteobacteria. J Bacteriol. 2010;192:2305–14.
Wu DY, Hugenholtz P, Mavromatis K, Pukall R, Dalin E, Ivanova NN, Kunin V, Goodwin L, Wu M, Tindall BJ, Hooper SD, Pati A, Lykidis A, Spring S, Anderson IJ, D’haeseleer P, Zemla A, Singer M, Lapidus A, Nolan M, Copeland A, Han C, Chen F, Cheng J-F, Lucas S, Kerfeld C, Lang E, Gronow S, Chain P, Bruce D, et al. A phylogeny-driven genomic encyclopaedia of Bacteria and Archaea. Nature. 2009;462:1056–60.
Fitch WM. Homology: a personal view on some of the problems. Trends Genet. 2000;16:227–31.
Koonin EV. Orthologs, paralogs and evolutionary genomics. Annu Rev Genet. 2005;39:309–38.
Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.
The UniProt Consortium. UniProt: a hub for protein information. Nucleic Acids Res. 2015;43:D204–12.
Newton ILG, Girguis PR, Cavanaugh CM. Comparative genomics of vesicomyid clam (Bivalvia: Mollusca) chemosynthetic symbionts. BMC Genom. 2008;9:585.
Barker D, Ferrier DEK, Holland PWH, Mitchell JBO, Plaisier H, Ritchie MG, Smart SD. 4273π: bioinformatics education on low cost ARM hardware. BMC Bioinform. 2013;14:243.
Cunningham F, Amode MR, Barrell D, Beal K, Billis K, Brent S, Carvalho-Silva D, Clapham P, Coates G, Fitzgerald S, Gil L, García Girón C, Gordon L, Hourlier T, Hunt SE, Janacek SH, Johnson N, Juettemann T, Kähäri AK, Keenan S, Martin FJ, Maurel T, McLaren W, Murphy DN, Nag R, Overduin B, Parker A, Patricio M, Perry E, Pignatelli M, et al. Ensembl 2015. Nucleic Acids Res. 2015;43:D662–9.
McGill S, Barker D. Additional files for ‘Comparison of the protein-coding genomes of three deep-sea, sulfur-oxidising bacteria: “Candidatus Ruthia magnifica”, “Candidatus Vesicomyosocius okutanii” and Thiomicrospira crunogena’[dataset]. Edinburgh Datashare. 2017. http://dx.doi.org/10.7488/ds/2065.
van Dongen S. A cluster algorithm for graphs. Rep Inform Syst. 2000;10:1–40.
Robson JF, Barker D. Comparison of the protein-coding gene content of Chlamydia trachomatis and Protochlamydia amoebophila using a Raspberry Pi computer. BMC Res Notes. 2015;8:561.
Johnson KS, Beehler CL, Sakamoto-Arnold CM, Childress JJ. In situ measurements of chemical distributions in a deep-sea hydrothermal vent field. Science. 1986;231:1139–41.
Lathe WC, Bork P. Evolution of tuf genes: ancient duplication, differential loss and gene conversion. FEBS Lett. 2001;502:113–6.
Isabel S, Leblanc E, Boissinot M, Boudreau DK, Grondin M, Picard FJ, Martel EA, Parham NJ, Chain PSG, Bader DE, Mulvey MR, Bryden L, Roy PH, Ouellette M, Bergeron MG. Divergence among genes encoding the elongation factor Tu of Yersinia species. J Bacteriol. 2008;190:7548–58.
Kondrashov FA, Gurbich TA, Vlasov PK. Selection for functional uniformity of tuf duplicates in γ-proteobacteria. Trends Genet. 2007;23:215–8.
Jack R, De Zamaroczy M, Merrick M. The signal transduction protein GlnK is required for NifL-dependent nitrogen control of nif gene expression in Klebsiella pneumoniae. J Bacteriol. 1999;181:1156–62.
Coutts G, Thomas G, Blakey D, Merrick M. Membrane sequestration of the signal transduction protein GlnK by the ammonium transporter AmtB. EMBO J. 2002;21:536–45.
Murray RG, Stackebrandt E. Taxonomic note: implementation of the provisional status candidatus for incompletely described procaryotes. Int J Syst Bacteriol. 1995;45:186–7.
Wreggelsworth KM, Barker D. A comparison of the protein-coding genomes of two green sulphur bacteria, Chlorobium tepidum TLS and Pelodictyon phaeoclathratiforme BU-1. BMC Res Notes. 2015;8:565.
Swingley WD, Blankenship RE, Raymond J. Integrating Markov clustering and molecular phylogenetics to reconstruct the cyanobacterial species tree from conserved protein families. Mol Biol Evol. 2008;25:643–54.
Latysheva N, Junker VL, Palmer WJ, Codd GA, Barker D. The evolution of nitrogen fixation in cyanobacteria. Bioinformatics. 2012;28:603–6.
Emms DM, Kelly S. OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015;16:157.
Nowell RW, Green S, Laue BE, Sharp PM. The extent of genome flux and its role in the differentiation of bacterial lineages. Genome Biol Evol. 2014;6:1514–29.
Lapage SP, Sneath PHA, Lessel EF, Skerman VBD, Seeliger HPR, Clark WA. International code of nomenclature of bacteria: bacteriological code, 1990 revision. Washington, DC: ASM Press; 1992.
Stewart FJ, Young CR, Cavanaugh CM. Lateral symbiont acquisition in a maternally transmitted chemosynthetic clam endosymbiosis. Mol Biol Evol. 2008;25:673–87.
Darby CA, Stolzer M, Ropp PJ, Barker D, Durand D. Xenolog classification. Bioinformatics. 2017;33:640–9.
This paper is based on work submitted by SEM as coursework for the module BL4273 Bioinformatics for Biologists, coordinated by DB at the University of St Andrews. SEM carried out the analyses. SEM and DB wrote the manuscript. Both authors read and approved the final manuscript.
We thank Hannah Currant and Kevin Kiesworo for providing their unpublished scripts, for purposes of verification.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and analysed during the current study are available in the University of Edinburgh Datashare repository, http://dx.doi.org/10.7488/ds/2065.
No specific funding was received for this research. The University of St Andrews provided funding for the open access charge.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
McGill, S.E., Barker, D. Comparison of the protein-coding genomes of three deep-sea, sulfur-oxidising bacteria: “Candidatus Ruthia magnifica”, “Candidatus Vesicomyosocius okutanii” and Thiomicrospira crunogena . BMC Res Notes 10, 296 (2017). https://doi.org/10.1186/s13104-017-2598-5