Defining bacterial species in the genomic era: insights from the genus Acinetobacter

Characterization of the diversity of type IV secretion system-encoding plasmids in Acinetobacter

The multi-drug resistant pathogen Acinetobacter baumannii has gained global attention as an important clinical challenge. Owing to its ability to survive on surfaces, its capacity for horizontal gene transfer, and its resistance to front-line antibiotics, A. baumannii has established itself as a successful pathogen. Bacterial conjugation is a central mechanism for pathogen evolution. The epidemic multidrug-resistant A. baumannii ACICU harbours a plasmid encoding a Type IV Secretion System (T4SS) with homology to the E. coli F-plasmid, and plasmids with homologous gene clusters have been identified in several A. baumannii sequence types. However the genetic and host strain diversity, global distribution, and functional ability of this group of plasmids is not fully understood. Using systematic analysis, we show that pACICU2 belongs to a group of almost 120 T4SS-encoding plasmids within four different species of Acinetobacter and one strain of Klebsiella pneumoniae from human and environmental origin, and globally distributed across 20 countries spanning 4 continents. Genetic diversity was observed both outside and within the T4SS-encoding cluster, and 47% of plasmids harboured resistance determinants, with two plasmids harbouring eleven. Conjugation studies with an extensively drug-resistant (XDR) strain showed that the XDR plasmid could be successfully transferred to a more divergent A. baumanii, and transconjugants exhibited the resistance phenotype of the plasmid. Collectively, this demonstrates that these T4SS-encoding plasmids are globally distributed and more widespread among Acinetobacter than previously thought, and that they represent an important potential reservoir for future clinical concern.

Embolic necrosuppurative pneumonia in domestic cats induced by a novel Neisseria species

Three cats, aged 2 to 11 years, presented to the University of Minnesota Veterinary Diagnostic Laboratory over a 3-year period following euthanasia or death due to respiratory distress. Thoracic radiographs revealed nodular, soft tissue opacities throughout the lung fields in all cases. On postmortem examination, approximately 60% to 80% of the lung parenchyma were expanded by multifocal to coalescing, well-demarcated, beige, semi-firm nodules. Histologically, large numbers of neutrophils, fewer macrophages, fibrin, and cellular and karyorrhectic debris effaced the pulmonary parenchyma. The inflammatory foci contained aggregates of gram-negative cocci. 16s rRNA Sanger sequencing and whole-genome sequencing identified the bacteria isolated from the lung of all cats under aerobic conditions as a novel Neisseria spp. Based on whole-genome sequence analysis, all 3 sequences shared 92.71% and 92.67% average nucleotide identity with closely related Neisseria animaloris NZ LR134440T and Neisseria animaloris GCA 002108605T, respectively. The in silico DNA-DNA hybridization identity compared to our isolates was 46.6% and 33.8% with strain DSM Neisseria zoodegmatis 21642 and strain DSM 21643, respectively. All 3 sequences have less than 95% average nucleotide identity and less than 70% DNA-DNA hybridization identity, suggesting that the 3 isolates are a novel species of the genus Neisseria. Infection with Neisseria spp. induces an embolic pneumonia in cats that radiographically and pathologically resembles a metastatic neoplastic process and should be considered among the etiologic differential diagnoses in cases of infectious pulmonary disease with a disseminated, nodular lung pattern.

First report on in-depth genome and comparative genome analysis of a metal-resistant bacterium Acinetobacter pittii S-30, isolated from environmental sample

A newly isolated bacterium Acinetobacter pittii S-30 was recovered from waste-contaminated soil in Ranchi, India. The isolated bacterium belongs to the ESKAPE organisms which represent the major nosocomial pathogens that exhibit high antibiotic resistance. Furthermore, average nucleotide identity (ANI) analysis also showed its closest match (>95%) to other A. pittii genomes. The isolate showed metal-resistant behavior and was able to survive up to 5 mM of ZnSO4. Whole genome sequencing and annotations revealed the occurrence of various genes involved in stress protection, motility, and metabolism of aromatic compounds. Moreover, genome annotation identified the gene clusters involved in secondary metabolite production (biosynthetic gene clusters) such as arylpolyene, acinetobactin like NRP-metallophore, betalactone, and hserlactone-NRPS cluster. The metabolic potential of A. pittii S-30 based on cluster of orthologous, and Kyoto Encyclopedia of Genes and Genomes indicated a high number of genes related to stress protection, metal resistance, and multiple drug-efflux systems etc., which is relatively rare in A. pittii strains. Additionally, the presence of various carbohydrate-active enzymes such as glycoside hydrolases (GHs), glycosyltransferases (GTs), and other genes associated with lignocellulose breakdown suggests that strain S-30 has strong biomass degradation potential. Furthermore, an analysis of genetic diversity and recombination in A. pittii strains was performed to understand the population expansion hypothesis of A. pittii strains. To our knowledge, this is the first report demonstrating the detailed genomic characterization of a heavy metal-resistant bacterium belonging to A. pittii. Therefore, the A. pittii S-30 could be a good candidate for the promotion of plant growth and other biotechnological applications.

Comparative Genomics Reveals Genetic Diversity and Variation in Metabolic Traits in Fructilactobacillus sanfranciscensis Strains

Fructilactobacillus sanfranciscensis is a significant and dominant bacterial species of sourdough microbiota from ecological and functional perspectives. Despite the remarkable prevalence of different strains of this species in sourdoughs worldwide, the drivers behind the genetic diversity of this species needed to be clarified. In this research, 14 F. sanfranciscensis strains were isolated from sourdough samples to evaluate the genetic diversity and variation in metabolic traits. These 14 and 31 other strains (obtained from the NCBI database) genomes were compared. The values for genome size and GC content, on average, turned out to 1.31 Mbp and 34.25%, respectively. In 45 F. sanfranciscensis strains, there were 162 core genes and 0 to 51 unique genes present in each strain. The primary functions of core genes were related to nucleotide, lipid transport, and amino acid, as well as carbohydrate metabolism. The size of core genes accounted for 41.18% of the pan-genome size in 14 F. sanfranciscensis strains, i.e., 0.70 Mbp of 1.70 Mbp. There were genetic variations among the 14 strains involved in carbohydrate utilization and antibiotic resistance. Moreover, exopolysaccharides biosynthesis-related genes were annotated, including epsABD, wxz, wzy. The Type IIA & IE CRISPR-Cas systems, pediocin PA-1 and Lacticin_3147_A1 bacteriocins operons were also discovered in F. sanfranciscensis. These findings can help to select desirable F. sanfranciscensis strains to develop standardized starter culture for sourdough fermentation, and expect to provide traditional fermented pasta with a higher quality and nutritional value for the consumers.

Taxonomic Diversity and Functional Traits of Soil Bacterial Communities under Radioactive Contamination: A Review

Studies investigating the taxonomic diversity and structure of soil bacteria in areas with enhanced radioactive backgrounds have been ongoing for three decades. An analysis of data published from 1996 to 2024 reveals changes in the taxonomic structure of radioactively contaminated soils compared to the reference, showing that these changes are not exclusively dependent on contamination rates or pollutant compositions. High levels of radioactive exposure from external irradiation and a high radionuclide content lead to a decrease in the alpha diversity of soil bacterial communities, both in laboratory settings and environmental conditions. The effects of low or moderate exposure are not consistently pronounced or unidirectional. Functional differences among taxonomic groups that dominate in contaminated soil indicate a variety of adaptation strategies. Bacteria identified as multiple-stress tolerant; exhibiting tolerance to metals and antibiotics; producing antioxidant enzymes, low-molecular antioxidants, and radioprotectors; participating in redox reactions; and possessing thermophilic characteristics play a significant role. Changes in the taxonomic and functional structure, resulting from increased soil radionuclide content, are influenced by the combined effects of ionizing radiation, the chemical toxicity of radionuclides and co-contaminants, as well as the physical and chemical properties of the soil and the initial bacterial community composition. Currently, the quantification of the differential contributions of these factors based on the existing published studies presents a challenge.

Two novel Bartonella (sub)species isolated from edible dormice (Glis glis): hints of cultivation stress-induced genomic changes

Bartonelloses are neglected emerging infectious diseases caused by facultatively intracellular bacteria transmitted between vertebrate hosts by various arthropod vectors. The highest diversity of Bartonella species has been identified in rodents. Within this study we focused on the edible dormouse (Glis glis), a rodent with unique life-history traits that often enters households and whose possible role in the epidemiology of Bartonella infections had been previously unknown. We identified and cultivated two distinct Bartonella sub(species) significantly diverging from previously described species, which were characterized using growth characteristics, biochemical tests, and various molecular techniques including also proteomics. Two novel (sub)species were described: Bartonella grahamii subsp. shimonis subsp. nov. and Bartonella gliris sp. nov. We sequenced two individual strains per each described (sub)species. During exploratory genomic analyses comparing two genotypes ultimately belonging to the same species, both factually and most importantly even spatiotemporally, we noticed unexpectedly significant structural variation between them. We found that most of the detected structural variants could be explained either by prophage excision or integration. Based on a detailed study of one such event, we argue that prophage deletion represents the most probable explanation of the observed phenomena. Moreover, in one strain of Bartonella grahamii subsp. shimonis subsp. nov. we identified a deletion related to Bartonella Adhesin A, a major pathogenicity factor that modulates bacteria-host interactions. Altogether, our results suggest that even a limited number of passages induced sufficient selective pressure to promote significant changes at the level of the genome.

Assessing the role of the gut microbiome in methylmercury demethylation and elimination in humans and gnotobiotic mice

The risk of methylmercury (MeHg) toxicity following ingestion of contaminated foodstuffs (e.g., fish) is directly related to the kinetics of MeHg elimination among individuals. Yet, the factors driving the wide range of inter-individual variability in MeHg elimination within a population are poorly understood. Here, we investigated the relationship between MeHg elimination, gut microbiome demethylation activity, and gut microbiome composition using a coordinated human clinical trial and gnotobiotic mouse modeling approach together with metagenomic sequence analysis. We first observed MeHg elimination half-lives (t1/2) ranging from 28 to 90 days across 27 volunteers. Subsequently, we found that ingestion of a prebiotic induced changes in the gut microbiome and mixed effects (increased, decrease, and no effect) on elimination in these same individuals. Nonetheless, elimination rates were found to correlate with MeHg demethylation activity in cultured stool samples. In mice, attempts to remove the microbiome via generation of germ-free (GF) animals or through antibiotic (Abx) treatment both diminished MeHg demethylation to a similar extent. While both conditions substantially slowed elimination, Abx treatment resulted in significantly slower elimination than the GF condition, indicating an additional role for host-derived factors in supporting elimination. Human fecal microbiomes transplanted to GF mice restored elimination rates to that seen in control mice. Metagenomic sequence analysis of human fecal DNA did not identify genes encoding proteins typically involved in demethylation (e.g., merB, organomercury lyase). However, the abundance of several anaerobic taxa, notably Alistipes onderdonkii, were positively correlated with MeHg elimination. Surprisingly, mono-colonization of GF free mice with A. onderdonkii did not restore MeHg elimination to control levels. Collectively, our findings indicate the human gut microbiome uses a non-conventional pathway of demethylation to increase MeHg elimination that relies on yet to be resolved functions encoded by the gut microbes and the hostClinical Trial NCT04060212, prospectively registered 10/1/2019.

Comparative genomic analysis of the Dietzia genus: an insight into genomic diversity, and adaptation

Dietzia strains are widely distributed in the environment, presenting an opportunistic role, and some species have undetermined taxonomic characteristics. Here, we propose the existence of errors in the classification of species in this genus using comparative genomics. We performed ANI, dDDH, pangenome and genomic plasticity analyses better to elucidate the phylogenomic relationships between Dietzia strains. For this, we used 55 genomes of Dietzia downloaded from public databases that were combined with a newly sequenced. Sequence analysis of a phylogenetic tree based on genome similarity comparisons and dDDH, ANI analyses supported grouping different Dietzia species into four distinct groups. The pangenome analysis corroborated the classification of these groups, supporting the idea that some species of Dietzia could be reassigned in a possible classification into three distinct species, each containing less variability than that found within the global pangenome of all strains. Additionally, analysis of genomic plasticity based on groups containing Dietzia strains found differences in the presence and absence of symbiotic Islands and pathogenic islands related to their isolation site. We propose that the comparison of pangenome subsets together with phylogenomic approaches can be used as an alternative for the classification and differentiation of new species of the genus Dietzia.

Phylogenomic analysis of Citrobacter sp. strain AAK_AS5 and its metabolic capabilities to support nitrogen removal behavior

Despite the ubiquity of the genus Citrobacter in clinical, industrial, and environmental scenarios, a large number of Citrobacter strains have not been explored at the genome-scale level. In this study, accurate taxonomic assignment of strain AAK_AS5 isolated from activated sludge was achieved by in-silico genomic comparison using Overall Genome-based Relatedness Indices (ANI(OAT): 97.55%, ANIb:97.28%, and ANIm: 97.83%) that indicated its closest identity to the related strain Citrobacter portucalensis A60^T . Results were consistent with a digital DNA-DNA hybridization value of 80% with C. portucalensis A60^T which was greater than the species boundary value >70% for delineating closely related bacterial species. Gene mining through Kyoto Encyclopedia of Genes and Genomes (KEGG), and annotation using rapid annotation subsystem technology (RAST) revealed the notable gene contents for nitrogen metabolism and other pathways associated with nitrate/nitrite ammonification (28 genes), ammonia assimilation (22 genes), and denitrification pathways (14 genes). Furthermore, the strain AAK_AS5 also exhibited a high soluble chemical oxygen demand (sCOD), NH₄ ⁺ -N, and NO₃ ^- -N removal efficiency of 91.4%, 90%, and 93.6%, respectively thus validating its genetic capability for utilizing both (NH₄ )₂ SO₄ and KNO₃ as the nitrogen source. The study provided deeper insights into the phylogenomics and the genetic potential of Citrobacter, sp. strain AAK AS5 associated with nitrogen metabolism thus signifying the potential application of the isolate for treating nitrogen-rich wastewaters.

Comparative genomics and informational content analysis uncovered internal regions of the core genes rpoD, pepN and gltX for an MLSA with genome-level resolving power within the genus Pseudomonas

In the field of prokaryotic taxonomy, there has been a recent transition towards phylogenomics as the gold standard approach. However, genome-based phylogenetics is still restrictive for its cost when managing large amounts of isolates. Fast, cheap, and taxonomically competent alternatives, like multilocus sequence analysis (MLSA) are thus recommendable. Nevertheless, the criteria for selecting the conserved genes for MLSA have not been explicit for different bacterial taxa, including the broadly diverse Pseudomonas genus. Here, we have carried out an unbiased and rational workflow to select internal sequence regions of Pseudomonas core genes (CG) for a MLSA with the best phylogenetic power, and with a resolution comparable to the genome-based ANI approach. A computational workflow was established to inspect 126 complete genomes of representatives from over 60 Pseudomonas species and subspecies, in order to identify the most informative CG internal regions and determine which combinations in sets of three partial CG sequences have comparable phylogenetic resolution to that of the current ANI standard. We found that the rpoD^346-1196-pepN^1711-2571-gltX^86-909 concatenated sequences were the best performing in terms of phylogenetic robustness and resulted highly sensitive and specific when contrasted with ANI. The rpoD-pepN-gltX MLSA was validated in silico and in vitro. Altogether, the results presented here supports the proposal of the rpoD-pepN-gltX MLSA as a fast, affordable, and robust phylogenetic tool for members of the Pseudomonas genus.

Halocatena salina sp. nov., a filamentous halophilic archaeon isolated from Aiding Salt Lake

A filamentous cell-shaped halophilic archaeon (strain AD-1^T) was isolated from Aiding Salt Lake, PR China. Its colonies on HCM7 agar plates were pinkish white, 1-4 mm (diameter), elevated and round. The optimum conditions for growth were observed at 42 °C, 4.3 M NaCl, 0.01 M MgCl₂ and pH 7. Strain AD-1^T could hydrolyse Tween 60, Tween 80, starch and gelatin. Phylogenetic analysis based on 16S rRNA gene, rpoB' and the concatenated 484 single-copy orthologous proteins revealed that strain AD-1^T formed a clade with Halocatena pleomorpha SPP-AMP-1^T. The average nucleotide identity and in silico DNA-DNA hybridization values between strain AD-1^T and Halocatena pleomorpha SPP-AMP-1^T were both below the species delineation thresholds (95~96 and 70 %, respectively). The major phospholipids of strain AD-1^T were phosphatidic acid, phosphatidylglycerol and phosphatidylglycerol phosphate methyl ester, while the major glycolipids were sulphated galactosyl mannosyl glucosyl diether, galactosyl mannosyl glucosyl diether and glucosyl mannosyl glucosyl diether. The phenotypic, phylogenetic and genome-based analyses suggested that strain AD-1^T (=CGMCC 1.13724^T=JCM 32960^T) represents a novel species, for which the name Halocatena salina sp. nov. is proposed.

Acinetobacter sedimenti sp. nov., isolated from beach sediment

A Gram-stain-negative, aerobic, non-motile, non-haemolytic, oxidase-negative, catalase-positive bacillus strain (A3.8^T) was isolated from beach sediment from Zhairuo Island, PR China. The strain grew at pH 6.0-9.0 (optimum, 7.0), with 0-4.5 % NaCl (optimum, 2 %) and at 10-35 °C (optimum, 30 °C). Its whole-genome sequence was 2.5 Mb in size, with a DNA G+C content of 41.6 mol%. On the basis of the results of core genome phylogenetic analysis, A3.8^T represents a separate branch within the clade formed by five species of the genus Acinetobacter with 'Acinetobacter marinus' as the most closely related species. The average nucleotide identity compared with the closely related species of the genus Acinetobacter was below 83.66 % and digital DNA-DNA hybridization values were less than 28.80 %. The predominant fatty acids included C_18 : 1ω9c, C_16 : 0 and summed feature 3 (C_16 : 1ω7c and/or C_16 : 1ω6c). Q-9 was the major respiratory quinone. The polar lipids are mainly composed of diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, two phospholipids, an aminolipid and four unknown lipids. A3.8^T cannot assimilate dl-lactate and weakly utilizes l-glutamate, l-leucine, l-phenylalanine and l-tartrate, which distinguishes it from other species of the genus Acinetobacter. On the basis of the genotype, phenotype and biochemical data, strain A3.8^T represents a novel species of the genus Acinetobacter, for which the name Acinetobacter sedimenti sp. nov. is proposed. The type strain is A3.8^T (=MCCC 1K07161^T=LMG 32568^T).

Phenotype-genotype analysis of Latilactobacills curvatus from different niches: Carbohydrate metabolism, antibiotic resistance, bacteriocin, phage fragments and linkages with CRISPR-Cas systems

The potential probiotic function of Latilactobacills curvatus has attracted the attention of researchers. To explore the differences in the genomes of L. curvatus, nine strains were isolated from various sources, including feces and fermented vegetables and compared with 25 strains from the NCBI database. The findings indicated that the average genome size, GC content, and CDS of L. curvatus were 1.94 MB, 41.9%, and 1825, respectively. Its core genome is associated with transcription, translation, carbohydrate transport and metabolism, and defense functions. The pan-genome of L. curvatus was in a closed state. The genetic diversity of L. curatus is mainly manifested in its ability to use carbohydrates, antibiotic resistance, bacteriocin operon, and polymeric regularly interspaced short palindromic repeats (CRISPR)-Cas for bacterial immunity. The CRISPR system of 34 strains of L. curvatus was predominantly found to be of the IIA type with a few IIC and IE types. These findings will contribute to a better understanding of this species.

Development of the SeqCode: A proposed nomenclatural code for uncultivated prokaryotes with DNA sequences as type

Over the last fifteen years, genomics has become fully integrated into prokaryotic systematics. The genomes of most type strains have been sequenced, genome sequence similarity is widely used for delineation of species, and phylogenomic methods are commonly used for classification of higher taxonomic ranks. Additionally, environmental genomics has revealed a vast diversity of as-yet-uncultivated taxa. In response to these developments, a new code of nomenclature, the Code of Nomenclature of Prokaryotes Described from Sequence Data (SeqCode), has been developed over the last two years to allow naming of Archaea and Bacteria using DNA sequences as the nomenclatural types. The SeqCode also allows naming of cultured organisms, including fastidious prokaryotes that cannot be deposited into culture collections. Several simplifications relative to the International Code of Nomenclature of Prokaryotes (ICNP) are implemented to make nomenclature more accessible, easier to apply and more readily communicated. By simplifying nomenclature with the goal of a unified classification, inclusive of both cultured and uncultured taxa, the SeqCode will facilitate the naming of taxa in every biome on Earth, encourage the isolation and characterization of as-yet-uncultivated taxa, and promote synergies between the ecological, environmental, physiological, biochemical, and molecular biological disciplines to more fully describe prokaryotes.

Characterization of symbiotic and nitrogen fixing bacteria

Symbiotic nitrogen fixing bacteria comprise of diverse species associated with the root nodules of leguminous plants. Using an appropriate taxonomic method to confirm the identity of superior and elite strains to fix nitrogen in legume crops can improve sustainable global food and nutrition security. The current review describes taxonomic methods preferred and commonly used to characterize symbiotic bacteria in the rhizosphere. Peer reviewed, published and unpublished articles on techniques used for detection, classification and identification of symbiotic bacteria were evaluated by exploring their advantages and limitations. The findings showed that phenotypic and cultural techniques are still affordable and remain the primary basis of species classification despite their challenges. Development of new, robust and informative taxonomic techniques has really improved characterization and identification of symbiotic bacteria and discovery of novel and new species that are effective in biological nitrogen fixation (BNF) in diverse conditions and environments.

Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies

BACKGROUND

The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82].

RESULTS

At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%).

CONCLUSIONS

Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.

New Insights into the Taxonomy of Bacteria in the Genomic Era and a Case Study with Rhizobia

Since early studies, the history of prokaryotes taxonomy has dealt with many changes driven by the development of new and more robust technologies. As a result, the number of new taxa descriptions is exponentially increasing, while an increasing number of others has been subject of reclassification, demanding from the taxonomists more effort to maintain an organized hierarchical system. However, expectations are that the taxonomy of prokaryotes will acquire a more stable status with the genomic era. Other analyses may continue to be necessary to determine microbial features, but the use of genomic data might be sufficient to provide reliable taxa delineation, helping taxonomy to reach the goal of correct classification and identification. Here we describe the evolution of prokaryotes' taxonomy until the genomic era, emphasizing bacteria and taking as an example the history of rhizobia taxonomy. This example was chosen because of the importance of the symbiotic nitrogen fixation of legumes with rhizobia to the nitrogen input to both natural ecosystems and agricultural crops. This case study reports the technological advances and the methodologies used to classify and identify bacterial species and indicates the actual rules required for an accurate description of new taxa.

Whole-Genome Shotgun Metagenomic Sequencing Reveals Distinct Gut Microbiome Signatures of Obese Cats

Overweight and obesity are growing health problems in domestic cats, increasing the risks of insulin resistance, lipid dyscrasias, neoplasia, cardiovascular disease, and decreasing longevity. The signature of obesity in the feline gut microbiota has not been studied at the whole-genome metagenomic level. We performed whole-genome shotgun metagenomic sequencing in the fecal samples of eight overweight/obese and eight normal cats housed in the same research environment. We obtained 271 Gbp of sequences and generated a 961-Mbp de novo reference contig assembly, with 1.14 million annotated microbial genes. In the obese cat microbiome, we discovered a significant reduction in microbial diversity (P < 0.01) and Firmicutes abundance (P = 0.005), as well as decreased Firmicutes/Bacteroidetes ratios (P = 0.02), which is the inverse of obese human/mouse microbiota. Linear discriminant analysis and quantitative PCR (qPCR) validation revealed significant increases of Bifidobacterium sp., Olsenella provencensis, Dialister sp.CAG:486, and Campylobacter upsaliensis as the hallmark of obese microbiota among 400 enriched species, whereas 1,525 bacterial species have decreased abundance in the obese microbiome. Phascolarctobacterium succinatutens and an uncharacterized Erysipelotrichaceae bacterium are highly abundant (>0.05%) in the normal gut with over 400-fold depletion in the obese microbiome. Fatty acid synthesis-related pathways are significantly overrepresented in the obese compared with the normal cat microbiome. In conclusion, we discovered dramatically decreased microbial diversity in obese cat gut microbiota, suggesting potential dysbiosis. A panel of seven significantly altered, highly abundant species can serve as a microbiome indicator of obesity. Our findings in the obese cat microbiome composition, abundance, and functional capacities provide new insights into feline obesity.

IMPORTANCE Obesity affects around 45% of domestic cats, and licensed drugs for treating feline obesity are lacking. Physical exercise and calorie restrictions are commonly used for weight loss but with limited efficacy. Through comprehensive analyses of normal and obese cat gut bacteria flora, we identified dramatic shifts in the obese gut microbiome, including four bacterial species significantly enriched and two species depleted in the obese cats. The key bacterial community and functional capacity alterations discovered from this study will inform new weight management strategies for obese cats, such as evaluations of specific diet formulas that alter the microbiome composition, and the development of prebiotics and probiotics that promote the increase of beneficial species and the depletion of obesity-associated species. Interestingly, these bacteria identified in our study were also reported to affect the weight loss success in human patients, suggesting translational potential in human obesity.

In-depth genome analysis of Bacillus sp. BH32, a salt stress-tolerant endophyte obtained from a halophyte in a semiarid region

Endophytic strains belonging to the Bacillus cereus group were isolated from the halophytes Atriplex halimus L. (Amaranthaceae) and Tamarix aphylla L. (Tamaricaceae) from costal and continental regions in Algeria. Based on their salt tolerance (up to 5%), the strains were tested for their ability to alleviate salt stress in tomato and wheat. Bacillus sp. strain BH32 showed the highest potential to reduce salinity stress (up to + 50% and + 58% of dry weight improvement, in tomato and wheat, respectively, compared to the control). To determine putative mechanisms involved in salt tolerance and plant growth promotion, the whole genome of Bacillus sp. BH32 was sequenced, annotated, and used for comparative genomics against the genomes of closely related strains. The pangenome of Bacillus sp. BH32 and its closest relative was further analyzed. The phylogenomic analyses confirmed its taxonomic position, a member of the Bacillus cereus group, with intergenomic distances (GBDP analysis) pinpointing to a new taxon (digital DNA-DNA hybridization, dDDH < 70%). Genome mining unveiled several genes involved in stress tolerance, production of anti-oxidants and genes involved in plant growth promotion as well as in the production of secondary metabolites. KEY POINTS : • Bacillus sp. BH32 and other bacterial endophytes were isolated from halophytes, to be tested on tomato and wheat and to limit salt stress adverse effects. • The strain with the highest potential was then studied at the genomic level to highlight numerous genes linked to plant growth promotion and stress tolerance. • Pangenome approaches suggest that the strain belongs to a new taxon within the Bacillus cereus group.

Comparative Genomic Analysis Reveals Intestinal Habitat Adaptation of Ligilactobacillus equi Rich in Prophage and Degrading Cellulase

Ligilactobacillus equi is common in the horse intestine, alleviates the infection of Salmonella, and regulates intestinal flora. Despite this, there have been no genomic studies on this species. Here, we provide the genomic basis for adaptation to the intestinal habitat of this species. We sequenced the genome of L. equi IMAU81196, compared this with published genome information from three strains in NCBI, and analyzed genome characteristics, phylogenetic relationships, and functional genes. The mean genome size of L. equi strains was 2.08 ± 0.09 Mbp, and the mean GC content was 39.17% ± 0.19%. The genome size of L. equi IMAU81196 was 1.95 Mbp, and the GC content was 39.48%. The phylogenetic tree for L. equi based on 1454 core genes showed that the independent branch of strain IMAU81196 was far from the other three strains. In terms of genomic characteristics, single-nucleotide polymorphism (SNP) sites, rapid annotation using subsystem technology (RAST), carbohydrate activity enzymes (CAZy), and predictions of prophage, we showed that strain L. equi JCM 10991T and strain DSM 15833T are not equivalent strains.It is worth mentioning thatthestrain of L. equi has numerous enzymes related to cellulose degradation, and each L. equi strain investigated contained at least one protophage. We speculate that this is the reason why these strains are adapted to the intestinal environment of horses. These results provide new research directions for the future.

Draft Genome Sequence of Oceanobacillus jordanicus Strain GSFE11, a Halotolerant Plant Growth-Promoting Bacterial Endophyte Isolated From the Jordan Valley

The bacterium Oceanobacillus jordanicus strain GSFE11 is a halotolerant endophyte isolated from sterilized roots of Durum wheat (Triticum turgidum ssp. Durum) growing in hot and arid environments of Ghor Safi area in the Jordan Valley. The draft genome sequence and annotation of this plant growth-promoting endophytic bacterium are reported in this study. The draft genome sequence of Oceanobacillus jordanicus strain GSFE11 has 3 839 208 bp with a G + C content of 39.09%. A total of 3893 protein-coding genes and 68 RNA coding genes were predicted. Several putative genes that are involved in secretion and delivery systems, transport, adhesion, motility, membrane proteins, plant cell wall modification, and detoxification were identified, some are characteristics of endophytes lifestyle including genes that are involved in metabolism of carbohydrate, genes for xylose, fructose and chitin utilization, quinone cofactors biosynthesis, genes associated with nitrogen, sulfur, phosphate and iron acquisition, in addition to genes involved in the biosynthesis of plant hormone auxin. This study highlights the importance of using genome analysis and phylogenomic analysis to resolve the differences between closely related species, such analysis showed Oceanobacillus jordanicus strain GSFE11 to be a new species closely related to Oceanobacillus picturae (genome size 3.67 Mb), Oceanobacillus jordanicus has higher a number of predicted genes compared with Oceanobacillus picturae (3961 genes vs 3823 genes).

Virulence Potential and Treatment Options of Multidrug-Resistant (MDR) Acinetobacter baumannii

Acinetobacter baumannii is an opportunistic pathogen which is undoubtedly known for a high rate of morbidity and mortality in hospital-acquired infections. A. baumannii causes life-threatening infections, including; ventilator-associated pneumonia (VAP), meningitis, bacteremia, and wound and urinary tract infections (UTI). In 2017, the World Health Organization listed A. baumannii as a priority-1 pathogen. The prevalence of A. baumannii infections and outbreaks emphasizes the direct need for the use of effective therapeutic agents for treating such infections. Available antimicrobials, such as; carbapenems, tigecycline, and colistins have insufficient effectiveness due to the appearance of multidrug-resistant strains, accentuating the need for alternative and novel therapeutic remedies. To understand and overcome this menace, the knowledge of recent discoveries on the virulence factors of A. baumannii is needed. Herein, we summarized the role of various virulence factors, including; outer membrane proteins, efflux pumps, biofilm, penicillin-binding proteins, and siderophores/iron acquisition systems. We reviewed the recent scientific literature on different A. baumannii virulence factors and the effective antimicrobial agents for the treatment and management of bacterial infections.

An evaluation of the species and subspecies of the genus Salmonella with whole genome sequence data: Proposal of type strains and epithets for novel S. enterica subspecies VII, VIII, IX, X and XI

Species and subspecies within the Salmonella genus have been defined for public health purposes by biochemical properties; however, reference laboratories have increasingly adopted sequence-based, and especially whole genome sequence (WGS), methods for surveillance and routine identification. This leads to potential disparities in subspecies definitions, routine typing, and the ability to detect novel subspecies. A large-scale analysis of WGS data from the routine sequencing of clinical isolates was employed to define and characterise Salmonella subspecies population structure, demonstrating that the Salmonella species and subspecies were genetically distinct, including those previously identified through phylogenetic approaches, namely: S. enterica subspecies londinensis (VII), subspecies brasiliensis (VIII), subspecies hibernicus (IX) and subspecies essexiensis (X). The analysis also identified an additional novel subspecies, reptilium (XI). Further, these analyses indicated that S. enterica subspecies arizonae (IIIa) isolates were divergent from the other S. enterica subspecies, which clustered together and, on the basis of ANI analysis, subspecies IIIa was sufficiently distinct to be classified as a separate species, S. arizonae. Multiple phylogenetic and statistical approaches generated congruent results, suggesting that the proposed species and subspecies structure was sufficiently biologically robust for routine application. Biochemical analyses demonstrated that not all subspecies were distinguishable by these means and that biochemical approaches did not capture the genomic diversity of the genus. We recommend the adoption of standardised genomic definitions of species and subspecies and a genome sequence-based approach to routine typing for the identification and definition of novel subspecies.

•

A large-scale analysis of genomic data demonstrate Salmonella species and subspecies are genetically distinct.

•

Biochemical analysis does not capture the genomic diversity of the Salmonella genus but routine species and subspecies identification can be achieved with rMLST

•

Average Nucleotide Identify (ANI) with a 95% criteria was suitable to distinguish species and 98% to distinguish subspecies.

•

Five novel S. enteric subspecies (VII-XI) type strains are defined.

•

Reclassification of S. arizonae as a separate species is recommended.

The genetic and ecophysiological diversity of Microcystis

Microcystis is a cyanobacterium that forms toxic blooms in freshwater ecosystems around the world. Biological variation among taxa within the genus is apparent through genetic and phenotypic differences between strains and via the spatial and temporal distribution of strains in the environment, and this fine-scale diversity exerts strong influence over bloom toxicity. Yet we do not know how varying traits of Microcystis strains govern their environmental distribution, the tradeoffs and links between these traits, or how they are encoded at the genomic level. Here we synthesize current knowledge on the importance of diversity within Microcystis and on the genes and traits that likely underpin ecological differentiation of taxa. We briefly review spatial and environmental patterns of Microcystis diversity in the field and genetic evidence for cohesive groups within Microcystis. We then compile data on strain-level diversity regarding growth responses to environmental conditions and explore evidence for variation of community interactions across Microcystis strains. Potential links and tradeoffs between traits are identified and discussed. The resulting picture, while incomplete, highlights key knowledge gaps that need to be filled to enable new models for predicting strain-level dynamics, which influence the development, toxicity and cosmopolitan nature of Microcystis blooms.

Ammonia oxidation at pH 2.5 by a new gammaproteobacterial ammonia-oxidizing bacterium

Ammonia oxidation was considered impossible under highly acidic conditions, as the protonation of ammonia leads to decreased substrate availability and formation of toxic nitrogenous compounds. Recently, some studies described archaeal and bacterial ammonia oxidizers growing at pH as low as 4, while environmental studies observed nitrification at even lower pH values. In this work, we report on the discovery, cultivation, and physiological, genomic, and transcriptomic characterization of a novel gammaproteobacterial ammonia-oxidizing bacterium enriched via continuous bioreactor cultivation from an acidic air biofilter that was able to grow and oxidize ammonia at pH 2.5. This microorganism has a chemolithoautotrophic lifestyle, using ammonia as energy source. The observed growth rate on ammonia was 0.196 day^-1, with a doubling time of 3.5 days. The strain also displayed ureolytic activity and cultivation with urea as ammonia source resulted in a growth rate of 0.104 day^-1 and a doubling time of 6.7 days. A high ammonia affinity (K_m(app) = 147 ± 14 nM) and high tolerance to toxic nitric oxide could represent an adaptation to acidic environments. Electron microscopic analysis showed coccoid cell morphology with a large amount of intracytoplasmic membrane stacks, typical of gammaproteobacterial ammonia oxidizers. Furthermore, genome and transcriptome analysis showed the presence and expression of diagnostic genes for nitrifiers (amoCAB, hao, nor, ure, cbbLS), but no nirK was identified. Phylogenetic analysis revealed that this strain belonged to a novel bacterial genus, for which we propose the name "Candidatus Nitrosacidococcus tergens" sp. RJ19.

Phylogenomic Framework for Taxonomic Delineation of Paracoccus spp. and Exploration of Core-Pan Genome

The taxonomic classification of metabolically versatile Paracoccus spp. has been so far performed using polyphasic approach. The topology of single gene phylogenies, however, has highlighted ambiguous species assignments. In the present study, genome based multi-gene phylogenies and overall genome related index were used for species threshold assessment. Comprehensive phylogenomic analysis of Paracoccus genomes (n = 103) showed concordant clustering of strains across multi-gene marker set phylogenies (nMC = 0.08-0.14); as compared to 16S rDNA phylogeny (nMC = 0.37-0.42) suggesting robustness of multi gene phylogenies in drawing phylogenetic inferences. Functional gene content distribution across the genus showed that only 1.7% gene content constitutes the core genome highlighting the significance of extensive genomic variability in the evolution of Paracoccus spp. Further, genome metrics were used to validate characterized strains, identifying classification anomalies (n = 13), and based on this, genome derived taxonomic amendments were notified in present study. Conclusively, validated metric tools can be employed on whole genome sequences, including draft assemblies, for the assessment and assignment of uncharacterized strains and species level ascription of newly isolated Paracoccus strains in future.

The genus bifidobacterium: From genomics to functionality of an important component of the mammalian gut microbiota running title: Bifidobacterial adaptation to and interaction with the host

Members of the genus Bifidobacterium are dominant and symbiotic inhabitants of the mammalian gastrointestinal tract. Being vertically transmitted, bifidobacterial host colonization commences immediately after birth and leads to a phase of host infancy during which bifidobacteria are highly prevalent and abundant to then transit to a reduced, yet stable abundance phase during host adulthood. However, in order to reach and stably colonize their elective niche, i.e. the large intestine, bifidobacteria have to cope with a multitude of oxidative, osmotic and bile salt/acid stress challenges that occur along the gastrointestinal tract (GIT). Concurrently, bifidobacteria not only have to compete with the myriad of other gut commensals for nutrient acquisition, but they also require protection against bacterial viruses. In this context, Next-Generation Sequencing (NGS) techniques, allowing large-scale comparative and functional genome analyses have helped to identify the genetic strategies that bifidobacteria have developed in order to colonize, survive and adopt to the highly competitive mammalian gastrointestinal environment. The current review is aimed at providing a comprehensive overview concerning the molecular strategies on which bifidobacteria rely to stably and successfully colonize the mammalian gut.

Grimontia sedimenti sp. nov., isolated from benthic sediments near coral reefs south of Kuwait

A Gram-stain-negative, rod and rod-curved shaped motile bacterium designated strain S25^T was obtained from benthic sediment collected near the Kubbar Island coral reefs south of Kuwait. Phenotypic analysis revealed that strain S25^T was slightly halophilic, mesophilic and facultative anaerobic, fermenting d-glucose, d-ribose, d-mannose, d-mannitol, maltose, fructose, gentiobiose, cellobiose, melibiose, trehalose and sucrose. It was positive for oxidase and indole production and negative for arginine dihydrolase and lysine and ornithine decarboxylases. It contained C_16 : 1  ω7c/C_16 : 1  ω6c (summed feature 3), C_18 : 1 ω7c (summed feature 8) and C_16 : 0 as the major fatty acids. Strain S25^T grew optimally at 30 °C and pH 8 in the presence of 3 % (w/v) NaCl. Phylogenetic analysis based on 16S rRNA sequences revealed that strain S25^T is related to species of the genus Grimontia, having 99.15 % similarity to 'Grimontia indica' AK16^T, 99.08 % to Grimontia celer 96-237^T and 98.66 % to Grimontia marina IMCC 5001^T. The DNA G+C content was 48.8 mol% and the full genome analysis for the strain S25^T showed that the bacterium has a genome size of 5 158 621 bp and contains 4730 predicted protein-encoding genes. The average nucleotide identity values between the S25^T genome and the genomes of its nearest matches ranged between 81.39 and 94.16 %. The strain was distinguishable from the phylogenetically related genera through differences in several phenotypic properties. On the basis of the phenotypic, phylogenetic and genetic data, strain S25^T represents a novel species in the genus Grimontia, for which the name Grimontia sedimenti sp. nov. is proposed. The type strain of Grimontia sedimenti is S25^T (=DSM 28878^T=LMG 28315^T).

Divergence in environmental adaptation between terrestrial clades of the Bacillus cereus group

The Bacillus cereus group encompasses beneficial and harmful species in diverse niches and has a much debated taxonomy. Investigating whether selection has led to ecological divergence between phylogenetic clades can help understand the basis of speciation, and has implications for predicting biological safety across this group. Using three most terrestrial species in this group (B. cereus, Bacillus thuringiensis and Bacillus mycoides) we charactererized ecological specialization in terms of resource use, thermal adaptation and fitness in different environmental conditions and tested whether taxonomic species or phylogenetic clade best explained phenotypic variation. All isolates grew vigorously in protein rich media and insect cadavers, but exploitation of soil or plant derived nutrients was similarly weak for all. For B. thuringiensis and B. mycoides, clade and taxonomic species were important predictors of relative fitness in insect infections. Fully psychrotolerant isolates could outcompete B. thuringiensis in insects at low temperature, although psychrotolerance predicted growth in artificial media better than clade. In contrast to predictions, isolates in the Bacillus anthracis clade had sub-optimal growth at 37°C. The common ecological niche in these terrestrial B. cereus species is the ability to exploit protein rich resources such as cadavers. However, selection has led to different phylogenetic groups developing different strategies for accessing this resource. Thus, clades, as well as traditional taxonomic phenotypes, predict biologically important traits.

Characterization and complete genome analysis of Bacillus velezensis CB6 revealed ATP synthase subunit α against foodborne pathogens

Given the serious threat of foodborne multidrug-resistant bacteria to animals and humans, finding an effective antibacterial compound has always been an important topic for scientists. Here, from the soil of Changbaishan, we have identified a bacterium that can inhibit the growth of Staphylococcus aureus. Nr genome database analysis and phylogenetic analysis showed that strain CB6 belongs to Bacillus velezensis. We found that the crude extract of strain CB6 has broad-spectrum antibacterial activity against foodborne pathogens. In addition, we showed that the crude extract loses antibacterial activity after treatment with papain. Next, strain CB6 was purified using ammonium sulfate precipitation, a Sephadex G-75 gel filtration column and high-performance liquid chromatography system (HPLC). Liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis indicated that the antibacterial compound was the protein ATP synthase subunit α (ATP-1), with a molecular weight of 55.397 KDa. Moreover, we reported the complete genome sequence of strain CB6, which is composed of a unique circular 3,963,507 bp chromosome with 3749 coding genes and a G + C content of 46.53%. The genome contained 12 gene clusters with antibacterial functions, which constituted over 20.947% of the complete genome. Of note, the amino acid sequence encoding the ATP-1 protein in the strain CB6 genome was identified. In addition to these findings, we speculate that the ATP-1 protein may provide energy for secondary metabolites, which in turn will improve the antibacterial activity of the secondary metabolites. All the above important features make the ATP-1 as a potential candidate for the development of new antibacterial drugs and food preservatives in the future.

Towards Genomic Criteria for Delineating Fungal Species

The discussion of fungal species delineation has yet to reach a consensus, despite the advancements in technology, which helped modernise traditional approaches. In particular, the phylogenetic species concept was one of the tools that has been used with considerable success across the fungal kingdom. The fast rise of fungal genomics provides an unprecedented opportunity to expand measuring the relatedness of fungal strains to the level of whole genomes. However, the use of genomic information in taxonomy has only just begun, and few methodological guidelines have been suggested so far. Here, a simple approach of computationally measuring genomic distances and their use as a standard for species delineation is investigated. A fixed threshold genomic distance calculated by the quick and easy-to-use tools Mash and Dashing proved to be an unexpectedly widely applicable and robust criterion for determining whether two genomes belong to the same or to different species. The accuracy of species delineation in an uncurated dataset of GenBank fungal genomes was close to 90%—and exceeded 90% with minimal curation. As expected, the discriminative power of this approach was lower at higher taxonomic ranks, but still significantly larger than zero. Simple instructions for calculation of a genomic distance between two genomes and species similarity thresholds at different k-mer sizes are suggested. The calculation of genomic distance is identified as a powerful approach for delineating fungal species and is proposed—not as the only criterion—but as an additional tool in the versatile toolbox of fungal taxonomy.

Methylococcus geothermalis sp. nov., a methanotroph isolated from a geothermal field in the Republic of Korea

A Gram-stain-negative, aerobic, non-motile and coccoid methanotroph, strain IM1^T, was isolated from hot spring soil. Cells of strain IM1^T were catalase-negative, oxidase-positive and displayed a characteristic intracytoplasmic membrane arrangement of type I methanotrophs. The strain possessed genes encoding both membrane-bound and soluble methane monooxygenases and grew only on methane or methanol. The strain was capable of growth at temperatures between 15 and 48 °C (optimum, 30-45 °C) and pH values between pH 4.8 and 8.2 (optimum, pH 6.2-7.0). Based on phylogenetic analysis of 16S rRNA gene and PmoA sequences, strain IM1^T was demonstrated to be affiliated to the genus Methylococcus. The 16S rRNA gene sequence of this strain was most closely related to the sequences of an uncultured bacterium clone FD09 (100 %) and a partially described cultured Methylococcus sp. GDS2.4 (99.78 %). The most closely related taxonomically described strains were Methylococcus capsulatus Texas^T (97.92 %), Methylococcus capsulatus Bath (97.86 %) and Methyloterricola oryzae 73a^T (94.21 %). Strain IM1^T shared average nucleotide identity values of 85.93 and 85.62 % with Methylococcus capsulatus strains Texas^T and Bath, respectively. The digital DNA-DNA hybridization value with the closest type strain was 29.90 %. The DNA G+C content of strain IM1^T was 63.3 mol% and the major cellular fatty acids were C_16 : 0 (39.0 %), C_16 : 1  ω7c (24.0 %), C_16 : 1  ω6c (13.6 %) and C_16 : 1  ω5c (12.0 %). The major ubiquinone was methylene-ubiquinone-8. On the basis of phenotypic, genetic and phylogenetic data, strain IM1^T represents a novel species of the genus Methylococcus for which the name Methylococcus geothermalis sp. nov. is proposed, with strain IM1^T (=JCM 33941^T=KCTC 72677^T) as the type strain.

Bradyrhizobium campsiandrae sp. nov., a nitrogen-fixing bacterial strain isolated from a native leguminous tree from the Amazon adapted to flooded conditions

The nitrogen-fixing bacterial strain UFLA 01-1174^T was isolated from nodules of Campsiandra laurilifolia Benth. originating from the Amazon region, Brazil. Its taxonomic position was defined using a polyphasic approach. Analysis of the 16S rRNA gene placed the strain in the Bradyrhizobium genus, the closest species being B. guangdongense CCBAU 51649^T and B. guangzhouense CCBAU 51670^T, both with 99.8% similarity. Multilocus sequence analysis (MLSA) of recA, gyrB, glnII, rpoB, atpD, and dnaK indicated that UFLA 01-1174^T is a new species, most closely related to B. stylosanthis BR 446^T (94.4%) and B. manausense BR 3351^T (93.7%). Average nucleotide identity (ANI) differentiated UFLA 01-1174^T from the closest species with values lower than 90%. The G + C content in the DNA of UFLA 01-1174^T is 63.6 mol%. Based on this data, we conclude that the strain represents a new species. The name proposed is Bradyrhizobium campsiandrae, with UFLA 01-1174^T (= INPA 394B^T = LMG 10099^T) as type strain.

Exploring the Diversity Within the Genus Francisella - An Integrated Pan-Genome and Genome-Mining Approach

Pan-genome analysis is a powerful method to explore genomic heterogeneity and diversity of bacterial species. Here we present a pan-genome analysis of the genus Francisella, comprising a dataset of 63 genomes and encompassing clinical as well as environmental isolates from distinct geographic locations. To determine the evolutionary relationship within the genus, we performed phylogenetic whole-genome studies utilizing the average nucleotide identity, average amino acid identity, core genes and non-recombinant loci markers. Based on the analyses, the phylogenetic trees obtained identified two distinct clades, A and B and a diverse cluster designated C. The sizes of the pan-, core-, cloud-, and shell-genomes of Francisella were estimated and compared to those of two other facultative intracellular pathogens, Legionella and Piscirickettsia. Francisella had the smallest core-genome, 692 genes, compared to 886 and 1,732 genes for Legionella and Piscirickettsia respectively, while the pan-genome of Legionella was more than twice the size of that of the other two genera. Also, the composition of the Francisella Type VI secretion system (T6SS) was analyzed. Distinct differences in the gene content of the T6SS were identified. In silico approaches performed to identify putative substrates of these systems revealed potential effectors targeting the cell wall, inner membrane, cellular nucleic acids as well as proteins, thus constituting attractive targets for site-directed mutagenesis. The comparative analysis performed here provides a comprehensive basis for the assessment of the phylogenomic relationship of members of the genus Francisella and for the identification of putative T6SS virulence traits.

New Insight into Antimicrobial Compounds from Food and Marine-Sourced Carnobacterium Species through Phenotype and Genome Analyses

Carnobacteriummaltaromaticum and Carnobacteriumdivergens, isolated from food products, are lactic acid bacteria known to produce active and efficient bacteriocins. Other species, particularly those originating from marine sources, are less studied. The aim of the study is to select promising strains with antimicrobial potential by combining genomic and phenotypic approaches on large datasets comprising 12 Carnobacterium species. The biosynthetic gene cluster (BGCs) diversity of 39 publicly available Carnobacterium spp. genomes revealed 67 BGCs, distributed according to the species and ecological niches. From zero to six BGCs were predicted per strain and classified into four classes: terpene, NRPS (non-ribosomal peptide synthetase), NRPS-PKS (hybrid non-ribosomal peptide synthetase-polyketide synthase), RiPP (ribosomally synthesized and post-translationally modified peptide). In parallel, the antimicrobial activity of 260 strains from seafood products was evaluated. Among the 60% of active strains, three genomes were sequenced and submitted to a dereplication process. C. inhibens MIP2551 produced a high amountof H₂O₂, probably thanks to the presence of four oxidase-encoding genes. C. maltaromaticum EBP3019 and SF668 strains were highly efficient against Listeria monocytogenes. A new extracellular 16 kDa unmodified bacteriocin in the EBP3019 strain and five different bacteriocins in SF668 were highlighted. In this study, the overview of antimicrobial BGC and inhibitory activities of Carnobacterium spp. allowed the prediction of potential innovative natural products that could be relevant for biotechnological applications.

Multidrug Resistance (MDR) and Collateral Sensitivity in Bacteria, with Special Attention to Genetic and Evolutionary Aspects and to the Perspectives of Antimicrobial Peptides-A Review

Antibiotic poly-resistance (multidrug-, extreme-, and pan-drug resistance) is controlled by adaptive evolution. Darwinian and Lamarckian interpretations of resistance evolution are discussed. Arguments for, and against, pessimistic forecasts on a fatal “post-antibiotic era” are evaluated. In commensal niches, the appearance of a new antibiotic resistance often reduces fitness, but compensatory mutations may counteract this tendency. The appearance of new antibiotic resistance is frequently accompanied by a collateral sensitivity to other resistances. Organisms with an expanding open pan-genome, such as Acinetobacter baumannii, Pseudomonas aeruginosa, and Klebsiella pneumoniae, can withstand an increased number of resistances by exploiting their evolutionary plasticity and disseminating clonally or poly-clonally. Multidrug-resistant pathogen clones can become predominant under antibiotic stress conditions but, under the influence of negative frequency-dependent selection, are prevented from rising to dominance in a population in a commensal niche. Antimicrobial peptides have a great potential to combat multidrug resistance, since antibiotic-resistant bacteria have shown a high frequency of collateral sensitivity to antimicrobial peptides. In addition, the mobility patterns of antibiotic resistance, and antimicrobial peptide resistance, genes are completely different. The integron trade in commensal niches is fortunately limited by the species-specificity of resistance genes. Hence, we theorize that the suggested post-antibiotic era has not yet come, and indeed might never come.

Genome Comparison Identifies Different Bacillus Species in a Bast Fibre-Retting Bacterial Consortium and Provides Insights into Pectin Degrading Genes

Retting of bast fibres requires removal of pectin, hemicellulose and other non-cellulosic materials from plant stem tissues by a complex microbial community. A microbial retting consortium with high-efficiency pectinolytic bacterial strains is effective in reducing retting-time and enhancing fibre quality. We report comprehensive genomic analyses of three bacterial strains (PJRB 1, 2 and 3) of the consortium and resolve their taxonomic status, genomic features, variations, and pan-genome dynamics. The genome sizes of the strains are ~3.8 Mb with 3729 to 4002 protein-coding genes. Detailed annotations of the protein-coding genes revealed different carbohydrate-degrading CAZy classes viz. PL1, PL9, GH28, CE8, and CE12. Phylogeny and structural features of pectate lyase proteins of PJRB strains divulge their functional uniqueness and evolutionary convergence with closely related Bacillus strains. Genome-wide prediction of genomic variations revealed 12461 to 67381 SNPs, and notably many unique SNPs were localized within the important pectin metabolism genes. The variations in the pectate lyase genes possibly contribute to their specialized pectinolytic function during the retting process. These findings encompass a strong foundation for fundamental and evolutionary studies on this unique microbial degradation of decaying plant material with immense industrial significance. These have preponderant implications in plant biomass research and food industry, and also posit application in the reclamation of water pollution from plant materials.

Ontologically simple theories do not indicate the true nature of complex biological systems: three test cases

A longstanding philosophical premise perceives simplicity as a desirable attribute of scientific theories. One of several raised justifications for this notion is that simple theories are more likely to indicate the true makeup of natural systems. Qualitatively parsimonious hypotheses and theories keep to a minimum the number of different postulated entities within a system. Formulation of such ontologically simple working hypotheses proved to be useful in the experimental probing of narrowly defined bio systems. It is less certain, however, whether qualitatively parsimonious theories are effective indicators of the true nature of complex biological systems. This paper assesses the success of ontologically simple theories in envisaging the makeup of three complex systems in bacteriology, immunology, and molecular biology. Evidence shows that parsimonious theories completely misconstrued the actual ontologically complex constitutions of the three examined systems. Since evolution and selective pressures typically produce ontologically intricate rather than simple bio systems, qualitatively parsimonious theories are mostly inapt indicators of the true nature of complex biological systems.

Enterobacter wuhouensis sp. nov. and Enterobacter quasihormaechei sp. nov. recovered from human sputum

Two novel strains of members of the genus Enterobacter, WCHEs120002^T and WCHEs120003^T, were recovered from the sputum of two patients at a hospital in PR China in 2017. The strains were Gram-stain-negative, facultatively anaerobic, motile and non-spore-forming. The two strains were subjected to whole-genome sequencing. Phylogenetic analysis based on core genes of type strains of species of the family Enterobacteriaceae revealed that the two strains belonged to the genus Enterobacter but were distinct from any previously known species of the genus. Both average nucleotide identity and in silico DNA-DNA hybridization values between strains WCHEs120002^T and WCHEs120003^T and type strains of all known species of the genus Enterobacter were lower than the recommended thresholds of 95 and 70  %, respectively, for species delineation. The major fatty acids of the two strains were C_16 : 0, C_17 : 0 cyclo and C_18:1ω7c, which are similar to those of other species of the genus Enterobacter. Genomic DNA G+C contents of strains WCHEs120002^T and WCHEs120003^T were 56.09 and 55.91 mol%, respectively. WCHEs120002^T ferments melibiose and sucrose but is negative for d-sorbitol and methyl-α-d-mannopyranoside reactions, which distinguish it from all other species of the genus Enterobacter. WCHEs120003^T can be differentiated from other species of the genus Enterobacter by its ability to ferment potassium gluconate and its negative reactions for d-sorbitol and l-fucose. Genotypic and phenotypic characteristics indicate that strains WCHEs120002^T and WCHEs120003^T represent two novel species of the genus Enterobacter, for which the names Enterobacter wuhouensis sp. nov. and Enterobacter quasihormaechei sp. nov. are proposed, respectively. The type strain of E. wuhouensis sp. nov. is WCHEs120002^T (=GDMCC1.1569^T=NCTC 14273^T) and the type strain of E. quasihormaechei sp. nov. is WCHEs120003^T (=GDMCC1.1568^T=NCTC 14274^T).

Complete Genome Sequence of Pseudomonas psychrotolerans CS51, a Plant Growth-Promoting Bacterium, Under Heavy Metal Stress Conditions

In the current study, we aimed to elucidate the plant growth-promoting characteristics of Pseudomonas psychrotolerans CS51 under heavy metal stress conditions (Zn, Cu, and Cd) and determine the genetic makeup of the CS51 genome using the single-molecule real-time (SMRT) sequencing technology of Pacific Biosciences. The results revealed that inoculation with CS51 induced endogenous indole-3-acetic acid (IAA) and gibberellins (GAs), which significantly enhanced cucumber growth (root shoot length) and increased the heavy metal tolerance of cucumber plants. Moreover, genomic analysis revealed that the CS51 genome consisted of a circular chromosome of 5,364,174 base pairs with an average G+C content of 64.71%. There were around 4774 predicted protein-coding sequences (CDSs) in 4859 genes, 15 rRNA genes, and 67 tRNA genes. Around 3950 protein-coding genes with function prediction and 733 genes without function prediction were identified. Furthermore, functional analyses predicted that the CS51 genome could encode genes required for auxin biosynthesis, nitrate and nitrite ammonification, the phosphate-specific transport system, and the sulfate transport system, which are beneficial for plant growth promotion. The heavy metal resistance of CS51 was confirmed by the presence of genes responsible for cobalt-zinc-cadmium resistance, nickel transport, and copper homeostasis in the CS51 genome. The extrapolation of the curve showed that the core genome contained a minimum of 2122 genes (95% confidence interval = 2034.24 to 2080.215). Our findings indicated that the genome sequence of CS51 may be used as an eco-friendly bioresource to promote plant growth in heavy metal-contaminated areas.

Spatially distinct physiology of Bacteroides fragilis within the proximal colon of gnotobiotic mice

A complex microbiota inhabits various microenvironments of the gut, with some symbiotic bacteria having evolved traits to invade the epithelial mucus layer and reside deep within intestinal tissue of animals. Whether these distinct bacterial communities across gut biogeographies exhibit divergent behaviors remains largely unknown. Global transcriptomic analysis to investigate microbial physiology in specific mucosal niches has been hampered technically by overabundance of host RNA. Herein, we employed hybrid selection RNA sequencing (hsRNA-Seq) to enable detailed spatial transcriptomic profiling of a prominent human commensal as it colonizes the colonic lumen, mucus or epithelial tissue of mice. Compared to conventional RNA-Seq, hsRNA-Seq increased reads mapping to the Bacteroides fragilis genome by 48- and 154-fold in mucus and tissue, respectively, allowing for high fidelity comparisons across biogeographic sites. Near the epithelium, B. fragilis up-regulated numerous genes involved in protein synthesis, indicating that bacteria inhabiting the mucosal niche are metabolically active. Further, a specific sulfatase (BF3086) and glycosyl hydrolase (BF3134) were highly induced in mucus and tissue compared to bacteria in the lumen. In-frame deletion of these genes impaired in vitro growth on mucus as a carbon source, as well as mucosal colonization of mice. Mutants in either B. fragilis gene displayed a fitness defect in competing for colonization against bacterial challenge, revealing the importance of site-specific gene expression for robust host-microbial symbiosis. As a versatile tool, hsRNA-Seq can be deployed to explore the in vivo spatial physiology of numerous bacterial pathogens or commensals.

Genome sequence-based criteria for demarcation and definition of species in the genus Rickettsia

Over recent years, genomic information has increasingly been used for prokaryotic species definition and classification. Genome sequence-based alternatives to the gold standard DNA–DNA hybridization (DDH) relatedness have been developed, notably average nucleotide identity (ANI), which is one of the most useful measurements for species delineation in the genomic era. However, the strictly intracellar lifestyle, the few measurable phenotypic properties and the low level of genetic heterogeneity made the current standard genomic criteria for bacterial species definition inapplicable to Rickettsia species. We evaluated a range of whole genome sequence (WGS)-based taxonomic parameters to develop guidelines for the classification of Rickettsia isolates at genus and species levels. By comparing the degree of similarity of 74 WGSs from 31 Rickettsia species and 61 WGSs from members of three closely related genera also belonging to the order Rickettsiales ( Orientia , 11 genomes; Ehrlichia , 22 genomes; and Anaplasma , 28 genomes) using digital DDH (dDDh) and ANI by orthology (OrthoANI) parameters, we demonstrated that WGS-based taxonomic information, which is easy to obtain and use, can serve for reliable classification of isolates within the Rickettsia genus and species. To be classified as a member of the genus Rickettsia , a bacterial isolate should exhibit OrthoANI values with any Rickettsia species with a validly published name of ≥83.63 %. To be classified as a new Rickettsia species, an isolate should not exhibit more than any of the following degrees of genomic relatedness levels with the most closely related species: >92.30 and >99.19 % for the dDDH and OrthoANI values, respectively. When applied to four rickettsial isolates of uncertain status, the above-described thresholds enabled their classification as new species in one case. Thus, we propose WGS-based guidelines to efficiently delineate Rickettsia species, with OrthoANI and dDDH being the most accurate for classification at the genus and species levels, respectively.

Paenibacillus odorifer, the Predominant Paenibacillus Species Isolated from Milk in the United States, Demonstrates Genetic and Phenotypic Conservation of Psychrotolerance but Clade-Associated Differences in Nitrogen Metabolic Pathways

Although Paenibacillus species isolates are frequently isolated from pasteurized fluid milk, the link between the genetic diversity and phenotypic characteristics of these isolates was not well understood, especially as some Bacillales isolated from milk are unable to grow at refrigeration temperatures. Our data demonstrate that Paenibacillus spp. isolated from fluid milk represent tremendous interspecies diversity, with P. odorifer being the predominant Paenibacillus sp. isolated. Furthermore, genetic and phenotypic data support that P. odorifer is well suited to transition from a soil-dwelling environment, where nitrogen fixation (and other nitrate/nitrite reduction pathways present only in clade A) may facilitate growth, to fluid milk, where its multiple cold shock-associated adaptations enable it to grow at refrigeration temperatures throughout the storage of milk. Therefore, efforts to reduce bacterial contamination of milk will require a systematic approach to reduce P. odorifer contamination of raw milk.

Paenibacillus is a spore-forming bacterial genus that is frequently isolated from fluid milk and is proposed to play a role in spoilage. To characterize the genetic and phenotypic diversity of Paenibacillus spp., we first used rpoB allelic typing data for a preexisting collection of 1,228 Paenibacillus species isolates collected from raw and processed milk, milk products, and dairy environmental sources. Whole-genome sequencing (WGS) and average nucleotide identity by BLAST (ANIb) analyses performed for a subset of 58 isolates representing unique and overrepresented rpoB allelic types in the collection revealed that these isolates represent 21 different Paenibacillus spp., with P. odorifer being the predominant species. Further genomic characterization of P. odorifer isolates identified two distinct phylogenetic clades, clades A and B, which showed significant overrepresentation of 172 and 164 ortholog clusters and 94 and 52 gene ontology (GO) terms, respectively. While nitrogen fixation genes were found in both clades, multiple genes associated with nitrate and nitrite reduction were overrepresented in clade A isolates; additional phenotypic testing demonstrated that nitrate reduction is specific to isolates in clade A. Hidden Markov models detected 9 to 10 different classes of cold shock-associated genetic elements in all P. odorifer isolates. Phenotypic testing revealed that all isolates tested here can grow in skim milk broth at 6°C, suggesting that psychrotolerance is conserved in P. odorifer. Overall, our data suggest that Paenibacillus spp. isolated from milk in the United States represent broad genetic diversity, which may provide challenges for targeted-control strategies aimed at reducing fluid milk spoilage.

IMPORTANCE Although Paenibacillus species isolates are frequently isolated from pasteurized fluid milk, the link between the genetic diversity and phenotypic characteristics of these isolates was not well understood, especially as some Bacillales isolated from milk are unable to grow at refrigeration temperatures. Our data demonstrate that Paenibacillus spp. isolated from fluid milk represent tremendous interspecies diversity, with P. odorifer being the predominant Paenibacillus sp. isolated. Furthermore, genetic and phenotypic data support that P. odorifer is well suited to transition from a soil-dwelling environment, where nitrogen fixation (and other nitrate/nitrite reduction pathways present only in clade A) may facilitate growth, to fluid milk, where its multiple cold shock-associated adaptations enable it to grow at refrigeration temperatures throughout the storage of milk. Therefore, efforts to reduce bacterial contamination of milk will require a systematic approach to reduce P. odorifer contamination of raw milk.

Genome-resolved metagenomics analysis provides insights into the ecological role of Thaumarchaeota in the Amazon River and its plume

Background

Thaumarchaeota are abundant in the Amazon River, where they are the only ammonia-oxidizing archaea. Despite the importance of Thaumarchaeota, little is known about their physiology, mainly because few isolates are available for study. Therefore, information about Thaumarchaeota was obtained primarily from genomic studies. The aim of this study was to investigate the ecological roles of Thaumarchaeota in the Amazon River and the Amazon River plume.

Results

The archaeal community of the shallow in Amazon River and its plume is dominated by Thaumarchaeota lineages from group 1.1a, which are mainly affiliated to Candidatus Nitrosotenuis uzonensis, members of order Nitrosopumilales, Candidatus Nitrosoarchaeum, and Candidatus Nitrosopelagicus sp. While Thaumarchaeota sequences have decreased their relative abundance in the plume, Candidatus Nitrosopelagicus has increased. One genome was recovered from metagenomic data of the Amazon River (ThauR71 [1.05 Mpb]), and two from metagenomic data of the Amazon River plume (ThauP25 [0.94 Mpb] and ThauP41 [1.26 Mpb]). Phylogenetic analysis placed all three Amazon genome bins in Thaumarchaeota Group 1.1a. The annotation revealed that most genes are assigned to the COG subcategory coenzyme transport and metabolism. All three genomes contain genes involved in the hydroxypropionate/hydroxybutyrate cycle, glycolysis, tricarboxylic acid cycle, oxidative phosphorylation. However, ammonia-monooxygenase genes were detected only in ThauP41 and ThauR71. Glycoside hydrolases and auxiliary activities genes were detected only in ThauP25.

Conclusions

Our data indicate that Amazon River is a source of Thaumarchaeota, where these organisms are important for primary production, vitamin production, and nitrification.

Comparative Genomics Analysis of Lactobacillus mucosae from Different Niches

The potential probiotic benefits of Lactobacillus mucosae have received increasing attention. To investigate the genetic diversity of L. mucosae, comparative genomic analyses of 93 strains isolated from different niches (human and animal gut, human vagina, etc.) and eight strains of published genomes were conducted. The results showed that the core genome of L. mucosae mainly encoded translation and transcription, amino acid biosynthesis, sugar metabolism, and defense function while the pan-genomic curve tended to be close. The genetic diversity of L. mucosae mainly reflected in carbohydrate metabolism and immune/competitive-related factors, such as exopolysaccharide (EPS), enterolysin A, and clustered regularly interspaced short palindromic repeats (CRISPR)-Cas. It was worth noting that this research firstly predicted the complete EPS operon shared among L. mucosae. Additionally, the type IIIA CRISPR-Cas system was discovered in L. mucosae for the first time. This work provided new ideas for the study of this species.

Characterization of Bradyrhizobium strains indigenous to Western Australia and South Africa indicates remarkable genetic diversity and reveals putative new species

Bradyrhizobium are N₂-fixing microsymbionts of legumes with relevant applications in agricultural sustainability, and we investigated the phylogenetic relationships of conserved and symbiotic genes of 21 bradyrhizobial strains. The study included strains from Western Australia (WA), isolated from nodules of Glycine spp. the country is one genetic center for the genus and from nodules of other indigenous legumes grown in WA, and strains isolated from forage Glycine sp. grown in South Africa. The 16S rRNA phylogeny divided the strains in two superclades, of B. japonicum and B. elkanii, but with low discrimination among the species. The multilocus sequence analysis (MLSA) with four protein-coding housekeeping genes (dnaK, glnII, gyrB and recA) pointed out seven groups as putative new species, two within the B. japonicum, and five within the B. elkanii superclades. The remaining eleven strains showed higher similarity with six species, B. lupini, B. liaoningense, B. yuanmingense, B. subterraneum, B. brasilense and B. retamae. Phylogenetic analysis of the nodC symbiotic gene clustered 13 strains in three different symbiovars (sv. vignae, sv. genistearum and sv. retamae), while seven others might compose new symbiovars. The genetic profiles of the strains evaluated by BOX-PCR revealed high intra- and interspecific diversity. The results point out the high level of diversity still to be explored within the Bradyrhizobium genus, and further studies might confirm new species and symbiovars.

Beneficial Endophytic Bacteria-Serendipita indica Interaction for Crop Enhancement and Resistance to Phytopathogens

Serendipita (=Piriformospora) indica is a fungal endophytic symbiont with the capabilities to enhance plant growth and confer resistance to different stresses. However, the application of this fungus in the field has led to inconsistent results, perhaps due to antagonism with other microbes. Here, we studied the impact of individual bacterial isolates from the endophytic bacterial community on the in vitro growth of S. indica. We further analyzed how combinations of bacteria and S. indica influence plant growth and protection against the phytopathogens Fusarium oxysporum and Rhizoctonia solani. Bacterial strains of the genera Bacillus, Enterobacter and Burkholderia negatively affected S. indica growth on plates, whereas Mycolicibacterium, Rhizobium, Paenibacillus strains and several other bacteria from different taxa stimulated fungal growth. To further explore the potential of bacteria positively interacting with S. indica, four of the most promising strains belonging to the genus Mycolicibacterium were selected for further experiments. Some dual inoculations of S. indica and Mycolicibacterium strains boosted the beneficial effects triggered by S. indica, further enhancing the growth of tomato plants, and alleviating the symptoms caused by the phytopathogens F. oxysporum and R. solani. However, some combinations of S. indica and bacteria were less effective than individual inoculations. By analyzing the genomes of the Mycolicibacterium strains, we revealed that these bacteria encode several genes predicted to be involved in the stimulation of S. indica growth, plant development and tolerance to abiotic and biotic stresses. Particularly, a high number of genes related to vitamin and nitrogen metabolism were detected. Taking into consideration multiple interactions on and inside plants, we showed in this study that some bacterial strains may induce beneficial effects on S. indica and could have an outstanding influence on the plant-fungus symbiosis.

Evolutionary insights into adaptation of Staphylococcus haemolyticus to human and non-human niches

Staphylococcus haemolyticus is a well-known member of human skin microbiome and an emerging opportunistic human pathogen. Presently, evolutionary studies are limited to human isolates even though it is reported from plants with beneficial properties and in environmental settings. In the present study, we report isolation of novel S. haemolyticus strains from surface sterilized rice seeds and compare their genome to other isolates from diverse niches available in public domain. The study showed expanding nature of pan-genome and revealed set of genes with putative functions related to its adaptability. This is seen by presence of type II lanthipeptide cluster in rice isolates, metal homeostasis genes in an isolate from copper coin and gene encoding methicillin resistance in human isolates. The present study on differential genome dynamics and role of horizontal gene transfers has provided novel insights into capability for ecological diversification of a bacterium of significance to human health.