Previous Article | Next Article ![]()
Journal of Clinical Microbiology, May 2007, p. 1483-1490, Vol. 45, No. 5
0095-1137/07/$08.00+0 doi:10.1128/JCM.02191-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

DST/NRF Centre of Excellence for Biomedical Tuberculosis Research/MRC Centre for Molecular and Cellular Biology, Division of Molecular Biology and Human Genetics, Department of Biomedical Sciences,1 Desmond Tutu Tuberculosis Centre, Department of Paediatrics and Child Health,2 Faculty of Health Sciences, Stellenbosch University, Tygerberg, South Africa, and Centre for Statistical Consultation, Stellenbosch University, Tygerberg, South Africa3
Received 26 October 2006/ Returned for modification 5 February 2007/ Accepted 6 March 2007
|
|
|---|
|
|
|---|
M. tuberculosis strains with a Beijing genotype are members of principal genetic group 1 and are characterized by a distinct spoligotype pattern, closely related IS6110 DNA fingerprints, identical variable-number tandem repeat sequences, and an IS6110 insertion in the NTF 1 region (Rv0001:Rv0002) (10). Comparative genomics have shown that Beijing strains have evolved through mechanisms including IS6110 transposition (2), SNPs in the mismatch repair genes (15), deletion of chromosomal domains (regions of difference [RDs]) (19), and synonymous SNPs (sSNPs) (9). An evolutionary scenario based on mutations in the mismatch repair genes mutT2 and mutT4 and ogt genes has suggested the absence of these mutations in the ancestral Beijing strains (15). Subsequent mutations in the mutT4 gene (codon 48) and ogt gene (codon 37) or ogt and mutT2 genes (codons 12 and 58, respectively) generated two independently evolving subgroups (15). Mapping of the RD domains showed that the Beijing strain family is a monophyletic clade, with four distinct evolutionary branches defined by the RD105, RD181, RD150, and RD142 deletions, respectively (19). Phylogenetic analysis according to sSNPs supported the monophyletic clade structure and suggested 11 Beijing sequence types (STs) (9). However, the order in which all of the above evolutionary events occurred remained to be determined.
In this study we aimed to reconstruct the evolutionary history of Beijing strains from the Western Cape, South Africa, through the analysis of the genomic structure of 40 different loci. This phylogeny was used to test the hypothesis that evolution has influenced the ability of the Beijing strains within the different Beijing sublineages to spread and cause disease.
|
|
|---|
Rural setting. Sputum isolates were routinely collected at diagnosis from all high-risk tuberculosis patients (patients whose sputum failed to convert after the intensive phase of therapy or patients who had previously been treated with antituberculosis drugs) who were attending primary health care clinics and who were resident in the Boland, Overberg, Karoo, and Southern Cape region of the Western Cape, South Africa, during the period February 2003 to August 2004. This region is located 150 to 500 km northeast of the epidemiological field site in Cape Town. Sputum smear microscopy (fluorescent staining) and/or culture (BACTEC 460, MGIT 960, or Löwenstein-Jensen medium) was performed by the National Health Laboratory Service or the Stellenbosch University laboratories to confirm the presence of M. tuberculosis. Clinical and demographic data including previous history of tuberculosis, gender and age, smear positivity, and drug susceptibility test results were recorded in a database. A survey of 384 tuberculosis retreatment patients (2003 to 2004) in this study setting showed that 13% were HIV coinfected. This study setting was included to enable a comparison of the M. tuberculosis strain population structures in two independent communities.
This study was approved by the Ethics Committee (Institutional Review Board) of Stellenbosch University.
Drug susceptibility testing. Drug susceptibility testing was done by the National Health Laboratory Service, using the indirect proportion method on Middlebrook medium containing critical concentrations of 0.2 µg/ml isoniazid and 30 µg/ml rifampin. In this study drug resistance was defined as resistance to either isoniazid or rifampin or a combination of isoniazid and rifampin (multidrug-resistant [MDR] tuberculosis). Drug susceptibility testing for other resistance markers was not routinely done.
Genotypic classification. Sputum isolates from patients resident in the urban epidemiological field site were subcultured on Löwenstein-Jensen medium and genotyped by IS6110 DNA fingerprinting (22) and spoligotyping according to internationally standardized methods (11). Sputum isolates from patients resident in the rural setting were cultured on BACTEC 12B medium, and a boiled aliquot was genotyped by spoligotyping according to internationally standardized methods (11). Subcultures of the isolates from the rural setting were not available for IS6110 DNA fingerprinting.
Isolates were classified as members of the Beijing genotype if they had the characteristic Beijing spoligotype (10). Only the first M. tuberculosis isolate from each case with a Beijing strain was included for subsequent genetic analysis.
Genomic comparison. Each Beijing strain from each patient was subjected to PCR amplification in a reaction mixture containing 0.2 µg DNA template, 5 µl Q buffer, 2.5 µl 10x buffer, 2 µl 25 mM MgCl2, 4 µl 10 mM deoxynucleoside triphosphates, 1 µl of each primer (50 pmol/µl) (Table 1), and 0.125 µl HotStarTaq DNA polymerase (QIAGEN, Germany) and made up to 25 µl with H2O. Primers 1 to 22 (Table 1) were used in combination with the universal forward primer and the internal control primers to determine the presence or absence of IS6110 elements at specific chromosomal loci (2). Primer set 23 to 26 (Table 1) was used to determine the presence or absence of SNPs in the mutT2 (codon 58), mutT4 (codon 48), and ogt (codons 12 and 37) genes (15) by the amplification refractory mutation system method. Primer set 27 to 30 (Table 1) was used to determine whether RDs (RD105, RD181, RD150, and RD142) (19) were present or absent. Primer set 31 to 40 (Table 1) was used to determine the presence or absence of sSNPs at chromosomal positions 797736, 909166, 1477596, 1548149, 1692069, 1892017, 2376135, 2532616, 2825581, and 4137829 (9) relative to the H37Rv whole-genome sequence (6) by the amplification refractory mutation system method. Amplification was initiated by incubation at 95°C for 15 min, followed by 35 to 45 cycles at 94°C for 1 min, annealing temperature (Table 1) for 1 min, and 72°C for 1 min. After the last cycle, the samples were incubated at 72°C for 10 min. PCR amplification products were electrophoretically fractionated in 3.0% agarose in 1x Tris-buffered EDTA, pH 8.3, at 6 V/cm for 4 h and visualized by staining with ethidium bromide. The existence of a mutational event was determined by the presence or absence and the size of the respective PCR product.
|
View this table: [in a new window] |
TABLE 1. Primer sequences
|
Statistical analysis. The z test for proportions was used for testing frequency of occurrence of cases in the respective Beijing sublineages. The Fisher exact test was used to determine the association between Beijing sublineage and gender, new versus retreatment cases, smear positivity within 2 months of diagnosis, drug resistance (for patients resident in either the urban or the rural setting), and IS6110 clustering (urban setting only). In all cases, a P value of 0.05 was used as the cutoff level for significance.
|
|
|---|
Phylogenetic analysis according to the neighbor joining algorithm grouped the Beijing strains from the two settings into four STs, ST11, ST19, ST22, and a novel ST (Beijing sublineages 6 and 7) (Fig. 1, left). The evolutionary history of these STs was further resolved according to SNPs in the mismatch repair genes, mutations induced by IS6110 insertion, and the structure of the RDs (Fig. 1, right). This phylogenetic tree was statistically robust, and the predicted evolutionary order was congruent with the tree based only on sSNPs (compare left and right sides of Fig. 1). This phylogeny demonstrates that the Beijing genotype is a monophyletic clade, defined by the RD105 deletion; IS6110 insertions at chromosome positions 1592, 1986638, 3127932, 3549199, and 3844681; and an sSNP at position 1548149 (Fig. 1, right). Subsequent evolution divided the Beijing strains into seven Beijing sublineages (Fig. 1, right) according to RDs (19), sSNPs (9), SNPs in the mismatch repair genes (15), and IS6110 insertions (2). The evolutionary order suggests that the RD181 deletion and the mutT4 (codon 48) SNP occurred before ST19, while the SNPs in the mutT2 and ogt mismatch repair genes occurred after ST19 and before ST10. Beijing sublineages 3 to 6 were distinguished by a combination of specific IS6110 insertions and/or sSNPs, while Beijing sublineage 7 was characterized by the RD150 deletion. The overall topology of the tree suggests divergent evolution, and this is supported by the increasing complexity in the IS6110 banding patterns between Beijing sublineage 1 (closest to the ancestor) and Beijing sublineage 7 (most distant from the ancestor) isolates from the urban setting (Fig. 2).
![]() View larger version (14K): [in a new window] |
FIG. 1. The evolutionary history of M. tuberculosis Beijing strains from an urban and a rural setting in the Western Cape, South Africa. (Left) Phylogenetic reconstruction based on sSNPs (9) using the bootstrapping and neighbor joining algorithms (PAUP* 4.0). The tree was rooted to M. bovis. Bootstrap values are given at internal nodes. The chromosomal positions of the predominant sSNPs are given at the nodes where they occur. All branches with a zero length were collapsed. (Right) Phylogenetic reconstruction based on sSNPs (9), SNPs in mismatch repair genes (15), IS6110 insertion points (2), and RDs (19) using the bootstrapping and neighbor joining algorithms (PAUP* 4.0). The tree was rooted to M. bovis. Bootstrap values are given at internal nodes. The respective evolutionary events are indicated at nodes where they occur (the order in which these events occurred is unknown). All branches with a zero length were collapsed. The scale indicates the number of steps per unit length.
|
![]() View larger version (29K): [in a new window] |
FIG. 2. IS6110 DNA fingerprints of M. tuberculosis Beijing strains from the urban setting. M. tuberculosis isolates were classified by IS6110 DNA fingerprinting according to the internationally standardized method (22). The DNA fingerprints from representative isolates are grouped according to their Beijing sublineage. Only a single representative of each IS6110 DNA fingerprint is shown.
|
In order to determine whether the evolutionary process had altered the ability of drug-sensitive Beijing strains to spread and cause disease, the frequency of occurrence of strains from the more recently evolved Beijing sublineage 7 was compared to that of strains belonging to the more distantly evolved Beijing sublineages 2 to 6. This analysis showed an overabundance of Beijing sublineage 7 strains in the urban population (Table 2) (P < 0.001 from the z test for the hypothesis that proportion of sublineage 7 cases = 0.14) as well as in the rural population (Table 3) (P < 0.001 from the z test for the hypothesis that proportion of sublineage 7 cases = 0.14). In addition, strain clustering according to IS6110 DNA fingerprinting in the urban community was strongly associated with the Beijing sublineage 7 (Table 2) (odds ratio, 5.82; 95% confidence interval, 3.13 to 10.82 [P < 0.001]). (A similar analysis could not be done for isolates from the rural setting as cultures were not available for DNA fingerprinting.) Together, this suggests that strains from sublineage 7 have evolved a phenotype which allows them to transmit and cause disease more frequently than strains from sublineages 2 to 6.
|
View this table: [in a new window] |
TABLE 2. Clinical, demographic, and molecular epidemiological data from patients resident in the urban epidemiological field site in Cape Town, South Africa
|
|
View this table: [in a new window] |
TABLE 3. Clinical and demographic data from patients resident in the rural setting of the Boland, Overberg, Karoo, and Southern Cape region of the Western Cape, South Africa
|
To determine whether evolution has altered the ability of the Beijing sublineages to acquire drug resistance, phenotypic drug susceptibility data were compared according to the respective strains' Beijing sublineage. The results from the 50 patients resident in the urban setting showed no association between the Beijing sublineage and the number of strains with drug resistance (Table 2). However, a high level of transmission of resistant strains was seen in Beijing sublineage 5, which was associated with MDR tuberculosis (Table 2). Similarly, the high number of drug-resistant tuberculosis cases seen with Beijing sublineages 6 and 7 in the rural setting was associated with MDR tuberculosis (Table 3).
|
|
|---|
In this study, cases with a sublineage 7 strain were overrepresented, suggesting that the evolutionary process has had a positive influence on the strain's ability to spread and cause disease. This concept was further supported by an association between sublineage 7 and transmissibility. Similarly, in a previous study, a significant association was observed between extrathoracic tuberculosis and a subset of Beijing strains characterized by either RD142 or RD150 deletion (12). Together these data support the hypothesis that phenotypic differences may exist within the Beijing strain family (7, 19). Molecular epidemiological studies also support the notion that the more recently evolved strains (termed typical Beijing strains) are adapted to spread and cause disease, given their frequency of occurrence, in comparison to distantly evolved strains (termed atypical Beijing strains) (14, 18). However, factors including social behavior which could lead to the spread of sublineage 7 strains were not investigated in this study. An alternative explanation for the high frequency of occurrence of the Beijing sublineage 7 strains in South Africa could be a strong founder effect. However, comparison between Beijing genotypes from South Africa and East Asia showed that at least nine different Beijing founder strains were introduced into South Africa (data not shown). The spectrum of founder strains is likely to reflect that in East Asia at the time of importation into South Africa since the East Asia immigrants were derived from many different locations throughout the region. Thus, it is probable that the subsequent spread of sublineage 7 in South Africa represents evolutionary selection of an individual sublineage rather than a founder effect.
The observation of a changing epidemiological picture concurrent with the accumulation of evolutionary events cannot be linked to the sSNPs analyzed as these events are thought to be largely neutral (9). Conversely, mutation in the mismatch repair genes (15), IS6110 integration (2), and chromosomal deletions (19) may influence the strains' ability to spread and cause disease by altering or eliminating gene function. Mutation in the mismatch repair genes has been hypothesized to enhance the ability of the Beijing strains to acquire drug resistance (15). IS6110 insertion is thought to influence the characteristics of the Beijing strains by altering the function of genes by either causing knockouts or disrupting or enhancing promoter activity (2). The accumulation of chromosomal deletions has been suggested to contribute to an altered virulence in the Beijing strain family (19), and in a previous study, the analysis of 100 M. tuberculosis strains suggested an association between deletion and an altered phenotype (20). In our study the most pathogenic strains, as measured by their frequency of occurrence, were associated with the accumulation of deletions, including RD150. The true extent of chromosomal evolution within the Beijing strain family almost certainly extends beyond the evolutionary events tested in this study, and the influence of such events on the phenotype of different strains remains to be determined.
Our study design differs from previous studies as we analyzed intrastrain family phenotypic-genotypic correlations compared to interstrain family phenotypic-genotypic correlations (5, 7, 16, 17, 21). In our study, the analysis of patient demographic and clinical data failed to demonstrate an association with any of the Beijing strain sublineages. This implies that the appearance of secondary cases occurs independently of the proportion of infectious (smear-positive) cases in the different sublineages, thereby suggesting that the characteristics of the strains from the different sublineages could be the key factors influencing the ability to spread and cause disease.
We noted the presence of drug resistance in all of the Beijing sublineages and noted that the frequency of occurrence of the drug resistance phenotype was not associated with the most recently evolved Beijing sublineage. The absence of spread of drug resistance in the most recently evolved Beijing sublineage could be explained by the fitness cost as a result of the acquisition of mutations in the genes conferring drug resistance (1). The high number of resistant cases in sublineage 5 (urban setting) and sublineages 6 and 7 (rural settings) was associated with MDR tuberculosis and clearly indicates that MDR tuberculosis can spread within communities if not managed and treated appropriately.
In summary, we suggest that the pathogenic characteristics of M. tuberculosis strains are not conserved within defined strain families but rather that individual lineages within strain families have evolved unique pathogenic characteristics.
We thank the nurses and data analysts for the collection and processing of clinical data. We thank I. Toms, Department of Health, City of Cape Town, and are indebted to the residents of the epidemiological field site.
Published ahead of print on 14 March 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»