- Research article
- Open Access
Comparative genomics of Salmonella entericaserovars Derby and Mbandaka, two prevalent serovars associated with different livestock species in the UK
BMC Genomics volume 14, Article number: 365 (2013)
Despite the frequent isolation of Salmonella enterica sub. enterica serovars Derby and Mbandaka from livestock in the UK and USA little is known about the biological processes maintaining their prevalence. Statistics for Salmonella isolations from livestock production in the UK show that S. Derby is most commonly associated with pigs and turkeys and S. Mbandaka with cattle and chickens. Here we compare the first sequenced genomes of S. Derby and S. Mbandaka as a basis for further analysis of the potential host adaptations that contribute to their distinct host species distributions.
Comparative functional genomics using the RAST annotation system showed that predominantly mechanisms that relate to metabolite utilisation, in vivo and ex vivo persistence and pathogenesis distinguish S. Derby from S. Mbandaka. Alignment of the genome nucleotide sequences of S. Derby D1 and D2 and S. Mbandaka M1 and M2 with Salmonella pathogenicity islands (SPI) identified unique complements of genes associated with host adaptation. We also describe a new genomic island with a putative role in pathogenesis, SPI-23. SPI-23 is present in several S. enterica serovars, including S. Agona, S. Dublin and S. Gallinarum, it is absent in its entirety from S. Mbandaka.
We discovered a new 37 Kb genomic island, SPI-23, in the chromosome sequence of S. Derby, encoding 42 ORFS, ten of which are putative TTSS effector proteins. We infer from full-genome synonymous SNP analysis that these two serovars diverged, between 182kya and 625kya coinciding with the divergence of domestic pigs. The differences between the genomes of these serovars suggest they have been exposed to different stresses including, phage, transposons and prolonged externalisation. The two serovars possess distinct complements of metabolic genes; many of which cluster into pathways for catabolism of carbon sources.
Salmonella enterica subspecies enterica is an important zoonotic pathogen of warm-blooded vertebrates, with both a broad host species range and geographical distribution. The subspecies is divided into over 1530 serovars based on the different epitopes of two surface antigens, the O lipopolysaccharide, and H flagellum of which there are commonly two phases .
Some serovars display association with a particular set of hosts that may be stable over many decades and large geographical distances suggesting a level of adaptation or restriction . With regard to serovars S. Derby and S. Mbandaka, both serovars are isolated with similar frequency in the UK and USA. Annually compiled statistics from several sources [3, 4] (HPA personal communication) showed that, whilst both serovars can readily cause disease in people, incidences in livestock show differing host associations. In the UK, for example, approximately 50% and 40% of incidences of S. Derby are in turkeys and pigs, respectively, and approximately 20% and 65% of incidences of S. Mbandaka are from cattle and chickens, respectively . In the USA approximately 80% of isolation of S. Derby are from pigs, while only 3% of isolations were from turkeys, 27% and 25% of S. Mbandaka isolations are from cattle and chickens. Unlike in the UK in the USA S. Mbandaka is isolated from pigs comprising 14% of the total . These host distributions have been maintained for over a decade and on two continents which gives rise to at least two hypotheses. First, is it possible that the differences in host association may relate to production systems and that these serotypes posses similar functional capabilities. Second, is it possible that the differences in host association reflect functional differences between serovars or genovars therein, whereby there exist bacterially encoded mechanisms that maintain these patterns. As a starting point to tackle these opposing hypotheses, we present the first full chromosome sequence of two UK isolates of both S. Derby and S. Mbandaka. We use functional genomics to describe genome features and to identify genes that are unique with a view to gaining insights into potential genetic components that contribute to the species distributions described above.
Results and discussion
The chromosomes of two strains of S. Derby and S. Mbandaka were sequenced and compared with the goal of identifying potential mechanistic differences between the two serovars that could explain their skewed isolation frequencies from subsets of livestock species in the UK. Strains were obtained from background monitoring performed by the Animal Health and Veterinary Laboratories Agency (AHVLA) in the UK between 2000 and 2010. In total 28 strains were selected spanning the decade and from differing geographic points of isolation across the UK (locations not shown due to sensitivity of data). The hosts of isolation of the selected strains were chosen to reflect the two most common hosts of each serovar, for S. Derby these were pigs and turkeys and for S. Mbandaka cows and chickens. Two isolates of each serovar isolated from separate geographical locations, with the same host species, and identical MLST sequence types (S. Derby strains ST40 and S. Mbandaka strains ST900) were chosen for full genome sequencing. We recognised that in the absence of information regarding the pan-genome of the population, that by comparing just two isolates of each serovar, we could potentially infer, incorrectly, that differences in gene complement between isolates of the same serovar isolated from different hosts were adaptations to these different hosts. The selection was therefore made with the aim of better understanding the genomic differences between strains which would typically be considered clonal. S. Derby strains D1 and D2 were both isolated in 2008 from porcine hosts. S. Mbandaka M1 and M2 were isolated from cattle in 2008 and 2009 respectively. No research has previously been performed on these strains.
General genome features of S. Derby D1 and D2 and S.Mbandaka M1 and M2
S. Derby strains D1 and D2 possessed chromosomes of 4.86 Mb nucleotides in length with a GC skew of 51.16% and 51.46% respectively. The RAST annotation system predicted that the chromosome sequence of S. Derby D1 encodes 4720 genes and the sequence of D2 4717 genes. The chromosome of S. Mbandaka strains M1 and M2 were both 4.72 MB nucleotides in length with a GC skew of 51.91% and 52.01% respectively. These were predicted to encode 4616 and 4619 genes respectively. Interestingly all four chromosomes contain different numbers of RNA coding sequences, D1 contains 69, D2 contains 73, M1 contains 74 and M2 contains 75 (Figure 1). RNA sequences are frequently sights for integration of horizontally acquired DNA sequences, in some cases leading to duplication of the RNA [5–7]. The difference in the number of RNAs in each genome could reflect a difference in evolutionary potential of each chromosome.
S. Mbandaka contains a large sequence inversion
S. Mbandaka contains a 860Kb sequence inversion between a mobile element protein and tRNA-ser-GGA (located between base 1086415 and 1947250 of M1, and 1132370 and 1992477 of M2) which was also found in S. Choleraesuis SC-B67, and was absent from S. Derby (Figure 1) and other sequenced S. enterica serovars including S. Agona SL483, S. Dublin CT02021853, S. Enteritidis P125109, S. Gallinarum 28791 and S. Typhimurium LT2 and SL1344. This region codes for 909 genes identified by the RAST gene caller. Large sequence inversions have a significant impact on the transcript composition of the cell during replication, as those genes closer to the origin of replication are present in duplicate for a longer period of time than those genes closer to the terminus of replication [8, 9]. The effects of increased gene dosage during replication are most noticeable when bacteria are growing at an optimal rate . In Escherichia coli DNA replication from the origin of replication to terminus of replication takes 22 minutes during a 40 minute cell cycle when grown in LB broth at 37°C . If we apply this duration to the inversion found in S. Mbandaka M1 and M2, where almost a quarter of the chromosome is in a different orientation to S. Derby D1 and D2, then there is an 8.6 minute difference between gene duplication events of the genes adjacent to the sites of inversion. These genes are therefore in duplicate and the other genes in singlet for 21% of the cell cycle. In S. Derby the ten genes closest to the mobile genetic element signifying the start of the inverted sequences do not pertain to a common mechanism. Though interestingly, amongst these ten genes is a permease of the drug/metabolite transporter (DMT) superfamily which in S. Mbandaka occupies the very furthest gene in the inversion. The ten genes at the terminus of the inversion in the chromosome of S. Derby D1 and D2 comprise of two operons, the Csg-curli operon (four genes) and Ycd-swarming operon (five genes). The most interesting aspect of these two operons is that they are associated with two diametrically opposed mechanisms; the curli operon is associated with biofilm development in a sessile population and the swarming operon is associated with directed movement of the bacterial population, both using quorum sensing [12, 13]. Both are population scale emergent phenotypes of gene regulation at a single cell level. csgD found in the curli operon regulates both curli expression and cellulose secretion, the main components of biofilms in Salmonella enterica. The ten genes at the centre of inversion and therefore in similar dosage throughout replication in S. Derby and S. Mbandaka are genes for formate dehydrogenase alpha, beta and gamma subunits which form a single transmembrane enzyme . Also contained within this region is a permease of the drug/metabolite transporter (DMT) superfamily. The ten genes at either end and in the centre of the sequence inversion can be found listed in the supplementary materials [Additional file 1].
Summary of functional annotation
The chromosome of S. Derby is 140 kb longer than that of S. Mbandaka, coding for 100 additional genes. RAST annotation was performed on 9/10/12 and achieved 67% coverage with FIGfam subsystems of the S. Derby D1 and D2 chromosomes and 68% of the S. Mbandaka M1 and M2 chromosomes . FIGfam clusters genes based on protein sequence similarity. From this the function of a novel gene may be inferred. These genes are then clustered into hierarchical subsystems that display increasing functional breadth . As the database has developed the subsystem coverage of the four genomes presented here has markedly increased by 8-9% over the same annotation performed on the 9/10/11. The most recent annotation is available through the following RAST IDs: D1 [RAST: 28144.16], D2 [RAST: 28144.17], M1 [RAST: 192954.16] and M2 [RAST: 192954.17] so that the chromosome can be re-annotated as the RAST databases are updated. The number of hypothetical genes between the annotations remained constant for each chromosome. In all cases almost a quarter of the genome annotation was found to be of hypothetical gene status. S. Derby contains 96 hypothetical/putative proteins which share less than 90% amino acid sequence homology with open reading frames in S. Mbandaka. S. Mbandaka contains 155 unique hypothetical/putative proteins.
Intra-serovar differences in the complement of functionally unique genes
The majority of the diversity between strains of the same serovars was in the complement of phage associated genes. All intra-serovar differences in gene complement can be found listed in the supplementary materials [Additional file 1]. Between the S. Derby isolates, D1 contains a single unique gene for an aconitate hydratase 2 (EC 18.104.22.168) associated with glyoxylate bypass. D2 contains 11 unique genes, of which five are associated with phage, the remaining are associated with metabolism. A single gene which is associated with the ribosome at stationary growth phase is absent from D1. There is less diversity between isolates of S. Mbandaka. M1 contains a single additional gene which encodes a phage tail fibre protein. M2 contains six additional genes, two for cytochrome-c biosynthesis, two phage genes and a gene encoding a 2,5-diketo-D-gluconic acid reductase B (EC 22.214.171.1244).
Inter-serovar differences of functionally unique genes
Genes pertaining to metabolite utilisation, prophage, CRISPR spacers and Salmonella pathogenicity islands will be dealt with separately. The following summarises the genes that do not fit into these categories. All inter-serovar differences in gene complement can be found listed in the supplementary materials [Additional file 1].
S. Derby contains 16 genes that are functionally unique. Of these 16 genes 12 are distributed between two operons. One is the mer operon conveying mercury resistance. This consists of five genes including merC and merT transport proteins, which actively take up toxic mercuric cations (Hg+2) for subsequent reduction to non-toxic metallic mercury (Hg0) . Mecuric cations enter the food chain from several sources, including fish, poultry and meat. Animals feed is frequently supplemented in the UK with fish meal which is high in mercury. Fish meal contains the second highest concentration of mercury per kg in animal feed/ pet food in Europe, with fish oil containing the highest concentration . The second operon, is a CRISPR operon made up of seven genes, one of putative status and cse1-cse4, cas1, cas2 and cas5e. Of the remaining four genes one is the gene for an UDP-galactopyranose mutase (EC 126.96.36.199) which is associated with the biosynthesis of alpha-D-galactofuranose, a component of the O-antigen in Salmonella enterica groups B, C2, D and E .
S. Mbandaka contains a cluster of four type VII secretion system Yad fimbrial chaperone proteins. Two genes with the same function from the HtrE fimbrial cluster sit approximately 100 KB away. A further 2 MB away at 4.5U sits a cluster of three beta-fimbriae genes also associated with type VII secretion. Three sialic acid metabolism genes associated with capsule production, nanC, nanM and a hypothetical gene, are clustered around 1U. Two cell death toxin-antitoxin genes, phd and doc, are unique to M1 and M2, and may be involved in plasmid addiction systems. Two genes associated with reduction in mutation rate due to exposure to bile salts are absent from S. Derby. These genes, umuC and umuD are part of the SOS DNA repair response and form DNA polymerase V. It has been shown in E. coli that in the absence of umuC genomic lesions are not repaired correctly by DNA polymerase III and can leave frame shift mutations which lead to pseudogene formation. DNA polymerase V has a higher rate of single nucleotide mutations than DNA polymerase III [21, 22]. This could lead to a higher rate of pseudogene formation in S. Mbandaka strains and SNP formation in S. Derby strains. However, this would need to be confirmed through further analyses.
There are only seventeen genes that are unique in function to either S. Derby or S. Mbandaka that are not clustered. Of these seventeen genes S. Mbandaka contains seven unique genes related to biogenesis of cytochrome-c, specifically the maturation of the molecule, and are spread across the chromosome. The genes ccmB, ccmC and ccmD convey the heme-b group to the product of CcmE, a monotopic membrane protein . The products of ccmF, ccmG and ccmH complex with CcmE to convey the heme-b group to the apocytochrome-c precursor of cytochrome-C [24, 25]. Though these genes are ubiquitous amongst Gram negative bacteria, strains of E. coli have been discovered that lack the ccm operon and yet are able to synthesis cytochrome-c containing heme-b .
Differences in metabolic gene complement between S. Derby and S. Mbandaka
Fourteen genes were identified by RAST subsystem annotations as being involved in primary or secondary metabolism which were found to differ between S. Derby and S. Mbandaka. Six of these genes belong to S. Mbandaka are associated with D-galactonate catabolism, this includes uptake, regulation and processing into central carbon metabolism. S. Derby contains six genes for the uptake and catabolism of six different carbon sources, this comprises an asparagine synthetase (EC 188.8.131.52), a hydroxyaromatic non-oxidative decarboxylase protein D (EC 4.1.1.-), a protein fumarylacetoacetate of the hydrolase family, phosphatase NagD predicted to act inN-acetylglucosamine utilization subsystem, an aconitate hydratase 2 (EC 184.108.40.206), a galactose-specific IIA component (EC 220.127.116.11) and the large subunit of a glycerol dehydratase reactivation factor.
The biological significance of the differences in metabolic genes was elaborated through construction of metabolic models from the genome sequences using SEEDmodel . These differences were then elaborated in context of the surrounding reactions. Metabolic reconstructions curated with phenotypic data are underway to better understand the effect of secondary metabolism on the optimal growth rate of S. Derby D1 and S. Mbandaka M1.
Alanine, aspartate and glutamate metabolism map 00250 created 1/6/12
S. Derby lacks a single gene, an aspartate—ammonia ligase (EC18.104.22.168) for the conversion of L-aspartate to L-asparagine. The same reaction is achievable through two additional reactions utilising an asparaginase/glutaminase (EC22.214.171.124) and an L-asparaginase (EC126.96.36.199) which are also present in S. Mbandaka.
Galactose metabolism map 00052 created 31/5/12
The three genes encoding products needed to feed D-galactonate into glycolysis (EC 188.8.131.52, EC 184.108.40.206 and EC 4.1.21) by conversion to D-glyceraldehyde-3P are present on the chromosome of S. Mbandaka and absent from that of S. Derby. There are no alternative routes from D-galactonate to glycolysis.
Nitrogen metabolism map 00910 created 21/8/12
A gene coding for the enzyme L-glutamine amido-ligase that converts L-glutamine to L-glutamate using one molecule of H2O in the process (EC 220.127.116.11) is missing from the chromosome of S. Derby D1. All strains contain a gene that catalyses the same reaction but with the requirement of a molecule of NADP+ as opposed to one of H2O (EC 18.104.22.168).
Starch and sucrose metabolism map 00500 created 9/7/12
A single reaction is missing from S. Mbandaka in this map for the conversion of alpha-D-Glucose-1-P to CDP-glucose (EC 22.214.171.124); there is no route to this compound other than this on the map. The CDP-glucose then leads into amino sugar and nucleotide sugar metabolism map 00520 created 19/1/10. In this map there is an additional reaction from CDP-glucose leading to CDP-4-keto-6-deoxy-D-Glucose missing in S. Mbandaka. This reaction is catalysed by the enzyme RfbG, a CDP-glucose 4,6-dehydratase (EC 126.96.36.199) which is found in Salmonella enterica groups A, B, C2, C3, D1 and D2 and required for binding of the O antigen to the core oligosaccharide [27, 28]. S. Mbandaka is a member of S. enterica group C1.
Streptomycin biosynthesis pathway map 00521 created 27/12/10
Two steps from D-glucose-1-P are present in both serovars (EC 188.8.131.52 and EC 184.108.40.206), following on from the terminal product of this reaction, two additional steps that lead to dTDP-L-rhamnose are missing in S. Mbandaka (EC 220.127.116.11 and EC 18.104.22.168). dDTP-L-Rhamnose feeds directly into novobiocin biosynthesis, diverted out of the streptomycin biosynthesis pathway. S. Mbandaka is left with a product which feeds into polyketide sugar unit biosynthesis (Pathway 00523, created 14/3/12, polyketide sugar unit biosynthesis).
The chromosome of Salmonella enterica comprises largely of a core sequence punctuated with horizontally acquired sequences . The complement of genomic islands within the chromosome of Salmonella enterica can vary amongst isolates of the same serovar [30, 31]. It has been postulated that the acquisition of horizontally acquired genes into a Salmonella pathogenicity island (SPI) led to the divergence of Salmonella from Escherichia coli[32, 33]. Salmonella pathogenicity island 1 (SPI-1) is found in all serovars of S. enterica (with the exception of S. Seftenberg and S. Litchfield) and is highly conserved [32, 34, 35]. There are currently 22 published Salmonella pathogenicity islands identified from the genomes of Salmonella enterica and Salmonella bongori. The gene content of some of these islands is highly plastic, as exemplified by the different gene complement of SPI-3 found in S. Dublin CT02021853 and S. Typhimurium LT2 . The Salmonella pathogenicity islands are well characterised in terms of genetic composition and putative function but less so, with notable exceptions, for their role in pathogenicity [38, 39]. Hence differences in SPI complement and gene content of D1, D2, M1 and M2 chromosomes may hint at mechanisms that maintain their respective host species range.
Complete or absent Salmonellapathogenicity islands
SPIs 2 and 4 found in the genome of S. Choloreaesuis SC-B67 and SPI-18 from S. Typhi CT18 are complete in the genomes of S. Derby D1 and D2, and S. Mbandaka M1 and M2. SPI-7, 8, 10, 15, 16, 17, 19, 20, 21 and 22 were absent from both S. Derby D1 and D2, and S. Mbandaka M1 and M2 genomes.
Variation in SPI-1of S. Derby and S.Mbandaka
SPI-1 in S. Mbandaka M1 and M2 shares 100% nucleotide sequence identity with S. Typhumirum LT2 with the addition of two ORFs coding for hypothetical proteins found in the SPI-1 of S. Choleraesuis SC-B67, SC2837 and SC2838 which are absent in S. Derby D1 and D2. S. Derby D1 and D2 lack three genes from SPI-1 of S. Typhimurium LT2, STM2901, STM2902 and STM2903 (Table 1). SIEVE an online server for the prediction of TTSS effector proteins, found that the S. Mbandaka M1 and M2 contained an ORF with 98% amino acid sequence homology with SC2837 from S. Choleraesuis SC-B67, is a likely candidate for an effector protein with a p-value of 0.003. With reference to well-characterised effector proteins, all four isolates contain intact versions of sopB and sopE. The two putative cytoplasmic proteins found in SPI-1 of S. Typhimurium LT2, STM2901 and STM2902 and here in S. Mbandaka M1 and M2 and not D1 and D2 are unlikely candidates for effector proteins with p-values of 0.142.
Variation in SPI-3 between other serovars and S. Derby and S.Mbandaka
SPI-3 is highly variable, between S. Typhimurium 14028 and S. Choleraesuis SC-B67 the only region of homology is the insertion sequence tRNA-selC. SPI-3 from S. Derby D1 and D2 is an amalgamation of 19 SPI-3 genes from S. Typhimurium 14028, S. Dublin, S. Choleroaeasuis SC-B67 and S. Typhi CT18. S. Mbandaka M1 and M2 also contain a unique SPI-3 gene complement, containing 12 genes found in S. Typhimurium 14028, S. Choleraesuis SC-B67 and S. Typhi CT18. Unlike S. Derby D1 and D2, S. Mbandaka M1 and M2 have no SPI-3 genes in common with S. Dublin. STY4039 previously unique to S. Typhi CT18 is present in S. Mbandaka M1 and M2 and absent from S. Derby D1 and D2 (Table 1). The main region of variation between S. Derby D1 and D2 and S. Mbandaka M1 and M2 SPI-3 is at the start of the island where the complete S. Dublin SPI3 is present, this was shown previously for S. Derby 9813031, 0010160 and 0010158 [38, 39]. This region contains seven genes relating to the adhesion structures, pili and fimbriae. S. Mbandaka M1 and M2 contains rhuM found in the SPI-3 of S. Typhimurium 14028; this sequence is absent from the SPI-3 of S. Derby D1 and D2. S. Derby D1 and D2 and S. Mbandaka M1 and M2 share 10 SPI-3 genes in common; this complement of genes is unique to these two serovars. Both serovars contain five virulence genes present in the SPI-3 of S. Typhimurium 14028 and S. Choleraesuis SC-B67. Both serovars lack the virulence gene mgtC which is present in S. Typhimurium 14028, S. Choloraesuis SC-B67 and S. Typhi CT18. In S. Typhimurium LT2 and 14028 mgtC was shown to be essential for intra-macrophage survival .
Variation in SPI-5 between other serovars and S. Derby and S.Mbandaka
It has previously been shown that SPI-5 from S. Derby 9813031, 0010160, 0010158 and S. Ohio 9815932, 9714920, 9714922 contain an additional unnamed ORF, this ORF was present in both S. Derby D1 and D2 and S. Mbandaka M1 and M2 .
Variation in SPI-6 between S. Derby and S.Mbandaka
SPI-6 is found in S. Typhi CT18, is 57 Kb in length and contains 59 genes. Between S. Derby D1 and D2 and S. Mbandaka M1 and M2, 24 genes that were not found in other islands on PAI-DB were identified by Glimmer3 (Table 1) . The annotations here were taken from NCBI BLASTn results, many of which were hypothetical or putative in description. SPI-6 also shows the largest variation between S. Derby D1 and D2 and S. Mbandaka M1 and M2 outside of prophage and SPI-23 nucleotide sequences. SPI-6 in D2 had 8 unique genes at the C terminus of the positive strand that were not found in the other isolates. This contains an AP like endonuclease gene related to egg albumin resistance, yafD. S. Derby has been isolated from inside of eggs while S. Mbandaka has been shown to grow slower than other S. enterica serovars in albumin . The remainder of the island showed no variation amongst isolates of the same serovar. Seven S. Typhi CT18 genes were absent from both serovars, these were STY0300-STY303, STY0342, STY0350 and STY03351. S. Derby D1 and D2 SPI-6 contained 8 genes from S. Typhi CT18 that were absent from S. Mbandaka M1 and M2. S. Mbandaka M1 and M2 SPI-6 contained 20 genes from S. Typhi CT18 that were absent from S. Derby D1 and D2.
The variation in the gene complement of SPI-6 in S. Derby and S. Mbandaka is of particular interest with regards to host adaptation. S. Derby possess the gene sirA which corresponds with ORF STY0300 in S. Typhi CT18, that codes for a transcription factor linked with regulation of the TTSS encoding SPI-1 when inside a mammalian host [8, 9]. Interestingly mutants for sirA in S. Typhimurium LT2 were attenuated in a bovine gastro-enteritis model, but were still proficient at causing typhoid fever in a mouse model [47, 48, 50, 64]. The SPI-6 of S. Mbandaka also contains a gene sciN which is an outer membrane lipoporotein essential for biofilm formation in E. coli which is absent from S. Derby .
Variation in SPI-9 from S.Typhi CT18
The alignment between SPI-9 of S. Derby D1 and D2 and S. Mbandaka M1 and M2 showed 100% sequence homology. SPI-9 from S. Typhi CT18 contains four genes as do the islands in S. Derby D1 and D2 and S. Mbandaka M1 and M2, though there is a difference in ORF length. STY2875 is at the start of the island and is 10.8 kb in length, in both S. Derby D1 and D2 and S. Mbandaka M1 and M2 an additional region of 595 bp is found between bases 3056 and 3057. The other three ORFs are truncated at the beginning of each sequence by 162 bp.
Variation in SPI-11 from S. Derby and S.Mbandaka
The same eight genes from SPI-11 of S. Choleroeaesuis SC-B67 are absent in S. Derby D1 and D2 and S. Mbandaka M1 and M2. One of these genes is the effector protein sopB, which has been implicated in fluid secretion in calf ileal loops and is essential for enteropathogenicity of S. Dublin [65, 66] although, as previously mentioned, a homolog to this gene was found elsewhere on the chromosome. SPI-11 also encodes the gene pagC, an envelope protein which increases survival within mouse macrophage .
Variation in SPI-12 between S. Derby and S.Mbandaka
SPI-12 is an 11 Kb island first identified in S. Choleraesuis SC-B67. The island is inserted at a tRNA-Pro. The insertion sequence was present in both S. Derby D1 and D2 and S. Mbandaka M1 and M2, though no genes were adjacent to this site. Alignment of the whole SPI-12 island from S. Choleroeaesuis SC-B67 with D1, D2, M1 and M2 identified homologs for all the genes in each sequence, not in a single unit, but spread across the chromosome.
Variation in CS54 between S. Derby and S.Mbandaka
CS54 identified from S. Typhimurium 14028 is associated with virulence and shows variation between isolates of the same serovar. All isolates lack the virulence genes shdA and ratB, and the untested gene ratC. M2 lacks the whole island with only the insertion sequences present (Table 1). This region in D1, D2 and M1 contains three genes ratA, sivI and sivH previously identified in CS54 of S. Typhimurium 14028. CS54 was previously described in S. Derby strain De1, in this instance ratB was also found, a gene essential for the colonisation of the cecum in BALB/c mice by S. Typhimurium IR715 .
A novel Salmonellapathogenicity island designated SPI-23
A new genomic island with a putative role in pathogenesis, SPI-23, was discovered in this study on the chromosome of D1 and D2 between bases 2027348–2065972 and 2052685–2089962 respectively flanked by tRNA-asn (GTT) and a hypothetical protein, docB (Figure 2). SPI-23 is composed of 42 ORFS with an overall GC composition of 38% differing largely from the 51% of the S. Derby genome. SPI-23 is completely missing from S. Mbandaka. Of the 42 ORFs 28 were of hypothetical status, of which, 17 contained no homology with an entry in the NCBI nucleotide database (accessed on 1/10/12). The island contains two genes, potR and talN, both implicated in type IV secretion and the production of pili. There is a single gene, zomB, predicted here to encode a lipoprotein. We also find five DNA binding proteins, furB, lamE, halF, mstR and numT and two putative membrane protein bigM and putM. SPI-23 contains a single NUDIX hydrolase, a very ubiquitous protein family involved in a multitude of regulator processes .
SIEVE effector protein predictor identified ten ORFs (sanA, janE, chlE, yuaM, genE, shaU, dumE, sadZ, tinY and docB) in SPI-23 of S. Derby D1 and D2 with a p-value of 0.05 or lower corresponding to a Z-Score of 1.5 or higher (Table 2). docB, encoding a putative effector protein was identified by RAST as a putative endoprotease and was found here to be conserved in S. Derby D1 and D2, S. Mbandaka M1 and M2, S. Agona SL483, S. Dublin ct02021853, S. Gallinarum SGG1, S. Enteritidis P125109, S. Newport SL254, and S. Typhimurium LT2. The functional prediction of docB fits with the function of other type III secretion effector proteins which have a cysteine protease activity . The high number of potential type III secretion system effector proteins makes SPI-23 a strong candidate for classification as a pathogenicity island. The acquisition of this sequence could be responsible for the modulation of the host’s cell, cytoskeleton, immune response and intracellular signalling [71, 72]. Though it is not possible to determine here if SPI-23 plays a role in defining the host range of S. Derby, the high number of potential effector proteins has identified it as a very interesting region for future experimental study of host adaptation.
Comparison of SPI-23 from S. Agona SL483, S. Dublin ct02021853 and S. Gallinarum SGG1 with SPI-23 found in S.Derby D1 and D2
There were no genes between gooN and docB in S. Enteritidis P125109 even though NCBI BLASTn showed 100% sequence homology with 17 genes from SPI-23 of S. Derby D1 and D2. A four way comparison between SPI-23 excised from the genomes of S. Agona SL483, S. Dublin CT02021853, S. Gallinarum SGG1 and S. Derby D1 was performed (Figure 3). The differences in SPI-23 between S. Agona SL483 and S. Derby D1 and D2 are dispersed across the island in four sections (Table 2). S. Agona SL483 contains seventeen unique genes and lacks twenty two genes when compared to the SPI-23 of S. Derby D1 and D2. All of the genes unique to S. Agona SL483 with the exception of three are of hypothetical status, relA a GDP/GTP pyrophosphokinase, a putative component of the TonB system and a P4 type intergrase. SPI-23 in S. Agona SL483 contains only four genes that are likely candidates for type III secretion system effector proteins (sanA, kayT, mstR and docB). All four genes are identical in nucleotide sequence to that of S. Derby D1. Serovars S. Dublin CT02021853 and S. Gallinarum SGG1 have identical sequences for SPI-23. There are fifteen genes in the SPI-23 of these two serovars that are not found in either S. Derby D1 and D2 or S. Agona SL483 and fifteen which are found in all five serovars. Only two of the hypothetical genes found in S. Agona SL483 and not S. Derby D1 and D2 are found in the SPI-23 of S. Dublin CT02021853 and S. Gallinarum SGG1. The SPI-23 of S. Dublin CT02021853 and S. Gallinarum SGG1 contains four unique genes that are not of hypothetical status. This comprises two pilV-like proteins, a DNA –binding protein HNS and a threonine operon leader protein. Both S. Dublin CT02021853 and S. Gallinarum SGG1 contain eight putative type III secretion system effector proteins, three of these are unique genes to these two sequences and are absent from S. Derby D1, D2 and S. Agona SL483. Two putative effector proteins are different between the SPI-23 sequences of S. Dublin CT02021853 and S. Gallinarum SGG1, sanA is present in S. Dublin CT02021853 but not identified as a putative effector protein and similarly a hypothetical gene in S. Dublin CT02021853 and S. Gallinarum SGG1. Interestingly SIEVE predicts the gene kayT as a type III secretion system effector protein from the amino acid sequences of S. Agona SL483, S. Dublin CT02021853 and S. Gallinarum SGG1 but not that of S. Derby D1 or D2. Similarly sanA is identified as a candidate effector protein in all sequences with the exception of S. Dublin CT02021853.
Bacteriophages are viruses that infect bacteria, integrating into the bacterial genome in order to replicate; in this form they are known as prophage. As a result of phage insertion the genome gains a substantial amount of foreign sequence, much of which encodes phage structural proteins. However, some phage carry cargo genes which convey a pathological advantage to the recipient . The process of lysogenic conversion prevents the prophage from destroying the host through maturation of progeny. The cargo genes and prophage remnants are therefore retained within the bacterial lineage, undergoing genetic mutation, drift and selection .
PHAST identified distinct complements of intact prophage and remnant prophage regions between S. Derby and S. Mbandaka . All isolates contain four phage regions, sharing only the remnants of a BcepMu phage in common. This remnant is identical in all strains, suggesting that the integration and degradation of this phage predates the split between S. Derby and S. Mbandaka. S. Mbandaka isolates contain the same prophage regions in the same locations along the chromosome. These comprise one intact prophage, resembling phage P2, two questionable prophage, similar to L413c and Epsilon34 and one incomplete prophage BcepMu. S. Derby isolates differed on the location of the prophage within the chromosome and the number of genes in all four phage regions. No ambiguous bases were identified in these regions. The partial prophage resembling SFV contains one additional ORF in D1 than in D2 and occupies the same region that the complete prophage of SFV occupies in D2. Whereas the complete copy of SFV in D1 occupies the position of the complete prophage in D2 and contains one fewer ORF. The BcepMu partial in D1 contains two additional ORFs than that found in D2. In D1 the intact prophage resembling ST64B comprises three additional ORFs than that found in D2, they occupy the same chromosomal region. ST64B is of particular interest as its homolog in S. Typhimurium SL1344 contains a gene with homology to a type III secreted effector protein Sske2, mutants of which have shown to have reduced pathogenicity in a bovine model . S. Derby contains an intact version of IN0, a transposon identified from Pseudomonas aeruginosa.
S. Derby and S.Mbandaka contain unique CRISPR spacer sequences
CRISPR operons convey an adaptive immunity against plasmids and bacteriophage to a broad range of archaeal and bacterial species. This is achieved through integration of unique regions of foreign DNA into the prokaryotic chromosome. Subsequent expressions of these fragments interfere with foreign nucleic acid, through complementation [76, 77]. The spacer sequences within a CRISPR operon reflect the historical interaction between the lineage of a strain and foreign DNA elements. The efficacy of invasion and ecological distribution of bacteriophage, transposons and plasmids have been found to associate with particular hosts and environments [78, 79]. Hence the different genomic complement of prophage and CRISPR operon elements in S. Derby and S. Mbandaka could reflect their particular niche or even define their niche within a specific group of livestock species.
S. Derby D1 and D2 contain four CRISPR operons each, with 34 and 35 spacers respectively. S. Mbandaka M1 contained two CRISPR operons with 25 spacers. M2 contains three CRISPR operons with 27 spacers. With the exception of two spacers, the sequences are completely unique to each serovar. S. Derby isolates contain four CRISPR spacer operons, the smallest contains only one sequence with the largest containing 25 spacers. D2 contains two additional spacer sequences and half of a much larger spacer than D1. S. Mbandaka isolates differ on the number of spacers they each contain; M1 contains two operons while M2 contains three. The majority of spacers are homologous between the isolates, with M2 containing four additional spacers. M2 CRISPR operon 2 and 3 contain all of the spacer sequences in M1 CRISPR operon 1. All spacer sequences can be found in the supplementary materials [Additional file 1].
We have already shown that S. Derby strains contain seven functionally unique CRISPR operon proteins. The lack of functional homolog in S. Mbandaka, leaves it without a functioning CRISPR operon. CRISPRdb shows here that S. Mbandaka strains contain arrays of CRISPR spacer sequences; these may be remnant from when S. Mbandaka had a fully functioning CRISPR operon. S. Mbandaka may now be susceptible to those phage and plasmid for which it once had resistance; this could reflect the loss of positive selection pressure on the operon from the surrounding environment.
Estimating the time since the divergence of S. Derby and S. Mbandaka
Whole genome alignment and SNP calling across CDS nucleotide sequences was used to estimate the years since divergence of D1, D2, M1 and M2. Interestingly the time since the divergence of S. Derby and S. Mbandaka is estimated at between 182,291 and 625,000 years, based on an average of the Ks values for all four pair-wise comparisons of S. Derby and S. Mbandaka isolates ranging between 0.015 and 0.019. The divergence of S. Derby and S. Mbandaka coincides with the estimated time of the divergence of all domesticated pig species, approximately 500,000 years ago [80–83]. The time since the split between D1 and D2 was estimated at between 350 and 1200 years ago based on 31 synonymous SNPs spread across 923506 synonymous positions. The isolates M1 and M2 are estimated to have diverged between 1271 to 4357 years ago based on 118 synonymous SNPs spread across 965114 synonymous positions.
We estimate here that S. Derby D1 and D2 diverged from S. Mbandaka M1 and M2 between 182kya and 625kya, during this period these serovars appear to have adapted towards two distinct ranges of host species. Comparative functional genomics has alluded to several mechanisms that could contribute towards distinct host adaptations of S. Derby D1 and D2 and S. Mbandaka M1 and M2. Most noteworthy of these differences are the diversity in SPI-6 gene complement and the discovery in the chromosome sequence of S. Derby, of a new 37 Kb genomic island, SPI-23 encoding 42 ORFS, ten of which are putative TTSS effector proteins. The absence of functional homologs to several CRISPR operon genes in the chromosome sequences of S. Mbandaka may reduce the fitness of the serovar in environments laden with actively integrative foreign genetic elements. The increased gene dosage of the Csg-biofilm operon and the Ycd-swarming operon in S. Mbandaka could make the implementation of these two behaviours more readily achievable. Both of these behaviours are considered stress responses. S. Mbandaka also possesses an operon pertaining to the uptake and metabolism of D-galactonate into glycolysis which is absent from the chromosome of S. Derby.
The genetic background in which the function of the genes discussed here have been characterised is non-isogenic to the chromosome of S. Derby D1 and D2 and S. Mbandaka M1 and M2. Due to the different context these genes are found in, firm conclusions on the function of these genes in these specific serovars can only be formed through further biological experimentation.
Bacterial strains and culturing
The original isolates were stored at RT on Dorset egg slopes from which bead stocks were made with HIB + 30% glycerol, samples were frozen to −80°C in 2010 and remained frozen throughout the study. Unless stated otherwise, strains were grown for 16 hours aerobically either on LB agar plates or in liquid broth vigorously agitated at 220 rpm.
DNA extraction, genome sequencing and assembly
DNA was extracted from 3 ml overnight cultures as per manufacturer’s instructions (Invitrogen EasyDNA kit). Sequencing was performed by the AHVLA Central Sequencing Unit, Weybridge. A Roche GSFLX titanium 454 pyrosequencer was used to produce rapid and paired-end libraries for whole genome DNA preparations of S. Derby D1 and D2 and S. Mbandaka M1 and M2. Roche protocols were used in all stages of sequencing. Paired-end library inserts were between 4 Kb to 9 Kb, containing 20,000 to 89,000 reads each. The rapid libraries contained between 69,000 and 173,000 reads each. Sequences were assembled de novo using Newbler v2.5. Scaffolds were reordered in ACT v9.0 in reference to a DoubleACT v2 comparison file of each genome with D1; D1 was chosen as the assembly consisted of a single scaffold [84, 85]. The final sequences were then formatted so as to begin at the gene thrL, in line with other published Salmonella enterica genomes.
Automated annotation, metabolic model construction and comparative genomics
Genomes were annotated using the RAST annotation system performed on 9/10/12, backfilling of gaps and automatic error fixing were enabled. Functional comparisons were implemented using the SEED genome viewer v2 [86, 87]. An automated metabolic reconstruction was also produced from the complete genome sequence using the ModelSEED server v1.0 . Differences in S. Derby and S. Mbandaka models were identified through gene overlays on top of KEGG maps . Reciprocal BLASTing was implemented in SEED genome viewer for each ORF that differed between isolates to identify functional homologs. The genomes were also compared through sequence homology. The population of “hypothetical” and “putative” genes were aligned with a cut off of 90% bi-directional amino acid sequence homology.
Mobile genetic elements
SPIs were identified from the genomes of S. Derby and S. Mbandaka through alignment of the insertion sequences with the newly acquired genomes. SPIs for the serovars S. Choleraesuis B67 (SPI-1, 2, 3, 4, 11, 12) and 1240 (SPI-11), S. Derby (SPI-5) isolate not specified, S. Gallinarum SGG1 (SPI-13) and SG8 (SPI-14), S. Dublin isolate not specified (SPI-3), S. Typhi CT18 (SPI-3, 6, 7, 8, 9, 10) S. Typhimurium LT2 (SPI-1, 5, 18) and 14028 (CS54, SPI-3 ) were acquired from PAI-DB website , these were aligned using DoubleACT with the newly isolated islands. From previously published annotated genomes, SPI-22 from S. Bongori  and SPI-15, 16 and 17 from S. Typhi CT18 were excised from the genome . SPI-19, 20 and 21 were excised from the genome of S. enterica subspecies arizonae (IIIa) serotype 62:z4,z23:- . In most cases the SPI from S. Derby and S. Mbandaka could be completely annotated through alignment using DoubleACT with the existing SPI. Where gaps were present the BLAST facility was first used on PAI-DB, when no results were obtained, the RAST annotation and NCBI BLASTn were used to annotate the genes, and extensive literature research was used to assign a putative role in pathogenesis. Prophage were identified and categorised as intact, questionable and partial using PHAST . CRISPR spacers were identified using CRISPRfinder [75, 91]. CRISPRdb BLAST facility was used to see if the spacers found in the newly sequenced genome were found within other bacterial species . Spacer sets of the newly sequenced strains were also cross compared to elicit the historical differences in exposure to phage that has occurred since their divergence. Hypothetical proteins found in SPI-1 and SPI-23 were tested for potential roles as TTSS effector through implementation of SIEVE-SVM based TTSS effector protein predictor . A Z-score above 1.5 was taken to reflect a good indicator of a type III effector protein. SPI-23 was identified and extracted from the publicly available genomes for S. Agona SL483, S. Dublin CT02021853 and S. Gallinarum RKS5078. The sequences were annotated using RAST and SEIVE. Sequences were compared using DoubleACT and ACT.
Estimation of years since S. Derby and S. Mbandaka diverged
For each genome (D1, D2, M1 and M2), nucleotide sequences of the CDS identified in the RAST annotation were converted into a single concatenated FASTA file using Artemis . Sequences were aligned in Mauve genome aligner . Aligned sequence blocks were reassembled from the Mauve alignment. A multiFASTA file was made for each combination of the four genomes. DNAsp was used to identify the synonymous and non-synonymous positions and SNPs [16, 95]. The years since the isolates diverged was estimated as described by Foster et al. 2009 using the following formula :
Where X is the years since divergence, Ks is the proportion of synonymous SNPs to synonymous sites, Z is the mutation rate per generation and Y the number of generations per year. The mutation rate for S. Typhi has been estimated at 1.6 x 10-10 mutations per nucleotide per generation, calculated over 20,000 generations . This is very close to the calculation for the rate of mutation in E. coli of 1.4 x 10-10 per nucleotide per generation estimated over 20,000 generations . The rate of mutation has been shown to be highly variable and dependent on the mutation rate phenotype of the lineage . Here we use both rates calculated between S. Typhi generations as the lower limit and E. coli generations as the upper limit, with the assumption that the true rate sits somewhere between these two. This is based on the assumption that the variation in the mutation rate correlates with the phylogenetic relationship between the strains [100, 101]. The number of generations per year (Y) is taken from the estimate produced for a wild population of E. coli of between 100 and 300. The denominator is multiplied by two as it applies to the number of SNPs between two genomes . To estimate the time of divergence between S. Derby and S. Mbandaka an average was taken of the four possible Ks values for each of the four pair wise comparisons.
Visualisation of sequence architecture
Genomic maps were constructed using CIRCOS circular visualization of data tool v 0.56 . A program for calculating GC skew in R v2.11.0 using the library SeqinR v3.0-6 was modified from R graphical manual example “fragment of E. coli chromosome” [104–106]. The GC skew was calculated under a 1 kb window at a 200 bp interval. The RAST annotation files were de-constructed into four tracks, forward and reverse coding DNA and RNA. The SPI-23 comparison maps were constructed from modified DoubleACT outputs for each combination of S. Derby D1 SPI-23 and the genomes of S. Agona SL483. S. Dublin CT02021853 and S. Gallinarum RKS5078, with a 100 bp cut-off for width between non-homologus sequences.
Grimont and Weill: Antigenic Formulae of the Salmonella Serovars. 2007, 9
Kingsley R, Bäumler J: Host adaptation and the emergence of infectious disease: the Salmonella paradigm. Mol Microbiol. 2000, 36: 1006-1014. 10.1046/j.1365-2958.2000.01907.x.
AHVLA D: Salmonella in livestock production, Annual report. 2011, Weybridge UK: AHVLA
CDC PHLIS: Salmonella annual survey. 2009, Atlanta, USA: CDC, 2010
Pickard D, Wain J, Baker S: Composition, acquisition, and distribution of the Vi exopolysaccharide-encoding Salmonella enterica pathogenicity island SPI-7. J Bacteriol. 2003, 185: 5055-5065. 10.1128/JB.185.17.5055-5065.2003.
Amavisit P, Lightfoot D, Browning GF, Markham PF: Variation between Pathogenic Serovars within Salmonella Pathogenicity Islands. J Bacteriol. 2003, 185: 3624-3635. 10.1128/JB.185.12.3624-3635.2003.
Lu S, Killoran PB, Riley LW: Association of Salmonella enterica serovar Enteritidis YafD with resistance to chicken egg albumen. Infect Immun. 2003, 71: 6734-6741. 10.1128/IAI.71.12.6734-6741.2003.
Messens W, Duboccage L, Grijspeerdt K, Heyndrickx M, Herman L: Growth of Salmonella serovars in hens’ egg albumen as affected by storage prior to inoculation. Food Microbiol. 2004, 21: 25-32. 10.1016/S0740-0020(03)00045-5.
Betancor L, Pereira M, Martinez A, Giossa G, Fookes M, Flores K, Barrios P, Repiso V, Vignoli R, Cordeiro N, Algorta G, Thomson N, Maskell D, Schelotto F, Chabalgoity J: Prevalence of Salmonella enterica in poultry and eggs in Uruguay during an epidemic due to Salmonella enterica serovar Enteritidis. J Clin Microbiol. 2010, 48: 2413-2423. 10.1128/JCM.02137-09.
Kingsley RA, Humphries AD, Weening EH, De Zoete MR, Winter S, Papaconstantinopoulou A, Dougan G, Baumler AJ: Molecular and phenotypic analysis of the CS54 island of Salmonella enterica serotype Typhimurium: identification of intestinal colonization and persistence determinants. Infect Immun. 2003, 71: 629-640. 10.1128/IAI.71.2.629-640.2003.
Shah DH, Lee M, Park J, Lee J, Eo S, Kwon J, Chae J: Identification of Salmonella gallinarum virulence genes in a chicken infection model using PCR-based signature-tagged mutagenesis. Microbiology. 2005, 151: 3957-3968. 10.1099/mic.0.28126-0.
Hammer BK, Bassler BL: Quorum sensing controls biofilm formation in Vibrio cholerae. Mol Microbiol. 2003, 50: 101-104. 10.1046/j.1365-2958.2003.03688.x.
Toguchi A, Siano M, Burkart M, Harshey RM: Genetics of swarming motility in Salmonella enterica Serovar Typhimurium: critical role for lipopolysaccharide. J Bacteriol. 2000, 182: 6308-6321. 10.1128/JB.182.22.6308-6321.2000.
Yoon SH, Park YK, Lee S, Choi D, Oh TK, Hur CG, Kim JF: Towards pathogenomics: a web-based resource for pathogenicity islands. Nucleic Acids Res. 2007, 35: 395-400. 10.1093/nar/gkl790.
Jormakka M, Byrne B, Iwata S: Formate dehydrogenase – a versatile enzyme in changing environments. Curr Opin Struct Biol. 2003, 13: 418-423. 10.1016/S0959-440X(03)00098-8.
Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, Formsma K, Gerdes S, Glass EM, Kubal M, Meyer F, Olsen GJ, Olson R, Osterman AL, Overbeek RA, McNeil LK, Paarmann D, Paczian T, Parrello B, Pusch GD, Reich C, Stevens R, Vassieva O, Vonstein V, Wilke A, Zagnitko O: The RAST server: rapid annotations using subsystems technology. BMC Genomics. 2008, 9: 75-90. 10.1186/1471-2164-9-75.
Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang HY, Cohoon M, De Crecy-Lagard V, Diaz N, Disz T, Edwards R, Fonstein M, Frank ED, Gerdes S, Glass EM, Goesmann A, Hanson A, Iwata-Reuyl D, Jensen R, Jamshidi N, Krause L, Kubal M, Larsen N, Linke B, McHardy AC, Meyer F, Neuweger H, Olsen G, Olson R, Osterman A, Portnoy V: The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res. 2005, 33: 5691-5702. 10.1093/nar/gki866.
Barkay T, Miller SM, Summers AO: Bacterial mercury resistance from atoms to ecosystems. FEMS Microbiol Rev. 2003, 27: 355-384. 10.1016/S0168-6445(03)00046-9.
Alexander J, Auðunsson GA, Benford D, Cockburn A, Dogliotti E, Di DA, Férnandez-cruz ML, Fink-gremmels J, Galli CL, Grandjean P, Gzyl J, Heinemeyer G, Johansson N, Mutti A, Schlatter J, Van LR: Mercury as undesirable substance in animal feed 1 scientific opinion of the panel on contaminants in the food chain adopted on 20 February 2008. The EFSA journal. 2008, 654: 1-76.
Stevenson G, Neal B, Liu DAN, Hobbs M, Packer NH, Batley M, Redmond JW, Lindquist L, Reeves P: Structure of the O Antigen of Escherichia coli K-12 and the Sequence of Its rIb Gene Cluster. J Bacteriol. 1994, 176: 4144-4156.
Reuven NB, Tomer G, Livneh Z: The mutagenesis proteins UmuD and UmuC prevent lethal frameshifts while increasing base substitution mutations. Mol Cell. 1998, 2: 191-199. 10.1016/S1097-2765(00)80129-X.
Merritt ME, Donaldson JR: Effect of bile salts on the DNA and membrane integrity of enteric bacteria. J Med Microbiol. 2009, 58: 1533-1541. 10.1099/jmm.0.014092-0.
Sanders C, Turkarslan S, Lee DW, Daldal F: Cytochrome c biogenesis: the Ccm system. Trends Microbiol. 2010, 18: 266-274. 10.1016/j.tim.2010.03.006.
Reid E, Cole J, Eaves D: The Escherichia coli CcmG protein fulfils a specific role in cytochrome c assembly. Biochem J. 2001, 58: 51-58.
Ren Q, Ahuja U, Thöny-Meyer L: A bacterial cytochrome c heme lyase. CcmF forms a complex with the heme chaperone CcmE and CcmH but not with apocytochrome c. J Biol Chem. 2002, 277: 7657-7663. 10.1074/jbc.M110979200.
Sinha N, Ferguson SJ: An Escherichia coli ccm (cytochrome c maturation) deletion strain substantially expresses Hydrogenobacter thermophilus cytochrome c552 in the cytoplasm: availability of haem influences cytochrome c552 maturation. FEMS Microbiol Lett. 1998, 161: 1-6. 10.1111/j.1574-6968.1998.tb12921.x.
Manning PA, Stroeher UH, Karageorgos LE, Morona R: Putative O-antigen transport genes within the rfb region of Vibrio cholerae O1 are homologous to those for capsule transport. Gene. 1995, 158: 1-7. 10.1016/0378-1119(95)00124-O.
Xiang SH, Reeves PR, Haase a M: Variation of the rfb gene clusters in Salmonella enterica. J Bacteriol. 1993, 175: 4877-4884.
Chan K, Baker S, Kim CC, Detweiler CS, Dougan G, Falkow S: Genomic comparison of Salmonella enterica serovars and Salmonella bongori by use of an S. enterica serovar typhimurium DNA microarray. J Bacteriol. 2003, 185: 553-563. 10.1128/JB.185.2.553-563.2003.
Saroj SD, Shashidhar R, Karani M, Bandekar JR: Distribution of Salmonella pathogenicity island (SPI)-8 and SPI-10 among different serotypes of Salmonella. J Med Microbiol. 2008, 57: 424-427. 10.1099/jmm.0.47630-0.
Seth-Smith HMB: SPI-7: Salmonella’s Vi-Encoding Pathogenicity Island. J Infect Dev Ctries. 2008, 2: 267-271.
Ochman H, Groisman EA: Distribution of pathogenicity islands in Salmonella spp. Infect Immun. 1996, 64: 5410-5412.
Ochman H, Groisman E a: How Salmonella became a pathogen. Trends Microbiol. 1997, 5: 343-349. 10.1016/S0966-842X(97)01099-8.
Matsushita S, Yamada S, Sekiguchi K, Kusunoki J, Ohta K, Kudoh Y: Serovar-distribution and drug-resistance of Salmonella strains isolated from domestic and imported cases in 1990–1994 in Tokyo. JJAID. 1996, 70: 42-50.
Hu Q, Coburn B, Deng W, Li Y, Shi X, Lan Q, Wang B, Coombes BK, Finlay BB: Salmonella enterica serovar Senftenberg human clinical isolates lacking SPI-1. J Clin Microbiol. 2008, 46: 1330-1336. 10.1128/JCM.01255-07.
Fookes M, Schroeder GN, Langridge GC, Blondel CJ, Mammina C, Connor TR, Seth-Smith H, Vernikos GS, Robinson KS, Sanders M, Petty NK, Kingsley RA, Bäumler AJ, Nuccio SP, Contreras I, Santiviago CA, Maskell D, Barrow P, Humphrey T, Nastasi A, Roberts M, Frankel G, Parkhill J, Dougan G, Thomson NR: Salmonella bongori provides insights into the evolution of the Salmonellae. PLoS Pathog. 2011, 7: e1002191-10.1371/journal.ppat.1002191.
Blanc-Potard AB, Groisman EA: The Salmonella selC locus contains a pathogenicity island mediating intramacrophage survival. EMBO J. 1997, 16: 5376-5385. 10.1093/emboj/16.17.5376.
Marcus SL, Brumell JH, Pfeifer CG, Finlay BB: Salmonella pathogenicity islands: big virulence in small packages. Microbes Infect. 2000, 2: 145-156. 10.1016/S1286-4579(00)00273-2.
Hensel M: Evolution of pathogenicity islands of Salmonella enterica. Int J Med Microbiol. 2004, 294: 95-102. 10.1016/j.ijmm.2004.06.025.
McClelland M, Sanderson KE, Spieth J, Clifton SW, Latreille P, Courtney L, Porwollik S, Ali J, Dante M, Du F, Hou S, Layman D, Leonard S, Nguyen C, Scott K, Holmes A, Grewal N, Mulvaney E, Ryan E, Sun H, Florea L, Miller W, Stoneking T, Nhan M, Waterston R, Wilson RK: Complete genome sequence of Salmonella enterica serovar Typhimurium LT2. Nature. 2001, 413: 852-856. 10.1038/35101614.
Chiu C-H, Tang P, Chu C, Hu S, Bao Q, Yu J, Chou Y-Y, Wang H-S, Lee Y-S: The genome sequence of Salmonella enterica serovar Choleraesuis, a highly invasive and resistant zoonotic pathogen. Nucleic Acids Res. 2005, 33: 1690-1698. 10.1093/nar/gki297.
Amavisit P, Lightfoot D: Variation between pathogenic serovars within Salmonella pathogenicity islands. J Bacteriol. 2003, 185: 3624-3635. 10.1128/JB.185.12.3624-3635.2003.
Han X, Dorsey-Oresto A, Malik M, Wang JY, Drlica K, Zhao X, Lu T: Escherichia coli genes that reduce the lethal effects of stress. BMC Microbiol. 2010, 10: 35-44. 10.1186/1471-2180-10-35.
Raina S, Missiakas D, Baird L, Kumar S, Georgopoulos C: Identification and transcriptional analysis of the Escherichia coli htrE operon which is homologous to pap and related pilin operons. J Bacteriol. 1993, 175: 5009-5021.
Tenor JL, McCormick BA, Ausubel FM, Aballay A: Caenorhabditis elegans-based screen identifies Salmonella virulence factors required for conserved host-pathogen interactions. Curr Biol. 2004, 14: 1018-1024. 10.1016/j.cub.2004.05.050.
Parkhill J, Dougan G, James KD, Thomson NR, Pickard D, Wain J, Churcher C, Mungall KL, Bentley SD, Holden MT, Sebaihia M, Baker S, Basham D, Brooks K, Chillingworth T, Connerton P, Cronin A, Davis P, Davies RM, Dowd L, White N, Farrar J, Feltwell T, Hamlin N, Haque A, Hien TT, Holroyd S, Jagels K, Krogh A, Larsen TS: Complete genome sequence of a multiple drug resistant Salmonella enterica serovar Typhi CT18. Nature. 2001, 413: 848-852. 10.1038/35101607.
Johnston C, Pegues DA, Hueck CJ, Lee A, Miller SI: Transcriptional activation of Salmonella typhimurium invasion genes by a member of the phosphorylated response-regulator superfamily. Mol Microbiol. 1996, 22: 715-727. 10.1046/j.1365-2958.1996.d01-1719.x.
Teplitski M, Goodier RI, Ahmer BMM: Pathways leading from BarA/SirA to motility and virulence gene expression in Salmonella. J Bacteriol. 2003, 185: 7257-7265. 10.1128/JB.185.24.7257-7265.2003.
Klumpp J, Fuchs TM: Identification of novel genes in genomic islands that contribute to Salmonella typhimurium replication in macrophages. Microbiology. 2007, 153: 1207-1220. 10.1099/mic.0. 2006/004747-0.
Ong SY, Ng FL, Badai SS, Yuryev A, Alam M: Analysis and construction of pathogenicity island regulatory pathways in Salmonella enterica serovar Typhi. J Integr Bioinform. 2010, 7: 145-179.
Aschtgen M-S, Gavioli M, Dessen A, Lloubes R, Cascales E: The SciZ protein anchors the enteroaggrative Eschericia coli type VI secretion system to the cell wall. J Molecular Microbiology. 2010, 75: 886-899. 10.1111/j.1365-2958.2009.07028.x.
Folkesson A, Advani A, Sukupolvi S, Pfeifer JD, Normark S, Löfdahl S: Multiple insertions of fimbrial operons correlate with the evolution of Salmonella serovars responsible for human disease. Mol Microbiol. 1999, 33: 612-622. 10.1046/j.1365-2958.1999.01508.x.
Groisman EA, Sturmoski MA, Solomon FR, Lin R, Ochman H: Molecular, functional, and evolutionary analysis of sequences specific to Salmonella. Proc Natl Acad Sci USA. 1993, 90: 1033-1037. 10.1073/pnas.90.3.1033.
Waterfield NR, Daborn PJ, Ffrench-Constant RH: Genomic islands in Photorhabdus. Trends Microbiol. 2002, 10: 541-545. 10.1016/S0966-842X(02)02463-0.
Bishop A, House D, Perkins T, Baker S, Kingsley RA, Dougan G: Interaction of Salmonella enterica serovar Typhi with cultured epithelial cells: roles of surface structures in adhesion and invasion. Microbiology. 2008, 154: 1914-1926. 10.1099/mic.0.2008/016998-0.
Harris JB, Baresch-Bernal A, Rollins SM, Alam A, LaRocque RC, Bikowski M, Peppercorn AF, Handfield M, Hillman JD, Qadri F, Calderwood SB, Hohmann E, Breiman RF, Brooks WA, Ryan ET: Identification of In Vivo-Induced Bacterial Protein Antigens during Human Infection with Salmonella enterica Serovar Typhi. Infect Immun. 2006, 74: 5161-5168. 10.1128/IAI.00488-06.
Ghosh S, Chakraborty K, Nagaraja T, Basak S, Koley H, Dutta S, Mitra U, Das S: An adhesion protein of Salmonella enterica serovar Typhi is required for pathogenesis and potential target for vaccine development. Proc Natl Acad Sci. 2011, 108: 3348-3353. 10.1073/pnas.1016180108.
Barbe V, Vallenet D, Fonknechten N, Kreimeyer A, Oztas S, Labarre L, Cruveiller S, Robert C, Duprat S, Wincker P, Ornston LN, Weissenbach J, Marlière P, Cohen GN, Médigue C: Unique features revealed by the genome sequence of Acinetobacter sp. ADP1, a versatile and naturally transformation competent bacterium. Nucleic Acids Res. 2004, 32: 5766-5779. 10.1093/nar/gkh910.
Xu Z, Wang Y, Han Y, Chen J, Zhang XH: Mutation of a novel virulence-related gene mltD in Vibrio anguillarum enhances lethality in zebra fish. Res Microbiol. 2011, 162: 144-150. 10.1016/j.resmic.2010.08.003.
Poon WW, Davis DE, Ha HT, Jonassen T, Rather PN, Clarke CF: Identification of Escherichia coli ubiB, a gene required for the first monooxygenase step in ubiquinone biosynthesis. J Bacteriol. 2000, 182: 5139-5146. 10.1128/JB.182.18.5139-5146.2000.
Shabala L, Bowman J, Brown J, Ross T, McMeekin T, Shabala S: Ion transport and osmotic adjustment in Escherichia coli in response to ionic and non-ionic osmotica. Environ Microbiol. 2009, 11: 137-148. 10.1111/j.1462-2920.2008.01748.x.
Vernikos GS, Parkhill J: Interpolated variable order motifs for identification of horizontally acquired DNA: revisiting the Salmonella pathogenicity islands. Bioinformatics. 2006, 22: 2196-2203. 10.1093/bioinformatics/btl369.
Delcher AL, Harmon D, Kasif S, White O, Salzberg SL: Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999, 27: 4636-4641. 10.1093/nar/27.23.4636.
Ahmer BM, Tran M, Heffron F: The virulence plasmid of Salmonella typhimurium is self-transmissible. J Bacteriol. 1999, 181: 1364-1368.
Norris FA, Wilson MP, Wallis TS, Galyov EE, Majerus PW: SopB, a protein required for virulence of Salmonella dublin, is an inositol phosphate phosphatase. Proc Natl Acad Sci. 1998, 95: 14057-14059. 10.1073/pnas.95.24.14057.
Galyov EE, Wood MW, Rosqvist R, Mullan PB, Watson PR, Hedges S, Wallis TS: A secreted effector protein of Salmonella dublin is translocated into eukaryotic cells and mediates inflammation and fluid secretion in infected ileal mucosa. Mol Microbiol. 1997, 25: 903-912. 10.1111/j.1365-2958.1997.mmi525.x.
Gunn JS, Alpuche-Aranda CM, Loomis WP, Belden WJ, Miller SI: Characterization of the Salmonella typhimurium pagC/pagD chromosomal region. J Bacteriol. 1995, 177: 5040-5047.
Aschtgen M-S, Bernard CS, De Bentzmann S, Lloubès R, Cascales E: SciN is an outer membrane lipoprotein required for type VI secretion in enteroaggregative Escherichia coli. J Bacteriol. 2008, 190: 7523-7531. 10.1128/JB.00945-08.
McLennan AG: The Nudix hydrolase superfamily. Cell Mol Life Sci. 2006, 63: 123-143. 10.1007/s00018-005-5386-7.
Samudrala R, Heffron F, McDermott JE: Accurate prediction of secreted substrates and identification of a conserved putative secretion signal for type III secretion systems. PLoS Pathog. 2009, 5: e1000375-10.1371/journal.ppat.1000375.
Dean P: Functional domains and motifs of bacterial type III effector proteins and their roles in infection. FEMS Microbiol Rev. 2011, 35: 1100-1125. 10.1111/j.1574-6976.2011.00271.x.
Shao F, Merritt PM, Bao Z, Innes RW, Dixon JE: A Yersinia effector and a Pseudomonas avirulence protein define a family of cysteine proteases functioning in bacterial pathogenesis. Cell. 2002, 109: 575-588. 10.1016/S0092-8674(02)00766-3.
Boyd EF, Carpenter MR, Chowdhury N: Mobile effector proteins on phage genomes. Bacteriophage. 2012, 2: 139-148.
Canchaya C, Proux C, Fournous G, Bruttin A, Brüssow H: Prophage genomics. Microbiol Mol Biol Rev. 2003, 67: 238-276. 10.1128/MMBR.67.2.238-276.2003.
Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS: PHAST: a fast phage search tool. Nucleic Acids Res. 2011, 39: 347-352. 10.1093/nar/gkq749.
Miao EA, Miller SI: Bacteriophages in the evolution of pathogen-host interactions. Proc Natl Acad Sci. 1999, 96: 9452-9454. 10.1073/pnas.96.17.9452.
Brown NF, Coombes BK, Bishop JL, Wickham ME, Lowden MJ, Gal-Mor O, Goode DL, Boyle EC, Sanderson KL, Finlay BB: Salmonella phage ST64B encodes a member of the SseK/NleB effector family. PLoS One. 2011, 6: e17824-10.1371/journal.pone.0017824.
Fricke WF, Mammel MK, McDermott PF, Tartera C, White DG, Leclerc JE, Ravel J, Cebula TA: Comparative genomics of 28 Salmonella enterica isolates: evidence for CRISPR-mediated adaptive sublineage evolution. J Bacteriol. 2011, 193: 3556-3568. 10.1128/JB.00297-11.
Sorek R, Kunin V, Hugenholtz P: CRISPR–a widespread system that provides acquired resistance against phages in bacteria and archaea. Nat Rev Microbiol. 2008, 6: 181-186. 10.1038/nrmicro1793.
Clokie MR, Millard AD, Letarov AV, Heaphy S: Phages in nature. Bacteriophage. 2011, 1: 31-45. 10.4161/bact.1.1.14942.
Qu A, Brulc JM, Wilson MK, Law BF, Theoret JR, Joens LA, Konkel ME, Angly F, Dinsdale EA, Edwards RA, Nelson KE, White BA: Comparative metagenomics reveals host specific metavirulomes and horizontal gene transfer elements in the chicken cecum microbiome. PLoS One. 2008, 3: e2945-10.1371/journal.pone.0002945.
Corpet DE: Ecological factors influencing the transfer of plasmids in vitro and in vivo. J Antimicrob Chemother. 1986, 18: 127-132.
Rotger R, Casadesús J: The virulence plasmids of Salmonella. Int Microbiol. 1999, 2: 177-184.
Giuffra E, Kijas JMH, Amarger V, Carlborg O, Jeon JT, Andersson L: The origin of the domestic pig: independent domestication and subsequent introgression. Genetics. 2000, 154: 1785-1791.
Chaisson MJ, Pevzner PA: Short read fragment assembly of bacterial genomes. Genome Res. 2008, 18: 324-330. 10.1101/gr.7088808.
Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J: ACT: the Artemis comparison tool. Bioinformatics. 2005, 21: 3422-3423. 10.1093/bioinformatics/bti553.
Henry CS, DeJongh M, Best AA, Frybarger PM, Linsay B, Stevens RL: High-throughput generation, optimization and analysis of genome-scale metabolic models. Nat Biotechnol. 2010, 28: 977-982. 10.1038/nbt.1672.
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 1999, 27: 29-34. 10.1093/nar/27.1.29.
Blondel C, Jimenez J, Contreras I, Santiviago C: Comparative genomic analysis uncovers 3 novel loci encoding type six secretion systems differentially distributed in Salmonella serotypes. BMC Genomics. 2009, 10: 354-371. 10.1186/1471-2164-10-354.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.
Grissa I, Vergnaud G, Pourcel C: CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 2007, 35: 52-57.
Grissa I, Vergnaud G, Pourcel C: The CRISPRdb database and tools to display CRISPRs and to generate dictionaries of spacers and repeats. BMC Bioinformatics. 2007, 8: 172-182. 10.1186/1471-2105-8-172.
Darling ACE, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence With rearrangements. Genome Res. 2004, 14: 1394-1403. 10.1101/gr.2289704.
Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA, Barrell B: Artemis: sequence visualization and annotation. Bioinformatics. 2000, 16: 944-945. 10.1093/bioinformatics/16.10.944.
Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.
Foster JT, Beckstrom-Sternberg SM, Pearson T, Beckstrom-Sternberg JS, Chain PSG, Roberto FF, Hnath J, Brettin T, Keim P: Whole-genome-based phylogeny and divergence of the genus Brucella. J Bacteriol. 2009, 191: 2864-2870. 10.1128/JB.01581-08.
Barrick JE, Yu DS, Yoon SH, Jeong H, Oh TK, Schneider D, Lenski RE, Kim JF: Genome evolution and adaptation in a long-term experiment with Escherichia coli. Nature. 2009, 461: 1243-1247. 10.1038/nature08480.
Lenski RE, Winkworth CL, Riley MA: Rates of DNA sequence evolution in experimental populations of Escherichia coli during 20,000 generations. J Mol Evol. 2003, 56: 498-508. 10.1007/s00239-002-2423-0.
Van Cuyck H, Farbos-Granger A, Leroy P, Yith V, Guillard B, Sarthou JL, Koeck JL, Kruy SL: MLVA polymorphism of Salmonella enterica subspecies isolated from humans, animals, and food in Cambodia. BMC Res Notes. 2011, 4: 306-314. 10.1186/1756-0500-4-306.
Fukushima M, Kakinuma K, Kawaguchi R: Phylogenetic analysis of Salmonella, Shigella, and Escherichia coli strains on the basis of the gyrB gene sequence. J Clin Microbiol. 2002, 40: 2779-2785. 10.1128/JCM.40.8.2779-2785.2002.
Ochman H, Elwyn S, Moran NA: Calibrating bacterial evolution. Proc Natl Acad Sci. 1999, 96: 12638-12643. 10.1073/pnas.96.22.12638.
Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics. Genome Res. 2009, 19: 1639-1645. 10.1101/gr.092759.109.
fragment of E. coli chromosome. http://rgm3.lab.nig.ac.jp/RGM/r_function?p=seqinr&f=m16j,
R Core Team: R: A language and environment for statistical computing. 2012
Charif D, Lobry J: SeqinR 1.0-2: a contributed package to the R project for statistical computing devoted to biological sequences retrieval and analysis. Structural Approaches to Sequence Evolution. Edited by: Bastolla U, Porto M, Roman HE, Vendruscolo M. 2007, Springer Berlin Heidelberg, 207-232.
The project was funded through a four year EPSRC studentship to MRH and the AHVLA seedcorn internal investment fund. We gratefully acknowledge funding provided to VAAJ through grant nr. BB/I004548/1 by the Biotechnology and Biological Sciences Research Council partners of the ERASysBio+ initiative supported under the EU ERA-NET Plus scheme in FP7. With special thanks to colleagues in BAC2 in the Department of Bacteriology and Food Safety at AHVLA (Weybridge).
Genome annotations and sequences available under guest log in on the RAST website, accessible from: http://rast.nmpdr.org/ using the following RAST IDs: D1 [RAST: 28144.16], D2 [RAST: 28144.17], M1 [RAST: 192954.16] and M2 [RAST: 192954.17].
The authors declare that they have no competing interests.
MRH, VAAJ and MJW conceived the study and wrote the manuscript. DNA extractions, bioinformatics and analysis were performed by MRH. All authors read and approved the final draft.
Electronic supplementary material
Additional file 1: Differences in gene and phage complement, the genes flanking the inversion, and the CRISPR spacer sequences, between isolates D1, D2, M1 and M2.(XLS 103 KB)
About this article
Cite this article
Hayward, M.R., Jansen, V.A. & Woodward, M.J. Comparative genomics of Salmonella entericaserovars Derby and Mbandaka, two prevalent serovars associated with different livestock species in the UK. BMC Genomics 14, 365 (2013). https://doi.org/10.1186/1471-2164-14-365
- S. Derby
- S. Mbandaka
- Functional genomics
- Host adaptation