VRprofile: gene-cluster-detection-based profiling of virulence and antibiotic resistance traits encoded within genome sequences of pathogenic bacteria

Overexpression of β-lactamase genes (blaKPC, blaSHV) and novel CirA deficiencies contribute to decreased cefiderocol susceptibility in carbapenem-resistant Klebsiella pneumoniae before its approval in China

Cefiderocol (FDC) is an effective antibiotic that is used to treat severe infections caused by carbapenem-resistant Klebsiella pneumoniae (CRKP). The mechanisms underlying FDC resistance and molecular epidemiology in China remain unclear. We collected 477 non-duplicate CRKP clinical isolates in central China and characterized their susceptibility to FDC, virulence genes, and sequence typing. The overall FDC susceptibility rate of CRKP was 99.2% in central China, which was higher than that in North America and Europe (96.1%), with MIC50/90 values of 1/2 mg/L. The decrease in FDC susceptibility in central China was concentrated in the ST11 CRKP-carrying virulence plasmids. Whole-genome sequencing (WGS) and quantitative reverse transcription PCR (qRT-PCR) experiments showed that serine β-lactamases, especially highly expressed KPC and SHV, substantially decreased FDC susceptibility in four FDC non-susceptible isolates (two resistant and two intermediate isolates). Notably, different CirA deficiencies, p.E450GfsTer16 and p.E133Ter, were found in both of the resistant isolates. In contrast, global WGS data indicate that the resistance mechanisms in North America and Europe were primarily associated with NDM and KPC variants, predominantly found in ST307 and ST147. Overall, FDC exhibits excellent activity against CRKP in central China, with resistance mechanisms primarily related to high KPC and SHV expression, along with deficiencies in CirA, frequently observed in ST11. This is remarkably different from the situation in North America and Europe and will directly impact the choice of clinical interventions. Additionally, the surveillance of FDC resistance in China is imperative.

Potential Convergence to Accommodate Pathogenicity Determinants and Antibiotic Resistance Revealed in Salmonella Mbandaka

Salmonella species are causal pathogens instrumental in human food-borne diseases. The pandemic survey related to multidrug resistant (MDR) Salmonella genomics enables the prevention and control of their dissemination. Currently, serotype Mbandaka is notorious as a multiple host-adapted non-typhoid Salmonella. However, its epidemic and MDR properties are still obscure, especially its genetic determinants accounting for virulence and MD resistance. Here, we aim to characterize the genetic features of a strain SMEH pertaining to Salmonella Mbandaka (S. Mbandaka), isolated from the patient's hydropericardium, using cell infections, a mouse model, antibiotic susceptibility test and comparative genomics. The antibiotic susceptibility testing showed that it could tolerate four antibiotics, including chloramphenicol, tetracycline, fisiopen and doxycycline by Kirby-Bauer (K-B) testing interpreted according to the Clinical and Laboratory Standards Institute (CLSI). Both the reproducibility in RAW 264.7 macrophages and invasion ability to infect HeLa cells with strain SMEH were higher than those of S. Typhimurium strain 14028S. In contrast, its attenuated virulence was determined in the survival assay using a mouse model. As a result, the candidate genetic determinants responsible for antimicrobial resistance, colonization/adaptability and their transferability were comparatively investigated, such as bacterial secretion systems and pathogenicity islands (SPI-1, SPI-2 and SPI-6). Moreover, collective efforts were made to reveal a potential role of the plasmid architectures in S. Mbandaka as the genetic reservoir to transfer or accommodate drug-resistance genes. Our findings highlight the essentiality of antibiotic resistance and risk assessment in S. Mbandaka. In addition, genomic surveillance is an efficient method to detect pathogens and monitor drug resistance. The genetic determinants accounting for virulence and antimicrobial resistance underscore the increasing clinical challenge of emerging MDR Mbandaka isolates, and provide insights into their prevention and treatment.

Comparative genomic analysis of Aeromonas dhakensis and Aeromonas hydrophila from diseased striped catfish fingerlings cultured in Vietnam

Introduction

Motile Aeromonas septicemia (MAS) is a burden for striped catfish (Pangasius hypophthalmus) farmers in Vietnam. MAS can be caused by several species of Aeromonas but Aeromonas hydrophila is seen as the leading cause of MAS in aquaculture, but recent reports suggest that A. dhakensis is also causing MAS.

Methods

Here we investigated the bacterial etiology of MAS and compared the genomic features of A. hydrophila and A. dhakensis. We collected 86 isolates from diseased striped catfish fingerlings over 5 years from eight provinces in Vietnam. Species identification was done using PCR, MALDI-TOF and whole genome sequence (WGS). The MICs of commonly used antimicrobials was established. Thirty presumed A. hydrophila isolates were sequenced for species confirmation and genomic comparison. A phylogenetic analysis was conducted using publicly available sequences and sequences from this study.

Results

A total of 25/30 isolates were A. dhakensis sequence type (ST) 656 and 5/30 isolates were A. hydrophila ST 251. Our isolates and all publicly available A. hydrophila isolates from Vietnam belonged to ST 251 and differed with <200 single nucleotide polymorphisms (SNP). Similarly, all A. dhakensis isolates from Vietnam belonged to ST 656 and differed with <100 SNPs. The tet(A) gene was found in 1/5 A. hydrophila and 19/25 A. dhakensis. All A. hydrophila had an MIC ≤2 mg/L while 19/25 A. dhakensis had MIC ≥8 mg/L for oxytetracycline. The floR gene was only found in A. dhakensis (14/25) which showed a MIC ≥8 mg/L for florfenicol. Key virulence genes, i.e., aerA/act, ahh1 and hlyA were present in all genomes, while ast was only present in A. dhakensis.

Discussion

This study confirms previous findings where A. dhakensis was the dominating pathogen causing MAS and that the importance of A. hydrophila has likely been overestimated. The differences in antimicrobial susceptibility between the two species could indicate a need for targeted antimicrobial treatment plans. The lipopolysaccharide regions and outer membrane proteins did not significantly differ in their immunogenic potentials, but it remains to be determined with in vivo experiments whether there is a difference in the efficacy of available vaccines against A. hydrophila and A. dhakensis.

Comprehensive Genomic Analysis Reveals Extensive Diversity of Type I and Type IV Secretion Systems in Klebsiella pneumoniae

The diversity and distribution of secretion systems in Klebsiella pneumoniae are unclear. In this study, the six common secretion systems (T1SS-T6SS) were comprehensively investigated in the genomes of 952 K. pneumoniae strains. T1SS, T2SS, type T subtype of T4SS, T5SS, and subtype T6SSi of T6SS were found. The findings indicated fewer types of secretion systems in K. pneumoniae than reported in Enterobacteriaceae, such as Escherichia coli. One conserved T2SS, one conserved T5SS, and two conserved T6SS were detected in more than 90% of the strains. In contrast, the strains displayed extensive diversity of T1SS and T4SS. Notably, T1SS and T4SS were enriched in the hypervirulent and classical multidrug resistance pathotypes of K. pneumoniae, respectively. The results expand the epidemiological knowledge of the virulence and transmissibility of pathogenic K. pneumoniae and contribute to identify the potential strains for safe applications.

Identification of Salmonella Pullorum Factors Affecting Immune Reaction in Macrophages from the Avian Host

The host-specific Salmonella serovar S. Pullorum (SP) modulates the chicken immune response to a Th2-biased response associated with persistent infection. This is different from the Th1-biased immune response induced by the genetically close serovar, S. Enteritidis (SE). Based on core genome differences between SP and SE, we used three complementary bioinformatics approaches to identify SP genes, which may be important for stimulation of the immune response. Defined mutants were constructed in selected genes, and the infection potential and ability of mutants to stimulate cytokine production in avian derived HD11 macrophages were determined. Deletion of large genomic regions unique to SP did not change infection potential nor immune stimulation significantly. Mutants in genes with conserved single nucleotide polymorphisms (SNPs) between the two serovars in the region 100 bp upstream of the start codon (conserved upstream SNPs [CuSNPs]) such as sseE, osmB, tolQ, a putative immune antigen, and a putative persistent infection factor, exhibited differences in induction of inflammatory cytokines compared to wild-type SP, suggesting a possible role of these CuSNPs in immune regulation. Single nucleotide SP mutants correcting for the CuSNP difference were constructed in the upstream region of sifA and pipA. The SNP corrected pipA mutant expressed pipA at a higher level than the wild-type SP strain, and the mutant differentially caused upregulation of proinflammatory cytokines. It suggests that this CuSNP is important for the suppression of proinflammatory responses. In conclusion, this study has identified putative immune stimulating factors of relevance to the difference in infection dynamics between SP and SE in avian macrophages. IMPORTANCE Salmonella Pullorum is host specific to avian species, where it causes life-threatening infection in young birds. It is unknown why it is host restricted and causes systemic disease, rather than gastroenteritis normally seen with Salmonella. In the present study, we identified genes and single nucleotide polymorphisms (SNPs; relative to the broad-host-range type Salmonella Enteritidis), which affected survival and immune induction in macrophages from hens suggesting a role in development of the host specific infection. Further studies of such genes may enable understanding of which genetic factors determine the development of host specific infection by S. Pullorum. In this study, we developed an in silico approach to predict candidate genes and SNPs for development of the host-specific infection and the specific induction of immunity associated with this infection. This study flow can be used in similar studies in other clades of bacteria.

Genomic insights into antibiotic resistance and mobilome of lactic acid bacteria and bifidobacteria

Lactic acid bacteria (LAB) and Bifidobacterium sp. (bifidobacteria) can carry antimicrobial resistance genes (ARGs), yet data on resistance mechanisms in these bacteria are limited. The aim of our study was to identify the underlying genetic mechanisms of phenotypic resistance in 103 LAB and bifidobacteria using whole-genome sequencing. Sequencing data not only confirmed the presence of 36 acquired ARGs in genomes of 18 strains, but also revealed wide dissemination of intrinsic ARGs. The presence of acquired ARGs on known and novel mobile genetic elements raises the possibility of their horizontal spread. In addition, our data suggest that mutations may be a common mechanism of resistance. Several novel candidate resistance mechanisms were uncovered, providing a basis for further in vitro studies. Overall, 1,314 minimum inhibitory concentrations matched with genotypes in 92.4% of the cases; however, prediction of phenotype based on genotypic data was only partially efficient, especially with respect to aminoglycosides and chloramphenicol. Our study sheds light on resistance mechanisms and their transferability potential in LAB and bifidobacteria, which will be useful for risk assessment analysis.

Genomic analysis of a new type VI secretion system in Vibrio parahaemolyticus and its implications for environmental adaptation in shrimp ponds

The type VI secretion system (T6SS) in Vibrio spp. is often used to kill heteroclonal neighbors by direct injection of toxic effectors, but its strategies in aquacultural environments receive limited attention. In this study, we conducted genomic analysis for a T6SS-harboring plasmid in V. parahaemolyticus strain VP157. Coculture assays were further conducted to verify its antibacterial function. The results showed that strain VP157 harbored a 132-kb plasmid, pVP157-1, which consists of two fragments: an 87.8-kb fragment identical to plasmid pTJ114-1 and a 44.2-kb T6SS gene cluster with only 4% DNA identity to T6SS1 in the V. parahaemolyticus reference genome. Gene-by-gene analysis of six genes representing core T6SS components suggested that each gene has distinct evolutionary origins. In vitro experimental evolution revealed that pVP157-1 can excise from the VP157 genome with an excision rate of 4%. A coculture assay suggested that strain VP157 had significantly higher antibacterial activity against Bacillus pumilus and V. cholerae than the strain without pVP157-1(VP157^∆T6SS). In contrast, a rapid decline was observed for the proportion of VP157^{∆ T6SS} in a mock microbial community, which decreased from 10.7% to 2.1% in 5 days. The results highlighted that the acquisition of T6SS fostered the fitness of V . parahaemolyticus in a complex environment.

Comparative genome analysis of the genus Shewanella unravels the association of key genetic traits with known and potential pathogenic lineages

Shewanella spp. are Gram-negative rods widely disseminated in aquatic niches that can also be found in human-associated environments. In recent years, reports of infections caused by these bacteria have increased significantly. Mobilome and resistome analysis of a few species showed that they are versatile; however, comprehensive comparative studies in the genus are lacking. Here, we analyzed the genetic traits of 144 genomes from Shewanella spp. isolates focusing on the mobilome, resistome, and virulome to establish their evolutionary relationship and detect unique features based on their genome content and habitat. Shewanella spp. showed a great diversity of mobile genetic elements (MGEs), most of them associated with monophyletic lineages of clinical isolates. Furthermore, 79/144 genomes encoded at least one antimicrobial resistant gene with their highest occurrence in clinical-related lineages. CRISPR-Cas systems, which confer immunity against MGEs, were found in 41 genomes being I-E and I-F the more frequent ones. Virulome analysis showed that all Shewanella spp. encoded different virulence genes (motility, quorum sensing, biofilm, adherence, etc.) that may confer adaptive advantages for survival against hosts. Our data revealed that key accessory genes are frequently found in two major clinical-related groups, which encompass the opportunistic pathogens Shewanella algae and Shewanella xiamenensis together with several other species. This work highlights the evolutionary nature of Shewanella spp. genomes, capable of acquiring different key genetic traits that contribute to their adaptation to different niches and facilitate the emergence of more resistant and virulent isolates that impact directly on human and animal health.

Whole genome sequence analysis of Ligilactobacillus agilis C7 isolated from pig feces revealed three bacteriocin gene clusters

We here report the whole genome sequence of Ligilactobacillus agilis C7 with anti-listerial activity, which was isolated from pig feces. The genome size of L. agilis C7 (~ 3.0 Mb) is relatively larger compared with other L. agilis strains. L. agilis C7 carries three bacteriocin gene clusters encoding garvicin Q, salivaricin A, and Blp family class II bacteriocin. Garvicin Q and salivaricin A are reported to be active against Listeria monocytogenes and Micrococcus luteus, respectively, as well as against other Gram-positive bacteria. Meanwhile, the bacteriocin encoded in the blp cassette was shown to be active against pneumococci, mediating intraspecies competition. This report highlights the potential of L. agilis C7 for the production of bacteriocins inhibiting pathogenic bacteria.

Transmission of Nonconjugative Virulence or Resistance Plasmids Mediated by a Self-Transferable IncN3 Plasmid from Carbapenem-Resistant Klebsiella pneumoniae

Klebsiella pneumoniae poses a critical challenge to clinical and public health. Along with conjugative plasmids, nonconjugative resistance or virulence plasmids associated with carbapenem-resistant K. pneumoniae (CRKP), hypervirulent K. pneumoniae (hvKP), and even carbapenem-resistant and hypervirulent K. pneumoniae (CR-hvKP) strains have been spreading globally. In this study, a clinical CRKP strain KP2648 was isolated, and the transferability of its plasmids was assessed using conjugation experiments. The transconjugants were characterized by polymerase chain reaction (PCR) detection, XbaI and S1-pulsed-field gel electrophoresis (PFGE), and/or whole-genome sequencing. Genetically modified IncN3 plasmids were employed to elucidate the self-transferability and the mobilization mechanisms. KP2648 has three natural plasmids: a nonconjugative IncFIB/IncHI3B virulence plasmid, a nonconjugative IncFII/IncR carbapenem-resistant plasmid, and a self-transferable IncN3 plasmid with a high conjugation frequency (7.54 ± 1.06) × 10-1. The IncN3 plasmid could mobilize the coexisting nonconjugative virulence/resistance plasmids either directly or by employing intermediate E. coli with two forms: a hybrid plasmid fused with IncN3 or a cotransfer with the helper plasmid, IncN3. Various mobile genetic elements, including ISKpn74, ISKpn14, IS26, ISShes11, ISAba11, and Tn3, are involved in the genetic transposition of diverse hybrid plasmids and the cotransfer process during the intra/interspecies transmission. IMPORTANCE Nowadays, the underlying mobilization mechanism and evolutionary processes of nonconjugative virulence or resistance plasmids in Klebsiella pneumoniae remain poorly understood. Our study revealed the high conjugation ability of IncN3 plasmid isolated from carbapenem-resistant K. pneumoniae and confirmed its capability to mobilize the nonconjugative virulence or resistance plasmids. The self-transferable IncN3 plasmid could facilitate the transmission of pathogenicity and genetic evolution of carbapenem-resistant and hypervirulent K. pneumoniae, including hv-CRKP (virulence plasmid obtained by carbapenem-resistant K. pneumoniae) and CR-hvKP (resistance plasmid obtained by hypervirulent K. pneumoniae), warranting further monitoring.

The Notable Achievements and the Prospects of Bacterial Pathogen Genomics

Throughout the entirety of human history, bacterial pathogens have played an important role and even shaped the fate of civilizations. The application of genomics within the last 27 years has radically changed the way we understand the biology and evolution of these pathogens. In this review, we discuss how the short- (Illumina) and long-read (PacBio, Oxford Nanopore) sequencing technologies have shaped the discipline of bacterial pathogen genomics, in terms of fundamental research (i.e., evolution of pathogenicity), forensics, food safety, and routine clinical microbiology. We have mined and discuss some of the most prominent data/bioinformatics resources such as NCBI pathogens, PATRIC, and Pathogenwatch. Based on this mining, we present some of the most popular sequencing technologies, hybrid approaches, assemblers, and annotation pipelines. A small number of bacterial pathogens are of very high importance, and we also present the wealth of the genomic data for these species (i.e., which ones they are, the number of antimicrobial resistance genes per genome, the number of virulence factors). Finally, we discuss how this discipline will probably be transformed in the near future, especially by transitioning into metagenome-assembled genomes (MAGs), thanks to long-read sequencing.

Salmonella Salamae and S. Waycross isolated from Nile perch in Lake Victoria show limited human pathogenic potential

Non-enterica subspecies of Salmonella enterica are rarely associated with human infections. Paradoxically, food safety legislations consider the entire genus Salmonella as pathogenic to humans. Globally, large amounts of seafoods are rejected and wasted due to findings of Salmonella. To inform better food safety decisions, we investigated the pathogenicity of Salmonella Salamae 42:r- and Salmonella Waycross isolated from Nile perch from Lake Victoria. Genome-wide analysis revealed absence of significant virulence determinants including on key Salmonella pathogenicity islands in both serovars. In epithelial cells, S. Salamae showed a weak invasion ability that was lower than the invH mutant of S. Typhimiurium used as negative control. Similarly, S. Salamae could not replicate inside macrophages. Moreover, intracellular replication in S. Waycross strains was significantly lower compared to the wild type S. Typhimurium. Our findings suggest a low pathogenicity of S. Salamae reinforcing the existing literature that non-enterica subspecies are avirulent. We propose that food legislations and actions taken on findings of Salmonella are revisited to avoid wasting valuable sea- and other foods.

A review: antimicrobial resistance data mining models and prediction methods study for pathogenic bacteria

Antimicrobials have paved the way for medical and social development over the last century and are indispensable for treating infections in humans and animals. The dramatic spread and diversity of antibiotic-resistant pathogens have significantly reduced the efficacy of essentially all antibiotic classes and is a global problem affecting human and animal health. Antimicrobial resistance is influenced by complex factors such as resistance genes and dosing, which are highly nonlinear, time-lagged and multivariate coupled, and the amount of resistance data is large and redundant, making it difficult to predict and analyze. Based on machine learning methods and data mining techniques, this paper reviews (1) antimicrobial resistance data storage and analysis techniques, (2) antimicrobial resistance assessment methods and the associated risk assessment methods for antimicrobial resistance, and (3) antimicrobial resistance prediction methods. Finally, the current research results on antimicrobial resistance and the development trend are summarized to provide a systematic and comprehensive reference for the research on antimicrobial resistance.

Colistin and Carbapenem-Resistant Acinetobacter baumannii Aci46 in Thailand: Genome Analysis and Antibiotic Resistance Profiling

Resistance to the last-line antibiotics against invasive Gram-negative bacterial infection is a rising concern in public health. Multidrug resistant (MDR) Acinetobacter baumannii Aci46 can resist colistin and carbapenems with a minimum inhibitory concentration of 512 µg/mL as determined by microdilution method and shows no zone of inhibition by disk diffusion method. These phenotypic characteristics prompted us to further investigate the genotypic characteristics of Aci46. Next generation sequencing was applied in this study to obtain whole genome data. We determined that Aci46 belongs to Pasture ST2 and is phylogenetically clustered with international clone (IC) II as the predominant strain in Thailand. Interestingly, Aci46 is identical to Oxford ST1962 that previously has never been isolated in Thailand. Two plasmids were identified (pAci46a and pAci46b), neither of which harbors any antibiotic resistance genes but pAci46a carries a conjugational system (type 4 secretion system or T4SS). Comparative genomics with other polymyxin and carbapenem-resistant A. baumannii strains (AC30 and R14) identified shared features such as CzcCBA, encoding a cobalt/zinc/cadmium efflux RND transporter, as well as a drug transporter with a possible role in colistin and/or carbapenem resistance in A. baumannii. Single nucleotide polymorphism (SNP) analyses against MDR ACICU strain showed three novel mutations i.e., Glu229Asp, Pro200Leu, and Ala138Thr, in the polymyxin resistance component, PmrB. Overall, this study focused on Aci46 whole genome data analysis, its correlation with antibiotic resistance phenotypes, and the presence of potential virulence associated factors.

Mobilization of the nonconjugative virulence plasmid from hypervirulent Klebsiella pneumoniae

Background

Klebsiella pneumoniae, as a global priority pathogen, is well known for its capability of acquiring mobile genetic elements that carry resistance and/or virulence genes. Its virulence plasmid, previously deemed nonconjugative and restricted within hypervirulent K. pneumoniae (hvKP), has disseminated into classic K. pneumoniae (cKP), particularly carbapenem-resistant K. pneumoniae (CRKP), which poses alarming challenges to public health. However, the mechanism underlying its transfer from hvKP to CRKP is unclear.

Methods

A total of 28 sequence type (ST) 11 bloodstream infection-causing CRKP strains were collected from Ruijin Hospital in Shanghai, China, and used as recipients in conjugation assays. Transconjugants obtained from conjugation assays were confirmed by XbaI and S1 nuclease pulsed-field gel electrophoresis, PCR detection and/or whole-genome sequencing. The plasmid stability of the transconjugants was evaluated by serial culture. Genetically modified strains and constructed mimic virulence plasmids were employed to investigate the mechanisms underlying mobilization. The level of extracellular polysaccharides was measured by mucoviscosity assays and uronic acid quantification. An in silico analysis of 2608 plasmids derived from 814 completely sequenced K. pneumoniae strains available in GenBank was performed to investigate the distribution of putative helper plasmids and mobilizable virulence plasmids.

Results

A nonconjugative virulence plasmid was mobilized by the conjugative plasmid belonging to incompatibility group F (IncF) from the hvKP strain into ST11 CRKP strains under low extracellular polysaccharide-producing conditions or by employing intermediate E. coli strains. The virulence plasmid was mobilized via four modes: transfer alone, cotransfer with the conjugative IncF plasmid, hybrid plasmid formation due to two rounds of single-strand exchanges at specific 28-bp fusion sites or homologous recombination. According to the in silico analysis, 31.8% (242) of the putative helper plasmids and 98.8% (84/85) of the virulence plasmids carry the 28-bp fusion site. All virulence plasmids carry the origin of the transfer site.

Conclusions

The nonconjugative virulence plasmid in ST11 CRKP strains is putatively mobilized from hvKP or E. coli intermediates with the help of conjugative IncF plasmids. Our findings emphasize the importance of raising public awareness of the rapid dissemination of virulence plasmids and the consistent emergence of hypervirulent carbapenem-resistant K. pneumoniae (hv-CRKP) strains.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00936-5.

The global emergence of a novel Streptococcus suis clade associated with human infections

Streptococcus suis, a ubiquitous bacterial colonizer in pigs, has recently extended host range to humans, leading to a global surge of deadly human infections and three large outbreaks since 1998. To better understand the mechanisms for the emergence of cross-species transmission and virulence in human, we have sequenced 366 S. suis human and pig isolates from 2005 to 2016 and performed a large-scale phylogenomic analysis on 1,634 isolates from 14 countries over 36 years. We show the formation of a novel human-associated clade (HAC) diversified from swine S. suis isolates. Phylogeographic analysis identified Europe as the origin of HAC, coinciding with the exportation of European swine breeds between 1960s and 1970s. HAC is composed of three sub-lineages and contains several healthy-pig isolates that display high virulence in experimental infections, suggesting healthy-pig carriers as a potential source for human infection. New HAC-specific genes are identified as promising markers for pathogen detection and surveillance. Our discovery of a human-associated S. suis clade provides insights into the evolution of this emerging human pathogen and extend our understanding of S. suis epidemics worldwide.

Pseudoprevotella muciniphila gen. nov., sp. nov., a mucin-degrading bacterium attached to the bovine rumen epithelium

A Gram-negative, strictly anaerobic mucin-degrading bacterium, which we designated strain E39^T, was isolated from the rumen epithelium of Korean cattle. The cells were non-motile and had a coccus morphology. Growth of strain E39^T was observed at 30–45°C (optimum, 39°C), pH 6.5–8.5 (optimum, pH 7.5), and in the presence of 0.0–1.0% (w/v) NaCl (optimum, 0.0–0.5%). Strain E39^T contained C_16:0, C_18:0, C_18:1ω9c, iso-C_15:0, and anteiso-C_15:0 as the major fatty acids. The major polar lipids were phosphatidylethanolamine, unidentified aminophospholipid, and unidentified lipids. The major respiratory isoprenoid quinones were MK-8 and MK-9. The major fermented end-products of mucin were acetate and succinate. The G+C content of the genomic DNA was 46.4 mol%. Strain E39^T was most closely related to Alloprevotella rava 81/4-12^T with an 87.3% 16S rRNA gene sequence similarity. On the basis of phenotypic, chemotaxonomic, and molecular properties, strain E39^T represents a novel genus of the family Prevotellaceae; as such, the name Pseudoprevotella muciniphila gen. nov., sp. nov. is proposed. A functional annotation of the whole genome sequences of P. muciniphila E39^T revealed that this bacterium has a putative mucin-degrading pathway and biosynthetic pathways of extracellular polymeric substances and virulence factors which enable bacteria to adhere to the epithelial cells and avoid the host’s immune responses.

Bacteriophages: from isolation to application

Bacteriophages are considered as a potential alternative to fight pathogenic bacteria during the antibiotic resistance era. With their high specificity, they are being widely used in various applications: medicine, food industry, agriculture, animal farms, biotechnology, diagnosis, etc. Many techniques have been designed by different researchers for phage isolation, purification, and amplification, each of which has strengths and weaknesses. However, all aim at having a reasonably pure phage sample that can be further characterized. Phages can be characterized based on their physiological, morphological or inactivation tests. Microscopy, in particular, has opened a wide gate not only for visualizing phage morphological structure, but also for monitoring biochemistry and behavior. Meanwhile, computational analysis of phage genomes provides more details about phage history, lifestyle, and potential for toxigenic or lysogenic conversion, which translate to safety in biocontrol and phage therapy applications. This review summarizes phage application pipelines at different levels and addresses specific restrictions and knowledge gaps in the field. Recently developed computational approaches, which are used in phage genome analysis, are critically assessed. We hope that this assessment provides researchers with useful insights for selection of suitable approaches for Phage-related research aims and applications.

Genome analysis provides insights into the biocontrol ability of Mitsuaria sp. strain TWR114

Mitsuaria sp. TWR114 is a biocontrol agent against tomato bacterial wilt (TBW). We aimed to gain genomic insights relevant to the biocontrol mechanisms and colonization ability of this strain. The draft genome size was found to be 5,632,523 bp, with a GC content of 69.5%, assembled into 1144 scaffolds. Genome annotation predicted a total of 4675 protein coding sequences (CDSs), 914 pseudogenes, 49 transfer RNAs, 3 noncoding RNAs, and 2 ribosomal RNAs. Genome analysis identified multiple CDSs associated with various pathways for the metabolism and transport of amino acids and carbohydrates, motility and chemotactic capacities, protection against stresses (oxidative, antibiotic, and phage), production of secondary metabolites, peptidases, quorum-quenching enzymes, and indole-3-acetic acid, as well as protein secretion systems and their related appendages. The genome resource will extend our understanding of the genomic features related to TWR114's biocontrol and colonization abilities and facilitate its development as a new biopesticide against TBW.

Microbial Genomics as a Catalyst for Targeted Antivirulence Therapeutics

Virulence arresting drugs (VAD) are an expanding class of antimicrobial treatment that act to “disarm” rather than kill bacteria. Despite an increasing number of VAD being registered for clinical use, uptake is hampered by the lack of methods that can identify patients who are most likely to benefit from these new agents. The application of pathogen genomics can facilitate the rational utilization of advanced therapeutics for infectious diseases. The development of genomic assessment of VAD targets is essential to support the early stages of VAD diffusion into infectious disease management. Genomic identification and characterization of VAD targets in clinical isolates can augment antimicrobial stewardship and pharmacovigilance. Personalized genomics guided use of VAD will provide crucial policy guidance to regulating agencies, assist hospitals to optimize the use of these expensive medicines and create market opportunities for biotech companies and diagnostic laboratories.

A Type VI Secretion System Facilitates Fitness, Homeostasis, and Competitive Advantages for Environmental Adaptability and Efficient Nicotine Biodegradation

Gram-negative bacteria employ secretion systems to translocate proteinaceous effectors from the cytoplasm to the extracellular milieu, thus interacting with the surrounding environment or microniche. It is known that bacteria can benefit from the type VI secretion system (T6SS) by transporting ions to combat reactive oxygen species (ROS). Here, we report that T6SS activities conferred tolerance to nicotine-induced oxidative stress in Pseudomonas sp. strain JY-Q, a highly active nicotine degradation strain isolated from tobacco waste extract. AA098_13375 was identified to encode a dual-functional effector with antimicrobial and anti-ROS activities. Wild-type strain JY-Q grew better than the AA098_13375 deletion mutant in nicotine-containing medium by antagonizing increased intracellular ROS levels. It was, therefore, tentatively designated TseN (type VI secretion system effector for nicotine tolerance), homologs of which were observed to be broadly ubiquitous in Pseudomonas species. TseN was identified as a Tse6-like bacteriostatic toxin via monitoring intracellular NAD+ TseN presented potential antagonism against ROS to fine tune the heavy traffic of nicotine metabolism in strain JY-Q. It is feasible that the dynamic tuning of NAD+ driven by TseN could satisfy demands from nicotine degradation with less cytotoxicity. In this scenario, T6SS involves a fascinating accommodation cascade that prompts constitutive biotransformation of N-heterocyclic aromatics by improving bacterial robustness/growth. In summary, the T6SS in JY-Q mediated resistance to oxidative stress and promoted bacterial fitness via a contact-independent growth competitive advantage, in addition to the well-studied T6SS-dependent antimicrobial activities.IMPORTANCE Mixtures of various pollutants and the coexistence of numerous species of organisms are usually found in adverse environments. Concerning biodegradation of nitrogen-heterocyclic contaminants, the scientific community has commonly focused on screening functional enzymes that transform pollutants into intermediates of attenuated toxicity or for primary metabolism. Here, we identified dual roles of the T6SS effector TseN in Pseudomonas sp. strain JY-Q, which is capable of degrading nicotine. The T6SS in strain JY-Q is able to deliver TseN to kill competitors and provide a growth advantage by a contact-independent pattern. TseN could monitor the intracellular NAD+ level by its hydrolase activity, causing cytotoxicity in competitive rivals but metabolic homeostasis on JY-Q. Moreover, JY-Q could be protected from TseN toxicity by the immunity protein TsiN. In conclusion, we found that TseN with cytotoxicity to bacterial competitors facilitated the nicotine tolerance of JY-Q. We therefore reveal a working model between T6SS and nicotine metabolism. This finding indicates that multiple diversified weapons have been evolved by bacteria for their growth and robustness.

Hidden Markov Model: a shortest unique representative approach to detect the protein toxins, virulence factors and antibiotic resistance genes

Objective

Currently, next generation sequencing (NGS) is widely used to decode potential novel or variant pathogens both in emergent outbreaks and in routine clinical practice. However, the efficient identification of novel or diverged pathogenomic compositions remains a big challenge. It is especially true for short DNA sequence fragments from NGS, since sequence similarity searching is vulnerable to false negatives or false positives, as is mismatching or matching with unrelated proteins. Therefore, this study aimed to establish a bioinformatics approach that can generate unique motif sequences for profiling searching, resulting in high specificity and sensitivity.

Results

In this study, we introduced a Shortest Unique Representative Hidden Markov Model (HMM) approach to identify bacterial toxin, virulence factor (VF), and antimicrobial resistance (AR) in short sequence reads. We first construct unique representative domain sequences of toxin genes, VFs, and ARs to avoid potential false positives, and then to use HMM models to accurately identify potential toxin, VF, and AR fragments. The benchmark shows this approach can achieve relatively high specificity and sensitivity if the appropriate cutoff value is applied. Our approach can be used to recognize the protein sequences of known toxins and pathogens, identifies their common characteristics and then searches for similar sequences in other organisms.

Additional Role of Nicotinic Acid Hydroxylase for the Transformation of 3-Succinoyl-Pyridine by Pseudomonas sp. Strain JY-Q

Nicotine and nicotinic acid (NA) are both considered to be representatives of N-heterocyclic aromatic compounds, and their degradation pathways have been revealed in Pseudomonas species. However, the cooccurrence of these two pathways has only been observed in Pseudomonas sp. strain JY-Q. The nicotine pyrrolidine catabolism pathway of strain JY-Q consists of the functional modules Nic1, Spm, and Nic2. The module enzyme, 3-succinoylpyridine monooxygenase (Spm), catalyzes transformation of 3-succinoyl-pyridine (SP) to 6-hydroxy-3-succinoyl-pyridine (HSP). There exist two homologous but not identical Spm enzymes (namely, Spm1 and Spm2) in JY-Q. However, when spm1 and spm2 were both in-frame deleted, the mutant still grew well in basic salt medium (BSM) supplemented with nicotine as the sole carbon/nitrogen nutrition, suggesting that there exists an alternative pathway responsible for SP catabolism in JY-Q. NicAB, an enzyme accounting for NA hydroxylation, contains reorganized domains similar to those of Spm. When the JY-Q_nicAB gene (nicAB in strain JY-Q) was introduced into another Pseudomonas strain, one that is unable to degrade NA, the resultant recombinant strain exhibited the ability to transform SP to HSP, but without the ability to metabolize NA. Here, we conclude that NicAB in strain JY-Q exhibits an additional role in SP transformation. The other genes in the NA cluster, NicXDFE (Nic2 homolog), then also exhibit a role in subsequent HSP metabolism for energy yield. This finding also suggests that the cooccurrence of nicotine and NA degradation genes in strain JY-Q represents an advantage for JY-Q, making it more effective and flexible for the degradation of nicotine.IMPORTANCE 3-Succinoyl-pyridine (SP) and 6-hydroxy-3-succinoyl-pyridine (HSP) are both valuable chemical precursors to produce insecticides and hypotensive agents. SP and HSP could be renewable through the nicotine microbial degradation pathway, in which 3-succinoylpyridine monooxygenases (Spm) account for transforming SP into HSP in Pseudomonas sp. strain JY-Q. However, when two homologous Spm genes (spm1 and spm2) were knocked out, the mutant retained the ability to degrade nicotine. Thus, in addition to Spm, JY-Q should have an alternative pathway for SP conversion. In this research, we showed that JY-Q_NicAB was responsible for this alternative SP conversion. Both of the primary functions for nicotinic acid dehydrogenation and the additional function for SP metabolism were detected in a recombinant strain harboring JY-Q_NicAB. As a result, both nicotinic acid and nicotine degradation pathways in JY-Q contribute to its remarkable nicotine tolerance and nicotine degradation availability. These findings also provide one more metabolic engineering strategy for accumulation for value-added intermediates.

Comparative genomics reveals broad genetic diversity, extensive recombination and nascent ecological adaptation in Micrococcus luteus

Background

Micrococcus luteus is a group of actinobacteria that is widely used in biotechnology and is being thought as an emerging nosocomial pathogen. With one of the smallest genomes of free-living actinobacteria, it is found in a wide range of environments, but intraspecies genetic diversity and adaptation strategies to various environments remain unclear. Here, comparative genomics, phylogenomics, and genome-wide association studies were used to investigate the genomic diversity, evolutionary history, and the potential ecological differentiation of the species.

Results

High-quality genomes of 66 M. luteus strains were downloaded from the NCBI GenBank database and core and pan-genome analysis revealed a considerable intraspecies heterogeneity. Phylogenomic analysis, gene content comparison, and average nucleotide identity calculation consistently indicated that the species has diverged into three well-differentiated clades. Population structure analysis further suggested the existence of an unknown ancestor or the fourth, yet unsampled, clade. Reconstruction of gene gain/loss events along the evolutionary history revealed both early events that contributed to the inter-clade divergence and recent events leading to the intra-clade diversity. We also found convincing evidence that recombination has played a key role in the evolutionary process of the species, with upto two-thirds of the core genes having been affected by recombination. Furthermore, distribution of mammal-associated strains (including pathogens) on the phylogenetic tree suggested that the last common ancestor had a free-living lifestyle, and a few recently diverged lineages have developed a mammal-associated lifestyle separately. Consistently, genome-wide association analysis revealed that mammal-associated strains from different lineages shared genes functionally relevant to the host-associated lifestyle, indicating a recent ecological adaption to the new host-associated habitats.

Conclusions

These results revealed high intraspecies genomic diversity of M. luteus and highlighted that gene gain/loss events and extensive recombination events played key roles in the genome evolution. Our study also indicated that, as a free-living species, some lineages have recently developed or are developing a mammal-associated lifestyle. This study provides insights into the mechanisms that drive the genome evolution and adaption to various environments of a bacterial species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07432-5.

Higher genome variability within metabolism genes associates with recurrent Clostridium difficile infection

BACKGROUND

Clostridium difficile (C. difficile) is a major source of healthcare-associated infection with a high risk of recurrence, attributable to many factors such as usage of antibiotics, older age and immunocompromised status of the patients. C. difficile has also a highly diverse genome, which may contribute to its high virulence. Herein we examined whether the genome conservation, measured as non-synonymous to synonymous mutations ratio (dN/dS) in core genes, presence of single genes, plasmids and prophages increased the risk of reinfection in a subset of 134 C. difficile isolates from our previous study in a singly hemato-oncology ward.

METHODS

C. difficile isolates were subjected to whole-genome sequencing (WGS) on Ion Torrent PGM sequencer. Genomes were assembled with MIRA5 and annotated with prokka and VRprofile. Logistic regression was used to asses the relationship between single gene presence and the odds of infection recurrence. DN/dS ratios were computed with codeml. Functional annotation was conducted with eggNOG-Mapper.

RESULTS

We have found that the presence of certain genes, associated with carbon metabolism and oxidative phosphorylation, increased the odds of infection recurrence. More core genes were under positive selective pressure in recurrent disease isolates - they were mostly associated with the metabolism of aminoacids. Finally, prophage elements were more prevalent in single infection isolates and plasmids did not influence the odds of recurrence.

CONCLUSIONS

Our findings suggest higher genetic plasticity in isolates causing recurrent infection, associated mainly with metabolism. On the other hand, the presence of prophages seems to reduce the isolates' virulence.

Dynamics and Microevolution of Vibrio parahaemolyticus Populations in Shellfish Farms

Globally, V. parahaemolyticus-related gastroenteritis outbreaks caused by seafood consumption represent an increasing threat to human health. Despite advances in our understanding of the global epidemiology of pandemic V. parahaemolyticus, fundamental questions about the key driving forces for the spread of V. parahaemolyticus at regional and national scales remain unanswered.

Vibrio parahaemolyticus is becoming the leading cause of acute bacterial gastroenteritis, but its population dynamics in aquafarms have received limited attention. To address this research gap, we selected three shellfish farms to examine the impacts of ocean currents and the transport of live aquatic animals on the transmission and microevolution of V. parahaemolyticus by using multilocus sequence typing (MLST) and whole-genome sequencing. MLST and genomic analysis revealed that the community structure of V. parahaemolyticus in Dalian and Donggang was relatively stable in the presence of ocean currents; however, horizontal gene transfer of mobile genetic elements (MGEs) between Dalian and Donggang was very common. Further analysis indicated that the transport of live aquatic animals from Dalian to Xiamen not only introduced new V. parahaemolyticus populations but also allowed the exchange of genetic material between the two sites. More interestingly, Dalian-originated strain ST722 was introduced to Xiamen farms, resulting in one MLST allele change and the acquisition of two genomic islands from indigenous isolates in Xiamen within 8 months; such alterations are thought to promote the adaptation of V. parahaemolyticus. These results provide direct observations of how ocean currents and the transport of live aquatic animals contribute to the dissemination and genetic mixture of V. parahaemolyticus, which provides insights into the dynamics and microevolution of V. parahaemolyticus in aquacultural environments.

IMPORTANCE Globally, V. parahaemolyticus-related gastroenteritis outbreaks caused by seafood consumption represent an increasing threat to human health. Despite advances in our understanding of the global epidemiology of pandemic V. parahaemolyticus, fundamental questions about the key driving forces for the spread of V. parahaemolyticus at regional and national scales remain unanswered. This study revealed that the transregional transport of aquatic animals and the movement of ocean currents both contributed to the mixing of V. parahaemolyticus populations. More importantly, this study demonstrated how genetic mixture occurred between introduced and endemic V. parahaemolyticus populations via the transport of aquatic animals, which accelerated bacterial adaptation by transferring ecologically important functions. These results suggest that human activities entail a risk of the emergence of new virulent populations for both aquatic animals and humans by horizontal gene transfer and provide important insights into the microevolution and population mixing of V. parahaemolyticus.

Closed Genome and Plasmid Sequences of Legionella pneumophila AW-13-4, Isolated from a Hot Water Loop System of a Large Occupational Building

Unused water in unoccupied buildings can become stagnant, with reductions in temperature and levels of disinfectant resulting in increased microbial growth. We report the closed and complete genome and plasmid of Legionella pneumophila strain AW-13-4 (serogroup 1), which was isolated from a hot water loop system of a large building.

Analysis of the core genome and pangenome of Clostridium butyricum

Clostridium butyricum is an anaerobic bacterium that inhabits broad niches. Clostridium butyricum is known for its production of butyrate, 1,3-propanediol, and hydrogen. This study aimed to present a comparative pangenome analysis of 24 strains isolated from different niches. We sequenced and annotated the genome of C. butyricum 3-3 isolated from the Chinese baijiu ecosystem. The pangenome of C. butyricum was open. The core genome, accessory genome, and strain-specific genes comprised 1011, 4543, and 1473 genes, respectively. In the core genome, Carbohydrate metabolism was the largest category, and genes in the biosynthetic pathway of butyrate and glycerol metabolism were conserved (in the core or soft-core genome). Furthermore, the 1,3-propanediol operon existed in 20 strains. In the accessory genome, numerous mobile genetic elements belonging to the Replication, recombination, and repair (L) category were identified. In addition, genome islands were identified in all 24 strains, ranging from 2 (strain KNU-L09) to 53 (strain SU1), and phage sequences were found in 17 of the 24 strains. This study provides an important genomic framework that could pave the way for the exploration of C. butyricum and future studies on the genetic diversification of C. butyricum.

Forensic genomics of a novel Klebsiella quasipneumoniae type from a neonatal intensive care unit in China reveals patterns of colonization, evolution and epidemiology

During March 2017, a neonatal patient with severe diarrhoea subsequently developed septicaemia and died, with Klebsiella isolated as the causative microorganism. In keeping with infection control protocols, the coincident illness of an attending staff member and three other neonates with Klebsiella infection triggered an outbreak response, leading to microbiological assessment of isolates collected from the staff member and all 21 co-housed neonates. Multilocus sequence typing and genomic sequencing identified that the isolates from the 21 neonates were of a new Klebsiella sequence type, ST2727, and taxonomically belonged to K. quasipneumoniae subsp. similipneumoniae (formerly referred to as KpIIB). Genomic characterization showed that the isolated ST2727 strains had diverged from other K. quasipneumoniae subsp. similipneumoniae strains at least 90 years ago, whereas the neonatal samples were highly similar with a genomic divergence of 3.6 months. There was no relationship to the Klebsiella isolate from the staff member. This demonstrates that no transmission occurred from staff to patient or between patients. Rather, the data suggest that ST2727 colonized each neonate from a common hospital source. Sequence-based analysis of the genomes revealed several genes for antimicrobial resistance and some virulence features, but suggest that ST2727 is neither extremely-drug resistant nor hypervirulent. Our results highlight the clinical significance and genomic properties of ST2727 and urge genome-based measures be implemented for diagnostics and surveillance within hospital environments. Additionally, the present study demonstrates the need to scale the power of genomic analysis in retrospective studies where relatively few samples are available.

Pathogenicity Factors of Genomic Islands in Intestinal and Extraintestinal Escherichia coli

Escherichia coli is a versatile bacterial species that includes both harmless commensal strains and pathogenic strains found in the gastrointestinal tract in humans and warm-blooded animals. The growing amount of DNA sequence information generated in the era of "genomics" has helped to increase our understanding of the factors and mechanisms involved in the diversification of this bacterial species. The pathogenic side of E. coli that is afforded through horizontal transfers of genes encoding virulence factors enables this bacterium to become a highly diverse and adapted pathogen that is responsible for intestinal or extraintestinal diseases in humans and animals. Many of the accessory genes acquired by horizontal transfers form syntenic blocks and are recognized as genomic islands (GIs). These genomic regions contribute to the rapid evolution, diversification and adaptation of E. coli variants because they are frequently subject to rearrangements, excision and transfer, as well as to further acquisition of additional DNA. Here, we review a subgroup of GIs from E. coli termed pathogenicity islands (PAIs), a concept defined in the late 1980s by Jörg Hacker and colleagues in Werner Goebel's group at the University of Würzburg, Würzburg, Germany. As with other GIs, the PAIs comprise large genomic regions that differ from the rest of the genome by their G + C content, by their typical insertion within transfer RNA genes, and by their harboring of direct repeats (at their ends), integrase determinants, or other mobility loci. The hallmark of PAIs is their contribution to the emergence of virulent bacteria and to the development of intestinal and extraintestinal diseases. This review summarizes the current knowledge on the structure and functional features of PAIs, on PAI-encoded E. coli pathogenicity factors and on the role of PAIs in host-pathogen interactions.

Experimental approaches to tracking mobile genetic elements in microbial communities

Horizontal gene transfer is an important mechanism of microbial evolution and is often driven by the movement of mobile genetic elements between cells. Due to the fact that microbes live within communities, various mechanisms of horizontal gene transfer and types of mobile elements can co-occur. However, the ways in which horizontal gene transfer impacts and is impacted by communities containing diverse mobile elements has been challenging to address. Thus, the field would benefit from incorporating community-level information and novel approaches alongside existing methods. Emerging technologies for tracking mobile elements and assigning them to host organisms provide promise for understanding the web of potential DNA transfers in diverse microbial communities more comprehensively. Compared to existing experimental approaches, chromosome conformation capture and methylome analyses have the potential to simultaneously study various types of mobile elements and their associated hosts. We also briefly discuss how fermented food microbiomes, given their experimental tractability and moderate species complexity, make ideal models to which to apply the techniques discussed herein and how they can be used to address outstanding questions in the field of horizontal gene transfer in microbial communities.

Inside environmental Clostridium perfringens genomes: antibiotic resistance genes, virulence factors and genomic features

Until recently, research has focused on Clostridium perfringens in clinical settings without considering environmental isolates. In this study, environmental genomes were used to investigate possible antibiotic resistance and the presence of virulence traits in C. perfringens strains from raw surface water. In silico assembly of three C. perfringens strains, DNA generated almost complete genomes setting their length ranging from 3.4 to 3.6 Mbp with GC content of 28.18%. An average of 3,175 open reading frames was identified, with the majority associated with carbohydrate and protein metabolisms. The genomes harboured several antibiotic resistance genes for glycopeptides, macrolide-lincosamide-streptogramin B, β-lactam, trimethoprim, tetracycline and aminoglycosides and also the presence of several genes encoding for polypeptides and multidrug resistance efflux pumps and 35 virulence genes. Some of these encode for haemolysins, sialidase, hyaluronidase, collagenase, perfringolysin O and phospholipase C. All three genomes contained sequences indicating phage, antibiotic resistance and pathogenic islands integration sites. A genomic comparison of these three strains confirmed high similarity and shared core genes with clinical C. perfringens strains, highlighting their health security risks. This study provides a genomic insight into the potential pathogenicity of C. perfringens present in the environment and emphasises the importance of monitoring this niche in the future.

Microbial genomic island discovery, visualization and analysis

Horizontal gene transfer (also called lateral gene transfer) is a major mechanism for microbial genome evolution, enabling rapid adaptation and survival in specific niches. Genomic islands (GIs), commonly defined as clusters of bacterial or archaeal genes of probable horizontal origin, are of particular medical, environmental and/or industrial interest, as they disproportionately encode virulence factors and some antimicrobial resistance genes and may harbor entire metabolic pathways that confer a specific adaptation (solvent resistance, symbiosis properties, etc). As large-scale analyses of microbial genomes increases, such as for genomic epidemiology investigations of infectious disease outbreaks in public health, there is increased appreciation of the need to accurately predict and track GIs. Over the past decade, numerous computational tools have been developed to tackle the challenges inherent in accurate GI prediction. We review here the main types of GI prediction methods and discuss their advantages and limitations for a routine analysis of microbial genomes in this era of rapid whole-genome sequencing. An assessment is provided of 20 GI prediction software methods that use sequence-composition bias to identify the GIs, using a reference GI data set from 104 genomes obtained using an independent comparative genomics approach. Finally, we present guidelines to assist researchers in effectively identifying these key genomic regions.

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.

Genomic Analysis of Bovine Staphylococcus aureus Isolates from Milk To Elucidate Diversity and Determine the Distributions of Antimicrobial and Virulence Genes and Their Association with Mastitis

Staphylococcus aureus causes persistent clinical and subclinical bovine intramammary infections (IMI) worldwide. However, there is a lack of comprehensive information regarding genetic diversity, the presence of antimicrobial resistance (AMR), and virulence genes for S. aureus in bovine milk in Canada. Here, we performed whole-genome sequencing (WGS) of 119 Canadian bovine milk S. aureus isolates and determined they belonged to 8 sequence types (ST151, ST352, ST351, ST2187, ST2270, ST126, ST133, and ST8), 5 clonal complexes (CC151, CC97, CC126, CC133, and CC8), and 18 distinct Spa types. Pan-, core, and accessory genomes were composed of 6,340, 1,279, and 2,431 genes, respectively. Based on phenotypic screening for AMR, resistance was common against beta-lactams (19% of isolates) and sulfonamides (7% of isolates), whereas resistance against pirlimycin, tetracycline, ceftiofur, and erythromycin and to the combination of penicillin and novobiocin was uncommon (3, 3, 3, 2, and 2% of all isolates, respectively). We also determined distributions of 191 virulence factors (VFs) in 119 S. aureus isolates after classifying them into 5 functional categories (adherence [n = 28], exoenzymes [n = 21], immune evasion [n = 20], iron metabolism [n = 29], and toxins [n = 93]). Additionally, we calculated the pathogenic potential of distinct CCs and STs and determined that CC151 (ST151 and ST351) had the highest pathogenic potential (calculated by subtracting core-VFs from total VFs), followed by CC97 (ST352 and ST2187) and CC126 (ST126 and ST2270), potentially linked to their higher prevalence in bovine IMI worldwide. However, there was no statistically significant link between the presence of VF genes and mastitis.IMPORTANCE Staphylococcus aureus is a major cause of bovine intramammary infections, leading to significant economic losses to dairy industry in Canada and worldwide. There is a lack of knowledge regarding genetic diversity, the presence of antimicrobial resistance (AMR), and virulence genes for S. aureus isolated from bovine milk in Canada. Based on whole-genome sequencing and genomic analysis, we have determined the phylogeny and diversity of S. aureus in bovine milk and concluded that it had a large accessory genome, limited distribution of AMR genes, variable VF gene profiles and sequence types (ST), and clonal complex (CC)-specific pathogenic potentials. Comprehensive information on the population structure, as well as the virulence and resistance characteristics of S. aureus from bovine milk, will allow for source attribution, risk assessment, and improved therapeutic approaches in cattle.

ICEberg 2.0: an updated database of bacterial integrative and conjugative elements

ICEberg 2.0 (http://db-mml.sjtu.edu.cn/ICEberg/) is an updated database that provides comprehensive information about bacterial integrative and conjugative elements (ICEs). Compared with the previous version, three major improvements were made. First, with the aid of text mining and manual curation, it now recorded the details of 1032 ICEs, including 270 with experimental supports and 762 from bioinformatics prediction. Second, as increasing evidence has shown that ICEs frequently mobilize the so-called ‘hitchhikers’, such as integrative and mobilizable elements (IMEs) and cis-mobilizable elements (CIMEs), 83 known transfer interactions between 49 IMEs and 7 CIMEs with 19 ICEs taken from the literature were included and illustrated with visually intuitive directed graphs. An expanded collection of 260 chromosome-borne IMEs and 235 CIMEs was also added. At last, ICEberg 2.0 provides an online tool ICEfinder to predict ICEs or IMEs in bacterial genome sequences. It combines a similarity search for the integrase, relaxase and/or type IV secretion system and the co-localization of these corresponding homologous genes. With the recent updates, ICEberg 2.0 might provide better support for understanding the biological traits of ICEs, especially as their interaction with cognate mobilizable elements may further promote horizontal gene flow.

VFDB 2019: a comparative pathogenomic platform with an interactive web interface

The virulence factor database (VFDB, http://www.mgc.ac.cn/VFs/) is devoted to providing the scientific community with a comprehensive warehouse and online platform for deciphering bacterial pathogenesis. The various combinations, organizations and expressions of virulence factors (VFs) are responsible for the diverse clinical symptoms of pathogen infections. Currently, whole-genome sequencing is widely used to decode potential novel or variant pathogens both in emergent outbreaks and in routine clinical practice. However, the efficient characterization of pathogenomic compositions remains a challenge for microbiologists or physicians with limited bioinformatics skills. Therefore, we introduced to VFDB an integrated and automatic pipeline, VFanalyzer, to systematically identify known/potential VFs in complete/draft bacterial genomes. VFanalyzer first constructs orthologous groups within the query genome and preanalyzed reference genomes from VFDB to avoid potential false positives due to paralogs. Then, it conducts iterative and exhaustive sequence similarity searches among the hierarchical prebuilt datasets of VFDB to accurately identify potential untypical/strain-specific VFs. Finally, via a context-based data refinement process for VFs encoded by gene clusters, VFanalyzer can achieve relatively high specificity and sensitivity without manual curation. In addition, a thoroughly optimized interactive web interface is introduced to present VFanalyzer reports in comparative pathogenomic style for easy online analysis.

PHAST, PHASTER and PHASTEST: Tools for finding prophage in bacterial genomes

PHAST (PHAge Search Tool) and its successor PHASTER (PHAge Search Tool - Enhanced Release) have become two of the most widely used web servers for identifying putative prophages in bacterial genomes. Here we review the main capabilities of these web resources, provide some practical guidance regarding their use and discuss possible future improvements. PHAST, which was first described in 2011, made its debut just as whole bacterial genome sequencing and was becoming inexpensive and relatively routine. PHAST quickly gained popularity among bacterial genome researchers because of its web accessibility, its ease of use along with its enhanced accuracy and rapid processing times. PHASTER, which appeared in 2016, provided a number of much-needed enhancements to the PHAST server, including greater processing speed (to cope with very large submission volumes), increased database sizes, a more modern user interface, improved graphical displays and support for metagenomic submissions. Continuing developments in the field, along with increased interest in automated phage and prophage finding, have already led to several improvements to the PHASTER server and will soon lead to the development of a successor to PHASTER (to be called PHASTEST).

PRAP: Pan Resistome analysis pipeline

Background

Antibiotic resistance genes (ARGs) can spread among pathogens via horizontal gene transfer, resulting in imparities in their distribution even within the same species. Therefore, a pan-genome approach to analyzing resistomes is necessary for thoroughly characterizing patterns of ARGs distribution within particular pathogen populations. Software tools are readily available for either ARGs identification or pan-genome analysis, but few exist to combine the two functions.

Results

We developed Pan Resistome Analysis Pipeline (PRAP) for the rapid identification of antibiotic resistance genes from various formats of whole genome sequences based on the CARD or ResFinder databases. Detailed annotations were used to analyze pan-resistome features and characterize distributions of ARGs. The contribution of different alleles to antibiotic resistance was predicted by a random forest classifier. Results of analysis were presented in browsable files along with a variety of visualization options. We demonstrated the performance of PRAP by analyzing the genomes of 26 Salmonella enterica isolates from Shanghai, China.

Conclusions

PRAP was effective for identifying ARGs and visualizing pan-resistome features, therefore facilitating pan-genomic investigation of ARGs. This tool has the ability to further excavate potential relationships between antibiotic resistance genes and their phenotypic traits.

Emergence of an NDM-5-Producing Hypervirulent Klebsiella pneumoniae Sequence Type 35 Strain with Chromosomal Integration of an Integrative and Conjugative Element, ICEKp1

Here, we report an NDM-5-producing sequence type 35 (ST35) hypervirulent Klebsiella pneumoniae strain, isolated from the blood of a male patient. It showed a remarkable resistance to serum killing and neutrophil phagocytosis and high virulence in a mouse peritonitis infection model. Instead of carrying a pLVPK-like virulence plasmid, chromosomal integration of ICEKp1 (∼76 kb) was identified in a specific asparagine-tRNA gene, harboring the iron acquisition system salmochelin genes (iroBCDN), a yersiniabactin gene, and a variant of the rmpA gene.

Towards Exploring Toxin-Antitoxin Systems in Geobacillus: A Screen for Type II Toxin-Antitoxin System Families in a Thermophilic Genus

The toxin-antitoxin (TA) systems have been attracting attention due to their role in regulating stress responses in prokaryotes and their biotechnological potential. Much recognition has been given to type II TA system of mesophiles, while thermophiles have received merely limited attention. Here, we are presenting the putative type II TA families encoded on the genomes of four Geobacillus strains. We employed the TA finder tool to mine for TA-coding genes and manually curated the results using protein domain analysis tools. We also used the NCBI BLAST, Operon Mapper, ProOpDB, and sequence alignment tools to reveal the geobacilli TA features. We identified 28 putative TA pairs, distributed over eight TA families. Among the identified TAs, 15 represent putative novel toxins and antitoxins, belonging to the MazEF, MNT-HEPN, ParDE, RelBE, and XRE-COG2856 TA families. We also identified a potentially new TA composite, AbrB-ParE. Furthermore, we are suggesting the Geobacillus acetyltransferase TA (GacTA) family, which potentially represents one of the unique TA families with a reverse gene order. Moreover, we are proposing a hypothesis on the xre-cog2856 gene expression regulation, which seems to involve the c-di-AMP. This study aims for highlighting the significance of studying TAs in Geobacillus and facilitating future experimental research.

Staphylococcus arlettae Genomics: Novel Insights on Candidate Antibiotic Resistance and Virulence Genes in an Emerging Opportunistic Pathogen

Coagulase Negative Staphylococci (CoNS) are becoming increasingly recognized as an important cause of human and animal infections. Notwithstanding their clinical relevance, annotation of genes potentially involved in pathogenicity and/or antibiotic resistance in the CoNS species Staphylococcus arlettae (SAR) is currently very limited. In the current work we describe the genome of a novel methicillin resistant isolate of SAR, which we named Bari, and present a comprehensive analysis of predicted antibiotic resistance profiles and virulence determinants for all the 22 currently available SAR genomes. By comparing predicted antibiotic resistance and virulence-associated genes with those obtained from a manual selection of 148 bacterial strains belonging to 14 different species of staphylococci and to two “outgroup” species, Bacillus subtilis (BS) and Macrococcus caseoliticus (MC), we derived some interesting observations concerning the types and number of antibiotic resistance-related and virulence-like genes in SAR. Interestingly, almost 50% of the putative antibiotic resistance determinants identified in this work, which include the clinically relevant mec, van, and cls genes, were shared among all the SAR strains herein considered (Bari included). Moreover, comparison of predicted antibiotic resistance profiles suggest that SAR is closely related to well-known pathogenic Staphylococcus species, such as Staphylococcus aureus (SA) and Staphylococcus epidermidis (SE). A similar analysis of predicted virulence factors, revealed that several genes associated with pathogenesis (including, for example, ica, nuc, and ssp), which are commonly found in the genomes of pathogenic staphylococci such as Staphylococcus haemolyticus (SH) and Staphylococcus saprophyticus (SS), are observed also in the SAR strains for which a genomic sequence is available. All in all, we believe that the analyses presented in the current study, by providing a consistent and comprehensive annotation of virulence and antibiotic resistance-related genes in SAR, can constitute a valuable resource for the study of molecular mechanisms of opportunistic pathogenicity in this species.

Identification and characterization of OmpT-like proteases in uropathogenic Escherichia coli clinical isolates

Bacterial colonization of the urogenital tract is limited by innate defenses, including the production of antimicrobial peptides (AMPs). Uropathogenic Escherichia coli (UPEC) resist AMP‐killing to cause a range of urinary tract infections (UTIs) including asymptomatic bacteriuria, cystitis, pyelonephritis, and sepsis. UPEC strains have high genomic diversity and encode numerous virulence factors that differentiate them from non‐UTI‐causing strains, including ompT. As OmpT homologs cleave and inactivate AMPs, we hypothesized that UPEC strains from patients with symptomatic UTIs have high OmpT protease activity. Therefore, we measured OmpT activity in 58 clinical E. coli isolates. While heterogeneous OmpT activities were observed, OmpT activity was significantly greater in UPEC strains isolated from patients with symptomatic infections. Unexpectedly, UPEC strains exhibiting the greatest protease activities harbored an additional ompT‐like gene called arlC (ompTp). The presence of two OmpT‐like proteases in some UPEC isolates led us to compare the substrate specificities of OmpT‐like proteases found in E. coli. While all three cleaved AMPs, cleavage efficiency varied on the basis of AMP size and secondary structure. Our findings suggest the presence of ArlC and OmpT in the same UPEC isolate may confer a fitness advantage by expanding the range of target substrates.

oriTfinder: a web-based tool for the identification of origin of transfers in DNA sequences of bacterial mobile genetic elements

oriTfinder is a web server that facilitates the rapid identification of the origin of transfer site (oriT) of a conjugative plasmid or chromosome-borne integrative and conjugative element. The utilized back-end database oriTDB was built upon more than one thousand known oriT regions of bacterial mobile genetic elements (MGEs) as well as the known MGE-encoding relaxases and type IV coupling proteins (T4CP). With a combination of similarity searches for the oriTDB-archived oriT nucleotide sequences and the co-localization of the flanking relaxase homologous genes, the oriTfinder can predict the oriT region with high accuracy in the DNA sequence of a bacterial plasmid or chromosome in minutes. The server also detects the other transfer-related modules, including the potential relaxase gene, T4CP gene and the type IV secretion system gene cluster, and the putative genes coding for virulence factors and acquired antibiotic resistance determinants. oriTfinder may contribute to meeting the increasing demands of re-annotations for bacterial conjugative, mobilizable or non-transferable elements and aid in the rapid risk accession of disease-relevant trait dissemination in pathogenic bacteria of interest. oriTfinder is freely available to all users without any login requirement at http://bioinfo-mml.sjtu.edu.cn/oriTfinder.

TADB 2.0: an updated database of bacterial type II toxin-antitoxin loci

TADB2.0 (http://bioinfo-mml.sjtu.edu.cn/TADB2/) is an updated database that provides comprehensive information about bacterial type II toxin–antitoxin (TA) loci. Compared with the previous version, the database refined and the new data schema is employed. With the aid of text mining and manual curation, it recorded 6193 type II TA loci in 870 replicons of bacteria and archaea, including 105 experimentally validated TA loci. In addition, the newly developed tool TAfinder combines the homolog searches and the operon structure detection, allowing the prediction for type II TA pairs in bacterial genome sequences. It also helps to investigate the genomic context of predicted TA loci for putative virulence factors, antimicrobial resistance determinants and mobile genetic elements via alignments to the specific public databases. Additionally, the module TAfinder-Compare allows comparing the presence of the given TA loci across the close relative genomes. With the recent updates, TADB2.0 might provide better support for understanding the important roles of type II TA systems in the prokaryotic life activities.

Insights into Klebsiella pneumoniae type VI secretion system transcriptional regulation

Background

Klebsiella pneumoniae (KP) is an opportunistic pathogen that mainly causes respiratory and urinary tract infections. The frequent occurrence of simultaneously virulent and multiple drug-resistant isolates led WHO to include this species in the list of top priorities for research and development of therapeutic alternatives. The comprehensive knowledge of the molecular mechanisms underlying KP virulence may lead to the proposal of more efficient and specific drugs. One of its virulence factors is the Type VI Secretion System (T6SS), which contributes to bacterial competition, cell invasion and in vivo colonisation. Despite the few studies showing the involvement of T6SS in KP pathogenesis, little is known concerning the regulation of its expression. The understanding of regulatory mechanisms may give more clues about the function of the system and the possibilities of future interference in this process. This work aimed to standardise the annotation of T6SS genes in KP strains and identify mechanisms of their transcriptional regulation through computational predictions.

Results

We analyzed the genomes of Kp52.145, HS11286 and NTUH-K2044 strains to perform a broad prediction and re-annotation of T6SS genes through similarity searches, comparative and linear discriminant analysis. 38 genes were found in Kp52.145, while 29 in HS11286 and 30 in NTUH-K2044. Genes coding for iron uptake systems are encoded in adjacencies of T6SS, suggesting that KP T6SS might also play a role in ion import.

Some of the T6SS genes are comprised in syntenic regions. 17 sigma 70-dependent promoter regions were identified in Kp52.145, 12 in HS11286 and 12 in NTUH-K2044. Using VirtualFootprint algorithm, binding sites for 13 transcriptional regulators were found in Kp52.145 and 9 in HS11286 and 17 in NTUH-K2044. Six of them are common to the 3 strains: OxyR, H-NS, RcsAB, GcvA, Fis, and OmpR.

Conclusions

The data presented herein are derived from computational analysis. Although future experimental studies are required to confirm those predictions, they suggest that KP T6SS might be regulated in response to environmental signals that are indeed sensed by the bacteria inside the human host: temperature (H-NS), nutrition-limitation (GcvA and Fis), oxidative stress (OxyR) and osmolarity (RscAB and OmpR).

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5885-9) contains supplementary material, which is available to authorized users.

Staphylococcus debuckii sp. nov., a coagulase-negative species from bovine milk

A novel type strain, designated SDB 2975^T (=CECT 9737^T=DSM 105892^T), of the novel species Staphylococcus debuckii sp. nov. isolated from bovine milk is described. The novel species belongs to the genus Staphylococcus and showed resistance to tetracycline and was oxidase- and coagulase-negative, catalase-positive, and Gram-stain-positive. Phylogenetic relationships of Staphylococcus debuckii SDB 2975^T to other staphylococcal species were inferred from 16S rRNA gene and whole-genome-based phylogenetic reconstruction. The 16S rRNA gene comparisons showed that the strain is closely related to Staphylococcus condimenti (99.73 %), Staphylococcus piscifermentans (99.66 %), Staphylococcus carnosus (99.59 %) and Staphylococcus simulans (98.03 %). Average nucleotide identity (ANI) values between S.taphylococcus debuckii SDB 2975^T and its closely related Staphylococcus species were 83.96, 94.5, 84.03 and 78.09 %, respectively, and digital DNA-DNA hybridization (dDDH) values were 27.70, 58.02, 27.70 and 22.00 %, respectively. The genome of Staphylococcus debuckii SDB 2975^T was sequenced with PacBio and Illumina technologies and is 2 691 850 bp long, has a G+C content of 36.6 mol% and contains 2678 genes and 80 RNAs, including six copies of each5S rRNA, 16S rRNA and 23S rRNA genes. Biochemical profiling and a newly developed PCR assay enabled differentiation of Staphylococcus debuckii SDB 2975^T and three other SDB strains from its closest staphylococcal species. Differentiation was also achieved by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF). Genes unique to Staphylococcus debuckii were identified and a PCR-based assay was developed to differentiate Staphylococcus debuckii from other staphylococcal species. In conclusion, the results of phylogenetic analysis along with the ANI values <95 %, and dDDH values <70 % from closely related species along with the phenotypic and biochemical characteristics and specific MALDI-TOF profiles demonstrated that Staphylococcus debuckii SDB 2975^T represents a novel species within the genus Staphylococcus, named Staphylococcus debuckii sp. nov. (SDB 2975^T=CECT 9737^T=DSM 105892^T).

Comprehensive Virulence Gene Profiling of Bovine Non-aureus Staphylococci Based on Whole-Genome Sequencing Data

Non-aureus staphylococci (NAS) are the most frequently isolated pathogens from milk in dairy cattle worldwide. The virulence factors (VFs) and mechanisms by which these bacteria cause udder infection are not fully known. We determined the distribution and associations of 191 VFs in 25 NAS species and investigated the relationship between VFs and disease. Although the overall number of VFs was not associated with disease severity, increasing numbers of toxin and host immune evasion genes specifically were associated with more severe disease outcomes. These findings suggest that the development of disease and the interactions of VFs with the host are complex and determined by the interplay of genes rather than just the presence of virulence genes. Together, our results provide foundational genetic knowledge to other researchers to design and conduct further experiments, focusing on understanding the synergy between VFs and roles of individual NAS species in IMI and characterizing species-specific effects on udder health.

Non-aureus staphylococci (NAS) are the most frequently isolated pathogens from intramammary infection (IMI) in dairy cattle. Virulence factors (VFs) and mechanisms by which NAS cause IMI are not fully known. Herein, we analyzed the distribution of 191 VFs in 441 genomes of 25 NAS species, after classifying VFs into functional categories: adherence (n = 28), exoenzymes (n = 21), immune evasion (n = 20), iron metabolism (n = 29), and toxins (n = 93). In addition to establishing VF gene profiles, associations of VF genes between and among functional categories were computed, revealing distinctive patterns of association among VFs for various NAS species. Associations were also computed for low, medium, and high somatic cell count (SCC) and clinical mastitis (CM) isolates, demonstrating distinctive patterns of associations for low SCC and CM isolates, but no differences between high SCC and CM isolates. To determine whether VF distributions had any association with SCC or CM, various clustering approaches, including complete linkages, Ward clustering, and t-distributed stochastic neighbor embedding, were applied. However, no clustering of isolates representing low SCC, medium SCC, or high SCC or CM was identified. Regression analysis to test for associations with individual VF functional categories demonstrated that each additional toxin and host immune evasion gene increased the odds of having high SCC or CM, although an overall increase in the number of VFs was not associated with increased SCC or occurrence of CM. In conclusion, we established comprehensive VF gene profiling, determined VF gene distributions and associations, calculated pathogenic potentials of all NAS species, and detected no clear link between VF genes and mastitis.

IMPORTANCE Non-aureus staphylococci (NAS) are the most frequently isolated pathogens from milk in dairy cattle worldwide. The virulence factors (VFs) and mechanisms by which these bacteria cause udder infection are not fully known. We determined the distribution and associations of 191 VFs in 25 NAS species and investigated the relationship between VFs and disease. Although the overall number of VFs was not associated with disease severity, increasing numbers of toxin and host immune evasion genes specifically were associated with more severe disease outcomes. These findings suggest that the development of disease and the interactions of VFs with the host are complex and determined by the interplay of genes rather than just the presence of virulence genes. Together, our results provide foundational genetic knowledge to other researchers to design and conduct further experiments, focusing on understanding the synergy between VFs and roles of individual NAS species in IMI and characterizing species-specific effects on udder health.