G. D. van der Spuy

A Recently Evolved Sublineage of the Mycobacterium tuberculosis Beijing Strain Family Is Associated with an Increased Ability to Spread and Cause Disease

ABSTRACT This study aimed to reconstruct the evolutionary history of Beijing strains of Mycobacterium tuberculosis and to test the hypothesis that evolution has influenced the ability of the Beijing strains within the different Beijing sublineages to spread and cause disease. A PCR-based method was used to analyze the genome structure of 40 different loci in 325 Beijing isolates collected from new and retreatment tuberculosis patients from an urban setting and 270 Beijing isolates collected from high-risk tuberculosis patients from a rural setting in the Western Cape, South Africa. The resulting data were subjected to phylogenetic analysis using the neighbor joining algorithm. Phylogenetic reconstructions were highly congruent with the “gold standard” phylogenetic tree based on synonymous single-nucleotide polymorphisms, thereby allowing a prediction of the order in which the evolutionary events had occurred. A total of seven independently evolving Beijing sublineages were identified. Analysis of epidemiological data in relation to the Beijing sublineage suggested an association between recent evolutionary change and frequency of occurrence in an urban population (P < 0.001) as well as in the rural population (P < 0.001). This concept was further supported by an association between more recently evolved Beijing strains and an increased ability to transmit and to cause disease (odds ratio, 5.82; 95% confidence interval, 3.13 to 10.82 [P < 0.001]). An association between Beijing sublineage and demographic and clinical parameters and drug resistance could not be demonstrated. From these data, we suggest that the pathogenic characteristics of Beijing strains are not conserved but rather that strains within individual lineages have evolved unique pathogenic characteristics.

Microevolution of the direct repeat region of Mycobacterium tuberculosis: Implications for interpretation of spoligotyping data

ABSTRACT The direct repeat (DR) region has been determined to be an important chromosomal domain for studying the evolution of Mycobacterium tuberculosis. Despite this, very little is known about microevolutionary events associated with clonal expansion and how such events influence the interpretation of both restriction fragment length polymorphism (RFLP) and spoligotype data. This study examined the structure of the DR region in three independently evolving lineages of M. tuberculosis with a combination of DR-RFLP, spoligotyping, and partial DNA sequencing. The results show that the duplication of direct variable repeat (DVR) sequences and single-nucleotide polymorphisms is rare; conversely, the deletion of DVR sequences and IS6110-mediated mutation is observed frequently. Deletion of either single or contiguous DVR sequences was observed. The deletion of adjacent DVR sequences occurred in a dependent manner rather than as an accumulation of independent events. Insertion of IS6110 into either the direct repeat or spacer sequences influenced the spoligotype pattern, resulting in apparent deletion of DVR sequences. Homologous recombination between adjacent IS6110 elements led to extensive deletion in the DR region, again demonstrating a dependent evolutionary mechanism. Different isolates from the same strain family and isolates from different strain families were observed to converge to the same spoligotype pattern. In conclusion, the binary data of the spoligotype are unable to provide sufficient information to accurately establish genotypic relationships between certain clinical isolates of M. tuberculosis. This has important implications for molecular epidemiologic strain tracking and for the application of spoligotype data to phylogenetic analysis of M. tuberculosis isolates.

Spoligotype Signatures in the Mycobacterium tuberculosis Complex

ABSTRACT Evolution of the direct repeat region in Mycobacterium tuberculosis has created unique spoligotype signatures specifically associated with IS6110-defined strain families. Spoligotyping signatures may enable the analysis of the strain population structure in different settings and will enable the rapid identification of strain families that acquire drug resistance or escape protective immunity in drug and vaccine trials.

Changing Mycobacterium tuberculosis population highlights clade-specific pathogenic characteristics.

Mycobacterium tuberculosis strains can be classified into a number of major clades according to defined evolutionary markers. It is hypothesised that strains comprising these clades have evolved different properties which may influence a local strain population structure. To investigate this, we analysed the incidence of tuberculosis caused by the predominant clades (Beijing, Haarlem, LAM, Quebec and the Low-Copy Clade) found in a community within the Cape Town metropole in South Africa over a 12-year period. We found that while the incidence of cases infected with strains of the Haarlem, LAM, Quebec and the Low-Copy Clades remained relatively stable, that of cases of the Beijing clade increased exponentially over time, with a doubling time of 4.86 years (P=0.018). This growth was exclusively attributable to drug-susceptible strains. Although drug-resistant Beijing cases remained constant in number, non-Beijing drug-resistant cases declined over time (P=0.007). Drug-susceptible Beijing-infected cases had a greater proportion of smear-positive sputa than their non-Beijing counterparts (P=0.013) and were less likely to be successfully treated (retreatment cases) (P=0.026). Recent evidence suggests that these differences likely reflect enhanced pathogenicity rather than transmissibility. The rapid emergence of Beijing strains demonstrates adaptation to conditions within the study community and poses a grave challenge to future TB control.

Discordance between Mycobacterial Interspersed Repetitive-Unit-Variable-Number Tandem-Repeat Typing and IS6110 Restriction Fragment Length Polymorphism Genotyping for Analysis of Mycobacterium tuberculosis Beijing Strains in a Setting of High Incidence of Tuberculosis

ABSTRACT IS6110 restriction fragment length polymorphism (RFLP) genotyping is the most widely used genotyping method to study the epidemiology of Mycobacterium tuberculosis. However, due to the complexity of the IS6110 RFLP genotyping technique, and the interpretation of RFLP data, mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) genotyping has been proposed as the new genotyping standard. This study aimed to determine the discriminatory power of different MIRU-VNTR locus combinations relative to IS6110 RFLP genotyping, using a collection of Beijing genotype M. tuberculosis strains with a well-established phylogenetic history. Clustering, diversity index, clustering concordance, concordance among unique genotypes, and divergent and convergent evolution were calculated for seven combinations of 27 different MIRU-VNTR loci and compared to IS6110 RFLP results. Our results confirmed previous findings that MIRU-VNTR genotyping can be used to estimate the extent of recent or ongoing transmission. However, molecular epidemiological linking of cases varied significantly depending on the genotyping method used. We conclude that IS6110 RFLP and MIRU-VNTR loci evolve independently and at different rates, which leads to discordance between transmission chains predicted by the respective genotyping methods. Concordance between the two genotyping methods could be improved by the inclusion of genetic distance (GD) into the clustering formulae for some of the MIRU-VNTR loci combinations. In summary, our findings differ from previous reports, which may be explained by the fact that in settings of low tuberculosis incidence, the genetic distance between epidemiologically unrelated isolates was sufficient to define a strain using either marker, whereas in settings of high incidence, continuous evolution and persistence of strains revealed the weaknesses inherent to these markers.

Evidence that the Spread of Mycobacterium tuberculosis Strains with the Beijing Genotype Is Human Population Dependent

ABSTRACT This study describes a comparative analysis of the Beijing mycobacterial interspersed repetitive unit types of Mycobacterium tuberculosis isolates from Cape Town, South Africa, and East Asia. The results show a significant association between the frequency of occurrence of strains from defined Beijing sublineages and the human population from whom they were cultured (P < 0.0001).

Genotypic and Phenotypic Characterization of Drug-Resistant Mycobacterium tuberculosis Isolates from Rural Districts of the Western Cape Province of South Africa

ABSTRACT Genotypic and phenotypic analysis of drug-resistant Mycobacterium tuberculosis isolates from the Western Cape Province of South Africa showed that drug resistance is widespread and recently transmitted. Multidrug-resistant (MDR) isolates comprise 40% of this collection, and a large pool of isoniazid monoresistance may be a future source of MDR tuberculosis.

Calculation of the Stability of the IS6110 Banding Pattern in Patients with Persistent Mycobacterium tuberculosis Disease

ABSTRACT The interpretation of molecular epidemiologic data of Mycobacterium tuberculosis infection is dependent on the understanding of the stability and evolutionary characteristics of the DNA fingerprinting marker used to classify clinical isolates. This study investigated the stability of the IS6110 banding pattern in serial tuberculosis isolates collected from patients resident in an area with a high incidence of tuberculosis. Evolutionary changes were observed in 4% of the strains, and a half-life (t1/2) of 8.74 years was calculated, assuming a constant rate of change over time. This rate may be composed of a high rate of change seen during the early disease phase (t1/2 = 0.57 years) and a low rate of change seen in the late disease phase (t1/2 = 10.69 years). The early rate probably reflects change occurring during active growth prior to therapy, while the low late rate may reflect change occurring during or after treatment. We demonstrate that the calculation of these rates is strongly influenced by the time interval between onset of disease and sputum sampling. These calculations are further complicated by partial replacement of the original strain population, resulting in the sporadic appearance of clonal variants in sputum specimens. Therefore, the true extent of genetic diversity may be underestimated within each host, thereby influencing molecular epidemiological data used to establish transmission chains.

Evolution of the IS6110-Based Restriction Fragment Length Polymorphism Pattern during the Transmission of Mycobacterium tuberculosis

ABSTRACT Interpretation of the molecular epidemiological data of Mycobacterium tuberculosis is dependent on the validity of the assumptions that have been made. It is assumed that the IS6110 banding pattern is sufficiently stable to define epidemiological events representing ongoing transmission. However, molecular epidemiological data also support the observation that the IS6110 banding pattern may change over time. Factors affecting this rate may include the nature and duration of disease in a host and the opportunity to experience different host environments during the transmission cycle. To estimate the rate of IS6110 change occurring during the process of transmission, M. tuberculosis isolates from epidemiologically linked patients were genotypically characterized by restriction fragment length polymorphism (RFLP) analysis. The identification of IS6110 banding pattern changes during ongoing transmission suggested that a rate could be estimated. IS6110 change was significantly associated with strains with >5 IS6110 elements (P = 0.013) and was not observed in low-copy-number isolates. The minimum rate of appearance of variant strains was calculated to be 0.14 variant cases per source-case per year. This data suggest that clustering of isolates based on identical RFLP patterns is expected to underestimate transmission in patients infected with high-copy-number isolates. A model based on the rate of appearance of both variant and invariant strains demonstrates that the genotypically defined population structure may change by 18.6% during the study period of approximately 6.5 years. The implications for the use of RFLP data for epidemiologic study are discussed.

Mapping of IS6110 flanking regions in clinical isolates of Mycobacterium tuberculosis demonstrates genome plasticity.

Southern hybridization was used in combination with IS6110 insertion‐locus‐specific probes in a comparative study to determine the structure of chromosomal domains flanking IS6110 elements in clinical isolates of Mycobacterium tuberculosis. The resulting restriction fragment length polymorphism (RFLP) data demonstrated three mutational mechanisms responsible for the polymorphisms observed: IS6110 insertion, chromosomal mutation and deletion. The frequency of IS6110 insertion within many of the chromosomal regions demonstrates that preferential integration regions are common in M. tuberculosis. Mapping the IS6110 insertion positions and chromosomal deletions in relation to the M. tuberculosis H37Rv and M. bovis BCG genome sequences reveals numerous disruptions of predicted open reading frames (ORFs). A phylogenetic tree, based on the mutational data, showed a number of independently evolving lineages of M. tuberculosis, while analysis of the mutational events occurring at each branch point suggests both divergent and convergent evolution. A significant positive correlation was demonstrated between the mutation rate and the frequency of occurrence of different isolates in families of strains, suggesting that evolution may impact on strain ‘fitness’ or that strain proliferation may increase the chance of mutation. We conclude that the genome of clinical isolates of M. tuberculosis continues to evolve.