Using polymerase chain reaction with primers based on the plcB-plcC intergenic region to detect Mycobacterium tuberculosis in clinical samples

Abstracts

OBJETIVO: Desenvolver um sistema para o diagnóstico molecular da tuberculose por reação em cadeia da polimerase, do inglês polymerase chain reaction (PCR), pela construção de primers baseados na diferença da organização de uma região intergênica da fosfolipase (phospholipase) C (região plcB-plcC), que diferencia Mycobacterium tuberculosis das outras micobactérias. MÉTODOS: Um produto de PCR com o tamanho esperado (432 pb) foi obtido somente de M. tuberculosis e M. africanum. Um total de 33 isolados micobacterianos e 273 amostras clínicas de pacientes com suspeita de tuberculose foram examinados. Estes foram submetidos ao estudo comparativo da técnica de PCR contra o cultivo. RESULTADOS: Os dados mostraram 93,8% de exatidão para PCR contra o cultivo (p < 0,0001), 93,1% de sensibilidade (IC: 88,7-96,0) e especificidade de 96,4% (IC: 96,1-99,4). O valor de Kappa foi de 0,82, demonstrando um alinhamento perfeito para a verificação do grau de concordância entre os testes. CONCLUSÃO: O uso da região plcB-plcC para a amplificação por PCR é mostrado como uma ferramenta importante e de confiança para o diagnóstico específico de tuberculose nas amostras clínicas analisadas.

Reação em cadeia da polimerase; Diagnóstico; Tuberculose; Mycobacterium tuberculosis

ARTIGO ORIGINAL

Using polymerase chain reaction with primers based on the plcB-plcC intergenic region to detect Mycobacterium tuberculosis in clinical samples^*,**

Hermides Pinto Júnior^I; Claudia Giuliano Bica^II; Moisés Palaci^III; Reynaldo Dietze^IV; Luiz Augusto Basso^V; Diógenes Santiago Santos^VI

^IDiretor de Projetos e Desenvolvimento do Amplicon, Porto Alegre (RS) Brasil

^IIProfessora Assistente da Fundação Faculdade Federal de Ciências Médicas de Porto Alegre – FFFCMPA – Porto Alegre (RS) Brasil

^IIIProfessor Adjunto de Microbiologia da Universidade Federal do Espírito Santo – UFES – Espírito Santo (ES) Brasil

^IVCoordenador do Núcleo de Doenças Infecciosas da Universidade Federal do Espírito Santo – UFES – Espírito Santo (ES) Brasil

^VDoutor em Ciências na área de Enzimologia. Pontíficia Universidade Católica do Rio Grande do Sul – PUCRS – Porto Alegre (RS) Brasil

^VIDoutor em Ciências na área de Biologia Molecular. Pontíficia Universidade Católica do Rio Grande do Sul – PUCRS – Porto Alegre (RS) Brasil

^{Endereço para correspondência}

ABSTRACT

OBJECTIVE: To develop a system for the molecular diagnosis of tuberculosis by polymerase chain reaction (PCR), constructing primers based on the difference in gene organization of the intergenic region of phospholipase C (plcB-plcC region), which differentiates Mycobacterium tuberculosis from other mycobacteria.

METHODS: A PCR product of the expected size (432 bp) was obtained from M. tuberculosis and M. africanum only. A total of 33 mycobacterial isolates and 273 clinical samples from patients suspected of having tuberculosis were examined. These were used in the comparative study of the PCR technique versus culture.

RESULTS: For PCR versus culture, the data showed 93.8% accuracy (p < 0.0001), 93.1% sensitivity (CI: 88.7-96.0), and 96.4% specificity (CI: 96.1-99.4). The Kappa value (0.82) shows that there was a near-perfect concordance between the two tests.

CONCLUSION: The use of the plcB-plcC region in PCR amplification was found to be an important and reliable tool for the specific diagnosis of tuberculosis in the samples analyzed.

Keywords: Polymerase chain reaction; Diagnosis; Tuberculosis; Mycobacterium tuberculosis.

Introduction

Fast and accurate diagnosis is an important element of global health measures to control tuberculosis (TB). Moreover, since 90% of TB cases occur in developing countries,⁽¹⁾ a cost-effective diagnostic tool for TB that could be routinely used would be of great benefit. A number of methods based on nucleic acid amplification, including polymerase chain reaction-restriction enzyme analysis (PCR-REA),⁽²⁾ commercial DNA assays,⁽³⁾ PCR-based sequencing technologies,⁽⁴⁾ and, more recently, real-time PCR assays,⁽⁵⁾ have been developed for use in the identification of mycobacteria. A shortcoming of the PCR-REA method is its inability to distinguish the small differences in band size between some species of mycobacteria,⁽²⁾ whereas a major drawback of commercial assays is their high cost. In addition, it has been reported that a commercial direct amplification test designed for the in vitro diagnostic detection of Mycobacterium tuberculosis complex rRNA can produce false-positive results for M. tuberculosis in AIDS patients infected with M. kansasii or M. avium, as well as in HIV-negative patients with chronic lung disease caused by M. avium.⁽⁶⁾ Another commercial test has such a low sensitivity for smear-negative specimens that the United States Food and Drug Administration has approved it only for direct detection of M. tuberculosis in acid-fast bacilli smear-positive respiratory specimens.^(7,8) In addition, the ligase chain reaction test has not demonstrated a sufficient degree of accuracy with smear-negative specimens to warrant its use on a routine basis.⁽⁹⁾ Although the use of nucleic acid sequencing is highly accurate, it is labor-intensive and requires expensive equipment. Therefore, various laboratories use in-house methods for the detection of M. tuberculosis. These assays have a wide variety of nucleic acid targets, such as insertion sequence (IS)6110, ribosomal ribonucleic acid genes, 65-kDa antigen, and heat shock protein 65 gene.⁽⁸⁾ In such in-house assays, IS6110 is probably the target most widely used in the molecular detection of M. tuberculosis.^(7,10) However, M. tuberculosis strains that do not contain IS6110 have been reported,^(5,11,12) and IS6110 homologs have been found in M. fortuitum, M. avium-M. intracellulare complex, M. kansasii, M. xenopi, M. malmoense, and M. chelonei clinical isolates,^(13,14) as well as in M. bovis.⁽¹⁵⁾ According to data from the Pan American Health Organization and World Health Organization, 7000 new cases of human TB caused by M. bovis are reported in South America each year.⁽¹⁶⁾ Differentiation between M. tuberculosis and M. bovis is particularly important, since the latter is resistant to pyrazinamide, a drug commonly used in the treatment of TB.⁽¹⁷⁾ However, current methodologies to distinguish between M. tuberculosis and M. bovis have been largely unsuccessful due to their lack of resolution,⁽¹⁸⁾ as well as to the fact that skilled laboratory personnel and relatively expensive consumable reagents are required.^(19,20)

Upon performing comparative hybridization experiments, eleven regions present in M. tuberculosis H37Rv were found to be absent from one or more virulent strains of M. bovis⁽²¹⁾; among them, a cluster of three phospholipase C genes (plcA, plcB, and plcC).⁽²¹⁾ Accordingly, we have attempted to develop a rapid, low-cost, direct PCR method based on the amplification of a DNA fragment including the plcB-plcC intergenic region for detection of M. tuberculosis, which will allow M. tuberculosis to be the differentiated from M. bovis in clinical samples.

Methods

A total of 33 mycobacterial isolates and 273 clinical samples were tested. All mycobacterial strains (M. tuberculosis H37Rv, M. africanum, M. bovis, M. abscessus, M. aichiense, M. avium, M. asiaticum, M. aurum, M. chitae, M. chubuense, M. duvalii, M. flavescens, M. fortuitum, M. gastri, M. gordonae, M. intracellulare, M. kansasii, M. lentiflavum, M. mageritense, M. marinum, M. microti, M. neoaurum, M. nonchromogenicum, M. parafortuitum, M. peregrinum, M. porcinum, M. scrofulaceum, M. simiae, M. smegmatis, M. szulgai, M. triviale, M. vaccae, and M. xenopi) were identified by conventional culturing and biochemical tests in the Molecular Biology Mycobacteria Laboratory of the Oswaldo Cruz Foundation (Rio de Janeiro, Brazil). The clinical samples were obtained patients treated at the Sanatório Partenon Hospital (Porto Alegre, Brazil) or the Espírito Santo State Center for Infectious Diseases (Vitória, Brazil). Sputum specimens were collected from the respiratory tract of patients. One aliquot was inoculated in Löwenstein-Jensen medium, and 200 µl were used for DNA purification/PCR amplification. All clinical isolates in this study were identified on the basis of conventional techniques, including the determination of microbiological characteristics and biochemical tests.

The mycobacterial strains and clinical samples were inactivated at 95 ºC for 10 min. The DNA was extracted using QIAamp Tissue DNA Kit (Qiagen, Hilden, Germany), according to the manufacturer's instructions. Purified DNA (10 µl) was added to 40 µl of reaction mixture containing 20 mM Tris-HCl (pH 8.4), 50 mM KCl, 2.0 mM MgCl₂, 1.0 U Platinum Taq DNA polymerase (Invitrogen Life Technologies, Gaithersburg, MD, USA), and 0.2 mM dNTPs (Invitrogen Life Technologies), together with 10 pmol each of the primers BC3 (5'-TCGACGGCAGAAAGCGTGCC-3') and BC5 (5'-GCGGCTCAATGCGCTTCCG-3'). The primers BC3 (complement of bases 200266 to 200285) and BC5 (bases 200679 to 200697), which are based on the M. tuberculosis H37Rv sequence (EMBL accession number BX842579), amplified a 432-bp DNA fragment comprising the 3' end of the plcB gene, the intergenic region, and the 5' end of the plcC gene. PCR amplification was performed as follows: initial denaturation at 94 ºC for 5 min; 35 cycles of denaturation at 94 ºC for 30 s, annealing at 65 ºC for 30 s, extension at 72 ºC for 30 s; and a final extension step at 72 ºC for 5 min. The PCR products were analyzed by agarose gel electrophoresis (2% w/v) and stained with ethidium bromide. The clinical samples were also amplified with specific primers for human b-globin gene GH20 and PC04⁽²²⁾ in order to detect inhibitory samples.

Direct sequencing of PCR fragments was performed in a commercial laboratory (CATG, Porto Alegre, Brazil) using the ABI PRISM BigDye Terminator kit (PE Applied Biosystems, Foster City, CA, USA), and the output was analyzed with an ABI 3700 DNA sequencer. The PCR amplification primers were also used as sequencing primers, and at least a single base overlap from two directions for each was usually achieved. The FASTA program⁽²³⁾ was used to compare the amplified sequence with a DNA sequence from the GenBank database.

The collected data were analyzed using SPSS software (version 1.0 for Windows). Pearson's chi-square was used to assess relationships between categorical variables, and the data were interpreted using the Kappa table.

Results

In order to optimize the PCR, we tested amplifications with different concentrations of MgCl₂ (1.5-2.5 mM). The best amplification efficiency was achieved with MgCl₂ at a concentration of 2 mM. A strong, 432-bp product, without nonspecific fragments and without abundant primer dimers, was observed. A PCR product of the expected size was obtained for M. tuberculosis and M. africanum. Particularly encouraging was the absence of any PCR product observed for M. bovis (Figure 1). In addition, no amplification was observed using DNA from other organisms, such as humans, Escherichia coli, Streptococcus pneumoniae, and Staphylococcus aureus. When we tested amplification with MgCl₂ concentration of 2.5 mM, a faint 445-bp band was amplified with M. marinum, whereas all other mycobacterial strains tested resulted either in PCR products with sizes different from those expected or in no observable PCR product (data not shown). All further amplifications were then performed with 2 mM MgCl₂. The 432-bp DNA fragment was isolated and sequenced. It was found to be identical to the 200266-200697 nucleotide sequence of M. tuberculosis H37Rv deposited in the GenBank (BX842579).

Of the 273 clinical samples analyzed by PCR (Table 1), 93.8% were in agreement with the culture results (p < 0.001). Only 2 of the 55 culture-negative samples were identified as positive in the PCR. However, 15 of the 203 culture-positive samples were identified as negative in the PCR. Comparing the PCR results with the culture results, the sensitivity, specificity, positive predictive value, and negative predictive value were, respectively, 93.1% (CI: 88.7-96.0), 96.4% (CI: 96.1-99.4), 99% (CI: 96.1-99.8), and 77.9% (CI: 65.9-86.7). The Kappa value (0.82) showed that there was a near-perfect concordance between the two tests. These results compare favorably with those compiled in a recent review of the literature.⁽⁸⁾

Discussion

The rapid identification of M. tuberculosis is crucial to optimizing patient recovery. Species identification can take several weeks to complete, during which time the patient might receive inappropriate antimycobacterial agents. This is particularly worrisome if the patient is infected with M. bovis rather than M. tuberculosis, since M. bovis is resistant to pyrazinamide, a first-line antitubercular drug.⁽¹⁹⁾ In the present study, we designed and tested the primers BC3 and BC5, which differentiate M. tuberculosis from M. bovis. The genomic fragment designated mtp40 has been tested as a species-specific DNA target for the identification of M. tuberculosis.⁽²⁴⁾ However, the mtp40 sequence has been shown to be absent from a number of M. tuberculosis strains.^(18,25) Interestingly, the mtp40 sequence has been shown to be part of the plcA gene.⁽²⁵⁾ In the present study, we used primers flanking the junction of the 3' end of the plcB gene and the 5' end of the plcC gene. In addition, our attempts to amplify a DNA fragment comprising the 3' end of the plcA gene, the intergenic region, as well as one comprising the 5' end of the plcB gene, failed (data not shown).

The absence of any observed PCR product for M. bovis was expected, since the deletion of three plc-encoding genes (plcA, plcB, and plcC) in M. bovis (virulent strains and BCG strains) had already been reported by some authors.⁽²¹⁾

Our PCR assay cannot differentiate M. tuberculosis from M. africanum, but both are susceptible to the same chemotherapeutic agents. Some authors⁽²⁶⁾ developed a PCR-based method to differentiate the subspecies of the M. tuberculosis complex on the basis of genomic deletions. Seven primer pairs were used in order to construct the M. tuberculosis complex PCR typing. However, the authors found that their method was unable to differentiate M. tuberculosis from M. africanum. The PCR method described in the present study was also unable to differentiate between M. tuberculosis and M. africanum. However, we used only one pair of primers, and the interpretation of results was straightforward. Real-time PCR assays for the detection, differentiation, and quantitation of M. tuberculosis have recently been developed.⁽⁵⁾ However, none of them can differentiate M. tuberculosis from M. africanum or from M. bovis.

The PCR results are in good agreement with the culture data, except for the 15 samples that were identified as negative in the PCR and positive in culture (out of 203 culture-positive samples) and the 2 samples that were identified as positive in the PCR and negative in culture (out of 55 culture-negative samples). In the case of the false-negatives, sensitivity might have been increased by the fact that more than one sample was collected from each patient.⁽¹⁰⁾ In contrast, the false-positives might be attributable to the inclusion of paucibacillary patients presenting symptoms indicating infection with TB.

Within the M. tuberculosis complex, bacterial plc-encoding genes have been reported to play a role in the pathogenesis of many bacteria. Genetic variations (IS6110 insertions and deletions) have been observed in the plcA-plcB-plcC and plcD regions.^(27,28) At M. tuberculosis plcC positions 19849 and 19668,⁽²⁷⁾ as well as 19589, 19645, and 19848,⁽²⁹⁾ IS6110 insertions have been observed. These findings indicate that this domain is a preferential integration region, defined as a <500-bp chromosomal domain.⁽²⁷⁾ One group of authors⁽³⁰⁾ reported that, in M. tuberculosis, there was a much higher frequency of IS6110 insertion and deletion in the plcD gene than in the plcA, plcB, and plcC genes. However, these alterations did not affect the analyses of the clinical samples studied. The BC3/BC5 primer pair amplified a 432-bp fragment from plcB-plcC at position 19384-19403 (GenBank accession number Z83860). Therefore, it is out of the chromosomal domain.

Herein, we have describes a PCR method based on amplification of the plcB-plcC intergenic sequence. This sequence appears to be a promising target for use in differentiating M. tuberculosis from M. bovis in clinical samples.

Acknowledgments

We would like to thank the Partenon Hospital (Porto Alegre, Brazil) and Philip Noel Suffys, respectively, for providing some of the clinical samples and mycobacterial strains used in this study. This study received financial support in the form of grants from the Conselho Nacional de Desenvolvimento Científico e Tecnológico/Ministério da Ciência e Tecnologia (CNPq/MCT, National Council for Scientific and Technological Development/Ministry of Science and Technology; grant no. 62.00555/01-4 PADCT III/MILÊNIO), the Financiadora de Estudos e Projetos (FINEP, Funding Body for Studies and Projects), and the Programa de Núcleos de Excelência (PRONEX; grant no. 661028/1998-4).

Referências

Recebido para publicação em 16/5/2006

Aprovado, após revisão, em 28/11/2006

• 1. Schroeder EK, de Souza N, Santos DS, Blanchard JS, Basso LA. Drugs that inhibit mycolic acid biosynthesis in Mycobacterium tuberculosis. Curr Pharm Biotechnol. 2002;3(3):197-225.
• 2. Roth A, Reischl U, Streubel A, Naumann L, Kroppenstedt RM, Habicht M, et al. 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-104.
• 3. Wang SX, Tay L. Evaluation of three nucleic acid amplification methods for direct detection of Mycobacterium tuberculosis complex in respiratory specimens. J Clin Microbiol. 1999;37(6):1932-4.
• 4. Cloud JL, Neal H, Rosenberry R, Turenne CY, Jama M, Hillyard DR, et al. Identification of Mycobacterium spp. by using a commercial 16S ribosomal DNA sequencing kit and additional sequencing libraries. J Clin Microbiol. 2002;40(2):400-6.
• 5. Broccolo F, Scarpellini P, Locatelli G, Zingale A, Brambilla AM, Cichero P, et al. Rapid diagnosis of mycobacterial infections and quantitation of Mycobacterium tuberculosis load by two real-time calibrated PCR assays. J Clin Microbiol. 2003;41(10):4565-72
• 6. Jorgensen JH, Salinas JR, Paxson R, Magnon K, Patterson JE, Patterson TF. False-positive Gen-Probe direct Mycobacterium tuberculosis amplification test results for patients with pulmonary M. kansasii and M. avium infections. J Clin Microbiol. 1999;37(1):175-8.
• 7. Eing BR, Becker A, Sohns A, Ringelmann R. Comparison of Roche Cobas Amplicor Mycobacterium tuberculosis assay with in-house PCR and culture for detection of M. tuberculosis. J Clin Microbiol. 1998;36(7):2023-9.
• 8. Kaul KL. Molecular detection of Mycobacterium tuberculosis: impact on patient care. Clin Chem. 2001;47(8):1553-8.
• 9. Garrino MG, Glupczynski Y, Degraux J, Nizet H, Delmée M. Evaluation of the Abbott LCx Mycobacterium tuberculosis assay for direct detection of Mycobacterium tuberculosis Complex in Human Samples. J Clin Microbiol. 1999;37(1):229-32.
• 10. Rodrigues VF, Queiroz-Mello FC, Ribeiro MO, Fonseca L, Kritski AL, Rossetti ML, et al. Detection of Mycobacterium avium in blood samples of patients with AIDS by using PCR. J Clin Microbiol. 2002;40(6):2297-9.
• 11. Sreevatsan S, Pan X, Stockbauer KE, Connell ND, Kreiswirth BN, Whittam TS, et al. Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination. Proc Natl Acad Sci USA. 1997;94(18):9869-74.
• 12. Agasino CB, Ponce de Leon A, Jasmer RM, Small PM. Epidemiology of Mycobacterium tuberculosis strains in San Francisco that do not contain IS6110. Int J Tuberc and Lung Dis. 1998;2(6):518-20.
• 13. Kent L, McHugh TD, Billington O, Dale JW, Gillespie SH. Demonstration of homology between IS6110 of Mycobacterium tuberculosis and DNAs of other Mycobacterium spp.? J Clin Microbiol. 1995;33(9):2290-3.
• 14. McHugh TD, Newport LE, Gillespie SH. IS6110 homologs are present in multiple copies in mycobacteria other than tuberculosis-causing mycobacteria. J Clin Microbiol. 1997;35(7):1769-71.
• 15. Plikaytis BB, Eisenach KD, Crawford JT, Shinnick TM. Differentiation of Mycobacterium tuberculosis and Mycobacterium bovis BCG by a polymerase chain reaction assay. Mol Cell Probes.1991;5(3):215-9.
• ¹⁶
  Pan American Health Organization. Health conditions in the Americas, 1992 edition. Washington, D.C.: Pan American Health Organization, Pan American Sanitary Bureau, Regional Office of the World Health Organization;1992.
• 17. Kataria Y P. Observations on human infection with Mycobacterium bovis. Tubercle. 1969;50(1):14-21.
• 18. Weil A, Plikaytis BB, Butler WR, Woodley CL, Shinnick TM. The mtp40 gene is not present in all strains of Mycobacterium tuberculosis. J Clin Microbiol. 1996;34(9):2309-11.
• 19. Niemann S, Richter E, Rüsch-Gerdes S. Differentiation among members of the Mycobacterium tuberculosis complex by molecular and biochemical features: evidence for two pyrazinamide-susceptible subtypes of M. bovis. J Clin Microbiol. 2000;38(1):152-7.
• 20. Caws M, Drobniewski FA. Molecular techniques in the diagnosis of Mycobacterium tuberculosis and the detection of drug resistance. Ann N Y Acad Sci. 2001;(953):138-45.
• 21. Behr MA, Wilson MA, Gill WP, Salamon H, Schoolnik GK, Rane S, et al. Comparative genomics of BCG vaccines by whole-genome DNA microarray. Science. 1999;284(5419):1520-3.
• 22. Bauer HM, Manos MM. PCR Detection of Genital Human Papillomavirus. In. Persing DH, Smith TF, Tenover FC, White TJ, editors. Diagnostic molecular microbiology: principals and applications. Washington, D.C.: American Society for Microbiology; 1993. p. 407-13.
• 23. Pearson WR, Lipman DJ. Improved tools for biological sequence comparison. Proc Natl Acad Sci U S A. 1988;85(8):2444-8.
• 24. Liébana E, Aranaz A, Francis B, Cousins D. Assessment of genetic markers for species differentiation within the Mycobacterium tuberculosis complex. J Clin Microbiol. 1996;34(4):933-8.
• 25. Vera-Cabrera L, Howard ST, Laszlo A, Johnson WM. Analysis of genetic polymorphism in the phospholipase region of Mycobacterium tuberculosis. J Clin Microbiol. 1997;35(5):1190-5.
• 26. Huard RC, de Oliveira Lazzarini LC, Butler WR, van Soolingen D, Ho JL. PCR-based method to differentiate the subspecies of the Mycobacterium tuberculosis Complex on the basis of genomic deletions. J Clin Microbiol. 2003;41(4):1637-50.
• 27. Viana-Niero C, de Haas PE, van Soolingen D, Leão SC. Analysis of genetic polymorphisms affecting the four phospholipase C (plc) genes in Mycobacterium tuberculosis complex clinical isolates. Microbiology. 2004;150(Pt 4):967-78.
• 28. Kong Y, Cave MD, Yang D, Zhang L, Marrs CF, Foxman B, et al. Distribution of insertion- and deletion-associated genetic polymorphisms among four mycobacterium tuberculosis phospholipase c genes and associations with extrathoracic tuberculosis: a population-based study. J. Clin. Microbiol. 2005;43(12):6048-53.
• 29. Viana-Niero C, Gutierrez C, Sola C, Filliol I, Boulahbal F, Vincent V, et al. Genetic diversity of Mycobacterium africanum clinical isolates based on IS6110-restriction fragment length polymorphism analysis, spoligotyping, and variable number of tandem DNA repeats. J Clin Microbiol. 2001;39(1):57-65.
• 30. Talarico S, Durmaz R, Yang Z. Insertion- and deletion-associated genetic diversity of mycobacterium tuberculosis phospholipase c-encoding genes among 106 clinical isolates from turkey. J Clin Microbiol. 2005;43(2):533-8.

*

Trabalho realizado no Centro de Pesquisa em Biologia Molecular e Funcional, Instituto de Pesquisas Biomédicas, Pontifícia Universidade Católica do Rio Grande do Sul – PUCRS – Porto Alegre (RS) Brasil.

**

A versão completa em português deste artigo está disponível em

www.jornaldepneumologia.com.br

Publication Dates

• Publication in this collection
  24 Oct 2007
• Date of issue
  Aug 2007

History

• Received
  16 May 2006
• Accepted
  28 Nov 2006