Use of a multilocus variable-number tandem repeat analysis method for molecular subtyping and phylogenetic analysis of Neisseria meningitidis isolates

Demographic Features of Invasive Meningococcal Disease in Taiwan, 1993 to 2020, and Genetic Characteristics of Neisseria meningitidis Isolates, 2003 to 2020

We present the demographic features of invasive meningococcal disease (IMD) in Taiwan between 1993 and 2020 and the genetic characteristics of Neisseria meningitidis isolates recovered from 2003 to 2020. IMD was rare in Taiwan between 1993 and 2020, with an annual incidence ranging from 0.009 to 0.204 per 100,000 people. The case fatality rate (CFR) declined from 18.1% for patients in 1993 to 2002 to 9.8% in 2003 to 2020. Infants less than 12 months were most susceptible to the disease. N. meningitidis serogroup B (NmB) was most predominant, responsible for 81.2% (134/165) of the IMD cases in 2003 to 2020. The majority of the isolates recovered from 2003 to 2020 belonged to 4 worldwide-spread hyperinvasive clonal complexes (cc), cc4821 (30.3%), cc32 (19.4%), cc41/44 (12.7%), cc23 (7.3%), and also a newly assigned clonal complex, cc3439 (10.3%). Core genome multilocus sequence typing (cgMLST) profile comparisons revealed that the cc4821 isolates with a T-to-I substitution at position 91 in gyrA were closely related to those originating from China. Of the 165 isolates, 20.0% and 53.3% were predicted to be covered by the Bexsero and Trumenba vaccines, respectively, whereas, 77.0% and 46.7% remained indeterminate. In conclusion, N. meningitidis isolates recovered in Taiwan between 2003 and 2020 were mostly highly diverse. Most IMD cases appeared sporadically and were caused by localized strains, although some patients were infected by recently introduced strains. cgMLST is a powerful tool for the rapid comparison of genetic relatedness among a large number of isolates. cgMLST profiling, based on 1,241 core genes, and strain tracking can be performed on the website of cgMLST@Taiwan (http://rdvd.cdc.gov.tw/cgMLST/). IMPORTANCE N. meningitidis can cause life-threatening invasive meningococcal disease (IMD), including meningitis and sepsis, resulting in a high CFR and long-term sequelae in survivors. Here, we report the demographic features of IMD in Taiwan over a 28-year period (1993 to 2020) and the genetic characteristics of N. meningitidis isolates recovered from patients with IMD over an 18-year period (2003 to 2020). We conducted a whole-genome sequence analysis to characterize the genetic features of the isolates and developed a cgMLST scheme for epidemiological investigation and strain tracking. The findings can be beneficial in understanding the epidemiology of IMD in Taiwan, the genetic characteristics of the bacterial strains, and the distribution of vaccine antigens for vaccine development and implementation.

Multilocus variable-number tandem-repeat analysis of Neisseria meningitidis serogroup C in China

This study characterized Neisseria meningitidis serogroup C strains in China in order to establish their genetic relatedness and describe the use of multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) to provide useful epidemiological information. A total of 215 N. meningitidis serogroup C strains, obtained from 2003 to 2012 in China, were characterized by MLVA with different published schemes as well as multilocus sequence typing. (i) Based on the MLVA scheme with a combination of five highly variable loci, 203 genotypes were identified; this level of discrimination supports its use for resolving closely related isolates. (ii) Based on a combination of ten low variable loci, clear phylogenetic relationships were established within sequence type complexes. In addition, there was evidence of microevolution of VNTR loci over the decade as strain lineages spread from Anhui to other provinces, the more distant the provinces from Anhui, the higher the genetic variation.

Application of multiple-locus variable number tandem repeat analysis to identify outbreak-associated Neisseria meningitides serogroup C sequence type 4821 in China

Neisseria meningitidis (N. meningitidis) serogroup C sequence type (ST)-4821 caused an outbreak in 2010 in Shandong province of China. Twenty-one non-outbreak-associated strains were isolated, along with twenty-eight N. meningitides serogroup C ST-4821 isolates. Therefore, it’s essential to identify and clarify characterization of the real outbreak-associated strains with a rapid method during an outbreak investigation. In this study, multiple-locus variable number tandem repeat analysis (MLVA) was applied to analyze 84 N. meningitidis strains, among which 58 were recovered from two outbreaks and 26 were sporadic isolates. Three MLVA schemes with different combination of VNTR loci were tested, and two of them were suitable for isolates from China: scheme 2 with six loci was found to separate ST into finer resolution, and scheme 3 with five loci can be used to identify outbreak-associated isolates from the same outbreak that caused by N. meningitidis serogroup C ST-4821.

Genetic diversity of Neisseria meningitidis serogroup C ST-4821 in China based on multiple-locus variable number tandem repeat analysis

Neisseria meningitidis sequence type (ST)-4821 was first reported in China in 2003, and a new hyper-virulent lineage has been designated as the ST-4821 complex. A large number of N. meningitidis ST-4821 strains have been identified in China since 2003; however, the microevolution characteristics of this complex are unclear. Different combinations of variable number of tandem repeats (VNTR) loci were used in multiple-locus VNTR analysis (MLVA) to analyze 118 N. meningitidis serogroup C ST-4821 strains isolated from seventeen provinces between 2003 and 2012. Additionally, MLVA with five VNTR loci was performed due to its high discriminatory power. One hundred and eighteen isolates were found to comprise 112 subtypes based on MLVA, and 16 outbreak-associated strains were clustered into one group. These data indicate a high level of diversity for N. meningitidis ST-4821 due to microevolution in the last decade. In addition, the results revealed high similarity between isolates from the same geographic origins, which is helpful when monitoring the spread of N. meningitidis serogroup C ST-4821 and will provide valuable information for the control and prevention of bacterial meningitis in China.

Multiple-locus variable number tandem repeat analysis of Neisseria gonorrhoeae isolates in Russia

In the present study, a multiple-locus variable number tandem repeat analysis (MLVA) was used to assess the molecular epidemiology of Neisseria gonorrhoeae clinical isolates originating from different regions of Russia. MLVA, based on seven loci, was performed on 218 isolates that were previously tested for susceptibility to penicillin, tetracycline and ciprofloxacin and for the presence of certain genetic determinants of drug resistance. In total, 83 MLVA types were identified, indicating that MLVA is a highly discriminatory technique with a Hunter-Gaston discriminatory index of 0.963 (95% CI, 0.950-0.977). MLVA type 16 was shown to be the most prevalent type and is undoubtedly associated with a multidrug resistant phenotype. The spread of MLVA type 16 from Moscow to Irkutsk suggests that this type has a highly successful transmission rate. Hierarchical cluster analysis of the MLVA profiles classified 208 isolates (95%) into six large groups (containing more than 10 isolates). Clusters differed in geographical characteristics and susceptibility profiles. MLVA cluster A comprised in total 34 isolates and was unambiguously associated with multidrug resistance. Most isolates in cluster A carried mutations in penA, ponA, rpsJ, mtrR, gyrA, and parC genes. MLVA cluster D was associated with resistance to penicillin and with mutations in ponA and rpsJ genes and the presence of plasmid-borne bla(TEM-1) gene. MLVA clusters B, C and E were associated with susceptibility to ciprofloxacin and had a lack of mutations in ponA, rpsJ, gyrA, and parC genes. We conclude that MLVA will be a useful tool for N. gonorrhoeae epidemiological studies.

Global Wolbachia prevalence, titer fluctuations and their potential of causing cytoplasmic incompatibilities in tsetse flies and hybrids of Glossina morsitans subgroup species

We demonstrate the high applicability of a novel VNTR-based (Variable-Number-Tandem-Repeat) molecular screening tool for fingerprinting Wolbachia-infections in tsetse flies. The VNTR-141 locus provides reliable and concise differentiation between Wolbachia strains deriving from Glossina morsitans morsitans, Glossina morsitans centralis, and Glossina brevipalpis. Moreover, we show that certain Wolbachia-infections in Glossina spp. are capable of escaping standard PCR screening methods by 'hiding' as low-titer infections below the detection threshold. By applying a highly sensitive PCR-blot technique to our Glossina specimen, we were able to enhance the symbiont detection limit substantially and, consequently, trace unequivocally Wolbachia-infections at high prevalence in laboratory-reared G. swynnertoni individuals. To our knowledge, Wolbachia-persistence was reported exclusively for field-collected samples, and at low prevalence only. Finally, we highlight the substantially higher Wolbachia titer levels found in hybrid Glossina compared to non-hybrid hosts and the possible impact of these titers on hybrid host fitness that potentially trigger incipient speciation in tsetse flies.

Molecular characterization of Neisseria gonorrhoeae isolates in Almaty, Kazakhstan, by VNTR analysis, Opa-typing and NG-MAST

In the present study, new variable number tandem repeats (VNTR) loci in the Neisseria gonorrhoeae genome were identified in silico. VNTR analysis scheme using PCR and agarose or polyacrylamide gel electrophoresis was developed based on nine VNTR loci with various degrees of polymorphism. The method was used to genotype a collection of 48 isolates, obtained from patients with gonorrhea in Almaty, Kazakhstan during the period from December 2008 to November 2009. This collection of isolates was also characterized by the opa-typing and multiantigen sequence typing (NG-MAST). The discriminatory power of the VNTR analysis translated by Hunter-Gaston Discrimination Index (HGDI) was similar to that of opa typing (HGDI=0.98 versus 0.97) and slightly higher than that of NG-MAST (HDGI=0.95). The adjusted Rand (AR) coefficients and Wallace coefficients showed that the overall concordance between the typing methods was not high. VNTR analysis described here is simple, inexpensive, easy to interpret, and it would be reliable for the comparison of data obtained in different laboratories. The proposed VNTR loci might be used for epidemiological studies of gonococcal infections.

Comparison of PCR ribotyping and multilocus variable-number tandem-repeat analysis (MLVA) for improved detection of Clostridium difficile

Background

Polymerase chain reaction (PCR) ribotyping is one of the globally accepted techniques for defining epidemic clones of Clostridium difficile and tracing virulence-related strains. However, the ambiguous data generated by this technique makes it difficult to compare data attained from different laboratories; therefore, a portable technique that could supersede or supplement PCR ribotyping should be developed. The current study attempted to use a new multilocus variable-number tandem-repeat analysis (MLVA) panel to detect PCR-ribotype groups. In addition, various MLVA panels using different numbers of variable-number tandem-repeat (VNTR) loci were evaluated for their power to discriminate C. difficile clinical isolates.

Results

At first, 40 VNTR loci from the C. difficile genome were used to screen for the most suitable MLVA panel. MLVA and PCR ribotyping were implemented to identify 142 C. difficile isolates. Groupings of serial MLVA panels with different allelic diversity were compared with 47 PCR-ribotype groups. A MLVA panel using ten VNTR loci with limited allelic diversity (0.54-0.83), designated MLVA10, generated groups highly congruent (98%) with the PCR-ribotype groups. For comparison of discriminatory power, a MLVA panel using only four highly variable VNTR loci (allelic diversity: 0.94-0.96), designated MLVA4, was found to be the simplest MLVA panel that retained high discriminatory power. The MLVA10 and MLVA4 were combined and used to detect genetically closely related C. difficile strains.

Conclusions

For the epidemiological investigations of C. difficile, we recommend that MLVA10 be used in coordination with the PCR-ribotype groups to detect epidemic clones, and that the MLVA4 could be used to detect outbreak strains. MLVA10 and MLVA4 could be combined in four multiplex PCR reactions to save time and obtain distinguishable data.

Multilocus variable-number tandem repeat analysis for molecular typing and phylogenetic analysis of Shigella flexneri

Background

Shigella flexneri is one of the causative agents of shigellosis, a major cause of childhood mortality in developing countries. Multilocus variable-number tandem repeat (VNTR) analysis (MLVA) is a prominent subtyping method to resolve closely related bacterial isolates for investigation of disease outbreaks and provide information for establishing phylogenetic patterns among isolates. The present study aimed to develop an MLVA method for S. flexneri and the VNTR loci identified were tested on 242 S. flexneri isolates to evaluate their variability in various serotypes. The isolates were also analyzed by pulsed-field gel electrophoresis (PFGE) to compare the discriminatory power and to evaluate the usefulness of MLVA as a tool for phylogenetic analysis of S. flexneri.

Results

Thirty-six VNTR loci were identified by exploring the repeat sequence loci in genomic sequences of Shigella species and by testing the loci on nine isolates of different subserotypes. The VNTR loci in different serotype groups differed greatly in their variability. The discriminatory power of an MLVA assay based on four most variable VNTR loci was higher, though not significantly, than PFGE for the total isolates, a panel of 2a isolates, which were relatively diverse, and a panel of 4a/Y isolates, which were closely-related. Phylogenetic groupings based on PFGE patterns and MLVA profiles were considerably concordant. The genetic relationships among the isolates were correlated with serotypes. The phylogenetic trees constructed using PFGE patterns and MLVA profiles presented two distinct clusters for the isolates of serotype 3 and one distinct cluster for each of the serotype groups, 1a/1b/NT, 2a/2b/X/NT, 4a/Y, and 6. Isolates that had different serotypes but had closer genetic relatedness than those with the same serotype were observed between serotype Y and subserotype 4a, serotype X and subserotype 2b, subserotype 1a and 1b, and subserotype 3a and 3b.

Conclusions

The 36 VNTR loci identified exhibited considerably different degrees of variability among S. flexneri serotype groups. VNTR locus could be highly variable in a serotype but invariable in others. MLVA assay based on four highly variable loci could display a comparable resolving power to PFGE in discriminating isolates. MLVA is also a prominent molecular tool for phylogenetic analysis of S. flexneri; the resulting data are beneficial to establish clear clonal patterns among different serotype groups and to discern clonal groups among isolates within the same serotype. As highly variable VNTR loci could be serotype-specific, a common MLVA protocol that consists of only a small set of loci, for example four to eight loci, and that provides high resolving power to all S. flexneri serotypes may not be obtainable.

Genetic diversity of O157:H7 and non-O157 verocytotoxigenic Escherichia coli from Argentina inferred from multiple-locus variable-number tandem repeat analysis (MLVA)

Although serotype O157:H7 has been implicated in most cases of haemolytic-uraemic syndrome (HUS), there is growing concern about non-O157 serotypes of verocytotoxigenic Escherichia coli (VTEC). Multiple-locus variable-number tandem repeat analysis (MLVA) has been focused on the specific typing of O157:H7 isolates, but recently, a generic MLVA assay for E. coli and Shigella has been developed. We performed a study of the polymorphism in 7 generic VNTR loci both in VTEC O157:H7 and non-O157 isolates from Argentina, in order to asses the ability of the method to type this group of isolates and to get insight into their genetic diversity. Sixty-four isolates from cattle, patients with diarrhoea, and contaminated food belonging to 8 different serotypes were studied. All of them could be typed by this method and revealed 41 different MLVA genotypes. The MLVA dendrogram showed 2 main clusters which corresponded to O157:H7 and non-O157, respectively. Our results confirm the suitability of this MLVA method for analyzing VTEC isolates belonging to several serotypes, both O157:H7 as well as non-O157, highlight the genetic variability of the O157:H7 serotype and the need of additional research in order to find more VNTR loci that could allow a higher discrimination among non-O157 VTEC.

Utility of multilocus variable-number tandem-repeat analysis as a molecular tool for phylogenetic analysis of Shigella sonnei

A panel of 916 isolates, including 703 closely related IST1 isolates, were characterized by inter-IS1 spacer typing (IST), pulsed-field gel electrophoresis (PFGE), and multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) to evaluate the utility of MLVA as a molecular tool for the phylogenetic analysis of Shigella sonnei. The global phylogenetic patterns determined by IST, PFGE, and MLVA were concordant. MLVA was carried out using 26 VNTR loci with a range of degrees of variability. MLVA data for the 703 IST1 isolates revealed that diversification among the closely related isolates was attributed mainly to four highly variable loci. The phylogenetic pattern for the closely related isolates determined using MLVA profiles of 8 highly variable loci was in agreement with that determined using the 26-locus profiles. A clustering analysis using the profiles of 18 loci with limited variability established clear phylogenetic relationships among IST clonal groups. Accordingly, MLVA is a useful tool for the phylogenetic analysis of S. sonnei. Combined VNTR loci with higher variability are useful markers for resolving closely related isolates, whereas combined loci with lower variability are suitable for establishing clear phylogenetic relationships between strains or clones that have evolved over a longer timescale.

Multiple locus variable number of tandem repeats analysis

Genotyping of bacteria through typing of loci containing a variable number of tandem repeats (VNTR) might become the gold standard for many pathogens. The development of genome sequencing has shown that such sequences were present in every species analyzed, and that polymorphism exists in at least a fraction of them. The length of these repetitions can vary from a single nucleotide to a few hundreds. This has implications for both the techniques used to measure the repeat number and the level of variability. In addition, tandem repeats can be part of coding regions or be intergenic and may play a direct role in the adaptation to the environment, thus having different observed evolution rates. For these reasons the choice of VNTR when setting a multiple-loci VNTR analysis (MLVA) assay is important. Although reasonable discrimination can be achieved with the typing of six to eight markers, in particular in species with high genomic diversity, it may be necessary to type 20 to 40 markers in monomorphic species or if an evolutionary meaningful assay is needed. Homoplasy (in the present context, two alleles containing the same repeat copy number in spite of a different history) is then compensated by the analysis of multiple markers. Finally, even if the underlying principles are relatively simple, quality standards must be implemented before this approach is widely accepted, and technology issues must be resolved to further lower the typing costs.

Variable-number-of-tandem-repeats analysis of genetic diversity in Pasteuria ramosa

Variable-number-of-tandem-repeats (VNTR) markers are increasingly being used in population genetic studies of bacteria. They were recently developed for Pasteuria ramosa, an endobacterium that infects Daphnia species. In the present study, we genotyped P. ramosa in 18 infected hosts from the United Kingdom, Belgium, and two lakes in the United States using seven VNTR markers. Two Daphnia species were collected: D. magna and D. dentifera. Six loci showed length polymorphism, with as many as five alleles identified for a single locus. Similarity coefficient calculations showed that the extent of genetic variation between pairs of isolates within populations differed according to the population, but it was always less than the genetic distances among populations. Analysis of the genetic distances performed using principal component analysis revealed strong clustering by location of origin, but not by host Daphnia species. Our study demonstrated that the VNTR markers available for P. ramosa are informative in revealing genetic differences within and among populations and may therefore become an important tool for providing detailed analysis of population genetics and epidemiology.

Multilocus variable-number tandem-repeat analysis for molecular typing of Shigella sonnei

A multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) method was developed and evaluated for the subtyping of Shigella sonnei isolates. A total of 26 VNTR loci were identified by exploring the repeat sequence loci in the genomic sequences of S. sonnei strains Ss046 and 53G and by testing 536 isolates that had previously been characterized by pulsed-field gel electrophoresis (PFGE). The discriminatory power of MLVA (Simpson's index of diversity [D], 0.9524; 95% confidence interval [CI], 0.9373 to 0.9564) for the 536 isolates was significantly higher than that of PFGE (D, 0.8882; CI, 0.8667 to 0.9097). MLVA typing with the four and eight most variable loci had D values of 0.9468 and 0.9481, respectively, results approaching that of 26 loci. The usefulness of MLVA for outbreak investigation was evaluated using 151 isolates from 10 shigellosis outbreaks and 22 PFGE-indistinguishable isolates collected from nine epidemiologically unrelated events in five different countries. The evaluations indicated that MLVA was a powerful typing tool to distinguish isolates for outbreak investigation and that it exhibited a good discrimination of the 22 PFGE-indistinguishable isolates. Single-locus variants did occur during the outbreak; therefore, S. sonnei isolates with MLVA profiles differing at no more than a single locus should be considered part of the same outbreak. The present study suggests that MLVA has the potential to replace PFGE as a standard method of typing S. sonnei isolates for disease surveillance and outbreak investigation.

The current MLVA typing scheme for Enterococcus faecium 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).

Results

All 51 but one hospital isolates from 1996–2006 were assigned to the clonal complex (CC) of epidemic-virulent, hospital-adapted lineages (MLST CC-17; MLVA CC-1) and differed from isolates of sporadic infections and colonizations (n = 7; 1991–1995) and other non-hospital origins (n = 27). Typing of all 58 hospital VRE revealed MLVA as the least discriminatory method (Simpson's diversity index 0.847) when compared to MLST (0.911) and PFGE (0.976). The two most common MLVA types MT-1 (n = 16) and MT-159 (n = 14) combined isolates of several MLST types including also major epidemic, hospital-adapted, clonal types (MT-1: ST-17, ST-18, ST-280, ST-282; MT-159: ST-78, ST-192, ST-203). These data clearly indicate that non-related E. faecium could possess an identical MLVA type being especially critical when MLVA is used to elucidate supposed outbreaks with E. faecium within a single or among different hospitals. Stability of a given MLVA profile MT-12 (ST-117) during an outbreak over a period of five years was also shown.

Conclusion

MLVA is a suitable method to assign isolates of E. faecium into distinct clonal complexes. To investigate outbreaks the current MLVA typing scheme for E. faecium does not discriminate enough and cannot be recommended as a standard superior to PFGE.

VNTRDB: a bacterial variable number tandem repeat locus database

Variable number tandem repeat-PCR (VNTR-PCR) is a novel method developed for molecular typing of microorganisms. This method has proven useful in epidemiological studies in medical microbiology. Although hundreds of bacterial genomes have been sequenced, variable number tandem repeats (TRs) derived from comparative genome analyses are scarce. This may hamper their application to the surveillance of bacteria in molecular epidemiology. Here, we present a freely accessible variable number tandem repeat database (VNTRDB) that is intended to be a resource for helping in the discovery of putatively polymorphic tandem repeat loci and to aid with assay design by providing the flanking sequences that can be used in subsequent PCR primer design. In order to reveal possible polymorphism, each TR locus was obtained by comparing the sequences between different sets of bacterial genera, species or strains. Through this comparison, TRs which are unique to a genus can also be identified. Moreover, a visualization tool is provided to ensure that the copy number and locus length of repeats are correct. The VNTRDB is available at .