Analysis of the hypervariable region of the Salmonella enterica genome associated with tRNA(leuX)

Structure classification of the proteins from Salmonella enterica pangenome revealed novel potential pathogenicity islands

Salmonella enterica is a pathogenic bacterium known for causing severe typhoid fever in humans, making it important to study due to its potential health risks and significant impact on public health. This study provides evolutionary classification of proteins from Salmonella enterica pangenome. We classified 17,238 domains from 13,147 proteins from 79,758 Salmonella enterica strains and studied in detail domains of 272 proteins from 14 characterized Salmonella pathogenicity islands (SPIs). Among SPIs-related proteins, 90 proteins function in the secretion machinery. 41% domains of SPI proteins have no previous sequence annotation. By comparing clinical and environmental isolates, we identified 3682 proteins that are overrepresented in clinical group that we consider as potentially pathogenic. Among domains of potentially pathogenic proteins only 50% domains were annotated by sequence methods previously. Moreover, 36% (1330 out of 3682) of potentially pathogenic proteins cannot be classified into Evolutionary Classification of Protein Domains database (ECOD). Among classified domains of potentially pathogenic proteins the most populated homology groups include helix-turn-helix (HTH), Immunoglobulin-related, and P-loop domains-related. Functional analysis revealed overrepresentation of these protein in biological processes related to viral entry into host cell, antibiotic biosynthesis, DNA metabolism and conformation change, and underrepresentation in translational processes. Analysis of the potentially pathogenic proteins indicates that they form 119 clusters or novel potential pathogenicity islands (NPPIs) within the Salmonella genome, suggesting their potential contribution to the bacterium's virulence. One of the NPPIs revealed significant overrepresentation of potentially pathogenic proteins. Overall, our analysis revealed that identified potentially pathogenic proteins are poorly studied.

The DNA Phosphorothioation Restriction-Modification System Influences the Antimicrobial Resistance of Pathogenic Bacteria

Bacterial defense barriers, such as DNA methylation-associated restriction-modification (R-M) and the CRISPR-Cas system, play an important role in bacterial antimicrobial resistance (AMR). Recently, a novel R-M system based on DNA phosphorothioate (PT) modification has been shown to be widespread in the kingdom of Bacteria as well as Archaea. However, the potential role of the PT R-M system in bacterial AMR remains unclear. In this study, we explored the role of PT R-Ms in AMR with a series of common clinical pathogenic bacteria. By analyzing the distribution of AMR genes related to mobile genetic elements (MGEs), it was shown that the presence of PT R-M effectively reduced the distribution of horizontal gene transfer (HGT)-derived AMR genes in the genome, even in the bacteria that did not tend to acquire AMR genes by HGT. In addition, unique gene variation analysis based on pangenome analysis and MGE prediction revealed that the presence of PT R-M could suppress HGT frequency. Thus, this is the first report showing that the PT R-M system has the potential to repress HGT-derived AMR gene acquisition by reducing the HGT frequency. IMPORTANCE In this study, we demonstrated the effect of DNA PT modification-based R-M systems on horizontal gene transfer of AMR genes in pathogenic bacteria. We show that there is no apparent association between the genetic background of the strains harboring PT R-Ms and the number of AMR genes or the kinds of gene families. The strains equipped with PT R-M harbor fewer plasmid-derived, prophage-derived, or integrating mobile genetic element (iMGE)-related AMR genes and have a lower HGT frequency, but the degree of inhibition varies among different bacteria. In addition, compared with Salmonella enterica and Escherichia coli, Klebsiella pneumoniae prefers to acquire MGE-derived AMR genes, and there is no coevolution between PT R-M clusters and bacterial core genes.

Reassembling a cannon in the DNA defense arsenal: Genetics of StySA, a BREX phage exclusion system in Salmonella lab strains

Understanding mechanisms that shape horizontal exchange in prokaryotes is a key problem in biology. A major limit on DNA entry is imposed by restriction-modification (RM) processes that depend on the pattern of DNA modification at host-specified sites. In classical RM, endonucleolytic DNA cleavage follows detection of unprotected sites on entering DNA. Recent investigation has uncovered BREX (BacteRiophage EXclusion) systems. These RM-like activities employ host protection by DNA modification, but immediate replication arrest occurs without evident of nuclease action on unmodified phage DNA. Here we show that the historical stySA RM locus of Salmonella enterica sv Typhimurium is a variant BREX system. A laboratory strain disabled for both the restriction and methylation activity of StySA nevertheless has wild type sequence in pglX, the modification gene homolog. Instead, flanking genes pglZ and brxC each carry multiple mutations (μ) in their C-terminal domains. We further investigate this system in situ, replacing the mutated pglZμ and brxCμ genes with the WT counterpart. PglZ-WT supports methylation in the presence of either BrxCμ or BrxC-WT but not in the presence of a deletion/insertion allele, ΔbrxC::cat. Restriction requires both BrxC-WT and PglZ-WT, implicating the BrxC C-terminus specifically in restriction activity. These results suggests that while BrxC, PglZ and PglX are principal components of the BREX modification activity, BrxL is required for restriction only. Furthermore, we show that a partial disruption of brxL disrupts transcription globally.

Horizontal gene transfer is a major driver of evolution and adaptation in bacteria. Genes from outside may be beneficial or dangerous to the receiving cell. Benefits include new food sources such as sugars, or new homes by adhesion, or new resistances, as to antibiotics. Dangers are posed by bacteriophages--viruses that take over the cell machinery, multiply, and release progeny to kill sister cells. Host-dependent restriction-modification systems enable defense that distinguishes relatives from strangers: using a modification pattern (M) carried by DNA bases added by the host cell to prevent restriction (R). Sisters and cousin cells will have the same protective pattern on DNA, while DNA of foreign origin will have the wrong M pattern and be restricted (R, rejected). Typically, restriction involves nuclease digestion. Here we address the enigmatic StySA RM system, one of the earliest to be genetically characterized. It is a variant of the newly recognized defense mechanism, BREX. BREX systems also track DNA history via modification pattern, but restrict by a novel, uncharacterized mechanism. Like other BREX family systems, StySA-BREX modification requires multiple components. When StySA-BREX transcription is unbalanced, we find global disruption of gene transcription. The disruption pattern does not suggest SOS-inducing damage to DNA.

The Invasin and Complement-Resistance Protein Rck of Salmonella is More Widely Distributed than Previously Expected

The rck open reading frame (ORF) on the pefI-srgC operon encodes an outer membrane protein responsible for invasion of nonphagocytic cell lines and resistance to complement-mediated killing. Until now, the rck ORF was only detected on the virulence plasmids of three serovars of Salmonella subsp. enterica (i.e., Bovismorbificans, Enteritidis, and Typhimurium). The increasing number of Salmonella genome sequences allowed us to use a combination of reference sequences and whole-genome multilocus sequence typing (wgMLST) data analysis to probe the presence of the operon and of rck in a wide array of isolates belonging to all Salmonella species and subspecies. We established the presence of partial or complete operons in 61 subsp. enterica serovars as well as in 4 other subspecies with various syntenies and frequencies. The rck ORF itself was retrieved in 36 subsp. enterica serovars and in two subspecies with either chromosomal or plasmid-borne localization. It displays high conservation of its sequence within the genus, and we demonstrated that most of the allelic variations identified did not alter the virulence properties of the protein. However, we demonstrated the importance of the residue at position 38 (at the level of the first extracellular loop of the protein) in the invasin function of Rck. Altogether, our results highlight that rck is not restricted to the three formerly identified serovars and could therefore have a more important role in virulence than previously expected. Moreover, this work raises questions about the mechanisms involved in rck acquisition and about virulence plasmid distribution and evolution.

IMPORTANCE The foodborne pathogen Salmonella is responsible for a wide variety of pathologies depending on the infected host, the infecting serovars, and its set of virulence factors. However, the implication of each of these virulence factors and their role in the specific host-pathogen interplay are not fully understood. The significance of our research is in determining the distribution of one of these factors, the virulence plasmid-encoded invasin and resistance to complement killing protein Rck. In addition to providing elements of reflection concerning the mechanisms of acquisition of specific virulence genes in certain serotypes, this work will help to understand the role of Rck in the pathogenesis of Salmonella.

Embracing Diversity: Differences in Virulence Mechanisms, Disease Severity, and Host Adaptations Contribute to the Success of Nontyphoidal Salmonella as a Foodborne Pathogen

Not all Salmonella enterica serovars cause the same disease. S. enterica represents an incredibly diverse species comprising >2,600 unique serovars. While some S. enterica serovars are host-restricted, others infect a wide range of hosts. The diseases that nontyphoidal Salmonella (NTS) serovars cause vary considerably, with some serovars being significantly more likely to cause invasive disease in humans than others. Furthermore, while genomic analyses have advanced our understanding of the genetic diversity of these serovars, they have not been able to fully account for the observed clinical differences. One overarching challenge is that much of what is known about Salmonella's general biology and virulence strategies is concluded from studies examining a select few serovars, especially serovar Typhimurium. As targeted control strategies have been implemented to control select serovars, an increasing number of foodborne outbreaks involving serovars that are less frequently associated with human clinical illness are being detected. Harnessing what is known about the diversity of NTS serovars represents an important factor in achieving the ultimate goal of reducing salmonellosis-associated morbidity and mortality worldwide. In this review we summarize the current understanding of the differences and similarities among NTS serovars, highlighting the virulence mechanisms, genetic differences, and sources that characterize S. enterica diversity and contribute to its success as a foodborne pathogen.

In silico Identification of Serovar-Specific Genes for Salmonella Serotyping

Salmonella enterica subspecies enterica is a highly diverse subspecies with more than 1500 serovars and the ability to distinguish serovars within this group is vital for surveillance. With the development of whole-genome sequencing technology, serovar prediction by traditional serotyping is being replaced by molecular serotyping. Existing in silico serovar prediction approaches utilize surface antigen encoding genes, core genome MLST and serovar-specific gene markers or DNA fragments for serotyping. However, these serovar-specific gene markers or DNA fragments only distinguished a small number of serovars. In this study, we compared 2258 Salmonella accessory genomes to identify 414 candidate serovar-specific or lineage-specific gene markers for 106 serovars which includes 24 polyphyletic serovars and the paraphyletic serovar Enteritidis. A combination of several lineage-specific gene markers can be used for the clear identification of the polyphyletic serovars and the paraphyletic serovar. We designed and evaluated an in silico serovar prediction approach by screening 1089 genomes representing 106 serovars against a set of 131 serovar-specific gene markers. The presence or absence of one or more serovar-specific gene markers was used to predict the serovar of an isolate from genomic data. We show that serovar-specific gene markers have comparable accuracy to other in silico serotyping methods with 84.8% of isolates assigned to the correct serovar with no false positives (FP) and false negatives (FN) and 10.5% of isolates assigned to a small subset of serovars containing the correct serovar with varied FP. Combined, 95.3% of genomes were correctly assigned to a serovar. This approach would be useful as diagnosis moves to culture-independent and metagenomic methods as well as providing a third alternative to confirm other genome-based analyses. The identification of a set of gene markers may also be useful in the development of more cost-effective molecular assays designed to detect specific gene markers of the all major serovars in a region. These assays would be useful in serotyping isolates where cultures are no longer obtained and traditional serotyping is therefore impossible.

Virulotyping of Salmonella enterica serovar Typhi isolates from Pakistan: Absence of complete SPI-10 in Vi negative isolates

The pathogenesis of Salmonella enterica serovar Typhi (S. Typhi), the cause of typhoid fever in humans, is mainly attributed to the acquisition of horizontally acquired DNA elements. Salmonella pathogenicity islands (SPIs) are indubitably the most important form of horizontally acquired DNA with respect to pathogenesis of this bacterium. The insertion or deletion of any of these transferrable SPIs may have impact on the virulence potential of S. Typhi. In this study, the virulence potential and genetic relatedness of 35 S. Typhi isolates, collected from 2004 to 2013 was determined by identification of SPI and non-SPI virulence factors through a combination of techniques including virulotyping, Whole Genome Sequencing (WGS), and Variable Number of Tandem Repeats (VNTR) profiling. In order to determine the virulence potential of local S. Typhi isolates, 56 virulence related genes were studied by PCR. These genes are located in the core as well as accessory genome (SPIs and plasmid). Major variations among studied virulence determinants were found in case of SPI-7 and SPI-10 associated genes. On the basis of presence of virulence related genes, the studied S. Typhi isolates from Pakistan were clustered into two virulotypes Vi-positive and Vi-negative. Interestingly, SPI-7 and SPI-10 were collectively absent or present in Vi-negative and Vi-positive strains, respectively. Two Vi-negative and 11 Vi-positive S. Typhi strains were also analyzed by whole genome sequencing (WGS) and their results supported the PCR results. Genetic diversity was tested by VNTR-based molecular typing. All 35 isolates were clustered into five groups. Overall, all Vi-negative isolates were placed in a single group (T5) whereas Vi-positive isolates were grouped into four types. Vi-negative and Vi-positive isolates were mutually exclusive. This is the first report on the comparative distribution of SPI and non-SPI related virulence genes in Vi-negative and Vi-positive S. Typhi isolates with an important finding that SPI-10 is absent in all Vi-negative isolates.

The distribution of virulence factors in S. Typhi can vary in isolates from different geographical regions and can have significant effect on the disease control. In this study, we have checked the distribution of 56 reported virulence associated factors in 35 local isolates of S. Typhi to identify any variations that can help in designing effective control strategies for typhoid. We have identified four naturally occurring variants which are simultaneously lacking SPI-7 and SPI-10, two adjacently located pathogenicity islands on S. Typhi chromosome. These isolates are not producing Vi capsular antigen hence the Vi based vaccines will not be effective against them. These findings highlight the need to develop typhoid vaccines specifically effective in Pakistan.

The Homolog of the Gene bstA of the BTP1 Phage from Salmonella enterica Serovar Typhimurium ST313 Is an Antivirulence Gene in Salmonella enterica Serovar Dublin

In a previous study, a novel virulence gene, bstA, identified in a Salmonella enterica serovar Typhimurium sequence type 313 (ST313) strain was found to be conserved in all published Salmonella enterica serovar Dublin genomes. In order to analyze the role of this gene in the host-pathogen interaction in S Dublin, a mutant where this gene was deleted (S Dublin ΔbstA) and a mutant which was further genetically complemented with bstA (S Dublin 3246-C) were constructed and tested in models of in vitro and in vivo infection as well as during growth competition assays in M9 medium, Luria-Bertani broth, and cattle blood. In contrast to the results obtained for a strain of S Typhimurium ST313, the lack of bstA was found to be associated with increased virulence in S Dublin. Thus, S Dublin ΔbstA showed higher levels of uptake than the wild-type strain during infection of mouse and cattle macrophages and higher net replication within human THP-1 cells. Furthermore, during mouse infections, S Dublin ΔbstA was more virulent than the wild type following a single intraperitoneal infection and showed an increased competitive index during competitive infection assays. Deletion of bstA did not affect either the amount of cytokines released by THP-1 macrophages or the cytotoxicity toward these cells. The histology of the livers and spleens of mice infected with the wild-type strain and the S Dublin ΔbstA mutant revealed similar levels of inflammation between the two groups. The gene was not important for adherence to or invasion of human epithelial cells and did not influence bacterial growth in rich medium, minimal medium, or cattle blood. In conclusion, a lack of bstA affects the pathogenicity of S Dublin by decreasing its virulence. Therefore, it might be regarded as an antivirulence gene in this serovar.

Tet(C) Gene Transfer between Chlamydia suis Strains Occurs by Homologous Recombination after Co-infection: Implications for Spread of Tetracycline-Resistance among Chlamydiaceae

Chlamydia suis is a swine pathogen that has also recently been found to cause zoonotic infections of the human eye, pharynx, and gastrointestinal tract. Many strains contain a tetracycline class C gene [tet(C)] cassette that confers tetracycline resistance. The cassette was likely originally acquired by horizontal gene transfer from a Gram-negative donor after the introduction of tetracycline into animal feed in the 1950s. Various research groups have described the capacity for different Chlamydia species to exchange DNA by homologous recombination. Since over 90% of C. suis strains are tetracycline resistant, they represent a potential source for antibiotic-resistance spread within and between Chlamydiaceae species. Here, we examined the genetics of tet(C)-transfer among C. suis strains. Tetracycline-sensitive C. suis strain S45 was simultaneously or sequentially co-infected with tetracycline-resistant C. suis strains in McCoy cells. Potential recombinants were clonally purified by a harvest assay derived from the classic plaque assay. C. suis strain Rogers132, lacking transposases IS200 and IS605, was the most efficient donor, producing two unique recombinants detected in three of the 56 (5.4%) clones screened. Recombinants were found to have a minimal inhibitory concentration (MIC) of 8-16 μg/mL for tetracycline. Resistance remained stable over 10 passages as long as recombinants were initially grown in tetracycline at twice the MIC of S45 (0.032 μg/mL). Genomic analysis revealed that tet(C) had integrated into the S45 genome by homologous recombination at two unique sites depending on the recombinant: a 55 kb exchange between nrqF and pckG, and a 175 kb exchange between kdsA and cysQ. Neither site was associated with inverted repeats or motifs associated with recombination hotspots. Our findings show that cassette transfer into S45 has low frequency, does not require IS200/IS605 transposases, is stable if initially grown in tetracycline, and results in multiple genomic configurations. We provide a model for stable cassette transfer to better understand the capability for cassette acquisition by Chlamydiaceae species that infect humans, a matter of public health importance.

Two Inducible Prophages of an Antarctic Pseudomonas sp. ANT_H14 Use the Same Capsid for Packaging Their Genomes - Characterization of a Novel Phage Helper-Satellite System

Two novel prophages ФAH14a and ФAH14b of a psychrotolerant Antarctic bacterium Pseudomonas sp. ANT_H14 have been characterized. They were simultaneously induced with mitomycin C and packed into capsids of the same size and protein composition. The genome sequences of ФAH14a and ФAH14b have been determined. ФAH14b, the phage with a smaller genome (16,812 bp) seems to parasitize ФAH14a (55,060 bp) and utilizes its capsids, as only the latter encodes a complete set of structural proteins. Both viruses probably constitute a phage helper-satellite system, analogous to the P2-P4 duo. This study describes the architecture and function of the ФAH14a and ФAH14b genomes. Moreover, a functional analysis of a ФAH14a-encoded lytic enzyme and a DNA methyltransferase was performed. In silico analysis revealed the presence of the homologs of ФAH14a and ФAH14b in other Pseudomonas genomes, which may suggest that helper-satellite systems related to the one described in this work are common in pseudomonads.

Novel genomic island modifies DNA with 7-deazaguanine derivatives

The discovery of ∼20-kb gene clusters containing a family of paralogs of tRNA guanosine transglycosylase genes, called tgtA5, alongside 7-cyano-7-deazaguanine (preQ0) synthesis and DNA metabolism genes, led to the hypothesis that 7-deazaguanine derivatives are inserted in DNA. This was established by detecting 2'-deoxy-preQ0 and 2'-deoxy-7-amido-7-deazaguanosine in enzymatic hydrolysates of DNA extracted from the pathogenic, Gram-negative bacteria Salmonella enterica serovar Montevideo. These modifications were absent in the closely related S. enterica serovar Typhimurium LT2 and from a mutant of S Montevideo, each lacking the gene cluster. This led us to rename the genes of the S. Montevideo cluster as dpdA-K for 7-deazapurine in DNA. Similar gene clusters were analyzed in ∼150 phylogenetically diverse bacteria, and the modifications were detected in DNA from other organisms containing these clusters, including Kineococcus radiotolerans, Comamonas testosteroni, and Sphingopyxis alaskensis Comparative genomic analysis shows that, in Enterobacteriaceae, the cluster is a genomic island integrated at the leuX locus, and the phylogenetic analysis of the TgtA5 family is consistent with widespread horizontal gene transfer. Comparison of transformation efficiencies of modified or unmodified plasmids into isogenic S. Montevideo strains containing or lacking the cluster strongly suggests a restriction-modification role for the cluster in Enterobacteriaceae. Another preQ0 derivative, 2'-deoxy-7-formamidino-7-deazaguanosine, was found in the Escherichia coli bacteriophage 9 g, as predicted from the presence of homologs of genes involved in the synthesis of the archaeosine tRNA modification. These results illustrate a deep and unexpected evolutionary connection between DNA and tRNA metabolism.

Genomic characterization of Salmonella Cerro ST367, an emerging Salmonella subtype in cattle in the United States

Background

Within the last decade, Salmonella enterica subsp. enterica serovar Cerro (S. Cerro) has become one of the most common serovars isolated from cattle and dairy farm environments in the northeastern US. The fact that this serovar is commonly isolated from subclinically infected cattle and is rarely associated with human disease, despite its frequent isolation from cattle, has led to the hypothesis that this emerging serovar may be characterized by reduced virulence. We applied comparative and population genomic approaches to (i) characterize the evolution of this recently emerged serovar and to (ii) gain a better understanding of genomic features that could explain some of the unique epidemiological features associated with this serovar.

Results

In addition to generating a de novo draft genome for one Salmonella Cerro strain, we also generated whole genome sequence data for 26 additional S. Cerro isolates, including 16 from cattle operations in New York (NY) state, 2 from human clinical cases from NY in 2008, and 8 from diverse animal sources (7 from Washington state and 1 from Florida). All isolates sequenced in this study represent sequence type ST367. Population genomic analysis showed that isolates from the NY cattle operations form a well-supported clade within S. Cerro ST367 (designated here “NY bovine clade”), distinct from isolates from Washington state, Florida and the human clinical cases. A molecular clock analysis indicates that the most recent common ancestor of the NY bovine clade dates back to 1998, supporting the recent emergence of this clone.

Comparative genomic analyses revealed several relevant genomic features of S. Cerro ST367, that may be responsible for reduced virulence of S. Cerro, including an insertion creating a premature stop codon in sopA. In addition, patterns of gene deletion in S. Cerro ST367 further support adaptation of this clone to a unique ecological or host related niche.

Conclusions

Our results indicate that the increase in prevalence of S. Cerro ST367 is caused by a highly clonal subpopulation and that S. Cerro ST367 is characterized by unique genomic deletions that may indicate adaptation to specific ecological niches and possibly reduced virulence in some hosts.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-427) contains supplementary material, which is available to authorized users.

Evidence for horizontal gene transfer, gene duplication and genetic variation as driving forces of the diversity of haemolytic phenotypes in Photobacterium damselae subsp. damselae

Photobacterium damselae subsp. damselae, a marine bacterium that causes infections in marine animals and in humans, produces up to three different haemolysins involved in virulence, which include the pPHDD1 plasmid-encoded damselysin (Dly) and HlyApl , and the chromosome-encoded HlyAch . We screened 45 isolates from different origins, and found a correlation between their haemolytic phenotypes and the differential haemolysin gene content. All highly and medium haemolytic strains harboured pPHDD1, with amino acid substitutions in HlyApl and HlyAch being the cause of the medium haemolytic phenotypes in some pPHDD1-harbouring strains. Weakly haemolytic strains contained only hlyAch , whereas nonhaemolytic isolates, in addition to lacking pPHDD1, either lacked hlyAch or contained a hlyAch pseudogene. Sequence analysis of the genomic context of hlyAch uncovered an unexpected genetic diversity, suggesting that hlyAch is located in an unstable chromosomal region. Phylogenetic analysis suggested that hlyApl and hlyAch originated by gene duplication within P. damselae subsp. damselae following acquisition by horizontal transfer. These observations together with the differential distribution of pPHDD1 plasmid among strains suggest that horizontal gene transfer has played a main role in shaping the haemolysin gene baggage in this pathogen.

Genome analysis and CRISPR typing of Salmonella enterica serovar Virchow

Background

Salmonella enterica subsp. enterica serovar Virchow has been recognized as a significant health burden in Asia, Australia and Europe. In addition to its global distribution, S. Virchow is clinically significant due to the frequency at which it causes invasive infections and its association with outbreaks arising from food-borne transmission. Here, we examine the genome of an invasive isolate of S. Virchow SVQ1 (phage type 8) from an outbreak in southeast Queensland, Australia. In addition to identifying new potential genotyping targets that could be used for discriminating between S. Virchow strains in outbreak scenarios, we also aimed to carry out a comprehensive comparative analysis of the S. Virchow genomes.

Results

Genome comparisons between S. Virchow SVQ1 and S. Virchow SL491, a previously published strain, identified a high degree of genomic similarity between the two strains with fewer than 200 single nucleotide differences. Clustered Regularly Interspaced Palindromic Repeats (CRISPR) regions were identified as a highly variable region that could be used to discriminate between S. Virchow isolates. We amplified and sequenced the CRISPR regions of fifteen S. Virchow isolates collected from seven different outbreaks across Australia. We observed three allelic types of the CRISPR region from these isolates based on the presence/absence of the spacers and were able to discriminate S. Virchow phage type 8 isolates originating from different outbreaks. A comparison with 27 published Salmonella genomes found that the S. Virchow SVQ1 genome encodes 11 previously described Salmonella Pathogenicity Islands (SPI), as well as additional genomic islands including a remnant integrative conjugative element that is distinct from SPI-7. In addition, the S. Virchow genome possesses a novel prophage that encodes the Type III secretion system effector protein SopE, a key Salmonella virulence factor. The prophage shares very little similarity to the SopE prophages found in other Salmonella serovars suggesting an independent acquisition of sopE.

Conclusions

The availability of this genome will serve as a genome template and facilitate further studies on understanding the virulence and global distribution of the S. Virchow serovar, as well as the development of genotyping methods for outbreak investigations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-389) contains supplementary material, which is available to authorized users.

Identification and characterization of novel Salmonella mobile elements involved in the dissemination of genes linked to virulence and transmission

The genetic diversity represented by >2,500 different Salmonella serovars provides a yet largely uncharacterized reservoir of mobile elements that can contribute to the frequent emergence of new pathogenic strains of this important zoonotic pathogen. Currently, our understanding of Salmonella mobile elements is skewed by the fact that most studies have focused on highly virulent or common serovars. To gain a more global picture of mobile elements in Salmonella, we used prediction algorithms to screen for mobile elements in 16 sequenced Salmonella genomes representing serovars for which no prior genome scale mobile element data were available. From these results, selected mobile elements underwent further analyses in the form of validation studies, comparative analyses, and PCR-based population screens. Through this analysis we identified a novel plasmid that has two cointegrated replicons (IncI1-IncFIB); this plasmid type was found in four genomes representing different Salmonella serovars and contained a virulence gene array that had not been previously identified. A Salmonella Montevideo isolate contained an IncHI and an IncN2 plasmid, which both encoded antimicrobial resistance genes. We also identified two novel genomic islands (SGI2 and SGI3), and 42 prophages with mosaic architecture, seven of them harboring known virulence genes. Finally, we identified a novel integrative conjugative element (ICE) encoding a type IVb pilus operon in three non-typhoidal Salmonella serovars. Our analyses not only identified a considerable number of mobile elements that have not been previously reported in Salmonella, but also found evidence that these elements facilitate transfer of genes that were previously thought to be limited in their distribution among Salmonella serovars. The abundance of mobile elements encoding pathogenic properties may facilitate the emergence of strains with novel combinations of pathogenic traits.

So similar, yet so different: uncovering distinctive features in the genomes of Salmonella enterica serovars Typhimurium and Typhi

Salmonella enterica represents a major human and animal pathogen. Many S. enterica genomes have been completed and many more genome sequencing projects are underway, constituting an excellent resource for comparative genome analysis studies leading to a better understanding of bacterial evolution and pathogenesis. Salmonella enterica serovar Typhimurium and Typhi are the best-characterized serovars, with the first being involved in localized gastroenteritis in many hosts and the latter causing a systemic human-specific disease. Here, we summarize the major genetic differences between the two different serovars. We detail the divergent repertoires of the virulence factors responsible for the pathogenesis of the organisms and that ultimately result in the distinct clinical outcomes of infection. This comparative genomic overview highlights hypotheses for future investigations on S. enterica pathogenesis and the basis of host specificity.

Viruses are essential agents within the roots and stem of the tree of life

In contrast with former definitions of life limited to membrane-bound cellular life forms which feed, grow, metabolise and replicate (i) a role of viruses as genetic symbionts, (ii) along with peripheral phenomena such as cryptobiosis and (iii) the horizontal nature of genetic information acquisition and processing broaden our view of the tree of life. Some researchers insist on the traditional textbook conviction of what is part of the community of life. In a recent review [Moreira, D., Lopez-Garcia, P., 2009. Ten reasons to exclude viruses from the tree of life. Nat. Rev. Microbiol. 7, 306-311.] they assemble four main arguments which should exclude viruses from the tree of life because of their inability to self-sustain and self-replicate, their polyphyly, the cellular origin of their cell-like genes and the volatility of their genomes. In this article we will show that these features are not coherent with current knowledge about viruses but that viral agents play key roles within the roots and stem of the tree of life.

Complete genome sequence of the chemolithoautotrophic marine magnetotactic coccus strain MC-1

The marine bacterium strain MC-1 is a member of the alpha subgroup of the proteobacteria that contains the magnetotactic cocci and was the first member of this group to be cultured axenically. The magnetotactic cocci are not closely related to any other known alphaproteobacteria and are only distantly related to other magnetotactic bacteria. The genome of MC-1 contains an extensive (102 kb) magnetosome island that includes numerous genes that are conserved among all known magnetotactic bacteria, as well as some genes that are unique. Interestingly, certain genes that encode proteins considered to be important in magnetosome assembly (mamJ and mamW) are absent from the genome of MC-1. Magnetotactic cocci exhibit polar magneto-aerotaxis, and the MC-1 genome contains a relatively large number of identified chemotaxis genes. Although MC-1 is capable of both autotrophic and heterotrophic growth, it does not appear to be metabolically versatile, with heterotrophic growth confined to the utilization of acetate. Central carbon metabolism is encoded by genes for the citric acid cycle (oxidative and reductive), glycolysis, and gluconeogenesis. The genome also reveals the presence or absence of specific genes involved in the nitrogen, sulfur, iron, and phosphate metabolism of MC-1, allowing us to infer the presence or absence of specific biochemical pathways in strain MC-1. The pathways inferred from the MC-1 genome provide important information regarding central metabolism in this strain that could provide insights useful for the isolation and cultivation of new magnetotactic bacterial strains, in particular strains of other magnetotactic cocci.

Identification of Salmonella enterica serovar Dublin-specific sequences by subtractive hybridization and analysis of their role in intestinal colonization and systemic translocation in cattle

Salmonella enterica serovar Dublin is a host-restricted serovar associated with typhoidal disease in cattle. In contrast, the fowl-associated serovar S. enterica serovar Gallinarum is avirulent in calves, yet it invades ileal mucosa and induces enteritis at levels comparable to those induced by S. enterica serovar Dublin. Suppression subtractive hybridization was employed to identify S. enterica serovar Dublin strain SD3246 genes absent from S. enterica serovar Gallinarum strain SG9. Forty-one S. enterica serovar Dublin fragments were cloned and sequenced. Among these, 24 mobile-element-associated genes were identified, and 12 clones exhibited similarity with sequences of known or predicted function in other serovars. Three S. enterica serovar Dublin-specific regions were homologous to regions from the genome of Enterobacter sp. strain 638. Sequencing of fragments adjacent to these three sequences revealed the presence of a 21-kb genomic island, designated S. enterica serovar Dublin island 1 (SDI-1). PCR analysis and Southern blotting showed that SDI-1 is highly conserved within S. enterica serovar Dublin isolates but rarely found in other serovars. To probe the role of genes identified by subtractive hybridization in vivo, 24 signature-tagged S. enterica serovar Dublin SD3246 mutants lacking loci not present in Salmonella serovar Gallinarum SG9 were created and screened by oral challenge of cattle. Though attenuation of tagged SG9 and SD3246 Salmonella pathogenicity island-1 (SPI-1) and SPI-2 mutant strains was detected, no obvious defects of these 24 mutants were detected. Subsequently, a DeltaSDI-1 mutant was found to exhibit weak but significant attenuation compared with the parent strain in coinfection of calves. SDI-1 mutation did not impair invasion, intramacrophage survival, or virulence in mice, implying that SDI-1 does not influence fitness per se and may act in a host-specific manner.

Comparative genome analysis of Salmonella Enteritidis PT4 and Salmonella Gallinarum 287/91 provides insights into evolutionary and host adaptation pathways

We have determined the complete genome sequences of a host-promiscuous Salmonella enterica serovar Enteritidis PT4 isolate P125109 and a chicken-restricted Salmonella enterica serovar Gallinarum isolate 287/91. Genome comparisons between these and other Salmonella isolates indicate that S. Gallinarum 287/91 is a recently evolved descendent of S. Enteritidis. Significantly, the genome of S. Gallinarum has undergone extensive degradation through deletion and pseudogene formation. Comparison of the pseudogenes in S. Gallinarum with those identified previously in other host-adapted bacteria reveals the loss of many common functional traits and provides insights into possible mechanisms of host and tissue adaptation. We propose that experimental analysis in chickens and mice of S. Enteritidis-harboring mutations in functional homologs of the pseudogenes present in S. Gallinarum could provide an experimentally tractable route toward unraveling the genetic basis of host adaptation in S. enterica.

Genetic diversity and composition of a plasmid metagenome from a wastewater treatment plant

Plasmid metagenome nucleotide sequence data were recently obtained from wastewater treatment plant (WWTP) bacteria with reduced susceptibility to selected antimicrobial drugs by applying the ultrafast 454-sequencing technology. The sequence dataset comprising 36,071,493 bases (346,427 reads with an average read length of 104 bases) was analysed for genetic diversity and composition by using a newly developed bioinformatic pipeline based on assignment of environmental gene tags (EGTs) to protein families stored in the Pfam database. Short amino acid sequences deduced from the plasmid metagenome sequence reads were compared to profile hidden Markov models underlying Pfam. Obtained matches evidenced that many reads represent genes having predicted functions in plasmid replication, stability and plasmid mobility which indicates that WWTP bacteria harbour genetically stabilised and mobile plasmids. Moreover, the data confirm a high diversity of plasmids residing in WWTP bacteria. The mobile organic peroxide resistance plasmid pMAC from Acinetobacter baumannii was identified as reference plasmid for the most abundant replication module type in the sequenced sample. Accessory plasmid modules encode different transposons, insertion sequences, integrons, resistance and virulence determinants. Most of the matches to Transposase protein families were identified for transposases similar to the one of the chromate resistance transposon Tn5719. Noticeable are hits to beta-lactamase protein families which suggests that plasmids from WWTP bacteria encode different enzymes possessing beta-lactam-hydrolysing activity. Some of the sequence reads correspond to antibiotic resistance genes that were only recently identified in clinical isolates of human pathogens. EGT analysis thus proofed to be a very valuable method to explore genetic diversity and composition of the present plasmid metagenome dataset.

The prpZ gene cluster encoding eukaryotic-type Ser/Thr protein kinases and phosphatases is repressed by oxidative stress and involved in Salmonella enterica serovar Typhi survival in human macrophages

The prpZ gene cluster consists of three ORFs coding for proteins with homology to eukaryotic-type Ser/Thr protein phosphatases 2C (prpZ) and Ser/Thr protein kinases (prkY and prkX). This cluster is present in the sequenced genomes of Salmonella enterica serovar Typhi (S. Typhi) strains Ty2 and CT18. This study investigated the genetic organization of this gene cluster, its regulation and its putative involvement in virulence. The three genes are transcribed as a polycistronic mRNA as demonstrated by reverse transcriptase (RT)-PCR. Analysis of a prpZ::lacZ transcriptional fusion showed that the prpZ locus is expressed throughout the growth phase. LacZ activity and real-time RT-PCR showed that transcription of the mRNA is negatively regulated upon exposure of cells to HOCl and, to a lesser extent, hydrogen peroxide. A deletion mutant of the prpZ gene cluster showed a significantly lower level of survival than the parental strain Ty2 in human macrophages at 48 h postinfection. Together these data suggest that prpZ, prkY and prkX are virulence genes that may be part of a signaling pathway controlling long-term survival of S. Typhi in host cells.

Prophage sequences defining hot spots of genome variation in Salmonella enterica serovar Typhimurium can be used to discriminate between field isolates

Sixty-one Salmonella enterica serovar Typhimurium isolates of animal and human origin, matched by phage type, antimicrobial resistance pattern, and place of isolation, were analyzed by microbiological and molecular techniques, including pulsed-field gel electrophoresis (PFGE) and plasmid profiling. PFGE identified 10 profiles that clustered by phage type and antibiotic resistance pattern with human and animal isolates distributed among different PFGE profiles. Genomic DNA was purified from 23 representative strains and hybridized to the composite Salmonella DNA microarray, and specific genomic regions that exhibited significant variation between isolates were identified. Bioinformatic analysis showed that variable regions of DNA were associated with prophage-like elements. Subsequently, simple multiplex PCR assays were designed on the basis of these variable regions that could be used to discriminate between S. enterica serovar Typhimurium isolates from the same geographical region. These multiplex PCR assays, based on prophage-like elements and Salmonella genomic island 1, provide a simple method for identifying new variants of S. enterica serovar Typhimurium in the field.

Evolution of Chlamydia trachomatis diversity occurs by widespread interstrain recombination involving hotspots

Chlamydia trachomatis is an obligate intracellular bacterium of major public health significance, infecting over one-tenth of the world's population and causing blindness and infertility in millions. Mounting evidence supports recombination as a key source of genetic diversity among free-living bacteria. Previous research shows that intracellular bacteria such as Chlamydiaceae may also undergo recombination but whether this plays a significant evolutionary role has not been determined. Here, we examine multiple loci dispersed throughout the chromosome to determine the extent and significance of recombination among 19 laboratory reference strains and 10 present-day ocular and urogenital clinical isolates using phylogenetic reconstructions, compatibility matrices, and statistically based recombination programs. Recombination is widespread; all clinical isolates are recombinant at multiple loci with no two belonging to the same clonal lineage. Several reference strains show nonconcordant phylogenies across loci; one strain is unambiguously identified as recombinantly derived from other reference strain lineages. Frequent recombination contrasts with a low level of point substitution; novel substitutions relative to reference strains occur less than one per kilobase. Hotspots for recombination are identified downstream from ompA, which encodes the major outer membrane protein. This widespread recombination, unexpected for an intracellular bacterium, explains why strain-typing using one or two genes, such as ompA, does not correlate with clinical phenotypes. Our results do not point to specific events that are responsible for different pathogenicities but, instead, suggest a new approach to dissect the genetic basis for clinical strain pathology with implications for evolution, host cell adaptation, and emergence of new chlamydial diseases.

ISCR elements: novel gene-capturing systems of the 21st century?

"Common regions" (CRs), such as Orf513, are being increasingly linked to mega-antibiotic-resistant regions. While their overall nucleotide sequences show little identity to other mobile elements, amino acid alignments indicate that they possess the key motifs of IS91-like elements, which have been linked to the mobility ent plasmids in pathogenic Escherichia coli. Further inspection reveals that they possess an IS91-like origin of replication and termination sites (terIS), and therefore CRs probably transpose via a rolling-circle replication mechanism. Accordingly, in this review we have renamed CRs as ISCRs to give a more accurate reflection of their functional properties. The genetic context surrounding ISCRs indicates that they can procure 5' sequences via misreading of the cognate terIS, i.e., "unchecked transposition." Clinically, the most worrying aspect of ISCRs is that they are increasingly being linked with more potent examples of resistance, i.e., metallo-beta-lactamases in Pseudomonas aeruginosa and co-trimoxazole resistance in Stenotrophomonas maltophilia. Furthermore, if ISCR elements do move via "unchecked RC transposition," as has been speculated for ISCR1, then this mechanism provides antibiotic resistance genes with a highly mobile genetic vehicle that could greatly exceed the effects of previously reported mobile genetic mechanisms. It has been hypothesized that bacteria will surprise us by extending their "genetic construction kit" to procure and evince additional DNA and, therefore, antibiotic resistance genes. It appears that ISCR elements have now firmly established themselves within that regimen.

Expression of phage P4 integrase is regulated negatively by both Int and Vis

Phage P4 int gene encodes the integrase responsible for phage integration into and excision from the Escherichia coli chromosome. Here, the data showing that P4 int expression is regulated in a complex manner at different levels are presented. First of all, the Pint promoter is regulated negatively by both Int and Vis, the P4 excisionase. The N-terminal portion of Int appears to be sufficient for such a negative autoregulation, suggesting that the Int N terminus is implicated in DNA binding. Second, full-length transcripts covering the entire int gene could be detected only upon P4 infection, whereas in P4 lysogens only short 5′-end covering transcripts were detectable. On the other hand, transcripts covering the 5′-end of int were also very abundant upon infection. It thus appears that premature transcription termination and/or mRNA degradation play a role in Int-negative regulation both on the basal prophage transcription and upon infection. Finally, comparison between Pint–lacZ transcriptional and translational fusions suggests that Vis regulates Int expression post-transcriptionally. The findings that Vis is also an RNA-binding protein and that Int may be translated from two different start codons have implications on possible regulation models of Int expression.

A hypervariable 130-kilobase genomic region of Magnetospirillum gryphiswaldense comprises a magnetosome island which undergoes frequent rearrangements during stationary growth

Genes involved in magnetite biomineralization are clustered in the genome of the magnetotactic bacterium Magnetospirillum gryphiswaldense. We analyzed a 482-kb genomic fragment, in which we identified an approximately 130-kb region representing a putative genomic "magnetosome island" (MAI). In addition to all known magnetosome genes, the MAI contains genes putatively involved in magnetosome biomineralization and numerous genes with unknown functions, as well as pseudogenes, and it is particularly rich in insertion elements. Substantial sequence polymorphism of clones from different subcultures indicated that this region undergoes frequent rearrangements during serial subcultivation in the laboratory. Spontaneous mutants affected in magnetosome formation arise at a frequency of up to 10(-2) after prolonged storage of cells at 4 degrees C or exposure to oxidative stress. All nonmagnetic mutants exhibited extended and multiple deletions in the MAI and had lost either parts of or the entire mms and mam gene clusters encoding magnetosome proteins. The mutations were polymorphic with respect to the sites and extents of deletions, but all mutations were found to be associated with the loss of various copies of insertion elements, as revealed by Southern hybridization and PCR analysis. Insertions and deletions in the MAI were also found in different magnetosome-producing clones, indicating that parts of this region are not essential for the magnetic phenotype. Our data suggest that the genomic MAI undergoes frequent transposition events, which lead to subsequent deletion by homologous recombination under physiological stress conditions. This can be interpreted in terms of adaptation to physiological stress and might contribute to the genetic plasticity and mobilization of the magnetosome island.