Multilocus Variable-Number Tandem-Repeat Method for Typing Salmonella enterica Serovar Newport

In recent years, the proportion of Salmonella enterica infections represented by S. enterica serovar Newport has increased markedly among humans and animals. Multilocus variable-number tandem-repeat analysis (MLVA) has proven to be useful in discriminating other highly clonal Salmonella serovars. Here, we report on the development of a highly discriminatory MLVA for Salmonella serovar Newport.Multilocus variable-number tandem-repeat analyses (MLVA) have been developed for subtyping bacterial pathogens, including Salmonella enterica (12, 13, 16, 20). MLVA detects polymorphisms at genomic loci that are composed of tandemly repeated DNA sequences. Pulsed-field gel electrophoresis (PFGE) has been used for many years to subtype Salmonella enterica with satisfactory results (7), but MLVA has the advantages of a shorter turnaround time, additional discriminatory power, and the possibility of providing phylogenetic information to investigators. We report here the development and application of an MLVA protocol for Salmonella enterica serovar Newport, a clinically and epidemiologically important serovar characterized by multidrug resistance (5, 6, 9, 11, 14, 19).Bovine isolates from the northwestern United States were obtained from previous field studies conducted at Washington State University (Pullman) and from clinical submissions to the Washington Animal Disease Diagnostic Laboratory (Pullman). Isolates from the northeastern United States originated from several independent cattle herd outbreaks and from epidemiologically independent sources (humans and wild birds) (see Table S1 in the supplemental material). Veterinary isolates were serotyped by the National Veterinary Services Laboratory (Ames, IA), and human isolates were serotyped by the state public health reference laboratories in Washington and New York State by using a standard protocol (15).Fourteen S. Newport isolates from diverse sources and with diverse PFGE profiles were initially assayed using a previously published protocol for S. Typhimurium (13). Two of the published loci, the STTR5 and STTR6 loci, were polymorphic (see Table S2 in the supplemental material). Tandem-repeat finder software (3) was used to identify potential variable-number tandem-repeat (VNTR) loci in the S. Newport genome (GenBank accession number {"type":"entrez-nucleotide","attrs":{"text":"NC_011080","term_id":"194442203","term_text":"NC_011080"}}NC_011080). Sixty-seven repeat loci were identified, including the STTR3, STTR5, STTR6, and STTR7 loci, previously reported by Lindstedt et al. for serovar Typhimurium (12, 13). Published primers for STTR7 (12) did not amplify a product, and STTR3 products were identical between S. Newport isolates. Among the remaining repeat loci, those with the highest number of repeat copies and the fewest mismatches in the repeated sequences were investigated further. Fourteen of the repeat loci met our criteria, and primers for these loci were designed using Primer3 (18), targeting an amplicon size between 150 and 600 bp (see Table S3 in the supplemental material). PCRs that produced products with visible fragment size differences between several strains of S. Newport (assessed using gel electrophoresis) were repeated with dye-labeled primers, and those products were submitted for fragment sizing by capillary electrophoresis (Laboratory for Biotechnology and Bioanalysis, Washington State University, Pullman). The Newport-A, Newport-B, Newport-L, and Newport-M loci were polymorphic in our set of test isolates, generated the expected product sizes, and could be assayed simultaneously (Table ​(Table11).

TABLE 1.

Primer sequences used for six VNTR target loci from Salmonella enterica serovar Newport

Target	Forward primera	Reverse primer	Repeat size (bp)b	Amplicon size (bp)
PCR 1
STTR6	5′-FAM-TCGGGCATGCGTTGAAA	5′-CTGGTGGGGAGAATGACTGG	6	397
Newport-A	5′-PET-ACTGAAAGGAAGGGGAGAGC	5′-GTCAGGGTGGAATAGAATGC	9	429
Newport-L	5′-VIC-GAAGTACCGAAGTGGGTGAT	5′-CGTCCGTTAGAGGAACGTAT	51	529
PCR 2
STTR5	5′-PET-ATGGCGAGGCGAGCAGCAGT	5′-GGTCAGGCCGAATAGCAGGAT	6	564
Newport-B	5′-VIC-GGCCGATATAGCTCAGTTGG	5′-GAACCTCGCTTAGGGTTGTG	12	350
Newport-M	5′-FAM-GGTCATAGAGGGTCTGCAT	5′-ATGGAGCACAGACCACTAAC	36	378

Open in a separate window^a6-Carboxyfluorescein (FAM), PET, and VIC are 5′-conjugated fluorescent dyes (Applied Biosystems).^bExpected repeat size for Salmonella enterica serovar Newport strain SL254 (GenBank accession number {"type":"entrez-nucleotide","attrs":{"text":"NC_011080","term_id":"194442203","term_text":"NC_011080"}}NC_011080).To avoid overlap in combinations of fragment sizes and dye colors, two separate triplex reactions were conducted using an iCycler thermal cycler (Bio-Rad, Hercules, CA) in 25-μl volumes (Table ​(Table1).1). All primers were purchased from Applied Biosystems (Foster City, CA). A colony boiled-lysate suspension was used as the PCR template, and cycling conditions for both reactions included an initial denaturation at 94°C for 15 min followed by 25 cycles of 94°C for 30 s, 55°C for 1 min, and 72°C for 1.5 min, with a final extension step at 72°C for 10 min. A size standard (LIZ600; Applied Biosystems) (0.125 μl) and Hi-Di formamide (19.375 μl) were added to the PCR products (0.5 μl) (total volume, 20 μl for capillary electrophoresis). Capillary electrophoresis was conducted using a 3730 DNA analyzer with a POP-7 polymer (Applied Biosystems). The resulting electropherograms were analyzed using GeneMarker software (Softgenetics LLC, State College, PA). Fragment sizes for each of the VNTR loci were analyzed as character data by BioNumerics software. Cluster analysis of the categorical coefficient was performed using the unweighted-pair group method using arithmetic averages (UPGMA) algorithm. The Simpson diversity index (SDI) and 95% confidence intervals (CI) around the SDI were calculated as previously described (8, 10).Isolates were also assayed using PFGE according to a previously published protocol (17). Bacterial DNA was digested in agarose plugs with the restriction enzyme XbaI, and the resulting DNA fragments were gel separated using a CHEF-DRII (Bio-Rad, Hercules, CA) apparatus. PFGE profiles were analyzed with BioNumerics software (Applied Maths, Sint-Martens-Latem, Belgium) by clustering Dice similarity coefficients, using the UPGMA algorithm with a 1% position tolerance and 0% optimization. Pearson''s correlation coefficient for the similarity matrices generated by MLVA and PFGE was 59.8%, as calculated using the “congruence of experiments” feature in BioNumerics.To determine the repeatability of MLVA results, 20 systematically chosen isolates were assayed twice. Among these, one replicate reaction produced an additional repeat length (STTR5 locus) compared with the initial PCR amplification result. Thus, for 120 individual allele determinations (20 isolates × 6 loci), 119 (99.2%) were stably repeated.MLVA divided 132 S. Newport isolates into 40 different types, where a type is defined as a set of alleles that is different from other isolates at one or more loci and where a null amplification product is considered a distinct allele after confirmation by a repeat assay (see Table S1 in the supplemental material). The number of isolates within each type ranged from 1 to 26, with an SDI of 0.91 (95% CI, 0.88 to 0.93). The SDI was 0.90 (95% CI, 0.86 to 0.95) for XbaI PFGE of the same group of isolates. The discriminatory power of MLVA relative to that of PFGE varied by region: for the northeastern isolates (Fig. ​(Fig.1),1), the SDIs (95% CI) for MLVA and PFGE were 0.90 (0.85 to 0.95) and 0.92 (0.87 to 0.97), respectively, while those for the northwestern isolates (Fig. ​(Fig.2)2) were 0.88 (0.83 to 0.93) for MLVA and 0.71 (0.61 to 0.82) for PFGE.Open in a separate windowFIG. 1.Dendrogram from cluster analysis of northeastern U.S. S. Newport isolates. The UPGMA algorithm was used to cluster categorical coefficients of MLVA allele data. STs were assigned by Alcaine et al. (1).Open in a separate windowFIG. 2.Dendrogram from cluster analysis of northwestern U.S. bovine S. Newport isolates. The UPGMA algorithm was used to cluster categorical coefficients of MLVA allele data.In the MLVA-generated cluster analysis of northeastern U.S. isolates, most (20 of 24) isolates from herd outbreaks of salmonellosis among cattle clustered with others from the same outbreak. The few exceptions, including isolates 16073 and 16055 from farm 510 and isolates 16071 and 16081 from farms 488 and 152, respectively, may have been misclassified by MLVA, but it is equally possible that they represent true dissemination of strains. MLVA clustering was consistent with multilocus sequence type (ST) data available from previous investigations (1, 2). Cluster 1 includes primarily ST 11 isolates, and cluster 2 includes human and avian isolates with diverse multiple MLVA types, PFGE types, and STs (Fig. ​(Fig.1).1). The cluster analysis of northwestern isolates from bovine sources provides additional support for the close correlation between our MLVA protocol and known epidemiologic data. As for the northeastern isolates, most isolates from the same farm or operation clustered together with few exceptions, and PFGE types corresponded well with MLVA types (Fig. ​(Fig.22).The MLVA protocol for S. Newport presented here demonstrated excellent repeatability and close congruence with PFGE and multilocus sequence typing, as well as having strong correspondence with epidemiological data. Our finding of occasional allele instability with this method is consistent with the instability of VNTR loci reported for Salmonella enterica serovar Typhimurium (4), and discrepancies between PFGE and MLVA may in part be explained by such instability (for example, isolates 10834 and 10835 in Fig. ​Fig.22 differed at a single locus and by one repeat length). However, given its relative stability and concordance with PFGE, the higher processing capacity and more-rapid turnaround of MLVA make it a strong complementary or alternative technique for regional surveillance and outbreak investigation.