Multiple-Locus Variable-Number Tandem Repeat Analysis of Vibrio cholerae in Comparison with Pulsed Field Gel Electrophoresis and Virulotyping

Molecular analysis of Malaysian Vibrio cholerae was carried out using a multiple-locus variable-number tandem repeat analysis (MLVA) assay based on 7 loci of V. cholerae. The discriminatory ability of the assay was compared with pulsed-field gel electrophoresis (PFGE) using 43 Malaysian V. cholerae isolated from various sources. In addition, the virulotypes of the strains were determined. Based on MLVA, 38 allelic profiles were obtained (F = 0.63) while PFGE generated 35 pulsotypes (F = 0.71). Simpson's index of diversity for different VNTR loci ranged from 0.59 to 0.92. The combined loci increased the discriminatory index to 0.99 which was comparable with PFGE (D = 0.99). Most of the environmental non-O1/non-O139 strains harbored rtxA, rstR, toxR, and hlyA only, and the virulotype of this serogroup was significantly different (P < .01) from clinical/environmental O1 and environmental O139 strains. In conclusion, the MLVA assay developed in this study was a useful genotyping tool with comparable discriminatory power with PFGE. In addition, the combination of the two approaches can further distinguish the strains from different sources and geographical regions of isolation.


Introduction
Vibrio cholerae is a gram negative bacterium which lives freely in aquatic environment and causes cholera [1]. Cholera is endemic in many parts of the world, especially the countries which lack proper sanitation managements.
In Malaysia, cholera outbreaks due to the V. cholerae O1 serotype which occurs periodically [2]. The ratio of distribution of V. cholerae O139 to O1 serogroups isolated from seafood from 1998 to 1999 was 14 : 1. Non-O1/non-O139 V. cholerae is also frequently isolated from seafood and water sources but has not been implicated in any major outbreaks [3,4]. Although non-O1/non-O139 V. cholerae is not associated with any major outbreak, it has been reported to be responsible for sporadic cases of diarrhea [5][6][7].
The well-known genes associated with colonization are ctxA and tcpA. These genes are commonly found in O1 and O139 serogroups. However, the mechanism of colonization for the non-O1/non-O139 strains remains relatively unknown as they do not harbor ctxA and tcpA [5]. Olivier et al. [8] had reported that accessory toxins such as hemolysin and multifunctional autoprocessing RTX toxin in El Tor V. cholerae are involved in prolonged colonization without cholera toxin (CT) or toxin-coregulated pili (TCP). As these accessory virulence genes are commonly found in all serogroups of V. cholerae, it is of interest to investigate the involvement of these accessory virulence genes for prolonged colonization in other serogroups of V. cholerae.
Molecular subtyping of pathogen is important for tracing a new or previously found virulent or multidrug-resistant clone. Genomic variation and epidemiological study for different serogroups of V. cholerae have been carried out using many DNA-fingerprinting tools. PFGE is the most common subtyping tool to define strains from outbreaks and from sporadic cases of cholera as it has the highest discriminatory ability [4]. However, PFGE is time consuming and requires strict adherence to standardize protocols for interlaboratory comparison.
Multiple-locus variable-number of tandem repeat analysis (MLVA), a method based on the tandem repeats in multiple loci, was developed and has been popularly adopted for differentiation of bacterial pathogens since 2006 [9]. Recently, Olsen et al. [10] evaluated the usefulness of MLVA on V. cholerae isolated worldwide based on polymorphism in 6 VNTR loci and demonstrated a high level of polymorphism among the strains tested.
The objective of the study was to develop a novel MLVA assay for genotypic study and to determine the presence of virulence-associated genes among Malaysian V. cholerae strains of different serogroups.

Bacteria Strains and DNA Extraction for Polymerase
Chain Reaction (PCR). Isolation of V. cholerae strains was carried out from June 2008 to April 2009. The identification and confirmation of the cultures were previously described in [11]. All together 43 V. cholerae strains were isolated from human stools (n = 15, from various public hospitals in Kuala Lumpur, Labuan, Kota Kinabalu, Kota Bharu, Tumpat, and Alor Setar), sea water (n = 13, from different locations in Klang, Banting, Morib, and Bachok), household wastewater (n = 6, from different locations in Petaling Jaya, Klang, Banting, Serdang, and Kajang), algae (n = 2, from costal water in Klang and Morib), and various raw seafood (n = 7) sold in wet markets in Petaling Jaya and used in this study. Crude DNA for each strain was prepared by suspending a loopful of bacterial colonies in 50 μl of deionized water and boiled at 99 • C for 5 minutes. After a quick centrifugation, 5 μl of the supernatant (∼50 ng of DNA template) was used for PCR.

Variable Number Tandem Repeats (VNTRs) Search and Primers Design.
A search for potential VNTR candidates in the genome of V. cholerae O1 biotype El Tor, N16961, was carried out (http://minisatellites.u-psud.fr/). Primers were designed for each VNTR based on a ∼500 bp sequence covering the flanking and tandem repeats region retrieved from the database. Primer 3 software [12] was used to design the primers. The specificity of the primers and VNTR polymorphisms were tested with in silico PCR (http://insilico.ehu.es/) and monoplex PCR using DNA templates from 5 Vibrio strains of different serogroups.
The PCR product was then diluted in 1 : 100 ratio with ddH 2 O. One microlitre of the dilution was mixed with 10 μl HiDi formamide and 0.2 μl GeneScan 500 LIZ size standard (Applied Biosystem). The mixture was boiled at 95 • C for 5 minutes and analyzed with ABI prism 3130xl Genetic Analyzer (Applied Biosystem). The amplicons were identified based on the peak height and colour while the size was determined with the GeneScan size standard using GeneMapper software (Applied Biosystem).

2.4.
Data Analysis for MLVA. The peak tables obtained from Gene Mapper were imported into BioNumerics 6.0 (Applied Maths, Belgium). The character values for each VNTR were determined and mapped. A Minimum Spanning Tree (MST) was also constructed using the categorical coefficient and WARD algorithm. The discrimination power of each VNTR was calculated individually and combined using Simpson's index of diversity [13].

Pulsed-Field Gel Electrophoresis (PFGE)
. PFGE was performed according to established protocol from PulseNet USA [14]. Briefly, equal volumes of standardized cell suspension (OD 610 = 0.8) of 1% (w/v) Seakem Gold agarose (Cambrex Bio Science Rockland, Inc, USA) were mixed to form plugs. The plugs were lysed in cell lysis buffer (50 mM Tris; 50 mM EDTA [pH 8.0], 1% Sarcosyl, and 1 mg/ml proteinase K) and incubated at 54 • C for 3 hours. The plugs were then washed thoroughly with sterile deionised water (twice) and TE buffer (6 times). A slice of DNA plug was digested overnight with 10 U restriction enzyme Not I (5 -GCGGCCGC-3 ) (Promega, Madison, Wis., USA) at 37 • C. The digested chromosomal DNA was subjected to PFGE on a 1.0 % (w/v) agarose gel (Sigma Type 1, St. Louis, Mo) with the CHEF DRIII (Bio-Rad, Hercules, CA) using the following conditions: 200 V for 13 hours with pulse times of 2-10 s and 6 hours with pulse time 20-25 s. Gel image was captured by using Gel Doc XR after staining with ethidium bromide (0.5 μg/ml).
PFGE banding patterns were analyzed with BioNumerics 6.0 (Applied Maths, Kortrijk, Belgium). The quantitative differences among the banding patterns were defined by the Dice coefficient; F. Cluster analysis was carried out based on the unweighted pair group method with arithmetic averages (UPGMA) using the position tolerance of 0.15. Simpson's index of diversity was also calculated [13].

Virulence Genes Detection.
Detections of ctxA, zot, rtxA, rstR, toxR, tcpA, tcpI, and hlyA genes were carried out via PCR using published primers [6,[15][16][17]. The genes selected were either major virulence genes or virulence-associated genes. The PCR amplifications were carried out in a Master cycler (Eppendorf, USA). The PCR mixture in a total volume of 25 μl contained 50 ng of DNA template, 1X PCR buffer, 2 mM MgCl 2 , 200 μM of each dNTP, 0.3 μM of each primer, Journal of Biomedicine and Biotechnology 3 and 2.5 U of Taq DNA polymerase (Promega, USA). The cycling conditions were set at 94 • C for 3 minutes (1 cycle), 94 • C for 30 s, 60 • C for 30 s, 72 • C for 1 minute (30 cycles), and 72 • C for 10 minutes (1 cycle). The products were then analyzed on 1.5% (w/v) agarose gel and run at 90 V for 1 hour. Gel images were captured and analyzed using Gel Doc XR (Bio Rad, USA). Representative amplicons of each gene were purified and sent to a commercial facility for sequencing to confirm the identity of amplicons.

MLVA.
Seven VNTR regions were chosen, and primers were designed based on the sequences retrieved from GenBank. Initially, the primers specificity was determined individually with monoplex PCR using 2 O1, 1 O139, and 2 non-O1/non-O139 strains. All 7 VNTR regions showed polymorphisms for those tested strains.
By using capillary electrophoresis, the sizes of DNA fragments could be determined more accurately based on the peaks. Overall, the size distributions of the 43 V. Thirty-eight MLVA profiles (F = 0.63) were obtained based on the character types (allelic profiles) ( Figure 3). The 43 strains were grouped into 3 major clusters in the MST (Figure 1). The first cluster (I) comprised of 14 O1 strains (8 clinical, 6 water) (90.6% similarity). In this cluster, 5 of the O1 strains shared the same profile (MLVA2: 09-08-09-12-05-02-23). These strains (VC4, VC52, VC63, VC65, and VC66) were isolated during a cholera outbreak in Kota Bahru, West Malaysia (n = 2) and Labuan, East Malaysia (n = 3) in 2008. Most of the environmental O1 strains were in this cluster. One environmental O1 strain, VC70, was highly similar (95 % similarity) with the O1 outbreak strains.
Simpson's index of diversity for each individual VCTR ranged from 0.59 to 0.92 (Table 1) while the index for the combined VCTR loci was 0.99.
When the strains were grouped according to serogroups and sources (clinical O1, environmental O1, environmental O139, and environmental non-O1/non-O139), significant difference in the aggregation was observed (one way-MANOVA F, Wilkes' λ = 0.0874, F 6,76 = 30.1700, and P < .0001). Further Post Hoc Tukey HSD test identified that there were significant differences between non-O1/non-O139 environmental Vibrio cholerae and the other 3 groups (P < .01). This was pertinent to the fact that all non-O1/non-139 strains lacked 4 virulent genes (ctxA, zot, tcpA, and tcpI). There was no significant difference in the virulotypes among the clinical O1, environmental O1, and environmental O139.

Discussion
The MLVA assay was developed based on the 7 loci in chromosomes I and II of V. cholerae. Among these loci, VCTR4 demonstrated the highest discriminatory ability (D = 0.92), followed by VCTR7 (D = 0.89). Both loci are situated in the small chromosome of V. cholerae genome (chromosome II) and suggest that there is a higher genetic variation in chromosome II of V. cholerae genome compared with chromosome I. Danin-Poleg et al. [18] and Olsen et al. [10] reported that the tandem repeat locus that showed the highest discriminating ability is physically located at 187 kb in chromosome II (such as VCTR7 in this study). However in our study, VCTR4 which is located at 303 kb in chromosome II showed even higher genetic variation (17 alleles, D = 0.92) among the 43 strains tested.
Both the PFGE and MLVA data showed that the V. cholerae O139 strain (VC4370) was closely related with an O1 strain with similarity of more than 90%. Strains VCGB and VCSW, isolated from river water and algae, respectively, were indistinguishable by both PFGE and MLVA. Both strains harbored the same virulence genes (rtxA + , rstR + , toxR + , and hlyA E+ ) although they were isolated from different locations. VCGB was isolated from river water in Banting, Selangor (August 2008) while VCSW was isolated from algae From the analysis of MST, those strains which belonged to the same PFGE profile (N11, N13, N14, and N16) were distinguishable by MLVA. However, some strains with same MLVA profiles could be further distinguished by PFGE.  Hence, a combination of PFGE and MLVA analysis may yield more information about the clonality of bacterial pathogens. For example, 3 environmental strains (VC87, VC88, and VC90) and 1 clinical strain (VC66) had identical pulsotype but were distinguishable by MLVA. Five strains (VC4, VC52, VC63, VC65, and VC66) were indistinguishable by MLVA as they shared identical allelic profile. However, PFGE was able to separate 2 of Kota Bahru outbreak strains (VC4 and VC65) (West Malaysia) from the Labuan outbreak strains in East Malaysia.
Cholera cases caused by O1 serogroup occur occasionally in Malaysia and the neighboring countries. PFGE is the most commonly used subtyping method to determine the epidemiological relatedness of the strains. In this study, the 23 O1 strains were subtyped into 18 pulsotypes. Based on visual comparison of the published pulsotypes [19], our pulsotypes N22 and N11 appeared to be similar to those pulsotypes of Hong Kong strains isolated during an outbreak in 2001. Pulsotype N13 which comprised of 3 environmental strains (VC87, VC88, and VC90) and 1 clinical strain (VC66) was similar to the pulsotype of a Philippine strain obtained during an outbreak in 2001. The Malaysian strains isolated from sea water had similar pulsotype with those Indonesian outbreak strains in 1995 [19]. Overall, comparison of the pulsotype distributions among V. cholerae in different neighboring countries would provide important information and further facilitate the epidemiological study of cholera in Southeast Asia. However, as different PFGE conditions were used by different researchers in the region, direct comparison was difficult. Adoption of a standardized PFGE protocol such as the PulseNet PFGE protocol proposed by CDC PulseNet, USA would greatly enhance interlaboratory comparison and improve tracking of V. cholerae strains among the endemic countries in the region.
In this study, 8 environmental strains harbored all virulence-associated genes, indicating the toxigenic potential and risk of infection from environmental sources [1]. The presence of accessory toxigenic genes in non-O1/non-O139 strains indicates the relative importance of these strains as they have the potential to cause infections [20].
VC 35 was isolated from a patient during a cholera outbreak in Kota Kinabalu, Sabah. This strain did not harbor the major virulence genes (ctxA, zot, tcpA, and tcpI) which are essential for colonization. Hence, the finding was in agreement with the finding of Olivier et al. [8] that secreted accessory toxins such as hemolysin and RTX toxin in El Tor V. cholerae would modify the host environment and prolong colonization without cholera toxin (CT) or toxin-coregulated pili (TCP). Overall, the MLVA assay developed in this study was a promising method to study the relatedness of Malaysian V. cholera, and the discriminatory ability was comparable with PFGE. In addition, further combination of MLVA and PFGE analysis may provide more information of the relatedness of strains and distinguish strains from different sources and geographical regions.