- Research article
- Open Access
Genomes of “Spiribacter”, a streamlined, successful halophilic bacterium
BMC Genomics volume 14, Article number: 787 (2013)
Thalassosaline waters produced by the concentration of seawater are widespread and common extreme aquatic habitats. Their salinity varies from that of sea water (ca. 3.5%) to saturation for NaCl (ca. 37%). Obviously the microbiota varies dramatically throughout this range. Recent metagenomic analysis of intermediate salinity waters (19%) indicated the presence of an abundant and yet undescribed gamma-proteobacterium. Two strains belonging to this group have been isolated from saltern ponds of intermediate salinity in two Spanish salterns and were named “Spiribacter”.
The genomes of two isolates of “Spiribacter” have been fully sequenced and assembled. The analysis of metagenomic datasets indicates that microbes of this genus are widespread worldwide in medium salinity habitats representing the first ecologically defined moderate halophile. The genomes indicate that the two isolates belong to different species within the same genus. Both genomes are streamlined with high coding densities, have few regulatory mechanisms and no motility or chemotactic behavior. Metabolically they are heterotrophs with a subgroup II xanthorhodopsin as an additional energy source when light is available.
This is the first bacterium that has been proven by culture independent approaches to be prevalent in hypersaline habitats of intermediate salinity (half a way between the sea and NaCl saturation). Predictions from the proteome and analysis of transporter genes, together with a complete ectoine biosynthesis gene cluster are consistent with these microbes having the salt-out-organic-compatible solutes type of osmoregulation. All these features are also consistent with a well-adapted fully planktonic microbe while other halophiles with more complex genomes such as Salinibacter ruber might have particle associated microniches.
The development of techniques to directly study genes in the environment, first by PCR and cloning of 16S rRNA genes and later by metagenomics has provided a more realistic view of the microbial community structure of natural environments. Historically, a gap between cultured microbes and the real natural diversity was recognized and expressed as the “great plate anomaly” . We are now in a better position to assess the selective bias introduced by nutrient-rich laboratory media. “Classical” marine bacteria, such as Vibrio species, easily isolated in pure culture, are now recognized as rare inhabitants selected by laboratory conditions . On the other hand, the most abundant marine microbes such as Candidatus Pelagibacter, are extremely difficult to retrieve in pure culture and are still largely studied by direct sequencing of the biomass or, at best, of a few isolated strains . Rather than such unexpected findings remaining restricted to the marine habitat, they rapidly turned out to be the norm, especially in the case of aquatic systems, e.g. the high abundances of low GC Actinobacteria in freshwater  and even in the ocean  and many other groups that are just starting to be known and studied [6, 7]. A similar discovery was made regarding hypersaline environments. The classical extreme halophiles that have been studied for decades such as Halobacterium or Haloarcula were barely detectable by molecular approaches in saturated brines [8, 9]. Instead, slow growing and extremely fastidious to isolate microbes, such as Haloquadratum walsbyi, have shown to be dominant worldwide in neutral saturated brines [10, 11]. Recently, a new group of extremely small archaea (the Nanohaloarchaea) were discovered by metagenomic approaches to be also relatively abundant in saturated brines, but no culture of these microbes is yet available [12–15].
In a large metagenomic study carried out in a solar saltern in which different salinities were studied by direct 454 pyrosequencing , it came as no surprise when in intermediate salinities (19%, about half a way between seawater salinity, 3.5% and saturation, 37%) the classical moderately halophilic bacteria such as Salinivibrio or Halomonas turned out to be nearly absent. The few studies carried out by PCR and cloning of 16S rRNA  had already indicated the presence of bacteria distantly related to Nitrococcus mobilis. Metagenomic assembly provided many contigs that clearly belonged to gamma-proteobacterial cells that gave consistent hits to the genomes of both Nitrococcus mobilis Nb-231 and Alkalilimnicola ehrlichii MLHE-1. Assembly of the metagenome allowed to reconstitute large genomic fragments and indicated that this was a new microbe (or group of microbes) distantly related of the previously mentioned.
In an attempt to retrieve this new microbe, oligotrophic media were designed to isolate bacteria from a ca. 20% salinity saltern ponds in Isla Cristina, Huelva (Spain) and from the large number of colonies screened, a microbe with 95% identity of the 16S rRNA gene to Alkalilimnicola could be isolated. Independently, a sample from the intermediate salinity saltern in Santa Pola also gave a similar isolate. Both microbes, like most environmentally successful aquatic microbes, are extremely slow growing and fastidious to maintain in the laboratory. Their genomes have now been sequenced and assembled and they have proven to be representatives of the bacteria that were dominant in the 19% saltern pond metagenome. Actually they seem to be very abundant at all medium salinity (15-25%) neutral hypersaline waters for which metagenomes are available.
Results and Discussion
Two strains, M19-40 and UAH-SP71, were isolated from saltern ponds of intermediate salinity (see methods) in South west and South east Spain (Isla Cristina and Santa Pola) respectively. Both microbes are representatives of a novel group of Ectothiorhodospiraceae that by metagenomics appeared to be the dominant bacterial group in saltern ponds of intermediate salinity . By 16S rRNA sequencing they were tentatively assigned to a new genus “Spiribacter” and two separate species “Spiribacter salinus” M19-40 and “Spiribacter” sp. UAH-SP71. A complete taxonomic description is being prepared (Leon et al., in preparation). The genomes of the two isolates have now been fully sequenced and assembled into a single contig (Table 1).
General characteristics and comparison with other Ectothiorhodospiraceae
The phylogenomic relationships of “Spiribacter” with other members of the family Ectothiorhodospiraceae that have been completely sequenced are shown in Figure 1. Allochromatium vinosum DSM 180 and Thioflavicoccus mobilis 8321, two purple sulfur marine bacteria which belong to the family Chromatiaceae, were used as outgroup. The family Ectothiorhodospiraceae is composed of purple sulfur photosynthetic bacteria that are in most cases halophilic and also often alkaliphilic . Many of the most classical isolates come from alkaline hypersaline lakes . The results confirm the initial phylogeny described by 16S rRNA gene sequence comparison (Additional file 1: Figure S1). Both strains appear monophyletic and only distantly related to the closest neighbor Arhodomonas aquaeolei DSM 8974, an aerobic chemoheterotroph , and Nitrococcus mobilis Nb-231, an obligate chemolitotrophic bacterium . The average nucleotide identity (ANI) between the two strains was only 77.3%, what fits with what is expected of different species of the same genus . Both microbes have genomes close to 2 Mbp (Table 1), making them the smallest genomes described within the Ecthiorhodospiraceae and for any halophilic bacterium. The genomes are very streamlined with a median intergenic spacer of only 14–19 nucleotides, also the smallest of any member of this family (Figure 2). The GC content is relatively high as seems to be the case with all the Ectothiorhodospiraceae sequenced till now. The genomes were composed of a single circular replicon with only one rRNA operon and with a high level of synteny, considering the relatively low ANI between the two genomes (Figure 3). It is clear from the comparison with the closest available complete genome (Alkalilimnicola ehrlichii MLHE-1) how “Spiribacter” species are simplified in their metabolic versatility, missing the chemolitotrophic and carbon fixation pathways. The small number of IS and other mobile genetic elements and absence of CRISPR system and flagellum was also remarkable. All these are frequently observed characteristics of oligotrophic microbes with streamlined genomes that reach high population densities in aquatic environments, such as Ca. Pelagibacter ubique in the ocean.
In a previous work, metagenomic sequencing indicated that the communities at 19% salinity had already a marked halophilic profile with the hyperhalophilic haloarchaeon Haloquadratum predominating in numbers . However, a major difference with the NaCl saturated ponds was the presence of the gammaproteobacterium associated to Alkalilimnicola-Nitrococcus that we now know corresponds largely (if not completely) to “Spiribacter” representatives. To confirm this, we have analyzed the recruitment of the two genomes from the available metagenomes of hypersaline waters. Besides the 19% metagenome we also used data from a San Diego saltern metagenome  also representing intermediate salinity ponds. As shown in Figure 4, M19-40 was the known bacterial genome recruiting most in San Diego 12-14% and in Santa Pola 19%, only two halophilic archaea, H. walsbyi and Ca. Haloredivivus, recruited more. The data confirms that “Spiribacter” is a very abundant microbe in intermediate salinities but its abundance decreases sharply at both high and low salinities (Figure 4A). In this sense, it can be considered a bona fide moderate halophile that prefers to inhabit intermediate salinities. Former definitions of moderate halophiles  were based on laboratory studies of growth rates at different salinities, but are often misleading as shown clearly by the lack of representation of most moderate halophiles defined this way in intermediate salinities. Besides, the growth salinity range in the laboratory is often very wide, with very similar growth rates over most intermediate salinities. Some of the assembled metagenomic contigs described before were now clearly associated with “Spiribacter” (Additional file 1: Figure S2). Actually, these were the largest contigs obtained from the 19% metagenome confirming the abundance of “Spiribacter” cells in this environment. However, many of these contigs had lower similarities while still being syntenic to “Spiribacter” genomes (Additional file 1: Figure S2) indicating that there might be other “Spiribacter” species present in significant amounts in this specific sample.
The recruitment of M19-40 (Figure 4B) shows the typical metagenomic islands  that are found when a well-represented microbe is compared to the metagenome of its habitat. These islands are a reflection of the diversity of co-existing clonal lineages that diverge in genomic regions delineated by the islands. The most remarkable feature in the case of strain M19-40 is the small number and size of these islands concordant with the highly streamlined genome and relatively restricted metabolic versatility. Genomic island (GI) 1 accumulates 4 of the 7 IS elements of the genome but no specific function could be attributed to the few other genes present. GI2 was the largest and contains a large urease cluster. Ureases are not present in strain UAH-SP71 and, although the cluster in M19-40 recruits very little at high identity, it showed many hits below 97% identity indicating the presence of other lineages that carry ureases. The rest of GI2 contains many toxin-antitoxin systems as is typical of flexible genomic islands of bacteria in general. GI3 is the cluster coding for the O-chain of the lipopolysaccharide which is a typical genomic island in all Gram-negative bacteria and its variability has been attributed to its major role as a phage recognition target [23, 24].
To assess the presence of “Spiribacter” associated microbes worldwide, we have searched the complete Ribosomal Database Project for 16S rRNA sequences highly similar (>97%) to that of either strain. The results described in Additional file 1: Figure S3 indicate that they are found worldwide (as already indicated by the recruitment in Spanish and Californian salterns) and always (when the data were available) in intermediate salinity hypersaline waters.
Both strains of “Spiribacter” were isolated on media for heterotrophic aerobic microbes (M19-40 on pyruvate and yeast extract and UAH-SP71 on pyruvate, dextrose and peptone). The genomes confirm that these microbes are devoid of any photosynthetic or chemolithotrophic capabilities, what sets them apart from most members of the Ectothiorhodospiraceae. No carbon fixing pathways were detected either. The only other member of this family with a purely heterotrophic metabolism is Arhodomonas aquaeolei, a halophilic aerobic heterotroph isolated from a petroleum reservoir . In concordance with their heterotrophic nature, transporters formed a large fraction of the genes in both genomes. TransporterDB (see methods) predicts 217 and 268 transport-associated genes for M19-40 and UAH-SP71, respectively. In comparison, in S. ruber, which has a much larger genome than either of them (3.5 Mbp), a similar number of genes (223) was detected by the same method.
One of the main carbon and energy sources for halophilic microorganisms in the salterns is glycerol. Massive amounts of this compound are produced by the unicellular green algae Dunaliella, the main primary producer in these hypersaline habitats, as a compatible solute to provide osmotic balance . In its catabolic degradation, glycerol is converted to dihydroxyacetone phosphate (DHAP), a glycolysis intermediate, by two different pathways. The first pathway, involving glycerol kinase and glycerol 3-phosphate dehydrogenase , was found as a complete cluster in both “Spiribacter” strains. These clusters included also two genes encoding a glycerol-3-phosphate ABC transporter and a glycerol-3-phosphate regulon repressor. In strain UAH-SP71, three other different glycerol 3-phosphate dehydrogenase genes were also found together in a separate location of the genome. In the second pathway, glycerol is metabolized to DHA first and then phosphorylated by a dihydroxyacetone kinase into DHAP. In strain M19-40 this second route is totally absent while in UAH-SP71 a DHA kinase was found, indicating the ability to use both pathways. Besides, located near this DHA kinase gene, another gene coding for a phosphoenolpyruvate dependent DHA phosphotransferase system was found. Although apparently redundant, this second enzyme seems to be involved in the transport of DHA rather than in catabolism of this compound . The same DHA uptake system has been found before in the haloarchaeon H. walsbyi. The authors of that study suggested that H. walsbyi can use DHA as carbon and energy source, which is known to be released by S. ruber (and likely other halophilic microbes) as a byproduct of glycerol metabolism.
Glycine betaine (trimethylated glycine) is another important osmoprotectant that could be an abundant carbon source in hypersaline environments. Along these lines, strain M19-40 has two paralogous operons of four sarcosine oxidase genes located in separate regions of the genome. Sarcosine oxidase is a heterotetrameric enzyme that catalyzes the oxidative demethylation of sarcosine (N-methylglycine) to yield glycine that can be further catabolized and used as a source of carbon and energy . Strain UAH-SP71 did not have these genes. Another difference indicating different lifestyles for both isolates was the presence of a complete cluster encoding a cytochrome c oxidase (cbb3) complex involved in microaerobic respiration in strain UAH-SP71. This heme-copper oxidase provides a better adaptation to respiration under microoxic conditions  indicating that strain UAH-SP71 might be specialized in comparatively less oxygenated microenvironments.
Regarding nitrogen metabolism, neither “Spiribacter” genomes have genes that code for nitrogen fixation or nitrate/nitrite assimilation. Both genomes have the gene coding for a high-affinity ammonium transporter Amt, indicating that nitrogen uptake occurs in its most reduced form, ammonia. Accordingly, the essential genes involved in ammonia assimilation and amino acid metabolism, glutamine synthetase and glutamate synthase, were found. Besides, several ABC transporters for exogenous nitrogen-rich organic compounds such as basic amino acids, putrescine/spermidine and polyamines were also found. In addition, each genome contained other non-shared genes that can provide nitrogen rich compounds. For example, strain M19-40 contained a complete urease gene cluster (see above), including the genes for the transport and conversion to ammonia, while strain UAH-SP71 carries an ABC transporter for the nitrogen rich osmolyte taurine (TauABC) . Inorganic phosphate appears to be mainly transported in both microbes by the high-affinity phosphate transport system, PstSCAB. Both genomes contained also the ABC transporter PhnCDE, a high affinity uptake system for phosphonates. Phosphonates are organophosphorus compounds characterized by the presence of a carbon- phosphorus bond and could be used also as a nutritional source of phosphorus in response to phosphate starvation .
As expected from the evident streamlining of these genomes, regulatory mechanisms were very scarce. Only two sigma factors involved in response to environmental stresses were found, σ32 a heat shock sigma factor and σ54 that has been linked to the regulation of nitrogen metabolism [33, 34]. Along similar lines, only two two-component regulatory systems were found, one involved on survival under nitrogen limited growth conditions (GlnLG) and PhoBR that regulates the response to variations in the level of extracellular inorganic phosphate [35, 36].
Despite the small genome size, both “Spiribacter” strains appear to have all the typical salt-out osmoregulatory mechanisms. The salt-out strategy is based in keeping most of the inorganic salts out and using organic osmolytes to balance the high salinity of the environment . These “organic compatible solutes” include amino acids and derivatives such as glycine betaine and ectoine . In the case of “Spiribacter”, judging by the number of betaine transporters in these genomes (six different ABC-type glycine betaine transport systems were found in both microbes) it seems that this compound has an important role. Besides, a choline transporter betH was also present. Choline itself is not an osmoprotectant but we found in both genomes the key enzymes for the synthesis of glycine betaine from choline, choline dehydrogenase (betA), glycine betaine aldehyde dehydrogenase (betB) and the transcriptional regulator (betI) . Ectoine is another widely distributed compatible solute commonly synthesized by halophilic bacteria . Ectoine biosynthesis is a separate branch along the pathway for the synthesis of the amino acid aspartate. In both “Spiribacter” we found the complete ectABC gene cluster, diaminobutyric acid (DABA) acetyltransferase (EctA), DABA aminotransferase (EctB) and ectoine synthase (EctC) .
It has been suggested that an important adaptation to saline environments may be an increase in the number of acidic amino acids to reduce their surface hydrophobicity. We have characterized the two strains M19-40 and UAH-SP71 proteomes by virtual 2D-gels and to analyze the differences in protein acidity we have compared them to those from a salt-in halophilic archaeon (Halobacterium sp. NRC-1), a salt in halophilic bacterium (Salinibacter ruber DSM 13855), a marine bacterium (Alteromonas macleodii DE1) and a freshwater bacterium (Fluviicola taffensis RW262T) (Figure 5). The whole proteome pI plots showed the change from the bimodal pattern of the freshwater marine bacterium, F. taffensis, to a single peak around 4.5 for Halobacterium sp. NRC-1. This figure shows how the pI distribution of the second peak representing basic proteins, that contains the transmembrane proteins, has evolved toward acidity in response to increase in salinity of the environment. In a manner similar to Halobacterium sp. NRC-1, both “Spiribacter” showed reduction of the amount of basic amino acids resulting in a single peak at around a pI 5.0. In any case, recent evidence indicates that the low pI peak that characterizes many halophiles is more a condition for a salt-in strategy than a reflection of that mechanism really taking place in the cells [42, 43].
Another way to gauge the salt adaptation of a microbe is the reliance on the Na+ gradient to actively transport nutrients into the cell. Both “Spiribacter” strains have secondary transporters that catalyze the translocation of solutes across the cytoplasmic membrane using electrochemical ion gradients, for example, the tripartite ATP-independent periplasmic family transporters . Thirteen and fifteen of these transporters (in M19-40 and UAH-SP71, respectively) that catalyze uptake of C4-dicarboxylates, sugars and metabolites using the Na+ gradient were found. Both microbes also possess several genes encoding Na+ symporters for proline, bicarbonate or anions and Na+/H+ and Na+/Ca2+ antiporters, consistent with adaptations to a high-salt environment . The genomes contain also a cluster of six genes coding for the multi-subunit Na+/H+ antiporter Mrp , which has been suggested to be the main mechanism to maintain pH homeostasis. In response to osmotic stress, bacteria can also accumulate K+ as an osmoregulatory solute and pH regulator [46, 47]. The uptake of K+ is catalyzed by multiple uptake systems . Interestingly, both “Spiribacter” strains only harbour the gene cluster trkAH that codes for the Trk transport system, while in S. ruber genome four copies of a trkA and two copies of a trkH gene and also the K+ efflux system KefB were present .
XRs are rhodopsin proton pumps that sometimes use a carotenoid pigment antenna to harvest light, a mechanism that is convergent with that of the chlorophylls of photosynthetic microbes and widens the light absorption spectrum . XR was first discovered in the extremely halophilic bacterium S. ruber but later on they have been found to be widespread in marine and halophilic microbes. Recently, a separation of xanthorhodopsins in two families has been proposed on the grounds of the gene sequence and predicted biophysical properties . One of these families, Group II XRs, lack some essential keto-carotenoid binding sites and have been described as being devoid of the carotenoid antenna characteristic of group I . Both “Spiribacter” genomes were found to contain rhodopsin-coding genes that, after sequence comparisons, cluster clearly with the XR genes found in other microbes (Additional file 1: Figure S4). The predicted proteins contained all the characteristic functional groups of proton-pumping rhodopsins (Additional file 1: Figure S5). These results suggest that both XRs are functional. Furthermore, the phylogenetic analysis (Additional file 1: Figure S4) and the characteristic genes found in the cluster (Figure 6) indicate that both closely related genes of “Spiribacter” encode XRs affiliated to subgroup II. This XR subgroup, was found in a heterogeneous group of different alpha and gamma proteobacteria, mostly of marine origins, and in marine eukaryotes. Actually, group II XRs have already been found by other authors in metagenomes in medium salinity, hypersaline and fresh water habitats .
Regardless of whether they have carotenoid antennas or not, all rhodopsins require the chromophore retinal. The gene clusters including the XR gene of group II, always contain four of the genes involved in the synthesis of the retinal precursor β-carotene from isopentenyl diphosphate (crtE, crtI, crtB and crtY) (Figure 6). In addition, UAH-SP71 has a paralogous gene (crtB2, 44% amino acid identity) coding for CrtB (phytoene synthase) located elsewhere in the genome. The presence of two phytoene synthase paralogs has been already reported in the genome of the Gram-positive bacterium Corynebacterium glutamicum. The authors reported that both phytoene synthases were functional leading to an overproduction of carotenoids. Another difference between the two strains involves the location of the β-carotene 15,15'-monooxygenase (Brp/Blh) gene. In M19-40 this gene was found within the cluster while in UAH-SP71 this gene appeared at a distant locus in the genome (Figure 6).
In recent years it has become apparent that planktonic aquatic habitats are often dominated by microbes that have highly streamlined genomes [2, 5]. They tend to be sessile and, in the case of heterotrophs in euphotic waters, often have rhodopsins. In high-salinity hypersaline habitats, two archaeal groups appear to occupy this niche, Haloquadratum and the Nanohaloarchaea. In intermediate salinities there appears to be more diversity and the “Spiribacter” representatives fill a significant gap. The widespread distribution of this microbe and its natural range indicates that it is present in 10 to 25% salinities worldwide. In other words, it is an ecological moderate halophile in the same manner Haloquadratum is an ecological extreme halophile. Thus the classical classification based on growth rate range in the laboratory can be replaced by a more realistic classification based on ecological distribution. The classical microbes classified as moderate halophiles, such as Halomonas or Salinivibrio, have been shown to be almost absent from hypersaline habitats and, at least in the case of Halomonas, appear to be just salt tolerant opportunistic marine bacteria that happen to survive higher salinities and be thus isolated often from moderately saline waters.
“Spiribacter” shows many characteristics in common with dominant aquatic microbes such as Ca. Pelagibacter or Ca. Actinomarina. The cells are small, although not as small as the marine counterparts, and so is the genome. Motility is absent as in most genuinely planktonic prokaryotes. Flagellar motility only has bearing when there are significant gradients of nutrients to traverse. This is unlikely to be the case in the water column of large aquatic habitats such as the ocean. Besides, the genes involved in motility and chemotaxis impose an additional genomic burden. A major difference between the planktonic marine habitat and its hypersaline counterpart (aside from salinity itself) is the nutrient status. Hypersaline waters tend to be eutrophic and organic matter is supposed to be plentiful for heterotrophs. However, the energetic efficiency of streamlined cells such as “Spiribacter” has to be a significant advantage even under relatively nutrient rich conditions. H. walsbyi genome was not considered to be streamlined. It is a relatively large genome (ca. 3.0 Mbp), has very low coding density (76.5%) and contains several IS, repetitive elements  and genomic islands [54, 55]. However, the abundance of IS elements seems to be a general feature of haloarchaea and could be related to their proven level of promiscuity . Furthermore, many regions identified as non-coding have revealed a rich content of small RNAs  that might be very important for the biology of this microbe, i.e. the apparent wastefulness of the Haloquadratum genome might just reflect our ignorance of archaeal biology.
The comparison of both “Spiribacter” with the only other bacterium that seems to be naturally abundant in hypersaline waters, the bacteroidete S. ruber, is also interesting. Both bacteria are carriers of similar xanthorhodopsins, although the type present in “Spiribacter” seems to operate without the associated carotenoid antenna that was originally characterized in this type of rhodopsins . Salinibacter however is motile and also carries sensory rhodopsins and chemotaxis machinery. Actually, the real relevance of S. ruber in the planktonic brines is not very clear. In metagenomes the available genomes do not recruit very extensively and, although its numbers are significant, they are always very far from those of the archaeal members of the population. Salinibacter together with many high GC haloarchaea, that are often isolated from salterns, could live attached to particles, a niche that is often overlooked but can explain the need for motility to commute between particles or get to the most nutrient rich areas.
Isolation of the strains
The two strains used in this study were isolated from water samples collected from two intermediate salinity ponds of the salterns “Braç del Port” located in Santa Pola (Alicante, Spain, 38º13´ N - 0º35´ W) and Isla Cristina (Huelva, Spain, 37°13´ N – 7°19´ W). “Spiribacter” sp. UAH-SP71 isolated from the Santa Pola salterns was grown in medium R2 Agar (Microkit, ref. DMT215) contained (in %): 0.0525 peptone, 0.035 yeast extract, 0.035 dextrose, 0.035 starch, 0.021 K2HPO4, 0.021 sodium pyruvate, 0.0175 tryptone, 0.00168 MgSO4 and supplemented with 10% of sea salts , pH 8. The medium for “Spiribacter salinus” M19-40 contained in % 15 sea salts, 0.5 casein digest, 0.1 yeast extract, 0.1 dextrose and 0.1 pyruvic acid, pH 7.5.
DNA sequencing and assembly
DNA was extracted by phenol-chloroform  and checked for quality on a 1% agarose gel. Whole genome sequencing for both isolates was performed using the Illumina HiSeq 2000 (100-bp paired-end read) sequencing platform. In addition, Pacific Bioscience 3- to 5-Kb reads were generated for “Spiribacter salinus” M19-40. Error correction of the PacBio reads was performed using the Illumina reads as described before . The reads were assembled de novo using MIRA  and Geneious Pro 5.0.1. Both results were compared for equal assemblies. Finally, oligonucleotides designed from the sequence of the ends of assembled contigs were used to obtain single closed contigs .
The genomes were annotated using the NCBI PGAAP annotation pipeline (http://www.ncbi.nlm.nih.gov/genome/annotation_prok/). The predicted protein sequences were also compared using BLASTP to the National Center for Biotechnology Information nr protein database (e value: 1e-5). BioEdit was used to manipulate the sequences . GC content was calculated using the EMBOSS tool geecee . For comparative analyses, reciprocal BLASTN and TBLASTXs searches between the genomes were carried out, leading to the identification of regions of identity, insertions and rearrangements. To allow the interactive visualization of genomic fragment comparisons Artemis v.12 , Artemis Comparison Tool ACTv.9  were used to compare the genomes. Average nucleotide identity (ANI) was calculated as defined before . Transporters were annotated using the TransporterDB database .
To determine the exact phylogenetic relationship of the new isolates within the family, phylogenetic analysis of 16S rRNA gene sequences for all the Ectothiorhodospiraceae members were carried out. Maximum likelihood tree was created using MEGA (version 4.0.2). The tree was rooted using Allochromatium vinosum DSM 180 and Thioflavicoccus mobilis 8321 as outgroups. For creating the phylogenomic tree in Figure 1, the complete genomes of Ectothiorhodospiraceae were analyzed using TIGRfams and 277 proteins conserved in all genomes were identified and concatenated. The concatenated proteins were aligned using Kalign  and a maximum likelihood tree was made using FastTree  using a JTT + CAT model and a gamma approximation.
Genomes recruitments from metagenomic datasets derived from some available marine habitats with different salinity [12, 22] were carried out via BLASTN. A restrictive cut-off of 95% of identity in at least 50 bp was established to guarantee that only similarities at the level of nearly identical microbes were counted. The number of hits was normalized against the genomes and the database sizes. As controls, similar recruitment experiments were carried out for other typical halophilic microorganism in this environment, three haloarchaea: Haloquadratum walsbyi, Halonotius pteroides, Ca. Haloredivivus and four bacteria Salinibacter ruber, Salinivibrio costicola, Chromohalobacter salexigens and Halomonas elongata. Recruitment of the genome of strain M19-40 against the metagenome from Santa Pola (19%) (Figure 4B), 95% of identity in 50% of the length of the metagenomic read was used as a cut-off to construct the plot.
The genomes have been deposited in NCBI GenBank and can be accessed with the following accession numbers: “Spiribacter salinus” M19-40 (CP005963) and “Spiribacter” sp. UAH-SP71 (CP005990).
Staley JT, Konopka A: Measurement of in situ activities of nonphotosynthetic microorganisms in aquatic and terrestrial habitats. Annu Rev Microbiol. 1985, 39: 321-346. 10.1146/annurev.mi.39.100185.001541.
Swan BK, Tupper B, Sczyrba A, Lauro FM, Martinez-Garcia M, González JM, Luo H, Wright JJ, Landry ZC, Hanson NW: Prevalent genome streamlining and latitudinal divergence of planktonic bacteria in the surface ocean. Proc Natl Acad Sci. 2013, 110: 11463-11468. 10.1073/pnas.1304246110.
Grote J, Thrash JC, Huggett MJ, Landry ZC, Carini P, Giovannoni SJ, Rappé MS: Streamlining and core genome conservation among highly divergent members of the SAR11 clade. mBio. 2012, 3: 1434-1437.
Ghai R, McMahon KD, Rodriguez-Valera F: Breaking a paradigm: cosmopolitan and abundant freshwater actinobacteria are low GC. Environ Microbiol Rep. 2012, 4: 29-35. 10.1111/j.1758-2229.2011.00274.x.
Ghai R, Mizuno CM, Picazo A, Camacho A, Rodriguez-Valera F: Metagenomics uncovers a new group of low GC and ultra-small marine actinobacteria. Sci Rep. 2013, 3: 2471-
Albertsen M, Hugenholtz P, Skarshewski A, Nielsen KL, Tyson GW, Nielsen PH: Genome sequences of rare, uncultured bacteria obtained by differential coverage binning of multiple metagenomes. Nat Biotech. 2013, 31: 533-538. 10.1038/nbt.2579.
Rinke C, Schwientek P, Sczyrba A, Ivanova NN, Anderson IJ, Cheng J-F, Darling A, Malfatti S, Swan BK, Gies EA: Insights into the phylogeny and coding potential of microbial dark matter. Nature. 2013, 499: 431-437. 10.1038/nature12352.
Benlloch S, Martínez-Murcia AJ, Rodríguez-Valera F: Sequencing of bacterial and archaeal 16S rRNA genes directly amplified from a hypersaline environment. Syst Appl Microbiol. 1995, 18: 574-581. 10.1016/S0723-2020(11)80418-2.
Benlloch S, López-López A, Casamayor EO, Øvreås L, Goddard V, Daae FL, Smerdon G, Massana R, Joint I, Thingstad F: Prokaryotic genetic diversity throughout the salinity gradient of a coastal solar saltern. Environ Microbiol. 2002, 4: 349-360. 10.1046/j.1462-2920.2002.00306.x.
Antón J, Llobet-Brossa E, Rodríguez-Valera F, Amann R: Fluorescence in situ hybridization analysis of the prokaryotic community inhabiting crystallizer ponds. Environ Microbiol. 1999, 1: 517-523. 10.1046/j.1462-2920.1999.00065.x.
Bolhuis H, Poele EM, Rodriguez-Valera F: Isolation and cultivation of Walsby's square archaeon. Environ Microbiol. 2004, 6: 1287-1291. 10.1111/j.1462-2920.2004.00692.x.
Ghai R, Pašić L, Fernández AB, Martin-Cuadrado A-B, Mizuno CM, McMahon KD, Papke RT, Stepanauskas R, Rodriguez-Brito B, Rohwer F: New abundant microbial groups in aquatic hypersaline environments. Sci Rep. 2011, 1: 135-
Narasingarao P, Podell S, Ugalde JA, Brochier-Armanet C, Emerson JB, Brocks JJ, Heidelberg KB, Banfield JF, Allen EE: De novo metagenomic assembly reveals abundant novel major lineage of Archaea in hypersaline microbial communities. ISME J. 2012, 6: 81-93. 10.1038/ismej.2011.78.
Podell S, Ugalde JA, Narasingarao P, Banfield JF, Heidelberg KB, Allen EE: Assembly-driven community genomics of a hypersaline microbial ecosystem. PLoS ONE. 2013, 8: e61692-10.1371/journal.pone.0061692.
Casanueva A, Galada N, Baker G, Grant W, Heaphy S, Jones B, Yanhe M, Ventosa A, Blamey J, Cowan D: Nanoarchaeal 16S rRNA gene sequences are widely dispersed in hyperthermophilic and mesophilic halophilic environments. Extremophiles. 2008, 12: 651-656. 10.1007/s00792-008-0170-x.
Ventosa A, Nieto JJ, Oren A: Biology of moderately halophilic aerobic bacteria. Microbiol Mol Biol Rev. 1998, 62: 504-544.
Tourova TP, Spiridonova EM, Berg IA, Slobodova NV, Boulygina ES, Sorokin DY: Phylogeny and evolution of the family Ectothiorhodospiraceae based on comparison of 16S rRNA, cbbL and nifH gene sequences. Int J Syst Evol Microbiol. 2007, 57: 2387-2398. 10.1099/ijs.0.65041-0.
Oremland RS, Hoeft SE, Santini JM, Bano N, Hollibaugh RA, Hollibaugh JT: Anaerobic oxidation of arsenite in mono lake water and by a facultative, arsenite-oxidizing chemoautotroph, strain MLHE-1. Appl Environ Microbiol. 2002, 68: 4795-4802. 10.1128/AEM.68.10.4795-4802.2002.
Adkins JP, Madigan MT, Mandelco L, Woese CR, Tanner RS: Arhodomonas aquaeolei gen. nov., sp. nov., an aerobic, halophilic bacterium isolated from a Subterranean Brine. Int J Syst Bacteriol. 1993, 43: 514-520. 10.1099/00207713-43-3-514.
Watson S, Waterbury J: Characteristics of two marine nitrite oxidizing bacteria, Nitrospina gracilis nov. gen. nov. sp. and Nitrococcus mobilis nov. gen. nov. sp. Archiv Mikrobiol. 1971, 77: 203-230. 10.1007/BF00408114.
Konstantinidis KT, Tiedje JM: Genomic insights that advance the species definition for prokaryotes. Proc Natl Acad Sci U S A. 2005, 102: 2567-2572. 10.1073/pnas.0409727102.
Rodriguez-Brito B, Li L, Wegley L, Furlan M, Angly F, Breitbart M, Buchanan J, Desnues C, Dinsdale E, Edwards R: Viral and microbial community dynamics in four aquatic environments. ISME J. 2010, 4: 739-751. 10.1038/ismej.2010.1.
Rodriguez-Valera F, Martin-Cuadrado A-B, Rodriguez-Brito B, Pasic L, Thingstad TF, Rohwer F, Mira A: Explaining microbial population genomics through phage predation. Nat Rev Micro. 2009, 7: 828-836. 10.1038/nrmicro2235.
Rodriguez-Valera F, Ussery DW: Is the pan-genome also a pan-selectome?. F1000Research. 2012, 1: 16-
Elevi Bardavid R, Khristo P, Oren A: Interrelationships between Dunaliella and halophilic prokaryotes in saltern crystallizer ponds. Extremophiles. 2008, 12: 5-14. 10.1007/s00792-006-0053-y.
Sherwood KE, Cano DJ, Maupin-Furlow JA: Glycerol-mediated repression of glucose metabolism and glycerol kinase as the sole route of glycerol catabolism in the haloarchaeon Haloferax volcanii. J Bacteriol. 2009, 191: 4307-4315. 10.1128/JB.00131-09.
Kotrba P, Inui M, Yukawa H: Bacterial phosphotransferase system (PTS) in carbohydrate uptake and control of carbon metabolism. J Biosci Bioeng. 2001, 92: 502-517.
Bolhuis H, Palm P, Wende A, Falb M, Rampp M, Rodriguez-Valera F, Pfeiffer F, Oesterhelt D: The genome of the square archaeon Haloquadratum walsbyi: life at the limits of water activity. BMC Genomics. 2006, 7: 169-10.1186/1471-2164-7-169.
Kvalnes-Krick K, Jorns MS: Bacterial sarcosine oxidase: comparison of two multisubunit enzymes containing both covalent and noncovalent flavin. Biochemistry. 1986, 25: 6061-6069. 10.1021/bi00368a034.
Cosseau C, Batut J: Genomics of the ccoNOQP-encoded cbb 3 oxidase complex in bacteria. Arch Microbiol. 2004, 181: 89-96. 10.1007/s00203-003-0641-5.
Denger K, Smits THM, Cook AM: Genome-enabled analysis of the utilization of taurine as sole source of carbon or of nitrogen by Rhodobacter sphaeroides 2.4.1. Microbiology. 2006, 152: 3197-3206. 10.1099/mic.0.29195-0.
Villarreal-Chiu JF, Quinn JP, McGrath JW: The genes and enzymes of phosphonate metabolism by bacteria, and their distribution in the marine environment. Front Microbiol. 2012, 3: 19-
Yura T, Nagai H, Mori H: Regulation of the heat-shock response in bacteria. Annu Rev Microbiol. 1993, 47: 321-350. 10.1146/annurev.mi.47.100193.001541.
Albright LM, Huala E, Ausubel FM: Prokaryotic signal transduction mediated by sensor and regulator protein pairs. Annu Rev Genet. 1989, 23: 311-336. 10.1146/annurev.ge.23.120189.001523.
Wanner BL: Phosphorus assimilation and control of the phosphate regulon. Escherichia coli and Salmonella typhimurium: cellular and molecular biology. Edited by: Neidhardt FC, Curtiss RIII, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schaechter M, Umbarger HE. 1996, Washington, D. C: ASM Press, 1357-1381.
Leigh JA, Dodsworth JA: Nitrogen regulation in bacteria and archaea. Annu Rev Microbiol. 2007, 61: 349-377. 10.1146/annurev.micro.61.080706.093409.
Empadinhas N, da Costa M: Osmoadaptation mechanisms in prokaryotes: distribution of compatible solutes. Int Microbiol. 2008, 11: 151-161.
Saum S, Muller V: Regulation of osmoadaptation in the moderate halophile Halobacillus halophilus: chloride, glutamate and switching osmolyte strategies. Saline Systems. 2008, 4: 4-10.1186/1746-1448-4-4.
Cánovas D, Vargas C, Kneip S, Morón Ma J, Ventosa A, Bremer E, Nieto JJ: Genes for the synthesis of the osmoprotectant glycine betaine from choline in the moderately halophilic bacterium Halomonas elongata DSM 3043. Microbiology. 2000, 146: 455-463.
Reshetnikov A, Khmelenina V, Trotsenko Y: Characterization of the ectoine biosynthesis genes of haloalkalotolerant obligate methanotroph “Methylomicrobium alcaliphilum 20Z”. Arch Microbiol. 2006, 184: 286-297. 10.1007/s00203-005-0042-z.
Cánovas D, Vargas C, Calderón MI, Ventosa A, Nieto JJ: Characterization of the genes for the biosynthesis of the compatible solute ectoine in the moderately halophilic bacterium Halomonas elongata DSM 3043. Syst Appl Microbiol. 1998, 21: 487-497. 10.1016/S0723-2020(98)80060-X.
Elevi Bardavid R, Oren A: Acid-shifted isoelectric point profiles of the proteins in a hypersaline microbial mat: an adaptation to life at high salt concentrations?. Extremophiles. 2012, 16: 787-792. 10.1007/s00792-012-0476-6.
Deole R, Challacombe J, Raiford DW, Hoff WD: An extremely halophilic proteobacterium combines a highly acidic proteome with a low cytoplasmic potassium content. J Biol Chem. 2013, 288: 581-588. 10.1074/jbc.M112.420505.
Mulligan C, Fischer M, Thomas GH: Tripartite ATP-independent periplasmic (TRAP) transporters in bacteria and archaea. FEMS Microbiol Rev. 2011, 35: 68-86. 10.1111/j.1574-6976.2010.00236.x.
Swartz T, Ikewada S, Ishikawa O, Ito M, Krulwich T: The Mrp system: a giant among monovalent cation/proton antiporters?. Extremophiles. 2005, 9: 345-354. 10.1007/s00792-005-0451-6.
Brooth IR: Regulation of cytoplasmic pH in bacteria. Microbiol Rev. 1985, 49: 359-378.
Epstein W: Osmoregulation by potassium transport in Escherichia coli. FEMS Microbiol Lett. 1986, 39: 73-78. 10.1111/j.1574-6968.1986.tb01845.x.
Corratgé-Faillie C, Jabnoune M, Zimmermann S, Véry AA, Fizames C, Sentenac H: Potassium and sodium transport in non-animal cells: the Trk/Ktr/HKT transporter family. Cell Mol Life Sci. 2010, 67: 2511-2532. 10.1007/s00018-010-0317-7.
Mongodin EF, Nelson KE, Daugherty S, DeBoy RT, Wister J, Khouri H, Weidman J, Walsh DA, Papke RT, Sanchez Perez G: The genome of Salinibacter ruber: convergence and gene exchange among hyperhalophilic bacteria and archaea. Proc Natl Acad Sci USA. 2005, 102: 18147-18152. 10.1073/pnas.0509073102.
Spudich JL, Jung K-H: Microbial rhodopsins: phylogenetic and functional diversity. Handbook of Photosensory Receptors. 2005, Weinheim, FRG: Wiley-VCH Verlag GmbH & Co. KGaA, 1-23.
Vollmers J, Voget S, Dietrich S, Gollnow K, Smits M, Meyer K, Brinkhoff T, Simon M, Daniel R: Poles apart: arctic and antarctic Octadecabacter strains share high genome plasticity and a New type of xanthorhodopsin. PLoS ONE. 2013, 8: e63422-10.1371/journal.pone.0063422.
Heider S, Peters-Wendisch P, Wendisch V: Carotenoid biosynthesis and overproduction in Corynebacterium glutamicum. BMC Microbiol. 2012, 12: 198-10.1186/1471-2180-12-198.
Chen Y, Zhou F, Li G, Xu Y: MUST: a system for identification of miniature inverted-repeat transposable elements and applications to Anabaena variabilis and Haloquadratum walsbyi. Gene. 2009, 436: 1-7. 10.1016/j.gene.2009.01.019.
Cuadros-Orellana S, Martin-Cuadrado A-B, Legault B, D'Auria G, Zhaxybayeva O, Papke RT, Rodriguez-Valera F: Genomic plasticity in prokaryotes: the case of the square haloarchaeon. ISME J. 2007, 1: 235-245. 10.1038/ismej.2007.35.
Dyall-Smith ML, Pfeiffer F, Klee K, Palm P, Gross K, Schuster SC, Rampp M, Oesterhelt D: Haloquadratum walsbyi: limited diversity in a global pond. PLoS ONE. 2011, 6: e20968-10.1371/journal.pone.0020968.
Nelson-Sathi S, Dagan T, Landan G, Janssen A, Steel M, McInerney JO, Deppenmeier U, Martin WF: Acquisition of 1,000 eubacterial genes physiologically transformed a methanogen at the origin of Haloarchaea. Proc Natl Acad Sci. 2012, 109: 20537-20542. 10.1073/pnas.1209119109.
Rodriguez-Valera F, Ruiz-Berraquero F, Ramos-Cormenzana A: Isolation of extreme halophiles from seawater. Appl Environ Microbiol. 1979, 38: 164-165.
Neumann B, Pospiech A, Schairer HU: Rapid isolation of genomic DNA from Gram-negative bacteria. Trends Genet. 1992, 8: 332-333.
Koren S, Schatz MC, Walenz BP, Martin J, Howard JT, Ganapathy G, Wang Z, Rasko DA, McCombie WR, Jarvis ED, Phillippy AM: Hybrid error correction and de novo assembly of single-molecule sequencing reads. Nat Biotech. 2012, 30: 693-700. 10.1038/nbt.2280.
Chevreux B: MIRA: an automated genome and EST assembler. 2005, Heidelberg, Germany: Ph. D. thesis. University R-K, 45-46.
Leon MJ, Ghai R, Fernandez AB, Sanchez-Porro C, Rodriguez-Valera F, Ventosa A: Draft genome of Spiribacter salinus M19-40, an abundant gammaproteobacterium in aquatic hypersaline environments. Genome Announcements. 2013, 1: 1-
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.
Rice P, Longden I, Bleasby A: EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.
Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation. Bioinformatics. 2000, 16: 944-945. 10.1093/bioinformatics/16.10.944.
Carver TJ, Rutherford KM, Berriman M, Rajandream M-A, Barrell BG, Parkhill J: ACT: the Artemis comparison tool. Bioinformatics. 2005, 21: 3422-3423. 10.1093/bioinformatics/bti553.
Ren Q, Chen K, Paulsen IT: TransportDB: a comprehensive database resource for cytoplasmic membrane transport systems and outer membrane channels. Nucleic Acids Res. 2007, 35: D274-D279. 10.1093/nar/gkl925.
Lassmann T, Sonnhammer E: Kalign - an accurate and fast multiple sequence alignment algorithm. BMC Bioinforma. 2005, 6: 298-10.1186/1471-2105-6-298.
Price MN, Dehal PS, Arkin AP: FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010, 5: e9490-10.1371/journal.pone.0009490.
Imasheva ES, Balashov SP, Choi AR, Jung K-H, Lanyi JK: Reconstitution of Gloeobacter violaceus rhodopsin with a light-harvesting carotenoid antenna. Biochemistry. 2009, 48: 10948-10955. 10.1021/bi901552x.
This work was supported by projects MAGYK (BIO2008-02444), MICROGEN (Programa CONSOLIDER-INGENIO 2010 CDS2009-00006), CGL2009-12651-C02-01, CGL2010-19303, CSD2009-00006 and BIO2011-12879-E from the Spanish Ministerio de Ciencia e Innovación, DIMEGEN (PROMETEO/2010/089), ACOMP/2009/155 from the Generalitat Valenciana, P10-CVI-6226 from the Junta de Andalucía, MaCuMBA from European Community (EC) (Ref. FP7-KBBE-2012-6-311975) and by FEDER funds. RG was supported by a Juan de la Cierva scholarship from the Spanish Ministerio de Ciencia e Innovación.
The authors declare that they have no competing interests.
FRV, AV and JS conceived the study, and participated in its design and coordination. ML, ARO, JLCP and CSP isolated the strains. Assembly and genome analysis was carried out by MLP, RG and MJL. FRV, MLP and RG wrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional 1: 16S rRNA phylogeny. Maximum likelihood phylogenetic tree based on the comparison of 16S rRNA gene sequences of the Ectothiorhodospiraceae. Allochromatium vinosum DSM 180 and Thioflavicoccus mobilis 8320, belonging to the Chromatiaceae were used as outgroup and are shown in red. Bootstrap values are indicated at the nodes. Figure S2. BLASTN comparisons of metagenomic contigs from 19% Santa Pola to “Spiribacter” genomes. The metagenomic contigs are shown in the middle. A color key for the similarity is shown on the top right. Figure S3. Global distribution of Spiribacter 16S rRNA gene sequences. Locations where 16S rRNA gene sequences from the Ribosomal Database Project were found (>97% identical, >300 bp) are indicated by colored boxes. The color code indicates the number of sequences found at each location. The map is a modified version of a freely available map from http://www.naturalearthdata.com. Figure S4. Phylogenetic tree of the two xanthorhodopsins found in both “Spiribacter” with all the xanthorhodopsins available. Taxonomy and origin of isolation of each strain are also shown. Figure S5. Xanthorhodopsin amino acid sequence alignment. Multiple alignments of all the predicted aminoacid sequences of the two xanthorhodopsin subgroups. Rectangles over the sequence indicate predicted transmembrane regions. Proton acceptor and donor and the conserved lysine linked to the cofactor retinal are marked by a rectangle with a yellow line. Yellow rectangles with red line indicated the residues that interact with the keto-carotenoid identified by . Maintaining nomenclature, the letters c, g, k and r, indicate contact with the chain, glucoside, keto group and ring of the carotenoid, respectively. (PDF 355 KB)
Authors’ original submitted files for images
About this article
Cite this article
López-Pérez, M., Ghai, R., Leon, M.J. et al. Genomes of “Spiribacter”, a streamlined, successful halophilic bacterium. BMC Genomics 14, 787 (2013). https://doi.org/10.1186/1471-2164-14-787
- Halophilic bacteria
- Moderate halophile