Previous Article | Next Article ![]()
Journal of Clinical Microbiology, November 2002, p. 4072-4076, Vol. 40, No. 11
0095-1137/02/$04.00+0 DOI: 10.1128/JCM.40.11.4072-4076.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Department of Biosciences, Faculty of Natural Science, University of Hertfordshire, Hatfield, Hertfordshire AL10 9AB,1 Molecular Biology Unit, Central Public Health Laboratory, London NW9 5HT, United Kingdom2
Received 23 May 2002/ Returned for modification 31 July 2002/ Accepted 22 August 2002
|
|
|---|
|
|
|---|
M. tuberculosis isolates cannot be resolved by phenotypic strain typing, and their genomes lack heterogeneity. IS6110 restriction fragment length polymorphism typing has been used to distinguish isolates of the M. tuberculosis complex (22). The method is based on the presence of different numbers of copies (varying from 0 to 25) of an insertion sequence, IS6110, in strains of M. tuberculosis and the variability of their positions in the genome. The technique is of limited value, as it cannot easily differentiate between M. tuberculosis isolates that contain only a few or no copies of the insertion sequence, as has been found in countries with a high incidence of tuberculosis, such as India (40% of isolates), and to a lesser degree (8% of isolates) in the United Kingdom (9, 15).
Amplified fragment length polymorphism (AFLP) analysis is a PCR-based technique involving the restriction of bacterial genomic DNA with two restriction enzymes, the ligation of adapters to the restriction fragments, and selective amplification of sets of the restriction fragments with adapter-specific primers. The amplified fragments are then sized by gel electrophoresis (24). The addition of a fluorescent label to one of the PCR primers allows the amplified fragments to be detected with an automated DNA sequencer.
AFLP analysis has been used in the typing of strains of different pathogens (10, 13, 19, 20) and also for plant and animal genetic mapping, medical diagnostics, and phylogenetic studies (18). The whole-genome fingerprinting technique fluorescent AFLP (FAFLP) analysis has also been used to identify polymorphisms in different strains of M. tuberculosis (12). FAFLP offers improved resolution compared with restriction fragment length polymorphism typing of isolates with a single copy of IS6110.
Using FAFLP, Goulding et al. (12) differentiated among 65 clinical M. tuberculosis isolates of known IS6110 type by analyzing DNA restriction fragments from 50 to 500 bp in size. Thirty-eight discriminatory fragments were identified, allowing differentiation of the isolates. Most of the epidemiologically related groups clustered in the same way as with IS6110 profiling, but a group of seven epidemiologically unlinked isolates with identical single-copy IS6110 profiles was split into four clusters by FAFLP analysis (12).
The present study examines whether including DNA fragments of up to 1,000 bp in the analysis increases the resolving power of FAFLP and, by comparing the results with the predicted in silico results obtained from the two sequenced strains of M. tuberculosis (8; http://www.tigr.org/CMR2), investigates whether a G+C-rich genome can be sized and identified accurately.
|
|
|---|
Experimental FAFLP. We used the method described previously by Goulding et al. (12). Briefly, 500 ng of genomic DNA was digested in a total volume of 20 µl consisting of 5 U of MseI (New England Biolabs, Hitchin, Hertfordshire, United Kingdom), 2 µl of 10x MseI buffer, 0.2 µl of 10x bovine serum albumin, and 1.0 µl of DNase-free RNase A (10 µg/µl) for 1 h at 37°C. To this digest was added 5 U (1.0 µl) of EcoRI (Invitrogen, Paisley, United Kingdom), 1.68 µl of 0.5 M Tris HCl (pH 7.6), and 2.1 µl of 0.5 M NaCl (total volume, 25 µl), and the reaction mixtures were incubated for a further hour at 37°C. Slight variations of this method were carried out to optimize the reaction, including increasing the digest times from 1 to 2 and 3 h and adding known PCR-enhancing agents (5% [vol/vol] dimethyl sulfoxide [DMSO], 10% [vol/vol] glycerol, 60 mM tetramethylammonium chloride [TMAC], or 1% [vol/vol] deionized formamide) to the digest or to the PCR mixture. Endonucleases were inactivated at 65°C for 10 min prior to ligation. To the double-digested DNA was added 25 µl of a solution containing 40 U of T4 DNA ligase (New England Biolabs), 10 pmol of EcoRI adapter, 100 pmol of MseI adapter, and 5 µl of 10x T4 ligase buffer. This reaction mixture was incubated at 12°C for 17 h, heated at 65°C for 10 min to inactivate the ligase, and stored at -20°C.
PCRs were performed in 20-µl volumes containing 2 µl of ligated DNA, 0.1 µM 5-carboxy-fluorescein-labeled nonselective primer specific to the EcoRI adapter, 0.25 µM MseI nonselective primer specific to the MseI adapter, 2 µl of 10x Taq polymerase buffer, 200 µM each deoxynucleoside triphosphate, 1.5 mM MgCl2, and 0.5 U of Taq polymerase. Touchdown PCR cycling conditions were used for amplifying the fragments: a 2-min denaturation step at 94°C (one cycle) followed by 30 cycles of denaturation at 94°C for 20 s, a 30-s annealing step (see below), and a 2-min extension step at 72°C. The annealing temperature for the first cycle was 66°C; for the next nine cycles, the temperature was decreased by 1°C at each cycle. The annealing temperature for the remaining 20 cycles was 56°C. This was followed by a final extension at 60°C for 30 min. PCR was performed in a PE-9600 thermocycler (Perkin-Elmer Corp., Norwalk, Conn.).
The amplification products were separated on a 5% denaturing polyacrylamide gel (FMC SingGel) on an ABI 377 DNA sequencer (Perkin-Elmer Corp.). The sample (1 µl) was added to 2.5 µl of loading dye, which was a mixture containing 5 µl of formamide, 1 µl of dextran blue-50 mM EDTA loading solution, and 0.5 µl of the internal lane standard GeneScan 2500 labeled with the fluorophore ROX (PE Biosystems). The sample mixture was heated at 95°C for 2 min, cooled on ice, and immediately loaded onto the gel. The electrophoresis conditions were 2 kV at 51°C for 16 h, using 1x Tris-borate-EDTA as the buffer.
GeneScan and GenoTyper software (PE Biosystems) were used to ascertain the presence or absence of precisely sized fragments.
Reproducibility. To check the reproducibility of the FAFLP analysis, at least two different sets of digestions and ligations were carried out on the DNAs of the clinical isolates. For each digestion-ligation, at least three different PCRs were carried out separately, and each was run on a different gel.
In silico FAFLP. The complete genome sequences of M. tuberculosis strains H37Rv and CDC1551 were analyzed with the Restriction Digest Tool of The Institute for Genomic Research (http://www.tigr.org/tigr-scripts/CMR2/restrict_display.pl) using EcoRI and MseI. Data concerning the sizes and predicted number of EcoRI-MseI fragments were imported into a spreadsheet, and fragment sizes were adjusted to allow for the addition of primer sequences during PCR, enabling direct comparison between the predicted in silico and in vitro results.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Comparison of in silico-predicted number of M. tuberculosis MseI/EcoRI fragments with actual number of fragments detected by FAFLP analysisa
|
Only 5 of the 11 predicted fragments (sizes, 195, 250, 252, 434, and 821 bp) were present in the experimental results of the MseI+TA-EcoRI FAFLP analysis of H37Rv (Table 2). These five fragments were within ±2 bp of their predicted sizes. Four fragments of other sizes were detected that were not predicted in silico (Table 2).
|
View this table: [in a new window] |
TABLE 2. Fragments predicted in silico and generated by FAFLP analysis of M. tuberculosis strain H37Rv using MseI + TA selective primer and EcoRI nonselective primer
|
|
View this table: [in a new window] |
TABLE 3. Sizes of discriminatory fragments generated by FAFLP analysis of 10 IS6110-typed M. tuberculosis strainsa
|
![]() View larger version (39K): [in a new window] |
FIG. 1. FAFLP profiles of M. tuberculosis single-copy IS6110 isolates 157 and 202. The plot shows the fragment sizes in the range from 965 to 1,010 bp. The numbered boxes show the size in base pairs of each fragment. A fragment at 996 bp is present in isolate 202 and absent in isolate 157. The vertical scale measures the efficiency of PCR amplification of each fragment in fluorescent units.
|
![]() View larger version (33K): [in a new window] |
FIG. 2. FAFLP reproducibility study for M. tuberculosis clinical isolate 157. The four plots show the FAFLP profiles for the same isolate across a 40-bp window (372 to 412 bp). The four reactions were carried out separately and run on different gels. The vertical scale measures the efficiency of PCR amplification of each fragment in fluorescent units.
|
|
|
|---|
Although DMSO, glycerol, and TMAC did not improve the method, preliminary results suggested that for FAFLP analysis of G+C-rich genomes, the use of deionized formamide in the digest step reduced the problem of unexpected fragments being produced and/or visualized. The identification and characterization of the four unpredicted restriction fragments generated by FAFLP analysis with the MseI+TA and EcoRI primers, and their locations in the genome, may help to elucidate the differences that occur between the in silico and experimental data.
If the failure to detect expected fragments is not due to incomplete digestion but to cut fragments not being amplified because of secondary structure, the addition of dITP in place of dGTP during the PCR step may improve the FAFLP profile. The effect of the dITP may be to increase the number of fragments seen, but it may also have an effect on their relative peak heights, as there appears to be some variability in the relative peak heights among different M. tuberculosis strains. However, this may be because amplification appears not to be uniform due to other factors, for example, if peaks with higher signal are comprised of multiple fragments of the same or similar sizes. Other potential causes for this lack of parity between in vitro and in silico data might be that there are errors in the sequences of H37Rv and CDC1551 or that passage of different stocks has led to differences between the isolates studied and the isolates sequenced. Mispriming of PCR may also have an effect on successful generation of expected fragments, but touchdown PCR is employed in FAFLP and should minimize this.
Another way of improving the FAFLP method for M. tuberculosis would be the use of different sets of enzymes, for example, XhoI and HhaI, which cut at IS6110 as well as throughout the rest of the genome. This combination should give FAFLP profiles corresponding to the IS6110 elements in a genome, as well as sampling throughout that genome.
The four-primer combination used in the study by Goulding et al. (12) of 65 IS6110-typed clinical isolates generated 38 discriminatory fragments under 500 bp in size. Of these fragments, only nine discriminated among the 10 (of 65) randomly selected isolates when studied as a subgroup. However, when nonselective primers were used, these 10 isolates gave a further 17 discriminatory fragments between 800 and 1,000 bp in size (Table 1). Two of these 17 polymorphic fragments (800 to 1,000 bp) differentiated between isolates 157 and 202, which both contain only a single copy of IS6110 (Fig. 1), making this a useful supplementary technique for M. tuberculosis isolates containing a single copy of IS6110. The profiles of the three isolates known to be epidemiologically linked (251, 252, and 253) were identical apart from one fragment observed in isolate 251 only. The reason for this is unclear, but it may be a reflection of real changes in the genome of isolate 251 compared with 252 and 253. The reason for the disproportionate number of polymorphisms in the higher size range is not understood. It is not due to the lack of reproducibility of amplification of larger-sized fragments.
In conclusion, experimental FAFLP does not exactly match in silico prediction, even when known PCR enhancers, such as TMAC (5) and formamide (17), are added to the digestion and/or PCR. The FAFLP technique is not as accurate for the G+C-rich genome of M. tuberculosis as for E. coli (1) and Campylobacter (10), as it does not strictly conform to the predicted model, but the fragments were precisely sized (±2 bp). The results were reproducible (Fig. 2) but not as clean as for genomes with low G+C content, such as E. coli, which again may be due to the G+C-rich content of the M. tuberculosis genome. Fragments below a certain signal strength were therefore ignored. FAFLP is capable of resolving strains that cannot be differentiated using IS6110 typing due to their low IS6110 copy numbers, as demonstrated by the two polymorphic fragments (898 and 996 bp) that were present in isolate 202 but not in 157 (both of which contain single IS6110 copies). This difference is significant, as up to 40% of M. tuberculosis isolates in certain countries have only one copy of IS6110 (9). Our study also shows that FAFLP has a greater resolving power when fragments between 501 and 1,000 bp are included, as well as those between 100 and 500 bp.
|
|
|---|
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»