New Multilocus Variable-Number Tandem-Repeat Analysis Tool for Surveillance and Local Epidemiology of Bacterial Leaf Blight and Bacterial Leaf Streak of Rice Caused by Xanthomonas oryzae

Multilocus variable-number tandem-repeat analysis (MLVA) is efficient for routine typing and for investigating the genetic structures of natural microbial populations. Two distinct pathovars of Xanthomonas oryzae can cause significant crop losses in tropical and temperate rice-growing countries. Bacterial leaf streak is caused by X. oryzae pv. oryzicola, and bacterial leaf blight is caused by X. oryzae pv. oryzae. For the latter, two genetic lineages have been described in the literature. We developed a universal MLVA typing tool both for the identification of the three X. oryzae genetic lineages and for epidemiological analyses. Sixteen candidate variable-number tandem-repeat (VNTR) loci were selected according to their presence and polymorphism in 10 draft or complete genome sequences of the three X. oryzae lineages and by VNTR sequencing of a subset of loci of interest in 20 strains per lineage. The MLVA-16 scheme was then applied to 338 strains of X. oryzae representing different pathovars and geographical locations. Linkage disequilibrium between MLVA loci was calculated by index association on different scales, and the 16 loci showed linear Mantel correlation with MLSA data on 56 X. oryzae strains, suggesting that they provide a good phylogenetic signal. Furthermore, analyses of sets of strains for different lineages indicated the possibility of using the scheme for deeper epidemiological investigation on small spatial scales.     

Features of Variable Number of Tandem Repeats in Yersinia pestis and the Development of a Hierarchical Genotyping Scheme

Background

Variable number of tandem repeats (VNTRs) that are widely distributed in the genome of Yersinia pestis proved to be useful markers for the genotyping and source-tracing of this notorious pathogen. In this study, we probed into the features of VNTRs in the Y. pestis genome and developed a simple hierarchical genotyping system based on optimized VNTR loci.

Methodology/Principal Findings

Capillary electrophoresis was used in this study for multi-locus VNTR analysis (MLVA) in 956 Y. pestis strains. The general features and genetic diversities of 88 VNTR loci in Y. pestis were analyzed with BioNumerics, and a “14+12” loci-based hierarchical genotyping system, which is compatible with single nucleotide polymorphism-based phylogenic analysis, was established.

Conclusions/Significance

Appropriate selection of target loci reduces the impact of homoplasies caused by the rapid mutation rates of VNTR loci. The optimized “14+12” loci are highly discriminative in genotyping and source-tracing Y. pestis for molecular epidemiological or microbial forensic investigations with less time and lower cost. An MLVA genotyping datasets of representative strains will improve future research on the source-tracing and microevolution of Y. pestis.     

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.     

A new typing technique for the Rickettsiales Ehrlichia ruminantium: multiple-locus variable number tandem repeat analysis

Ehrlichia ruminantium (ER) is a member of the order Rickettsiales transmitted by Amblyomma ticks. This obligatory intracellular bacterium is the causative agent of a fatal disease in ruminants, named heartwater. It represents a constraint on breeding development in sub-Saharan Africa and in the Caribbean. The genetic diversity of the strains of ER, which could be a limiting factor to obtain effective vaccines, needs to be better characterized. For this purpose, we developed a molecular typing technique based on the polymorphism of variable number tandem repeat (VNTR) sequences, MLVA (multiple locus VNTR analysis).Eight (out of 21) VNTR candidates were validated using 17 samples representing a panel of ER strains from different geographical origins from West, South Africa, and Caribbean areas and in ER infected ticks and goat tissues. This result demonstrated the ability of these VNTRs to type a wide range of strains. The stability of the selected VNTR markers was very good, at the time scale needed for epidemiological purposes: in particular, no difference in the VNTR profiles was observed between virulent and attenuated strains (for Gardel and Senegal strains) and between strains (Gardel and Blonde strains) isolated in the same area 19 years apart. We validated the strong discriminatory power of MLVA for ER and found a high level of polymorphism between the available strains, with 10 different profiles out of 13 ER strains.The MLVA scheme described in this study is a rapid and efficient molecular typing tool for ER, which allows rapid and direct typing of this intracellular pathogen without preliminary culture and gives reliable results that can be used for further epidemiological studies.     

Differentiation of Strains of Xylella fastidiosa by a Variable Number of Tandem Repeat Analysis

Short sequence repeats (SSRs) with a potential variable number of tandem repeat (VNTR) loci were identified in the genome of the citrus pathogen Xylella fastidiosa and used for typing studies. Although mono- and dinucleotide repeats were absent, we found several intermediate-length 7-, 8-, and 9-nucleotide repeats, which we examined for allelic polymorphisms using PCR. Five genuine VNTR loci were highly polymorphic within a set of 27 X. fastidiosa strains from different hosts. The highest average Nei's measure of genetic diversity (H) estimated for VNTR loci was 0.51, compared to 0.17 derived from randomly amplified polymorphic DNA (RAPD) analysis. For citrus X. fastidiosa strains, some specific VNTR loci had a H value of 0.83, while the maximum value given by specific RAPD loci was 0.12. Our approach using VNTR markers provides a high-resolution tool for epidemiological, genetic, and ecological analysis of citrus-specific X. fastidiosa strains.     

Highly diverse variable number tandem repeat loci in the E. coli O157:H7 and O55:H7 genomes for high-resolution molecular typing

AIM: Evaluation of the Escherichia coli genome for variable number tandem repeat (VNTR) loci in order to provide a subtyping tool with greater discrimination and more efficient capacity. METHODS AND RESULTS: Twenty-nine putative VNTR loci were identified from the E. coli genomic sequence. Their variability was validated by characterizing the number of repeats at each locus in a set of 56 E. coli O157:H7/HN and O55:H7 isolates. An optimized multiplex assay system was developed to facility high capacity analysis. Locus diversity values ranged from 0.23 to 0.95 while the number of alleles ranged from two to 29. This multiple-locus VNTR analysis (MLVA) data was used to describe genetic relationships among these isolates and was compared with PFGE (pulse field gel electrophoresis) data from a subset of the same strains. Genetic similarity values were highly correlated between the two approaches, through MLVA was capable of discrimination amongst closely related isolates when PFGE similar values were equal to 1.0. CONCLUSIONS: Highly variable VNTR loci exist in the E. coli O157:H7 genome and are excellent estimators of genetic relationships, in particular for closely related isolates. SIGNIFICANCE AND IMPACT OF THE STUDY: Escherichia coli O157:H7 MLVA offers a complimentary analysis to the more traditional PFGE approach. Application of MLVA to an outbreak cluster could generate superior molecular epidemiology and result in a more effective public health response.     

Subtyping of a Large Collection of Historical Listeria monocytogenes Strains from Ontario,Canada, by an Improved Multilocus Variable-Number Tandem-Repeat Analysis (MLVA)

Listeria monocytogenes is responsible for severe and often fatal food-borne infections in humans. A collection of 2,421 L. monocytogenes isolates originating from Ontario''s food chain between 1993 and 2010, along with Ontario clinical isolates collected from 2004 to 2010, was characterized using an improved multilocus variable-number tandem-repeat analysis (MLVA). The MLVA method was established based on eight primer pairs targeting seven variable-number tandem-repeat (VNTR) loci in two 4-plex fluorescent PCRs. Diversity indices and amplification rates of the individual VNTR loci ranged from 0.38 to 0.92 and from 0.64 to 0.99, respectively. MLVA types and pulsed-field gel electrophoresis (PFGE) patterns were compared using Comparative Partitions analysis involving 336 clinical and 99 food and environmental isolates. The analysis yielded Simpson''s diversity index values of 0.998 and 0.992 for MLVA and PFGE, respectively, and adjusted Wallace coefficients of 0.318 when MLVA was used as a primary subtyping method and 0.088 when PFGE was a primary typing method. Statistical data analysis using BioNumerics allowed for identification of at least 8 predominant and persistent L. monocytogenes MLVA types in Ontario''s food chain. The MLVA method correctly clustered epidemiologically related outbreak strains and separated unrelated strains in a subset analysis. An MLVA database was established for the 2,421 L. monocytogenes isolates, which allows for comparison of data among historical and new isolates of different sources. The subtyping method coupled with the MLVA database will help in effective monitoring/prevention approaches to identify environmental contamination by pathogenic strains of L. monocytogenes and investigation of outbreaks.     

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.     

Development of New Multilocus Variable Number of Tandem Repeat Analysis (MLVA) for Listeria innocua and Its Application in a Food Processing Plant

Listeria innocua is an important hygiene indicator bacterium in food industries because it behaves similar to Listeria monocytogenes, which is pathogenic to humans. PFGE is often used to characterize bacterial strains and to track contamination source. However, because PFGE is an expensive, complicated, time-consuming protocol, and poses difficulty in data sharing, development of a new typing method is necessary. MLVA is a technique that identifies bacterial strains on the basis of the number of tandem repeats present in the genome varies depending on the strains. MLVA has gained attention due to its high reproducibility and ease of data sharing. In this study, we developed a MLVA protocol to assess L. innocua and evaluated it by tracking the contamination source of L. innocua in an actual food manufacturing factory by typing the bacterial strains isolated from the factory. Three VNTR regions of the L. innocua genome were chosen for use in the MLVA. The number of repeat units in each VNTR region was calculated based on the results of PCR product analysis using capillary electrophoresis (CE). The calculated number of repetitions was compared with the results of the gene sequence analysis to demonstrate the accuracy of the CE repeat number analysis. The developed technique was evaluated using 60 L. innocua strains isolated from a food factory. These 60 strains were classified into 11 patterns using MLVA. Many of the strains were classified into ST-6, revealing that this MLVA strain type can contaminate each manufacturing process in the factory. The MLVA protocol developed in this study for L. innocua allowed rapid and easy analysis through the use of CE. This technique was found to be very useful in hygiene control in factories because it allowed us to track contamination sources and provided information regarding whether the bacteria were present in the factories.     

Origin of the Outbreak in France of Pseudomonas syringae pv. actinidiae Biovar 3, the Causal Agent of Bacterial Canker of Kiwifruit,Revealed by a Multilocus Variable-Number Tandem-Repeat Analysis

The first outbreaks of bacterial canker of kiwifruit caused by Pseudomonas syringae pv. actinidiae biovar 3 were detected in France in 2010. P. syringae pv. actinidiae causes leaf spots, dieback, and canker that sometimes lead to the death of the vine. P. syringae pv. actinidifoliorum, which is pathogenic on kiwi as well, causes only leaf spots. In order to conduct an epidemiological study to track the spread of the epidemics of these two pathogens in France, we developed a multilocus variable-number tandem-repeat (VNTR) analysis (MLVA). MLVA was conducted on 340 strains of P. syringae pv. actinidiae biovar 3 isolated in Chile, China, France, Italy, and New Zealand and on 39 strains of P. syringae pv. actinidifoliorum isolated in Australia, France, and New Zealand. Eleven polymorphic VNTR loci were identified in the genomes of P. syringae pv. actinidiae biovar 3 ICMP 18744 and of P. syringae pv. actinidifoliorum ICMP 18807. MLVA enabled the structuring of P. syringae pv. actinidiae biovar 3 and P. syringae pv. actinidifoliorum strains in 55 and 16 haplotypes, respectively. MLVA and discriminant analysis of principal components revealed that strains isolated in Chile, China, and New Zealand are genetically distinct from P. syringae pv. actinidiae strains isolated in France and in Italy, which appear to be closely related at the genetic level. In contrast, no structuring was observed for P. syringae pv. actinidifoliorum. We developed an MLVA scheme to explore the diversity within P. syringae pv. actinidiae biovar 3 and to trace the dispersal routes of epidemic P. syringae pv. actinidiae biovar 3 in Europe. We suggest using this MLVA scheme to trace the dispersal routes of P. syringae pv. actinidiae at a global level.     

Study of polymorphic variable-number of tandem repeats loci in the ECOR collection and in a set of pathogenic Escherichia coli and Shigella isolates for use in a genotyping assay

The Escherichia coli (E. coli) reference collection, ECOR, consists of 72 strains that are representative of the genotypic diversity, as indexed by multilocus enzyme electrophoresis (MLEE), in the species as a whole. MLEE revealed 4 main phylogenetic groups designated A, B1, B2 and D. We present a study of the relationship between the ECOR strains as determined by polymorphisms in seven variable-number of tandem repeats (VNTR) loci. Seven tandem repeats that were present in more than one of the fully sequenced E. coli strains were selected, and primers were constructed in order to amplify the targets in all species where the loci were present. The combined result for all VNTR loci was adapted as a multiple-locus variable-number tandem repeats analysis (MLVA) and showed that the ECOR collection was divided into 63 distinct genotypes. The ECOR phylogenetic groups defined by MLEE were not well conserved by MLVA. A set of 61 pathogenic isolates of both E. coli and Shigella spp. was then tested with the same set of VNTR loci, and revealed 54 distinct genotypes. In addition, the MLVA method was used to genotype isolates from patients and suspected sources in a recent outbreak of E. coli O103 in Norway. The pathogenic E. coli isolates contained the diarrhea causing categories EIEC, EAEC, STEC, ETEC and EPEC. Shigella isolates were of species S. flexneri, S. boydii, S. sonnei and S. dysenteriae. The MLVA method rapidly genotyped all isolates in the study at a Simpson's index of diversity of D=0.98.     

Rapid and high resolution genotyping of all Escherichia coli serotypes using 10 genomic repeat-containing loci

Our laboratory has previously published two multiple-locus variable-number tandem-repeats analysis (MLVA) methods for rapid genotyping of Escherichia coli (E. coli), which are now in routine use for surveillance and outbreak detection. The first assay developed was specific for E. coli O157:H7; however this assay was not suitable for genotyping other E. coli serotypes. A new generic MLVA-assay was then developed with the capability of genotyping all E. coli serotypes. This generic E. coli MLVA (GECM7) was based on polymorphism in seven variable number of tandem repeats (VNTR) loci. GECM7 worked well with the majority of E. coli serotypes; however we wanted to increase the resolution for this method based in part of comparison with PFGE typing of E. coli O26:H11, where PFGE appeared to display higher resolution. The GECM7 method was improved by adding three new repeat-loci to a total of ten (GECM10), and a considerable increase in resolution was observed (from 296 to 507 genotypes on the same set of strains).     

The use of MLVA for Corynebacterium diphtheriae genotyping--preliminary studies

The complete genome sequence of strain NCTC 13129 C. diphtheriae were investigated in order to identify tandem repeats (VNTR). From 75 VNTR loci identified in the genome 14 were selected. Primers were designed and PCR conditions were optimized for amplification of the selected VNTR markers. Preliminary studies of usefulness of selected VNTR markers were conducted using a group of 28 C. diphtheriae strains. From 14 markers 8 were regarded as potentially useful. The diversity of individual markers ranged from 1 to 6 alleles (Simpson index from 0 to 0,746). No diversity were observed for 3 VNTR markers but it could be a results of too small group of strains analyzed in the tests. Simpson diversity index calculated for all the markers tested on 28 strains was 0,87. Results of the preliminary studies showed usefulness of MLVA for C. diphtheriae genotyping. Nevertheless, confirmation of reliability of the method should be done using a large group of strains. Moreover, the method should be compared with other genotyping methods.     

Development of Multiple-Locus Variable-Number Tandem-Repeat Analysis for Molecular Subtyping of Campylobacter jejuni by Using Capillary Electrophoresis

Campylobacter jejuni is a common cause of the frequently reported food-borne diseases in developed and developing nations. This study describes the development of multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA) using capillary electrophoresis as a novel typing method for microbial source tracking and epidemiological investigation of C. jejuni. Among 36 tandem repeat loci detected by the Tandem Repeat Finder program, 7 VNTR loci were selected and used for characterizing 60 isolates recovered from chicken meat samples from retail shops, samples from chicken meat processing factory, and stool samples. The discrimination ability of MLVA was compared with that of multilocus sequence typing (MLST). MLVA (diversity index of 0.97 with 31 MLVA types) provided slightly higher discrimination than MLST (diversity index of 0.95 with 25 MLST types). The overall concordance between MLVA and MLST was estimated at 63% by adjusted Rand coefficient. MLVA predicted MLST type better than MLST predicted MLVA type, as reflected by Wallace coefficient (Wallace coefficient for MLVA to MLST versus MLST to MLVA, 86% versus 51%). MLVA is a useful tool and can be used for effective monitoring of C. jejuni and investigation of epidemics caused by C. jejuni.     

Genotyping and Phylogenetic Analysis of Yersinia pestis by MLVA: Insights into the Worldwide Expansion of Central Asia Plague Foci

Background

The species Yersinia pestis is commonly divided into three classical biovars, Antiqua, Medievalis, and Orientalis, belonging to subspecies pestis pathogenic for human and the (atypical) non-human pathogenic biovar Microtus (alias Pestoides) including several non-pestis subspecies. Recent progress in molecular typing methods enables large-scale investigations in the population structure of this species. It is now possible to test hypotheses about its evolution which were proposed decades ago. For instance the three classical biovars of different geographical distributions were suggested to originate from Central Asia. Most investigations so far have focused on the typical pestis subspecies representatives found outside of China, whereas the understanding of the emergence of this human pathogen requires the investigation of strains belonging to subspecies pestis from China and to the Microtus biovar.

Methodology/Principal Findings

Multi-locus VNTR analysis (MLVA) with 25 loci was performed on a collection of Y. pestis isolates originating from the majority of the known foci worldwide and including typical rhamnose-negative subspecies pestis as well as rhamnose-positive subspecies pestis and biovar Microtus. More than 500 isolates from China, the Former Soviet Union (FSU), Mongolia and a number of other foci around the world were characterized and resolved into 350 different genotypes. The data revealed very close relationships existing between some isolates from widely separated foci as well as very high diversity which can conversely be observed between nearby foci.

Conclusions/Significance

The results obtained are in full agreement with the view that the Y. pestis subsp. pestis pathogenic for humans emerged in the Central Asia region between China, Kazakhstan, Russia and Mongolia, only three clones of which spread out of Central Asia. The relationships among the strains in China, Central Asia and the rest of the world based on the MLVA25 assay provide an unprecedented view on the expansion and microevolution of Y. pestis.     

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.     

Multiple-locus variable-number of tandem-repeats analysis distinguishes Vibrio parahaemolyticus pandemic O3:K6 strains

A specific serotype of Vibrio parahaemolyticus, O3:K6, has recently been linked to epidemics of gastroenteritis in Southeast Asia, Japan, and North America. These pandemic O3:K6 strains appear to have recently spread across continents from a single origin to reach global coverage, based on profiling of strains by several molecular typing methods. In this study, variable-number tandem repeats (VNTR)-based fingerprinting was applied to clinical and environmental V. parahaemolyticus O3:K6 strains in an attempt to develop a molecular method with increased sensitivity for discriminating strains; the relative discriminatory powers were compared with ribotyping and pulsed-field gel electrophoresis (PFGE). All clinical strains tested were independent human isolates obtained from different outbreaks or from sporadic cases in Tokyo during the period from 1996 to 2003. Multiple-locus VNTR analysis (MLVA) was shown to have high resolution and reproducibility for typing of V. parahaemolyticus clones. MLVA analysis of 28 pandemic V. parahaemolyticus O3:K6 strains isolated from human cases produced 28 distinct VNTR patterns. The VNTR loci displayed between 2 and 15 alleles at each of eight loci with Nei's diversity index ranging from 0.35 and 0.91. These data demonstrated that MLVA is useful for individual strain typing of new O3:K6 strains, which appear to be closely related by other molecular methods.     

Use of the VNTR typing technique to determine the origin of Mycobacterium tuberculosis strains isolated from Filipino patients in Korea

With increasing international interchange of personnel, international monitoring is necessary to decrease tuberculosis incidence in the world. This study aims to develop a new tool to determine origin of Mycobacterium tuberculosis strains isolated from Filipino patients living in Korea. Thirty-two variable number tandem repeat (VNTR) loci were used for discrimination of 50 Filipino M. tuberculosis strains isolated in the Philippines, 317 Korean strains isolated in Korea, and 8 Filipino strains isolated in Korea. We found that the VNTR loci 0580, 0960, 2531, 2687, 2996, 0802, 2461, 2163a, 4052, 0424, 1955, 2074, 2347, 2401, 3171, 3690, 2372, 3232, and 4156 had different mode among copy numbers or exclusively distinct copy number in VNTR typing between Filipino and Korean M. tuberculosis strains. When these differences of the VNTR loci were applied to 8 Filipino M. tuberculosis strains isolated in Korea, 6 of them revealed Filipino type while 2 of them had Korean type. Using the differences of mode or repeated number of VNTR loci were very useful in distinguishing the Filipino strain from Korean strain.     

Genetic diversity of Francisella tularensis subspecies holarctica strains isolated in Japan

The recently developed MLVA has high discriminatory power for the typing of individual strains or isolates of Francisella tularensis. In the present study, MLVA was applied to 33 Japanese F. tularensis subspecies holarctica strains to examine the genetic diversity of F. tularensis isolated. Among the seven VNTR loci analyzed, Ft-M2, Ft-M10, and Ft-M20 loci showed high genetic polymorphism in Japanese strains, whereas Ft-M3 was most variable in non-Japanese strains. These results provide novel extended information about the genomic diversity among the strains of F. tularensis ssp. holarctica distributed in Japan and enable determination of whether a given isolate is indigenous to Japan by examining these loci using MLVA.     

Tandem repeat regions within the Burkholderia pseudomallei genome and their application for high resolution genotyping

Background

The facultative, intracellular bacterium Burkholderia pseudomallei is the causative agent of melioidosis, a serious infectious disease of humans and animals. We identified and categorized tandem repeat arrays and their distribution throughout the genome of B. pseudomallei strain K96243 in order to develop a genetic typing method for B. pseudomallei. We then screened 104 of the potentially polymorphic loci across a diverse panel of 31 isolates including B. pseudomallei, B. mallei and B. thailandensis in order to identify loci with varying degrees of polymorphism. A subset of these tandem repeat arrays were subsequently developed into a multiple-locus VNTR analysis to examine 66 B. pseudomallei and 21 B. mallei isolates from around the world, as well as 95 lineages from a serial transfer experiment encompassing ~18,000 generations.

Results

B. pseudomallei contains a preponderance of tandem repeat loci throughout its genome, many of which are duplicated elsewhere in the genome. The majority of these loci are composed of repeat motif lengths of 6 to 9 bp with 4 to 10 repeat units and are predominately located in intergenic regions of the genome. Across geographically diverse B. pseudomallei and B.mallei isolates, the 32 VNTR loci displayed between 7 and 28 alleles, with Nei's diversity values ranging from 0.47 and 0.94. Mutation rates for these loci are comparable (>10^-5 per locus per generation) to that of the most diverse tandemly repeated regions found in other less diverse bacteria.

Conclusion

The frequency, location and duplicate nature of tandemly repeated regions within the B. pseudomallei genome indicate that these tandem repeat regions may play a role in generating and maintaining adaptive genomic variation. Multiple-locus VNTR analysis revealed extensive diversity within the global isolate set containing B. pseudomallei and B. mallei, and it detected genotypic differences within clonal lineages of both species that were identical using previous typing methods. Given the health threat to humans and livestock and the potential for B. pseudomallei to be released intentionally, MLVA could prove to be an important tool for fine-scale epidemiological or forensic tracking of this increasingly important environmental pathogen.