ORIGINAL RESEARCH article
Cold-Active Shewanella glacialimarina TZS-4T nov. Features a Temperature-Dependent Fatty Acid Profile and Putative Sialic Acid Metabolism
- 1RNAcious Laboratory, Molecular and Integrative Biosciences Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland
- 2Doctoral Programme in Microbiology and Biotechnology, University of Helsinki, Helsinki, Finland
- 3Molecular and Integrative Biosciences Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland
- 4Helsinki University Lipidomics Unit HiLIPID, Helsinki Institute of Life Science HiLIFE and Biocenter Finland, Helsinki, Finland
- 5The Laboratory of Structural Biology, Helsinki Institute of Life Science HiLIFE, Helsinki, Finland
Species of genus Shewanella are among the most frequently identified psychrotrophic bacteria. Here, we have studied the cellular properties, growth dynamics, and stress conditions of cold-active Shewanella strain #4, which was previously isolated from Baltic Sea ice. The cells are rod-shaped of ~2μm in length and 0.5μm in diameter, and they grow between 0 and 25°C, with an optimum at 15°C. The bacterium grows at a wide range of conditions, including 0.5–5.5% w/v NaCl (optimum 0.5–2% w/v NaCl), pH 5.5–10 (optimum pH 7.0), and up to 1mM hydrogen peroxide. In keeping with its adaptation to cold habitats, some polyunsaturated fatty acids, such as stearidonic acid (18:4n-3), eicosatetraenoic acid (20:4n-3), and eicosapentaenoic acid (20:5n-3), are produced at a higher level at low temperature. The genome is 4,456kb in size and has a GC content of 41.12%. Uniquely, strain #4 possesses genes for sialic acid metabolism and utilizes N-acetyl neuraminic acid as a carbon source. Interestingly, it also encodes for cytochrome c3 genes, which are known to facilitate environmental adaptation, including elevated temperatures and exposure to UV radiation. Phylogenetic analysis based on a consensus sequence of the seven 16S rRNA genes indicated that strain #4 belongs to genus Shewanella, closely associated with Shewanella aestuarii with a ~97% similarity, but with a low DNA–DNA hybridization (DDH) level of ~21%. However, average nucleotide identity (ANI) analysis defines strain #4 as a separate Shewanella species (ANI score=76). Further phylogenetic analysis based on the 92 most conserved genes places Shewanella strain #4 into a distinct phylogenetic clade with other cold-active marine Shewanella species. Considering the phylogenetic, phenotypic, and molecular characterization, we conclude that Shewanella strain #4 is a novel species and name it Shewanella glacialimarina sp. nov. TZS-4T, where glacialimarina means sea ice. Consequently, S. glacialimarina TZS-4T constitutes a promising model for studying transcriptional and translational regulation of cold-active metabolism.
The genus Shewanella belongs to the order Alteromonadales of class Gammaproteobacteria. In 1985, MacDonell and Colwell (1985) first characterized members of this group comprising Gram stain negative, rod-shaped, and facultative anaerobic bacteria, and in 2004, they were classified into family Shewanellaceae (Ivanova et al., 2004). Members of this family have been isolated from diverse aquatic and marine habitats, including deep sea, ocean sediments, freshwater, and sea ice (Li et al., 2014; Luhtanen et al., 2014; Kim et al., 2016; Wang and Sun, 2016). Shewanella species are known to be active in a broad temperature range, comprising both psychrotrophic (Shewanella baltica, Shewanella frigidimarina, and Shewanella polaris) and mesophilic (Shewanella oneidensis, Shewanella algae, and Shewanella putrefaciens) bacteria (Rossello-Mora et al., 1995; Khashe and Janda, 1998; Venkateswaran et al., 1999; Bozal et al., 2002; Vogel et al., 2005; Cha et al., 2020). Shewanella bacteria are known for their versatile metabolic pathways and wide range of electron acceptors, including various oxidized metals, such as Mn(III), Mn(IV), Fe(III), U(VI) (Myers and Nealson, 1988; Nealson and Saffarini, 1994; Fredrickson et al., 2000), as well as nitrate, sulfite, thiosulfate, and elemental sulfur (Fredrickson et al., 2008), which are vital for bioremediation and biogeochemical cycles.
Shewanella spp. are saprophytic and form an integral part of the marine microflora. Certain Shewanella species, such as S. putrefaciens and S. baltica, have been associated with spoilage of fishery food products due to their psychrotrophic nature, which permits growth at low temperatures (Jørgensen and Huss, 1989; Vogel et al., 2005). These bacteria reduce trimethyl-amine-N-oxide (TMAO) to trimethylamine (TMA), which generates a pungent odor that alongside hydrogen sulfide (H2S) gas, produced by degradation of amino acids, further exacerbates the spoilage process (Gram et al., 1987). Recent emergence of antibiotic resistance and auxiliary metabolic pathways, such as protein and lipid degradation, may further benefit some food spoilage species of Shewanella, resulting in additional cold storage problems.
Shewanella strain #4 was isolated from Baltic Sea ice around Tvärminne Zoological Station in Hanko, Finland (Luhtanen et al., 2014). Phylogenetic analysis of Shewanella strain #4 revealed a distinct clade of cold-active Shewanella species, including S. frigidimarina, S. polaris, and S. arctica (Bozal et al., 2002; Kim et al., 2012; Cha et al., 2020). Its ability to thrive in cold, icy conditions and its association with marine life makes Shewanella strain #4, a suitable candidate to further study the metabolic and gene expression changes in bacteria adapted to cold environments. For example, we located a complete gene cluster associated with sialic acid metabolism and showed that strain #4 can metabolize sialic acid as a carbon source – a characteristic that is generally found in Vibrio choleora, Yersinia pestis, Clostridium perfringens, etc., and regarded as advantageous in commensal and pathogenic bacteria. Based on our biochemical, physiological, and phylogenetic results, we propose Shewanella strain #4 as a novel Shewanella species named Shewanella glacialimarina TZS-4T.
Materials and Methods
Bacterial Strains and Isolation
The bacterial strain #4 (hereafter referred to as S. glacialimarina TZS-4T) was previously isolated from a Baltic Sea ice sample (Luhtanen et al., 2014). Two other Shewanella strains, S. frigidimarina ACAM 591 and S. baltica LMG 2250, were purchased from the Leibniz Institute DSMZ-German collection of microorganisms under catalog number DSM-12253 and DSM-9439, respectively. All strains were grown on a solid medium containing rich 25% w/v marine broth [abbreviated as rMB; 7.5g peptone (Sigma-Aldrich), 1.5g yeast extract (Fisher Bioreagents), and 9.35g marine broth (BD-Difco) in 1L of ddH2O (Milli-Q, Merck)] agar (15g/L) and stored in 30% v/v glycerol stock at −80°C.
General Characterization of Bacteria
A glass slide fixed with bacterial cells was first stained with crystal violet dye solution (Sigma) (2% w/v crystal violet dye in 95% v/v ethanol) for Gram staining. These slides were flooded with Gram’s iodine [1g iodine (Riedel-de Haen) and 2g potassium iodide (Merck) in 300ml distilled water] and counterstained with safranin (Merck) (2.5% w/v in 95% v/v ethanol; Beveridge, 2001). Samples were observed under a microscope at 100× and 1,000× magnification.
The cells were cultured on rMB agar. After 48h of growth, a single colony of S. glacialimarina TZS-4T was transferred to 50ml of rMB and grown for 48h at 15°C with constant aeration at 200rpm. Fresh rMB was inoculated with the starter culture to an optical density at 600nm (OD600) of 0.2, and the cells were grown until the desired OD600 was reached. The turbidity was measured using an Eppendorf Biophotometer.
Hemolysis, Motility, and Hydrogen Sulfide Production Test
Blood agar plates (3% v/v sheep blood with agar base no.2) (Labema) were used to check hemolytic activity. Iron agar test was performed to determine bacterial motility and H2S production (Gram et al., 1987). Both blood agar and iron agar tubes were inoculated with the starter culture and incubated at 15°C for 72h.
Catalase and Oxidase Test
Catalase enzyme activity was confirmed with 3% v/v hydrogen peroxide (H2O2) (Fisher Chemical) (Taylor and Achanzar, 1972), and oxidase reagent (1% w/v tetramethyl-p-phenylenediamine dihydrochloride) (Acros Organics) was used for the oxidase test (Kovacs, 1956). Both catalase and oxidase test were done according to established protocols.
Temperature, pH, and Salinity Conditions
The cells were prepared in rMB media as described in section Growth Conditions, except for the starter culture to assess growth at 0°C, which was grown at 4°C for 48h (necessary for growth at 0°C). Various growth conditions, namely, temperature (0, 4, 15, and 25°C), pH (4.5–10.5), salinity (5–55g/L NaCl), and hydrogen peroxide (0.5–4mM)-induced oxidative stress, were evaluated as a time-course series by determining the change in OD600 values.
Carbon Assimilation Using Minimal Growth Media
Carbohydrate assimilation was determined by growing S. glacialimarina TZS-4T on M-9 minimal media [64.0g Na2HPO4, 15.0g KH2PO4 (Acros Organics), 2.5g NaCl (Fisher Chemical), and 5.0g NH4Cl (Fisher Chemical) in 1L ddH2O (Milli-Q) to make 5× stock] in the presence of a sole carbon source, namely, 200mM glucose (VWR Chemicals), fructose (Fisher Chemical), maltose (Alfa Aesar), or galactose (Fluka), respectively. The cells were grown, as described in section Growth Conditions, in 10ml of M-9 minimal media using starter cultures grown in rMB. Bacterial growth was monitored by measuring OD600. M-9 media without any carbon source was used as negative control.
Ammonia production was analyzed by colorimetric assays using Nessler’s reagent (VWR Chemicals). The cells were grown in 50ml rMB for 24 or 48h, as described in section Growth Conditions. Cells were removed by centrifugation for 5min at 3200g, 4°C, after which 10ml of supernatant was collected. Subsequently, 1ml of Nessler’s reagent was added to the supernatant and the formation of a dark blue color indicating the presence of ammonia.
Phenotypic Characterization Using Biolog GENIII Plates
Biolog GENIII microplates were used for initial phenotypic characterization, providing a standardized platform to test 94 different phenotypic characteristics, including 71 carbon utilization assays and 23 chemical sensitivity tests. Samples were prepared and tested according to the manufacturer’s instructions. All assays were performed in triplicates.
The bacteria were grown in both rMB media and rMB agar in the presence of different antibiotics, such as ampicillin (100μg/ml in ddH2O), kanamycin (50μg/ml in ddH2O), tetracycline (10μg/ml in ddH2O), and chloramphenicol (25μg/ml in 99.5% EtOH). The growth on antibiotic media was observed after 24h, and results were interpreted as resistant (+) or susceptible (−).
Genome Sequencing and Phylogenetic Analysis of Shewanella
The genomic sequence of S. glacialimarina TZS-4T was determined using the Pacific Biosciences PacBio Sequel sequencing technology at the DNA Sequencing and Genomics Laboratory (Helsinki Institute of Life Science, University of Helsinki). Phylogenetic trees based on maximum likelihood, neighbor-joining, and maximum parsimony methods using the consensus sequence of the seven gene copies of 16S rRNA were constructed on MEGA software with 100 replicates as bootstrap value. To infer genome-based phylogeny, we constructed a concatenated bacterial core-gene set phylogenetic tree by the up-to-date bacterial core-gene set (UBCG) method (Na et al., 2018). A model phylogenetic tree was depicted for species delineation based on average nucleotide identities (ANI) score using the fastANI (ver. 1.31) and OrthoANI tools (ver. 0.90; Lee et al., 2016; Jain et al., 2018). OrthoANI calculated the ANI score of query organism against a group of reference genomes and generated a distance matrix. The phylogenetic tree with its heat-map shows the relationship among the genome sets. We used AAI-profiler1 (Medlar et al., 2018) as an additional method to predict the taxonomic identity of the S. glacialimarina TZS-4T. Open reading frames were predicted using GenMark.hmm prokaryotic version 3.25 (Besemer, 2001) applying the same search parameters as used for S. frigidimarina NCIMB400. The genome sequence has been submitted to GenBank under the accession number CP041216.
Transmission Electron Microscopy
The size and shape of S. glacialimarina TZS-4T were determined by transmission electron microscopy (TEM). The cells were grown in rMB media at 15°C until they reached an OD600 of 0.6, which corresponds to approximately 108 cells/ml. Later, 2μl of the culture was incubated on carbon-coated Cu mesh grids for 1min and then negative stained with filtered, neutral 2% uranyl acetate for 15s. TEM examination was conducted using a JEOL JEM-1400 microscope (Jeol Ltd., Tokyo, Japan) operating at 80kV and equipped with Gatan Orius SC 1000B bottom-mounted CCD camera (Gatan Inc., United States).
Scanning Electron Microscopy
For scanning electron microscopy (SEM), cells were grown to OD600 of 0.6, as described in section Growth Conditions. Cells were harvested from 10ml of culture by centrifugation (5min, 3,200g, 4°C) and resuspended in the original volume with phosphate buffer saline (PBS). Cells were fixed with 2.5% v/v glutaraldehyde at 4°C for 20h. After fixation, the cells were washed twice with PBS using centrifugation and resuspended in half of the original volume of PBS. Cell suspension (100μl) was spread on concanavalin A-coated glasses, fixed with osmium tetroxide, and dehydrated through ethanol series (50, 70, and 96%). Dehydrated cells were incubated in hexamethyldisilazane overnight and coated with platinum. Finally, SEM examination was conducted using FEI Quanta 250 Field Emission Gun (FEG) Scanning Electron Microscope.
Fatty Acid Composition
Starter culture of S. glacialimarina TZS-4T was used to inoculate a 50ml culture as described in section Growth Conditions and grown until it reached OD600 of 0.8. Samples were prepared at different temperature conditions, i.e., 5, 15, and 25°C, and the pellet from 5ml culture was washed five times with PBS. Finally, the washed pellet was used for the analysis of fatty acid (FA) composition by gas chromatography. The bacterial pellets were subjected to transmethylation by heating with 1% v/v H2SO4 in methanol under nitrogen, and the formed FA methyl esters were extracted with hexane (Christie, 1993). The samples were dried with anhydrous Na2SO4 and concentrated. The FA methyl esters were identified based on their mass spectra recorded by GCMS-QP2010 Ultra (Shimadzu Scientific Instruments, Kyoto, Japan) and compared to the spectra of several authentic standard mixtures (including Supelco 47,080-U Bacterial Acid Methyl Ester BAME Mix) and published reference mass spectra (Christie, 2019). Quantitative composition was determined by using a Shimadzu GC-2010 Plus with flame ionization detector. The responses were corrected by using the theoretical correction factors for this detector (Ackman, 2007). Both gas chromatographs were equipped with a Zebron ZB-wax capillary column (30m, 0.25mm ID and film thickness 0.25μm; Phenomenex, Torrance, CA, United States). The FA compositions are expressed as mol% profiles, and the FAs are marked using the following abbreviations: (carbon number):(number of double bonds) n-(position of the first double bond calculated from the methyl end) (e.g., 16:1n-7 for palmitoleate). The identified iso and anteiso branches and the 3-OH substituents of the FAs were indicated. The FA components exceeding 0.1mol% at least in one sample were listed, and the remaining were summed as trace FAs. The fold change analysis of unsaturated FA relative concentrations was performed w.r.t 15°C samples as reference. SD was calculated using propagation of error.
RNA isolation was performed using acidic phenol:bromochloropropane (BCP; Chomczynski and Mackey, 1995; Gregorova et al., 2020). The cells were grown from the starter culture in 50ml of rMB until it reached OD600 of 0.8. The cells were harvested at 3200g, 4°C for 10min and stored in −80°C. The pellet was resuspended in 4ml of 0.9% w/v NaCl solution followed by 4ml of acidic phenol (Sigma) and 800μl of BCP (Acros organic). Glass beads were added to break the bacterial cells and vortexed for 10min at room temperature (RT). The lysate was centrifuged at 10000g for 10min at RT, and aqueous phase was collected. The aqueous phase was re-extracted twice with phenol:BCP (2ml of acidic phenol and 400μl of BCP). Total RNA from the aqueous phase was precipitated by adding 2.5 vol of 99.6% v/v EtOH, at −20°C overnight and pelleted by centrifugation at 10000g for 20min at 4°C. The RNA pellets were air-dried and resuspended in RNAase-free ddH2O water. RNA concentration was measured using a Nanodrop 2000 spectrophotometer (Thermo Scientific) and run on 2% w/v agarose, Tris-borate EDTA (TBE) gel with of Midori green (4μl per 100ml of 2% agarose) for quality assessment. The images were captured on Gel-doc XR (Biorad).
Reverse Transcription-Quantitative PCR
Primers were designed for each gene using the IDT PrimerQuest tool2 (Supplementary Table ST-1). The purified RNA (30μg) was DNAase treated with 7U of RQ1 RNAase-Free DNase (Promega). cDNA was synthesized using 200U of Maxima reverse transcriptase primed with 1μl random hexamers (0.2μg/μl) (Thermo Scientific) and 3μg of DNase-treated RNA. Both DNAase treatment and reverse transcription reactions were performed according to the manufacturers’ instructions. cDNA synthesis was validated by PCR (using primers listed in Supplementary Table ST-1), and the products were analyzed on a 2% agarose gel (Supplementary Figure SF-1). Quantitative PCR reactions were performed using Perfecta SYBR green FastMix, Low ROX (Quantabio) consisting of 5μl of 2× Perfecta mix, 0.4μl of 2μM primer mix (fwd and rev primers), adjusted to 10μl of final volume with ddH2O water. Control samples were included to check amplification arising from contaminating genomic DNA and from the primer-dimer formation. All the samples were run in technical triplicates on Quantstudio 3 Real-time PCR system (Thermo Fisher Scientific). The transcript was amplified using the following conditions: 95°C for 3min followed by 50cycles of 95°C for 30s, 62°C for 30s, and 72°C for 30s. A melting curve profile was generated to determine the formation of a single amplification product. Primer efficiency for all the genes was determined using the formula Efficiency (%)=100×(10-1/slope-1) with 5-fold dilution series of template cDNA. Comparison of fold change between different target genes was using ∆∆CT method, and gyrase A (gyrA) gene was selected as reference for normalization (Livak and Schmittgen, 2001). Data analysis and statistical tests for RT-qPCR were performed on GraphPad Prism 9.
PCR and Agarose Gel Electrophoresis
The RT-qPCR primers (Supplementary Table ST-1; housekeeping genes gyrA, 16S rRNA, and recA, and target genes nanA and nanH) were tested on cDNA using end-point PCR. The PCR reaction conditions were 95°C for 2min followed by 35cycles of 95°C for 30s, 62°C for 30s and 72°C for 30s. The amplification temperature of the primer pairs was confirmed using gradient PCR from 60 to 64°C (data not shown). The products of PCR reactions were analyzed on a 2% agarose TBE gel. The images were captured using a Gel-doc XR (Biorad).
Nucleotide Sequence Accession Numbers
GenBank accession numbers of Shewanella whole genome sequences used for fastANI and UBCG analysis were as follows: NZ_CP050313.1, NZ_CP047422.1, NZ_CP018456.1, NZ_CP033575.1, NZ_CP046378.1, NC_008700.1, NC_017571.1, NC_017579.1, NC_009052.1, NC_009665.1, NC_009997.1, NC_011663.1, NC_016901.1, NZ_CP028730.1, NZ_CP028355.1, NZ_LS483452.1, NZ_CP022358.1, NZ_CP045857.1, NC_007954.1, NZ_CP041783.1, NC_008345.1, NC_010334.1, NZ_CP020472.1, NZ_CP020373.1, NZ_CP034015.1, NC_009092.1, NZ_CP041153.1, NZ_CP022272.1, NZ_CP036200.1, NC_004347.2, NC_009901.1, NC_011566.1, NZ_CP041036.1, NZ_CP014782.1, NC_017566.1, NC_009438.1, NZ_CP046329.1, NZ_LR134321.1, NZ_LR134303.1, NZ_CP028435.1, NC_009831.1, NC_008577.1, NZ_CP048031.1, NZ_CP022089.2, NZ_CP039928.1, NC_008321.1, NC_008322.1, NZ_CP015194.1, NZ_CP041329.1, NZ_CP041151.1, NC_008750.1, NZ_CP023019.1, NZ_CP041614.1, NC_014012.1, and NC_010506.1.
Shewanella glacialimarina TZS-4T Grows Well in Enriched Marine Broth
Bacterial strain S. glacialimarina TZS-4T was originally isolated from the Baltic Sea ice outside of Tvärminne Zoological station, Hanko, Finland, and cultivated on Zobell media containing undefined Baltic Sea water (Luhtanen et al., 2014). To standardize the growth media and remove any influence of seasonal fluctuation in seawater composition, we investigated if marine broth (MB) could be used to cultivate S. glacialimarina TZS-4T. First, the salt concentration of the Baltic Sea water used was determined to be 6.9g/L, which is consistent with brackish coastal water. We then compared the growth in Zobell media and in rMB, which has a defined composition with a total salt concentration of ~7.8g/L. rMB was found to be well suited for S. glacialimarina TZS-4T and supports a faster growth rate and a higher cell density at stationary phase (Supplementary Figure SF-2). Based on the favorable growth characteristics, rMB was adopted for all subsequent experiments in this study.
Phylogenetic Analysis and Average Nucleotide Identity Distinguishes S. glacialimarina TZS-4T as a Distinct Shewanella Species
Genomic DNA sequencing of S. glacialimarina TZS-4T was performed on a PacBio Sequel platform and assembled using hierarchical genome-assembly process (HGAP) analysis (Chin et al., 2013), achieving a coverage depth of 90.54×, indicating a high read confidence. The number of PacBio sequence reads was 57,211 with a N50 read length of 16,906bp. Its genome comprises a single circular chromosome of 4,456kb with a GC content of 41.12%, featuring 3,906 predicted open reading frames (ORFs) with 3,725 protein-coding sequences (Figure 1).3 There are a total of 118 RNA coding genes, 91 of which encode for tRNA and 22 for rRNA, of which eight encode for 5S rRNA, seven for 16S rRNA, and seven for 23S rRNA. In addition, two complete and one incomplete prophage segments were identified (Figure 1) by phispy tool on the RAST server and PHASTER (Arndt et al., 2016) with hits to predicted proteins of, e.g., phage 1/44, which was previously reported to infect the cold-active Shewanella sp. #44 strain (Luhtanen et al., 2014; Senčilo et al., 2015).
Figure 1. Genome map of Shewanella glacialimarina TZS-4T. Circular map from outside to center; coding DNA sequence (CDS) forward strand and CDS reverse strand are shown in magenta and green, respectively. tRNA genes are depicted as orange triangles and rRNA genes as blue circles. Prophage gene segments are represented as grey arrows. The innermost two segments denote GC content (salmon/blue) and GC skew (purple/olive) around the chromosome. The map was produced using the DNAplotter tool on Artemis (Carver et al., 2012).
Next, we wanted to position S. glacialimarina TZS-4T in regard to other Shewanella species. To this end, we constructed phylogenetic trees based on the 16S rRNA gene by maximum likelihood, neighbor-joining, and maximum parsimony methods with 100 replicates as bootstrapping and compared the tree topologies. The trees included 29 Shewanella strains representing 13 different phenotypic sub-groups with Vibrio and Pseudomonas strains as out-groups. The phylogenetic analyses yielded equivalent results (Supplementary Figure SF-3), with maximum likelihood showing that S. glacialimarina TZS-4T forms a differentiated branch with S. aestuarii SC18 (97.25%) and S. denitrificans OS217 (96.53%) branching next to it (Figure 2A). Similar results were obtained using AAI-profiler which compares the predicted proteome of an organism to Uniprot database (Supplementary Figure SF-4; Medlar et al., 2018). Limiting the phylogenetic analysis to the 92 most conserved genes yields a similar distribution, with cold-active Antarctic and Arctic Sea bacteria, such as Shewanella sp. Arc9-LZ, S. aestuarii SC18, S. polaris, and S. frigidimarina clustering alongside S. glacialimarina TZS-4T (Supplementary Figure SF-5).
Figure 2. Phylogenetic analysis of S. glacialimarina TZS-4T. (A) Phylogenetic tree based on Maximum Likelihood tree construction of 16S rRNA sequences with 100 replicates as bootstrap. The taxonomic position of S. glacialimarina TZS-4T (bolded) alongside other closely related strains of genus Shewanella. Pseudomonas putida is shown as an outgroup at the end of the tree. (B) Phylogenetic tree constructed by the orthroANI tool showing the position of S. glacialimarina TZS-4T and nine closely related Shewanella species. ANI score matrix is represented as a heat map from blue (low score) to red (high score).
To obtain independent confirmation, two bioinformatic tools – fastANI and orthroANI – were used to calculate the average nucleotide identity (ANI) values for S. glacialimarina TZS-4T in relation to other Shewanella species. For fastANI analysis, the comparison was made against whole genome sequences (WGS) of 55 Shewanella strains available in NCBI. The result showed that S. aestuarii yielded the highest ANI score with 611 matching sequences S. glacialimarina TZS-4T, followed by S. Arc9-Lz and S. frigidimarina (Supplementary Table ST-2). Further comparison was done using orthroANI on a limited subset of nine closely related Shewanella species, which delineated S. glacialimarina TZS-4T (ANI=75.31) along with S. aestuarii (ANI=76.08; Figure 2B). The ANI scores derived from fastANI and orthroANI are both below 95%, which is currently considered as the benchmark value for classifying bacterial species (Richter and Rosselló-Móra, 2009). Taken together, this implies that S. glacialimarina TZS-4T is genetically divergent and constitutes a distinct Shewanella species, which clusters with other cold-active marine Shewanella species.
Phenotypic Characterization of S. glacialimarina TZS-4T Uncovers a Resilient, Cold-Active Environmental Bacterium
To determine the characteristics of S. glacialimarina TZS-4T, we grew the bacteria on rMB-agar where they form smooth, circular, and convex colonies with a diameter of 1–2mm and a milky white appearance. Gram staining confirmed that S. glacialimarina TZS-4T is a Gram-negative, rod-shaped bacterium. Interestingly, after 3days of growth the colonies turn reddish-brown. A similar phenomenon was observed when the bacteria were grown for up to 29h in rMB-broth at 15 and 25°C (Supplementary Figure SF-6). We also found that pigment formation is not affected by light intensity (data not shown), but the cells are less pigmented at low growth temperatures (Supplementary Figure SF-6). This suggests that the pigment may have a role in temperature or UV-radiation adaptation, as was previously proposed for a S. frigidimarina strain isolated from Antarctic glacier snow (Martin-Cerezo et al., 2015), where the red pigment in question was identified as cytochrome c3. Indeed, a closer inspection of the S. glacialimarina TZS-4T genome revealed genes encoding for a putative cytochrome c3 family protein.
Next, we examined the physiology of S. glacialimarina TZS-4T using transmission electron microscopy (TEM). This showed a rod-shaped bacterium 0.5μm in diameter and ~2μm in length (Figure 3A) with clearly visible long polar flagella (Figure 3B). We also observed various bead-shaped extensions surrounding the bacteria. To visualize these, SEM analysis was performed, exposing a granular structure extending from the outer membrane of the bacteria (Figure 3C; highlighted in the insert). These extensions resemble a previously identified nanowire-like structure in S. oneidensis MR-1 (Gorby et al., 2006). The nanowires function as channels for extracellular electron transport in bacterial communities, allowing electron transfer from external surfaces, such as oxidized metals (Pirbadian et al., 2014). Genes encoding for nanowire-like structures, such as omcA/mtrC, were also identified in the S. glacialimarina TZS-4T genome, but the potential function of these structures is unclear and remains to be determined.
Figure 3. Electron microscopy visualization of S. glacialimarina TZS-4T. (A,B) Cells visualized with TEM after negative staining with 2% uranyl acetate. (C) Cells visualized with SEM. The insert (white box at the top right) shows a close-up of the nanowire-like structures originating from the outer membrane of the bacteria.
To further study the phenotypic characteristics of S. glacialimarina TZS-4T, we performed a number of standardized panel tests, comparing the results to the related psychrophilic bacterium S. frigidimarina (type strain ACAM 591) and the mesophilic reference bacterium S. baltica (type strain LMG 2250; Table 1). Shewanella glacialimarina TZS-4T is a motile, facultatively anaerobic bacteria that utilizes glucose and maltose as its sole source of carbon. However, the bacterium does not produce H2S, which is a characteristic feature of many Shewanella strains associated with fish spoilage in cold storage. Furthermore, S. glacialimarina TZS-4T is also capable of cytochrome oxidase and catalase production.
Table 1. Phenotypic characteristics of S. glacialimarina TZS-4T alongside other Shewanella-type strains.
Next, to determine the growth rate, we explored different physiological conditions, such as temperature, pH, salinity, and hydrogen peroxide induced oxidative stress. We observed that S. glacialimarina TZS-4T is psychrophilic and grows from low to ambient temperatures in a range between 0 and 25°C, with an optimum at 15°C (Figure 4A; Supplementary Figure SF-2). However, growth at 0°C requires preconditioning (growth at 4°C) of the starter cultures, as those grown at 15°C do not cope with the cold-shock and yield no growth at freezing conditions. During the first 6h post-inoculation, the growth rate is the slowest at 0 and 4°C, whereas it is the fastest at 25°C. On the other hand, there is an earlier cessation of growth at high temperature, as the cell density recorded at 24h post-inoculation is lower at 25°C than at 15°C (Figure 4A). Accordingly, at 24h we also observed a lower number of viable cells at 25°C (~1 ×109cfu/ml) than at 15°C (~8 ×109 cfu/ml), and the cells lose their ability to grow altogether at 28°C (data not shown). Indeed, a significant loss of viability is also observed during the 24–48h post-inoculation window (Supplementary Figure SF-2), which is consistent with the culture entering the death phase. S. glacialimarina TZS-4T has a propensity for neutral pH conditions with an optimum at pH 7.0, although it tolerates pH ranging from 5.0 to 9.5 (Figure 4B). It also has a high NaCl tolerance, with 35 and 55g/L NaCl reducing growth but not abolishing it (Figure 4C). Since the oxygen solubility from the environment is increased at low temperature (Karbowiak et al., 2009), we analyzed the ability of the cells to tolerate reactive oxygen species by supplementing the growth media with different concentrations of hydrogen peroxide (H2O2; Figure 4D). The cells grow in the presence of 1mM H2O2, but 2mM and 4mM H2O2 concentrations resulted in the complete cessation of growth.
Figure 4. Characterization of growth conditions of S. glacialimarina TZS-4T bacteria. Growth curves at different (A) temperature, (B) pH, (C) NaCl, and (D) hydrogen peroxide concentration. The error bars denote SD.
Finally, we investigated the antibiotic resistance of S. glacialimarina TZS-4T and found it to be susceptible to chloramphenicol, kanamycin, and tetracycline, but resistant to ampicillin (Table 1). Indeed, sequence analysis revealed that S. glacialimarina encodes resistance genes for ampicillin, including class-D β-lactamase and metallo-β-lactamase. Other virulence-associated factors, such as putative T1SS-secreted agglutinin RTX (cytotoxin), putative collagenase genes, and clp proteases, were also identified from the genomic sequence.
Unsaturated Fatty Acid Production Responds Dynamically to Changes in Temperature
To further investigate the cold adaptation of S. glacialimarina TZS-4T, we analyzed the FA modulation as it is a well-known strategy of cold-active bacteria to maintain membrane fluidity and viability. Our analysis revealed that palmitoleic acid (16:1n-7) was the most prevalent FA component of S. glacialimarina TZS-4T, constituting 32–37mol% of the total FAs, with the highest values detected at 5°C. In general, iso-13:0 and iso-15:0 accounted for more than 10mol% and were followed by 12:0, 14:0, and 16:0 that constituted close to 5% (Supplementary Table ST-3). A significant increase of certain unsaturated FAs, including 7-tetradecenoic acid (14:1n-7), eicosatetraenoic acid (20:4n-3), stearidonic acid (18:4n-3), and eicosapentaenoic acid (20:5n-3), was observed at low temperatures (Figure 5). Indeed, at 5°C the ratio of 20:5n-3 had raised close to 6mol%. A previous study on S. electrodiphila MAR441 showed similar temperature-dependent regulation of the relative concentration of certain n-3 FAs that have been associated with improved membrane fluidity and increased cell viability at low temperatures (Zhang and Burgess, 2017).
Figure 5. Relative concentrations of unsaturated FAs in S. glacialimarina TZS-4T analyzed by GC-FID and described as fold change. Figu re depicts the FA mol% values of S. glacialimarina TZS-4T grown to exponential phase (OD600=0.8) at 5 and 25°C normalized to the values at 15°C. All FAs are abbreviated as (carbon number):(number of double bonds) n-(position of the first double bond calculated from the methyl end). Error bar shows the SD, Welch t-test with Bonferroni correction, two-tailed *** p ≤0.001 and ** p=0.002 w.r.t 15°C.
Sialic Acid Metabolic Genes Are Expressed in S. glacialimarina TZS-4T
Interestingly, S. glacialimarina TZS-4T features a complete sialic acid metabolism gene cluster in its chromosome, including nanA that encodes for the major catalytic enzyme N-acetylneuraminate lyase and the sialidase encoding gene nanH. Upon blastn search, two distantly related Shewanella species, S. pealeana and Shewanella sp. YLB-09, also contain a similar sialic acid metabolism gene cluster (Figure 6A). Phylogenetic comparison using nanA revealed that among Shewanella, only these species carry this gene cluster, which is otherwise commonly found among pathogenic bacteria, such as Vibrio cholerae and Yersinia pestis (Figure 6B). Since other Shewanella species, such as S. putrefaciens, have been reported to cause serious health disorders among freshwater fish (Paździor et al., 2019), we hypothesize that the sialic acid metabolism gene cluster may play a role as a putative pathogenicity factor in S. glacialimarina TZS-4T. To determine if the nanA and nanH genes are transcribed at normal growth conditions, we performed an RT-qPCR expression analysis at 5, 15, and 25°C. Our results show that nanA and nanH are stably expressed at 5 and 15°C, although a significant upregulation occurs at 25°C w.r.t 15°C (Figure 6C). Furthermore, nanA and nanH expression remains stable throughout the exponential growth phase and only a slight, albeit not statistically significant downregulation is observed at the stationary phase for cells grown at 5°C (Supplementary Figure SF-7). This shows that the sialic acid enzymes genes are transcribed in S. glacialimarina TZS-4T and there is a significant temperature-dependent expression profile. However, the upregulation seen at 25°C suggests that the sialic acid metabolism gene cluster most likely does not further cold-active growth or a role in food spoilage during cold storage.
Figure 6. (A) Sialic acid metabolism gene cluster in S. glacialimarina TZS-4T, Shewanella pealeana, and Shewanella sp. YLB-09. Genes involved in the sialic acid metabolism shown in different colors located around a central N-acetylneuraminate lysase (nanA) colored in red. (B) Phylogenetic tree of N-acetyl neuraminate lyase (nanA) from S. glacialimarina TZS-4T and closely related nanA genes from other bacterial species. The sequences were aligned using ClustalW tool, and tree was constructed using 100 replicates as bootstrapping value on MEGA X software. (C) Expression analysis by RT-qPCR of sialic acid catabolic enzymes genes in S. glacialimarina TZS-4T at different temperatures. SD is shown by the error bars, independent t-test, two-tailed *** p ≤0.001 w.r.t 15°C.
Bacteria of the genus Shewanella, family Shewanellaceae, are capable of inhabiting a wide range of aquatic environments due to their versatile physiological characteristics (Lemaire et al., 2020). The Baltic Sea is one of the largest brackish water basins in the world and is partly covered by ice during the winter months (Granskog et al., 2006). During ice formation, channels containing brine are formed (Granskog et al., 2006), which serve as habitats for primary producers, such as unicellular algae, and heterotrophic bacteria, including Shewanella (Kaartokallio et al., 2008). Here, we have further characterized Shewanella strain #4, which was previously isolated from a Baltic Sea ice sample (Luhtanen et al., 2014), by sequencing its genome and studying its physiological properties, thereby uncovering it to be a new species that we name Shewanella glacialimarina TZS-4T. Shewanella glacialimarina TZS-4T is cold-active, able to grow in freezing temperatures, and in high salinity (up to 6% NaCl; Table 1; Figure 4). Such conditions are fairly common in the Baltic Sea, where Shewanella are most often found at depths corresponding to the oxic/anoxic interface (Brettar et al., 2002; Rak and Wieczorek, 2012), which is mixed annually (Hordoir et al., 2015). Furthermore, the shift toward freezing temperature occurs gradually in the Baltic Sea. This gradual change is crucial for S. glacialimarina TZS-4T, as it enables the bacteria to adapt to the prevailing environmental conditions. Indeed, we have shown that growth of S. glacialimarina TZS-4T does not recover following a sudden decrease in temperature from 15 to 0°C, as only bacteria that are preconditioned, i.e., maintained at 4°C, can survive and grow at 0°C (Table 1). During the cold-adaptation phase, bacteria generally produce cold-shock proteins (Csp) and store solutes, such as glycerol, sucrose, or mannitol, which provide both cryoprotection and yield an additional carbon source (Tribelli and López, 2018).
Shewanella glacialimarina TZS-4T also has additional strategies to cope with the conditions native to its habitat. This is exemplified by the dynamic changes in unsaturated FA production in response to growth temperature. For the cells grown at low temperatures, we observed large relative increases for 14:1n-7, 20:4n-3, 18:4n-3, and 20:5n-3, whereas the profile was shifted toward relatively long saturated and monounsaturated FAs at high temperatures (Figure 5; Supplementary Table ST-3). Such a temperature-specific regulation has previously been reported to modulate lipid composition of bacterial cell membrane maintaining membrane fluidity and cell viability (Zhang and Burgess, 2017). On the other hand, temperature-regulated sialic acid metabolism also directs toward another subtle change in bacterial metabolism upon experiencing environmental change (Figure 6C). However, our initial hypothesis was that sialic acid metabolism might be part of the strategy for S. glacialimarina TZS-4T to survive in cold temperatures. In contrast, expression analysis of the nanA and nanH genes uncovered a statistically significant upregulation at 25°C instead of 5°C. This suggests that sialic acid metabolism might instead be linked to nutrient utilization. Phytoplankton and fungi tend to be prevalent in the Baltic Sea during summer times when the temperature is above 10°C (Kaikkonen et al., 2020). These eukaryotes have peripheral sialic acid that could act as a stimulus for sialic acid metabolising bacteria to exploit this resource. A previous study found Shewanella to be the third most abundant bacteria associated with Pseudonitzschia fraudulenta in the 2010 algal bloom in Monterey Bay, California (Sison-Mangus et al., 2016). Therefore, it is possible that certain bacterial communities possessing distinct catalytic enzymes, such as those found in S. glacialimarina TZS-4T, could have a role in degrading phytoplankton-derived organic matter.
Despite the insight that the nanA gene cluster is expressed in S. glacialimarina TZS-4T, we still lack a detailed understanding of the role and function of these gene products in the life cycle of the bacterium. Furthermore, the presence of some pathogenicity-associated traits, such as alpha hemolysis on blood agar, the type IV secretion system, and antibiotic resistance, may hint toward potential pathogenic properties. Notably, S. putrefaciens has been reported to cause disease among freshwater fish (Paździor et al., 2019) and S. algae is a potential human pathogen (Torri et al., 2018). However, S. glacialimarina TZS-4T does not harm mammalians as it cannot grow above 25°C, nor does the presence of virulence-associated gene clusters constitute sufficient evidence for a potential disease-causing role in fish. Hence, further studies are needed to determine the function and utilization of these genetic elements.
Here, we have further characterized a novel cold-active Shewanella isolate named S. glacialimarina TZS-4T. The cold-active nature of S. glacialimarina TZS-4T makes it an interesting model to study transcriptional and translational adaptations that facilitate metabolism at low temperatures. For instance, elucidating the RNA-based regulatory mechanisms that govern translation – such as the dynamics of post-transcriptional RNA modification, which is known to among other facilitate the translation of stress-specific gene transcripts (Koh and Sarin, 2018) – may offer further insights into the modulation of cellular responses to environmental stress. Moreover, S. glacialimarina TZS-4T could also provide a platform to study cold-active enzymes, such as sialidase and N-acetylneuramidase. Consequently, based on our phylogenetic, phenotypic, and genomic characterization, we present S. glacialimarina TZS-4T as a new species belonging to the genus Shewanella.
Data Availability Statement
The data presented in the study are deposited in the GenBank repository, accession number CP041216.
MSQ and LPS conceived the study. MSQ, ML, M-MH, and BG-Z performed the experimental work. MSQ, RK, and ER analyzed the data with input from DB and LPS. MSQ wrote the first draft of the manuscript. ML, ER, and LPS edited the first draft of the manuscript. All authors contributed to the article and approved the submitted version.
This work was supported by the Academy of Finland Academy Research Fellow program (grants #294917, #307215, and #327181 to LPS) and Academy Professor program (grants #255342, #256518, and #283072 to DB); the Novo Nordisk Foundation Emerging Investigator grant (#NNF18OC0054454 to LPS); and the Sigrid Jusélius Foundation Young Group Leader grant (to LPS). Open access fees are covered by the Helsinki University Library.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
The authors wish to thank Salla Kalaniemi and Sanna P. Sihvo for technical assistance, as well as Pavlína Gregorová for her valuable advice on RT-qPCR analysis. The authors also thank the Electron Microscopy Unit of the Institute of Biotechnology, University of Helsinki, for providing laboratory facilities. All members of the RNAcious laboratory are thanked for their critical feedback and supportive discussions. MSQ is a fellow of the Doctoral Program in Microbiology and Biotechnology, University of Helsinki.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.737641/full#supplementary-material
Ackman, R. G. (2007). “Application of gas—liquid chromatography to lipid separation and analysis: qualitative and quantitative analysis,” in Fatty Acids in Foods and their Health Implications. 3rd Edn. ed. C. K. Chow (Boca Raton, Florida US: CRC press), 47–65.
Arndt, D., Grant, J. R., Marcu, A., Sajed, T., Pon, A., Liang, Y., et al. (2016). PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 44, W16–W21. doi: 10.1093/nar/gkw387
Besemer, J. (2001). GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 29, 2607–2618. doi: 10.1093/nar/29.12.2607
Bozal, N., Montes, M. J., Tudela, E., Jiménez, F., and Guinea, J. (2002). Shewanella frigidimarina and Shewanella livingstonensis sp. nov. isolated from Antarctic coastal areas. Int. J. Syst. Evol. Microbiol. 52, 195–205. doi: 10.1099/00207713-52-1-195
Brettar, I., Christen, R., and Höfle, M. G. (2002). Shewanella denitrificans sp. nov., a vigorously denitrifying bacterium isolated from the oxic-anoxic interface of the Gotland deep in the Central Baltic Sea. Int. J. Syst. Evol. Microbiol. 52, 2211–2217. doi: 10.1099/00207713-52-6-2211
Carver, T., Harris, S. R., Berriman, M., Parkhill, J., and Mcquillan, J. A. (2012). Artemis: an integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinformatics 28, 464–469. doi: 10.1093/bioinformatics/btr703
Cha, Q.-Q., Ren, X.-B., Sun, Y.-Y., He, X.-Y., Su, H.-N., Chen, X.-L., et al. (2020). Shewanella polaris sp. nov., a psychrotolerant bacterium isolated from Arctic brown algae. Int. J. Syst. Evol. Microbiol. 70, 2096–2102. doi: 10.1099/ijsem.0.004022
Chin, C.-S., Alexander, D. H., Marks, P., Klammer, A. A., Drake, J., Heiner, C., et al. (2013). Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat. Methods 10, 563–569. doi: 10.1038/nmeth.2474
Christie, W. (2019). LIPID MAPS Lipidomics Gateway, The Lipid Web, Methyl Esters of Fatty Acids, Archive of Mass Spectra., Available at: https://lipidmaps.org/resources/lipidweb/index.php?page=ms/methesters/me-arch/index.htm (Accessed July 8, 2021).
Fredrickson, J. K., Romine, M. F., Beliaev, A. S., Auchtung, J. M., Driscoll, M. E., Gardner, T. S., et al. (2008). Towards environmental systems biology of Shewanella. Nat. Rev. Microbiol. 6, 592–603. doi: 10.1038/nrmicro1947
Fredrickson, J. K., Zachara, J. M., Kennedy, D. W., Duff, M. C., Gorby, Y. A., Li, S.-M. W., et al. (2000). Reduction of U(VI) in goethite (α-FeOOH) suspensions by a dissimilatory metal-reducing bacterium. Geochim. Cosmochim. Acta 64, 3085–3098. doi: 10.1016/S0016-7037(00)00397-5
Gorby, Y. A., Yanina, S., Mclean, J. S., Rosso, K. M., Moyles, D., Dohnalkova, A., et al. (2006). Electrically conductive bacterial nanowires produced by Shewanella oneidensis strain MR-1 and other microorganisms. Proc. Natl. Acad. Sci. 103, 11358–11363. doi: 10.1073/pnas.0604517103
Gram, L., Trolle, G., and Huss, H. H. (1987). Detection of specific spoilage bacteria from fish stored at low (0°C) and high (20°C) temperatures. Int. J. Food Microbiol. 4, 65–72. doi: 10.1016/0168-1605(87)90060-2
Gregorova, P., Sipari, N. H., and Sarin, L. P. (2020). Broad-range RNA modification analysis of complex biological samples using rapid C18-UPLC-MS. RNA Biol. doi: 10.1080/15476286.2020.1853385 [Epub ahead of print]
Hordoir, R., Axell, L., Löptien, U., Dietze, H., and Kuznetsov, I. (2015). Influence of sea level rise on the dynamics of salt inflows in the B altic S Ea. J. Geophys. Res. Oceans 120, 6653–6668. doi: 10.1002/2014jc010642
Ivanova, E. P., Flavier, S., and Christen, R. (2004). Phylogenetic relationships among marine Alteromonas-like proteobacteria: emended description of the family Alteromonadaceae and proposal of Pseudoalteromonadaceae fam. Nov., Colwelliaceae fam. Nov., Shewanellaceae fam. Nov., Moritellaceae fam. Nov., Ferri. Int. J. Syst. Evol. Microbiol. 54, 1773–1788. doi: 10.1099/ijs.0.02997-0
Jain, C., Rodriguez-R, L. M., Phillippy, A. M., Konstantinidis, K. T., and Aluru, S. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 9:5114. doi: 10.1038/s41467-018-07641-9
Kaartokallio, H., Tuomainen, J., Kuosa, H., Kuparinen, J., Martikainen, P. J., and Servomaa, K. (2008). Succession of sea-ice bacterial communities in the Baltic Sea fast ice. Polar Biol. 31, 783–793. doi: 10.1007/s00300-008-0416-1
Kaikkonen, L., Enberg, S., Blomster, J., Luhtanen, A.-M., Autio, R., and Rintala, J.-M. (2020). Autumn to spring microbial community in the northern Baltic Sea: temporal variability in bacterial, viral and nanoflagellate abundance during the cold-water season. Polar Biol. 43, 1193–1206. doi: 10.1007/s00300-020-02700-8
Karbowiak, T., Gougeon, R. D., Alinc, J.-B., Brachais, L., Debeaufort, F., Voilley, A., et al. (2009). Wine oxidation and the role of cork. Crit. Rev. Food Sci. Nutr. 50, 20–52. doi: 10.1080/10408390802248585
Kim, S.-J., Park, S.-J., Oh, Y.-S., Lee, S.-A., Shin, K.-S., Roh, D.-H., et al. (2012). Shewanella arctica sp. nov., an iron-reducing bacterium isolated from arctic marine sediment. Int. J. Syst. Evol. Microbiol. 62, 1128–1133. doi: 10.1099/ijs.0.031401-0
Kim, J.-Y., Yoo, H.-S., Lee, D.-H., Park, S.-H., Kim, Y.-J., and Oh, D.-C. (2016). Shewanella algicola sp. nov., a marine bacterium isolated from brown algae. Int. J. Syst. Evol. Microbiol. 66, 2218–2224. doi: 10.1099/ijsem.0.001014
Koh, C. S., and Sarin, L. P. (2018). Transfer RNA modification and infection: implications for pathogenicity and host responses. Biochim. Biophys. Acta Gene Regul. Mech. 1861, 419–432. doi: 10.1016/j.bbagrm.2018.01.015
Lee, I., Ouk Kim, Y., Park, S.-C., and Chun, J. (2016). OrthoANI: an improved algorithm and software for calculating average nucleotide identity. Int. J. Syst. Evol. Microbiol. 66, 1100–1103. doi: 10.1099/ijsem.0.000760
Lemaire, O. N., Méjean, V., and Iobbi-Nivol, C. (2020). The Shewanella genus: ubiquitous organisms sustaining and preserving aquatic ecosystems. FEMS Microbiol. Rev. 44, 155–170. doi: 10.1093/femsre/fuz031
Li, Z., Lin, S., Liu, X., Tan, J., Pan, J., and Yang, H. (2014). A freshwater bacterial strain, Shewanella sp. Lzh-2, isolated from Lake Taihu and its two algicidal active substances, hexahydropyrrolo[1,2-a]pyrazine-1,4-dione and 2, 3-indolinedione. Appl. Microbiol. Biotechnol. 98, 4737–4748. doi: 10.1007/s00253-014-5602-1
Luhtanen, A.-M., Eronen-Rasimus, E., Kaartokallio, H., Rintala, J.-M., Autio, R., and Roine, E. (2014). Isolation and characterization of phage–host systems from the Baltic Sea ice. Extremophiles 18, 121–130. doi: 10.1007/s00792-013-0604-y
MacDonell, M. T., and Colwell, R. R. (1985). Phylogeny of the vibrionaceae, and recommendation for two new genera, Listonella and Shewanella. Syst. Appl. Microbiol. 6, 171–182. doi: 10.1016/S0723-2020(85)80051-5
Martin-Cerezo, M. L., Garcia-Lopez, E., and Cid, C. (2015). Isolation and identification of a red pigment from the antarctic bacterium Shewanella; frigidimarina. Protein Pept. Lett. 22, 1076–1082. doi: 10.2174/0929866522666150915122247
Medlar, A. J., Törönen, P., and Holm, L. (2018). AAI-profiler: fast proteome-wide exploratory analysis reveals taxonomic identity, misclassification and contamination. Nucleic Acids Res. 46, W479–W485. doi: 10.1093/nar/gky359
Na, S.-I., Kim, Y. O., Yoon, S.-H., Ha, S.-M., Baek, I., and Chun, J. (2018). UBCG: up-to-date bacterial core gene set and pipeline for phylogenomic tree reconstruction. J. Microbiol. 56, 280–285. doi: 10.1007/s12275-018-8014-6
Nealson, K. H., and Saffarini, D. (1994). Iron and manganese in anaerobic respiration: environmental significance, physiology, and regulation. Annu. Rev. Microbiol. 48, 311–343. doi: 10.1146/annurev.mi.48.100194.001523
Paździor, E., Pękala-Safińska, A., and Wasyl, D. (2019). Phenotypic diversity and potential virulence factors of the Shewanella putrefaciens group isolated from freshwater fish. J. Vet. Res. 63, 321–332. doi: 10.2478/jvetres-2019-0046
Pirbadian, S., Barchinger, S. E., Leung, K. M., Byun, H. S., Jangir, Y., Bouhenni, R. A., et al. (2014). Shewanella oneidensis MR-1 nanowires are outer membrane and periplasmic extensions of the extracellular electron transport components. Proc. Natl. Acad. Sci. 111, 12883–12888. doi: 10.1073/pnas.1410551111
Rossello-Mora, R. A., Caccavo, F., Osterlehner, K., Springer, N., Spring, S., Schüler, D., et al. (1995). Isolation and taxonomic characterization of a halotolerant, facultatively iron-reducing bacterium. Syst. Appl. Microbiol. 17, 569–573. doi: 10.1016/s0723-2020(11)80078-0
Senčilo, A., Luhtanen, A.-M., Saarijärvi, M., Bamford, D. H., and Roine, E. (2015). Cold-active bacteriophages from the Baltic Sea ice have diverse genomes and virus-host interactions. Environ. Microbiol. 17, 3628–3641. doi: 10.1111/1462-2920.12611
Sison-Mangus, M. P., Jiang, S., Kudela, R. M., and Mehic, S. (2016). Phytoplankton-associated bacterial community composition and succession during toxic diatom bloom and non-bloom events. Front. Microbiol. 7:1433. doi: 10.3389/fmicb.2016.01433
Torri, A., Bertini, S., Schiavone, P., Congestrì, F., Matteucci, M., Sparacino, M., et al. (2018). Shewanella algae infection in Italy: report of 3 years' evaluation along the coast of the northern Adriatic Sea. New Microbes New Infect. 23, 39–43. doi: 10.1016/j.nmni.2018.01.002
Venkateswaran, K., Moser, D. P., Dollhopf, M. E., Lies, D. P., Saffarini, D. A., Macgregor, B. J., et al. (1999). Polyphasic taxonomy of the genus Shewanella and description of Shewanella oneidensis sp. nov. Int. J. Syst. Evol. Microbiol. 49, 705–724. doi: 10.1099/00207713-49-2-705
Vogel, B. F., Venkateswaran, K., Satomi, M., and Gram, L. (2005). Identification of Shewanella baltica as the most important H2S-producing species during iced storage of danish marine fish. Appl. Environ. Microbiol. 71, 6689–6697. doi: 10.1128/aem.71.11.6689-6697.2005
Keywords: Shewanella, cold-active bacteria, sea ice, sialic acid metabolism, polyunsaturated (essential) fatty acids
Citation: Qasim MS, Lampi M, Heinonen M-MK, Garrido-Zabala B, Bamford DH, Käkelä R, Roine E and Sarin LP (2021) Cold-Active Shewanella glacialimarina TZS-4T nov. Features a Temperature-Dependent Fatty Acid Profile and Putative Sialic Acid Metabolism. Front. Microbiol. 12:737641. doi: 10.3389/fmicb.2021.737641
Edited by:Richard Allen White III, University of North Carolina at Charlotte, United States
Reviewed by:Geeta Chhetri, Dongguk University Seoul, South Korea
Jinwei Zhang, University of Exeter, United Kingdom
Copyright © 2021 Qasim, Lampi, Heinonen, Garrido-Zabala, Bamford, Käkelä, Roine and Sarin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Leif Peter Sarin, firstname.lastname@example.org
†Present address: Berta Garrido-Zabala, Science for Life Laboratory, Department of Immunology, Genetics and Pathology, Rudbeck Laboratory, Uppsala University, Uppsala, Sweden