Comparison of repetitive-sequence-based polymerase chain reaction with random amplified polymorphic DNA analysis for rapid genotyping of nontuberculosis mycobacteria

Nontuberculosis mycobacteria (NTM) are an important cause of human disease and infections. Though less notorious than tuberculosis, these infections are clinically significant and have been associated with outbreaks in various settings. To accommodate outbreak investigations for the numerous species of NTM, we evaluated a DiversiLab repetitive-sequence-based PCR (rep-PCR) kit for genotyping of mycobacteria. This kit was used to genotype both rapidly and slowly growing mycobacteria and was compared with other PCR-based genotyping methods, including random amplified polymorphic DNA (RAPD) analysis, hsp65 gene sequencing, and mycobacterial interspersed repetitive unit?- variable number of tandem repeat (MIRU-VNTR) analysis. Compared with RAPD analysis, rep-PCR achieved better reproducibility in testing. When compared with hsp65 gene sequencing and MIRU-VNTR for Mycobacterium avium , rep-PCR provided results that agreed with these less discriminatory genotyping methods but provided a higher level of discrimination for situations such as outbreak investigations. We also evaluated the kit for its ability to identify closely related rapidly growing NTM. While rep-PCR was informative in some cases, a much larger library of isolates would be necessary to truly evaluate it as an identification tool. Overall, rep-PCR was able to provide improved reproducibility over RAPD and a discriminatory genotyping method for the isolates evaluated in this study.     

Trinucleotide repeat sequence-based PCR as a potential approach for genotyping Mycobacterium gordonae strains

Diseases that are caused by non-tuberculous mycobacteria (NTM) continue to pose difficult clinical problems, and the epidemiological aspect of NTM-caused diseases is of great importance. In the case of Mycobacterium gordonae there is no adequate genotyping scheme. Here we present a potential rapid and reproducible genetic assay that uses trinucleotide repeat sequence-based PCR (TRS-PCR) for genotyping M. gordonae. The proposed method constitutes a useful single-primer PCR screen for genotyping this species. Among 10 TRS-containing primers, after applying (CAC)₄-based PCR to 36 strains of M. gordonae, we found a discriminatory index of 0.975. The accuracy of this analysis was supported by a reasonable reproducibility of 92%. These results were compared with the Enterobacterial Repetitive Intergenic Consensus Sequences (ERIC)-PCR typing scheme which had lower discriminatory index of 0.93 and its reproducibility was only 86.3%.     

Genotyping of Mycobacterium tuberculosis isolates from a low-endemic setting in northwestern state of Paraná in Southern Brazil

The purpose of this study was to provide information about the genetic diversity and prevalent genotype of Mycobacterium tuberculosis in a low-endemic setting in northwestern state of Paraná in Southern Brazil. We employed spoligotyping and mycobacterial interspersed repetitive units-variable number tandem repeat (MIRU-VNTR) techniques to genotype M. tuberculos isisolates from patients with pulmonary tuberculosis (TB). The 93 isolates analyzed by spoligotyping were divided into 36 different patterns, 30 of which were described in the SITVIT database. Latin American and Mediterranean, Haarlem and T families were responsible for 26.9%, 17.2% and 11.8% of TB cases, respectively. From the 84 isolates analyzed by MIRU-VNTR, 58 shared a unique pattern and the remaining 26 belonged to nine clusters. The MIRU loci 40, 23, 10 and 16 were the most discriminatory. A combination of MIRU-VNTR and spoligotyping resulted in 85.7% discriminatory power (Hunter-Gaston index = 0.995). Thus, combining spoligotyping and MIRU-VNTR typing proved to be most useful for epidemiological study in this low-endemic setting in Southern Brazil. The current study demonstrated that there is significant diversity in circulating strains in the city of Maringá and the surrounding regions, with no single genotype of M. tuberculosis predominating.     

Suspicion of Mycobacterium avium subsp. paratuberculosis transmission between cattle and wild-living red deer (Cervus elaphus) by multitarget genotyping

Multitarget genotyping of the etiologic agent Mycobacterium avium subsp. paratuberculosis is necessary for epidemiological tracing of paratuberculosis (Johne's disease). The study was undertaken to assess the informative value of different typing techniques and individual genome markers by investigation of M. avium subsp. paratuberculosis transmission between wild-living red deer and farmed cattle with known shared habitats. Fifty-three M. avium subsp. paratuberculosis type II isolates were differentiated by short sequence repeat analysis (SSR; 4 loci), mycobacterial interspersed repetitive-unit-variable-number tandem-repeat analysis (MIRU-VNTR; 8 loci), and restriction fragment length polymorphism analysis based on IS900 (IS900-RFLP) using BstEII and PstI digestion. Isolates originated from free-living red deer (Cervus elaphus) from Eifel National Park (n = 13), six cattle herds living in the area of this park (n = 23), and five cattle herds without any contact with these red deer (n = 17). Data based on individual herds and genotypes verified that SSR G2 repeats did not exhibit sufficient stability for epidemiological studies. Two common SSR profiles (without G2 repeats), nine MIRU-VNTR patterns, and nine IS900-RFLP patterns were detected, resulting in 17 genotypes when combined. A high genetic variability was found for red deer and cattle isolates within and outside Eifel National Park, but it was revealed only by combination of different typing techniques. Results imply that within this restricted area, wild-living and farmed animals maintain a reservoir for specific M. avium subsp. paratuberculosis genotypes. No host relation of genotypes was obtained. Results suggested that four genotypes had been transmitted between and within species and that one genotype had been transmitted between cattle herds only. Use of multitarget genotyping for M. avium subsp. paratuberculosis type II strains and sufficiently stable genetic markers is essential for reliable interpretations of epidemiological studies on paratuberculosis.     

Hypervariable loci that enhance the discriminatory ability of newly proposed 15-loci and 24-loci variable-number tandem repeat typing method on Mycobacterium tuberculosis strains predominated by the Beijing family

The newly proposed 15- and 24-loci mycobacterial interspersed repetitive unit (MIRU)-variable-number tandem repeat (VNTR) typing method was evaluated for its ability to differentiate 181 Mycobacterium tuberculosis Beijing family strains. Compared with the original 12-loci MIRU-VNTR typing method, the 15-loci system dramatically improved the discriminatory power for Beijing strains; however, large clusters that could be further differentiated by IS6110 restriction fragment length polymorphism (RFLP) were still obtained. The clonal stability and allelic diversity of a total of 31 VNTR loci were evaluated. VNTRs 3232, 3820, and 4120 were identified as the effective hypervariable VNTR set for the second-line typing of clustered strains following the 15-loci based scheme. Consequently, the discriminatory power of the new scheme (18 loci) equaled that of IS6110 RFLP.     

Spoligotyping and variable number tandem repeat analysis of Mycobacterium bovis isolates from cattle in Brazil

We performed spoligotyping and 12-mycobacterial interspersed repetitive unit-variable number tandem repeats (MIRU-VNTRs) typing to characterise Mycobacterium bovis isolates collected from tissue samples of bovines with lesions suggestive for tuberculosis during slaughter inspection procedures in abattoirs in Brazil. High-quality genotypes were obtained with both procedures for 61 isolates that were obtained from 185 bovine tissue samples and all of these isolates were identified as M. bovis by conventional identification procedures. On the basis of the spoligotyping, 53 isolates were grouped into nine clusters and the remaining eight isolates were unique types, resulting in 17 spoligotypes. The majority of the Brazilian M. bovis isolates displayed spoligotype patterns that have been previously observed in strains isolated from cattle in other countries. MIRU-VNTR typing produced 16 distinct genotypes, with 53 isolates forming eight of the groups, and individual isolates with unique VNTR profiles forming the remaining eight groups. The allelic diversity of each VNTR locus was calculated and only two of the 12-MIRU-VNTR loci presented scores with either a moderate (0.4, MIRU16) or high (0.6, MIRU26) discriminatory index (h). Both typing methods produced similar discriminatory indexes (spoligotyping h = 0.85; MIRU-VNTR h = 0.86) and the combination of the two methods increased the h value to 0.94, resulting in 29 distinct patterns. These results confirm that spoligotyping and VNTR analysis are valuable tools for studying the molecular epidemiology of M. bovis infections in Brazil.     

Genetic diversity of Mycobacterium tuberculosis isolates from Tibetans in Tibet, China

Background

Tuberculosis (TB) is a serious health problem in Tibet where Tibetans are the major ethnic group. Although genotyping of Mycobacterium tuberculosis (M. tuberculosis) isolates is a valuable tool for TB control, our knowledge of population structure of M. tuberculosis circulating in Tibet is limited.

Methodology/Principal Findings

In our study, a total of 576 M. tuberculosis isolates from Tibetans in Tibet, China, were analyzed via spoligotyping and 24-locus MIRU-VNTR. The Beijing genotype was the most prevalent family (90.63%, n = 522). Shared-type (ST) 1 was the most dominant genotype (88.89%, n = 512). We found that there was no association between the Beijing genotype and sex, age and treatment status. In this sample collection, 7 of the 24 MIRU-VNTR loci were highly or moderately discriminative according to their Hunter-Gaston discriminatory index. An informative set of 12 loci had similar discriminatory power with 24 loci set.

Conclusions/Significance

The population structure of M. tuberculosis isolates in Tibetans is homogeneous and dominated by Beijing genotype. The analysis of 24-locus MIRU-VNTR data might be useful to select appropriate VNTR loci for the genotyping of M. tuberculosis.     

Interleukin-6 is used as a growth factor by virulent Mycobacterium avium: presence of specific receptors.

In this paper, we examined the contribution of the lymphokine interleukin-6 (IL-6) to the growth of four virulent strains of Mycobacterium avium and the nature of the binding moieties on the mycobacteria. First, we showed that human or mouse recombinant interleukin-6 are potent growth factors for four strains of virulent M. avium. This was shown to occur in tissue culture medium, which does not support maximal growth of M. avium. Bioactive IL-6 was required, inasmuch as heat-activating IL-6 or adding an antibody against IL-6 blocked this growth-enhancing ability. The rapid uptake of IL-6 by M. avium was indicated by the fact that the incubation of IL-6 with the four M. avium strains led to a rapid removal of the bioactivity from the culture medium and a rapid removal of radiolabeled IL-6. Scatchard analysis of receptor interaction showed that the M. avium strains had a single receptor species with a Kd of 50 nM and the number of receptor sites was approximately 15,000 bacterium. Blocking experiments showed that the binding of radiolabeled IL-6 was fully displaceable with cold IL-6, but not with other lymphokines. These data suggest that IL-6 may play an important role in the pathogenesis of M. avium infections, notably by promoting growth of M. avium, and that some virulent M. avium strains bind IL-6 in a specific manner.     

Mycobacterium avium subsp. paratuberculosis and M. avium subsp. avium are independently evolved pathogenic clones of a much broader group of M. avium organisms

Mycobacterium avium comprises organisms that share the same species designation despite considerable genomic and phenotypic variability. To determine the degree and nature of variability between subspecies and strains of M. avium, we used multilocus sequencing analysis, studying 56 genetically diverse strains of M. avium that included all described subspecies. In total, 8,064 bp of sequence from 10 gene loci were studied, with 205 (2.5%) representing variable positions. The majority (149/205) of these variations were found among M. avium subsp. hominissuis organisms. Recombination was also evident in this subspecies. In contrast, there was comparatively little variability and no evidence of recombination within the pathogenic subspecies, M. avium subsp. paratuberculosis, M. avium subsp. avium, and M. avium subsp. silvaticum. Phylogenetic analysis showed that M. avium subsp. avium and M. avium subsp. silvaticum strains clustered together on one branch, while a distinct branch defined M. avium subsp. paratuberculosis organisms. Despite the independent origin of these pathogenic subspecies, an analysis of their rates of nonsynonymous (dN) to synonymous (dS) substitutions showed increased dN/dS ratios for both: 0.67 for M. avium subsp. paratuberculosis and 0.50 for M. avium subsp. avium/M. avium subsp. silvaticum, while the value was 0.08 for M. avium subsp. hominissuis organisms. In conclusion, M. avium subsp. hominissuis represents a diverse group of organisms from which two pathogenic clones (M. avium subsp. paratuberculosis and M. avium subsp. avium/M. avium subsp. silvaticum) have evolved independently.     

Fluorescence in situ hybridization using peptide nucleic acid probes for rapid detection of Mycobacterium avium subsp. avium and Mycobacterium avium subsp. paratuberculosis in potable-water biofilms

Here, we present for the first time a high-affinity peptide nucleic acid (PNA) oligonucleotide sequence for detecting Mycobacterium avium bacteria, including the opportunistically pathogenic subspecies M. avium subsp. avium, M. avium subsp. paratuberculosis, and M. avium subsp. silvaticum, by the fluorescence in situ hybridization (FISH) method. There is evidence that M. avium subsp. avium especially is able to survive and grow in drinking-water biofilms and possibly transmit via drinking water. The designed PNA probe (MAV148) specificity was tested with several bacterial species, including other mycobacteria and mycolic acid-containing bacteria. From the range of bacterial strains tested, only M. avium subsp. avium and M. avium subsp. paratuberculosis strains were hybridized. The PNA FISH method was applied successfully to detect M. avium subsp. avium spiked in water samples and biofilm established within a Propella biofilm reactor fed with potable water from a distribution supply.     

Characterization of <Emphasis Type="Italic">Mycobacterium tuberculosis</Emphasis> Central Asian Strain1 using mycobacterial interspersed repetitive unit genotyping

Background  

The Central Asian Strain1 (CAS1) genogroup of Mycobacterium tuberculosis (MTB) is the most prevalent in Pakistan, India and Bangladesh. Mycobacterial interspersed repetitive units variable number tandem repeat (MIRU-VNTR) typing is a reliable and reproducible method for differentiation of MTB isolates. However, information of its utility in determining the diversity of CAS1 strain is limited. We performed standard 12 loci based MIRU-VNTR typing on previously spoligotyped CAS1 strains and 'unique' strains in order to evaluate its discriminatory power for these isolates.     

Genotypic characteristics of Mycobacterium avium subsp. hominissuis strains

Mycobacterium avium are typical environmental, non-tuberculosis microorganisms that occasionally cause mycobacteriosis, an infectious disease in wild and domestic animals, birds, and humans. Here, we report the results of the first study on the genetic diversity of the Russian population of M. avium. A total of 85 M. avium subsp. hominissuis (MAH) clinical strains were isolated from patients (including 30 HIV-positive individuals) with mycobacteriosis in St. Petersburg, 2008–2011. The identification of the microorganisms was carried out using biochemical tests and the PCR detection of the mobile elements IS901 and IS900, as well as of the polymorphism of restriction fragments of the hsp65 gene. The genetic diversity of the isolates was evaluated by VNTR typing based on eight variable-number tandem repeats (VNTRs) (292, X3, 25, 47, 3, 7, 10, and 32 [Thibault et al., 2007]). The MAH population studied was characterized by 15 VNTR types, including nine unique patterns and six clusters of isolates with identical eight-digit profiles. The largest clusters (22221128 and 24221128) included 45 (59.2%) and 15 (19.7%) isolates, respectively; the others contained 2–7 strains. The strains of the cluster 2533112’8 possessed a truncated TR10 locus (allele 2′). Taking into account the absence of the epidemiological links between the patients and the fact that the infection was presumably delivered from the environment, the high rate of clustering of MAH isolates can be explained by the low discriminatory power of the eight-locus VNTR-typing scheme (HGDI 0–0.61).     

Genetic diversity of the Mycobacterium tuberculosis Beijing family based on SNP and VNTR typing profiles in Asian countries

The Mycobacterium tuberculosis (MTB) Beijing strain is highly virulent, drug resistant, and endemic over Asia. To explore the genetic diversity of this family in several different regions of eastern Asia, 338 Beijing strains collected in Taiwan (Republic of China) were analyzed by mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing and compared with published MIRU-VNTR profiles and by the Hunter-Gaston diversity index (HGDI) of Beijing strains from Japan and South Korea. The results revealed that VNTR2163b (HGDI>0.6) and five other loci (VNTR424, VNTR4052, VNTR1955, VNTR4156 and VNTR 2996; HGDI>0.3) could be used to discriminate the Beijing strains in a given geographic region. Analysis based on the number of VNTR repeats showed three VNTRs (VNTR424, 3192, and 1955) to be phylogenetically informative loci. In addition, to determine the geographic variation of sequence types in MTB populations, we also compared sequence type (ST) data of our strains with published ST profiles of Beijing strains from Japan and Thailand. ST10, ST22, and ST19 were found to be prevalent in Taiwan (82%) and Thailand (92%). Furthermore, classification of Beijing sublineages as ancient or modern in Taiwan was found to depend on the repeat number of VNTR424. Finally, phylogenetic relationships of MTB isolates in Taiwan, South Korea, and Japan were revealed by a minimum spanning tree based on MIRU-VNTR genotyping. In this topology, the MIRU-VNTR genotypes of the respective clusters were tightly correlated to other genotypic characters. These results are consistent with the hypothesis that clonal evolution of these MTB lineages has occurred.     

Genetic diversity analysis of isolates belonging to Bordetella genus

In this study, RandomAmplified Polymorphic DNA (RAPD) method was used to track differences among human and animals isolates of eight known species of Bordetella genus. One hundred representative strains of these species from international and Polish bacterial collections were genotyped according to RAPD protocols using primer 1254 or 1247. Problems with reproducibility and discriminatory power, frequently discussed in literature, have been overcome by precise optimization procedure, which allowed to achieve reliable conditions for Bordetella species analysed. This study proved high potential of RAPD method using primers 1247/1254 for intra-species differentiation of B. bronchiseptica and B. hinzii isolates. Additionally RAPD method has been found to be useful for confirming B. avium and B. holmesii identification.     

A Novel Molecular Strategy for Surveillance of Multidrug Resistant Tuberculosis in High Burden Settings

Background

In South Africa and other high prevalence countries, transmission is a significant contributor to rising rates of multidrug resistant tuberculosis (MDR-TB). Thus, there is a need to develop an early detection system for transmission clusters suitable for high burden settings. We have evaluated the discriminatory power and clustering concordance of a novel and simple genotyping approach, combining spoligotyping with pncA sequencing (SpoNC), against two well-established methods: IS6110-RFLP and 24-loci MIRU-VNTR.

Methods

A total of 216 MDR-TB isolates collected from January to June 2010 from the NHLS Central TB referral laboratory in Braamfontein, Johannesburg, representing a diversity of strains from South Africa, were included. The isolates were submitted for genotyping, pncA sequencing and analysis to the Centre for Tuberculosis in South Africa and the Public Health Research Institute Tuberculosis Center at Rutgers University in the United States. Clustering rates, Hunter-Gaston Discriminatory Indexes (HGI) and Wallace coefficients were compared between the methods.

Results

Overall clustering rates were high by both IS6110-RFLP (52.8%) and MIRU-VNTR (45.8%), indicative of on-going transmission. Both 24-loci MIRU-VNTR and IS6110-RFLP had similar HGI (0.972 and 0.973, respectively), with close numbers of unique profiles (87 vs. 70), clustered isolates (129 vs. 146), and cluster sizes (2 to 26 vs. 2 to 25 isolates). Spoligotyping alone was the least discriminatory (80.1% clustering, HGI 0.903), with 28 unique types. However, the discriminatory power of spoligotyping was improved when combined with pncA sequencing using the SpoNC approach (61.8% clustering, HGI 0.958). A high proportion of MDR-TB isolates had mutations in pncA (68%, n = 145), and pncA mutations were significantly associated with clustering (p = 0.007 and p = 0.0013 by 24-loci MIRU-VNTR and IS6110-RFLP, respectively), suggesting high rates of resistance to pyrazinamide among all MDR-TB cases and particularly among clustered cases.

Conclusion

We conclude that SpoNC provides good discrimination for MDR-TB surveillance and early identification of outbreaks in South Africa, with 24-loci MIRU-VNTR applied for pncA wild-type strains as needed.     

Value of morphotyping for the characterization of Candida albicans clinical isolates

Until recently, morphotyping, a method evaluating fringe and surface characteristics of streak colonies grown on malt agar, has been recommended as a simple and unexpensive typing method for Candida albicans isolates. The discriminatory power and reproducibility of Hunter's modified scheme of Phongpaichit's morphotyping has been evaluated on 28 C. albicans isolates recovered from the oral cavity of asymptomatic human immunodeficiency virus-positive subjects, and compared to two molecular typing methods: randomly amplified polymorphic DNA (RAPD) fingerprinting, and contour clamped homogeneous electric field (CHEF) electrophoretic karyotyping. Morphological features of streak colonies allowed to distinguish 11 different morphotypes while RAPD fingerprinting yielded 25 different patterns and CHEF electrophoresis recognized 9 karyotypes. The discriminatory power calculated with the formula of Hunter and Gaston was 0.780 for morphotyping, 0.984 for RAPD fingerprinting, and 0.630 for karyotyping. Reproducibility was tested using 43 serial isolates from 15 subjects (2 to 6 isolates per subject) and by repeating the test after one year storage of the isolates. While genetic methods generally recognized a single type for all serial isolates from each of the subjects studied, morphotyping detected strain variations in five subjects in the absence of genetic confirmation. Poor reproducibility was demonstrated repeating morphotyping after one year storage of the isolates since differences in at least one character were detected in 92.9% of the strains.     

Mycobacterial Agglutination-Inhibition Test

The mycobacterial agglutination-inhibition test (MAIT) is described as a modification of the mycobacterial agglutination of Schaefer (MATS). It is an organized assemblage of absorption and agglutination tests which is especially useful in the identification of Mycobacterium avium isolates normally excluded from serological identification because of unstable cell suspension characteristics. The identity of 100 unstable M. avium and 3 unstable Mycobacterium intracellulare isolates was determined by the MAIT method and by animal pathogenicity tests. The comparability of the MAIT and MATS methods was demonstrated by testing 10 stable M. avium and 10 stable M. intracellulare isolates by both methods and obtaining the same serological identification for each culture. The reproducibility of the MAIT was confirmed when the same result was obtained in three consecutive tests on each of the 123 cultures used in this evaluation.     

Comparison of Variable-Number Tandem-Repeat Markers typing and IS1245 Restriction Fragment Length Polymorphism fingerprinting of Mycobacterium avium subsp. hominissuis from human and porcine origins

Background  

Animal mycobacterioses are regarded as a potential zoonotic risk and cause economical losses world wide. M. avium subsp. hominissuis is a slow-growing subspecies found in mycobacterial infected humans and pigs and therefore rapid and discriminatory typing methods are needed for epidemiological studies. The genetic similarity of M. avium subsp. hominissuis from human and porcine origins using two different typing methods have not been studied earlier. The objective of this study was to compare the IS1245 RFLP pattern and MIRU-VNTR typing to study the genetic relatedness of M. avium strains isolated from slaughter pigs and humans in Finland with regard to public health aspects.     

MIRU-VNTR genotyping of Mycobacterium tuberculosis strains from East Siberian: Beijing family versus Kilimanjaro family

Seventy-two M. tuberculosis strains isolated from patients in different regions of East Siberia were studied by MIRU-VNTR repeats. Phylogenetic trees of relationship between the studied strains and of those from other world regions were investigated by phylogenetic analysis, i.e. NJ, UPGMA and ML. Phylogenetic evidence is provided for the prevalence of Beijing and Beijing-like strains (71%) in the above region. Besides, there was found a group of strains (115%) with the robust genotype which was given the name of ASU (African-Siberian-Ural). A linkage between MIRU10 (5-10 repeat) and MIRU26 (single repeat) was established as a predominant ASU sign. According to preliminary results, apart from Siberia, ASU strains can be encountered in the Urals and Africa (Tanzania, Kilimanjaro genotype). MIRU-VNTR consensus of ASU for East Siberia is 22( > 6)225113233.     

Multi-virulence-locus sequence typing of Listeria monocytogenes

A multi-virulence-locus sequence typing (MVLST) scheme was developed for subtyping Listeria monocytogenes, and the results obtained using this scheme were compared to those of pulsed-field gel electrophoresis (PFGE) and the published results of other typing methods, including ribotyping (RT) and multilocus sequence typing (MLST). A set of 28 strains (eight different serotypes and three known genetic lineages) of L. monocytogenes was selected from a strain collection (n > 1,000 strains) to represent the genetic diversity of this species. Internal fragments (ca. 418 to 469 bp) of three virulence genes (prfA, inlB, and inlC) and three virulence-associated genes (dal, lisR, and clpP) were sequenced and analyzed. Multiple DNA sequence alignment identified 10 (prfA), 19 (inlB), 13 (dal), 10 (lisR), 17 (inlC), and 16 (clpP) allelic types and a total of 28 unique sequence types. Comparison of MVLST with automated EcoRI-RT and PFGE with ApaI enzymatic digestion showed that MVLST was able to differentiate strains that were indistinguishable by RT (13 ribotypes; discrimination index = 0.921) or PFGE (22 profiles; discrimination index = 0.970). Comparison of MVLST with housekeeping-gene-based MLST analysis showed that MVLST provided higher discriminatory power for serotype 1/2a and 4b strains than MLST. Cluster analysis based on the intragenic sequences of the selected virulence genes indicated a strain phylogeny closely related to serotypes and genetic lineages. In conclusion, MVLST may improve the discriminatory power of MLST and provide a convenient tool for studying the local epidemiology of L. monocytogenes.