Qing Liu
Hong-Can Liu1- 1China General Microbiological Culture Collection Center (CGMCC), Institute of Microbiology, Chinese Academy of Sciences, Beijing, China
- 2State Key Laboratory of Microbial Resources, Institute of Microbiology, Chinese Academy of Sciences, Beijing, China
Numerous mountain glaciers located on the Tibetan Plateau are inhabited by abundant microorganisms. The microorganisms on the glacier surface are exposed to the cold, barren, and high-ultraviolet radiation environments. Although the microbial community composition on glaciers has been revealed by high-throughput sequencing, little is known about the microevolution and adaptive strategy of certain bacterial populations. In this study, we used a polyphasic approach to determine the taxonomic status of 11 psychrophilic Flavobacterium strains isolated from glaciers on the Tibetan Plateau and performed a comparative genomic analysis. The phylogenetic tree based on the concatenated single-copy gene sequences showed the 11 strains clustered together, forming a distinct and novel clade in the genus Flavobacterium. The average nucleotide identity (ANI) values among these strains were higher than 96%. However, the values much lower than 90% between them and related species indicated that they represent a novel species and the name Flavobacterium bomense sp. nov. is proposed. The core and accessory genomes of strains in this new Flavobacterium species showed diverse distinct patterns of gene content and metabolism pathway. In order to infer the driving evolutionary forces of the core genomes, homologous recombination was found to contribute twice as much to nucleotide substitutions as mutations. A series of genes encoding proteins with known or predicted roles in cold adaptation were found in their genomes, for example, cold-shock protein, proteorhodopsin, osmoprotection, and membrane-related proteins. A comparative analysis of the group with optimal growth temperature (OGT) ≤ 20°C and the group with OGT > 20°C of the 32 Flavobacterium type strains and 11 new strains revealed multiple amino acid substitutions, including the decrease of the proline and glutamine content and the increase of the methionine and isoleucine content in the group with OGT ≤ 20°C, which may contribute to increased protein flexibility at low temperatures. Thus, this study discovered a novel Flavobacterium species in glaciers, which has high intraspecific diversity and multiple adaptation mechanisms that enable them to cope and thrive in extreme habitats.
Introduction
Glaciers harbor large number of microorganisms, including viruses, bacteria, archaea, and microeukaryotes, which have the ability to survive at low temperatures (Boetius et al., 2015). Low temperature is an important evolutionary force in the diversification of microorganisms. The collection of these cold-adapted microbes and analysis of their genetic characteristics is very important for us to understand how they can live in such cold environment and how the high species diversity is maintained. In recent years, numerous studies have been conducted to determine the microbial community structure on glaciers (Liu et al., 2009, 2015a,b, 2017; Franzetti et al., 2013). These studies revealed the most abundant bacterial genera on glaciers, including Flavobacterium, Arthrobacter, Hymenobacter, Deinococcus, Cryobacterium, Polaromonas, and Sphingomonas.
The genus Flavobacterium was proposed by Bergey et al. (1923) with Flavobacterium aquatile as the type species. It belongs to the phylum Bacteroidetes (formerly the Cytophaga-Flavobacterium-Bacteroides group), and its description was emended by Bernardet et al. (1996). Currently, it contains 150 species and is the largest genus within the family Flavobacteriaceae. Members of Flavobacterium produce carotenoids and flexirubin-type pigments, which make them yellowish or orange (Reichenbach et al., 1980; Bernardet and Bowman, 2015). Flavobacterium strains live in a variety of environments throughout temperate and polar regions, including terrenes, lakes, oceans, glaciers, plants, animals, and so on (Bernardet and Bowman, 2015). Some strains isolated from animals are pathogenic. Flavobacterium was found to be one of the most abundant genera in the Arctic and Antarctic sea ice (Boetius et al., 2015). In a survey of bacterial diversity on the surfaces of mountain glaciers in China, Flavobacterium was one of the three most abundant genera with more than 5% of average abundance (Liu et al., 2015b). Up to now, a total of seven psychrophilic/psychrotolerant Flavobacterium species with validly published names have been collected from glaciers, including F. collinsense (Zhang et al., 2016), F. glaciei (Zhang et al., 2006), F. noncentrifugens (Zhu et al., 2013), F. xinjiangense (Zhu et al., 2003), F. omnivorum (Zhu et al., 2003), F. sinopsychrotolerans (Xu et al., 2011), F. tiangeerense (Xin et al., 2009), F. xueshanense (Dong et al., 2012), and F. urumqiense (Dong et al., 2012).
With the development of modern “omics” technologies, some cold-adaptation strategies of microorganisms have been revealed. These specific physiological mechanisms are associated with cellular membrane fluidity (Casanueva et al., 2010), compatible solutes (Welsh, 2000), antifreeze protein (Chattopadhyay, 2006), ice-binding proteins (Raymond et al., 2007, 2008), anti-nucleating proteins (Kawahara, 2002), cold-shock proteins (Phadtare, 2004), cold acclimation protein (Phadtare, 2004), DEAD-box RNA helicase (Kuhn, 2012), cold-active enzymes (Feller and Gerday, 2003; Damico et al., 2006; Feller, 2010; Kasana and Gulati, 2011), energy generation and conservation (Amato and Christner, 2009), and genome plasticity (Casanueva et al., 2010).
Although the development of high-throughput sequencing techniques enables us to analyze bacterial diversity and community composition in cold environments, questions remain about the microevolution of certain psychrophilic/psychrotolerant groups. For instance, to survive for a long time at low temperature, which kinds of adaptations to low temperature have evolved in these lineages? How do these bacteria endure the harsh conditions and dwell in cold environments? During a survey of bacterial diversity on the surface of four glaciers located in Bome County, Tibetan Autonomous Region, P.R. China, we collected 11 psychrophilic Flavobacterium strains. Using a polyphasic approach, we determined the taxonomic status of these strains and identified a novel species, which we proposed to be named Flavobacterium bomense sp. nov. Furthermore, the intraspecies diversity, driving forces of microevolution, and cold-adaptation strategies of these glacier-inhabiting Flavobacterium strains were analyzed by comparative genomics.
Materials and Methods
Bacterial Strains and Culture Conditions
A total of 11 strains were isolated from melt water and ice samples from the surface of the Laigu, Zepu, Renlongba, and Gawalong glaciers in Bome County, Tibetan Autonomous Region, P.R. China (Table 1). All isolates were picked up as yellow, round colonies from PYG (peptone, yeast extract, and glucose) medium (Liu et al., 2015a) and 1/4 R2A (BD Difco, Becton, Dickinson and Company, Franklin Lakes, NJ, USA) agar plate. Six type strains of closely related species, namely F. urumqiense CGMCC 1.9230T, F. sinopsychrotolerans CGMCC 1.8704T, F. tiangeerense CGMCC 1.6847T, F. xueshanense CGMCC 1.9227T, F. omnivorum CGMCC 1.2747T, and F. frigidarium CGMCC 1.9172T, were used as reference strains for comparative analysis. All strains were routinely incubated in PYG medium at 15°C.
Table 1. Information of the Flavobacterium strains analyzed in this study and GenBank accession nos. of the 16S rRNA gene.
DNA Extraction, Amplification, and Sequencing
Genomic DNA was extracted using the Genomic DNA Rapid Isolation Kit for Bacterial Cell (BioDev-Tech, Co., Beijing, China) following the manufacturer’s instructions. The 16S rRNA gene was amplified and sequenced using the universal primers 27F and 1492R (Lane, 1991). Sequencing was performed using a PRISM 3730XL DNA analyzer (Applied Biosystems, Foster City, CA, USA) at SinoGenoMax Co. (Beijing, China). Whole genome sequencing of the 11 strains was performed using the Illumina HiSeq 4,000 platform (Illumina, San Diego, CA, USA) according to the manufacturer’s protocols. The assemblies of short reads were performed using the SPAdes 3.11 program with default parameters (Bankevich et al., 2012). The quality of the genomes was assessed by CheckM (Parks et al., 2015) and QUAST v5 (Mikheenko et al., 2018).
Phylogenetic Analysis Based on 16S rRNA Genes
For identification of the new strains, the 16S rRNA gene sequences were submitted to the EzBioCloud server (Yoon et al., 2017) to search for their closely phylogenetic neighbors. After multiple sequence alignment of the 16S rRNA gene sequences by ClustalW (Thompson et al., 1994), the neighbor-joining (NJ) and maximum-likelihood (ML) trees were built with 1,000 bootstrap replicates using the MEGA V. 5.2 software (Tamura et al., 2011). Kimura’s two-parameter model (K2P) was used to calculate the genetic distances (Kimura, 1980). The GTR + G + I was selected as the best nucleotide substitution model for ML tree construction.
Comparative Genomic Analysis
The draft genomes determined in this study were submitted for annotation using RASTkt (Brettin et al., 2015). To further confirm the taxonomic status of the new strains, a total of 32 whole genome sequences of related type strains were obtained from the NCBI genome database (Supplementary Table S1). For construction of species trees, single-copy orthologues were selected from genomic sequences using GET_HOMOLOGUES (Contreras-Moreira and Vinuesa, 2013) and GET_PHYLOMARKERS programs (Vinuesa et al., 2018). The alignments were generated by Clustal Omega software (Sievers and Higgins, 2014) implemented in GET_PHYLOMARKERS. Phylogenetic trees were generated by ML algorithms with a GTR + F + R5 model in the IQ-TREE software (Nguyen et al., 2015) based on the concatenated gene sequences with 1,000 bootstrap replicates. The average nucleotide identity (ANI) values were calculated by the GET_HOMOLOGUES program. The pan-genome analysis was performed using the BPGA tool with default parameters to determine the core and accessory genes (Chaudhari et al., 2016). Recombination analysis was performed with ClonalFrameML (Didelot and Wilson, 2015). Gene gain and loss rates were determined by the Count program using a birth-and-death model (Csuros, 2010) with the homologous table, which was inferred by using OMCL and bidirectional best hit (BDBH) methods implemented in GET_HOMOLOGUES package with the following parameters: 60% identity and 75% coverage.
Phenotypic Characterization of Novel Species
A polyphasic taxonomic analysis of the 11 strains was performed in this study. The morphology of the colonies was determined after culturing the 11 strains on PYG agar for 7 days. The cellular morphology of strain RB1N8T was examined by transmission electron microscopy using a JEM-1400 transmission electron microscope (JEOL Ltd., Tokyo, Japan). The growth at different temperatures (4, 15, 20, 22, 25, and 28°C), the tolerance to NaCl (0–4.0% (w/v) at 0.5% intervals), and pH values (ranging from pH 5.0 to 10.0 with 1 pH unit intervals) were tested in PYG broth. Hydrolyses of casein, starch, and Tween 80 were performed according to Smibert and Krieg (1994). The presence of flexirubin-type pigment was examined using 20% KOH (w/v). The utilization of sole carbon source, enzyme activities, and other phenotypic characteristics were tested using the API 20E, 20NE, ID 32 GN, and ZYM strips (bioMérieux, Marcy-l’Étoile, France) according to the manufacturer’s instructions.
For analysis of cellular fatty acid composition, the cells of the tested strains were harvested from colonies on the same sectors of the PYG plates after incubation at 15°C. The extraction of saponified and methylated fatty acids was performed according to the protocol of MIDI 6.0 system (Sasser, 1990). The samples were separated and identified on an Agilent 6,890 N gas chromatography system (Agilent Technologies, Santa Clara, CA, USA) using the TSBA6 database. Respiratory quinones and polar lipids were determined in cells of strain RB1N8T, which were harvested from PYG broth after incubation at 15°C for 5 days. The extracts of respiratory quinones and polar lipids were analyzed according to reported methods (Tindall et al., 2007).
Nucleotide Sequence Accession Numbers
The GenBank accession numbers for the 16S rRNA gene sequences of the 11 Flavobacterium strains are MK346152–MK346162 (Table 1). The Whole Genome Shotgun projects have been deposited at DDBJ/ENA/GenBank under the accession numbers: RYDG00000000, RYDH00000000, RYDI00000000, RYDJ00000000, YDK00000000, RYDL00000000, RYDM00000000, RYDN00000000, RYDO00000000, RYDP00000000, and RYDF00000000, respectively (Supplementary Table S2).
Results and Discussion
Phylogenetic Analysis Based on 16S rRNA Gene Sequences
The average evolutionary divergence (K2P) over all the 11 16S rRNA gene sequence pairs was 0.001%. Accordingly, the sequence of strain RB1N8T was selected as a representative to compare with the EzBioCloud database. The highest 16S rRNA gene sequence similarities were found between strain RB1N8T and F. xueshanense Sr22T (97.77%), F. psychrolimnae LMG 22018T (97.49%), F. fryxellicola DSM 16209T (97.49%), and F. tiangeerense 0563T (97.40%), is lower than the 98.65% threshold value for species delineation (Kim et al., 2014), supporting the notion that the new group represents a novel species of the genus Flavobacterium. Phylogenetic analysis revealed that the 11 strains formed an independent lineage with strong bootstrap support of 100% (Supplementary Figure S1). Although some strains shared 100% sequence identities, strain diversity was noticeable in this new clade. The tree topologies constructed with ML algorithms were similar to that of the NJ trees (Supplementary Figure S2).
Taxonomic Features and Phylogenetic Analysis Based on Genomes
Except for the genome sequence of strain GSN2, in which reads contamination was detected, the genome sequence assemblies of the other 10 strains are high quality (Supplementary Table S2). Thus, the genome sequence of strain GSN2 was ignored in the subsequent analysis, except for the single-copy gene phylogenetic reconstruction. These genomes were similar in size, ranging from 3.20 to 3.82 Mb. The DNA GC content was calculated to be from 34.91 to 35.04%, which is in agreement with the GC content of their closely related strains F. xueshanense Sr22T (34.1%), F. psychrolimnae LMG 22018T (34.2%), F. fryxellicola DSM 16209T (34.6%), and F. tiangeerense 0563T (33.6%). The ANI values among the 11 strains ranged from 96.89 to 99.17%, which are higher than the proposed cut-off value for species boundary (95–96%; Richter and Rosselló-Móra, 2009), indicating that they belong to one species. Furthermore, the ANI values between these strains and their closely related species with validly published names are much lower than 90% (Supplementary Table S1), revealing that this new Flavobacterium group represents a novel species.
A total of 230 single-copy core genes were extracted from the genome sequences of 32 Flavobacterium type strains and 11 new strains. After sequence alignment, the concatenated single-copy gene sequences were used to construct ML and NJ phylogenetic trees with 1,000 bootstrap replicates (Figure 1), and the trees reconstructed by different methods showed identical topologies. The species tree clearly showed that the 11 tested strains clustered together, forming an independent and novel clade in the genus Flavobacterium. This novel clade and 10 other type strains, almost all of which were isolated from glaciers or Antarctic lakes, formed a larger lineage, revealing the diversification of Flavobacterium inhabiting in cold environments.
Figure 1. Phylogeny of the new strains and related species constructed from an ML analysis of the 230 concatenated single-copy gene sequences. Bootstrap values ≥70% based on 1,000 replicates are indicated at branch points. Bar = 0.05 nt substitutions per site. PR, proteorhodopsin; OGT, optimal growth temperature.
Pan-Genome of the Novel Species
To gain a deeper understanding of the intraspecies genomic diversity of glacier-inhabiting Flavobacterium species, pan-genome analysis was performed. In this study, the pan-genome refers to the total number of orthologous gene families in the 10 strains isolated here with uncontaminated genome assemblies. The core pan-plot showed that the size of the pan-genome increased unlimitedly with the addition of new genomes (Supplementary Figure S3). The more genomes added, the more orthologous clusters produced. Thus, the pan-genome of the novel species was open. The core genome decreased as the genomes were added one by one and ultimately became relatively constant. A total of 2,269 genes formed the core gene pool of the novel species. Additionally, the total number of accessory genes comprised more than half of the pan-genome. Every strain contained a certain number of strain-specific genes (unique), and the number varied considerably (61–403) depending on individual strains, indicating that the ongoing genetic flow led to the generation of strain specificity (Supplementary Table S3).
Both the core genes and accessory genes were involved in all the categories of COG functions (Figure 2). However, most of the core genes was classified into the basic functions, such as “amino acid transport and metabolism [E],” “translation, ribosomal structure, and biogenesis [J],” “energy production and conversion [C],” “coenzyme transport and metabolism [H],” “lipid transport and metabolism [I],” and “nucleotide transport and metabolism [F].” Most of the accessory genes was related to functions of “cell wall/membrane/envelope biogenesis [M],” “transcription [K],” “carbohydrate transport and metabolism [G],” “replication, recombination, and repair [L],” and “signal transduction mechanisms [T],” which may contribute to adaptation to changing environments. In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis (Supplementary Figure S4) revealed that the unique genes were mainly involved in carbohydrate, amino acid, lipid, energy metabolism pathways, and signal transduction, especially carbohydrate metabolism, which corresponded with the variable ability of carbon source utilization tested by API ID 32 GN (Supplementary Table S4). These genomic characteristics indicated the diversity of metabolic pathways in different Flavobacterium strains. Thus, the pan-genome analysis showed the distinct patterns of gene content and metabolism pathway within this new Flavobacterium species.
Figure 2. COG distribution for the pan-genome of 10 glacier-inhabiting Flavobacterium strains determined using the BPGA tool with default parameters. The x axis represents the relative percentage of the number of genes in each functional category.
Microevolution of the Novel Species
In their evolutionary history, bacteria experience frequent gene gain and loss (Polz et al., 2013; Vos et al., 2015) or undergo homologous recombination and mutation within gene families (Didelot and Maiden, 2010), all of which may result in the variation of the pan-genome. These varied parts of the pan-genome are essential to the ability of the bacteria to survive in their individual habitat (McInerney et al., 2017). In order to investigate the evolutionary history of Flavobacterium strains, gene gain and loss rates were estimated (Figure 3). The results showed that the rates of gene gain and loss were about the same in every strain, except for strain RB1N8T, which gained more genes than it lost. The gene gain and loss of these Flavobacterium strains may have resulted from horizontal gene transfer (HGT). Bacteria could benefit from HGT, which would introduce new functions to adapt to the changing environment (Zhaxybayeva and Doolittle, 2011).
Figure 3. Hierarchical clustering of 10 strains based on the heatmap of the orthologous genes. The presence and absence of the orthologous genes for each strain are indicated in red and blue color, respectively.
The contribution of homologous recombination and mutation to microevolution of the tested strains was further investigated by ClonalFrameML using a single-copy gene subset. The ratio of recombination to mutation rate (ρ/θ) and the ratio of nucleotide substitutions due to homologous recombination or mutation (r/m) were calculated based on 33 type strains and 10 new strains. Within species, ρ/θ and r/m values were 0.5064 and 2.1862, respectively. The values of ρ/θ and r/m between Flavobacterium species were 0.0052 and 0.3358. Accordingly, the rate of intraspecies recombination was much higher than that of interspecies. Although the r/m value was much lower than the values within many bacterial species (Vos and Didelot, 2009), homologous recombination still introduced more than twice the number of substitutions than that of mutations within the glacier-inhabiting Flavobacterium species.
Therefore, gain and loss of genetic material, recombination, and mutation all would result in the variations of the genomes. Such evolutionary events led to the formation of intraspecific diversity and phylogenetic cohesion of the novel Flavobacterium species.
Insights of the Genome Sequences Linked to the Strategy of Living in Cold Environment
A number of known genes related to cold adaptation were identified in the novel species and compared with the genomes of two phylogenetic-related mesophiles, namely F. granuli DSM 19729T and F. crassostreae LPB0076T (Supplementary Table S5). The genes of proteorhodopsin, polysaccharide transporter, and ice-binding protein (IBPs) were absent in the two reference genomes, but the other adaptation-related genes were also found in the two mesophiles.
Proteorhodopsin (PR) is a light-dependent proton pump of bacteria, which can utilize solar irradiance to produce ATP (Béjà et al., 2001) and consequently can enhance the growth rate of some bacteria, but not all of them (Giovannoni et al., 2005; Gómez-Consarnau et al., 2007, 2016; Lami et al., 2009). Fuhrman et al. (2008) suggested that PR may contribute to the physiological adaptation of bacteria to stressful conditions. Comparative genome analysis revealed the presence of a gene encoding PR in all the new strains as well as in 12 closely related type strains, which were isolated from glaciers or the Antarctic. However, the PR gene was absent in the genomes of the 20 type strains of Flavobacterium analyzed in this study, most of which was mesophilic (Figure 1). Accordingly, the spread of the PR gene in cold environments informed our understanding of its ecological function. We inferred that harboring PR genes might be particularly advantageous to the survival of Flavobacterium bacteria in frigid and barren environments.
Extracellular polysaccharides (EPSs) serve as cryoprotectants in bacteria living in marine or other cold environments (Nichols et al., 2005; Reid et al., 2006). All the tested strains in this study contained one gene copy of lptA, lptB, lptC, lptF, and lptG, which are related to the transportation of polysaccharide out of the membrane, suggesting that these Flavobacterium strains could produce and export EPSs to protect the cells against the harm of low temperature.
Ice-binding proteins (IBPs) inhibit the growth of ice crystals inside and outside the cells. The production of IBPs could make the bacteria counteract deleterious effects under subzero conditions (Mangiagalli et al., 2017). One gene copy of IBP was observed in the 10 Flavobacterium strains. IBP genes have also been identified in other psychrophilic bacteria, such as Antarctic Colwellia sp. SLW05 (Raymond et al., 2007) and Psychroflexus torquis (Mangiagalli et al., 2017)
Sigma factors perform all initiation of transcription in bacteria. Except for the single essential housekeeping σ that promotes the transcription of thousands of genes, many σs promote the transcription of specialized genes in response to a particular stress or stimulus (Feklistov et al., 2014). The σ70 factors are commonly related to gene transcription, stress response, cell development, and auxiliary metabolism. Multiple copies of the σ70 gene (rpoD) have been found to increase resistance to cold stress in psychrophilic bacteria (Riley et al., 2008; Mykytczuk et al., 2013; Dsouza et al., 2015). In this study, the gene rpoD was also found in the novel Flavobacterium strains. Additionally, in their genomes, there were 9–10 copies of the rpoE gene, whose product (σ24) was associated with regulating cellular responses to heat-shock and other stresses on cellular membrane and periplasmic proteins (Dsouza et al., 2015).
In cold environments, the solubility of O2 increases at low temperatures and more reactive oxygen species (ROS) are formed (Chattopadhyay, 2002). Thus, for survival under oxidative stress conditions, psychrophilic bacteria could remove the ROS through proteins encoded by some functional genes, such as sodA (superoxide dismutase), katE and katG (catalase), peroxiredoxin bcp types (thiol peroxidases), osmC/ohr (organic hydroperoxide reductase), trxB (thioredoxin reductase), and trxA (thiol-disulfide isomerase and thioredoxins). Several copies of these genes were found in the genomes of the novel glacier-inhabiting Flavobacterium strains, and the proteins encoded by these genes may contribute to their survival on surfaces of glaciers with low temperature and high UV radiation.
When bacteria are suddenly exposed to a cold environment, a number of physiological changes occur in the cells and a set of small molecule proteins is expressed. Genes encoding these cold-shock-inducible proteins, including cold-shock protein (cspA), ribosome-binding factor A (rbfA), translation initiation factors, IF-1 and IF-2 (infA, infB), polynucleotide phosphorylase (pnp), and transcription termination protein A (nusA), were identified in all 10 Flavobacterium strains. In these strains, just one gene copy of cold-shock protein (cspA) was identified. However, several copies were observed in Antarctic Arthrobacter (Dsouza et al., 2015), Colwellia psychrerythraea (Methé et al., 2005), Psychrobacter arcticus (Ayaladelrío et al., 2010), and Shewanella oneidensis (Gao et al., 2006). Carotenoid could stabilize the cellular membrane at low temperature and therefore contribute to the adaptation of carotenoid-pigmented bacteria to a cold environment (Fong et al., 2001; Dieser et al., 2010). The intact pathway of carotenoid biosynthesis was found in the 10 genomes, including the genes idi (isopentenyl-diphosphate delta-isomerase), crtB (phytoene synthase), crtI (phytoene dehydrogenase), crtY (lycopene beta cyclase), and crtZ (beta-carotene hydroxylase). The presence of these genes was consistent with the yellow color of the colonies of these strains.
Glycogen is thought to help bacteria resist to stressful condition of low temperature (Bresolin et al., 2006; Dalmasso et al., 2012). Three enzymes, glycogen synthase (glgA), glycogen branching enzyme (glgB), and glucose-1-phosphate adenylytransferase (glgC), participate in the biosynthesis of glycogen (Cifuente et al., 2016). In this study, one to two copies of these genes were found in the genomes of 10 tested strains, which would enable them to accumulate carbon and energy reserves to cope with the cold and barren environments.
Proline is an important organic metabolite and has been proposed to be a protective osmolyte (Hoffmann and Bremer, 2011). Genes involved in proline biosynthesis and Na+/proline symporter were found in the 10 tested genomes, but genes of other compatible solutes such as glycine betaine were absent.
Maintaining the permeability and fluidity of the cellular membrane is important for bacteria living in cold environments, and this could be achieved through the synthesis of polyunsaturated fatty acids (Methé et al., 2005). The genome analysis in this study revealed the presence of genes for fatty acid desaturases (des), which would introduce the unsaturated double bond into the saturated fatty acid. Concerning the cellular fatty acid compositions (%) of Flavobacterium strains, several polyunsaturated fatty acids, including iso-C15:1 G, C15:1 ω6c, iso-C16:1 H, C17:1 ω6c, C18:1 ω5c, summed in Feature 3 (C16:1 ω6c and/or C16:1 ω7c), summed in Feature 4 (iso-C17:1 I/ anteiso-C17:1 B), and summed in Feature 9 (iso-C17:1 ω9c/10-methyl C16:0), were detected (Supplementary Table S6).
In order to increase the tRNA flexibility, psychrophilic bacteria can posttranscriptionally incorporate dihydrouridine in the tRNA (Dalluge et al., 1997). Similar to some other cold-adapted bacterial groups (Saunders et al., 2003; Wang et al., 2008; Qin et al., 2014), two copies of the gene encoding tRNA-dihydrouridine synthase (dusB) were found in the Flavobacterium strains, which could increase the conformational flexibility of their tRNAs to guarantee their survival in cold glaciers.
The system of Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and its associated proteins (Cas) serve as an adaptive immune system for protecting prokaryotes against viral predators and foreign invaders (Hille and Charpentier, 2016). CRISPR-associated proteins Cas1 and Cas2, CRISPR-associated endonuclease Cas9, and CRISPR repeats (35–36 bp) were found present in the genomes of these Flavobacterium strains, indicating that the CRISPR-Cas system of Flavobacterium strains belonged to Type II system. There were 9–45 CRISPR repeats depending on the particular strain, showing a variety of CRISPR in the Flavobacterium strains. Some researchers have isolated bacteriophages from polar regions, deep sea, permafrost regions, high latitude lakes, and glaciers (Sawstrom et al., 2007; Li et al., 2016). Thus, the CRISPR-Cas systems found in these genomes indicate that viruses are also an important part of the glacier habitat.
At low temperatures, the most important property for proteins is the maintenance of sufficient flexibility, which can increase their interactions with substrates and reduce the required activation energy. The analysis of amino acid components revealed that the decrease in the number of proline residues (Fields, 2001), the introduction of methionine residues in the place of other buried hydrophobic residues (Marshall, 1997), and the decrease in the number of glutamine residues (Damico et al., 2006) were the strategies to increase the flexibility in cold-adapted proteins. In order to analyze the amino acid substitutions in proteins of strains with different optimal growth temperatures (OGT), the amino acid frequencies of the single-copy gene alignments of the 32 Flavobacterium type strains and 11 new strains were calculated. The OGT information of the type strains was acquired from the online catalogues of DSMZ1, JCM2, and CGMCC3. The strains with OGT ≤ 20°C contain more methionine, while the strains with OGT >20°C contain more proline and glutamine. Additionally, more isoleucine was found in the group with OGT ≤ 20°C. Our results suggested that the Flavobacterium strains may employed multiple amino acid substitutions to decrease protein stability and increase protein flexibility, which make the proteins more active at low temperatures (Figure 4).
Figure 4. Difference of amino acid residue frequencies between the two groups of strains with OGT ≤ 20°C and OGT > 20°C in the genus Flavobacterium. The boxplot was produced in “ggplot2” package implemented in R. The significant differences were identified using Student’s t test.
Phenotypic and Chemotaxonomic Characteristics
The new species of Flavobacterium was analyzed using phenotypic and chemotaxonomic methods. These strains were found to be Gram-negative, catalase and oxidase-positive, non-motile rods with yellow, and round colonies. No spores or flagella were detected (Supplementary Figure S5). These strains could grow in a temperature range of 0–22°C, and the maximum growth temperature for most of them was below 20°C (Supplementary Table S7). They were therefore considered to be a group of psychrophiles. They showed growth in a NaCl concentration range of 0–2.0/2.5 and a pH value range of 6.0/7.0–8.0/9.0, varying with different strains. Flexirubin-type pigment was absent. These strains could be easily distinguished from the type strains of their closely related species by their physiological and biochemical characteristics (Supplementary Table S7).
The profiles of cellular fatty acids in novel strains and reference strains are shown in Supplementary Table S6. The cellular fatty acid composition of novel strains and related type strains were similar, although there were some minor differences in certain components. The main fatty acids of the novel strains were summed Feature 3 (C16:1 ω6c and/or C16:1 ω7c, 17.4–33.8%) and iso-C15:0 (13.9–17.4%). The major hydroxyl fatty acids were iso-C15:0 3-OH (5.9–8.3%), iso-C16:0 3-OH (3.7–6.2%), and iso-C17:0 3-OH (4.4–7.6%). Menaquinone 6 (MK-6) was the only isoprenoid quinone detected in strain RB1N8T, which is consistent with the other members of the genus Flavobacterium (Bernardet and Bowman, 2015). Polar lipids of RB1N8T were phosphatidylethanolamine (PE), three unidentified polar lipids, and four unidentified aminolipids (Supplementary Figure S6), which is similar to the polar lipid profiles of other species of Flavobacterium (Bernardet and Bowman, 2015).
Conclusion
In this study, 11 bacterial strains isolated from glaciers in China were classified as a novel Flavobacterium species using a polyphasic taxonomy method, for which the name Flavobacterium bomensis sp. nov. is proposed. A comparative genomic analysis of these strains was performed. Their pan-genome was open, and every genome contained some unique genes, which showed the diversity of gene content and metabolic pathways in this species. The driving forces of microevolution were also investigated. The results revealed that the accessory genomes of the new species gained and lost genes at certain rates, while homologous recombination showed a twofold higher contribution to nucleotide substitutions than mutation in their core genomes. In addition, all genomes harbored numerous genes related to adaptation to cold and harsh environments. These genes were involved in various adaptive aspects, such as carotenoid biosynthesis, carbon and energy reserves, the fluidity of cellular membrane, osmotic and oxidative stress, proteorhodopsin, cold-shock protein, and ice-binding protein. Additionally, in the Flavobacterium group with OGT ≤ 20°C, multiple amino acid residue substitutions were found, which were involved in decreasing protein stability and increasing protein flexibility. In brief, the comparison and analysis help us gain a better understanding of the microevolution and the adaptive strategies of the dominant group Flavobacterium in the cryosphere.
Description of Flavobacterium bomensis sp. nov.
Flavobacterium bomensis (bo.men.′sis N.L. adj. Bomensis referring to Bome County, Tibetan Autonomous Region, P.R. China, from which the strains were isolated).
Cells are aerobic, Gram-negative, non-spore-forming, non-motile and non-gliding, chemoorganotrophic rods, 0.32–0.38 mm wide, and 1.04–2.08 mm long. Colonies grown on PYG agar for 9 days are circular, yellow, convex with entire margins, and about 1.5 mm in diameter. Growth occurs at 0–22°C (optimum, 10–15°C). The pH value range for growth is from 6.0 to 9.0 (optimum, pH 7.0). The salinity range for growth is from 0 to 2.5% NaCl. Catalase and oxidase are positive. No flexirubin-type pigments produced on PYG agar. Do not reduce nitrates to nitrites. The hydrolysis of gelatin and casein is variable.
No acid is produced from D-glucose, D-mannitol, inositol, L-rhamnose, D-sucrose, D-melibiose, and L-arabinose; acid production from D-sorbitol and amygdalin is variable. All strains can utilize D-glucose, D-sucrose, D-maltose, glycogen, and L-proline as sole carbon source and cannot utilize the following carbohydrate as sole carbon source: D-mannitol, L-fucose, D-sorbitol, caprate, valerate, citrate, L-histidine, 2-ketogluconate, 3-hydroxy-butyrate, 4-hydroxy-benzoate, L-rhamnose, D-ribose, inositol, itaconate, suberate, malonate, acetate, lactate, L-alanine, 5-ketogluconate, 3-hydroxy-benzoate, and L-serine. The utilization of salicin, D-melibiose, L-arabinose, N-acetyl-D-glucosamine, and propionate is variable.
Positive for starch hydrolysis, alkaline phosphatase, leucine arylamidase, valine arylamidase, and naphthol-AS-BI-phosphohydrolase. Enzyme activities of acid phosphatase and N-acetyl-β-glucosaminidase are variable. Negative for Voges-Proskauer test, indole and H2S production, esculin hydrolysis, citrate utilization, glucose fermentation, arginine dihydrolase, lysine decarboxylase, ornithine decarboxylase, urease, tryptophan deaminase, esterase (C4), esterase lipase (C8), lipase (C14), cystine arylamidase trypsin, α-chymotrypsin, α-galactosidase, β-galactosidase, β-glucuronidase, α-glucosidase, β-glucosidase, α-mannosidase, and α-fucosidase.
The most predominant cellular fatty acids are iso-C15:0, iso-C15:0 3OH, anteiso-C15:0, and summed in Feature 3 (comprising C16:1 ω7c and/or C16:1 ω6c). The major polar lipid is phosphatidylethanolamine. The DNA GC content is 34.9–35.1 mol%.
Strains were isolated from melt water and ice samples on the surface of Laigu, Zepu, Renlongba, and Gawalong glaciers in Bome County, Tibetan Autonomous Region, P.R. China, in 2016. The type strain is RB1N8T (= CGMCC 1.23902 = NBRC 113662).
Author Contributions
QL and Y-HX designed the project and analyzed the data. QL, Y-HX, and Y-GZ collected the samples. QL and Y-HX purified the strains. QL performed the bioinformatic analysis of the genome sequences. H-CL performed the fatty acid analysis. QL and Y-HX wrote the manuscript.
Funding
This study was supported by grant nos. 31670003 and 31600007 from the National Natural Science Foundation of China (NSFC) and grant no. 2015FY110100 from the National Science and Technology Foundation project.
Conflict of Interest Statement
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.
Acknowledgments
We thank Jing-Nan Liang of the Institute of Microbiology, Chinese Academy of Sciences (IMCAS), for her help in the use of transmission electron microscopes. We thank Bing-Da Sun of IMCAS for his kind help in the sample collection. We are grateful to Dr. Yofre C. of University of Portsmouth for English editing of the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01069/full#supplementary-material
Footnotes
References
Amato, P., and Christner, B. C. (2009). Energy metabolism response to low temperature and frozen conditions in Psychrobacter cryohalolentis. Appl. Environ. Microbiol. 75, 711–718. doi: 10.1128/AEM.02193-08
Ayaladelrío, H. L., Chain, P. S., Grzymski, J. J., Ponder, M. A., Ivanova, N., Bergholz, P. W., et al. (2010). The genome sequence of Psychrobacter arcticus 273-4, a psychroactive Siberian permafrost bacterium, reveals mechanisms for adaptation to low-temperature growth. Appl. Environ. Microbiol. 76, 2304–2312. doi: 10.1128/AEM.02101-09
Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477. doi: 10.1089/cmb.2012.0021
Béjà, O., Spudich, E. N., Spudich, J. L., Leclerc, M., and DeLong, E. F. (2001). Proteorhodopsin phototrophy in the ocean. Nature 411, 786–789. doi: 10.1038/35081051
Bergey, D. H., Harrison, F. C., Breed, R. S., Hammer, B. W., and Huntoon, F. M. (1923). Bergey’s manual of determinative bacteriology. (Baltimore: The Williams & Wilkins Co).
Bernardet, J. F., and Bowman, J. P. (2015). “Flavobacterium” in Bergey's manual of systematics of archaea and bacteria. ed. W. Whitman (John Wiley & Sons, Inc.), 1–75. doi: 10.1002/9781118960608.gbm00312
Bernardet, J. F., Segers, P., Vancanneyt, M., Berthe, F., Kersters, K., and Vandamme, P. (1996). Cutting a Gordian knot: emended classification and description of the genus Flavobacterium, emended description of the family Flavobacteriaceae, and proposal of Flavobacterium hydatis nom. nov. (basonym, Cytophaga aquatilis Strohl and Tait 1978). Int. J. Syst. Bacteriol. 46, 128–148. doi: 10.1099/00207713-46-1-128
Boetius, A., Anesio, A. M., Deming, J. W., Mikucki, J. A., and Rapp, J. Z. (2015). Microbial ecology of the cryosphere: sea ice and glacial habitats. Nat. Rev. Microbiol. 13, 677–690. doi: 10.1038/nrmicro3522
Bresolin, G., Neuhaus, K., Scherer, S., and Fuchs, T. M. (2006). Transcriptional analysis of long-term adaptation of Yersinia enterocolitica to low-temperature growth. J. Bacteriol. 188, 2945–2958. doi: 10.1128/JB.188.8.2945-2958.2006
Brettin, T., Davis, J. J., Disz, T., Edwards, R. A., Gerdes, S., Olsen, G. J., et al. (2015). RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci. Rep. 5:8365. doi: 10.1038/srep08365
Casanueva, A., Tuffin, M., Cary, C., and Cowan, D. A. (2010). Molecular adaptations to psychrophily: the impact of ‘omic’ technologies. Trends Microbiol. 18, 374–381. doi: 10.1016/j.tim.2010.05.002
Chattopadhyay, M. K. (2006). Mechanism of bacterial adaptation to low temperature. J. Biosci. Bioeng. 31, 157–165. doi: 10.1007/BF02705244
Chaudhari, N. M., Gupta, V. K., and Dutta, C. (2016). BPGA-an ultra-fast pan-genome analysis pipeline. Sci. Rep. 6:24373. doi: 10.1038/srep24373
Cifuente, J. O., Comino, N., Madariagamarcos, J., Lopezfernandez, S., and Garciaalija, M. (2016). Structural basis of glycogen biosynthesis in bacteria. Structure 24, 1613–1622. doi: 10.1016/j.str.2016.06.023
Contreras-Moreira, B., and Vinuesa, P. (2013). GET_HOMOLOGUES, a versatile software package for scalable and robust microbial pangenome analysis. Appl. Environ. Microbiol. 79, 7696–7701. doi: 10.1128/AEM.02411-13
Csuros, M. (2010). Count: evolutionary analysis of phylogenetic profiles with parsimony and likelihood. Bioinformatics 26, 1910–1912. doi: 10.1093/bioinformatics/btq315
Dalluge, J. J., Hamamoto, T., Horikoshi, K., Morita, R. Y., Stetter, K. O., and Mccloskey, J. A. (1997). Posttranscriptional modification of tRNA in psychrophilic bacteria. J. Bacteriol. 179, 1918–1923. doi: 10.1128/jb.179.6.1918-1923.1997
Dalmasso, M., Aubert, J., Briardbion, V., Chuat, V., Deutsch, S. M., Even, S., et al. (2012). A temporal-omic study of Propionibacterium freudenreichii CIRM-BIA1T adaptation strategies in conditions mimicking cheese ripening in the cold. PLoS One 7:e29083. doi: 10.1371/journal.pone.0029083
Damico, S., Collins, T., Marx, J.-C., Feller, G., and Gerday, C. (2006). Psychrophilic microorganisms: challenges for life. EMBO Rep. 7, 385–389. doi: 10.1038/sj.embor.7400662
Didelot, X., and Maiden, M. C. (2010). Impact of recombination on bacterial evolution. Trends Microbiol. 18, 315–322. doi: 10.1016/j.tim.2010.04.002
Didelot, X., and Wilson, D. J. (2015). ClonalFrameML: efficient inference of recombination in whole bacterial genomes. PLoS Comput. Biol. 11:e1004041. doi: 10.1371/journal.pcbi.1004041
Dieser, M., Greenwood, M., and Foreman, C. M. (2010). Carotenoid pigmentation in Antarctic heterotrophic bacteria as a strategy to withstand environmental stresses. FEMS Microbiol. Ecol. 42, 396–405. doi: 10.1657/1938-4246-42.4.396
Dong, K., Liu, H., Zhang, J., Zhou, Y., and Xin, Y. (2012). Flavobacterium xueshanense sp. nov. and Flavobacterium urumqiense sp. nov., two psychrophilic bacteria isolated from glacier ice. Int. J. Syst. Evol. Microbiol. 62, 1151–1157. doi: 10.1099/ijs.0.030049-0
Dsouza, M., Taylor, M. W., Turner, S. J., and Aislabie, J. (2015). Genomic and phenotypic insights into the ecology of Arthrobacter from Antarctic soils. BMC Genomics 16, 36–36. doi: 10.1186/s12864-015-1220-2
Feklistov, A., Sharon, B. D., Darst, S. A., and Gross, A. C. (2014). Bacterial sigma factors: a historical, structural, and genomic perspective. Annu. Rev. Microbiol. 68, 357–376. doi: 10.1146/annurev-micro-092412-155737
Feller, G. (2010). Protein stability and enzyme activity at extreme biological temperatures. J. Phys-Condens. Mat. 22:323101. doi: 10.1088/0953-8984/22/32/323101
Feller, G., and Gerday, C. (2003). Psychrophilic enzymes: hot topics in cold adaption. Nat. Rev. Microbiol. 1, 200–208. doi: 10.1038/nrmicro773
Fields, P. A. (2001). Review: protein function at thermal extremes: balancing stability and flexibility. Comp. Biochem. Physiol., Part A Mol. Integr. Physiol. 129, 417–431. doi: 10.1016/S1095-6433(00)00359-7
Fong, N., Burgess, M., Barrow, K., and Glenn, D. (2001). Carotenoid accumulation in the psychrotrophic bacterium Arthrobacter agilis in response to thermal and salt stress. Appl. Microbiol. Biotechnol. 56, 750–756. doi: 10.1007/s002530100739
Franzetti, A., Tatangelo, V., Gandolfi, I., Bertolini, V., Bestetti, G., Diolaiuti, G., et al. (2013). Bacterial community structure on two alpine debris-covered glaciers and biogeography of Polaromonas phylotypes. ISME J. 7, 1483–1492. doi: 10.1038/ismej.2013.48
Fuhrman, J. A., Schwalbach, M. S., and Stingl, U. (2008). Proteorhodopsins: an array of physiological roles? Nat. Rev. Microbiol. 6, 488–494. doi: 10.1038/nrmicro1893
Gao, H., Yang, Z. K., Wu, L., Thompson, D. K., and Zhou, J. (2006). Global transcriptome analysis of the cold shock response of Shewanella oneidensis MR-1 and mutational analysis of its classical cold shock proteins. J. Bacteriol. 188, 4560–4569. doi: 10.1128/JB.01908-05
Giovannoni, S. J., Bibbs, L., Cho, J.-C., Stapels, M. D., Desiderio, R., Vergin, K. L., et al. (2005). Proteorhodopsin in the ubiquitous marine bacterium SAR11. Nature 438, 82–85. doi: 10.1038/nature04032
Gómez-Consarnau, L., Gonzalez, J. M., Collllado, M., Gourdon, P., Pascher, T., Neutze, R., et al. (2007). Light stimulates growth of proteorhodopsin-containing marine Flavobacteria. Nature 445, 210–213. doi: 10.1038/nature05381
Gómez-Consarnau, L., Gonzalez, J. M., Riedel, T., Jaenicke, S., Wagner-Dobler, I., Sanudo-Wilhelmy, S. A., et al. (2016). Proteorhodopsin light-enhanced growth linked to vitamin-B1 acquisition in marine Flavobacteria. ISME J. 10, 1102–1112. doi: 10.1038/ismej.2015.196
Hille, F., and Charpentier, E. (2016). CRISPR-Cas: biology, mechanisms and relevance. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 371. doi: 10.1098/rstb.2015.0496
Hoffmann, T., and Bremer, E. (2011). Protection of Bacillus subtilis against cold stress via compatible-solute acquisition. J. Bacteriol. 193, 1552–1562. doi: 10.1128/JB.01319-10
Kasana, R. C., and Gulati, A. (2011). Cellulases from psychrophilic microorganisms: a review. J. Basic Microbiol. 51, 572–579. doi: 10.1002/jobm.201000385
Kawahara, H. (2002). The structures and functions of ice crystal-controlling proteins from bacteria. J. Biosci. Bioeng. 94, 492–496. doi: 10.1016/S1389-1723(02)80185-2
Kim, M., Oh, H. S., Park, S. C., and Chun, J. (2014). Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. Int. J. Syst. Evol. Microbiol. 64, 346–351. doi: 10.1099/ijs.0.059774-0
Kimura, M. (1980). A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol. 16, 111–120. doi: 10.1007/BF01731581
Kuhn, E. (2012). Toward understanding life under subzero conditions: the significance of exploring psychrophilic “cold-shock” proteins. Astrobiology 12, 1078–1086. doi: 10.1089/ast.2012.0858
Lami, R., Cottrell, M. T., Campbell, B. J., and Kirchman, D. L. (2009). Light-dependent growth and proteorhodopsin expression by Flavobacteria and SAR11 in experiments with Delaware coastal waters. Environ. Microbiol. 11, 3201–3209. doi: 10.1111/j.1462-2920.2009.02028.x
Lane, D. J. (1991). “16S/23S rRNA sequencing” in Nucleic acid techniques in bacterial systematics. ed. E. Stackebrandt (New York: John Wiley and Sons), 115–175.
Li, M., Wang, J., Zhang, Q., Lin, L., Kuang, A., Materon, L., et al. (2016). Isolation and characterization of the lytic cold-active bacteriophage MYSP06 from the Mingyong glacier in China. Curr. Microbiol. 72, 120–127. doi: 10.1007/s00284-015-0926-3
Liu, Q., Liu, H. C., Zhang, J. L., Zhou, Y. G., and Xin, Y. H. (2015a). Sphingomonas psychrolutea sp. nov., a psychrotolerant bacterium isolated from glacier ice. Int. J. Syst. Evol. Microbiol. 65, 2955–2959. doi: 10.1099/ijs.0.000362
Liu, Y., Vick-Majors, T. J., Priscu, J. C., Yao, T., Kang, S., Liu, K., et al. (2017). Biogeography of cryoconite bacterial communities on glaciers of the Tibetan plateau. FEMS Microbiol. Ecol. 93, fix072. doi: 10.1093/femsec/fix072
Liu, Y., Yao, T., Jiao, N., Kang, S., Xu, B., Zeng, Y., et al. (2009). Bacterial diversity in the snow over Tibetan plateau glaciers. Extremophiles 13, 411–423. doi: 10.1007/s00792-009-0227-5
Liu, Q., Zhou, Y. G., and Xin, Y. H. (2015b). High diversity and distinctive community structure of bacteria on glaciers in China revealed by 454 pyrosequencing. Syst. Appl. Microbiol. 38, 578–585. doi: 10.1016/j.syapm.2015.09.005
Mangiagalli, M., Bardolev, M., Tedesco, P., Natalello, A., Kaleda, A., Brocca, S., et al. (2017). Cryo-protective effect of an ice-binding protein derived from Antarctic bacteria. FEBS J. 284, 163–177. doi: 10.1111/febs.13965
Marshall, C. J. (1997). Cold-adapted enzymes. Trends Biotechnol. 15, 359–364. doi: 10.1016/S0167-7799(97)01086-X
McInerney, J. O., McNally, A., and O'Connell, M. J. (2017). Why prokaryotes have pangenomes. Nat. Microbiol. 2:17040. doi: 10.1038/nmicrobiol.2017.40
Methé, B. A., Nelson, K. E., Deming, J. W., Momen, B., Melamud, E., Zhang, X., et al. (2005). The psychrophilic lifestyle as revealed by the genome sequence of Colwellia psychrerythraea 34H through genomic and proteomic analyses. Proc. Natl. Acad. Sci. USA 102, 10913–10918. doi: 10.1073/pnas.0504766102
Mikheenko, A., Prjibelski, A. D., Saveliev, V., Antipov, D., and Gurevich, A. (2018). Versatile genome assembly evaluation with QUAST-LG. Bioinformatics 34, i142–i150. doi: 10.1093/bioinformatics/bty266
Mykytczuk, N. C., Foote, S. J., Omelon, C. R., Southam, G., Greer, C. W., and Whyte, L. G. (2013). Bacterial growth at −15°C; molecular insights from the permafrost bacterium Planococcus halocryophilus Or1. ISME J. 7, 1211–1226. doi: 10.1038/ismej.2013.8
Nguyen, L. T., Schmidt, H. A., von Haeseler, A., and Minh, B. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274. doi: 10.1093/molbev/msu300
Nichols, C. A., Guezennec, J., and Bowman, J. P. (2005). Bacterial exopolysaccharides from extreme marine environments with special consideration of the southern ocean, sea ice, and deep-sea hydrothermal vents: a review. Mar. Biotechnol. 7, 253–271. doi: 10.1007/s10126-004-5118-2
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114
Phadtare, S. (2004). Recent developments in bacterial cold-shock response. Curr. Issues Mol. Biol. 6, 125–136.
Polz, M. F., Alm, E. J., and Hanage, W. P. (2013). Horizontal gene transfer and the evolution of bacterial and archaeal population structure. Trends Genet. 29, 170–175. doi: 10.1016/j.tig.2012.12.006
Qin, Q., Xie, B., Yu, Y., Shu, Y. L., Rong, J. C., Zhang, Y. J., et al. (2014). Comparative genomics of the marine bacterial genus Glaciecola reveals the high degree of genomic diversity and genomic characteristic for cold adaptation. Environ. Microbiol. 16, 1642–1653. doi: 10.1111/1462-2920.12318
Raymond, J. A., Christner, B. C., and Schuster, S. C. (2008). A bacterial ice-binding protein from the Vostok ice core. Extremophiles 12, 713–717. doi: 10.1007/s00792-008-0178-2
Raymond, J. A., Fritsen, C., and Shen, K. (2007). An ice-binding protein from an Antarctic Sea ice bacterium. FEMS Microbiol. Ecol. 61, 214–221. doi: 10.1111/j.1574-6941.2007.00345.x
Reichenbach, H., Kohl, W., Bttgervetter, A., and Achenbach, H. (1980). Flexirubin-type pigments in Flavobacterium. Arch. Microbiol. 126, 291–293. doi: 10.1007/BF00409934
Reid, I. N., Sparks, W. B., Lubow, S., Mcgrath, M., Livio, M., Valenti, J., et al. (2006). Terrestrial models for extraterrestrial life: methanogens and halophiles at Martian temperatures. Int. J. Aatrobiol. 5, 89–97. doi: 10.1017/S1473550406002916
Richter, M., and Rosselló-Móra, R. (2009). Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. U. S. A. 106, 19126–19131. doi: 10.1073/pnas.0906412106
Riley, M., Staley, J. T., Danchin, A., Wang, T. Z., Brettin, T. S., Hauser, L. J., et al. (2008). Genomics of an extreme psychrophile, Psychromonas ingrahamii. BMC Genomics 9:210. doi: 10.1186/1471-2164-9-210
Sasser, M. (1990). Identification of bacteria by gas chromatography of cellular fatty acids, MIDI technical note 101. Newark, DE: MIDI Inc. Placeholder Text.
Saunders, N. F., Thomas, T., Curmi, P. M., Mattick, J. S., Kuczek, E., Slade, R., et al. (2003). Mechanisms of thermal adaptation revealed from the genomes of the Antarctic Archaea Methanogenium frigidum and Methanococcoides burtonii. Genome Res. 13, 1580–1588. doi: 10.1101/gr.1180903
Sawstrom, C., Graneli, W., Laybournparry, J., and Anesio, A. M. (2007). High viral infection rates in Antarctic and Arctic bacterioplankton. Environ. Microbiol. 9, 250–255. doi: 10.1111/j.1462-2920.2006.01135.x
Sievers, F., and Higgins, D. G. (2014). Clustal omega, accurate alignment of very large numbers of sequences. Methods Mol. Biol. 1079, 105–116. doi: 10.1007/978-1-62703-646-7_6
Smibert, R. M., and Krieg, N. R. (1994). “Phenotypic characterization” in Methods for general and molecular bacteriology. eds. P. Gerhardt, R. G. E. Murray, W. A. Wood, and N. R. Krieg (Washington, DC: American Society for Microbiology), 607–654.
Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28, 2731–2739. doi: 10.1093/molbev/msr121
Thompson, J. D., Higgins, D. G., and Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673–4680. doi: 10.1093/nar/22.22.4673
Tindall, B. J., Sikorski, J., Smibert, R. M., and Kreig, N. R. (2007). “Phenotypic characterization and the principles of comparative systematics” in Methods for general and molecular microbiology. eds. C. A. Reddy, T. J. Beveridge, J. A. Breznak, G. Marzluf, and T. M. Schmidt (Washington, DC: American Society for Microbiology), 330–393.
Vinuesa, P., Ochoa-Sánchez, L. E., and Contreras-Moreira, B. (2018). GET_PHYLOMARKERS, a software package to select optimal orthologous clusters for phylogenomics and inferring pan-genome phylogenies, used for a critical geno-taxonomic revision of the genus Stenotrophomonas. Front. Microbial. 9, 771. doi: 10.3389/fmicb.2018.00771
Vos, M., and Didelot, X. (2009). A comparison of homologous recombination rates in bacteria and archaea. ISME J. 3, 199–208. doi: 10.1038/ismej.2008.93
Vos, M., Hesselman, M. C., te Beek, T. A., van Passel, M. W. J., and Eyre-Walker, A. (2015). Rates of lateral gene transfer in prokaryotes: high but why? Trends Microbiol. 23, 598–605. doi: 10.1016/j.tim.2015.07.006
Wang, F., Wang, J., Jian, H., Zhang, B., Li, S., Zeng, X., et al. (2008). Environmental adaptation: genomic analysis of the piezotolerant and psychrotolerant deep-sea iron reducing bacterium Shewanella piezotolerans WP3. PLoS One 3:e1937. doi: 10.1371/journal.pone.0001937
Welsh, D. T. (2000). Ecological significance of compatible solute accumulation by micro-organisms: from single cells to global climate. FEMS Microbiol. Rev. 24, 263–290. doi: 10.1111/j.1574-6976.2000.tb00542.x
Xin, Y., Liang, Z., Zhang, D., Liu, H., Zhang, J., Yu, Y., et al. (2009). Flavobacterium tiangeerense sp. nov., a cold-living bacterium isolated from a glacier. Int. J. Syst. Evol. Microbiol. 59, 2773–2777. doi: 10.1099/ijs.0.007906-0
Xu, M., Xin, Y., Tian, J., Dong, K., Yu, Y., Zhang, J., et al. (2011). Flavobacterium sinopsychrotolerans sp. nov., isolated from a glacier. Int. J. Syst. Evol. Microbiol. 61, 20–24. doi: 10.1099/ijs.0.014126-0
Yoon, S. H., Ha, S. M., Kwon, S., Lim, J., Kim, Y., Seo, H., et al. (2017). Introducing EzBioCloud: a taxonomically united database of 16S rRNA gene sequences and whole-genome assemblies. Int. J. Syst. Evol. Microbiol. 67, 1613–1617. doi: 10.1099/ijsem.0.001755
Zhang, Y., Jiang, F., Chang, X., Qiu, X., Ren, L., Qu, Z., et al. (2016). Flavobacterium collinsense sp. nov., isolated from a till sample of an Antarctic glacier. Int. J. Syst. Evol. Microbiol. 66, 172–177. doi: 10.1099/ijsem.0.000688
Zhang, D., Wang, H., Liu, H., Dong, X., and Zhou, P. (2006). Flavobacterium glaciei sp. nov., a novel psychrophilic bacterium isolated from the China no.1 glacier. Int. J. Syst. Evol. Microbiol. 56, 2921–2925. doi: 10.1099/ijs.0.64564-0
Zhaxybayeva, O., and Doolittle, W. F. (2011). Lateral gene transfer. Curr. Biol. 21, R242–R246. doi: 10.1016/j.cub.2011.01.045
Zhu, L., Liu, Q., Liu, H., Zhang, J., Dong, X., Zhou, Y., et al. (2013). Flavobacterium noncentrifugens sp. nov., a psychrotolerant bacterium isolated from glacier meltwater. Int. J. Syst. Evol. Microbiol. 63, 2032–2037. doi: 10.1099/ijs.0.045534-0
Keywords: Flavobacterium bomense, glacier-inhabiting bacteria, microevolution, cold adaptation, psychrophilic
Citation: Liu Q, Liu H-C, Zhou Y-G and Xin Y-H (2019) Microevolution and Adaptive Strategy of Psychrophilic Species Flavobacterium bomense sp. nov. Isolated From Glaciers. Front. Microbiol. 10:1069. doi: 10.3389/fmicb.2019.01069
Edited by:
Brian D. Lanoil, University of Alberta, CanadaReviewed by:
Alexandre Anesio, Aarhus University, DenmarkCamilla Lothe Nesbø, University of Alberta, Canada
Copyright © 2019 Liu, Liu, Zhou and Xin. 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: Yu-Guang Zhou, zhouyg@im.ac.cn
Yu-Hua Xin, xinyh@im.ac.cn