T346Hunter: a novel web-based tool for the prediction of type III, type IV and type VI secretion systems in bacterial genomes

Multidrug resistance among uropathogenic clonal group A E. Coli isolates from Pakistani women with uncomplicated urinary tract infections

OBJECTIVE

Multi-drug resistance (MDR) has notably increased in community acquired uropathogens causing urinary tract infections (UTIs), predominantly Escherichia coli. Uropathogenic E. coli causes 80% of uncomplicated community acquired UTIs, particularly in pre-menopausal women. Considering this high prevalence and the potential to spread antimicrobial resistant genes, the current study was conducted to investigate the presence of clinically important strains of E. coli in Pakistani women having uncomplicated cystitis and pyelonephritis. Women belonging to low-income groups were exclusively included in the study. Seventy-four isolates from urine samples were processed, phylotyped, and screened for the presence of two Single Nucleotide Polymorphisms (SNPs) particularly associated with a clinically important clonal group A of E. coli (CgA) followed by antibiotic susceptibility testing and genome sequence analysis.

RESULTS

Phylogroup B2 was most prevalent in patients and 44% of isolates were positive for the presence of CgA specific SNPs in Fumarate hydratase and DNA gyrase subunit B genes. Antibiotic susceptibility testing showed widespread resistance to trimethoprim-sulfamethoxazole and extended-spectrum beta-lactamase production. The infection analysis revealed the phylogroup B2 to be more pathogenic as compared to the other groups. The genome sequence of E. coli strain U17 revealed genes encoding virulence, multidrug resistance, and host colonization mechanisms.

CONCLUSIONS

Our research findings not only validate the significant occurrence of multidrug-resistant clonal group A E. coli (CgA) in premenopausal Pakistani women suffering from cystitis and pyelonephritis but also reveal the presence of genes associated withvirulence, and drug efflux pumps. The detection of highly pathogenic, antimicrobial-resistant phylogroup B2 and CgA E. coli strains is likely to help in understanding the epidemiology of the pathogen and may ultimately help to reduce the impact of these strains on human health. Furthermore, the findings of this study will particularly help to reduce the prevalence of uncomplicated UTIs and the cost associated with their treatment in women belonging to low-income groups.

T9GPred: A Comprehensive Computational Tool for the Prediction of Type 9 Secretion System, Gliding Motility, and the Associated Secreted Proteins

Type 9 secretion system (T9SS) is one of the least characterized secretion systems exclusively found in the Bacteroidetes phylum, which comprises various environmental and economically relevant bacteria. While T9SS plays a central role in bacterial movement termed gliding motility, survival, and pathogenicity, there is an unmet need for a comprehensive tool that predicts T9SS, gliding motility, and proteins secreted via T9SS. In this study, we develop such a computational tool, Type 9 secretion system and Gliding motility Prediction (T9GPred). To build this tool, we manually curated published experimental evidence and identified mandatory components for T9SS and gliding motility prediction. We also compiled experimentally characterized proteins secreted via T9SS and determined the presence of three unique types of C-terminal domain signals, and these insights were leveraged to predict proteins secreted via T9SS. Notably, using recently published experimental evidence, we show that T9GPred has high predictive power. Thus, we used T9GPred to predict the presence of T9SS, gliding motility, and associated secreted proteins across 693 completely sequenced Bacteroidetes strains. T9GPred predicted 402 strains to have T9SS, of which 327 strains are also predicted to exhibit gliding motility. Further, T9GPred also predicted putative secreted proteins for the 402 strains. In a nutshell, T9GPred is a novel computational tool for systems-level prediction of T9SS and streamlining future experimentation. The source code of the computational tool is available in our GitHub repository: https://github.com/asamallab/T9GPred. The tool and its predicted results are compiled in a web server available at: https://cb.imsc.res.in/t9gpred/.

Influence of mobile genetic elements and insertion sequences in long- and short-term adaptive processes of Acidithiobacillus ferrooxidans strains

The recent revision of the Acidithiobacillia class using genomic taxonomy methods has shown that, in addition to the existence of previously unrecognized genera and species, some species of the class harbor levels of divergence that are congruent with ongoing differentiation processes. In this study, we have performed a subspecies-level analysis of sequenced strains of Acidithiobacillus ferrooxidans to prove the existence of distinct sublineages and identify the discriminant genomic/genetic characteristics linked to these sublineages, and to shed light on the processes driving such differentiation. Differences in the genomic relatedness metrics, levels of synteny, gene content, and both integrated and episomal mobile genetic elements (MGE) repertoires support the existence of two subspecies-level taxa within A. ferrooxidans. While sublineage 2A harbors a small plasmid related to pTF5, this episomal MGE is absent in sublineage 2B strains. Likewise, clear differences in the occurrence, coverage and conservation of integrated MGEs are apparent between sublineages. Differential MGE-associated gene cargo pertained to the functional categories of energy metabolism, ion transport, cell surface modification, and defense mechanisms. Inferred functional differences have the potential to impact long-term adaptive processes and may underpin the basis of the subspecies-level differentiation uncovered within A. ferrooxidans. Genome resequencing of iron- and sulfur-adapted cultures of a selected 2A sublineage strain (CCM 4253) showed that both episomal and large integrated MGEs are conserved over twenty generations in either growth condition. In turn, active insertion sequences profoundly impact short-term adaptive processes. The ISAfe1 element was found to be highly active in sublineage 2A strain CCM 4253. Phenotypic mutations caused by the transposition of ISAfe1 into the pstC2 encoding phosphate-transport system permease protein were detected in sulfur-adapted cultures and shown to impair growth on ferrous iron upon the switch of electron donor. The phenotypic manifestation of the △pstC2 mutation, such as a loss of the ability to oxidize ferrous iron, is likely related to the inability of the mutant to secure the phosphorous availability for electron transport-linked phosphorylation coupled to iron oxidation. Depletion of the transpositional △pstC2 mutation occurred concomitantly with a shortening of the iron-oxidation lag phase at later transfers on a ferrous iron-containing medium. Therefore, the pstII operon appears to play an essential role in A. ferrooxidans when cells oxidize ferrous iron. Results highlight the influence of insertion sequences and both integrated and episomal mobile genetic elements in the short- and long-term adaptive processes of A. ferrooxidans strains under changing growth conditions.

Density-based binning of gene clusters to infer function or evolutionary history using GeneGrouper

MOTIVATION

Identifying variant forms of gene clusters of interest in phylogenetically proximate and distant taxa can help to infer their evolutionary histories and functions. Conserved gene clusters may differ by only a few genes, but these small differences can in turn induce substantial phenotypes, such as by the formation of pseudogenes or insertions interrupting regulation. Particularly as microbial genomes and metagenomic assemblies become increasingly abundant, unsupervised grouping of similar, but not necessarily identical, gene clusters into consistent bins can provide a population-level understanding of their gene content variation and functional homology.

RESULTS

We developed GeneGrouper, a command-line tool that uses a density-based clustering method to group gene clusters into bins. GeneGrouper demonstrated high recall and precision in benchmarks for the detection of the 23-gene Salmonella enterica LT2 Pdu gene cluster and four-gene Pseudomonas aeruginosa PAO1 Mex gene cluster among 435 genomes spanning mixed taxa. In a subsequent application investigating the diversity and impact of gene-complete and -incomplete LT2 Pdu gene clusters in 1130 S.enterica genomes, GeneGrouper identified a novel, frequently occurring pduN pseudogene. When investigated in vivo, introduction of the pduN pseudogene negatively impacted microcompartment formation. We next demonstrated the versatility of GeneGrouper by clustering distant homologous gene clusters and variable gene clusters found in integrative and conjugative elements.

AVAILABILITY AND IMPLEMENTATION

GeneGrouper software and code are publicly available at https://pypi.org/project/GeneGrouper/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Predicting Type III Effector Proteins Using the Effectidor Web Server

Various Gram-negative bacteria use secretion systems to secrete effector proteins that manipulate host biochemical pathways to their benefit. We and others have previously developed machine-learning algorithms to predict novel effectors. Specifically, given a set of known effectors and a set of known non-effectors, the machine-learning algorithm extracts features that distinguish these two protein groups. In the training phase, the machine learning learns how to best combine the features to separate the two groups. The trained machine learning is then applied to open reading frames (ORFs) with unknown functions, resulting in a score for each ORF, which is its likelihood to be an effector. We developed Effectidor, a web server for predicting type III effectors. In this book chapter, we provide a step-by-step introduction to the application of Effectidor, from selecting input data to analyzing the obtained predictions.

Preliminary studies on the different roles of T6SSs in pathogenicity of Pseudomonas plecoglossicida NB2011

Pseudomonas plecoglossicida, the causative agent of visceral granulomas in the large yellow croaker (Larimichthys crocea) in China, encodes three sets of type Ⅵ secretion systems (T6SS1-3). The purpose of this study was to characterize the different roles of T6SSs involved in infection. In-frame deletion of T6SSs was constructed, which resulted in 8 mutants. Competition against E. coli DH5α, virulence against the croaker and in vivo survival ability of the mutants were tested. The expression and secretion of Hcp by P. plecoglossicida NB2011 were investigated. The results showed T6SS2 mutant failed to inhibit the growth of E. coli, which is an indication of T6SS2 acting against environmental bacteria. The LD₅₀ value of T6SS1 mutant strongly increased; T6SS2 and T6SS3 mutants were similar to that of the wild type; and the virulence of double deletion or triple deletion mutant was drastically alleviated, indicating that T6SS1 being one of the major virulence factors, and T6SS2 and T6SS3 directly or indirectly being involved in the pathogenicity. T6SS1 mutant disappeared in the fish spleen in 3 days, while other strains kept increasing, indicating the T6SS1 stimulation bacteria replication in vivo. Hcp1 secreted at 12-28°C and Hcp2 secreted at 12-35°C, while Hcp3 secretion not detected in vitro. This study has thrown some insights on the understanding of pathogenicity mechanisms of this pathogen.

Pseudomonas savastanoi pv. mandevillae pv. nov., a Clonal Pathogen Causing an Emerging, Devastating Disease of the Ornamental Plant Mandevilla spp

Commercial production of the ornamental plant dipladenia (Mandevilla spp.) is threatened by dipladenia leaf and stem spot disease, caused by the bacterium Pseudomonas savastanoi. P. savastanoi includes four pathovars of woody hosts differentiated by a characteristic host range in olive, oleander, ash, and broom plants. However, isolates from dipladenia have not been ascribed to any particular lineage or P. savastanoi pathovar. Here we report that isolates from dipladenia represent a distinct, clonal lineage. First, dipladenia isolates display very similar plasmid profiles, including a plasmid encoding the iaaM gene for biosynthesis of indole-3-acetic acid. Second, multilocus sequence analysis and core genome single-nucleotide polymorphisms phylogenies showed a monophyletic origin for dipladenia isolates, which cluster with isolates from oleander (pathovar nerii) in a distinct clade well separated from other P. savastanoi strains. Metabolic profiling and cross-pathogenicity tests in olive, oleander, ash, broom, and dipladenia clearly distinguished dipladenia isolates from the four P. savastanoi pathovars. Comparative genomics of the draft genome sequence of the dipladenia strain Ph3 with the other four pathovars showed that Ph3 encodes very few strain-specific genes and a similar set of virulence genes to pv. nerii, including its repertoire of type III secretion system effectors. However, hierarchical clustering based on the catalog of effectors and their allelic variants clearly separated Ph3 from pv. nerii strains. Based on their distinctive pathogenicity profile, we propose a de novo pathovar for P. savastanoi isolates from dipladenia, P. savastanoi pv. mandevillae pv. nov., for which strain Ph3 (CFBP 8832^PT) has been designated as the pathotype strain.

Immunity proteins of dual nuclease T6SS effectors function as transcriptional repressors

Bacteria utilize type VI secretion system (T6SS) to deliver antibacterial toxins to target co-habiting bacteria. Here, we report that Burkholderia gladioli strain NGJ1 deploys certain T6SS effectors (TseTBg), having both DNase and RNase activities to kill target bacteria. RNase activity is prominent on NGJ1 as well as other bacterial RNA while DNase activity is pertinent to only other bacteria. The associated immunity (TsiTBg) proteins harbor non-canonical helix-turn-helix motifs and demonstrate transcriptional repression activity, similar to the antitoxins of type II toxin-antitoxin (TA) systems. Genome analysis reveals that homologs of TseTBg are either encoded as TA or T6SS effectors in diverse bacteria. Our results indicate that a new ORF (encoding a hypothetical protein) has evolved as a result of operonic fusion of TA type TseTBg homolog with certain T6SS-related genes by the action of IS3 transposable elements. This has potentially led to the conversion of a TA into T6SS effector in Burkholderia. Our study exemplifies that bacteria can recruit toxins of TA systems as T6SS weapons to diversify its arsenal to dominate during inter-bacterial competitions.

Computational prediction of secreted proteins in gram-negative bacteria

Gram-negative bacteria harness multiple protein secretion systems and secrete a large proportion of the proteome. Proteins can be exported to periplasmic space, integrated into membrane, transported into extracellular milieu, or translocated into cytoplasm of contacting cells. It is important for accurate, genome-wide annotation of the secreted proteins and their secretion pathways. In this review, we systematically classified the secreted proteins according to the types of secretion systems in Gram-negative bacteria, summarized the known features of these proteins, and reviewed the algorithms and tools for their prediction.

Learning transferable deep convolutional neural networks for the classification of bacterial virulence factors

MOTIVATION

Identification of virulence factors (VFs) is critical to the elucidation of bacterial pathogenesis and prevention of related infectious diseases. Current computational methods for VF prediction focus on binary classification or involve only several class(es) of VFs with sufficient samples. However, thousands of VF classes are present in real-world scenarios, and many of them only have a very limited number of samples available.

RESULTS

We first construct a large VF dataset, covering 3446 VF classes with 160 495 sequences, and then propose deep convolutional neural network models for VF classification. We show that (i) for common VF classes with sufficient samples, our models can achieve state-of-the-art performance with an overall accuracy of 0.9831 and an F1-score of 0.9803; (ii) for uncommon VF classes with limited samples, our models can learn transferable features from auxiliary data and achieve good performance with accuracy ranging from 0.9277 to 0.9512 and F1-score ranging from 0.9168 to 0.9446 when combined with different predefined features, outperforming traditional classifiers by 1-13% in accuracy and by 1-16% in F1-score.

AVAILABILITY AND IMPLEMENTATION

All of our datasets are made publicly available at http://www.mgc.ac.cn/VFNet/, and the source code of our models is publicly available at https://github.com/zhengdd0422/VFNet.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Delineating the key virulence factors and intraspecies divergence of Vibrio harveyi via whole-genome sequencing

Vibrio harveyi is one of the major pathogens in aquaculture. To identify the key virulence factors affecting pathogenesis of V. harveyi towards fish, we conducted a field investigation for three representative fish farms infected with V. harveyi. Multilocus sequence typing (MLST) and whole-genome sequencing were conducted to delineate the phylogenetic relationship and genetic divergence of V. harveyi. A total of 25 V. harveyi strains were isolated from the diseased fish and groundwater and were subtyped into 12 sequence types by MLST. Five virulence genes, mshB, pilA, hutR, ureB, and ureG, were variably presented in the sequenced strains. The virulence gene profiles strongly correlated with the distinct pathogenicity of V. harveyi strains, with a strain harboring all five genes exhibiting the highest virulence towards fish. Phenotype assay confirmed that reduced virulence correlated with decreased motility and biofilm formation ability. Additionally, three types of type VI secretion system, namely T6SS1, T6SS2, and T6SS3, were identified in V. harveyi strains, which can be classified into six, four, and 12 subtypes, respectively. In conclusion, the results indicated that the virulence level of V. harveyi is mainly determined by the above virulence genes, which may play vital roles in environmental adaptation for V. harveyi.

Genome-wide Analysis of Four Enterobacter cloacae complex type strains: Insights into Virulence and Niche Adaptation

Enterobacter cloacae complex (Ecc) species are widely distributed opportunistic pathogens mainly associated with humans and plants. In this study, the genomes of clinical isolates including E. hormaechei, E. kobei, and E. ludwigii and non-clinical isolate including E. nimipressuralis were analysed in combination with the genome of E. asburiae by using the reference strain E. cloacae subsp. cloacae ATCC 13047; the Ecc strains were tested on artificial sputum media (ASM), which mimics the host, to evaluate T6SS genes as a case study. All five Ecc strains were sequenced in our lab. Comparative genome analysis of the Ecc strains revealed that genes associated with the survival of Ecc strains, including genes of metal-requiring proteins, defence-associated genes and genes associated with general physiology, were highly conserved in the genomes. However, the genes involved in virulence and drug resistance, specifically those involved in bacterial secretion, host determination and colonization of different strains, were present in different genomic regions. For example, T6SS accessory and core components, T4SS, and multidrug resistance genes/efflux system genes seemed vital for the survival of Ecc strains in various environmental niches, such as humans and plants. Moreover, the ASM host-mimicking growth medium revealed significantly high expression of T6SS genes, including PrpC, which is a regulatory gene of the T6SS, in all tested Ecc strains compared to the control medium. The variations in T6SS gene expression in ASM vs. control showed that the ASM system represents a simple, reproducible and economical alternative to animal models for studies such as those aimed at understanding the divergence of Ecc populations. In summary, genome sequencing of clinical and environmental Ecc genomes will assist in understanding the epidemiology of Ecc strains, including the isolation, virulence characteristics, prevention and treatment of infectious disease caused by these broad-host-range niche-associated species.

Marine Sediments Illuminate Chlamydiae Diversity and Evolution

The bacterial phylum Chlamydiae is so far composed of obligate symbionts of eukaryotic hosts. Well known for Chlamydiaceae, pathogens of humans and other animals, Chlamydiae also include so-called environmental lineages that primarily infect microbial eukaryotes. Environmental surveys indicate that Chlamydiae are found in a wider range of environments than anticipated previously. However, the vast majority of this chlamydial diversity has been underexplored, biasing our current understanding of their biology, ecological importance, and evolution. Here, we report that previously undetected and active chlamydial lineages dominate microbial communities in deep anoxic marine sediments taken from the Arctic Mid-Ocean Ridge. Reaching relative abundances of up to 43% of the bacterial community, and a maximum diversity of 163 different species-level taxonomic units, these Chlamydiae represent important community members. Using genome-resolved metagenomics, we reconstructed 24 draft chlamydial genomes, expanding by over a third the known genomic diversity in this phylum. Phylogenomic analyses revealed several novel clades across the phylum, including a previously unknown sister lineage of the Chlamydiaceae, providing new insights into the origin of pathogenicity in this family. We were unable to identify putative eukaryotic hosts for these marine sediment chlamydiae, despite identifying genomic features that may be indicative of host-association. The high abundance and genomic diversity of Chlamydiae in these anoxic marine sediments indicate that some members could play an important, and thus far overlooked, ecological role in such environments and may indicate alternate lifestyle strategies.

Comparative Genomic Analysis Confirms Five Genetic Populations of the Select Agent, Rathayibacter toxicus

Rathayibacter toxicus is a Gram-positive, nematode-vectored bacterium that infects several grass species in the family Poaceae. Unique in its genus, R. toxicus has the smallest genome, possesses a complete CRISPR-Cas system, a vancomycin-resistance cassette, produces tunicamycin, a corynetoxin responsible for livestock deaths in Australia, and is designated a Select Agent in the United States. In-depth, genome-wide analyses performed in this study support the previously designated five genetic populations, with a core genome comprising approximately 80% of the genome for all populations. Results varied as a function of the type of analysis and when using different bioinformatics tools for the same analysis; e.g., some programs failed to identify specific genomic regions that were actually present. The software variance highlights the need to verify bioinformatics results by additional methods; e.g., PCR, mapping genes to genomes, use of multiple algorithms). These analyses suggest the following relationships among populations: RT-IV ↔ RT-I ↔ RT-II ↔ RT-III ↔ RT-V, with RT-IV and RT-V being the most unrelated. This is the most comprehensive analysis of R. toxicus that included populations RT-I and RT-V. Future studies require underrepresented populations and more recent isolates from varied hosts and geographic locations.

Isolation, Characterization, and Complete Genome Sequence of a Bradyrhizobium Strain Lb8 From Nodules of Peanut Utilizing Crack Entry Infection

In many legumes, the colonization of roots by rhizobia is via "root hair entry" and its molecular mechanisms have been extensively studied. However, the nodulation of peanuts (Arachis hypogaea L.) by Bradyrhizobium strains requires an intercellular colonization process called "crack entry," which is understudied. To understand the intercellular crack entry process, it is critical to develop the tools and resources related to the rhizobium in addition to focus on investigating the mechanisms of the plant host. In this study, we isolated a Bradyrhizobium sp. strain, Lb8 from peanut root nodules and sequenced it using PacBio long reads. The complete genome sequence was a circular chromosome of 8,718,147 base-pair (bp) with an average GC content of 63.14%. No plasmid sequence was detected in the sequenced DNA sample. A total of 8,433 potential protein-encoding genes, one rRNA cluster, and 51 tRNA genes were annotated. Fifty-eight percent of the predicted genes showed similarity to genes of known functions and were classified into 27 subsystems representing various biological processes. The genome shared 92% of the gene families with B. diazoefficens USDA 110^T. A presumptive symbiosis island of 778 Kb was detected, which included two clusters of nif and nod genes. A total of 711 putative protein-encoding genes were in this region, among which 455 genes have potential functions related to symbiotic nitrogen fixation and DNA transmission. Of 21 genes annotated as transposase, 16 were located in the symbiosis island. Lb8 possessed both Type III and Type IV protein secretion systems, and our work elucidated the association of flagellar Type III secretion systems in bradyrhizobia. These observations suggested that complex rearrangement, such as horizontal transfer and insertion of different DNA elements, might be responsible for the plasticity of the Bradyrhizobium genome.

Host Range Determinants of Pseudomonas savastanoi Pathovars of Woody Hosts Revealed by Comparative Genomics and Cross-Pathogenicity Tests

The study of host range determinants within the Pseudomonas syringae complex is gaining renewed attention due to its widespread distribution in non-agricultural environments, evidence of large variability in intra-pathovar host range, and the emergence of new epidemic diseases. This requires the establishment of appropriate model pathosystems facilitating integration of phenotypic, genomic and evolutionary data. Pseudomonas savastanoi pv. savastanoi is a model pathogen of the olive tree, and here we report a closed genome of strain NCPPB 3335, plus draft genome sequences of three strains isolated from oleander (pv. nerii), ash (pv. fraxini) and broom plants (pv. retacarpa). We then conducted a comparative genomic analysis of these four new genomes plus 16 publicly available genomes, representing 20 strains of these four P. savastanoi pathovars of woody hosts. Despite overlapping host ranges, cross-pathogenicity tests using four plant hosts clearly separated these pathovars and lead to pathovar reassignment of two strains. Critically, these functional assays were pivotal to reconcile phylogeny with host range and to define pathovar-specific genes repertoires. We report a pan-genome of 7,953 ortholog gene families and a total of 45 type III secretion system effector genes, including 24 core genes, four genes exclusive of pv. retacarpa and several genes encoding pathovar-specific truncations. Noticeably, the four pathovars corresponded with well-defined genetic lineages, with core genome phylogeny and hierarchical clustering of effector genes closely correlating with pathogenic specialization. Knot-inducing pathovars encode genes absent in the canker-inducing pv. fraxini, such as those related to indole acetic acid, cytokinins, rhizobitoxine, and a bacteriophytochrome. Other pathovar-exclusive genes encode type I, type II, type IV, and type VI secretion system proteins, the phytotoxine phevamine A, a siderophore, c-di-GMP-related proteins, methyl chemotaxis proteins, and a broad collection of transcriptional regulators and transporters of eight different superfamilies. Our combination of pathogenicity analyses and genomics tools allowed us to correctly assign strains to pathovars and to propose a repertoire of host range-related genes in the P. syringae complex.

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.

Plant-Pathogenic Agrobacterium tumefaciens Strains Have Diverse Type VI Effector-Immunity Pairs and Vary in In-Planta Competitiveness

The type VI secretion system (T6SS) is used by gram-negative bacteria to translocate effectors that can antagonize other bacterial cells. Models predict the variation in collections of effector and cognate immunity genes determine competitiveness and can affect the dynamics of populations and communities of bacteria. However, the outcomes of competition cannot be entirely explained by compatibility of effector-immunity (EI) pairs. Here, we characterized the diversity of T6SS loci of plant-pathogenic Agrobacterium tumefaciens and showed that factors other than EI pairs can impact interbacterial competition. All examined strains encode T6SS active in secretion and antagonism against Escherichia coli. The spectra of EI pairs as well as compositions of gene neighborhoods are diverse. Almost 30 in-planta competitions were tested between different genotypes of A. tumefaciens. Fifteen competitions between members of different species-level groups resulted in T6SS-dependent suppression in in-planta growth of prey genotypes. In contrast, ten competitions between members within species-level groups resulted in no significant effect on the growth of prey genotypes. One strain was an exceptional case and, despite encoding a functional T6SS and toxic effector protein, could not compromise the growth of the four tested prey genotypes. The data suggest T6SS-associated EI pairs can influence the competitiveness of strains of A. tumefaciens, but genetic features have a significant role on the efficacy of interbacterial antagonism.

Molecular epidemiological analysis and risk factors for acquisition of carbapenemase-producing Enterobacter cloacae complex in a Japanese university hospital

Background

To clarify the molecular epidemiology of carbapenem-resistant Enterobacter cloacae complex (CREC) and the risk factors for acquisition of carbapenemase-producing E. cloacae complex (CPEC).

Methods

Using clinical CREC isolates detected in a Japanese university hospital over 4 years, carbapenemase production was screened with phenotypic methods. Carbapenemase genes were analysed by PCR and sequencing. Molecular epidemiological analyses were conducted with repetitive extragenic palindromic (REP)-PCR and multilocus sequence typing (MLST). CRECs were identified to the subspecies level by hsp60 sequencing. Whole-genome sequencing of plasmids was conducted. A case-control study was performed to identify risk factors for acquisition of CPEC among patients with CREC.

Results

Thirty-nine CRECs including 20 CPECs carrying bla _IMP-1 were identified. Patients with CPEC had longer hospital stay before detection (26.5 days vs. 12 days, p = 0.008), a urinary catheter (odds ratio [OR], 5.36; 95% confidence interval [CI], 1.14-30.9; p = 0.023), or intubation (OR, 7.53; 95% CI, 1.47-53.8; p = 0.008) compared to patients without CPEC. Four genetically closely related CPEC clusters were observed, which showed that three of four CPEC clusters corresponded to E. asburiae (ST 53), E. hormaechei subsp. steigerwaltii (ST 113 and ST 1047) and E. cloacae subsp. cloacae (ST 513) by MLST and hsp60 sequencing. Seven representative plasmids shared structures with class I integron containing bla _IMP-1 and IncHI2A replicon type.

Conclusions

A longer hospital stay, presence of a urinary catheter, and intubation are risk factors for CPEC acquisition. In addition to horizontal transmission of genetically indistinguishable CPECs, IncHI2A plasmid carrying bla _IMP-1 appeared to be transferred among genetically different ECs.

An account of in silico identification tools of secreted effector proteins in bacteria and future challenges

Bacterial pathogens secrete numerous effector proteins via six secretion systems, type I to type VI secretion systems, to adapt to new environments or to promote virulence by bacterium-host interactions. Many computational approaches have been used in the identification of effector proteins before the subsequent experimental verification because they tolerate laborious biological procedures and are genome scale, automated and highly efficient. Prevalent examples include machine learning methods and statistical techniques. In this article, we summarize the computational progress toward predicting secreted effector proteins in bacteria, with an opening of an introduction of features that are used to discriminate effectors from non-effectors. The mechanism, contribution and deficiency of previous developed detection tools are presented, which are further benchmarked based on a curated testing data set. According to the results of benchmarking, potential improvements of the prediction performance are discussed, which include (1) more informative features for discriminating the effectors from non-effectors; (2) the construction of comprehensive training data set of the machine learning algorithms; (3) the advancement of reliable prediction methods and (4) a better interpretation of the mechanisms behind the molecular processes. The future of in silico identification of bacterial secreted effectors includes both opportunities and challenges.

SLING: a tool to search for linked genes in bacterial datasets

Gene arrays and operons that encode functionally linked proteins form the most basic unit of transcriptional regulation in bacteria. Rules that govern the order and orientation of genes in these systems have been defined; however, these were based on a small set of genomes that may not be representative. The growing availability of large genomic datasets presents an opportunity to test these rules, to define the full range and diversity of these systems, and to understand their evolution. Here we present SLING, a tool to Search for LINked Genes by searching for a single functionally essential gene, along with its neighbours in a rule-defined proximity (https://github.com/ghoresh11/sling/wiki). Examining this subset of genes enables us to understand the basic diversity of these genetic systems in large datasets. We demonstrate the utility of SLING on a clinical collection of enteropathogenic Escherichia coli for two relevant operons: toxin antitoxin (TA) systems and RND efflux pumps. By examining the diversity of these systems, we gain insight on distinct classes of operons which present variable levels of prevalence and ability to be lost or gained. The importance of this analysis is not limited to TA systems and RND pumps, and can be expanded to understand the diversity of many other relevant gene arrays.

Loss and Gain in the Evolution of the Salmonella enterica Serovar Gallinarum Biovar Pullorum Genome

Pullorum disease, an acute poultry septicemia caused by Salmonella Gallinarum biovar Pullorum, is fatal for young chickens and is a heavy burden on poultry industry. The pathogen is rare in most developed countries but still extremely difficult to eliminate in China. Efficient epidemiological surveillance necessitates clarifying the origin of the isolates from different regions and their phylogenic relationships. Genomic epidemiological analysis of 97 S. Pullorum strains was carried out to reconstruct the phylogeny and transmission history of S. Pullorum. Further analysis demonstrated that functional gene loss and acquisition occurred simultaneously throughout the evolution of S. Pullorum, both of which reflected adaptation to the changing environment. The result of our study will be helpful in surveillance and prevention of pullorum disease.

Salmonella enterica subspecies enterica serovar Gallinarum biovar Pullorum (S. Pullorum) is the etiological agent of pullorum disease, causing white diarrhea with high mortality in chickens. There are many unsolved issues surrounding the epidemiology of S. Pullorum, including its origin and transmission history as well as the discordance between its phenotypic heterogeneity and genetic monomorphism. In this paper, we report the results of whole-genome sequencing of a panel of 97 S. Pullorum strains isolated between 1962 and 2014 from four countries across three continents. We utilized 6,795 core genome single nucleotide polymorphisms (SNPs) to reconstruct a phylogenetic tree within a spatiotemporal Bayesian framework, estimating that the most recent common ancestor of S. Pullorum emerged in ∼914 CE (95% confidence interval [95%CI], 565 to 1273 CE). The extant S. Pullorum strains can be divided into four distinct lineages, each of which is significantly associated with geographical distribution. The intercontinental transmissions of lineages III and IV can be traced to the mid-19th century and are probably related to the “Hen Fever” prevalent at that time. Further genomic analysis indicated that the loss or pseudogenization of functional genes involved in metabolism and virulence in S. Pullorum has been ongoing since before and after divergence from the ancestor. In contrast, multiple prophages and plasmids have been acquired by S. Pullorum, and these have endowed it with new characteristics, especially the multidrug resistance conferred by two large plasmids in lineage I. The results of this study provide insight into the evolution of S. Pullorum and prove the efficiency of whole-genome sequencing in epidemiological surveillance of pullorum disease.

IMPORTANCE Pullorum disease, an acute poultry septicemia caused by Salmonella Gallinarum biovar Pullorum, is fatal for young chickens and is a heavy burden on poultry industry. The pathogen is rare in most developed countries but still extremely difficult to eliminate in China. Efficient epidemiological surveillance necessitates clarifying the origin of the isolates from different regions and their phylogenic relationships. Genomic epidemiological analysis of 97 S. Pullorum strains was carried out to reconstruct the phylogeny and transmission history of S. Pullorum. Further analysis demonstrated that functional gene loss and acquisition occurred simultaneously throughout the evolution of S. Pullorum, both of which reflected adaptation to the changing environment. The result of our study will be helpful in surveillance and prevention of pullorum disease.

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

VRprofile is a Web server that facilitates rapid investigation of virulence and antibiotic resistance genes, as well as extends these trait transfer-related genetic contexts, in newly sequenced pathogenic bacterial genomes. The used backend database MobilomeDB was firstly built on sets of known gene cluster loci of bacterial type III/IV/VI/VII secretion systems and mobile genetic elements, including integrative and conjugative elements, prophages, class I integrons, IS elements and pathogenicity/antibiotic resistance islands. VRprofile is thus able to co-localize the homologs of these conserved gene clusters using HMMer or BLASTp searches. With the integration of the homologous gene cluster search module with a sequence composition module, VRprofile has exhibited better performance for island-like region predictions than the other widely used methods. In addition, VRprofile also provides an integrated Web interface for aligning and visualizing identified gene clusters with MobilomeDB-archived gene clusters, or a variety set of bacterial genomes. VRprofile might contribute to meet the increasing demands of re-annotations of bacterial variable regions, and aid in the real-time definitions of disease-relevant gene clusters in pathogenic bacteria of interest. VRprofile is freely available at http://bioinfo-mml.sjtu.edu.cn/VRprofile.

Comparative genomics of Paraburkholderia kururiensis and its potential in bioremediation, biofertilization, and biocontrol of plant pathogens

Burkholderia harbors versatile Gram-negative species and is β-Proteobacteria. Recently, it was proposed to split the genus in two main branches: one of animal and plant pathogens and another, Paraburkholderia, harboring environmental and plant-beneficial species. Currently, Paraburkholderia comprises more than 70 species with ability to occupy very diverse environmental niches. Herein, we sequenced and analyzed the genome of Paraburkholderia kururiensis type strain KP23^T , and compared to P. kururiensis M130, isolated in Brazil, and P. kururiensis susbp. thiooxydans, from Korea. This study focused on the gene content of the three genomes with special emphasis on their potential of plant-association, biocontrol, and bioremediation. The comparative analyses revealed several genes related to plant benefits, including biosynthesis of IAA, ACC deaminase, multiple efflux pumps, dioxygenases, and degradation of aromatic compounds. Importantly, a range of genes for protein secretion systems (type III, IV, V, and VI) were characterized, potentially involved in P. kururiensis well documented ability to establish endophytic association with plants. These findings shed light onto bacteria-plant interaction mechanisms at molecular level, adding novel information that supports their potential application in bioremediation, biofertilization, and biocontrol of plant pathogens. P. kururiensis emerges as a promising model to investigate adaptation mechanisms in different ecological niches.

Whole-Genome Comparisons Among the Genus Shewanella Reveal the Enrichment of Genes Encoding Ankyrin-Repeats Containing Proteins in Sponge-Associated Bacteria

The bacterial members of the genus Shewanella are widely distributed and inhabit both freshwater and marine environments. Some members of Shewanella have gained considerable attention due to its ability to survive in redox-stratified environments. However, a gap of knowledge exists on the key genomic features of the sponge-associated Shewanella sp. involving the successful host-bacteria interaction, as sponge-symbiotic Shewanella are largely underrepresented in the public repositories. With the aim of identifying the genomic signatures of sponge-Shewanella association, we generated a high-quality genome data of a sponge-associated, Shewanella sp. OPT22, isolated from the intertidal marine sponge Ophlitaspongia papilla and performed comprehensive comparative analyses of 68 genome strains of the genus Shewanella including two previously reported genomes of sponge-associated bacteria, Shewanella spongiae KCTC 22492 and Shewanella sp. Alg231_23. The 16S rRNA-based phylogenetic reconstruction showed the well-supported affiliation of OPT22 and KCTC 22492 with previously reported sponge-associated bacteria, affirming the “sponge-specific” nature of these two bacterial strains isolated from different marine sponge species from the Atlantic and Pacific (East Sea) Oceans, respectively. The genome comparison of the 68 strains of Shewanella inhabiting different habitats revealed the unusual/previously unreported abundance of genes encoding for ankyrin-repeat containing proteins (ANKs) in the genomes of the two sponge-associated strains, OPT22 (ANKs; n = 45) and KCTC 22492 (ANKs; n = 52), which might be involved in sponge-Shewanella interactions. Focused analyses detected the syntenic organization of the gene cluster encoding major secretion system (type III/IV/VI) components and the presence of effector homologs in OPT22 and KCTC 22492 that seem to play a role in the virulence of the sponge bacteria. The genomic island (GI) of Shewanella sp. OPT22 was identified to localize a gene cluster encoding T4SS components and ANK (n = 1), whereas S. spongiae KCTC 22492 harbored a total of seven ANKs within multiple GIs. GIs may play a pivotal role in the dissemination of symbioses-related genes (ANKs) through the horizontal gene transfer, contributing to the diversification and adaptation of sponge-associated Shewanella. Overall, the genome analyses of Shewanella isolates from marine sponges revealed genomic repertoires that might be involved in establishing successful symbiotic relationships with the sponge hosts.

Tiny architects: biogenesis of intracellular replicative niches by bacterial pathogens

Co-evolution of bacterial pathogens with their hosts led to the emergence of a stunning variety of strategies aiming at the evasion of host defences, colonisation of host cells and tissues and, ultimately, the establishment of a successful infection. Pathogenic bacteria are typically classified as extracellular and intracellular; however, intracellular lifestyle comes in many different flavours: some microbes rapidly escape to the cytosol whereas other microbes remain within vacuolar compartments and harness membrane trafficking pathways to generate their host-derived, pathogen-specific replicative niche. Here we review the current knowledge on a variety of vacuolar lifestyles, the effector proteins used by bacteria as tools to take control of the host cell and the main membrane trafficking signalling pathways targeted by vacuolar pathogens as source of membranes and nutrients. Finally, we will also discuss how host cells have developed countermeasures to sense the biogenesis of the aberrant organelles harbouring bacteria. Understanding the dialogue between bacterial and eukaryotic proteins is the key to unravel the molecular mechanisms of infection and in turn, this may lead to the identification of new targets for the development of new antimicrobials.

Secretion Systems and Secreted Proteins in Gram-Negative Entomopathogenic Bacteria: Their Roles in Insect Virulence and Beyond

Many Gram-negative bacteria have evolved insect pathogenic lifestyles. In all cases, the ability to cause disease in insects involves specific bacterial proteins exported either to the surface, the extracellular environment, or the cytoplasm of the host cell. They also have several distinct mechanisms for secreting such proteins. In this review, we summarize the major protein secretion systems and discuss examples of secreted proteins that contribute to the virulence of a variety of Gram-negative entomopathogenic bacteria, including Photorhabdus, Xenorhabdus, Serratia, Yersinia, and Pseudomonas species. We also briefly summarize two classes of exported protein complexes, the PVC-like elements, and the Tc toxin complexes that were first described in entomopathogenic bacteria.

Typing and comparative genome analysis of Brucella melitensis isolated from Lebanon

Brucella melitensis is the main causative agent of the zoonotic disease brucellosis. This study aimed at typing and characterizing genetic variation in 33 Brucella isolates recovered from patients in Lebanon. Bruce-ladder multiplex PCR and PCR-RFLP of omp31, omp2a and omp2b were performed. Sixteen representative isolates were chosen for draft-genome sequencing and analyzed to determine variations in virulence, resistance, genomic islands, prophages and insertion sequences. Comparative whole-genome single nucleotide polymorphism analysis was also performed. The isolates were confirmed to be B. melitensis. Genome analysis revealed multiple virulence determinants and efflux pumps. Genome comparisons and single nucleotide polymorphisms divided the isolates based on geographical distribution but revealed high levels of similarity between the strains. Sequence divergence in B. melitensis was mainly due to lateral gene transfer of mobile elements. This is the first report of an in-depth genomic characterization of B. melitensis in Lebanon.

Spread of Plasmid-Encoded NDM-1 and GES-5 Carbapenemases among Extensively Drug-Resistant and Pandrug-Resistant Clinical Enterobacteriaceae in Durban, South Africa

Whole-genome sequence analyses revealed the presence of bla_NDM-1 (n = 31), bla_GES-5 (n = 8), bla_OXA-232 (n = 1), or bla_NDM-5 (n = 1) in extensively drug-resistant and pandrug-resistant Enterobacteriaceae organisms isolated from in-patients in 10 private hospitals (2012 to 2013) in Durban, South Africa. Two novel NDM-1-encoding plasmids from Klebsiella pneumoniae were circularized by PacBio sequencing. In p19-10_01 [IncFIB(K); 223.434 bp], bla_NDM-1 was part of a Tn1548-like structure (16.276 bp) delineated by IS26 The multireplicon plasmid p18-43_01 [IncR_1/IncFIB(pB171)/IncFII(Yp); 212.326 bp] shared an 80-kb region with p19-10_01, not including the bla_NDM-1-containing region. The two plasmids were used as references for tracing NDM-1-encoding plasmids in the other genome assemblies. The p19-10_01 sequence was detected in K. pneumoniae (n = 7) only, whereas p18-43_01 was tracked to K. pneumoniae (n = 4), Klebsiella michiganensis (n = 1), Serratia marcescens (n = 11), Enterobacter spp. (n = 7), and Citrobacter freundii (n = 1), revealing horizontal spread of this bla_NDM-1-bearing plasmid structure. Global phylogeny showed clustering of the K. pneumoniae (18/20) isolates together with closely related carbapenemase-negative ST101 isolates from other geographical origins. The South African isolates were divided into three phylogenetic subbranches, where each group had distinct resistance and replicon profiles, carrying either p19-10_01, p18-10_01, or pCHE-A1 (8,201 bp). The latter plasmid carried bla_GES-5 and aacA4 within an integron mobilization unit. Our findings imply independent plasmid acquisition followed by local dissemination. Additionally, we detected bla_OXA-232 carried by pPKPN4 in K. pneumoniae (ST14) and bla_NDM-5 contained by a pNDM-MGR194-like genetic structure in Escherichia coli (ST167), adding even more complexity to the multilayer molecular mechanisms behind nosocomial spread of carbapenem-resistant Enterobacteriaceae in Durban, South Africa.

The mutagenesis of a type IV secretion system locus of Piscirickettsia salmonis leads to the attenuation of the pathogen in Atlantic salmon, Salmo salar

Piscirickettsiosis is a threatening infectious disease for the salmon industry, due to it being responsible for significant economic losses. The control of outbreaks also poses considerable environmental challenges. Despite Piscirickettsia salmonis having been discovered as the aetiological agent of the disease more than 25 years ago, its pathogenicity remains poorly understood. Among virulence factors identified so far, type four secretion systems (T4SS) seem to play a key role during the infection caused by the bacterium. We report here the genetic manipulation of P. salmonis by means of the transference of plasmid DNA in mating assays. An insertion cassette was engineered for targeting the icmB gene, which encodes a putative T4SS-ATPase and is carried by one of the chromosomal T4SS clusters found within the genome of P. salmonis PM15972A1, a virulent representative of the EM-90-like strain. The molecular characterization of the resulting mutant strain demonstrated that the insertion interrupted the target gene. Further in vitro testing of the icmB mutant showed a dramatic drop in infectivity as tested in CHSE-214 cells, which is in agreement with its attenuated behaviour observed in vivo. Altogether, our results demonstrate that, similar to other facultative intracellular pathogens, P. salmonis' virulence relies on an intact T4SS.

Genome analysis of two novel Pseudomonas strains exhibiting differential hypersensitivity reactions on tobacco seedlings reveals differences in nonflagellar T3SS organization and predicted effector proteins

Multilocus sequence analysis (MLSA) of two new biological control strains (S1E40 and S3E12) of Pseudomonas was performed to assess their taxonomic position relative to close lineages, and comparative genomics employed to investigate whether these strains differ in key genetic features involved in hypersensitivity responses (HRs). Strain S3E12, at high concentration, incites HRs on tobacco and corn plantlets while S1E40 does not. Phylogenies based on individual genes and 16S rRNA-gyrB-rpoB-rpoD concatenated sequence data show strains S1E40 and S3E12 clustering in distinct groups. Strain S3E12 consistently clustered with Pseudomonas marginalis, a bacterium causing soft rots on plant tissues. MLSA data suggest that strains S1E40 and S3E12 are novel genotypes. This is consistent with the data of genome-based DNA-DNA homology values that are below the proposed cutoff species boundary. Comparative genomics analysis of the two strains revealed major differences in the type III secretion systems (T3SS) as well as the predicted T3SS secreted effector proteins (T3Es). One nonflagellar (NF-T3SS) and two flagellar T3SSs (F-T3SS) clusters were identified in both strains. While F-T3SS clusters in both strains were relatively conserved, the NF-T3SS clusters differed in the number of core components present. The predicted T3Es also differed in the type and number of CDSs with both strains having unique predicted protease-related effectors. In addition, the T1SS organization of the S3E12 genome has protein-coding sequences (CDSs) encoding for key factors such as T1SS secreted agglutinin repeats-toxins (a group of cytolysins and cytotoxins), a membrane fusion protein (LapC), a T1SS ATPase of LssB family (LapB), and T1SS-associated transglutaminase-like cysteine proteinase (LapP). In contrast, strain S1E40 has all CDSs for the seven-gene operon (pelA-pelG) required for Pel biosynthesis but not S3E12, suggesting that biofilm formation in these strains is modulated differently. The data presented here provide an insight of the genome organization of these two phytobacterial strains.

Indigenous Pseudomonas spp. Strains from the Olive (Olea europaea L.) Rhizosphere as Effective Biocontrol Agents against Verticillium dahliae: From the Host Roots to the Bacterial Genomes

The use of biological control agents (BCA), alone or in combination with other management measures, has gained attention over the past decades, driven by the need to seek for sustainable and eco-friendly alternatives to confront plant pathogens. The rhizosphere of olive (Olea europaea L.) plants is a source of bacteria with potential as biocontrol tools against Verticillium wilt of olive (VWO) caused by Verticillium dahliae Kleb. A collection of bacterial isolates from healthy nursery-produced olive (cultivar Picual, susceptible to VWO) plants was generated based on morphological, biochemical and metabolic characteristics, chemical sensitivities, and on their in vitro antagonistic activity against several olive pathogens. Three strains (PIC25, PIC105, and PICF141) showing high in vitro inhibition ability of pathogens' growth, particularly against V. dahliae, were eventually selected. Their effectiveness against VWO caused by the defoliating pathotype of V. dahliae was also demonstrated, strain PICF141 being the rhizobacteria showing the best performance as BCA. Genotypic and phenotypic traits traditionally associated with plant growth promotion and/or biocontrol abilities were evaluated as well (e.g., phytase, xylanase, catalase, cellulase, chitinase, glucanase activities, and siderophore and HCN production). Multi-locus sequence analyses of conserved genes enabled the identification of these strains as Pseudomonas spp. Strain PICF141 was affiliated to the “Pseudomonas mandelii subgroup,” within the “Pseudomonas fluorescens group,” Pseudomonas lini being the closest species. Strains PIC25 and PIC105 were affiliated to the “Pseudomonas aeruginosa group,” Pseudomonas indica being the closest relative. Moreover, we identified P. indica (PIC105) for the first time as a BCA. Genome sequencing and in silico analyses allowed the identification of traits commonly associated with plant-bacteria interactions. Finally, the root colonization ability of these olive rhizobacteria was assessed, providing valuable information for the future development of formulations based on these strains. A set of actions, from rhizosphere isolation to genome analysis, is proposed and discussed for selecting indigenous rhizobacteria as effective BCAs.

Pathogenic, phenotypic and molecular characterisation of Xanthomonas nasturtii sp. nov. and Xanthomonas floridensis sp. nov., new species of Xanthomonas associated with watercress production in Florida

We describe two new species of the genus Xanthomonas, represented by yellow mucoid bacterial strains isolated from diseased leaves of watercress (Nasturtium officinale) produced in Florida, USA. One strain was pathogenic on watercress, but not in other species including a range of brassicas; other strains were not pathogenic in any of the tested plants. Data from Biolog carbon source utilization tests and nucleotide sequence data from 16S and gyrB loci suggested that both pathogenic and non-pathogenic strains were related to, yet distinct from, previously described Xanthomonas species. Multilocus sequence analysis and whole genome-wide comparisons of the average nucleotide identity (ANI) of genomes of two strains from watercress showed that these are distinct and share less than 95 % ANI with all other known species; the non-pathogenic strain WHRI 8848 is close to Xanthomonascassavae (ANI of 93.72 %) whilst the pathogenic strain WHRI 8853 is close to a large clade of species that includes Xanthomonasvesicatoria (ANI ≤90.25 %). Based on these results, we propose that both strains represent new Xanthomonas species named Xanthomonas floridensis sp. nov. (type strain WHRI 8848=ATCC TSD-60=ICMP 21312=LMG 29665=NCPPB 4601) and Xanthomonas nasturtii sp. nov. (type strain WHRI 8853=ATCC TSD-61=ICMP 21313=LMG 29666=NCPPB 4600), respectively. The presence of non-pathogenic Xanthomonas strains in watercress and their interaction with pathogenic strains needs to be further investigated. Although the importance of the new pathogenic species is yet to be determined, the bacterial disease that it causes constitutes a threat to watercress production and its distribution should be monitored.

Complete Genome Sequences of Three Xanthomonas citri Strains from Texas

The complete genome sequences of three Xanthomonas citri strains isolated from lime trees in Texas were found to belong to the A^w group. All carried nearly identical large plasmids with similarity to those of a citrus canker strain from India and to xanthomonads from Africa and Colombia. All three strains harbored unusual pthA homologs.

Type VI Secretion Effectors: Methodologies and Biology

The type VI secretion system (T6SS) is a nanomachine deployed by many Gram-negative bacteria as a weapon against eukaryotic hosts or prokaryotic competitors. It assembles into a bacteriophage tail-like structure that can transport effector proteins into the environment or target cells for competitive survival or pathogenesis. T6SS effectors have been identified by a variety of approaches, including knowledge/hypothesis-dependent and discovery-driven approaches. Here, we review and discuss the methods that have been used to identify T6SS effectors and the biological and biochemical functions of known effectors. On the basis of the nature and transport mechanisms of T6SS effectors, we further propose potential strategies that may be applicable to identify new T6SS effectors.

Pan-Genomic Analysis Permits Differentiation of Virulent and Non-virulent Strains of Xanthomonas arboricola That Cohabit Prunus spp. and Elucidate Bacterial Virulence Factors

Xanthomonas arboricola is a plant-associated bacterial species that causes diseases on several plant hosts. One of the most virulent pathovars within this species is X. arboricola pv. pruni (Xap), the causal agent of bacterial spot disease of stone fruit trees and almond. Recently, a non-virulent Xap-look-a-like strain isolated from Prunus was characterized and its genome compared to pathogenic strains of Xap, revealing differences in the profile of virulence factors, such as the genes related to the type III secretion system (T3SS) and type III effectors (T3Es). The existence of this atypical strain arouses several questions associated with the abundance, the pathogenicity, and the evolutionary context of X. arboricola on Prunus hosts. After an initial characterization of a collection of Xanthomonas strains isolated from Prunus bacterial spot outbreaks in Spain during the past decade, six Xap-look-a-like strains, that did not clustered with the pathogenic strains of Xap according to a multi locus sequence analysis, were identified. Pathogenicity of these strains was analyzed and the genome sequences of two Xap-look-a-like strains, CITA 14 and CITA 124, non-virulent to Prunus spp., were obtained and compared to those available genomes of X. arboricola associated with this host plant. Differences were found among the genomes of the virulent and the Prunus non-virulent strains in several characters related to the pathogenesis process. Additionally, a pan-genomic analysis that included the available genomes of X. arboricola, revealed that the atypical strains associated with Prunus were related to a group of non-virulent or low virulent strains isolated from a wide host range. The repertoire of the genes related to T3SS and T3Es varied among the strains of this cluster and those strains related to the most virulent pathovars of the species, corylina, juglandis, and pruni. This variability provides information about the potential evolutionary process associated to the acquisition of pathogenicity and host specificity in X. arboricola. Finally, based in the genomic differences observed between the virulent and the non-virulent strains isolated from Prunus, a sensitive and specific real-time PCR protocol was designed to detect and identify Xap strains. This method avoids miss-identifications due to atypical strains of X. arboricola that can cohabit Prunus.

Pan-Genomic Analysis Permits Differentiation of Virulent and Non-virulent Strains of Xanthomonas arboricola That Cohabit Prunus spp. and Elucidate Bacterial Virulence Factors

Xanthomonas arboricola is a plant-associated bacterial species that causes diseases on several plant hosts. One of the most virulent pathovars within this species is X. arboricola pv. pruni (Xap), the causal agent of bacterial spot disease of stone fruit trees and almond. Recently, a non-virulent Xap-look-a-like strain isolated from Prunus was characterized and its genome compared to pathogenic strains of Xap, revealing differences in the profile of virulence factors, such as the genes related to the type III secretion system (T3SS) and type III effectors (T3Es). The existence of this atypical strain arouses several questions associated with the abundance, the pathogenicity, and the evolutionary context of X. arboricola on Prunus hosts. After an initial characterization of a collection of Xanthomonas strains isolated from Prunus bacterial spot outbreaks in Spain during the past decade, six Xap-look-a-like strains, that did not clustered with the pathogenic strains of Xap according to a multi locus sequence analysis, were identified. Pathogenicity of these strains was analyzed and the genome sequences of two Xap-look-a-like strains, CITA 14 and CITA 124, non-virulent to Prunus spp., were obtained and compared to those available genomes of X. arboricola associated with this host plant. Differences were found among the genomes of the virulent and the Prunus non-virulent strains in several characters related to the pathogenesis process. Additionally, a pan-genomic analysis that included the available genomes of X. arboricola, revealed that the atypical strains associated with Prunus were related to a group of non-virulent or low virulent strains isolated from a wide host range. The repertoire of the genes related to T3SS and T3Es varied among the strains of this cluster and those strains related to the most virulent pathovars of the species, corylina, juglandis, and pruni. This variability provides information about the potential evolutionary process associated to the acquisition of pathogenicity and host specificity in X. arboricola. Finally, based in the genomic differences observed between the virulent and the non-virulent strains isolated from Prunus, a sensitive and specific real-time PCR protocol was designed to detect and identify Xap strains. This method avoids miss-identifications due to atypical strains of X. arboricola that can cohabit Prunus.

SecretEPDB: a comprehensive web-based resource for secreted effector proteins of the bacterial types III, IV and VI secretion systems

Bacteria translocate effector molecules to host cells through highly evolved secretion systems. By definition, the function of these effector proteins is to manipulate host cell biology and the sequence, structural and functional annotations of these effector proteins will provide a better understanding of how bacterial secretion systems promote bacterial survival and virulence. Here we developed a knowledgebase, termed SecretEPDB (Bacterial Secreted Effector Protein DataBase), for effector proteins of type III secretion system (T3SS), type IV secretion system (T4SS) and type VI secretion system (T6SS). SecretEPDB provides enriched annotations of the aforementioned three classes of effector proteins by manually extracting and integrating structural and functional information from currently available databases and the literature. The database is conservative and strictly curated to ensure that every effector protein entry is supported by experimental evidence that demonstrates it is secreted by a T3SS, T4SS or T6SS. The annotations of effector proteins documented in SecretEPDB are provided in terms of protein characteristics, protein function, protein secondary structure, Pfam domains, metabolic pathway and evolutionary details. It is our hope that this integrated knowledgebase will serve as a useful resource for biological investigation and the generation of new hypotheses for research efforts aimed at bacterial secretion systems.

Comprehensive analysis of draft genomes of two closely related pseudomonas syringae phylogroup 2b strains infecting mono- and dicotyledon host plants

Background

In recent years, the damage caused by bacterial pathogens to major crops has been increasing worldwide. Pseudomonas syringae is a widespread bacterial species that infects almost all major crops. Different P. syringae strains use a wide range of biochemical mechanisms, including phytotoxins and effectors of the type III and type IV secretion systems, which determine the specific nature of the pathogen virulence.

Results

Strains 1845 (isolated from dicots) and 2507 (isolated from monocots) were selected for sequencing because they specialize on different groups of plants. We compared virulence factors in these and other available genomes of phylogroup 2 to find genes responsible for the specialization of bacteria. We showed that strain 1845 belongs to the clonal group that has been infecting monocots in Russia and USA for a long time (at least 50 years). Strain 1845 has relatively recently changed its host plant to dicots.

Conclusions

The results obtained by comparing the strain 1845 genome with the genomes of bacteria infecting monocots can help to identify the genes that define specific nature of the virulence of P. syringae strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3358-y) contains supplementary material, which is available to authorized users.

Draft Genome Sequence of Piscirickettsia litoralis, Isolated from Seawater

One species of Piscirickettsia, a pathogen of salmonid fish, has been described. The genome sequence of a putative second and free-living species may provide insights into the evolution of pathogenicity in the genus.

Classification of a Hypervirulent Aeromonas hydrophila Pathotype Responsible for Epidemic Outbreaks in Warm-Water Fishes

Lineages of hypervirulent Aeromonas hydrophila (vAh) are the cause of persistent outbreaks of motile Aeromonas septicemia in warm-water fishes worldwide. Over the last decade, this virulent lineage of A. hydrophila has resulted in annual losses of millions of tons of farmed carp and catfish in the People's Republic of China and the United States (US). Multiple lines of evidence indicate US catfish and Asian carp isolates of A. hydrophila affiliated with sequence type 251 (ST251) share a recent common ancestor. To address the genomic context for the putative intercontinental transfer and subsequent geographic spread of this pathogen, we conducted a core genome phylogenetic analysis on 61 Aeromonas spp. genomes, of which 40 were affiliated with A. hydrophila, with 26 identified as epidemic strains. Phylogenetic analyses indicate all ST251 strains form a coherent lineage affiliated with A. hydrophila. Within this lineage, conserved genetic loci unique to A. hydrophila were identified, with some genes present in consistently higher copy numbers than in non-epidemic A. hydrophila isolates. In addition, results from analyses of representative ST251 isolates support the conclusion that multiple lineages are present within US vAh isolated from Mississippi, whereas vAh isolated from Alabama appear clonal. This is the first report of genomic heterogeneity within US vAh isolates, with some Mississippi isolates showing closer affiliation with the Asian grass carp isolate ZC1 than other vAh isolated in the US. To evaluate the biological significance of the identified heterogeneity, comparative disease challenges were conducted with representatives of different vAh genotypes. These studies revealed that isolate ZC1 yielded significantly lower mortality in channel catfish, relative to Alabama and Mississippi vAh isolates. Like other Asian vAh isolates, the ZC1 lineage contains all core genes for a complete type VI secretion system (T6SS). In contrast, more virulent US isolates retain only remnants of the T6SS (clpB, hcp, vgrG, and vasH) which may have functional implications. Collectively, these results characterize a hypervirulent A. hydrophila pathotype that affects farmed fish on multiple continents.

Prediction of bacterial associations with plants using a supervised machine-learning approach

Recent scenarios of fresh produce contamination by human enteric pathogens have resulted in severe food-borne outbreaks, and a new paradigm has emerged stating that some human-associated bacteria can use plants as secondary hosts. As a consequence, there has been growing concern in the scientific community about these interactions that have not yet been elucidated. Since this is a relatively new area, there is a lack of strategies to address the problem of food-borne illnesses due to the ingestion of fruits and vegetables. In the present study, we performed specific genome annotations to train a supervised machine-learning model that allows for the identification of plant-associated bacteria with a precision of ∼93%. The application of our method to approximately 9500 genomes predicted several unknown interactions between well-known human pathogens and plants, and it also confirmed several cases for which evidence has been reported. We observed that factors involved in adhesion, the deconstruction of the plant cell wall and detoxifying activities were highlighted as the most predictive features. The application of our strategy to sequenced strains that are involved in food poisoning can be used as a primary screening tool to determine the possible causes of contaminations.

Host-Associated Genomic Features of the Novel Uncultured Intracellular Pathogen Ca. Ichthyocystis Revealed by Direct Sequencing of Epitheliocysts

Advances in single-cell and mini-metagenome sequencing have enabled important investigations into uncultured bacteria. In this study, we applied the mini-metagenome sequencing method to assemble genome drafts of the uncultured causative agents of epitheliocystis, an emerging infectious disease in the Mediterranean aquaculture species gilthead seabream. We sequenced multiple cyst samples and constructed 11 genome drafts from a novel beta-proteobacterial lineage, Candidatus Ichthyocystis. The draft genomes demonstrate features typical of pathogenic bacteria with an obligate intracellular lifestyle: a reduced genome of up to 2.6 Mb, reduced G + C content, and reduced metabolic capacity. Reconstruction of metabolic pathways reveals that Ca Ichthyocystis genomes lack all amino acid synthesis pathways, compelling them to scavenge from the fish host. All genomes encode type II, III, and IV secretion systems, a large repertoire of predicted effectors, and a type IV pilus. These are all considered to be virulence factors, required for adherence, invasion, and host manipulation. However, no evidence of lipopolysaccharide synthesis could be found. Beyond the core functions shared within the genus, alignments showed distinction into different species, characterized by alternative large gene families. These comprise up to a third of each genome, appear to have arisen through duplication and diversification, encode many effector proteins, and are seemingly critical for virulence. Thus, Ca Ichthyocystis represents a novel obligatory intracellular pathogenic beta-proteobacterial lineage. The methods used: mini-metagenome analysis and manual annotation, have generated important insights into the lifestyle and evolution of the novel, uncultured pathogens, elucidating many putative virulence factors including an unprecedented array of novel gene families.

Complete genome sequence of Vibrio alginolyticus ATCC 33787(T) isolated from seawater with three native megaplasmids

Vibrio alginolyticus, an opportunistic pathogen, is commonly associated with vibriosis in fish and shellfish and can also cause superficial and ear infections in humans. V. alginolyticus ATCC 33787(T) was originally isolated from seawater and has been used as one of the type strains for exploring the virulence factors of marine bacteria and for developing vaccine against vibriosis. Here we sequenced and assembled the whole genome of this strain, and identified three megaplasmids and three Type VI secretion systems, thus providing useful information for the study of virulence factors and for the development of vaccine for Vibrio.

The outer-membrane export signal of Porphyromonas gingivalis type IX secretion system (T9SS) is a conserved C-terminal β-sandwich domain

In the recently characterized Type IX Secretion System (T9SS), the conserved C-terminal domain (CTD) in secreted proteins functions as an outer membrane translocation signal for export of virulence factors to the cell surface in the Gram-negative Bacteroidetes phylum. In the periodontal pathogen Porphyromonas gingivalis, the CTD is cleaved off by PorU sortase in a sequence-independent manner, and anionic lipopolysaccharide (A-LPS) is attached to many translocated proteins, thus anchoring them to the bacterial surface. Here, we solved the atomic structure of the CTD of gingipain B (RgpB) from P. gingivalis, alone and together with a preceding immunoglobulin-superfamily domain (IgSF). The CTD was found to possess a typical Ig-like fold encompassing seven antiparallel β-strands organized in two β-sheets, packed into a β-sandwich structure that can spontaneously dimerise through C-terminal strand swapping. Small angle X-ray scattering (SAXS) revealed no fixed orientation of the CTD with respect to the IgSF. By introducing insertion or substitution of residues within the inter-domain linker in the native protein, we were able to show that despite the region being unstructured, it nevertheless is resistant to general proteolysis. These data suggest structural motifs located in the two adjacent Ig-like domains dictate the processing of CTDs by the T9SS secretion pathway.

Identification of protein secretion systems in bacterial genomes

Bacteria with two cell membranes (diderms) have evolved complex systems for protein secretion. These systems were extensively studied in some model bacteria, but the characterisation of their diversity has lagged behind due to lack of standard annotation tools. We built online and standalone computational tools to accurately predict protein secretion systems and related appendages in bacteria with LPS-containing outer membranes. They consist of models describing the systems’ components and genetic organization to be used with MacSyFinder to search for T1SS-T6SS, T9SS, flagella, Type IV pili and Tad pili. We identified ~10,000 candidate systems in bacterial genomes, where T1SS and T5SS were by far the most abundant and widespread. All these data are made available in a public database. The recently described T6SSⁱⁱⁱ and T9SS were restricted to Bacteroidetes, and T6SSⁱⁱ to Francisella. The T2SS, T3SS, and T4SS were frequently encoded in single-copy in one locus, whereas most T1SS were encoded in two loci. The secretion systems of diderm Firmicutes were similar to those found in other diderms. Novel systems may remain to be discovered, since some clades of environmental bacteria lacked all known protein secretion systems. Our models can be fully customized, which should facilitate the identification of novel systems.

Draft genome sequence for virulent and avirulent strains of Xanthomonas arboricola isolated from Prunus spp. in Spain

Xanthomonas arboricola is a species in genus Xanthomonas which is mainly comprised of plant pathogens. Among the members of this taxon, X. arboricola pv. pruni, the causal agent of bacterial spot disease of stone fruits and almond, is distributed worldwide although it is considered a quarantine pathogen in the European Union. Herein, we report the draft genome sequence, the classification, the annotation and the sequence analyses of a virulent strain, IVIA 2626.1, and an avirulent strain, CITA 44, of X. arboricola associated with Prunus spp. The draft genome sequence of IVIA 2626.1 consists of 5,027,671 bp, 4,720 protein coding genes and 50 RNA encoding genes. The draft genome sequence of strain CITA 44 consists of 4,760,482 bp, 4,250 protein coding genes and 56 RNA coding genes. Initial comparative analyses reveals differences in the presence of structural and regulatory components of the type IV pilus, the type III secretion system, the type III effectors as well as variations in the number of the type IV secretion systems. The genome sequence data for these strains will facilitate the development of molecular diagnostics protocols that differentiate virulent and avirulent strains. In addition, comparative genome analysis will provide insights into the plant-pathogen interaction during the bacterial spot disease process.

EffectiveDB--updates and novel features for a better annotation of bacterial secreted proteins and Type III, IV, VI secretion systems

Protein secretion systems play a key role in the interaction of bacteria and hosts. EffectiveDB (http://effectivedb.org) contains pre-calculated predictions of bacterial secreted proteins and of intact secretion systems. Here we describe a major update of the database, which was previously featured in the NAR Database Issue. EffectiveDB bundles various tools to recognize Type III secretion signals, conserved binding sites of Type III chaperones, Type IV secretion peptides, eukaryotic-like domains and subcellular targeting signals in the host. Beyond the analysis of arbitrary protein sequence collections, the new release of EffectiveDB also provides a ‘genome-mode’, in which protein sequences from nearly complete genomes or metagenomic bins can be screened for the presence of three important secretion systems (Type III, IV, VI). EffectiveDB contains pre-calculated predictions for currently 1677 bacterial genomes from the EggNOG 4.0 database and for additional bacterial genomes from NCBI RefSeq. The new, user-friendly and informative web portal offers a submission tool for running the EffectiveDB prediction tools on user-provided data.

Bioinformatics Analysis of the Complete Genome Sequence of the Mango Tree Pathogen Pseudomonas syringae pv. syringae UMAF0158 Reveals Traits Relevant to Virulence and Epiphytic Lifestyle

The genome sequence of more than 100 Pseudomonas syringae strains has been sequenced to date; however only few of them have been fully assembled, including P. syringae pv. syringae B728a. Different strains of pv. syringae cause different diseases and have different host specificities; so, UMAF0158 is a P. syringae pv. syringae strain related to B728a but instead of being a bean pathogen it causes apical necrosis of mango trees, and the two strains belong to different phylotypes of pv.syringae and clades of P. syringae. In this study we report the complete sequence and annotation of P. syringae pv. syringae UMAF0158 chromosome and plasmid pPSS158. A comparative analysis with the available sequenced genomes of other 25 P. syringae strains, both closed (the reference genomes DC3000, 1448A and B728a) and draft genomes was performed. The 5.8 Mb UMAF0158 chromosome has 59.3% GC content and comprises 5017 predicted protein-coding genes. Bioinformatics analysis revealed the presence of genes potentially implicated in the virulence and epiphytic fitness of this strain. We identified several genetic features, which are absent in B728a, that may explain the ability of UMAF0158 to colonize and infect mango trees: the mangotoxin biosynthetic operon mbo, a gene cluster for cellulose production, two different type III and two type VI secretion systems, and a particular T3SS effector repertoire. A mutant strain defective in the rhizobial-like T3SS Rhc showed no differences compared to wild-type during its interaction with host and non-host plants and worms. Here we report the first complete sequence of the chromosome of a pv. syringae strain pathogenic to a woody plant host. Our data also shed light on the genetic factors that possibly determine the pathogenic and epiphytic lifestyle of UMAF0158. This work provides the basis for further analysis on specific mechanisms that enable this strain to infect woody plants and for the functional analysis of host specificity in the P. syringae complex.