Relatedness of Vibrio cholerae O1/O139 Isolates from Patients and Their Household Contacts,Determined by Multilocus Variable-Number Tandem-Repeat Analysis

The genetic relatedness of Vibrio cholerae O1/O139 isolates obtained from 100 patients and 146 of their household contacts in Dhaka, Bangladesh, between 2002 and 2005 was assessed by multilocus variable-number tandem-repeat analysis. Isolate genotypes were analyzed at five loci containing tandem repeats. Across the population, as well as within households, isolates with identical genotypes were clustered in time. Isolates from individuals within the same household were more likely to have similar or identical genotypes than were isolates from different households, but even within a household, isolates from different individuals often had different genotypes. When household contacts were sampled regularly for 3 weeks after the illness of the household index patient, isolates with genotypes related to the index patient appeared in contacts, on average, ∼3 days after the index patient, while isolates with unrelated genotypes appeared in contacts ∼6 days after. Limited data revealed that multiple isolates from the same individual collected within days of each other or even from a single stool sample may have identical, similar, or unrelated genotypes as well. Our results demonstrate that genetically related V. cholerae strains cluster in local outbreaks but also suggest that multiple distinct strains of V. cholerae O1 may circulate simultaneously within a household.Vibrio cholerae is the etiologic agent of cholera, a secretory diarrheal disease with a high mortality rate in humans if untreated (25). Serogroups of V. cholerae, a motile, Gram-negative, curved rod, can be defined serologically by the O side chain of the lipopolysaccharide (LPS) component of the outer membrane (9). V. cholerae is found in a variety of forms in aquatic ecosystems (41, 42), and more than 200 different serogroups have been isolated, mostly from environmental sources (45). However, the vast majority of V. cholerae strains that cause the clinical disease cholera belong to serogroup O1 or O139 (37, 42). V. cholerae O1, the historical agent of epidemic and pandemic cholera and the current leading cause of cholera both globally and in Bangladesh (42), is classified into two major biotypes, classical and El Tor (44), and two major serotypes, Ogawa and Inaba (48). The current global pandemic is caused by V. cholerae O1 El Tor. A second pathogenic serogroup, O139, emerged in the Bengal region in 1992 by horizontal transfer of new LPS biosynthesis-encoding genes into the El Tor biotype (1, 4). This new serogroup continues to cocirculate with El Tor V. cholerae O1 serotypes Ogawa and Inaba as a cause of disease in humans, although it accounts for a smaller proportion of all cholera now than in its first years of circulation (16, 20). Recently, comparative genomics has revealed an extensive amount of lateral gene transfer between strains, suggesting that genomic classification may be an alternative to serogrouping for classifying pathogenic V. cholerae strains (11).Toxigenic V. cholerae may be present in environmental sources in regions of endemicity and emerge, often seasonally, to cause cholera in humans (12, 18). Once an outbreak has begun, organisms from one infected individual are more infectious for the next individual, a property termed hyperinfectivity, and these forms may be able to pass directly from human to human through fecal-oral contamination (35). However, because vibrio organisms are difficult to isolate from implicated environmental or domestic water sources (28, 29), little is known about the diversity of V. cholerae in inocula that cause human infection.Established laboratory methods for differentiating V. cholerae strains, apart from serogrouping and serotyping, include rRNA restriction fragment length polymorphism (ribotyping), pulsed-field gel electrophoresis (PFGE), and multilocus sequence typing (MLST). These methods, however, have a limited capacity to differentiate between pathogenic V. cholerae strains, as clinical isolates are relatively genetically monomorphic. For instance, V. cholerae O1 comprises approximately 30 ribotypes (39); however, only a few ribotypes are common in clinical isolates, ribotypes evolve slowly, and all isolates of a given pathogenic V. cholerae serotype in a local area over a period of multiple years often belong to a single ribotype (8, 14, 17). In a broad sampling of 154 V. cholerae isolates from Bangladesh and worldwide over several decades, only 15 ribotypes were identified, and of these, many were found in nonpathogenic environmental isolates only; only five ribotypes were associated with the V. cholerae O1 El Tor biotype that currently predominates as the cause of clinical disease, while pathogenic isolates of serogroup O139 were indistinguishable from each other by ribotype (19).PFGE, in which restriction endonuclease digestion of genomic DNA generates mutation-sensitive banding patterns, is often more sensitive than ribotyping in detecting strain variation (7, 34, 51) and detects extensive genetic variation within nonpathogenic V. cholerae serogroups (3, 46). However, PFGE types change slowly and are useful primarily for distinguishing between strains in different pandemics or between different continental branches of those pandemics. In an analysis of 180 mostly western-hemisphere isolates (7), PFGE differences had developed from a prior pandemic strain over the 30 years since its arrival in Latin America, but a new strain that had been causing disease for 2 years still had only a single PFGE type across the 64 isolates analyzed. Similarly, in a Japanese study (2), although 19 PFGE types were identified among O1 isolates, the majority of the domestic isolates, along with several imported isolates, belonged to a single PFGE type.Further differentiation between V. cholerae isolates is achievable by MLST, which characterizes isolates by internal DNA sequences in selected housekeeping genes (32). Nevertheless, epidemic strains also cluster tightly in this typing scheme (5, 32) and the method has been useful primarily for determining relationships between nontoxigenic strains (36) or for linking regional outbreaks (which typically appear monoclonal by these methods) with the pandemic strain responsible (5, 33).Although these methods have distinguished major pandemic clones from other nonpathogenic human and environmental isolates of V. cholerae, the near clonality of pathogenic O1 and O139 strains means that established methods may not provide sufficiently robust differentiation of these genetically similar pathogenic strains to answer important epidemiological questions. Therefore, there is a need for other methods that can distinguish among clinical O1 and O139 isolates and track the epidemiology of outbreaks in a restricted geographic area on a shorter time scale.Multilocus variable-number tandem-repeat (VNTR) analysis (MLVA) is one method that may be useful for differentiating between pathogenic V. cholerae O1 and O139 strains that would be indistinguishable by other techniques (15). This method examines short repeating DNA segments at various locations in the genome that can vary in number at each location and uses the number of repeats at each varying locus as a fingerprint to distinguish between isolates.Escherichia coli is the paradigm organism for demonstrating the value of the MLVA method. Noller et al. (38) showed that E. coli O157 isolates that were indistinguishable by MLST could be distinguished to some extent by PFGE but that MLVA distinguished between isolates that had the same PFGE type and did so in a manner consistent with the known epidemiology of the isolates (38a). In addition, machine-scored VNTR assays have been demonstrated to be robust and portable and to discriminate clearly between isolates by using relatively few loci, therefore limiting the effect of compounding genotyping errors (6).For V. cholerae, five VNTR loci have been identified (15), and the initial application of MLVA at those loci has demonstrated distinct populations of clinical isolates of V. cholerae in different geographic regions within Bangladesh and India (23, 47). Predominant isolates in each of two rural Bangladeshi regions varied gradually over a time scale of months to years (47), and isolates collected from India over a 15-year period varied widely, with individual MLVA types clustering in time and place—some with widespread dissemination and others with limited local occurrence only (23). MLVA has also been used to classify hybrid and altered V. cholerae variants and to demonstrate their genetic distance from the pandemic El Tor strain (10). Use of the MLVA method for epidemiologic study of cholera requires that V. cholerae VNTR alleles remain reasonably stable during bacterial replication in patients or in laboratory culture after isolation. Some degree of stability of two of the five loci used in V. cholerae MLVA has been demonstrated previously by serial passage in vitro through four overnight cultures (15). In this study, we used MLVA to examine V. cholerae O1 and O139 isolates obtained from infected patients and their household contacts—including multiple isolates from the same individual and isolates from multiple individuals within the same household—in a large city where cholera is endemic.