|
|
||||||||
1Epidemiology and Molecular Biology Unit, The Children's University Hospital, Temple Street, Dublin 1, Ireland 2Department of Clinical Microbiology, Royal College of Surgeons in Ireland, York Street, Dublin 2, Ireland
Correspondence Désirée E. Bennett Desiree.Bennett{at}tsch.ie
Received February 19, 2003
Accepted June 2, 2003
Invasive disease-associated strains of Neisseria meningitidis were analysed by multilocus restriction typing (MLRT), which involves the restriction fragment-length polymorphism analysis of PCR products generated from the seven loci of housekeeping genes used in MLST. Several different restriction patterns (alleles) were observed for each of the seven loci examined. Greater allelic variation was observed with the fumC and pgm loci than with the abcZ and adk loci, suggesting that the latter were more conserved. The alleles at each of the seven loci were combined to give an allelic profile or restriction type (RT). A good correlation between RT and serogroup, serotype and serosubtype was observed, as all C 2ap1.2,5 strains were contained in a single RT, as were all but one strain of B 4p1.4. In this study, MLRT proved to be an efficient, effective and relatively inexpensive method for N. meningitidis strain characterization.
| INTRODUCTION |
|---|
|
|
|---|
Multilocus enzyme electrophoresis (MLEE) and its recently developed DNA sequence-based counterpart, multilocus sequencing typing (MLST), are the most valuable techniques adapted for the study of the long-term or global epidemiology of N. meningitidis (Maiden et al., 1998; Selander et al., 1986). MLST offers many advantages, including speed and enhanced discrimination, over MLEE and has now almost entirely superseded MLEE as a strain typing method (Enright & Spratt, 1999; Spratt, 1999; Yakubu et al., 1999). The electronic portability of the nucleotide sequence data, enabling rapid global exchange of molecular typing information for epidemiological comparisons, is one of the biggest advantages of MLST over MLEE. However, while the technology for nucleotide sequencing has improved substantially in the last few years, with the advent of automated sequencers, sequencing of multiple housekeeping genes still remains time-consuming (Clarke et al., 2001a; Tondella et al., 1999). In addition, a further factor that must be taken into consideration with MLST is cost and, despite the additional commercial competition on the DNA sequencing market, automated DNA sequencing instruments remain expensive, as do fluorescently labelled dideoxynucleotide sequencing kits (Hookey & Arnold, 2001; Goulding et al., 2000; Olive & Bean, 1999). Initial start-up costs are high, ranging from approximately $100 000 to $300 000 (Diggle & Clarke, 2002), and it is generally accepted that high throughput is required to minimize nucleotide sequencing costs. Therefore, to achieve maximum cost- and time-effectiveness and efficiency, an automated sequencer must be used to maximum capacity every time, which, in most laboratories, may correspond to a single use per week. Even in the recent description of semi-automated MLST performed using a robotic liquid-handling system, it was reported that MLST was probably only cost-effective in national reference laboratories or large clinical laboratories where high throughput could justify the cost (Clarke et al., 2001a). These authors later report that the cost of consumables alone for performing MLST on seven genes (for both strands) arrives at $32 per bacterial isolate (Diggle & Clarke, 2002). Furthermore, manual DNA sequencing, even with an automated DNA sequencer, requires a high level of technical competence in carrying out the DNA sequencing reactions and interpreting the results and also in instrument maintenance. In short, in a standard laboratory without robotics and personnel dedicated to sequencing, it is too expensive, time-consuming, laborious and tedious to perform complete MLST on seven loci for a large number of isolates, on a routine basis (Clarke et al., 2001a; Olive & Bean, 1999; Shlush et al., 2002). A more user-friendly, relatively cheap, rapid and simple method based on the multilocus approach is required. A technique using the multilocus approach, based on cleavase fragment length polymorphism (CFLP) analysis of metabolic genes, was recently demonstrated to be reproducible, rapid and discriminatory for N. meningitidis epidemiology (Tondella et al., 1999). However, it has the drawback of being as complex and as expensive to perform as DNA sequencing (Olive & Bean, 1999). A further technique with phylogenetic significance, fluorescent amplified-fragment length polymorphism (FAFLP) genotyping, has been applied successfully to identify hypervirulent, hyperendemic lineages of N. meningitidis (Goulding et al., 2000; Hookey & Arnold, 2001). However, FAFLP analysis appears time-consuming and laborious and still requires the use of an automated sequencer. In a recent report, denaturing HPLC (DHPLC) was used as a simple and faster approach than DNA sequencing, for comparing the products obtained by MLST. The authors demonstrated 100 % correlation between MLST in its original format and their DHPLC-MLST method (Shlush et al., 2002).
The goal of the studies described in this paper was to evaluate a simple, rapid and reproducible procedure as an alternative to MLST for multilocus typing of N. meningitidis isolates. An alternative to MLST is multilocus restriction typing (MLRT; Müller-Graf et al., 1999), in which variation at several loci is indexed by restriction fragment-length polymorphism (RFLP) analysis of PCR-amplified genes. During the preparation of the present report, a paper was published documenting the use of MLRT as a novel tool for studying the global epidemiology of Burkholderia cepacia (Coenye & LiPuma, 2002). A further paper by the same researchers (Coenye et al., 2002) demonstrated a good agreement between MLRT and other commonly used genotyping methods for B. cepacia. Here, we describe MLRT analysis of the seven housekeeping genes used in MLST and, as an initial step in the assessment of MLRT, we compared data obtained by MLRT and serology for a large number of N. meningitidis isolates that were collected in Ireland.
| METHODS |
|---|
|
|
|---|
Gene targets, PCR primers and cycling conditions.
The seven gene targets (abcZ, adk, aroE, fumC, gdh, pdhC and pgm) and PCR primer pairs used in this study were those described previously for MLST analysis (Maiden et al., 1998). Touchdown PCR cycling conditions were used with annealing temperatures of 59 °C for the initial three cycles, 57 °C for the next three cycles followed by 52 °C for the final 30 cycles on a PTC 200 thermal cycler (MJ Research). These cycling conditions are slightly different from those described for MLST analysis by Clarke et al. (2001a). PCRs were performed in a final volume of 50 µl consisting of 1 µl extracted DNA as template, 5 µl 10x PCR buffer, 1.25 U Taq DNA polymerase (Gibco BRL) with 50 pmol of each primer and 0.2 mM of each of dATP, dCTP, TTP and dGTP. A MgCl2 concentration of 2.5 mM was used for all primer pairs except aroE primers, where 2.0 mM MgCl2 was used.
Restriction endonuclease typing.
Following an initial study with a panel of frequent-cutting restriction endonucleases (4 bp recognition sequence) and a serogroup B and C strain, a single restriction endonuclease that yielded at least three fragments for each locus and strain was chosen. MspI was used to analyse PCR products of abcZ, adk, fumC and gdh. MnlI was used to analyse the PCR products of aroE, pdhC and pgm. Restriction patterns were assessed on the basis of the number and size of restriction fragments obtained following digestion and electrophoretic separation (4 V cm-1) in a horizontal 3.5 % (w/v) agarose gel containing 0.5 µg ethidium bromide ml-1. The restriction patterns obtained with each locus were analysed using the band detection and molecular mass determination features of the 1D Advanced software (version 4.01.2) and compared using the 1D Database software (version 1.12) available from Phoretix International to remove any subjectivity or bias in their analysis. Different allele numbers were assigned arbitrarily to patterns that differed in either the number and/or size of any of the bands. As with MLST, the combination of alleles at each of the seven loci gave an allelic profile or restriction type (RT) and strains were said to be different if the allele at any of the seven loci differed.
Using the RT data, specific clonal complexes within the 83 strains were investigated with the BURST algorithm within the START software package, version 1.0.0 (K. Jolley, University of Oxford, UK) (Jolley & Urwin, 2001). In addition, a dendrogram was constructed using the percentage of alleles at which strains differed as the distance between them, as calculated using the unweighted pair group method with arithmetic mean (UPGMA) algorithm within the START software package.
| RESULTS |
|---|
|
|
|---|
|
|
Following analysis of the seven-loci RTs with the BURST program, the strains were grouped into two lineages, with a single singleton RT. Group 1 contained 26 RTs and 79 strains (95.2 %) and group 2 contained only two RTs or three strains (RTs 21 and 24). RT7 was not represented in either of the groups, indicating that this RT varied by three or more loci from every other RT.
MLRT association with serology
The 83 strains examined included representatives of all the main serogroups of N. meningitidis, with 38 serogroup B strains, 30 C strains, four each of W135 and NG, two each of Y and X and one each of A, 29E and Z (Fig. 2). Six RTs contained all the serogroup C strains examined. RTs 1 and 6 exclusively contained 25 (83.3 %) of these strains, all with either subtype 2a, NT or 2b and with serosubtype p1.2, p1.5, p1.2,5 or NT. The serogroup B strains examined yielded 14 different RTs. Eight RTs contained a single strain each and a further three RTs contained two strains each. A single RT, RT5, contained 14 serogroup B strains and, in conjunction with two other RTs, accounted for 24 (63.5 %) of all serogroup B strains. RT15 contained strains of serogroups B and C, two strains of C 4p1.10 and three strains each of B 15p1.7,16 and B 4p1.15. This was the only RT to contain strains of both serogroups. Among the four serogroup W135 strains analysed, three RT combinations were obtained, which differed from each other by only either the fumC or pdhC locus. All of these were designated serotype NT and were either serosubtype p1.3,6 or NT. In addition, both serogroup X strains examined yielded the same RT, RT21, but both were of the same serology combination, X 21p1.16. These strains, along with the ATCC 13101 reference strain, formed the basis of the second lineage as defined by BURST analysis. The two serogroup Y strains examined yielded different RTs, differing at four loci, but were of completely different subtype and serosubtype combinations. Each of the strains of serogroups A, 29E and Z yielded unique RTs; the serogroup A strain was RT25, the serogroup 29E strain was RT7 and the serogroup Z strain was RT27. The four NG strains examined yielded three different RTs, RT5, 8 and 9. The serogroup 29E strain yielding RT7 constituted the only singleton as defined by BURST analysis. Strains of all other serogroups, A, B, C (except ATCC 13101), W135, Y, Z and NG, yielded RTs that were all deemed to be part of the same lineage as defined by BURST analysis.
In general, there was a good correlation between RT and serogroup, serotype and serosubtype, apart from single strain exceptions. Strains of RT1 (n=21) were all serogroup C; 19 were serotype 2a and the serotypes of the remaining two were designated non-typable (NT) by the conventional serology typing scheme. Of these 21 strains, 10 were serosubtype p1.2,5, six were designated NT and the remaining five serosubtypes were p1.2, p1.5 (each n=2) and p1.10. All four C 2bp1.2,5 strains yielded RT6 and all but one of the 11 B 4p1.4 strains examined yielded RT5. The remaining B 4p1.4 strain yielded a unique RT, RT29, differing in the fumC and pgm alleles. However, RT5 also contained single B 4p1.7, B 1p1.4, B 4p1.3,6 and B 15p1.12,13 strains and two NG 4p1.4 strains. A further grouping of B NTp1.15 was also observed, with three of the five strains examined yielding the same RT, RT16. The other two B NTp1.15 strains were recovered from the same patient and differed at only one allele, the pgm locus, from RT16, producing RT17 (Fig. 2). RT15 contained a variety of strains with three distinct serogroup, subtype and serosubtype combinations, three strains of B 15p1.7,16, three strains of B 4p1.15 and two strains of C 4p1.10.
Of the 12 strains recovered from individuals without invasive meningococcal disease, seven were contained in RTs along with strains having similar serology combinations, whereas the other five, each with unique serology combinations, yielded distinct RTs. It is difficult to make any further conclusions regarding congruence between RT and serological typing, as only single strains of other serology combinations or strains, whose specific serology combinations were unknown, as they were designated non-typable and non-serosubtypable by conventional serology, were examined. Among the four NG strains examined, three different RTs were obtained. The two NG 4p1.4 strains were contained in RT5 with strains of serogroup B possessing the same serotypes and serosubtypes. The NG NTp1.16 strain yielded a unique RT, RT8, and the NG NTp1.5 strains had the same RT, RT9, as a serogroup Y NTp1.5 strain.
| DISCUSSION |
|---|
|
|
|---|
Two clonal complexes and one singleton RT were identified following BURST analysis. The primary clonal complex contained all but three RTs, representing 95.2 % of the strains examined and including all the serogroup A, B, C (except ATCC 13101), W135, Y, NG and Z strains analysed. The second group contained the two serogroup X strains and strain ATCC 13101, and the singleton RT was the serogroup 29E strain. To be included in a clonal complex/group, every strain must share at least five identical alleles with at least one other strain in that group. In the primary complex in this study, all 79 strains fulfil this criterion, but the five shared alleles are not the same five alleles in all of the 79 strains. Therefore, it is unlikely that these lineage delineations represent the true phylogenetic relationship among these strains.
An excellent association between RTs and serology was observed, apart from single strain exceptions. Strains of the same RT were mainly of the same serogroup and serotype combination, i.e. strains with serogroup C, serotype 2b and serosubtype p1.2,5 were all contained within RT6, and strains with RT1 were predominantly serogroup C, serotype 2a and serosubtype p1.2,5. However, two strains whose serotypes were non-typable and six strains whose serosubtypes were non-typable were also contained within RT1. It is possible that these strains were 2a serotype and p1.2 or p1.5 or both serosubtypes but did not express the genes. Similarly, strains with RT5 were mainly B 4p1.4, although a B 4p1.4 strain also formed a unique RT, RT29, and RT5 also contained single B 4p1.7, B 4p1.3,6, B 15p1.12,13 and B 1p1.4 strains. However, this is not surprising, as it has been documented previously that B 4p1.4 strains can be closely related genetically to strains of other diverse serological phenotypes (Van Looveren et al., 1998, 2001). Therefore, it could be concluded that serotype 4 could be important for inclusion in RT5, although other strains of serotype 4, three B 4p1.15 strains and two C 4p1.10 strains, yielded a different RT, RT15. RT15 also contained three strains with a B 15p1.7,16 phenotype. RT5 also contained two strains whose serogroup was nongroupable but whose serotypes were 4p1.4; this suggests that these strains were probably serogroup B strains that did not express the B capsule protein. Furthermore, another NG strain, NG NTp1.5, yielded the same RT as a serogroup Y strain with a similar serology combination, Y NTp1.5. This could suggest that the NG strain was, in fact, a non-expressing serogroup Y strain.
In almost all cases, RTs contained strains of only one serogroup (apart from NG strains). Two RTs, RT1 and RT6, exclusively contained 83.3 % of serogroup C strains examined and were either subtype 2a, NT or 2b with serosubtype p1.2, p1.5, NT or p1.2,5. These serotypes and serosubtypes represent 98 % of the types and subtypes of the serogroup C isolates recovered in Ireland during 2000 and approximately 88 % of serogroup C isolates recovered during 1999. Serogroup C strains were contained in only four other RTs. This suggests that the prevalent serogroup C strains associated with invasive meningococcal disease in Ireland are a genetically homogeneous population. In contrast, 14 RTs contained all the serogroup B strains examined. RTs 5, 15 and 16 contained 63 % of all serogroup B strains examined. These strains were predominantly subtypes 4, NT and 15 and serosubtypes p1.4, p1.7, p1.15, NT and p1.7,16. These serotype and serosubtype combinations account for approximately 65 % of all serogroup B isolates recovered during 1999 and 2000 in Ireland. RT15 contained the Norwegian serogroup B strain, H44/76, which has previously been designated as belonging to the clonal complex ST-32 complex/ ET-5 complex by MLST/MLEE characterization (http://www.mlst.net).
In this study, there was no evidence of strains with the same genetic background possessing different serogroups, suggestive of capsular switching, although there was evidence of poor or absent capsule expression. However, any serogroup B isolates yielding RT1 or RT6 or serogroup C isolates yielding RT5 could be indicative of capsular switching. The phenomenon of capsule switching has been described by a number of groups (Swartley et al., 1997; Vogel et al., 2000) and, in light of the recent introduction of MenC vaccination in Ireland (October 2000), it is important to monitor strains for capsular switching. The significance of RT15 containing both serogroup B and C strains is unknown, except that strains of both serogroups have been previously documented to belong to the clonal complex ST-32 complex/ET-5 complex, based on a search of the MLST database (http://www.mlst.net), although none with the C 4p1.10 serotype and serosubtype combination was documented.
MLRT is an application of proven concepts and a variation on the methods of MLST and is therefore suitable for long-term or global epidemiological purposes. Many of the advantages associated with MLST are also valid for MLRT, including 100 % typability, enhanced discrimination achievable through using a multilocus approach and potential for typing organisms directly from clinical material (Clarke et al., 2001b; Kriz et al., 2002). Other advantages of MLRT are simplicity, speed and relative cost and, from this limited study, it can be seen that MLRT also affords a high degree of discrimination, yielding an excellent congruence with serology.
The power of discrimination possible by MLST is undoubtedly much greater than could be expected from MLRT, even though larger gene fragments are examined in MLRT than are sequenced in MLST. Additional major advantages of MLST over MLRT include the portability of sequence data versus restriction patterns and the unambiguity of allele identity in cases where fragments with the same electrophoretic mobility may have completely different sequences. Nevertheless, the benefits of MLRT, combining high discriminatory power and acceptable turnaround time and lower cost compared with complete MLST analysis (Hookey & Arnold, 2001; Olive & Bean, 1999; Shlush et al., 2002), outweigh these advantages, especially for an in-house initial screening method. MLRT is easy to perform, without the need for a high level of technical expertise or the use of expensive specialized equipment; therefore, it could be performed on relatively large numbers of isolates in a routine molecular laboratory. However, as the number of isolates to be examined increases, potentially yielding several distinct DNA banding patterns that may be difficult to compare visually, computer-assisted image acquisition and DNA pattern analysis software may be required, at least to eliminate user bias and subjectivity in restriction pattern classification. An additional benefit of MLRT is its potential for adaptability to MLST, as the previously generated PCR product can be used as a sequencing template. It is these advantages that point towards the use of MLRT as an initial screening method to distinguish different isolates at a basal level without the need for sequencing. Isolates that are indistinguishable by MLRT can then be typed further by MLST without the need to repeat the PCR.
To conclude, in this study, MLRT has proved to be a discriminatory, reproducible, easy-to-perform, rapid and relatively inexpensive typing method to characterize N. meningitidis isolates genetically in a way that reflects their serological classification. However, the reliance of MLRT on comparisons of DNA patterns indicates its use either as an in-house alternative to MLST and/or as an initial screening method to distinguish different isolates prior to MLST analysis. We expect that MLRT will be of value during epidemiological investigations of potential meningococcal outbreak situations, especially in laboratories that do not have direct access to robotics, automated sequencing facilities and specifically trained personnel.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. H. Dyet and D. R. Martin Sequence Variation in the porB Gene from B:P1.4 Meningococci Causing New Zealand's Epidemic J. Clin. Microbiol., February 1, 2005; 43(2): 838 - 842. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. H. Dyet, R. S. Simmonds, and D. R. Martin Multilocus Restriction Typing Method To Predict the Sequence Type of Meningococci J. Clin. Microbiol., April 1, 2004; 42(4): 1742 - 1745. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | J MED MICROBIOL | MICROBIOLOGY | J GEN VIROL | ALL SGM JOURNALS |