Previous Article | Next Article ![]()
Journal of Clinical Microbiology, June 2007, p. 1728-1734, Vol. 45, No. 6
0095-1137/07/$08.00+0 doi:10.1128/JCM.00018-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Beijing Friendship HospitalAffiliate of Capital University of Medical Sciences, Beijing Tropical Medicine Research Institute, 95 Yong An Road, Beijing 100050, People's Republic of China,1 Department of Microbiology, Immunology, and Pathology, Colorado State University, Fort Collins, Colorado 805232
Received 4 January 2007/ Returned for modification 5 February 2007/ Accepted 28 March 2007
|
|
|---|
|
|
|---|
Wenshan Prefecture in Yunnan Province can be considered an area in China where leprosy is endemic because the new-case detection rate has decreased only moderately from 4.7 to 3.6 per 100,000 during 1985 to 2002. In particular, certain townships of Qiubei County in this prefecture have had detection rates as high as 9.7/100,000 in recent years (6). An epidemiological investigation conducted from 1998 to 2000 based on a serological test for antibodies against phenolic glycolipid I, a dominant M. leprae antigen, showed that over 30% of the general population tested positive. Furthermore, in 28% of household contacts of leprosy patients tested in two villages where leprosy is endemic, Tonghong and Nanqiu in Qiubei County, M. leprae DNA was detectable in nasal swabs (21). This study revealed a significant transmission of leprosy in this area. To further understand the transmission patterns, the present study takes advantage of novel molecular strain typing methods reported since the availability of the M. leprae genome sequence of a strain (TN) isolated from a patient in Tamil Nadu, India (2). Two types of genomic markers have been identified. Three single nucleotide polymorphisms (SNPs) discovered through the partial genome sequencing of an isolate from a patient from Brazil have been useful to track evolution and dispersion of M. leprae at global and species levels (11). On the other hand, the detection of multiple variable-number tandem-repeat (VNTR) loci with a wide range of allelic diversity should be more suitable for the differentiation of isolates within a community and thus useful for detecting recent transmission events (5, 25, 26). To describe M. leprae diversity and assess the value of multiple-locus VNTR analysis (MLVA), we present our findings for Qiubei County.
|
|
|---|
VNTR amplification and DNA sequencing of PCR products. The M. leprae VNTR locus nomenclature, primers for PCR amplification, and PCR conditions were described previously (5). The primers for the (AC)8a locus were described previously by Zhang et al. (26). During the course of the study, we implemented multiplex PCR in which certain primer sets were combined. The Multiplex PCR kit (QIAGEN) was utilized using 1 µl of a twofold-diluted template and cycled as follows: enzyme activation for 15 min at 95°C, followed by 40 cycles, each comprised of 30 s of denaturation at 94°C, 90 s of annealing at 60°C, and 90 s of extension at 72°C. A final extension step of 10 min at 72°C was performed prior to storage at 10°C. A 5-µl aliquot from the PCR amplification was applied to a 3% agarose gel (Invitrogen Corporation, Carlsbad, CA) in 1x Tris-borate-EDTA buffer and electrophoresed at 80 V for 90 min. PCR products were detected by staining the gel with ethidium bromide. If SYBR Safe DNA gel stain (Invitrogen) was used, it was added to the gel prior to electrophoresis. An EZ Load 20-bp molecular ruler (Bio-Rad, Hercules, CA) was run on lanes 1 and 8 of an eight-well gel. PCR-positive products (1 to 2 µl) were submitted to the BigDye Terminator v 3.1 cycle sequencing method using the forward primer at the Macromolecular Resources Core Facility at Colorado State University or at Boya Bio-Tech Company, Beijing, People's Republic of China. The VNTR allele at each locus was defined as the number of repeat units present for each product. The allele at each of the microsatellites and the 6-7 minisatellite locus were determined by direct sequencing of PCR products.
VNTR data analysis.
The Microsatellite Tool kit (available from the website http://oscar.gen.tcd.ie/
sdepark/ms-toolkit/) was used to calculate allele frequency. The allele diversity, H, for each locus was calculated by using the equation H = 1
fi2, where f is the frequency of the ith allele. To infer phylogeny, we performed parsimony analysis through the application of PAUP version 4.0 beta (17). All analyses were preformed by heuristic search procedures, with tree bisection-reconnection branch swapping and the MULTREES option. A 50% consensus tree was used to report tree topology. A weighted-step procedure was used based on the assumption that a loss or gain of a repeat consists of a single step. For example, if the number of repeats at a VNTR locus was 6 in one strain and was 9 in a second strain, this constituted three steps.
|
|
|---|
![]() View larger version (30K): [in a new window] |
FIG. 1. (A) Distribution of patients diagnosed from 2002 to 2005 in Qiubei County, Yunnan Province, People's Republic of China. The capital letters represent the names of the townships of this county. XD, Xindian; NJ, Nijiao; YZ, Yuezhe; BDS, Badaoshao; SP, Shupi; SLY, Shuanglongying; JP, Jinping; TX, Tianxing; PZ, Pingzai; WL, Wenliu; YJ, Yangjie. The inset depicts the Yunnan Province in southern China. Patients are coded by numbers. Multicase families (F1 to F11) are shown as rectangles, and the patient numbers are indicated. (B) The distribution of characteristic strain types A to E identified by MLVA (in Fig. 2) is depicted for some of the patients represented in panel A.
|
![]() View larger version (27K): [in a new window] |
FIG. 2. Phylogenetic analysis of M. leprae strains isolated from leprosy patients in Qiubei County, Yunnan Province, People's Republic of China. A 50% consensus tree is illustrated, which was obtained with VNTR alleles for eight loci for 57 isolates using parsimony principles as described in Materials and Methods. The numbers at the nodes of branches represent percentages of trees in which the nodes were found. To the right of the tree, the eight locus alleles (repeat number) for each patient are indicated. Alleles 1 to 9 are indicated by the number, while a letter code is used for alleles greater than 9. Alleles 22, 26, 31, 32, 33, 38, 39, 40, 41, and 42 were not found in the data set. Thus, the alleles could be accommodated using 26 letter codes as follows: A to L, 10 to 21, respectively; M to O, 23 to 25, respectively; P to S, 27 to 30, respectively; T to W, 34 to 37, respectively; X and Y, 43 and 44, respectively; and Z, all alleles 45. The township of residence is listed after the alleles using capital letters (as described in the legend of Fig. 1), followed by the patient number. The last column indicates the family number for multicase families. The rectangles marked with the letters A to C represent VNTR profiles grouped on the basis of similarity between the strains. The (GTA)9 allele is 9 in clusters A and B and 11 to 13 in cluster C. Strains grouped in clusters D and E are more variable but distinguished by larger (GTA)9 alleles: 15 to 25 in cluster D and >26 in cluster E.
|
|
View this table: [in a new window] |
TABLE 1. Allelic properties of the VNTR loci in M. leprae isolates in leprosy patients from Qiubei Countya
|
Distribution of strain types in Qiubei County. Combining the distribution of patients, as illustrated in Fig. 1, with the phylogenetic analysis (Fig. 2) revealed that more than one strain type can be seen in each township (Fig. 1B). Although the number of patients from each township is limited (1 to 10 patients), it can be seen that group A and group B isolates are distributed mainly in the north part of the county (townships GH, GZ, and YZ) (Fig. 1B). On the other hand, patients in the southern and eastern parts (townships Tianxing, Shupi, Jinping, Wenliu, and Badaoshao) harbor the C, D, and E types of strains.
Transmission scenarios in Qiubei County and the influence of multicase families. In Qiubei County, intrafamilial transmission appears to be common. In each of groups A to E, there is at least one multicase family. Within intrafamilial patient pairs (F1 [patients 237 and 135], F3, F5, F9, and F10), matching alleles were detected at five or more of the eight loci, indicating a common infectious source and recent transmission within the household. The allelic differences of one to two copies within each pair of the multicase families are seen mostly at the AT-rich loci. Allelic differences across several loci seen within patient pairs of families F6 and F7 suggest that the family members harbor different strains.
Furthermore, we have identified cases in the same or adjoining townships that harbor closely matching genotypes. For example, closely matching M. leprae VNTR patterns were seen in cases 239 and 245 in the GZ township, in the F1 multicase family (patients 237 and 135), and in case 240 in the adjacent township YZ. Similarly, in GH, one member of F6 (patient 277) shared VNTR pattern similarities with two other patients (patients 241 and 271) in the same township as well as with patient 244 in the neighboring GZ township (Fig. 2).
It is worth noting that there are three multicase families in SD: F3 (patients 136 and 246), F4 (patients 73 and 154), and F5 (patients 247 and 249). Families F3 and F4 are similar, both belonging to group C, while F5 is quite different (group E). F3 and F4 cluster with case 268 of the adjoining GZ township.
A trade fair held weekly in the SD township is a meeting place for people from GZ and YZ and serves as a main mode of interaction among the local people. The average distance from the outside villages to SD is 10 km. We postulate that the trade fair may be a means of transmission to other places, although other possibilities must be sought and investigated.
|
|
|---|
With regard to the identification of cluster-related genotypes and close tracking of transmissible isolates of an infectious disease, it has been widely recognized that genetic markers with sufficient heterogeneity within the population are necessary to distinguish multiple transmission clusters yet should be stable enough to link isolates belonging to common recent transmission events (15). Highly stable markers with limited allelic diversity in the study population such as SNPs are not suitable for the study of short-term transmission (7). Recently, three SNPs were used to develop a model of the origin and dispersion of M. leprae (11). The isolates from China belong to one (SNP type 3) of the four major types and represent an alternate eastward branch/route of transmission from the proposed origin in East Africa/Central Asia. In this study, four randomly selected isolates from different townships, 78 (Badaoshao), 79 (Pingzai), 80 (Nijiao), and 84 (Shupi), are all SNP type 3. Additional isolates from prior sample collections in China were also shown to be SNP type 3 (data not shown).
In recent years, short tandem-repeat sequences that generate polymorphisms by the mechanisms of strand slippage and recombination have been widely accepted as being useful markers for epidemiology and transmission surveillance of multiple pathogens. MLVA has been compared with other traditional methods such as restriction fragment length polymorphism, spoligotyping, and multilocus sequence typing and was found to be suitable based on discriminatory potential and convenience in laboratory applications (1, 3, 8, 13, 16). Thus, after the sequence of the M. leprae genome was published, a list of 44 candidate VNTR loci was identified, 9 of which showed polymorphisms in a small panel of reference strains (5). Subsequently, VNTRs were demonstrated in clinical isolates in these and other loci (25, 26). Stability in M. leprae microsatellite loci has been demonstrated during multiple passages in nude mice and armadillo hosts, with an occasional shift by one or two copies in some loci (10, 18).
Although VNTRs are accepted markers for strain typing, the optimal model for the quantification of allelic changes during expansion or contraction of repeats is not known. Furthermore, the process and patterns of mutations at different loci may differ from locus to locus, depending on the sequence and the length of the repeat motif and the array size at each locus (4, 12). We have assumed a stepwise model of mutation for the expansion or contraction of repeat numbers at VNTR loci in order to construe genetic relatedness between isolates. Recent studies of two bacterial species amenable to multiple in vitro passages for measurements of mutation rates have indeed demonstrated that single-repeat insertions and deletions are equally probable and that they constitute a significant proportion of the mutations (
80%) responsible for strain evolution (9, 19, 20). Furthermore, the mutation rate tends to depend on allele diversity.
A survey of four microsatellites and six minisatellites on 50 M. leprae isolates from China revealed that two of the microsatellites [(AC)8b and (GGT)5] and five minisatellites (12-5, 18-8, 21-3, 23-3, and 27-5) exhibited no or minimal polymorphisms (22). Therefore, in the present study, the polymorphic microsatellites (AC)9 and (GTA)9 were combined with one minisatellite locus (6-7) and five other candidate loci to resolve prevalent strains. In developing this approach, results from previous field studies were taken into consideration (10, 25). When only three microsatellite VNTR loci were applied to investigate the genetic diversity of M. leprae in India, only 2 out of 33 isolates were found to have the same VNTR pattern, with no known epidemiological link. However, within a known multicase family, these three loci were informative and identical (25). Therefore, it was clear that several loci are required to discern strain identity, but which loci would prove to be informative was not known. The study reported previously by Zhang et al. utilized isolates collected mainly from Japan and a few others from Thailand, Korea, and Indonesia (26). Those strains were collected at different times and maintained in mice by multiple passages. Although many loci were mapped, they were not ideally suited for finding relationships relevant to local transmission. However, striking strain identity was detected when multiple patients in households were studied, indicating infection with the same strain.
Our present study was therefore focused on newly diagnosed leprosy patients residing in the small area of Qiubei County over a 3-year period. MLVA was shown to be more meaningful than strain typing with a single locus or a few loci. We illustrate that if only loci with many alleles, such as (GAA)21 and (TA)18, are used as first-level strain-typing markers, it is likely that related strains may not be recognized. Furthermore, microevolution at microsatellite loci is evident in Qiubei, even in strains infecting patients in the same household. However, by including six other loci, even some with limited allelic diversity, major groupings were discernible. Therefore, in addition to having a sufficient number of loci, it may be necessary to pay attention to the allelic diversity of loci in the prevalent population.
The study depicts several characteristic features of M. leprae from patients in Qiubei County. First, the clustering of patients within families is severe in Qiubei County. Multicase families, sharing closely matching VNTR genotypes, indicate ongoing transmission prior to detection and treatment with MDT. One or two unit differences in copy number are expected in loci prone to high rates of mutation, such as the AT-rich regions, due to the long incubation period and unknown number of transmission cycles before the onset of disease. We propose that multicase families may constitute epidemic foci and a source of M. leprae in villages; therefore the predominant strain or cluster(s) tends to be that found in multicase families (Fig. 1 and 2). The high level of transmission involving families is indicated by the finding that three other multicase households (F12, F13, and F14) have been identified in the GZ township since early 2006. Of these, F13 and F14 arose from previously defined single-case families. Furthermore, new cases were found in existing multicase families (as in F1).
A long term-study is required to monitor the evolution and transmission of M. leprae in Qiubei and to verify the trends in data from 2002 to 2005 deduced via this MLVA molecular epidemiology approach. The mechanism of expansion of the (GTA)9 locus in isolates with large allele sizes as seen in group E deserves further study.
We acknowledge NIH/NIAID grants NO1 AI-25469, RO3-AI-59644, and RO1-AI-63457.
Published ahead of print on 11 April 2007. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»