Multiple-locus variable number tandem repeat analysis for Streptococcus pneumoniae: comparison with PFGE and MLST

In the era of pneumococcal conjugate vaccines, surveillance of pneumococcal disease and carriage remains of utmost importance as important changes may occur in the population. To monitor these alterations reliable genotyping methods are required for large-scale applications. We introduced a high throughput multiple-locus variable number tandem repeat analysis (MLVA) and compared this method with pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST). The MLVA described here is based on 8 BOX loci that are amplified in two multiplex PCRs. The labeled PCR products are sized on an automated DNA sequencer to accurately determine the number of tandem repeats. The composite of the number of repeats of the BOX loci makes up a numerical profile that is used for identification and clustering. In this study, MLVA was performed on 263 carriage isolates that were previously characterized by MLST and PFGE. MLVA, MLST and PFGE (cut-off of 80%) yielded 164, 120, and 87 types, respectively. The three typing methods had Simpson's diversity indices of 98.5% or higher. Congruence between MLST and MLVA was high. The Wallace of MLVA to MLST was 0.874, meaning that if two strains had the same MLVA type they had an 88% chance of having the same MLST type. Furthermore, the Wallace of MLVA to clonal complex of MLST was even higher: 99.5%. For some isolates belonging to a single MLST clonal complex although displaying different serotypes, MLVA was more discriminatory, generating groups according to serotype or serogroup. Overall, MLVA is a promising genotyping method that is easy to perform and a relatively cheap alternative to PFGE and MLST. In the companion paper published simultaneously in this issue we applied the MLVA to assess the pneumococcal population structure of isolates causing invasive disease in The Netherlands before the introduction of the 7-valent conjugate vaccine.     

Multiple-locus variable-number tandem repeat analysis of Salmonella Enteritidis isolates from human and non-human sources using a single multiplex PCR

Simplified multiple-locus variable-number tandem repeat analysis (MLVA) was developed using one-shot multiplex PCR for seven variable-number tandem repeats (VNTR) markers with high diversity capacity. MLVA, phage typing, and PFGE methods were applied on 34 diverse Salmonella Enteritidis isolates from human and non-human sources. MLVA detected allelic variations that helped to classify the S. Enteritidis isolates into more evenly distributed subtypes than other methods. MLVA-based S. Enteritidis clonal groups were largely associated with sources of the isolates. Nei's diversity indices for polymorphism ranged from 0.25 to 0.70 for seven VNTR loci markers. Based on Simpson's and Shannon's diversity indices, MLVA had a higher discriminatory power than pulsed field gel electrophoresis (PFGE), phage typing, or multilocus enzyme electrophoresis. Therefore, MLVA may be used along with PFGE to enhance the effectiveness of the molecular epidemiologic investigation of S. Enteritidis infections.     

Multiple-locus variable-number tandem repeat analysis for strain typing of Clostridium perfringens

Clostridium perfringens is ubiquitous in the environment and causes diseases in man and animals, with syndromes ranging from enteritis, enterotoxemia, and sudden death to food poisoning and gas gangrene. Understanding the epidemiology of these infections and of the evolution of virulence in C. perfringens necessitate an efficient, time and cost effective strain typing method. Multiple-locus variable-number tandem repeat analysis (MLVA) has been applied to typing of other pathogens and we describe here the development of a MLVA scheme for C. perfringens. We characterized five variable tandem repeat (VNTR) loci, four of which are contained within protein encoding genes and screened 112 C. perfringens isolates to evaluate typability, reproducibility, and discriminatory power of the scheme. All the isolates were assigned a MLVA genotype and the technique has excellent reproducibility, with a numerical index of discrimination for the five VNTR loci of 0.995. Thus MLVA is an efficient tool for C. perfringens strain typing, and being PCR based makes it rapid, easy, and cost effective. In addition, it can be employed in epidemiological, ecological, and evolutionary investigations of the organism.     

Multilocus variable number tandem repeat analysis as a tool to discern genetic relationships among strains of Yersinia enterocolitica biovar 1A

Aims:  To identify variable number tandem repeat (VNTR)-containing loci, and to use multilocus VNTR (MLVA) to discern genetic relationships among strains of Yersinia enterocolitica biovar 1A isolated from diverse sources.
Methods and Results:  The whole genome sequence of Y. enterocolitica 8081 was analysed and eight VNTR loci with repeat sizes between 4 and 9 bp, and each containing more than four repeat copies were selected for MLVA typing of 88 strains of Y. enterocolitica . Of these, four loci were polymorphic and generated 26 MLVA genotypes among 81 strains of Y. enterocolitica biovar 1A. MLVA was found to be quite discriminatory (DI = 0·87). Cluster analysis and population modelling using minimum spanning tree (MST) clearly clustered Y. enterocolitica biovar 1A into two major groups.
Conclusions:  The MLVA is easy to perform and can be used to discern clonal relationship among strains of Y. enterocolitica . Also the phylogenetic relationships obtained with MLVA genotypes were in good agreement with those established by other typing methods.
Significance and Impact of the Study:  The MLVA method reported is relatively more discriminatory than the other genotyping methods and has the potential to be used as an epidemiological tool for the study of strains of Y. enterocolitica biovar 1A.     

Four Genotyping Schemes for Phylogenetic Analysis of Pseudomonas aeruginosa: Comparison of Their Congruence with Multi-Locus Sequence Typing

Several molecular typing schemes have been proposed to differentiate among isolates and clonal groups, and hence establish epidemiological or phylogenetic links. It has been widely accepted that multi-locus sequence typing (MLST) is the gold standard for phylogenetic typing/long-term epidemiological surveillance, but other recently described methods may be easier to carry out, especially in settings with limited access to DNA sequencing. Comparing the performance of such techniques to MLST is therefore of relevance. A study was therefore carried out with a collection of P. aeruginosa strains (n = 133) typed by four typing schemes: MLST, multiple-locus variable number tandem repeat analysis (MLVA), pulsed-field gel electrophoresis (PFGE) and the commercial DiversiLab microbial typing system (DL). The aim of this study was to compare the results of each typing method with MLST. The Simpson''s indices of diversity were 0.989, 0.980, 0.961 and 0.906 respectively for PFGE, MLVA, DL and MLST. The congruence between techniques was measured by the adjusted Wallace index (W): this coefficient indicates the probability that a pair of isolates which is assigned to the same type by one typing method is also typed as identical by the other. In this context, the congruence between techniques was recorded as follow: MLVA-type to predict MLST-type (93%), PFGE to MLST (92%), DL to MLST (64.2%), PFGE to MLVA (63.5%) and PFGE to DL (61.7%). Conversely, for all above combinations, prediction was very poor. The congruence was increased at the clonal complex (CC) level. MLST is regarded the gold standard for phylogenetic classification of bacteria, but is rather laborious to carry out in many settings. Our data suggest that MLVA can predict the MLST-type with high accuracy, and even higher when studying the clonal complex level. Of the studied three techniques MLVA was therefore the best surrogate method to predict MLST.     

Within-host competition between Borrelia afzelii ospC strains in wild hosts as revealed by massively parallel amplicon sequencing

Infections frequently consist of more than one strain of a given pathogen. Experiments have shown that co-infecting strains often compete, so that the infection intensity of each strain in mixed infections is lower than in single strain infections. Such within-host competition can have important epidemiological and evolutionary consequences. However, the extent of competition has rarely been investigated in wild, naturally infected hosts, where there is noise in the form of varying inoculation doses, asynchronous infections and host heterogeneity, which can potentially alleviate or eliminate competition. Here, we investigated the extent of competition between Borrelia afzelii strains (as determined by ospC genotype) in three host species sampled in the wild. For this purpose, we developed a protocol for 454 amplicon sequencing of ospC, which allows both detection and quantification of each individual strain in an infection. Each host individual was infected with one to six ospC strains. The infection intensity of each strain was lower in mixed infections than in single ones, showing that there was competition. Rank-abundance plots revealed that there was typically one dominant strain, but that the evenness of the relative infection intensity of the different strains in an infection increased with the multiplicity of infection. We conclude that within-host competition can play an important role under natural conditions despite many potential sources of noise, and that quantification by next-generation amplicon sequencing offers new possibilities to dissect within-host interactions in naturally infected hosts.     

Multiple-locus variable-number tandem repeat analysis of Dutch Bordetella pertussis strains reveals rapid genetic changes with clonal expansion during the late 1990s

Bordetella pertussis, the causative agent of whooping cough, has remained endemic in The Netherlands despite extensive nationwide vaccination since 1953. In the 1990s, several epidemic periods have resulted in many cases of pertussis. We have proposed that strain variation has played a major role in the upsurges of this disease in The Netherlands. Therefore, molecular characterization of strains is important in identifying the causes of pertussis epidemiology. For this reason, we have developed a multiple-locus variable-number tandem repeat analysis (MLVA) typing system for B. pertussis. By combining the MLVA profile with the allelic profile based on multiple-antigen sequence typing, we were able to further differentiate strains. The relationships between the various genotypes were visualized by constructing a minimum spanning tree. MLVA of Dutch strains of B. pertussis revealed that the genotypes of the strains isolated in the prevaccination period were diverse and clearly distinct from the strains isolated in the 1990s. Furthermore, there was a decrease in diversity in the strains from the late 1990s, with a remarkable clonal expansion that coincided with the epidemic periods. Using this genotyping, we have been able to show that B. pertussis is much more dynamic than expected.     

Genotypes of pathogenic Leptospira spp isolated from rodents in Argentina

Leptospirosis is the most widespread zoonosis in the world and significant efforts have been made to determine and classify pathogenic Leptospira strains. This zoonosis is maintained in nature through chronic renal infections of carrier animals, with rodents and other small mammals serving as the most important reservoirs. Additionally, domestic animals, such as livestock and dogs, are significant sources of human infection. In this study, a multiple-locus variable-number tandem repeat analysis (MLVA) was applied to genotype 22 pathogenic Leptospira strains isolated from urban and periurban rodent populations from different regions of Argentina. Three MLVA profiles were identified in strains belonging to the species Leptospira interrogans (serovars Icterohaemorrhagiae and Canicola); one profile was observed in serovar Icterohaemorrhagiae and two MLVA profiles were observed in isolates of serovars Canicola and Portlandvere. All strains belonging to Leptospira borgpetersenii serovar Castellonis exhibited the same MLVA profile. Four different genotypes were isolated from urban populations of rodents, including both mice and rats and two different genotypes were isolated from periurban populations.     

Diversification and Distribution of Ruminant Chlamydia abortus Clones Assessed by MLST and MLVA

Chlamydia abortus, an obligate intracellular bacterium, is the most common infectious cause of abortion in small ruminants worldwide and has zoonotic potential. We applied multilocus sequence typing (MLST) together with multiple-locus variable-number tandem repeat analysis (MLVA) to genotype 94 ruminant C. abortus strains, field isolates and samples collected from 1950 to 2011 in diverse geographic locations, with the aim of delineating C. abortus lineages and clones. MLST revealed the previously identified sequence types (STs) ST19, ST25, ST29 and ST30, plus ST86, a recently-assigned type on the Chlamydiales MLST website and ST87, a novel type harbouring the hemN_21 allele, whereas MLVA recognized seven types (MT1 to MT7). Minimum-spanning-tree analysis suggested that all STs but one (ST30) belonged to a single clonal complex, possibly reflecting the short evolutionary timescale over which the predicted ancestor (ST19) has diversified into three single-locus variants (ST86, ST87 and ST29) and further, through ST86 diversification, into one double-locus variant (ST25). ST descendants have probably arisen through a point mutation evolution mode. Interestingly, MLVA showed that in the ST19 population there was a greater genetic diversity than in other STs, most of which exhibited the same MT over time and geographical distribution. However, the evolutionary pathways of C. abortus STs seem to be diverse across geographic distances with individual STs restricted to particular geographic locations. The ST30 singleton clone displaying geographic specificity and represented by the Greek strains LLG and POS was effectively distinguished from the clonal complex lineage, supporting the notion that possibly two separate host adaptations and hence independent bottlenecks of C. abortus have occurred through time. The combination of MLST and MLVA assays provides an additional level of C. abortus discrimination and may prove useful for the investigation and surveillance of emergent C. abortus clonal populations.     

Ruminant rhombencephalitis-associated Listeria monocytogenes alleles linked to a multilocus variable-number tandem-repeat analysis complex

Listeria monocytogenes is among the most important food-borne pathogens and is well adapted to persist in the environment. To gain insight into the genetic relatedness and potential virulence of L. monocytogenes strains causing central nervous system (CNS) infections, we used multilocus variable-number tandem-repeat analysis (MLVA) to subtype 183 L. monocytogenes isolates, most from ruminant rhombencephalitis and some from human patients, food, and the environment. Allelic-profile-based comparisons grouped L. monocytogenes strains mainly into three clonal complexes and linked single-locus variants (SLVs). Clonal complex A essentially consisted of isolates from human and ruminant brain samples. All but one rhombencephalitis isolate from cattle were located in clonal complex A. In contrast, food and environmental isolates mainly clustered into clonal complex C, and none was classified as clonal complex A. Isolates of the two main clonal complexes (A and C) obtained by MLVA were analyzed by PCR for the presence of 11 virulence-associated genes (prfA, actA, inlA, inlB, inlC, inlD, inlE, inlF, inlG, inlJ, and inlC2H). Virulence gene analysis revealed significant differences in the actA, inlF, inlG, and inlJ allelic profiles between clinical isolates (complex A) and nonclinical isolates (complex C). The association of particular alleles of actA, inlF, and newly described alleles of inlJ with isolates from CNS infections (particularly rhombencephalitis) suggests that these virulence genes participate in neurovirulence of L. monocytogenes. The overall absence of inlG in clinical complex A and its presence in complex C isolates suggests that the InlG protein is more relevant for the survival of L. monocytogenes in the environment.     

MLVA Based Classification of Mycobacterium tuberculosis Complex Lineages for a Robust Phylogeographic Snapshot of Its Worldwide Molecular Diversity

Multiple-locus variable-number tandem repeat analysis (MLVA) is useful to establish transmission routes and sources of infections for various microorganisms including Mycobacterium tuberculosis complex (MTC). The recently released SITVITWEB database contains 12-loci Mycobacterial Interspersed Repetitive Units – Variable Number of Tandem DNA Repeats (MIRU-VNTR) profiles and spoligotype patterns for thousands of MTC strains; it uses MIRU International Types (MIT) and Spoligotype International Types (SIT) to designate clustered patterns worldwide. Considering existing doubts on the ability of spoligotyping alone to reveal exact phylogenetic relationships between MTC strains, we developed a MLVA based classification for MTC genotypic lineages. We studied 6 different subsets of MTC isolates encompassing 7793 strains worldwide. Minimum spanning trees (MST) were constructed to identify major lineages, and the most common representative located as a central node was taken as the prototype defining different phylogenetic groups. A total of 7 major lineages with their respective prototypes were identified: Indo-Oceanic/MIT57, East Asian and African Indian/MIT17, Euro American/MIT116, West African-I/MIT934, West African-II/MIT664, M. bovis/MIT49, M.canettii/MIT60. Further MST subdivision identified an additional 34 sublineage MIT prototypes. The phylogenetic relationships among the 37 newly defined MIRU-VNTR lineages were inferred using a classification algorithm based on a bayesian approach. This information was used to construct an updated phylogenetic and phylogeographic snapshot of worldwide MTC diversity studied both at the regional, sub-regional, and country level according to the United Nations specifications. We also looked for IS6110 insertional events that are known to modify the results of the spoligotyping in specific circumstances, and showed that a fair portion of convergence leading to the currently observed bias in phylogenetic classification of strains may be traced back to the presence of IS6110. These results shed new light on the evolutionary history of the pathogen in relation to the history of peopling and human migration.     

Frequent occurrence of multidrug‐resistant CC17 Enterococcus faecium among clinical isolates in Sweden

Aims: To screen for the globally spread cluster of Enterococcus faecium, clonal complex 17 (CC17) and characterize the genetic profile of Swedish clinical Ent. faecium isolates. Methods: A total of 203 consecutive isolates collected from 2004 to 2007 from patients with bacteraemia in Sweden. All isolates were genotyped using multiple‐locus variable‐number tandem repeat analysis (MLVA) and 20 isolates representing different MLVA types (MT) were chosen for multilocus sequence typing (MLST). Minimal inhibitory concentrations against clinically relevant antibiotics were determined with agar dilution. Presence of the virulence genes esp and hyl was investigated using PCR. Results: A total of 65% (n = 109) of all isolates belonged to MT‐1, and the second most common MLVA type was MT‐159 (13%, n = 21). MLST analysis confirmed the presence of CC17 during the entire study period. The number of isolates resistant to gentamicin and vancomycin, as well as the presence of hyl, increased significantly during the investigation period. Conclusions: The present study demonstrates that nosocomial infections caused by Ent. faecium CC17 are commonly occurring in Sweden. Significance and Impact of the Study: This is the first report of CC17 Ent. faecium in Sweden. The increase of antibiotic resistance and virulence indicates that these strains are further adapting to the hospital environment.     

The current MLVA typing scheme for <Emphasis Type="Italic">Enterococcus faecium</Emphasis> is less discriminatory than MLST and PFGE for epidemic-virulent,hospital-adapted clonal types

Background  

MLVA (multiple-locus variable-number tandem repeat analysis) is a reliable typing technique introduced recently to differentiate also isolates of Enterococcus faecium. We used the established VNTR (variable number of tandem repeats) scheme to test its suitability to differentiate 58 E. faecium isolates representing mainly outbreaks and clusters of infections and colonizations among patients from 31 German hospitals. All isolates were vancomycin-resistant (vanA type). Typing results for MLVA are compared with results of macrorestriction analysis in PFGE (pulsed-field gel electrophoresis) and MLST (multi-locus sequence typing).     

新疆喀什维吾尔族结核病患者结核分枝杆菌5个VNTRs基因分型研究

目的:应用多位点数目可变串联重复序列分析(multiple loci VNTR analysis,MLVA)技术,对新疆喀什地区维吾尔族结核病患者结核分枝杆菌临床分离株进行基因分型,探讨5个数目可变串联重复序列(VNTR)基因型种类及其分布。方法:收集结核分枝杆菌,采用PCR和琼脂糖凝胶电泳技术,结合BioNumerics5.0软件,对其5个VNTR位点进行结果分析。结果:分离出58株结核分枝杆菌,分为4个基因群21个基因型,分别为Ⅰ群占19.1%,含7个基因型;Ⅱ群占3.4%,含2个基因型;Ⅲ群占67.2%,含9个基因型;Ⅳ群占10.3%,含5个基因型。结论:新疆喀什地区维吾尔族结核病患者的结核分枝杆菌存在明显的基因多态性,且存在主要流行菌群。     

Genetic Evolution of Mycobacterium bovis Causing Tuberculosis in Livestock and Wildlife in France since 1978

To study the dynamics of bovine tuberculosis (bTB) in France, 4,654 M. bovis strains isolated mainly from livestock and wildlife since 1978 were characterized by spoligotyping and MLVA based on MIRU-VNTR. In our study spoligotyping allowed the discrimination of 176 types although 3 spoligotypes are predominant and account for more than half of the total strain population: SB0120 (26%), SB0134 (11%) and SB0121 (6%). In addition, 11% of the isolates, principally from Southern France, showing close spoligotypes and MIRU-VNTR types have been gathered in a family designated as the “F4-family”. MLVA typing allowed extensive discrimination, particularly for strains with predominant spoligotypes, with a total of 498 genotypes, several of which were highly regionalized. The similarity of the strains’ genetic relationships based on spoligotyping and MIRU-VNTR markers supports the co-existence of different clonal populations within the French M. bovis population. A genetic evolution of the strains was observed both geographically and in time. Indeed, as a result of the reduction of bTB due to the national control campaigns, a large reduction of the strains’ genetic variability took place in the last ten years. However, in the regions were bTB is highly prevalent at present, cases in both livestock and in wildlife are due to the spread of unique local genotype profiles. Our results show that the highly discriminating genotyping tools used in this study for molecular studies of bTB are useful for addressing pending questions, which would lead to a better insight into the epidemiology of the disease, and for finding proper solutions for its sustainable control in France.     

Genetic characterization of the Guinea-Bissau family of Mycobacterium tuberculosis complex strains

In a previous study of Mycobacterium tuberculosis complex isolates from Guinea-Bissau in West Africa, we identified a unique group of strains, designated here as the Guinea-Bissau family of strains, which, although genotypically closely related, phenotypically demonstrated a considerable heterogeneity. We conducted here a detailed genotypic analysis of a subset (n = 35) of these isolates. Based on the data obtained, and by comparison of known corresponding genes in mycobacteria outside the M. tuberculosis complex, we propose that the Guinea-Bissau strains belong to a unique branch of the M. tuberculosis complex tree in between classical M. tuberculosis and classical M. bovis. These observations are discussed in their significance in M. tuberculosis complex classification.     

Optimization of variable number tandem repeat typing set for differentiating Mycobacterium tuberculosis strains in the Beijing family

Variable number tandem repeat (VNTR) typing of Mycobacterium tuberculosis, as presently used, has often failed to differentiate clonal strains especially in areas where the Beijing family genotype is predominant. We have evaluated mycobacterial interspersed repetitive units (MIRUs), Queen's University of Belfast (QUB) loci and exact tandem repeat (ETR) loci individually for their abilities to differentiate the Beijing and non-Beijing genotype families of M. tuberculosis. The best locus was QUB 3232 followed by QUB 11b. Available MIRU, QUB and ETR loci were then consecutively arranged in declining order of discriminatory power, and combined to form an optimal set of 12-loci VNTR typing (MIRU 10, 26, 39, 40; QUB 11a, 11b, 15, 18, 26, 3232, 3336; ETR A) that possessed comparable discriminatory ability to IS6110 restriction fragment length polymorphism. This optimized 12-loci VNTR typing set can become an important tool for tracking M. tuberculosis of the Beijing family.     

Molecular characterization of Korean Bacillus anthracis isolates by amplified fragment length polymorphism analysis and multilocus variable-number tandem repeat analysis

We analyzed the genetic relationships and molecular characteristics of 34 Bacillus anthracis isolates from soil and clinical samples in various regions of Korea and 17 related Bacillus species, using the amplified fragment length polymorphism (AFLP) and multilocus variable-number tandem repeat (MLVA) approaches. Triplicate AFLP profiles of these strains showed high reproducibility and identified 376 polymorphisms. AFLP phylogenetic analysis of B. anthracis isolates showed a high level of similarity, 0.93, and this monomorphic fragment profile proved to be useful to differentiate B. anthracis strains from other Bacillus species. The B. cereus group was separated from other Bacillus species at a level of similarity of 0.68. Among them, some B. cereus strains showed genetic interspersion with B. thuringiensis strains. The evolutionary pattern of nucleotide differences among B. anthracis strains with the eight MLVA markers showed nine MLVA types. Three MLVA types, M1 to M3, were pathogenic B. anthracis isolates and were assigned as new genotypes belonging to the A4 and B3 clusters, compared with 89 genotypes deduced from previous data. This indicates that differences in cluster prevalence and distribution may be influenced more by MLVA markers on two plasmids loci and human activity. Consequently, we suggest that the novel MLVA type may represent significant evidence for historic adaptation to environmental conditions of the Asian continent, particularly Korea. Therefore, MLVA techniques may be available for molecular monitoring on anthrax-release-related bioterrorism and further study is required for the continuous epidemiological study of variable anthrax collections.     

Multi-locus variable number tandem repeat analysis for investigation of the genetic association of Clostridium difficile isolates from food, food animals and humans

Clostridium difficile is the primary known cause of antibiotic-associated diarrhea. Diarrheal disease in food animals due to C. difficile infection has been well documented. Recently, reports of C. difficile infections in patients with no known risk factors for disease have raised concern of community acquisition through food animals and food. In this study, multi-locus variable number tandem repeat analysis (MLVA) was performed on a collection of 97C. difficile isolates of human, animal and food origin belonging to either the North American pulsed-field type (NAP) 1 or NAP7/NAP8. MLVA discriminated between NAP1 and NAP7/NAP8 populations. Three clusters of food, food animal and human NAP1 isolates were highly related by MLVA. These data suggest the possibility of either laboratory contamination or widespread distribution of clonal C. difficile populations. Community-associated NAP1 isolates were unrelated to NAP1 food and food animal isolates. Two MLVA loci were absent and 1 was invariant in all NAP7/NAP8 isolates. Therefore, MLVA discrimination was not sufficient to make assessments regarding the genetic associations among food, food animal and human isolates belonging to the NAP7/NAP8 pulsovar. Rigorous epidemiologic and laboratory investigations that employ highly discriminatory genotyping methods are necessary to compare C. difficile isolates from food and food animals to those from humans.