CVTree3 Web Server for Whole-genome-based and Alignment-free Prokaryotic Phylogeny and Taxonomy

CVTree: A Parallel Alignment-free Phylogeny and Taxonomy Tool based on Composition Vectors of Genomes

CVTree is an alignment-free algorithm to infer phylogenetic relationships from genome sequences. It had been successfully applied to study phylogeny and taxonomy of viruses, prokaryotes, and fungi based on the whole genomes, as well as chloroplasts, mitochondria, and metagenomes. Here we presented the standalone software for the CVTree algorithm. In the software, an extensible parallel workflow for the CVTree algorithm was designed. Based on the workflow, new alignment-free methods were also implemented. And by examining the phylogeny and taxonomy of 13,903 prokaryotes based on 16S rRNA sequences, we showed that CVTree software is an efficient and effective tool for the studying of phylogeny and taxonomy based on genome sequences. Code availability: https://github.com/ghzuo/cvtree.

Bioinformatic and Functional Characterization of Hsp70s in Myxococcus xanthus

Hsp70 proteins are among the most ubiquitous chaperones and play important roles in maintaining proteostasis and resisting environmental stress. Multiple copies of Hsp70s are widely present in eukaryotic cells with redundant and divergent functions, but they have been less well investigated in prokaryotes. Myxococcus xanthus DK1622 is annotated as having many hsp70 genes. In this study, we performed a bioinformatic analysis of Hsp70 proteins and investigated the functions of six hsp70 genes in DK1622, including two genes that encode proteins with the conserved PRK00290 domain (MXAN_3192 and MXAN_6671) and four genes that encode proteins with the cl35085 or cd10170 domain. We found that only MXAN_3192 is essential for cell survival and heat shock induction. MXAN_3192, compared with the other hsp70 genes, has a high transcriptional level, far exceeding that of any other hsp70 gene, which, however, is not the reason for its essentiality. Deletion of MXAN_6671 (sglK) led to multiple deficiencies in development, social motility, and oxidative resistance, while deletion of each of the other four hsp70 genes decreased sporulation and oxidative resistance. MXAN_3192 or sglK, but not the other genes, restored the growth deficiency of the E. colidnaK mutant. Our results demonstrated that the PRK00290 proteins play a central role in the complex cellular functions of M. xanthus, while the other diverse Hsp70 superfamily homologues probably evolved as helpers with some unknown specific functions.

IMPORTANCE Hsp70 proteins are highly conserved chaperones that occur in all kingdoms of life. Multiple copies of Hsp70s are often present in genome-sequenced prokaryotes, especially taxa with complex life cycles, such as myxobacteria. We investigated the functions of six hsp70 genes in Myxococcus xanthus DK1622 and demonstrated that the two Hsp70 proteins with the PRK00290 domain play a central role in complex cellular functions in M. xanthus, while other Hsp70 proteins probably evolved as helpers with some unknown specific functions.

Acidovorax antarcticus sp. nov., isolated from a soil sample of Collins Glacier front, Antarctica

A Gram-stain-negative, rod-shaped, and aerobic bacterium, strain 16-35-5^T, was isolated from Collins Glacier front soil from the Fildes Peninsula, Antarctica. The bacterium grew optimally at 28 °C, pH 7.0 and in the presence of 0-4.0 % (w/v) NaCl. On the basis of the results of 16S rRNA gene sequence phylogenetic analyses, it was concluded that 16-35-5^T represented a member of the genus Acidovorax and had the highest sequence similarities with Acidovorax anthurii CFBP 3232^T (96.48 %). The genome of 16-35-5^T is 4.2 Mb long with a DNA G+C content of 66.3 mol%. Average nucleotide identity (ANI) value between the genomes of 16-35-5^T and Acidovorax wautersii DSM 27981^T, was 85.29 %. Strain 16-35-5^T had ubiquinone-8 (Q-8) as the respiratory ubiquinone. The polar lipids of 16-35-5^T were consisted of phosphatidylethanolamine, diphosphatidylglycerol and phosphatidylglycerol. The main fatty acids were summed feature 3 (C_16 : 1ω7c and/or C_16 : 1ω6c, 25.2 %), summed feature 8 (C_18 : 1ω7c and/or C_18 : 1ω6c, 12.9 %), C_16 : 0 (35.2 %), and C_17 : 0 cyclo (19.0 %). On the basis of the evidence presented in this study, 16-35-5^T should be classified as representing a novel species of the genus Acidovorax, for which the name Acidovorax antarcticus sp. nov., is proposed, with the type strain 16-35-5^T (=CCTCC AB 2019325^T=KCTC 72915^T).

Kaistella flava sp. nov., isolated from Antarctic tundra soil, and emended descriptions of Kaistella yonginensis, Kaistella jeonii, Kaistella antarctica and Kaistella chaponensis

A rod-shaped, yellow-pigmented, Gram-stain-negative, non-motile and aerobic bacterium, designated 7-3A^T, was isolated from soil from King George Island, maritime Antarctica, and subjected to a polyphasic taxonomic study. Growth occurred at 4-37 °C (optimum, 20°C) and at pH 5.0-9.0 (optimum, pH 7.0-8.0). Tolerance to NaCl was up to 4 % (w/v) with optimum growth in the absence of NaCl. The results of phylogenetic analysis based on 16S rRNA gene sequences indicated that strain 7-3A^T represented a member of the family Flavobacteriaceae. Strain 7-3A^T showed the highest sequence similarities with Kaistella yonginensis HMD 1043^T (96.65 %), Kaistella carnis NCTC 13525^T (96.53 %), Kaistella chaponensis DSM 23145^T (96.27 %), Kaistella antarctica LMG 24720^T (96.13 %) and Kaistella jeonii DSM 17048^T (96.06 %). A whole genome-level comparison of 7-3A^T with K. jeonii DSM 17048^T, K. antarctica LMG 24720^T, K. chaponensis DSM 23145^T, and Kaistella palustris DSM 21579^T revealed average nucleotide identity (ANI) values of 79.03, 82.25, 78.12, and 74.42 %, respectively. The major respiratory isoprenoid quinone was identified as MK-6 and a few ubiquinones Q-10 were identified. In addition, flexirubin-type pigments were absent. The polar lipid profile of 7-3A^T was found to contain one phosphatidylethanolamine, six unidentified aminolipids (AL) and two unidentified lipids (L). The G+C content of the genomic DNA was determined to be 34.54 mol%. The main fatty acids were iso-C_15 : 0, summed feature 9 (comprising iso-C_17 : 1ω9c and/or C_16 : 0 10-methyl), anteiso-C_15 : 0, iso-C_13 : 0 and summed feature 3 (comprising C_16 : 1ω7c and/or C_16 : 1ω6c). On the basis of the evidence presented in this study, a novel species of the genus Kaistella, Kaistella flava sp. nov., is proposed, with the type strain 7-3A^T (=CCTCC AB 2016141^T= KCTC 52492^T). Emended descriptions of Kaistella yonginensis, Kaistella jeonii, Kaistella antarctica and Kaistella chaponensis are also given.

Aerophototrophica crusticola gen. nov., sp. nov., isolated from desert biocrusts

A pink-pigmented, Gram-stain-negative, rod-shaped, strictly aerobic bacterial strain MIMtkB3^T, was isolated from moss crusts in Hunshandake desert of China. Cells grew at 15-45 °C (optimum of 28 °C), at pH of 6.0-8.5 (optimum of 7.0) and with 0-1.0 % (w/v) NaCl (optimum of 0 %). The strain could biosynthesize the green-coloured pigment bacteriochlorophyll a (BChl a). The respiratory quinone was ubiquinone Q-10, while C_18 : 1 ω7c and C_18 : 1 2OH were the major fatty acids. Phosphatidylethanolamine, phosphatidylglycerol, diphosphatidylglycerol, an unidentified aminophospholipid, one unidentified phospholipid, three unidentified glycolipid and one unidentified lipid were the major polar lipids. Strain MIMtkB3^T was most closely related to Oleisolibacter albus NAU-10^T, Niveispirillum fermenti CC-LY736^T, and Rhodocista centenaria SW of the family Rhodospirillaceae with 16S rRNA gene similarities of 93.09, 92.02 and 91.73%, respectively. The genomic DNA G+C content calculated on complete genome sequencing was 69.3 mol%. The average nucleotide identity between strain MIMtkB3^T and its closely related type strains in Rhodospirillaceae was below 77.96 % and digital DNA-DNA hybridization lower than 24.70 %. Full light utilization pathway of aerobic anoxygenic phototrophic bacteria was identified in the genome. Based on phenotypic, chemotaxonomic and phylogenetic characteristics, strain MIMtkB3^T represents a novel genus of the family Rhodospirillaceae, for which the name Aerophototrophica crusticola gen. nov., sp. nov. is proposed. The type strain is MIMtkB3^T (=KCTC 42633^T=MCCC 1K00570^T).

Comparative genomics of the Pasteurella multocida toxin

Pasteurella multocida is a zoonotic pathogen whose genetic heterogeneity is well known. Five serogroups (A, B, D, E, and F) and 16 serotypes of P. multocida have been recognized thus far based on capsular polysaccharide typing and lipopolysaccharide typing, respectively. Progressive atrophic rhinitis in domestic pigs is caused by P. multocida strains containing toxA, which encodes a 146 kDa heat-labile toxin. Among the five serogroups, only some strains of serogroups A and D are toxigenic. In this study, by comparative analyses of the genomes of many strains, it has been shown that toxA is sparsely distributed in P. multocida. Furthermore, full-length homologs of P. multocida toxA were found only in two other bacterial species. It has also been shown that toxA is usually associated with a prophage, and that some strains contain an orthologous prophage but not toxA. Among the toxA-containing prophages that were compared, an operon putatively encoding a type II restriction-modification system was present only in strains LFB3, HN01, and HN06. These results indicate that the selection and maintenance of the heat-labile toxin and the type II restriction-modification system are evolutionarily less favorable among P. multocida strains. Phylogenetic analysis using the alignment- and parameter-free method CVTree3 showed that deduced proteome sequences can be used as effectively as whole/core genome single nucleotide polymorphisms to group P. multocida strains in relation to their serotypes and (or) genotypes. It remains to be determined if the toxA-containing prophages in strains HN01 and HN06 are inducible, and if they can be used for lysogenic transfer of toxA to other bacteria.

A comprehensive analysis of the phylogenetic signal in ramp sequences in 211 vertebrates

Ramp sequences increase translational speed and accuracy when rare, slowly-translated codons are found at the beginnings of genes. Here, the results of the first analysis of ramp sequences in a phylogenetic construct are presented. Ramp sequences were compared from 247 vertebrates (114 Mammalian and 133 non-mammalian), where the presence and absence of ramp sequences was analyzed as a binary character in a parsimony and maximum likelihood framework. Additionally, ramp sequences were mapped to the Open Tree of Life synthetic tree to determine the number of parallelisms and reversals that occurred, and those results were compared to random permutations. Parsimony and maximum likelihood analyses of the presence and absence of ramp sequences recovered phylogenies that are highly congruent with established phylogenies. Additionally, 81% of vertebrate mammalian ramps and 81.2% of other vertebrate ramps had less parallelisms and reversals than the mean from 1000 randomly permuted trees. A chi-square analysis of completely orthologous ramp sequences resulted in a p-value < 0.001 as compared to random chance. Ramp sequences recover comparable phylogenies as other phylogenomic methods. Although not all ramp sequences appear to have a phylogenetic signal, more ramp sequences track speciation than expected by random chance. Therefore, ramp sequences may be used in conjunction with other phylogenomic approaches if many orthologs are taken into account. However, phylogenomic methods utilizing few orthologs should be cautious in incorporating ramp sequences because individual ramp sequences may provide conflicting signals.

Biocontrol Using Bacillus amyloliquefaciens PP19 Against Litchi Downy Blight Caused by Peronophythora litchii

Bacillus amyloliquefaciens has been widely used in the agriculture, food, and medicine industries. Isolate PP19 was obtained from the litchi fruit carposphere and showed biocontrol efficacy against litchi downy blight (LDB) whether applied preharvest or postharvest. To further understand the underlying regulatory mechanisms, the genome of PP19 was sequenced and analyzed. The genome comprised a 3,847,565 bp circular chromosome containing 3990 protein-coding genes and 121 RNA genes. It has the smallest genome among 36 sequenced strains of B. amyloliquefaciens except for RD7-7. In whole genome phylogenetic analysis, PP19 was clustered into a group with known industrial applications, indicating that it may also produce high-yield metabolites that have yet to be identified. A large chromosome structural variation and large numbers of single nucleotide polymorphisms (SNPs) between PP19 (industrial strain) and UMAF6639 (plant-associated strain) were detected through comparative analysis, which may shed light on their functional differences. Preharvest treatment with PP19 enhanced resistance to LDB, by decreasing the plant H2O2 content and increasing the SOD activity. This is the first report of an industrial strain of B. amyloliquefaciens showing a plant-associated function and with major potential for the biocontrol of LDB.

Isolation and characterization of the novel Pseudomonas stutzeri bacteriophage 8P

Bacteriophage 8P was isolated with a Pseudomonas stutzeri strain isolated from an oil reservoir as its host bacterium. The phage genome comprises 63,753 base pairs with a G+C content of 64.35. The phage encodes 63 predicted proteins, and 27 of them were functionally assigned. No tRNA genes were found. Comparative genomics analysis showed that 8P displayed some relatedness to F116-like phages (78% identity, 20% query coverage). The genome has very low sequence similarity to the other phage genomes in the GenBank database and Viral Sequence Database. Based on whole-genome analysis and transmission electron microscopy imaging, 8P is proposed to be a member of a new species in the genus Hollowayvirus, family Podoviridae.

Phylogenomic analysis of Clostridioides difficile ribotype 106 strains reveals novel genetic islands and emergent phenotypes

Clostridioides difficile infection (CDI) is a major healthcare-associated diarrheal disease. Consistent with trends across the United States, C. difficile RT106 was the second-most prevalent molecular type in our surveillance in Arizona from 2015 to 2018. A representative RT106 strain displayed robust virulence and 100% lethality in the hamster model of acute CDI. We identified a unique 46 KB genomic island (GI1) in all RT106 strains sequenced to date, including those in public databases. GI1 was not found in its entirety in any other C. difficile clade, or indeed, in any other microbial genome; however, smaller segments were detected in Enterococcus faecium strains. Molecular clock analyses suggested that GI1 was horizontally acquired and sequentially assembled over time. GI1 encodes homologs of VanZ and a SrtB-anchored collagen-binding adhesin, and correspondingly, all tested RT106 strains had increased teicoplanin resistance, and a majority displayed collagen-dependent biofilm formation. Two additional genomic islands (GI2 and GI3) were also present in a subset of RT106 strains. All three islands are predicted to encode mobile genetic elements as well as virulence factors. Emergent phenotypes associated with these genetic islands may have contributed to the relatively rapid expansion of RT106 in US healthcare and community settings.

Pseudopuniceibacterium antarcticum sp. nov., isolated from an Antarctic marine sponge

A Gram-negative, aerobic, rod-shaped, non-motile bacterium, designated strain HQ09^T, was isolated from a marine sponge off the coast of Fields Peninsula, West Antarctica. Strain HQ09^T grew at 4-35 °C (optimum, 25 °C), pH 5-9 (optimum, pH 7.0), and with 1-10% NaCl (optimum, 2 %). Phylogenetic analysis based on the 16S rRNA gene sequences showed that strain HQ09^T was affiliated with the genus Pseudopuniceibacterium in the family Rhodobacteraceae, sharing 99.64 % identity with the type strain of Pseudopuniceibacterium sediminis, the only known species in the genus. However, the low digital DNA-DNA hybridization (dDDH) (27.2 %) and average nucleotide identity (ANI) (83.63 %) values between strain HQ09^T and the type strain of Pseudopuniceibacterium sediminis indicated that they did not belong to the same species. Strain HQ09^T could also be differentiated from Pseudopuniceibacterium sediminis by many phenotypic characteristics. The major fatty acids (>5 %) of strain HQ09^T were summed feature 8 (C_18 : 1  ω7c/C_18 : 1 ω6c), 11-methyl C_18 : 1  ω7c, C_16 : 0 and C19 : 0 cyclo ω8c. The polar lipids included phosphatidylglycerol, phosphatidylcholine, two unidentified aminolipids and one unidentified phospholipid. The predominant respiratory quinone was ubiquinone 10 (Q-10). The genomic DNA G+C content was 62.63 mol%. Four secondary metabolite biosynthetic gene clusters were detected in the genome, potentially producing ectoine and three types of unknown compounds. On the basis of the polyphasic evidences obtained in this study, strain HQ09^T represents a novel species of the genus Pseudopuniceibacterium, for which the name Pseudopuniceibacterium antarcticum sp. nov. is proposed, with the type strain being HQ09^T (=KCTC 52229^T=CGMCC 1.15538^T).

Noviherbaspirillum aerium sp. nov., isolated from air

A Gram-stain-negative, strictly aerobic, non-spore-forming, rod-shaped, motile with polar flagella and pale-orange bacterium, designated strain 122213-3^T, was isolated from air, collected at the foot of the Xiangshan Mountain, located in Beijing, PR China. Optimal growth occurred at 28 °C, at pH 7 and in the presence of 0-1 % (w/v) NaCl. Phylogenetic analyses based on 16S rRNA gene sequences revealed that 122213-3^T clustered with species of the genus Noviherbaspirillum and formed a distinct sublineage, showing highest similarities to Noviherbaspirillum malthae CC-AFH3^T (96.88 %), Noviherbaspirillum massiliense JC206^T (95.78 %) and Noviherbaspirillum aurantiacum SUEMI08^T (95.78 %). The predominant cellular fatty acids were summed feature 3 (C_16 : 1  ω6c and/or C_16 : 1  ω7c), summed feature 8 (C_18 : 1  ω7c and/or C_18 : 1  ω6c) and C_16 : 0. The predominant quinone was ubiquinone 8 (Q-8). The polar lipids comprised phosphatidylethanolamine, phosphatidylglycerol, unidentified phospholipid and two unidentified polar lipids. The polyamine pattern showed the presence of putrescine as the major polyamine, with minor amounts of 2-hydroxyputrescine. The DNA G+C content was 60.1 mol%. The phylogenetic analysis and physiological and biochemical data showed that strain 122213-3^T should be classified as representing a novel species in the genus Noviherbaspirillum, for which the name Noviherbaspirillum aerium sp. nov. is proposed. The type strain of N. aerium is 122213-3^T (=CFCC 14286^T=LMG 30131^T).

Genomic analysis of Ranavirus and exploring alternative genes for phylogenetics

Ranaviruses can infect both captive and wild cold-blooded vertebrates, leading to significant economic and environmental losses. With the cases of ranavirus infection increasing, many ranavirus genomic sequences were published, but little is known about ranavirus taxonomy on a whole-genome level. In this study, 44 ranaviruses core genes were identified in 32 ranaviruses genome sequences by using PanX. The neighbour-joining phylogenetic trees (NJ-tree) based on 44 ranaviruses core genes and 24 iridoviridae core genes and composition vector phylogenetic tree (CV-Tree) based on whole genome were constructed. The three of phylogenetic trees showed that 32 ranavirus isolates can be divided into 4 different subgroups including SGIV-like, EHNV-like, FV3-like and CMTV-like, and subgroups taxonomic position of three phylogenetic trees were consistent. However, the phylogenetic position of ToRV could not be determined if it belongs to FV3-like or CMTV-like group. Subsequently, we carried out dot plot analysis and confirmed that ToRV should belong to CMTV-like group. Based on dot plot analysis and phylogenetic trees, the taxonomic classification of ranaviruses was confirmed. Finally, four genes which are suitable for the construction of phylogenetic tree were selected from ranavirus core genes by recombination analysis, substitution saturation analysis and single-gene phylogenetic analysis. Phylogenetic tree based on concatenated sequences of the four selected genes showed that the classification of subgroups was identical with three of the phylogenetic trees. Conclusion: Our results confirmed taxonomic identification of ranaviruses; the four selected genes used in phylogenetic analysis will make taxonomic identification more convenient and accurate.

Streptomyces arboris sp. nov., isolated from Populus euphratica wetland soil

A novel actinobacterial strain (TRM 68085T) was isolated from soil ofPopulus euphraticawetland. A polyphasic approach was used to study the taxonomy of TRM 68085Tand the results showed a range of phylogenetic and chemotaxonomic properties consistent with those of the genusStreptomyces. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain TRM 68085Tshowed the highest similarity value toStreptomyces capitiformicae1H-SSA4T(98.6 %), and phylogenetically clustered withStreptomyces kanasensisZX01T(97.5 %) andStreptomyces ipomoeaeNBRC 13050T(97.4 %). The genomic DNA G+C content of strain TRM 68085Tbased on the genome sequence was 71.4 mol%. The levels of DNA–DNA relatedness between the genome of the isolate and its nearest phylogenetic neighbours,S. capitiformicae1H-SSA4T,S. kanasensisZX01TandS. ipomoeaeNBRC 13050T, were 19.2±0.4, 21.8±0.5 and 19.3±0.6 %, respectively. Chemotaxonomic data revealed that strain TRM 68085Tpossessed MK-9(H6) and MK-9(H8) as the predominant menaquinones.ll-Diaminopimelic acid and a small amount ofmeso-diaminopimelic acid were the diagnostic diamino acids. Ribose, xylose, glucose and galactose were the whole-cell sugars. The major cellular fatty acids were C16 : 0(25.4 %) and iso-C16 : 0(18.3 %). On the basis of these genotypic and phenotypic data, it is concluded that strain TRM 68085Trepresents a novel species of the genusStreptomyces, for which the nameStreptomyces arborissp. nov. is proposed. The type strain is TRM 68085T(=CCTCC AA2019031T=LMG 31492T).

CGIDLA:Developing the Web Server for CpG Island Related Density and LAUPs (Lineage-Associated Underrepresented Permutations) Study

It is well known that CpG island plays an important role in gene methylation. Since CpG island is closely related to human genetic characteristics such as TATA-box, tissue expression specificity, and LAUPs (Lineage-associated Underrepresented Permutations), it is important to investigate the sequence specificity of CpG island as well as the potential genetic characteristics related to CpG island to further understand the methylation related regulation mechanism. Therefore, this study develops such an online service website for CpG island related density and LAUPs analysis (CGIDLA, www.combio-lezhang.online/cgidla/index.html), that not only can investigate the relationship among the CpG island density, TATA-box feature, and expression breadth of human genes, but also deposit LAUPs of 32 representative species to help molecular biologists investigate the relationship between CpG island and LUAPs. Moreover, CGIDLA provides the source code download service and the related LAUPs counting functions.

Tabrizicola algicola sp. nov. isolated from culture of microalga Ettlia sp

A novel Gram-stain-negative, aerobic, non-spore-forming, non-motile, and rod-shaped bacterium, strain ETT8^T was isolated from a chemostat culture of microalga Ettlia sp. YC001. Optimal growth was with 0-2% NaCl and at 25-37 °C on R2A medium. Phylogenetic analysis based on the 16S rRNA gene and genome sequence showed that strain ETT8^T belongs to the genus Tabrizicola, with the close neighbours being T. sediminis DRYC-M-16^T (98.1 %), T. alkalilacus DJC^T (97.6 %), T. fusiformis SY72^T (96.9 %), T. piscis K13M18^T (96.8 %), and T. aquatica RCRI19^T (96.5 %). The genomic comparison of strain ETT8^T with type species in the genus Tabrizicola was analysed using the genome-to-genome distance calculator (GGDC), average nucleotide identity (ANI), and average amino acid identity (AAI) (values indicated ≤17.7, ≤75.4 and ≤71.9 %, respectively). The genomic DNA G+C content of strain ETT8^T was 64.4 %, plus C_18 : 1  ω6c and C_18 : 0-iso were the major fatty acids and Q-10 the major respiratory quinone. Strain ETT8^T contained diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, phosphatidylcholine aminolipid, and four unidentified lipids as the major polar lipids. Based on the chemotaxonomic, genotypic, and phenotype results, strain ETT8^T was recognized as a novel species of the genus Tabrizicola for which the name Tabrizicola algicola sp. nov. is proposed. The type strain is ETT8^T (=KCTC 72206^T=JCM 31893^T=MCC 4339^T).

Pelagovum pacificum gen. nov., sp. nov., a novel member of the family Rhodobacteraceae isolated from surface seawater of the Mariana Trench

A Gram-stain-negative, aerobic, ovoid-rod-shaped bacterium, designated strain SM1903^T, was isolated from surface seawater of the Mariana Trench. The strain grew at 15-37 °C (optimum, 35 °C) and with 1-15 % (optimum, 4 %) NaCl. It hydrolysed aesculin but did not reduce nitrate to nitrite and hydrolyse Tween 80. Phylogenetic analysis based on the 16S rRNA gene sequences revealed that strain SM1903^T formed a separate lineage within the family Rhodobacteraceae, sharing the highest 16S rRNA gene sequence similarity with type strains of Pseudooceanicola antarcticus (95.7 %) and Roseisalinus antarcticus (95.7 %). In phylogenetic trees based on single-copy OCs and whole proteins sequences, strain SM1903^T fell within a sub-cluster encompassed by Oceanicola granulosus, Roseisalinus antarcticus and Histidinibacterium lentulum and formed a branch adjacent to Oceanicola granulosus. The major cellular fatty acids were summed feature 8 (C_18 : 1  ω7c and/or C_18 : 1  ω6c), C_16 : 0 and 11-methyl-C_18 : 1  ω7c. The polar lipids mainly comprised phosphatidylglycerol, phosphatidylcholine, one unidentified lipid, one unidentified aminolipid, and one unidentified glycolipid. The solo respiratory quinone was ubiquinone-10. The genomic DNA G+C content of strain SM1903^T was 66.0 mol%. Based on the results of phenotypic, chemotaxonomic, and phylogenetic characterization for strain SM1903^T, it is considered to represent a novel species of a novel genus in the family Rhodobacteraceae, for which the name Pelagovum pacificum gen. nov., sp. nov. is proposed. The type strain is SM1903^T (=MCCC 1K03608^T=KCTC 72046^T).

Rufibacter hautae sp. nov., a red-pigmented bacterium from freshwater lake sediment, and proposal of Rufibacter quisquiliarum as a latter heterotypic synonym of Rufibacter ruber

A taxonomic identification using a polyphasic approach was performed on strain NBS58-1^T, which was isolated from the interfacial sediment of Taihu Lake in China. Strain NBS58-1^T was Gram-stain-negative, aerobic, non-spore-forming and catalase-positive. Phylogenetic analyses based on 16S rRNA gene and three housekeeping genes (rpoB, gyrB and dnaK) sequences supported the position that strain NBS58-1^T should be classified within the genus Rufibacter. The 16S rRNA gene sequence of strain NBS58-1^T possessed the highest similarity to Rufibacter sediminis H-1^T (96.60 %), followed by Rufibacter glacialis MDT1-10-3^T (96.17 %). And the ANI value between strain NBS58-1^T and R. glacialis MDT1-10-3^T was 79.3 %. The respiratory quinone was menaquinone 7 (MK-7). The major cellular fatty acids comprised iso-C_15 : 0 and summed feature 3. Phosphatidylethanolamine, two unidentified phospholipids and four unidentified lipids were the main polar lipids. The genomic DNA G+C content was 51.3 mol%. Based on phenotypic features and phylogenetic position, a novel species with the name Rufibacter hautae sp. nov. is proposed. The type strain is NBS58-1^T=(KACC 21309^T=MCCC 1K04037^T). We also proposed Rufibacter quisquiliarum as a latter heterotypic synonym of Rufibacter ruber.

Novosphingobium aquimarinum sp. nov., isolated from seawater

A novel Gram-stain-negative, aerobic, and rod-shaped bacterial strain, M24A2M^T, was isolated from seawater in the Republic of Korea. On the basis of the 16S rRNA gene phylogeny, strain M24A2M^T was found to be closely related to Novosphingobium pentaromativorans US6-1^T and Novosphingobium mathurense SM117^T with pair-wise sequence similarities of 97.4 and 96.9 %, respectively. Phylogenetic analysis of 16S rRNA sequences indicated that M24A2M^T formed a branch with Novosphingobium pentaromativorans US6-1^T and represented a member of the genus Novosphingobium. The predominant cellular fatty acids were C_14 : 0 2-OH, summed feature 3 (C_16 : 1 ω7c and/or C_16 : 1 ω6c), and summed feature 8 (C_18 : 1 ω7c and/or C_18 : 1 ω6c). The polar lipids of strain M24A2M^T consisted mainly of diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol, phosphatidylcholine, one unidentified phospholipid, and two unidentified lipids. The respiratory quinone was ubiquinone Q-10. The genomic DNA G+C content was 63.9 %. Given the phenotypic characteristics along with the phylogenetic distinctness and chemotaxonomic features, strain M24A2M^T is considered to represent a novel species within the genus Novosphingobium, for which the name Novosphingobium aquimarinum sp. nov. is proposed. The type strain of Novosphingobium aquimarinum sp. nov. is M24A2M^T (=KCTC 72894^T=JCM 33983^T).

Paucihalobacter ruber gen. nov., sp. nov., isolated from a haloalkaline lake sediment sample

A Gram-stain-negative, strictly aerobic, non-motile, rod-shaped bacterium, designated CWB-1^T, was isolated from a haloalkaline lake sediment sample collected from the bottom of Chaiwopu Lake, Urumchi, Xinjiang Province, PR China. Strain CWB-1^T grew at 4-40 °C (optimum, 30-35 °C), pH 6.5-9.0 (optimum, pH 6.5-7.0) and with 0.5-5.5 % (w/v) NaCl (optimum, 2.5-3.0 %). Phylogenetic analyses based on the 16S rRNA gene sequence and the whole genome sequence both revealed that strain CWB-1^T belonged to the family Flavobacteriaceae. The strain had the highest similarity of the 16S rRNA gene sequence to Psychroserpens jangbogonensis PAMC 27130^T (92.8 %). The genome of strain CWB-1^T was 3 548 011 bp long with 36.3 % DNA G+C content. The predominant fatty acids (>10 %) in the CWB-1^T cells were iso-C_15 : 0, iso-C_17 : 0 3-OH and summed feature 1 (iso-C_15 : 1 H/C_13 : 0 3-OH). The major respiratory quinone was menaquinone-6 and the major polar lipids were phosphatidylethanolamine, an unidentified aminolipid and two unidentified lipids. Based on the phylogenetic analyses, as well as the phenotypic characteristics, a novel genus and species of the family Flavobacteriaceae, Paucihalobacter ruber gen. nov., sp. nov., is proposed. The type strain is CWB-1^T (=KCTC 72450^T=CGMCC 1.17149^T).

Genome compartmentalization predates species divergence in the plant pathogen genus Zymoseptoria

Background

Antagonistic co-evolution can drive rapid adaptation in pathogens and shape genome architecture. Comparative genome analyses of several fungal pathogens revealed highly variable genomes, for many species characterized by specific repeat-rich genome compartments with exceptionally high sequence variability. Dynamic genome structure may enable fast adaptation to host genetics. The wheat pathogen Zymoseptoria tritici with its highly variable genome, has emerged as a model organism to study genome evolution of plant pathogens. Here, we compared genomes of Z. tritici isolates and of sister species infecting wild grasses to address the evolution of genome composition and structure.

Results

Using long-read technology, we sequenced and assembled genomes of Z. ardabiliae, Z. brevis, Z. pseudotritici and Z. passerinii, together with two isolates of Z. tritici. We report a high extent of genome collinearity among Zymoseptoria species and high conservation of genomic, transcriptomic and epigenomic signatures of compartmentalization. We identify high gene content variability both within and between species. In addition, such variability is mainly limited to the accessory chromosomes and accessory compartments. Despite strong host specificity and non-overlapping host-range between species, predicted effectors are mainly shared among Zymoseptoria species, yet exhibiting a high level of presence-absence polymorphism within Z. tritici. Using in planta transcriptomic data from Z. tritici, we suggest different roles for the shared orthologs and for the accessory genes during infection of their hosts.

Conclusion

Despite previous reports of high genomic plasticity in Z. tritici, we describe here a high level of conservation in genomic, epigenomic and transcriptomic composition and structure across the genus Zymoseptoria. The compartmentalized genome allows the maintenance of a functional core genome co-occurring with a highly variable accessory genome.

Halomicroarcula amylolytica sp. nov., a novel halophilic archaeon isolated from a salt mine

A halophilic archaeon named strain LR21^T was isolated from a salt mine in Yunnan Province, PR China. Cells were spherical, Gram-stain-negative and motile. Strain LR21^T grew at 20-50 °C (optimum, 42 °C), with 8-30 % (w/v) NaCl (optimum, 23 %) and at pH 5.5-9.5 (optimum, pH 7.5-8.5). Mg²⁺ was not required for growth. The major polar lipid profile comprised phosphatidylglycerol, phosphatidylglycerol phosphate methyl ester and phosphatidylglycerol sulfate. Strain LR21^T had two dissimilar 16S rRNA genes (rrnA and rrnB) and they were closely related to Halomicroarcula limicola YGHS32^T, Hma. pellucida BNERC31^T and Hma. salina YGHS18^T with sequence similarities of 95.3-99.0, 93.0-96.2 and 93.2-95.9 %, respectively, and much lower values to other members. The rpoB' gene sequence similarities between strain LR21^T and Hma. limicola YGHS32^T, Hma. pellucida BNERC31^T and Hma. salina YGHS18^T were 95.2, 91.2 and 91.2 % respectively. The values of average nucleotide identity (ANI) and average amino-acid identity (AAI) between strain LR21^T and Hma. limicola YGHS32^T, were 89.0 and 90.1 %, respectively. DNA relatedness between strains LR21^T and Hma. limicola YGHS32^T determined by in silico DNA-DNA hybridization was 36.8 %. Values of ANI and AAI between strain LR21^T and other members in the genus Halomicroarcula were far below 95 % and the DNA-DNA relatedness values between strain LR21^T and its close relatives were much lower than 70 %, which is far below the boundary for delineation of a new species prokaryote. The DNA G+C content of strain LR21^T was 62.0 mol% (genome). The results suggested that strain LR21^T represents a novel species of the genus Halomicroarcula, for which the name Halomicroarcula amylolytica sp. nov. is proposed. The type strain is LR21^T (=CGMCC 1.13611^T=NBRC 113588^T).

Chitinophaga solisilvae sp. nov., isolated from forest soil

A Gram-stain-negative, strictly aerobic, non-motile and rod-shaped bacterial strain, O9T, was isolated from a forest soil sample collected at Dai, Xishuangbanna, Yunnan Province, PR China. Strain O9T grew optimally at pH 7.0, at 28‒30 °C and in the absence of NaCl. 16S rRNA gene sequence analysis placed strain O9T within the genus Chitinophaga of the family Chitinophagaceae, with Chitinophaga terrae KP01T (97.8 %), Chitinophaga jiangningensis JN53T (97.7 %), Chitinophaga niastensis JS16-4T (97.4 %), Chitinophaga qingshengii JN246T (97.3 %) and Chitinophaga dinghuensis DHOC24T (97.3 %) as its closest relatives. Strain O9T hydrolysed casein, gelatin and Tween 80. It could not assimilate l-arabinose, l-rhamnose, sucrose, melibiose, gentiobiose or d-fructose as a carbon source. It was negative for esterase lipase (C8) and β-glucosidase. Phosphatidylethanolamine was the predominant polar lipid. The major respiratory quinone of strain O9T was MK-7. Its major fatty acids were iso-C15:0 (34.2 %), C16:1 ω5c (20.9 %) and iso-C17:0 3-OH (12.6 %). The genomic DNA G+C content of strain O9T was 49.0 mol% based on total genome calculations. The average nucleotide identity score between the genomic sequence of strain O9T and that of Chitinophaga terrae KP01T was 72.9%. The Genome-to-Genome Distance Calculator showed that DNA‒DNA hybridization values for strain O9T and Chitinophaga terrae KP01T were 13.6, 21.1 and 14.4%, respectively. Based on the polyphasic taxonomic data, strain O9T represents a novel species of the genus Chitinophaga , for which the name Chitinophaga solisilvae sp. nov. is proposed. The type strain is O9T (=CGMCC 1.12462T=KCTC 32404T).

Whole-Genome Comparative and Pathogenicity Analysis of Salmonella enterica subsp. enterica Serovar Rissen

Salmonella are a type of bacteria known to cause food-borne illness. Their host range varies widely, and their susceptibility to the host determines its pathogenicity. Salmonella enterica serovar Rissen (S. Rissen) is a widely distributed serotype; however, its virulence and pathogenicity are poorly understood. In this study, the pathogenicity and antibiotic resistance of a representative S. Rissen isolate were investigated. The cell model results showed that S. Rissen preferred to replicate in human macrophage cells U937 compared to murine macrophage cells RAW264.7, suggesting that it has a level of host adaptability. Genome sequencing and comparison analysis revealed that the distribution and nonsynonymous single nucleotide polymorphisms of virulence factors in S. Rissen were similar to those in S. Typhi rather than to those in S. Typhimurium. Taken together, our results suggest that although S. Rissen is a common serotype distributed in swine herds, pork and chicken products, it has strong ability to infect humans.

Genome-Wide Analysis Reveals Genetic Potential for Aromatic Compounds Biodegradation of Sphingopyxis

Members of genus Sphingopyxis are frequently found in diverse eco-environments worldwide and have been traditionally considered to play vital roles in the degradation of aromatic compounds. Over recent decades, many aromatic-degrading Sphingopyxis strains have been isolated and recorded, but little is known about their genetic nature related to aromatic compounds biodegradation. In this study, bacterial genomes of 19 Sphingopyxis strains were used for comparative analyses. Phylogeny showed an ambiguous relatedness between bacterial strains and their habitat specificity, while clustering based on Cluster of Orthologous Groups suggested the potential link of functional profile with substrate-specific traits. Pan-genome analysis revealed that 19 individuals were predicted to share 1,066 orthologous genes, indicating a high genetic homogeneity among Sphingopyxis strains. Notably, KEGG Automatic Annotation Server results suggested that most genes pertaining aromatic compounds biodegradation were predicted to be involved in benzoate, phenylalanine, and aminobenzoate metabolism. Among them, β-ketoadipate biodegradation might be the main pathway in Sphingopyxis strains. Further inspection showed that a number of mobile genetic elements varied in Sphingopyxis genomes, and plasmid-mediated gene transfer coupled with prophage- and transposon-mediated rearrangements might play prominent roles in the evolution of bacterial genomes. Collectively, our findings presented that Sphingopyxis isolates might be the promising candidates for biodegradation of aromatic compounds in pollution sites.

Comparative Genomics of Acinetobacter baumannii Clinical Strains From Brazil Reveals Polyclonal Dissemination and Selective Exchange of Mobile Genetic Elements Associated With Resistance Genes

Acinetobacter baumannii is an opportunistic bacterial pathogen infecting immunocompromised patients and has gained attention worldwide due to its increased antimicrobial resistance. Here, we report a comparative whole-genome sequencing and analysis coupled with an assessment of antibiotic resistance of 46 Acinetobacter strains (45 A. baumannii plus one Acinetobacter nosocomialis) originated from five hospitals from the city of Recife, Brazil, between 2010 and 2014. An average of 3,809 genes were identified per genome, although only 2,006 genes were single copy orthologs or core genes conserved across all sequenced strains, with an average of 42 new genes found per strain. We evaluated genetic distance through a phylogenetic analysis and MLST as well as the presence of antibiotic resistance genes, virulence markers and mobile genetic elements (MGE). The phylogenetic analysis recovered distinct monophyletic A. baumannii groups corresponding to five known (ST1, ST15, ST25, ST79, and ST113) and one novel ST (ST881, related to ST1). A large number of ST specific genes were found, with the ST79 strains having the largest number of genes in common that were missing from the other STs. Multiple genes associated with resistance to β-lactams, aminoglycosides and other antibiotics were found. Some of those were clearly mapped to defined MGEs and an analysis of those revealed known elements as well as a novel Tn7-Tn3 transposon with a clear ST specific distribution. An association of selected resistance/virulence markers with specific STs was indeed observed, as well as the recent spread of the OXA-253 carbapenemase encoding gene. Virulence genes associated with the synthesis of the capsular antigens were noticeably more variable in the ST113 and ST79 strains. Indeed, several resistance and virulence genes were common to the ST79 and ST113 strains only, despite a greater genetic distance between them, suggesting common means of genetic exchange. Our comparative analysis reveals the spread of multiple STs and the genomic plasticity of A. baumannii from different hospitals in a single metropolitan area. It also highlights differences in the spread of resistance markers and other MGEs between the investigated STs, impacting on the monitoring and treatment of Acinetobacter in the ongoing and future outbreaks.

Flavobacterium alkalisoli sp. nov., isolated from rhizosphere soil of Suaeda salsa

A Gram-negative, strictly aerobic, gliding motility, none-spore forming, yellow, rods bacterial strain, designated XS-5T, was isolated from rhizosphere soil of Suaeda salsa, in Tumd Right Banner, Inner Mongolia, PR China. A phylogenetic tree based on the 16S rRNA gene sequences and the phylogenomic tree both showed that strain XS-5T clustered with Flavobacterium beibuense F44-8T (shared 97.2 % of 16S rRNA gene similarity) and Flavobacterium rakeshii FCS-5T (97.6 %), and shared <96.0 % of 16S rRNA gene similarities with all other type strains. Strain XS-5T contained MK-6 as the major respiratory quinone. Its major polar lipids were phosphatidylethanolamine, an unidentified aminolipid and an unidentified lipid; and the major fatty acids were iso-C15 : 0, iso-C17 : 0 3-OH, C16 : 0, iso-C15 : 0 3-OH, Summed feature 3 (iso-C15 : 0 2-OH and/or C16 : 1  ω7c), and Summed feature 9 (iso-C17 : 1  ω9c and/or C16 : 0 10-methyl). The genome consisted of a 3 985 855 bp circular chromosome, with a G+C content of 37.9 mol%, predicting 3616 coding sequences genes, 45 tRNA genes and three rRNA operons. The average nucleotide identity, amino acid identity and digital DNA–DNA hybridization values of strain XS-5T to F. beibuense F44-8T and F. rakeshii FCS-5T were 79.2 and 79.2 %, 81.7 and 81.6 %, 22.3 and 22.2 %, respectively. The results of phylogenetic, physiological and biochemical tests allowed the discrimination of strain XS-5T from its phylogenetic relatives. Flavobacterium alkalisoli sp. nov. is therefore proposed with strain XS-5T (=CGMCC 1.17077T=KCTC 72459T) as the type strain.

Novel cultivated endophytic Verrucomicrobia reveal deep-rooting traits of bacteria to associate with plants

Despite the relevance of complex root microbial communities for plant health, growth and productivity, the molecular basis of these plant-microbe interactions is not well understood. Verrucomicrobia are cosmopolitans in the rhizosphere, nevertheless their adaptations and functions are enigmatic since the proportion of cultured members is low. Here we report four cultivated Verrucomicrobia isolated from rice, putatively representing four novel species, and a novel subdivision. The aerobic strains were isolated from roots or rhizomes of Oryza sativa and O. longistaminata. Two of them are the first cultivated endophytes of Verrucomicrobia, as validated by confocal laser scanning microscopy inside rice roots after re-infection under sterile conditions. This extended known verrucomicrobial niche spaces. Two strains were promoting root growth of rice. Discovery of root compartment-specific Verrucomicrobia permitted an across-phylum comparison of the genomic conformance to life in soil, rhizoplane or inside roots. Genome-wide protein domain comparison with niche-specific reference bacteria from distant phyla revealed signature protein domains which differentiated lifestyles in these microhabitats. Our study enabled us to shed light into the dark microbial matter of root Verrucomicrobia, to define genetic drivers for niche adaptation of bacteria to plant roots, and provides cultured strains for revealing causal relationships in plant-microbe interactions by reductionist approaches.

Codon Pairs are Phylogenetically Conserved: A comprehensive analysis of codon pairing conservation across the Tree of Life

Identical codon pairing and co-tRNA codon pairing increase translational efficiency within genes when two codons that encode the same amino acid are translated by the same tRNA before it diffuses from the ribosome. We examine the phylogenetic signal in both identical and co-tRNA codon pairing across 23 428 species using alignment-free and parsimony methods. We determined that conserved codon pairing typically has a smaller window size than the length of a ribosome, and codon pairing tracks phylogenies across various taxonomic groups. We report a comprehensive analysis of codon pairing, including the extent to which each codon pairs. Our parsimony method generally recovers phylogenies that are more congruent with the established phylogenies than our alignment-free method. However, four of the ten taxonomic groups did not have sufficient orthologous codon pairings and were therefore analyzed using only the alignment-free methods. Since the recovered phylogenies using only codon pairing largely match phylogenies from the Open Tree of Life and the NCBI taxonomy, and are comparable to trees recovered by other algorithms, we propose that codon pairing biases are phylogenetically conserved and should be considered in conjunction with other phylogenomic techniques.

Catellatospora sichuanensis sp. nov., a novel actinobacterium isolated from soil

A novel actinobacterium, designated strain H14505^T, was isolated from a soil sample collected in Hong Yuan, Sichuan, southwest PR China. The temperature, pH and NaCl ranges for growth were determined to be 15-35 °C (optimum, 28 °C), 6.0-8.0 (optimum, pH 7.0) and 0-2 % (w/v; optimum without NaCl), respectively. The polar lipdis detected for strain H14505^T were diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol, glycolipid and four unidentified lipids. The predominant menaquinones of strain H14505^T were MK-9(H₄) and MK-9(H₆), and the prevalent fatty acids (>10 %) were C_18 : 1 ω9c, C_17 : 1 ω8c, summed feature 5 (anteiso-C_18 : 0/ C_18 : 2 ω6,9c) and C_16 : 0. Phylogenetic analysis based on 16S rRNA gene and whole-genome sequences indicated that strain H14505^T showed high similarity to Catellatospora vulcania NEAU-JM1^T (99.0 %) and Catellatospora paridis NEAU-CL2^T (99.0 %), and formed a monophyletic clade within the the genus Catellatospora in the phylogenetic trees. However, the average nucleotide indentity and DNA-DNA hybridization values between strain H14505^T and closely related Catellatospora species showed that it belonged to a distinct species. Furthermore, the results of morphological, physiological and biochemical tests allowed further phenotypic differentiation of strain H14505^T from its closest relatives. Thus, it is proposed that strain H14505^T represents a novel species of the genus Catellatospora, for which the name Catellatospora sichuanensis sp. nov. is proposed. The type strain of Catellatospora sichuanensis is H14505^T (=JCM 32394^T=CICC 11042^T).

Chitinimonas arctica sp. nov., isolated from Arctic tundra soil

A Gram-stain-negative, rod-shaped, green-pigmented, aerobic and motile bacterium, strain R3-44^T, was isolated from Arctic tundra soil. Stain R3-44^T clustered closely with members of the genus Chitinimonas, which belongs to the family Burkholderiaceae, and showed the highest 16S rRNA sequence similarity to Chitinimonas naiadis AR2^T (96.10%). Strain R3-44^T grew optimally at pH 7.0, 28 °C and in the presence of 0-0.5 % (w/v) NaCl. The predominant respiratory isoprenoid quinone of strain R3-44^T was identified as ubiquinone Q-8. The polar lipids consisted of phosphatidylglycerol, phosphatidylethanolamine, unidentified aminolipid and unidentified phospholipid. The main fatty acids were summed feature 3 (comprising C_16 : 1  ω7c and/or C_16 : 1  ω6c, 40.6 %) and C_16 : 0 (29.3 %). The DNA G+C content of strain R3-44^T was 60.8 mol%. On the basis of the evidence presented in this study, strain R3-44^T represents a novel species of the genus Chitinimonas, for which the name Chitinimonas arctica sp. nov. is proposed, with the type strain R3-44^T (=CCTCC AB 2010422^T=KCTC 72602^T).

Decolorization and biodegradation of textile di-azo dye Acid Blue 113 by Pseudomonas stutzeri AK6

Textile industry is one of the anthropogenic activities that consume a large amount of water and pollute water bodies. It uses a massive amount of dyes, which is one of the main constituents of polluting textile effluent. In the present research, biodegradation of Acid Blue 113 dye, a commonly used textile di-azo dye, has been studied exploiting Pseudomonas stutzeri, strain AK6. The dye (300 ppm) was decolorized up to 86.2% within 96 h. The metabolites of Acid Blue 113 obtained after biodegradation were identified by various analytical techniques viz. HPLC (high-performance liquid chromatography) and GC-MS (gas chromatography-mass spectrometry). Genome analysis of isolate AK6 using IMG/M (Integrated Microbial Genomes and Microbiomes) system supported the role of azoreductase and laccase for the decolorization and degradation of azo dye. The ability of P. stutzeri AK6 to tolerate high amount of dye makes it a potential candidate for bioremediation and pre-processing to remove dyes from textile effluents.

Alignment- and reference-free phylogenomics with colored de Bruijn graphs

BACKGROUND

The increasing amount of available genome sequence data enables large-scale comparative studies. A common task is the inference of phylogenies-a challenging task if close reference sequences are not available, genome sequences are incompletely assembled, or the high number of genomes precludes multiple sequence alignment in reasonable time.

RESULTS

We present a new whole-genome based approach to infer phylogenies that is alignment- and reference-free. In contrast to other methods, it does not rely on pairwise comparisons to determine distances to infer edges in a tree. Instead, a colored de Bruijn graph is constructed, and information on common subsequences is extracted to infer phylogenetic splits.

CONCLUSIONS

The introduced new methodology for large-scale phylogenomics shows high potential. Application to different datasets confirms robustness of the approach. A comparison to other state-of-the-art whole-genome based methods indicates comparable or higher accuracy and efficiency.

Functional role of bloom-forming cyanobacterium Planktothrix in ecologically shaping aquatic environments

Diverse metabolic behaviors endow microorganisms with various ecological functions, and metabolic activities of microbial species may affect the environmental conditions of their habitats. In this study, genome-guided analysis of Planktothrix spp. first divided these strains into six distinct groups, and comparisons of Planktothrix genomes revealed the inter- and intra-species variation. Prediction of central metabolism showed the functional diversity with regard to uptake of carbon, nitrogen, and sulfur sources. As the carbon-fixing microorganisms, Planktothrix isolates played a critical role in transforming the atmospheric carbon into organic carbon-the waterbodies' pool of available carbon. Diazotrophic lifestyle in certain Planktothrix strains may provide valuable avenues for supporting the equilibrium community. Furthermore, genome mining supported the exploration of biosynthetic gene clusters dedicated to cyanobacterial natural products, mainly including non-ribosomal peptide, polyketide, cyanobactin, and microviridin. Notably, some Planktothrix strains had the potential to non-ribosomally synthesize the microcystin (MC), a potent cyclic heptapeptide toxin, and MC-mediated cycling might strengthen the association between MC-producing and MC-degrading microorganisms. In short, genome-wide study of Planktothrix strains advances our current understanding of their metabolic potential and especially ecological roles in shaping natural environments.

Insights into ecological roles and potential evolution of Mlr-dependent microcystin-degrading bacteria

Over decades many studies have focused on the biodegradation of microcystins (MCs), and some Mlr-dependent MC-degrading bacteria were recorded, but the ecological functions, metabolic traits, and potential evolution of these organisms remain poorly understood. In this study, 16S rRNA-based phylogeny unraveled a wide range of genetic diversity across bacterial lineage, accompanied by re-evaluation of taxonomic placement of some MC-degrading species. Genome-wide comparison showed that considerable genes unique in individual organisms were identified, suggesting genetic differentiation among these Mlr-dependent MC-degrading bacteria. Notably, analyses of metabolic profiles first revealed the presence of functional genes involved in phenylacetate biodegradation in the specialized genomic regions, and mlr gene cluster was located around the neighborhood. The identification of transposable elements further indicated that these genomic regions might undergo horizontal gene transfer events to recruit novel functionalities, suggesting an adaptive force driving genome evolution of these organisms. In short, phylogenetic and genetic content analyses of Mlr-dependent MC-degraders shed light on their metabolic potential, ecological roles, and bacterial evolution, and expand the understanding of ecological status of MCs biodegradation.

Lysobacter alkalisoli sp. nov., a chitin-degrading strain isolated from saline-alkaline soil

Strains of Lysobacter, thought to play vital roles in the environment for their high enzyme production capacity, are ubiquitous in various ecosystems. During an analysis of bacterial diversity in saline soil, a Gram-stain-negative, aerobic, chitin-degrading bacterial strain, designated SJ-36^T, was isolated from saline-alkaline soil sampled at Tumd Right Banner, Inner Mongolia, PR China. Strain SJ-36^T grew at 4-40 °C (optimum, 30 °C), pH 5.0-10.0 (optimum, pH 7.0-8.0) and 0-6 % NaCl (optimum, 1.0 %). Oxidase and catalase activities were positive. A phylogenetic tree based on 16S rRNA gene sequences and the phylogenomic tree both showed that strain SJ-36^T formed a tight clade with Lysobacter maris KMU-14^T (sharing 97.6 % 16S rRNA gene similarity) and Lysobacter aestuarii S2-C^T (97.8 %). The major polar lipids of strain SJ-36^T were phosphatidylethanolamine, diphosphatidylglycerol, phosphatidylglycerol, two unidentified lipids and one unidentified phospholipid. The major fatty acids were iso-C_15 : 0 (37.5 %), summed feature 9 (14.0 %; iso-C_17 : 1ω9c and/or C_16 : 0 10-methyl) and iso-C_11 : 0 (10.6 %). Q-8 was the predominant ubiquinone. Its genomic DNA G+C content was 66.6 mol%. The average nucleotide identity values of strain SJ-36^T to L. maris KMU-14^T, L. aestuarii S2-C^T and other type strains were 81.5, 79.1 and <79.0 %, respectively. The results of physiological, phenotypic and phylogenetic characterizations allowed the discrimination of strain SJ-36^T from its phylogenetic relatives. Lysobacter alkalisoli sp. nov. is therefore proposed with strain SJ-36^T (=CGMCC 1.16756^T=KCTC 43039^T) as the type strain.

Antarcticimicrobium sediminis gen. nov., sp. nov., isolated from Antarctic intertidal sediment, transfer of Ruegeria lutea to Antarcticimicrobium gen. nov. as Antarcticimicrobium luteum comb. nov

A Gram-stain-negative, aerobic, non-flagellated and rod- or ovoid-shaped bacterium, designated as strain S4J41^T, was isolated from Antarctic intertidal sediment. The isolate grew at 0-37 °C and with 0.5-10 % (w/v) NaCl. It reduced nitrate to nitrite and hydrolysed Tween 80 and gelatin. Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain S4J41^T constituted a distinct phylogenetic line within the family Rhodobacteraceae and was closely related with some species in the genera Ruegeria, Phaeobacter, Pseudopuniceibacterium, Sulfitobacter, Puniceibacterium and Poseidonocella with 98.6-95.7 % 16S rRNA gene sequence similarities. The major cellular fatty acids were C_16 : 0, summed feature 8 (C_18 : 1  ω7c and/or C_18 : 1  ω6c) and C_18 : 0 and the major polar lipids were phosphatidylglycerol, phosphatidylcholine, diphosphatidylglycerol, phosphatidylethanolamine and one unidentified aminolipid. The sole respiratory quinone was Q-10. The genomic DNA G+C content of strain S4J41^T was 60.3 mol%. Based on the phylogenetic, chemotaxonomic and phenotypic data obtained in this study, strain S4J41^T is considered to represent a novel species in a new genus within the family Rhodobacteraceae, for which the name Antarcticimicrobium sediminis gen. nov., sp. nov. is proposed. The type strain is S4J41^T (=MCCC 1K03508^T=KCTC 62793^T). Moreover, the transfer of Ruegeria lutea Kim et al. 2019 to Antarcticimicrobium gen. nov. as Antarcticimicrobium luteum comb. nov. (type strain 318-1^T=JCM 30927^T=KCTC 72105^T) is also proposed.

Roseovarius arcticus sp. nov., a bacterium isolated from Arctic marine sediment

An aerobic, Gram-stain-negative, motile, rod or long-rod-shaped bacterial isolate, strain MK6-18T, was isolated from a marine sediment sample from Kongsfjorden, Arctic. The bacterium grew optimally at 20 °C, pH 7.0 and in the presence of 1.0–2.0 % (w/v) NaCl. Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain MK6-18T belonged to the genus Roseovarius . Its closest phylogenetic neighbour was Roseovarius nanhaiticus NH52JT showing 96.97 % 16S rRNA gene sequence similarity. The genome of strain MK6-18T is 4.2 Mb long in size with a G+C content of 59.5 mol%. The average nucleotide identity value between the genomes of strain MK6-18T and Roseovarius nanhaiticus NH52JT, was 78.0 %. Similar to other species of the genus Roseovarius , strain MK6-18T had ubiquinone 10 as the predominant ubiquinone and C12 : 0, C16 : 0, summed feature 3 (C16 : 1ω7c/ω6c) and summed feature 8 (C18 : 1ω7c/ω6c) as the major fatty acids. The polar lipid pattern consisted of diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine and phosphatidylcholine; one unidentified polar lipid, one unidentified aminolipid and one unidentified lipid were also detected. This is the first time that a member of the genus Roseovarius has been isolated from the Arctic, which may promote the study of the distribution characteristics and environmental adaptability of this genus. On the basis of the data provided here, strain MK6-18T should be classed as representing a novel species of the genus Roseovarius , for which the name Roseovarius arcticus sp. nov. is proposed. The type strain is MK6-18T (=CCTCC AB 2018219T=KCTC 72187T).

Genomic Analysis of Bacillus megaterium NCT-2 Reveals Its Genetic Basis for the Bioremediation of Secondary Salinization Soil

Bacillus megaterium NCT-2 is a nitrate-uptake bacterial, which shows high bioremediation capacity in secondary salinization soil, including nitrate-reducing capacity, phosphate solubilization, and salinity adaptation. To gain insights into the bioremediation capacity at the genetic level, the complete genome sequence was obtained by using a multiplatform strategy involving HiSeq and PacBio sequencing. The NCT-2 genome consists of a circular chromosome of 5.19 Mbp and ten indigenous plasmids, totaling 5.88 Mbp with an average GC content of 37.87%. The chromosome encodes 5,606 genes, 142 tRNAs, and 53 rRNAs. Genes involved in the features of the bioremediation in secondary salinization soil and plant growth promotion were identified in the genome, such as nitrogen metabolism, phosphate uptake, the synthesis of organic acids and phosphatase for phosphate-solubilizing ability, and Trp-dependent IAA synthetic system. Furthermore, strain NCT-2 has great ability of adaption to environments due to the genes involved in cation transporters, osmotic stress, and oxidative stress. This study sheds light on understanding the molecular basis of using B. megaterium NCT-2 in bioremediation of the secondary salinization soils.

A comparative genomics study of 23 Aspergillus species from section Flavi

Section Flavi encompasses both harmful and beneficial Aspergillus species, such as Aspergillus oryzae, used in food fermentation and enzyme production, and Aspergillus flavus, food spoiler and mycotoxin producer. Here, we sequence 19 genomes spanning section Flavi and compare 31 fungal genomes including 23 Flavi species. We reassess their phylogenetic relationships and show that the closest relative of A. oryzae is not A. flavus, but A. minisclerotigenes or A. aflatoxiformans and identify high genome diversity, especially in sub-telomeric regions. We predict abundant CAZymes (598 per species) and prolific secondary metabolite gene clusters (73 per species) in section Flavi. However, the observed phenotypes (growth characteristics, polysaccharide degradation) do not necessarily correlate with inferences made from the predicted CAZyme content. Our work, including genomic analyses, phenotypic assays, and identification of secondary metabolites, highlights the genetic and metabolic diversity within section Flavi.

Phylogenetic Analysis of HIV-1 Genomes Based on the Position-Weighted K-mers Method

HIV-1 viruses, which are predominant in the family of HIV viruses, have strong pathogenicity and infectivity. They can evolve into many different variants in a very short time. In this study, we propose a new and effective alignment-free method for the phylogenetic analysis of HIV-1 viruses using complete genome sequences. Our method combines the position distribution information and the counts of the k-mers together. We also propose a metric to determine the optimal k value. We name our method the Position-Weighted k-mers (PWkmer) method. Validation and comparison with the Robinson-Foulds distance method and the modified bootstrap method on a benchmark dataset show that our method is reliable for the phylogenetic analysis of HIV-1 viruses. PWkmer can resolve within-group variations for different known subtypes of Group M of HIV-1 viruses. This method is simple and computationally fast for whole genome phylogenetic analysis.

An assessment of the genomics, comparative genomics and cellulose degradation potential of Mucilaginibacter polytrichastri strain RG4-7

In this study, whole genome sequencing and comparative genomic analyses were performed for Mucilaginibacter polytrichastri RG4-7 and its carboxymethyl cellulose degradation potential was assessed. The results showed that the genome of strain RG4-7 was 5.84 Mb and contained 5019 predicted genes, in which a high proportion of strain-specific genes were related to carbohydrate metabolism. The carboxymethyl cellulose (CMC) degradation and cellulase activity tests revealed the strong cellulose degradation ability, CMCase and β-glucosidase activity in strain RG4-7. Real-time RT-PCR testing of most cellulose degradation related glycoside hydrolase (GH) families showed that GH9 (OKS85969), GH1 (OKS85832), GH3 (OKS89331 and OKS85615) were significantly up-regulated when strain RG4-7 was inoculated with CMC-Na, which suggested that GH9, GH1 and GH3 might determine its cellulose degradation ability. Certainly, further research need to be done to elucidate cellulose degradation mechanisms in strain RG4-7 in order to develop its industrial application value in lignocellulosic biomass degradation and waste management.

Georhizobium profundi gen. nov., sp. nov., a piezotolerant bacterium isolated from a deep-sea sediment sample of the New Britain Trench

A novel alphaproteobacterium, strain WS11T, was isolated from a deep-sea sediment sample collected from the New Britain Trench. The full-length 16S rRNA gene of strain WS11T had the highest sequence similarity of 97.6 % to Rhizobium subbaraonis JC85T, followed by Mycoplana ramosa DSM 7292T (96.9 %) and Rhizobium azooxidifex Po 20/26T (96.8 %). Phylogenetic analysis of concatenated 16S rRNA, atpD and recA gene sequences showed that strain WS11T was deeply separated from the species within the family Rhizobiaceae. Phylogenomic analysis based on the whole-genome protein sequences showed that strain WS11T formed an independent monophyletic branch in the family Rhizobiaceae, paralleled with the species in the families Brucellaceae and Phyllobacteriaceae within the order Rhizobiales. Cells were Gram-stain-negative, oxidase- and catalase-positive, and aerobic short rods (1.5-2.4×0.9-1.0 µm). Growth was observed at salinities ranging from 0 to 5% (optimum, 1 %), from pH 6.5 to 9 (optimum, pH 7) and at temperatures between 20 and 30 °C (optimum, 28 °C). Strain WS11T was piezotolerant, growing optimally at 0.1 MPa (range 0.1-70 MPa). The main fatty acid was summed feature 8 (C18 : 1 ω7c/C18  : 1 ω 6c). The sole respiratory quinone was ubiquinone-10 (Q-10). The predominant polar lipids were phosphatidylcholine, two unidentified aminophospholipids and an unidentified phospholipid. The genome size was about 4.36 Mbp and the G+C content was 62.3 mol%. The combined genotypic and phenotypic data show that strain WS11T represents a novel species of a novel genus in the family Rhizobiaceae, for which the name Georhizobium profundi gen. nov., sp. nov. is proposed (type strain WS11T=MCCC 1K03498T=KCTC 62439T).

Dyadobacter flavalbus sp. nov., isolated from lake sediment

A novel bacterial strain, designated NS28^T, was isolated from interfacial sediment sampled at Taihu Lake, PR China. Cells were rod-shaped, Gram-negative, aerobic and non-motile on Reasoner's 2A medium. The results of phylogenetic analysis based on 16S rRNA gene sequences indicated that strain NS28^T was most closely related to species from the genus Dyadobacter, with 98.4 and 96.0 % 16S rRNA gene sequence similarity to its closest phylogenetic neighbours Dyadobacter sediminis CGMCC 1.12895^T and Dyadobacter luticola CCTCC AB 2017091^T, respectively. MK-7 was the only cellular menaquinone. The major fatty acids were summed feature 3 (C_16 : 1ω7c and/or C_16 : 1ω6c), iso-C_15 : 0 and C_16 : 1ω5c. The major polar lipids were phosphatidylethanolamine, one phospholipid, one aminolipid, one lipid and two unidentified lipids. Genomic analysis of strain NS28^T indicated that the total genome size was 6 477 094 bp with a G+C content of 44.8 mol%, 5380 protein-coding genes, 79 contigs and an N50 length of 299584 bp. On the basis of the genomic DNA sequence, the average nucleotide identity values were 90.5 and 74.1 % with D. sediminis CGMCC 1.12895^T and D. luticola CCTCC AB 2017091^T, respectively. Digital DNA-DNA hybridization results of strain NS28^T with D. sediminis CGMCC 1.12895^T and D. luticola CCTCC AB 2017091^T were 40.9 and 18.6 %, respectively. Based on the phenotypic, chemotaxonomic, phylogenetic and genome sequence data presented here, it is proposed that strain NS28^T represents a novel species of the genus Dyadobacter for which the name Dyadobacter flavalbus is proposed . The type strain is NS28^T (=NBRC 113854^T=MCCC 1K03764^T).

Physiological and genomic features of Paraoceanicella profunda gen. nov., sp. nov., a novel piezophile isolated from deep seawater of the Mariana Trench

A novel piezophilic alphaproteobacterium, strain D4M1^T, was isolated from deep seawater of the Mariana Trench. 16S rRNA gene analysis showed that strain D4M1^T was most closely related to Oceanicella actignis PRQ‐67^T (94.2%), Oceanibium sediminis O448^T (94.2%), and Thioclava electrotropha ElOx9^T (94.1%). Phylogenetic analyses based on both 16S rRNA gene and genome sequences showed that strain D4M1^T formed an independent monophyletic branch paralleled with the genus Oceanicella in the family Rhodobacteraceae. Cells were Gram‐stain‐negative, aerobic short rods, and grew optimally at 37°C, pH 6.5, and 3.0% (w/v) NaCl. Strain D4M1^T was piezophilic with the optimum pressure of 10 MPa. The principal fatty acids were C_18:1ω7c/C_18:1ω6c and C_16:0, major respiratory quinone was ubiquinone‐10, and predominant polar lipids were phosphatidylglycerol, phosphatidylethanolamine, and an unidentified aminophospholipid. The complete genome contained 5,468,583‐bp with a G + C content of 70.2 mol% and contained 4,855 protein‐coding genes and 78 RNA genes. Genomic analysis revealed abundant clues on bacterial high‐pressure adaptation and piezophilic lifestyle. The combined evidence shows that strain D4M1^T represents a novel species of a novel genus in the family Rhodobacteraceae, for which the name Paraoceanicella profunda gen. nov., sp. nov. is proposed (type strain D4M1^T = MCCC 1K03820^T = KCTC 72285^T).

Paraconexibacter algicola gen. nov., sp. nov., a novel actinobacterium isolated from a eutrophic lake during the end of cyanobacterial harmful algal blooms, and proposal of Paraconexibacteraceae fam. nov. in the order Solirubrobacterales

A novel bacterium, strain Seoho-28^T, was isolated from a shallow eutrophic lake during the end of cyanobacterial harmful algal blooms and was characterized taxonomically and phylogenetically. Strain Seoho-28^T was a Gram-stain-negative, aerobic, rod-shaped and non-motile bacterium. The strain grew optimally with 0 % NaCl and at 25-30 °C on Reasoner's 2A medium. The phylogenetic analysis based on 16S rRNA gene sequences positioned the novel strain among the order Solirubrobacterales, but sequence similarities to known species were less than 94.7 %. The genomic DNA G+C content of the strain Seoho-28^T was 74.2 mol%. Genomic comparisons of strain Seoho-28^T with families in the order Solirubrobacterales were made using the Genome-to-Genome Distance Calculator, average nucleotide identity and average amino acid identity analyses (values indicated ≤14.9, ≤73.5 and ≤57.8 %, respectively). Strain Seoho-28^T contained C_16 : 0-iso, C_18 : 1  ω9c and C_16 : 0 as major fatty acids and MK-7 (H₄) as the major quinone. Strain Seoho-28^T contained diphosphatidylglycerol, phosphatidylinositol and an unidentified phospholipid as major polar lipids. Meso- and ll-diaminopimelic acids were the diagnostic diamino acids in the cell-wall peptidoglycan. Based on the genotypic, chemotaxonomic and phenotypic results, strain Seoho-28^T represents a novel genus and species, Paraconexibacter algicola gen. nov., sp. nov., which belongs to a new family Paraconexibacteraceae in the order Solirubrobacterales and the class Thermoleophilia. The type strain is Seoho-28^T (=KCTC 39791^T=JCM 31881^T).