Reliable identification of mycobacterial species by PCR-restriction enzyme analysis (PRA)-hsp65 in a reference laboratory and elaboration of a sequence-based extended algorithm of PRA-hsp65 patterns

  • Erica Chimara1,

    Affiliated with

    • Lucilaine Ferrazoli1,

      Affiliated with

      • Suely Yoko Misuka Ueky1,

        Affiliated with

        • Maria Conceição Martins1,

          Affiliated with

          • Alan Mitchel Durham2,

            Affiliated with

            • Robert D Arbeit3 and

              Affiliated with

              • Sylvia Cardoso Leão4Email author

                Affiliated with

                BMC Microbiology20088:48

                DOI: 10.1186/1471-2180-8-48

                Received: 17 October 2007

                Accepted: 20 March 2008

                Published: 20 March 2008

                Abstract

                Background

                Identification of nontuberculous mycobacteria (NTM) based on phenotypic tests is time-consuming, labor-intensive, expensive and often provides erroneous or inconclusive results. In the molecular method referred to as PRA-hsp65, a fragment of the hsp65 gene is amplified by PCR and then analyzed by restriction digest; this rapid approach offers the promise of accurate, cost-effective species identification. The aim of this study was to determine whether species identification of NTM using PRA-hsp65 is sufficiently reliable to serve as the routine methodology in a reference laboratory.

                Results

                A total of 434 NTM isolates were obtained from 5019 cultures submitted to the Institute Adolpho Lutz, Sao Paulo Brazil, between January 2000 and January 2001. Species identification was performed for all isolates using conventional phenotypic methods and PRA-hsp65. For isolates for which these methods gave discordant results, definitive species identification was obtained by sequencing a 441 bp fragment of hsp65. Phenotypic evaluation and PRA-hsp65 were concordant for 321 (74%) isolates. These assignments were presumed to be correct. For the remaining 113 discordant isolates, definitive identification was based on sequencing a 441 bp fragment of hsp65. PRA-hsp65 identified 30 isolates with hsp65 alleles representing 13 previously unreported PRA-hsp65 patterns. Overall, species identification by PRA-hsp65 was significantly more accurate than by phenotype methods (392 (90.3%) vs. 338 (77.9%), respectively; p < .0001, Fisher's test). Among the 333 isolates representing the most common pathogenic species, PRA-hsp65 provided an incorrect result for only 1.2%.

                Conclusion

                PRA-hsp65 is a rapid and highly reliable method and deserves consideration by any clinical microbiology laboratory charged with performing species identification of NTM.

                Background

                The genus Mycobacterium comprises organisms that are heterogeneous in terms of metabolism, growth, environmental niche, epidemiology, pathogenicity, geographic distribution and disease association [1]. While there are notable pathogens such as Mycobacterium tuberculosis, Mycobacterium bovis and Mycobacterium leprae, most are environmental organisms typically acting as opportunistic pathogens. These species, often collectively called nontuberculous mycobacteria (NTM), have been associated with a variety of problems including pulmonary, lymph node, skin, soft tissue, skeletal, and disseminated infections as well as nosocomial outbreaks related to inadequate disinfection/sterilization of medical devices [2]. In recent years, infections due to the subset of rapidly growing NTM, including Mycobacterium fortuitum, Mycobacterium chelonae and Mycobacterium abscessus, have been reported as complications of numerous surgical procedures, particularly involving foreign bodies (e.g., augmentation mammaplasty), high risk sites (e.g., eye) and injections of natural products used as alternative medicines [38].

                In most laboratories, identification of mycobacterial species is based on in vitro growth and metabolic activities. Such phenotypic tests are labor-intensive and time-consuming to perform and may take several days to weeks to complete. Further, for many NTM species, the tests may be poorly reproducible [9], and consequently, the identifications may be ambiguous or erroneous [10].

                DNA-based methods offer the promise of rapid and accurate species identification. However, commercially available DNA probes are available only for a handful of mycobacterial species; moreover, reagents are quite costly. Nucleotide sequence analyses can be used to resolve essentially any bacterial species, but requires both amplification and sequencing.

                Telenti and coworkers described a DNA-based method for species identification of mycobacteria in which a portion of hsp65, the gene encoding the 65 kDa heat shock protein, was amplified by PCR and then analyzed by restriction digest [11]. This approach, referred to as PRA-hsp65, required only routine PCR and agarose gel electrophoresis equipment and could be completed within a few hours. The different species of mycobacteria yielded distinctly different patterns of restriction fragments and thus the species of an unknown isolate could be determined by comparing the fragments observed with published analyses of clinical isolates [1117] and of newly described species [4, 1824]. The availability of an on-line internet resource facilitates the process [25].

                Some studies have observed limitations to PRA-hsp65 which could, potentially, render the approach impractical for routine use. First, within commonly encountered species of clinical significance, such as Mycobacterium avium and Mycobacterium kansasii, as many as six distinct PRA-hsp65 patterns have been encountered [20, 2628]. Such variability could result in a high frequency of ambiguous or uninterpretable patterns. Second, validated protocols for electrophoresis and internal standards have not been defined [17, 29]. Lastly, published tables present patterns which differ within a range of 5–15 bp and lack patterns for recently described species [11, 14, 16]. The aim of this study was to determine whether PRA-hsp65 of mycobacterial isolates provides sufficiently reliable species identification to enable it to be used as the routine methodology in a reference laboratory.

                Results

                Species identification by phenotype and PRA-hsp65 considered separately

                Among the 434 isolates studied, biochemical and phenotypic evaluation alone assigned 371 (85.5%) isolates a species or complex; PRA-hsp65 assigned 404 (93%) isolates a species. Inconclusive results were obtained for 63 (14.5%) isolates by conventional methods compared with 30 (6.9%) isolates using the rapid DNA-based approach; these included nine isolates that could not be identified by either method.

                Species identification by phenotype and PRA-hsp65 compared to sequencing

                For 321 (74.0%) of the 434 isolates both methods gave the same species identification, i.e., the results were concordant (Table 1). Based on prior experience by the authors and others [26, 30], these identifications were presumed to be correct. The hsp65 genes of the remaining 113 (26.0%) isolates giving discordant or inconclusive results were sequenced. Among these, phenotypic testing had assigned 50 isolates to a species or a complex, but sequencing indicated that 33 (66%) of these assignments were incorrect (Table 2). For 63 isolates the phenotypic results were ambiguous and provided only a broad Runyon classification. Even among these, 19 (30.2%) were misclassified compared to conventional expectations [9, 31], including 12 with regard to rate of growth (i.e., slow vs. rapid) and 7 with regard to chromogen production (Table 2). Overall, phenotypic species identification was correct for only 17 (15%) of 113 isolates for which hsp65 sequencing was performed.
                Table 1

                Species identification of 321 isolates which had concordant results by both phenotypic and PRA-hsp65 methods.

                Phenotypic identification

                PRA-hsp65 b

                N (%)

                M. avium complex (146)a

                M. avium 1

                107 (33.5)

                 

                M. avium 2

                24 (7.5)

                 

                M. avium 3

                1 (0.3)

                 

                M. intracellulare 1

                13 (4.1)

                 

                M. intracellulare 4

                1 (0.3)

                M. kansasii (95)

                M. kansasii 1

                95 (29.7)

                M. gordonae (30)

                M. gordonae 1

                2 (0.6)

                 

                M. gordonae 3

                19 (6.0)

                 

                M. gordonae 4

                2 (0.6)

                 

                M. gordonae 5

                1 (0.3)

                 

                M. gordonae 7

                3 (0.9)

                 

                M. gordonae 8

                3 (0.9)

                M. fortuitum complex (24)

                M. fortuitum 1

                21 (6.6)

                 

                M. peregrinum 2

                1 (0.3)

                 

                M. peregrinum 3

                2 (0.6)

                M. chelonae complex (21)

                M. chelonae 1

                5 (1.6)

                 

                M. abscessus 1

                14 (4.4)

                 

                M. abscessus 2

                2 (0.6)

                M. marinum (2)

                M. marinum 1

                2 (0.6)

                M. terrae complex (2)

                M. terrae 1

                1 (0.3)

                 

                M. nonchromogenicum 2

                1 (0.3)

                M. szulgai (1)

                M. szulgai 1

                1 (0.3)

                a Number of isolates.

                b PRA-hsp65 designation; see text for details.

                Table 2

                Results for 96 NTM isolates for which phenotypic methods gave incorrect species identification as determined by hsp65 sequencing.

                Species

                Nb

                Phenotypic result

                M. abscessus (1)c

                1

                SGN

                M. arupense (5)

                1

                M. chelonae complex

                 

                2

                SGN

                 

                1

                SGS

                M. asiaticum (3)

                2

                M. avium complex

                 

                1

                M. gordonae

                M. avium (18)

                1

                M. chelonae complex

                 

                1

                M. fortuitum

                 

                2

                M. kansasii

                 

                10

                SGN

                 

                1

                RGN

                 

                1

                SGP

                 

                1

                SGS

                M. celatum (2)

                1

                M. xenopi

                 

                1

                SGN

                M. chelonae (2)

                2

                SGN

                M. cosmeticum (1)

                1

                M. chelonae

                M. farcinogenes (1)

                1

                M. chelonae complex

                M. flavescens (1)

                1

                RGS

                M. fortuitum (6)

                2

                M. chelonae complex

                 

                1

                RGN

                 

                1

                SGN

                M. genavense (1)

                1

                SGN

                M. gordonae (26)

                1

                RGP

                 

                1

                SGN

                 

                12

                SGS

                M. hassiacum (1)

                1

                RGS

                M. intracellulare (9)

                1

                M. chelonae complex,

                 

                1

                M. gordonae

                 

                7

                SGN

                M. kansasii (7)

                1

                M. nonchromogenicum

                 

                2

                RGP

                 

                1

                SGN

                 

                2

                SGP

                 

                1

                SGS

                M. lentiflavum (3)

                2

                M. avium complex

                 

                1

                M. gordonae

                M. mageritense (1)

                1

                M. fortuitum

                M. marinum (1)

                1

                M. kansasii

                M. mucogenicum (8)

                2

                M. chelonae complex

                 

                1

                M. fortuitum complex

                 

                1

                M. peregrinum

                 

                3

                SGN

                 

                1

                SGS

                M. nebraskense (1)

                1

                M. gordonae

                M. nonchromogenicum (2)

                1

                SGN

                M. peregrinum (4)

                4

                M. chelonae complex

                M. phlei (1)

                1

                RGS

                M. scrofulaceum (3)

                1

                M. avium complex

                 

                2

                SGN

                M. sherrisii (3)

                2

                M. avium complex

                 

                1

                SGN

                M. szulgai (1)

                1

                SGS

                M. terrae (2)

                1

                SGN

                a Species identification was determined by hsp65 sequencing for 113 isolates that had discordant results by PRA-hsp65 and phenotypic studies. For 17 isolates sequencing confirmed the species identification obtained by phenotypic methods.

                b Number of isolates for which the phenotypic identification shown was incorrect.

                c Total number of isolates of that species sequenced. SGS: slowly growing scotochromogen; SGN: slowly growing nonchromogen; SGP: slowly growing photochromogen; RGS: rapidly growing scotochromogen; RGN: rapidly growing nonchromogen; RGP: rapidly growing photochromogen.

                Among the 113 isolates with discordant or inconclusive results, PRA-hsp65 assigned 83 isolates to a species; 71 (85.5%) of these assignments were confirmed by hsp65 partial gene sequencing (Table 3). For most of the remaining isolates, the identifications resolved by PRA-hsp65 and sequencing were consistent with close evolutionary relationships (e.g., M. kansasii and Mycobacterium gastri, Mycobacterium intracellulare and M. avium) (Table 3).
                Table 3

                Results for 12 NTM isolates for which PRA-hsp65 gave incorrect species identification as determined by hsp65 sequencing.

                Species

                Nb

                PRA-hsp65 result

                M. avium (18)c

                1

                M. kansasii 1

                M. farcinogenes (1)

                1

                M. scrofulaceum 1

                M. intracellulare (9)

                1

                M. avium 3

                M. kansasii (7)

                1

                M. avium 2

                 

                1

                M. gastri 1

                M. mucogenicum (8)

                1

                M. chitae 1

                 

                1

                M. gordonae 1

                 

                1

                M. nonchromogenicum 1

                M. nebraskense (1)

                1

                M. avium 3

                M. scrofulaceum (3)

                2

                M. lentiflavum 3

                 

                1

                M. simiae 1

                a Species identification was determined by hsp65 sequencing for 113 isolates that had discordant results by PRA-hsp65 and phenotypic studies. For 71 isolates sequencing confirmed the species identification obtained by PRA-hsp65. For an additional 30 isolates, the PRA-hsp65 patterns obtained were previously unreported (see Table 4).

                b N, number of isolates for which the PRA-hsp65 identification shown was incorrect.

                c Total number of isolates of that species sequenced.

                There were 30 isolates representing 13 PRA-hsp65 patterns not in the available databases and the species was resolved by sequencing. The observed BstEII and HaeIII fragments for these new patterns (designated NP), the source of these isolates and the species identification based on sequencing are listed in Table 4; the observed phenotypes, including antimicrobial susceptibilities, are presented in Table 5. In four instances (NP1, NP11, NP14 and NP17, representing Mycobacterium gordonae Mycobacterium terrae, Mycobacterium sherrisii and Mycobacterium arupense, respectively) multiple isolates with the pattern were identified.
                Table 4

                BstEII and HaeIII fragment lengths (base pairs) for 30 isolates with new patterns by PRA-hsp65.

                Speciesa

                PRA-hsp65

                N

                Fragment BstEII

                Length (bp) HaeIII

                M. arupense

                NP17

                5

                320-115

                145-75-60

                M. avium

                NP10

                1

                320-115

                140-90-60

                M. cosmeticum

                NP6

                1

                320-115

                150-95-80

                M. fortuitum

                NP12

                1

                235-120-85

                140-120-100-55

                M. fortuitum

                NP19

                1

                235-120-100

                145-140-100-55

                M. gordonae

                NP1

                11

                235-120-100

                130-110-95

                M. gordonae

                NP3

                1

                320-130

                130-60

                M. gordonae

                NP13

                1

                235-120-85

                130-90

                M. gordonae

                NP22

                1

                235-130-85

                160-90-60

                M. mageritense

                NP5

                1

                240-130-85

                145-100-50

                M. nonchromogenicum

                NP4

                1

                235-120-85

                145-80-60

                M. sherrisii

                NP14

                3

                235-120-85

                145-130

                M. terrae

                NP11

                2

                235-210

                140-115-70

                a Species identification based on sequencing of hsp65 gene. Bold indicates sequences submitted to GenBank and patterns included in the updated PRA-hsp65 algorithm (see Figures 1, 2 and 3). GenBank accession numbers: NP1, EF601222; NP11, EF601223; NP14, AY365190 [23]; NP17, DQ168662 [18].

                All isolates with new PRA-hsp65 profiles were cultured from sputum, with the following exceptions: NP1: urine (2), feces, liver biopsy and unknown (one each); NP17: unknown (2).

                Table 5

                Phenotypic characteristics of isolates demonstrating previously unreported PRA-hsp65 patterns.

                Species

                PRA hsp65

                25°C

                37°C

                45°C

                pg

                TCH

                nit

                Tween

                NaCl

                Aryl3

                Aryl15

                ag

                pic

                β-gal

                LJ

                HA

                PNB

                INH

                RF

                EMB

                CIP

                OFL

                M. arupense

                NP17

                3

                3

                0

                N

                2–3

                0

                1

                0

                0–1

                0–2

                nd

                0

                0–1

                nd

                3

                2–3

                3

                0

                0

                0–1

                3

                M. avium

                NP10

                3

                3

                0

                N

                3

                3

                1

                1

                2

                3

                3

                0

                0

                3

                3

                1

                1

                1

                0

                0

                0

                M. cosmeticum

                NP6

                3

                3

                2

                N

                3

                0

                1

                0

                1

                2

                1

                2

                1

                nd

                2

                3

                2

                3

                2

                1

                2

                M. fortuitum

                NP12

                3

                3

                0

                S

                3

                3

                1

                3

                1

                3

                2

                1

                0

                3

                3

                3

                3

                3

                3

                0

                3

                M. fortuitum

                NP19

                3

                3

                0

                N

                3

                3

                0

                3

                3

                3

                3

                3

                1

                3

                3

                3

                3

                3

                3

                0

                2

                M. gordonae

                NP1

                2–3

                3

                0

                S

                3

                0

                1

                0

                1

                2

                0

                0

                0

                3

                1–3

                3

                0–3

                0–3

                0–1

                0–1

                1–3

                M. gordonae

                NP3

                3

                1

                0

                S

                3

                0

                2

                0

                0

                2

                0

                0

                0

                3

                1

                2

                0

                0

                0

                0

                0

                M. gordonae

                NP13

                1

                2

                0

                P

                3

                1

                0

                nd

                0

                1

                0

                0

                0

                3

                3

                3

                0

                1

                1

                2

                1

                M. gordonae

                NP22

                3

                3

                3

                S

                3

                3

                2

                0

                0

                0

                nd

                0

                0

                3

                0

                3

                1

                3

                3

                0

                nd

                M. mageritense

                NP5

                3

                3

                0

                N

                3

                2

                1

                0

                0

                2

                3

                0

                0

                3

                3

                3

                0

                3

                0

                0

                0

                M. nonchromogenicum

                NP4

                2

                3

                2

                N

                2

                1

                1

                0

                0

                0

                nd

                1

                1

                nd

                2

                2

                2

                2

                0

                0

                0

                M. sherrisii

                NP14

                1

                2

                0

                S

                1–2

                0

                0

                1–2

                0

                0

                1

                0

                1

                2

                1–2

                1–2

                1

                1

                2

                2

                1–2

                M. terrae

                NP11

                2

                3

                0

                N/S

                2–3

                0–3

                1, 2

                1

                0

                0–1

                0

                0

                0–1

                3

                3

                1, 2

                3

                0

                0

                0, 1

                3

                Phenotypes: 24°C, 36°C, 45°C: growth at temperature shown; pg: pigmentation (N, nonchromogen; P, photochromogen; S, scotochromogen); TCH: growth on thiophene-2-carboxylic acid hydrazide; nit: nitrate reduction; Tween: hydrolysis of Tween 80; NaCl: growth on 5% NaCl; Aryl 3. Aryl 15: arylsulfatase activity after 3 and 15 days of growth, respectively; ag: growth on nutrient agar; pic: growth on picric acid; β-gal: β-galactosidase activity; LJ: growth on Löwenstein-Jensen media; HA: growth on hydroxylamine 500 μg/ml; PNB: growth on p-nitrobenzoic acid; INH: isoniazid; RF: rifampicin; EMB: ethambutol; CIP: ciprofloxacin; OFL: ofloxacin. Responses are graded 0 (negative, no growth, no activity expressed) to 3 (positive, heavy growth, strong activity expressed); nd, not done. For patterns with multiple isolates, the result shown represents the most common phenotype(s) or the range of phenotypes observed.

                Overview of results

                The overall results of the two methods are summarized in Table 6. Among 434 NTM isolates, PRA-hsp65 provided correct species identification significantly more frequently than phenotypic/biochemical testing (392 (90.3%) vs 338 (77.9%), respectively; p < .0001, Fisher's exact test).
                Table 6

                Summary of concordance among species identification results obtained by PRA-hsp65, phenotypic evaluation and sequence analysis of the hsp65 gene.

                hsp65 sequence

                N

                PRA-hsp65

                Phenotypic identification

                  

                Concordant

                New Pattern

                Discordant

                Concordant

                Ambiguous

                Discordant

                Not donea

                321

                321

                --

                --

                321

                --

                --

                Done

                113

                71

                30

                12

                17

                63

                33

                Total

                434

                392 (90.3%)

                30 (6.9%)

                12 (2.8%)

                338 (77.9%)

                63 (14.5%)

                33 (7.6%)

                a Isolates for which species identification by PRA-hsp65 and phenotypic/biochemical evaluation were concordant were not sequenced. Based on prior reports by the authors and others, sequencing hsp65 in such isolates almost invariably confirms the species identification of the other methods.

                The four species or complex of NTM most commonly associated with clinically significant disease are M. avium complex, M. fortuitum complex, M. chelonae complex and M. kansasii. These represented 333 (76.7%) of the 434 isolates in this collection. PRA-hsp65 provided incorrect species identification for only 4 (1.2%) of these isolates and a new pattern for an additional 3 (0.9%). In contrast, phenotypic/biochemical testing provided incorrect assignments for 9 (2.7%) and ambiguous results for 31 (9.3%). Thus, the frequency of incorrect or uncertain species identification among these isolates of potential clinical importance was almost 6-fold higher for the phenotypic method than for PRA-hsp65 (40 (12.0%) vs. 7 (2.1%), respectively; p < .0001, Fisher's exact test).

                PRA-hsp65 algorithm

                Figures 1, 2 and 3 display an updated algorithm relating observed restriction fragments to particular species. We have included refinements of previously assigned fragment sizes based on our observations and analysis of available hsp65 sequences from validated mycobacterial species found online [32]. Sequences retrieved from GenBank [33] comprising the 441 bp Telenti fragment were analyzed using BioEdit, version 7.0.5.3. [34] and/or the DNASIS Max version 1 program (Hitashi Software Engineering Co., USA). BstEII restriction patterns were distributed in seven possible configurations: 440, 320-130, 320-120, 235-210, 235-130-85, 235-120-100, and 235-120-85. HaeIII fragment sizes were adjusted considering the nearest number multiple of 5, to facilitate interpretation of gel bands. These adjustments were performed based in our experience with analysis of more than 500 gels both visually and using the GelCompar program. HaeIII restriction fragments shorter than 50 bp were not taken in account as their discrimination in 4% agarose gels is often inaccurate. Different variants of PRA-hsp65 profiles from each species were numbered using Arabic numbers after the designation of the species, as reported in the PRASITE, except for M. avium, for which variants M. avium 1 and M. avium 2 were defined as reported in Leao et al. [20] and Smole et al. [27]. There were also PRA-hsp65 patterns frequently found in our routine work that had no sequence deposited. These patterns were included according to published data [1117] or the PRASITE [25]. Figures 2 and 3 also include the two new patterns we observed in two or more isolates (NP11 and NP1) and for which we propose PRA-hsp65 designations, M. terrae 4 and M. gordonae 10, respectively. The partial hsp65 gene sequences of these isolates have been deposited in GenBank [GenBank:EF601223 and GenBank:EF601222, respectively]. The figures also indicate the basic phenotypic characteristics (time for growth and pigment production) observed for each species.
                http://static-content.springer.com/image/art%3A10.1186%2F1471-2180-8-48/MediaObjects/12866_2007_Article_481_Fig1_HTML.jpg
                Figure 1

                Algorithm of PRA-hsp65patterns based on analysis of the 441 bp fragment of thehsp65gene. BstEII patterns: 440 bp, 320 bp/130 bp, 320 bp/115 bp. Columns 1 and 2: calculated BstEII and HaeIII fragment sizes in base pairs. Column 3: species names according to [32]. Column 4: PRA-hsp65 pattern type. Column 5: RGN: rapidly growing non-pigmented, RGS: rapidly growing scotochromogen, RGP: rapidly growing photochromogen, SGN: slowly growing non-pigmented, SGS: slowly growing scotochromogen, SGP: slowly growing photochromogen. Column 6: strain(s) used for hsp65 sequencing or reference of the publication describing this pattern.

                http://static-content.springer.com/image/art%3A10.1186%2F1471-2180-8-48/MediaObjects/12866_2007_Article_481_Fig2_HTML.jpg
                Figure 2

                Algorithm of PRA-hsp65 patterns based on analysis of the 441 bp fragment of the hsp65 gene. BstEII patterns: 235 bp/210 bp. Columns 1 and 2: calculated BstEII and HaeIII fragment sizes in base pairs. Column 3: species names according to [32]. Column 4: PRA-hsp65 pattern type. Column 5: RGN: rapidly growing non-pigmented, RGS: rapidly growing scotochromogen, RGP: rapidly growing photochromogen, SGN: slowly growing non-pigmented, SGS: slowly growing scotochromogen, SGP: slowly growing photochromogen. Column 6: strain(s) used for hsp65 sequencing or reference of the publication describing this pattern.

                http://static-content.springer.com/image/art%3A10.1186%2F1471-2180-8-48/MediaObjects/12866_2007_Article_481_Fig3_HTML.jpg
                Figure 3

                Algorithm of PRA-hsp65patterns based on analysis of the 441 bp fragment of thehsp65gene. BstEII patterns: 235 bp/130 bp/85 bp, 235 bp/120 bp/100 bp, 235 bp/120 bp/85 bp. Columns 1 and 2: calculated BstEII and HaeIII fragment sizes in base pairs. Column 3: species names according to [32]. Column 4: PRA-hsp65 pattern type. Column 5: RGN: rapidly growing non-pigmented, RGS: rapidly growing scotochromogen, RGP: rapidly growing photochromogen, SGN: slowly growing non-pigmented, SGS: slowly growing scotochromogen, SGP: slowly growing photochromogen. Column 6: strain(s) used for hsp65 sequencing or reference of the publication describing this pattern.

                Discussion

                The incidence of individual infections and outbreaks associated with NTM has risen dramatically over the past decade establishing these organisms as significant human pathogens. Traditionally, the identification of mycobacteria to the species level has relied upon biochemical tests, which require three to six weeks to complete. Biochemical identification, even when performed by skilled microbiologists, may yield uncertain or even misleading results because (a) the tests used are inherently poorly reproducible; (b) the expected phenotypes are not an absolute property of the species, but may exhibit substantial variability; and (c) the database of phenotypic characteristics is limited to common species [10].

                In recent years, DNA-based techniques have greatly facilitated identifying the species of NTM isolates and enabled a number of new species to be documented as infecting agents [3539]. These approaches can be applied to a single isolated colony and a definitive result can typically be obtained within a day. PRA-hsp65, first described by Telenti et al., is based on detection of restriction fragment polymorphisms in the hsp65 gene and thereby resolving the species of a mycobacterial isolate [11].

                In the present study, 434 NTM isolates from clinical specimens were analyzed by conventional phenotypic methods and by PRA-hsp65; further, those isolates for which the results from the two methods were discordant were analyzed using nucleotide sequencing of the hsp65 gene. For 63 (14.5%) isolates phenotypic methods could not provide a species identification and for almost a third of these isolates even the apparent Runyon classification proved inconsistent with conventional expectations. For an additional 33 (7.6%) isolates the phenotypic identification proved incorrect. Phenotypic variability among fresh clinical isolates has been observed in other studies [10, 40, 41].

                In contrast, PRA-hsp65 correctly identified over 90% of evaluable isolates using currently available databases of restriction digest patterns. For most of the remaining isolates, the PRA-hsp65 pattern observed was not previously reported. There were only 4 (1.2%) clinically significant isolates for which the current PRA algorithm indicated an incorrect species.

                PRA-hsp65 has proven similarly effective in other studies. Hafner et al. used 16S rDNA sequencing to analyze 126 isolates selected at random from a larger collection [17]. The hsp65 method correctly identified 120 (95.2%) of these isolates. They also sequenced 10 additional isolates from the larger collection that gave PRA-hsp65 patterns not previously reported. All these isolates represented environmental species rarely associated with clinically significant disease.

                Among our 434 isolates, 30 (6.9%) provided 13 PRA-hsp65 profiles not previously reported. Our series represents isolates cultured from varied clinical specimens collected in the metropolitan and surrounding areas of the city of Sao Paulo, Brazil. Most of the isolates with new PRA-hsp65 patterns were cultured from sputum. Many represented species typically considered non-pathogens; clinical correlation was not available and these isolates may reflect colonization by environmental organisms. Previous studies have similarly documented considerable species diversity as well as the genotypic diversity among mycobacteria isolates in Brazil [42, 43]. Sequence analysis confirmed that the new profiles were allelic variations within the species, consistent with previous studies [13, 17, 20]. Of interest, four profiles were represented by more than one isolate, suggesting that they are potentially prevalent lineages rather than singular mutation events.

                The most commonly identified new profile (designated NP1) was observed in 11 isolates, representing 20% of all M. gordonae in this collection. Comparison to the prototype M. gordonae sequence indicated two point mutations that resulted in the loss of two HaeIII sites and the addition of 95-bp fragment to the profile [GenBank:EF601222]. A similar profile was assigned to M. gordonae by da Silva Rocha et al. [13], although sequence confirmation was not reported. Hafner et al. also noted that M. gordonae is a particularly polymorphic species [17].

                The NP17 profile, demonstrated for five isolates, was identified by sequencing as M. arupense, a recently described species related to the M. terrae complex [18]. The NP14 profile, observed for three Mycobacterium simiae isolates, was similar profile to the M. simiae 3 pattern reported by Legrand et al. [22] as well as to the prototype M. simiae 1 pattern [11]. Sequencing confirmed that the nucleotide sequence is intermediate between those two strains. The sequence also matches that recently reported by Selvarangan et al., who proposed that their isolates represented a new species (M. sherrisii sp. nov) based on a distinct pattern of cellular fatty acids and a unique 16S rRNA gene [23]. The NP11 profile, represented by two isolates of M. terrae, was similar to a PRA-hsp65 pattern described by McNabb et al. [44] with the addition of a unique HaeIII restriction site [GenBank:EF601223].

                We would concur with Hafner et al. that additional work is required to define and standardize the most effective electrophoresis conditions for resolving hsp65 digests of mycobacteria [17]. In a recent multicenter study evaluating PRA-hsp65, variations related to gel preparation, running conditions and documentation tools all complicated the interpretation of digestion patterns [29].

                The ever-increasing amount of data available and the identification of new profiles make the analysis more complex. We present an updated PRA-hsp65 algorithm, which includes 174 patterns among 120 species and sub-species and have the basic cultural characteristics (rate of growth and pigment production). These core phenotypic traits can be readily determined and, as emphasized in a recent statement by the American Thoracic Society [45], can assist in confirming the molecular identification, detecting mixed cultures, and classifying species with indistinguishable PRA-hsp65 patterns.

                Despite the complexities noted above, PRA-hsp65 analysis proved both more rapid and more reliable than phenotypic methods; it was particularly effective at resolving the most common pathogenic species. Commercial DNA probes are available only for a very few species and their expense may be prohibitive in some settings. DNA sequencing is more definitive, but sequencing capability is not yet widely available in clinical laboratories.

                Conclusion

                Based on our extensive practical experience, we believe that PRA-hsp65 has the potential to provide clinicians with more timely, more accurate and, ultimately, more useful information and therefore deserves consideration by any clinical microbiology laboratory charged with performing species identification on NTM.

                Methods

                Mycobacterial isolates

                From January 2000 to January 2001, 5019 cultures were received at Institute Adolfo Lutz, São Paulo, Brazil for mycobacterial identification. M. tuberculosis complex was identified by direct observation of colony aspect and by Ziehl-Neelsen stained smears for presence/absence of cord formation. Cord-positive isolates with nonpigmented rough cultures were excluded from this study.

                A total of 439 isolates consistent with NTM were cultured from 435 (8.7%) specimens; five isolates were excluded because they could not be unambiguously resolved as NTM by the three methods used (phenotypic, PRA-hsp65 and sequencing), leaving a total of 434 isolates in the study. The specimens yielding NTM included sputum (280), blood (41), bronchial lavage (13), bone marrow (13), urine (7), skin biopsy (6), lymph node (5), feces (6), corneal scraping (4), pleural fluid (4), ascitic fluid (2), liver biopsy (2), liquor (1), gastric fluid (1), synovial fluid (1), abscess/secretion from unknown origin (11) and unknown (38). The majority (61.4%) of these specimens were from the Metropolitan Region of São Paulo, with 36.1% from elsewhere in São Paulo State and 2.5% from other States in Brazil.

                Conventional identification

                Isolates were identified based on phenotypic characteristics, including growth rate (fast/slow), pigment production, growth in different temperatures (26°C, 37°C and 45°C), biochemical tests (nitrate reduction, catalase activity, urease activity, tween 80 hydrolysis, arylsulfatase), specific chemicals (sodium chloride 5%, sodium salicylate), and growth in the presence of drugs (isoniazid 10 μg/ml, rifampicin 25 μg/ml, ethambutol 5 μg/ml, thiophen-2-carboxylic acid hydrazide 5 μg/ml, p-nitro-benzoic acid 0.5 μg/ml, cycloserine 30 μg/ml, ciprofloxacin 5 μg/ml, hydroxylamine 500 μg/ml, ofloxacin 2.5 μg/ml) [9, 31]. Some closely related mycobacterial species cannot be resolved by these biochemical tests. In such instances, isolates were designated as M. avium complex, M. terrae complex, M. chelonae complex or M. fortuitum complex, as appropriate.

                DNA extraction and PRA-hsp65 method

                For DNA extraction, a loop-full of organisms grown on Löwenstein-Jensen medium was suspended in 500 μl of ultrapure water, boiled for 10 min and frozen at -20°C for at least 18 h. Five microliters of DNA-containing supernatant were subjected to PCR amplification of the 441 bp of the gene hsp65 [11]. Separate aliquots of the PCR product were digested with BstEII and HaeIII, and the resulting restriction fragments separated by electrophoresis in a 4% agarose gel (Nusieve, FMC Bioproducts, Rockland, Maine USA) with 50 bp ladder as molecular size standard.

                Analysis of PRA-hsp65 results

                Gels were stained with ethidium bromide, photographed on a UV transilluminator, the images scanned, the restriction fragment sizes estimated using GelCompar II software, version 2.5 (AppliedMaths, St. Marten Latem, Belgium) and the patterns observed compared to the patterns reported on PRASITE [25], in publications [1117] or calculated in silico from sequences deposited in GenBank [33] using BioEdit, version 7.0.5.3 [34].

                hsp65 partial gene sequencing

                For those isolates for which conventional and PRA-hsp65 methods gave discordant or inconclusive results, the hsp65 amplicon was purified using Novagen Spin-prep Kit (Novagen, Canada) and then sequenced using BigDye terminator cycle sequencing reagents. Cycle sequencing was performed by using a Perkin-Elmer 9600 GeneAmp PCR system programmed for 25 cycles at 96°C for 20 s, 50°C for 10 s and 60°C for 4 min. Sequencing products were cleaned with CentriSep Spin Columns (Princeton Separations, Applied Biosystems) and then analyzed on a ABI Prism 377 sequencer (Perkin-Elmer).

                Sequence data analysis

                Data produced by the sequencer was automatically processed using the EGene platform [46]. The trace files were initially submitted to Phred [47] for base calling and quality assessment. Then, sequences were submitted to a quality filter that eliminated reads that did not present at least one window of 200 bases where 190 bases had phred quality above 15. After, low quality bases were trimmed from the sequence. For each sequence, the trimming procedure isolated a "good quality" subsequence. In this remaining subsequence, any window of 15 bases have at least 12 bases above the quality threshold of 15. After trimming, contaminant screening was performed using Blastn [48] against Homo sapiens, Salmonella typhimurium and Gallus gallus databases. Finally the clean isolates were identified by similarity using Blastn against a database of hsp65 genes. Sequences were considered a positive match when they presented a minimum similarity of 80 percent over a local alignment of at least 90 bases and ev-value of 1e-20. Species identification was confirmed if = 97% match was achieved, according to criteria proposed by McNabb et al. [44],. with any sequence deposited in databases and published.

                Declarations

                Acknowledgements

                Simone Guadagnucci Morillo is acknowledged for helping with DNA sequencing. This work was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Proc. N. 01/13389-6. Erica Chimara was the recipient of a fellowship from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES).

                Authors’ Affiliations

                (1)
                Instituto Adolfo Lutz
                (2)
                Instituto de Matemática e Estatística, Universidade de São Paulo
                (3)
                Division of Geographic Medicine and Infectious Diseases, Tufts University School of Medicine
                (4)
                Departamento de Microbiologia, Imunologia e Parasitologia, Universidade Federal de São Paulo

                References

                1. Roth A, Reischl U, Streubel A, Naumann L, Kroppenstedt RM, Habicht M, Fischer M, Mauch H: Novel diagnostic algorithm for identification of mycobacteria using genus-specific amplification of the 16S-23S rRNA gene spacer and restriction endonucleases. J Clin Microbiol 2000,38(3):1094–1104.PubMed
                2. Saiman L: The mycobacteriology of non-tuberculous mycobacteria. Paediatr Respir Rev 2004, 5 Suppl A:S221–3.View ArticlePubMed
                3. Sampaio JL, Chimara E, Ferrazoli L, da Silva Telles MA, Del Guercio VM, Jerico ZV, Miyashiro K, Fortaleza CM, Padoveze MC, Leao SC: Application of four molecular typing methods for analysis of Mycobacterium fortuitum group strains causing post-mammaplasty infections. Clin Microbiol Infect 2006,12(2):142–149.View ArticlePubMed
                4. Sampaio JL, Junior DN, de Freitas D, Hofling-Lima AL, Miyashiro K, Alberto FL, Leao SC: An outbreak of keratitis caused by Mycobacterium immunogenum. J Clin Microbiol 2006,44(9):3201–3207.View ArticlePubMed
                5. John T, Velotta E: Nontuberculous (atypical) mycobacterial keratitis after LASIK: current status and clinical implications. Cornea 2005,24(3):245–255.View ArticlePubMed
                6. Tiwari TS, Ray B, Jost KC Jr., Rathod MK, Zhang Y, Brown-Elliott BA, Hendricks K, Wallace RJ Jr.: Forty years of disinfectant failure: outbreak of postinjection Mycobacterium abscessus infection caused by contamination of benzalkonium chloride. Clin Infect Dis 2003,36(8):954–962.View ArticlePubMed
                7. Wallace RJ Jr., Brown BA, Griffith DE: Nosocomial outbreaks/pseudo-outbreaks caused by nontuberculous mycobacteria. Annu Rev Microbiol 1998, 52:453–490.View ArticlePubMed
                8. Freitas D, Alvarenga L, Sampaio J, Mannis M, Sato E, Sousa L, Vieira L, Yu MC, Martins MC, Hofling-Lima A, Belfort R Jr.: An outbreak of Mycobacterium chelonae infection after LASIK. Ophthalmology 2003,110(2):276–285.View ArticlePubMed
                9. Collins CH, Grange JM, Yates MD: Tuberculosis Bacteriology: Organization and Practice. 2nd edition Edition Oxford, Butterworth-Heinemann 1997.
                10. Springer B, Stockman L, Teschner K, Roberts GD, Bottger EC: Two-laboratory collaborative study on identification of mycobacteria: molecular versus phenotypic methods. J Clin Microbiol 1996,34(2):296–303.PubMed
                11. Telenti A, Marchesi F, Balz M, Bally F, Bottger EC, Bodmer T: Rapid identification of mycobacteria to the species level by polymerase chain reaction and restriction enzyme analysis. J Clin Microbiol 1993,31(2):175–178.PubMed
                12. Brunello F, Ligozzi M, Cristelli E, Bonora S, Tortoli E, Fontana R: Identification of 54 mycobacterial species by PCR-restriction fragment length polymorphism analysis of the hsp65 gene. J Clin Microbiol 2001,39(8):2799–2806.View ArticlePubMed
                13. da Silva Rocha A, Werneck Barreto AM, Dias Campos CE, Villas-Boas da Silva M, Fonseca L, Saad MH, Degrave WM, Suffys PN: Novel allelic variants of Mycobacteria isolated in Brazil as determined by PCR-restriction enzyme analysis of hsp65. J Clin Microbiol 2002,40(11):4191–4196.View ArticlePubMed
                14. Devallois A, Goh KS, Rastogi N: Rapid identification of mycobacteria to species level by PCR-restriction fragment length polymorphism analysis of the hsp65 gene and proposition of an algorithm to differentiate 34 mycobacterial species. J Clin Microbiol 1997,35(11):2969–2973.PubMed
                15. Leao SC, Martin A, Mejia GI, Palomino JC, Robledo J, Telles MAS, Portaels F: Practical handbook for the phenotypic and genotypic identification of mycobacteria. Brugges, Vanden BROELLE 2004.
                16. Rastogi N: An introduction to bacterial taxonomy, structure, drug resistance, and pathogenesis. Textbook-atlas of intestinal infections in AIDS (Edited by: Dioniosio D). Milano, Italia, Springer-Verlag 2003, 89–115.
                17. Hafner B, Haag H, Geiss HK, Nolte O: Different molecular methods for the identification of rarely isolated non-tuberculous mycobacteria and description of new hsp65 restriction fragment length polymorphism patterns. Mol Cell Probes 2004,18(1):59–65.View ArticlePubMed
                18. Cloud JL, Meyer JJ, Pounder JI, Jost KC Jr., Sweeney A, Carroll KC, Woods GL: Mycobacterium arupense sp. nov., a non-chromogenic bacterium isolated from clinical specimens. Int J Syst Evol Microbiol 2006,56(Pt 6):1413–1418.View ArticlePubMed
                19. Cooksey RC, de Waard JH, Yakrus MA, Rivera I, Chopite M, Toney SR, Morlock GP, Butler WR: Mycobacterium cosmeticum sp. nov., a novel rapidly growing species isolated from a cosmetic infection and from a nail salon. Int J Syst Evol Microbiol 2004,54(Pt 6):2385–2391.View ArticlePubMed
                20. Leao SC, Briones MR, Sircili MP, Balian SC, Mores N, Ferreira-Neto JS: Identification of two novel Mycobacterium avium allelic variants in pig and human isolates from Brazil by PCR-restriction enzyme analysis. J Clin Microbiol 1999,37(8):2592–2597.PubMed
                21. Leao SC, Sampaio JL, Martin A, Palomino JC, Portaels F: Profiling Mycobacterium ulcerans with hsp65. Emerg Infect Dis 2005,11(11):1795–1796.PubMed
                22. Legrand E, Goh KS, Sola C, Rastogi N: Description of a novel Mycobacterium simiae allelic variant isolated from Caribbean AIDS patients by PCR-restriction enzyme analysis and sequencing of hsp65 gene. Mol Cell Probes 2000,14(6):355–363.View ArticlePubMed
                23. Selvarangan R, Wu WK, Nguyen TT, Carlson LD, Wallis CK, Stiglich SK, Chen YC, Jost KC Jr., Prentice JL, Wallace RJ Jr., Barrett SL, Cookson BT, Coyle MB: Characterization of a novel group of mycobacteria and proposal of Mycobacterium sherrisii sp. nov. J Clin Microbiol 2004,42(1):52–59.View ArticlePubMed
                24. Turenne C, Chedore P, Wolfe J, Jamieson F, May K, Kabani A: Phenotypic and molecular characterization of clinical isolates of Mycobacterium elephantis from human specimens. J Clin Microbiol 2002,40(4):1230–1236.View ArticlePubMed
                25. PRASITE: Identification of Mycobacteria[http://​app.​chuv.​ch/​prasite]
                26. Chimara E, Giampaglia CM, Martins MC, Telles MA, Ueki SY, Ferrazoli L: Molecular characterization of Mycobacterium kansasii isolates in the State of Sao Paulo between 1995–1998. Mem Inst Oswaldo Cruz 2004,99(7):739–743.PubMed
                27. Smole SC, McAleese F, Ngampasutadol J, Von Reyn CF, Arbeit RD: Clinical and epidemiological correlates of genotypes within the Mycobacterium avium complex defined by restriction and sequence analysis of hsp65. J Clin Microbiol 2002,40(9):3374–3380.View ArticlePubMed
                28. Jimenez-Pajares MS, Herrera L, Valverde A, Saiz P, Saez-Nieto JA: [Phenotypic and genotypic characteristics of Mycobacterium kansasii strains isolated in Spain (2000–2003)]. Enferm Infecc Microbiol Clin 2005,23(5):254–258.View ArticlePubMed
                29. Leao SC, Bernardelli A, Cataldi A, Zumarraga M, Robledo J, Realpe T, Mejia GI, da Silva Telles MA, Chimara E, Velazco M, Fernandez J, Rodrigues PA, Guerrero MI, Leon CI, Porras TB, Rastogi N, Goh KS, Suffys P, da Silva Rocha A, dos Santos Netto D, Ritacco V, Lopez B, Barrera L, Palomino JC, Martin A, Portaels F: Multicenter evaluation of mycobacteria identification by PCR restriction enzyme analysis in laboratories from Latin America and the Caribbean. J Microbiol Methods 2005,61(2):193–199.View ArticlePubMed
                30. Prammananan T, Cheunoy W, Na-Ubol P, Tingtoy N, Srimuang S, Chaiprasert A: Evaluation of polymerase chain reaction and restriction enzyme analysis for routine identification of mycobacteria: accuracy, rapidity, and cost analysis. Southeast Asian J Trop Med Public Health 2005,36(5):1252–1260.PubMed
                31. Tsukamura M, Sekine K, Yokota A, Kuze A, Shibata M, Sato K: Lung infection due to Mycobacterium xenopi: report of the first case in Japan. Microbiol Immunol 1984,28(1):123–127.PubMed
                32. J.P. Euzéby: List of Prokaryotic names with Standing in Nomenclature[http://​www.​bacterio.​cict.​fr/​m/​mycobacterium.​html]
                33. GenBank[http://​www.​ncbi.​nlm.​nih.​gov]
                34. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 1999, 41:95–98.
                35. da Silva Rocha A, da Costa Leite C, Torres HM, de Miranda AB, Pires Lopes MQ, Degrave WM, Suffys PN: Use of PCR-restriction fragment length polymorphism analysis of the hsp65 gene for rapid identification of mycobacteria in Brazil. J Microbiol Methods 1999,37(3):223–229.View ArticlePubMed
                36. Mun HS, Kim HJ, Oh EJ, Kim H, Bai GH, Yu HK, Park YG, Cha CY, Kook YH, Kim BJ: Mycobacterium seoulense sp. nov., a slowly growing scotochromogenic species. Int J Syst Evol Microbiol 2007,57(Pt 3):594–599.View ArticlePubMed
                37. Reischl U, Melzl H, Kroppenstedt RM, Miethke T, Naumann L, Mariottini A, Mazzarelli G, Tortoli E: Mycobacterium monacense sp. nov. Int J Syst Evol Microbiol 2006,56(Pt 11):2575–2578.View ArticlePubMed
                38. Turenne CY, Cook VJ, Burdz TV, Pauls RJ, Thibert L, Wolfe JN, Kabani A: Mycobacterium parascrofulaceum sp. nov., novel slowly growing, scotochromogenic clinical isolates related to Mycobacterium simiae. Int J Syst Evol Microbiol 2004,54(Pt 5):1543–1551.View ArticlePubMed
                39. Wallace Jr RJ, Brown-Elliott BA, Brown J, Steigerwalt AG, Hall L, Woods G, Cloud J, Mann L, Wilson R, Crist C, Jost KC Jr., Byrer DE, Tang J, Cooper J, Stamenova E, Campbell B, Wolfe J, Turenne C: Polyphasic characterization reveals that the human pathogen Mycobacterium peregrinum type II belongs to the bovine pathogen species Mycobacterium senegalense. J Clin Microbiol 2005,43(12):5925–5935.View Article
                40. Kellogg JA, Bankert DA, Withers GS, Sweimler W, Kiehn TE, Pfyffer GE: Application of the Sherlock Mycobacteria Identification System using high-performance liquid chromatography in a clinical laboratory. J Clin Microbiol 2001,39(3):964–970.View ArticlePubMed
                41. Cook VJ, Turenne CY, Wolfe J, Pauls R, Kabani A: Conventional methods versus 16S ribosomal DNA sequencing for identification of nontuberculous mycobacteria: cost analysis. J Clin Microbiol 2003,41(3):1010–1015.View ArticlePubMed
                42. Barreto AMW, Campos CED: Micobacterias nao-tuberculosas no Brasil. Bol Pneum Sanit 2000,8(1):23–32.
                43. Ueki SYM, Martins MC, Telles MAS, Virgilio MC, Giampaglia CMS, Chimara E, Ferrazoli L: Nontuberculous mycobacteria: species diversity in Sao Paulo state, Brazil. J Bras Patol Med Lab 2005,41(1):1–8.View Article
                44. McNabb A, Eisler D, Adie K, Amos M, Rodrigues M, Stephens G, Black WA, Isaac-Renton J: Assessment of partial sequencing of the 65-kilodalton heat shock protein gene (hsp65) for routine identification of Mycobacterium species isolated from clinical sources. J Clin Microbiol 2004,42(7):3000–3011.View ArticlePubMed
                45. Griffith DE, Aksamit T, Brown-Elliott BA, Catanzaro A, Daley C, Gordin F, Holland SM, Horsburgh R, Huitt G, Iademarco MF, Iseman M, Olivier K, Ruoss S, von Reyn CF, Wallace RJ Jr., Winthrop K: An official ATS/IDSA statement: diagnosis, treatment, and prevention of nontuberculous mycobacterial diseases. Am J Respir Crit Care Med 2007,175(4):367–416.View ArticlePubMed
                46. Durham AM, Kashiwabara AY, Matsunaga FT, Ahagon PH, Rainone F, Varuzza L, Gruber A: EGene: a configurable pipeline generation system for automated sequence analysis. Bioinformatics 2005,21(12):2812–2813.View ArticlePubMed
                47. Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res 1998,8(3):175–185.PubMed
                48. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990,215(3):403–410.PubMed

                Copyright

                © Chimara et al. 2008

                This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

                Advertisement