- Research article
- Open Access
Identification of pathogenic Leptospira species by conventional or real-time PCR and sequencing of the DNA gyrase subunit B encoding gene
BMC Microbiologyvolume 6, Article number: 95 (2006)
Leptospira is the causative genus of the disease, leptospirosis. Species identification of pathogenic Leptospira in the past was generally performed by either DNA-DNA hybridisation or 16s rRNA gene sequencing. Both methods have inherent disadvantages such as the need for radio-labelled isotopes or significant homology between species. A conventional and real-time PCR amplification and sequencing method was developed for an alternate gene target: DNA gyrase subunit B (gyrB). Phylogenetic comparisons were undertaken between pathogenic Leptospira 16srRNA and gyrB genes using clustering and minimum evolution analysis. In addition 50 unidentified Leptospira isolates were characterised by gyrB sequencing and compared with conventional 16s rRNA sequencing.
A conventional and real-time PCR methodology was developed and optimised for the amplification of the gyrB from pathogenic Leptospira species. Non pathogenic and opportunistic Leptospira species such as L. fainei and L. broomi were not amplified. The gyrB gene shows greater nucleotide divergence (3.5% to 16.1%) than the 16s rRNA gene (0.1% to 1.4%). Minimum evolution analysis reveals that the gyrB has a different evolution topology for L. kirschneri and L. interrogans. When the two genes were compared for the identification of the 50 unknown isolates there was 100% agreement in the results.
This research has successfully developed a methodology for the identification of pathogenic Leptospira using an alternate gene to 16s rRNA. The gyrB encoding gene shows higher nucleotide/evolutionary divergence allowing for superior identification and also the potential for the development of DNA probe based identification.
Leptospirosis is the zoonotic disease caused by members of the genus, Leptospira. They are motile helical spirochaetes that metabolise long chain fatty acids as their carbon source. There are 17 species of Leptospira as determined by DNA-DNA hybridisation [1–4]. These species can be further divided into pathogenic, non-pathogenic and opportunistic/possibly pathogenic Leptospira with pathogenic species. The pathogenic Leptospira include; L. interrogans, L. kirschneri, L. santarosai, L. weilii, L. alexanderi, L. borgpetersenii, L. genomospecies 1 and L. noguchii. The non pathogenic Leptospira include: L. biflexa, L. meyeri, L. wolbachii, L. genomospecies 3, L. genomospecies 4, L. genomospecies 5 and opportunistic/intermediate pathogens Leptospira include L. broomi, L. fainei and L. inadai . The grouping of the last three species as opportunistic or possible pathogens is due to the lack of information on the pathogenicity of the species, different phenotypic characteristics compared to the pathogenic Leptospira and also the limited number of reports of these species involvement in human leptospirosis.
Before molecular techniques such as DNA-DNA hybridisation or 16s rRNA gene sequencing became available, speciation of the genus Leptospira was limited to the classifications of pathogenic (L. interrogans sensu lato) or saprophytic (L. biflexa sensu lato) and was performed using phenotypic tests such as growth at 13°C/30°C or growth in the presence of a chemical such as 8-Azaguanine . These tests can take up to 28 days to complete and the results can vary within a species . Since the introduction of molecular techniques, the identification of Leptospira species has generally been performed using either DNA-DNA hybridisation [2, 4] or 16s rRNA gene sequencing. Both of these methods have inherent disadvantages; DNA-DNA hybridisation is laborious and requires the use of radio-labelled isotopes  and the 16s rRNA gene has significant sequence homology between species which requires the majority of the gene to be sequenced for a definitive Leptospira identification. As an alternative target to 16s rRNA for species identification, the DNA Gyrase Subunit B gene (gyrB) has been successfully used for species identification in a wide variety of bacterial genera [7–13]. More recently the gyrB gene has been used for the identification of Leptospira borgpetersenii isolates from the Amami Islands  though this study used universal gyrB primers and only conducted limited phylogenetic analysis.
This paper reports the development of a molecular technique for the identification of pathogenic Leptospira species using conventional or real-time PCR amplification and sequencing of a partial fragment of the gyrB gene. The method was then used to ascertain gyrB sequences from representative reference strains of the eight pathogenic species. These sequences were used for phylogenetic and evolutionary comparisons between the species themselves and also between the gyrB gene and 16s rRNA gene. To highlight the potential value of the gyrB as an alternate identification gene, a blind trial was conducted between the two gene targets to identify previously uncharacterised clinical Leptospira isolates.
Results and Discussion
Using conventional and real-time PCR we were able to amplify a 504 bp product from the eight pathogenic Leptospira species: L. interrogans, L. borgpetersenii, L. weilii, L. santarosai, L. alexanderi, L. genomospecies 1, L. noguchii and L. kirschneri. No PCR products were amplified from the representatives of non-pathogenic species L. biflexa, L. meyeri or from the pathogenic/intermediate species such as L. inadai, L. fainei or L. broomi (Figure 1 and 2). The development of both conventional and real-time PCR methodologies for the amplification of the gyrB gene enables this method to be instituted at the majority of laboratories and is not dependant on having relatively expensive real-time PCR equipment. The advantage of using real-time PCR over conventional PCR is that it is quicker (amplification is completed in less than an hour) and there is no need to perform agarose gel electrophoresis or capture the gel image. Confirmation of the gyrB gene amplification was performed using the melting curve analysis on the LightCycler instrument. The Tm of the gyrB PCR product was found to be between 83.4°C and 84.8°C (Figure 2). Cycle sequencing of the gyrB PCR product and comparison of the DNA sequences enabled species specific identification. The gyrB DNA sequences of the reference strains were deposited on GenBank (Table 1).
The amplification of only the pathogenic species from the genus Leptospira has created an assay which has a wide potential in this field of research. For example, the gyrB PCR could be use to identify pathogenic Leptospira from cultures or identify Leptospira isolates that have been overgrown with bacteria or fungi. Additionally it would be possible to apply this method to clinical samples that contain high concentrations of Leptospira organisms such as kidney tissue. The non-culture identification of pathogenic Leptospira would be difficult without the use of specific gyrB primers as universal gyrB primers such as those used by Kawabata et al.  would amplify DNA from all bacteria present in a sample.
The potential of this assay must be balanced by three apparent limitations. Firstly the lack of sensitivity of conventional detection methodologies would not enable this test to be used in diagnosis of human infections where there are generally only low levels of Leptospira in the blood. Secondly, potentially pathogenic species such as L. fainei, L. inadai or L. broomi are not amplified, and therefore could be missed or excluded during molecular investigations. Finally, if a culture or sample contained two different Leptospira species then it would result in a mixed sequencing result, requiring use of DNA cloning and multiple sequencing reactions. The ultimate evolution of this method would be through the use of specific DNA probe detection either through Taqman/FRET probe in a real-time PCR guise or through a more conventional chemical-luminescent detection format.
Comparative phylogenetic analysis
Multiple alignments of the DNA sequences allowed phylogenetic comparisons between the species and between the gyrB gene and 16s rRNA gene to be performed. Global clustering (Figure 3) and similarity matrix analysis (Table 2 and 3) was performed using the aligned sequence data and the Unweighted Pair Group Method with Arithmetic Mean (UPGMA) algorithm for the 16s rRNA, gyrB and a consensus of both genes to examine the level of relatedness between the pathogenic species. 16s rRNA and gyrB genes show significant difference in total relatedness as shown by the percentage scale in Figure 3 and more accurately in the similarity matrix analysis (Table 2 and 3). The maximum nucleotide difference for the 16s rRNA gene ranges from 0.1% to 1.4% whilst for the gyrB gene it ranges from 3.5% to 16.1%. The differences in nucleotide divergence between the two genes is due to gyrB having a higher rate of base substitution (0.7–0.8% per 1 million years) when compared to 16s rRNA (1% per 50 million years) .
In addition to the cluster analysis, minimum evolution trees for 16s rRNA and gyrB were constructed using 1000 bootstrap replications (Figure 4). The trees have nearly identical topology except the evolution of L. kirschneri and L. interrogans in the gyrB gene is quite distinctly different to that of the evolution pattern in the 16s rRNA gene. The changes in topology signifies that the gyrB gene in pathogenic Leptospira has undergone significant evolutionary divergence compared to that of the 16s rRNA gene and this result is consisted with research conducted in other bacterial genera including the preliminary work conducted by Kawabata et al. with Leptospira [7–9, 11, 12, 14–18].
Identification of Leptospira clinical isolates using gyrB PCR
To validate the use of gyrB as an alternative target to 16s rRNA, a comparison study was performed using 50 unidentified clinical isolates all from human sources. The gyrB sequences of the unknown isolates were compared to those deposited on GenBank using a BLASTn search. Confirmation of the gyrB result was performed using 16s rRNA sequencing and BLAST analysis as described in the methods. There was 100% agreement between the 16s rRNA and gyrB gene sequencing for the identification of pathogenic Leptospira species. Within the 50 isolates tested there was found to be the following number of species; L. alexanderi (1 isolate), L. weilli (10 isolates), L. borgpetersenii (10 isolates) and L. interrogans (29 isolates). The advantage of using the gyrB gene over the 16s rRNA gene is that during the BLASTn searches on GenBank, the score values were generally higher and the E values generally lower for the predicted species when compared to the 16s rRNA gene BLASTn searches allowing for greater confidence in the final result.
We have developed and validated a conventional and real-time PCR method for the amplification of gyrB gene from pathogenic Leptospira. When compared to the 16s rRNA gene, the gyrB gene shows greater evolutionary divergence and an alternate evolutionary topology for L. kirschneri and L. interrogans using minimum evolution analysis. Additionally the greater divergence of the gyrB gene makes it more amendable to the identification of pathogenic Leptospira either through sequencing as shown in this study or in the future by Real-time PCR using DNA probe technology.
Leptospira strains and DNA extraction
In total, 37 reference strains from the eight pathogenic Leptospira species and 50 clinical Leptospira isolates from human sources were obtained from the WHO/FAO/OIE Collaborating Centre for Reference & Research on Leptospirosis, Brisbane, Australia. Genomic DNA was extracted by the following method: 500 μL of Ellinghausen McCullough Johnson Harris (EMJH) media containing actively growing Leptospira was centrifuged in a micro-centrifuge tube at 12,000 g for 5 min. The supernatant was removed and the pellet re-suspended in 400 μL of 1× TE Buffer (10 mM Tris, 1 mM EDTA, pH 8.0). This suspension was boiled for 10 min and then centrifuged at 12,000 g for 5 min.
gyrB amplification: Conventional PCR
PCR primers were developed from the two available Leptospira interrogans genome sequences: NC_005823[19, 20] and NC_004342 using Primer Premier 5.0 (Premier Biosoft) to amplify a 502 base pair (bp) fragment of the gyrB gene. PCR amplification was performed in a final volume of 25 μL using 1 × PCR buffer, 2.5 mM Magnesium Chloride (MgCl2), 200 μM dNTPs, 12.5 pmol of oligonucleotides; 2For and 504Rev (Table 4), one unit of AmpliTaq Gold, 2 μl of DNA extract and double distilled water (ddH2O) to make up the final volume. Thermal cycling was as following: Initial denaturation at 94°C for 10 min, followed by 35 cycles of 94°C for 30 s, 60°C for 30 s and 72°C for 30 s with a final extension at 72°C for 10 min. 5 μL of PCR product was electrophoresis on 1.5% agarose gel at 80 V for 60 min (Figure 1).
gyrB amplification: Real-time PCR
Real-time amplification of the gyrB gene was performed in a total volume of 20 μL containing 2 μL of Fast-Start Sybr green mix (Roche), 2.4 μL of 25 mM MgCl2, 10 pmol of oligonucleotides; 2For and 504Rev, 2 μL of template DNA and ddH2O to make up the final volume. Thermal cycling was performed on a LightCycler real-time thermalcycler (Roche) using the following program: initial denaturation at 95°C for 10 min and 40 cycles of 95°C for 10 s, 60°C for 20 s and 72°C for 20 s. Fluorescence readings were taken at the end of each extension cycle in the F1 (FAM/Sybr green) channel. Melting curve analysis was performed by heating the PCR product from 60°C to 95°C and monitoring the fluorescence change every 0.2°C. The melting temperature or Tm was calculated by calculated on the initial fluorescence curve (F1/I) by plotting the negative derivative of fluorescence over temperature versus temperature (-dF1/dT versus T). Amplified products were removed before sequencing from the capillaries by uncapping and inverting the capillary in a micro-centrifuge tube and centrifuging at 1,500 g for 10 s.
16s rRNA amplification
16s rRNA amplification was performed in a final volume of 25 μL containing 1× PCR buffer, 2.0 mM of MgCl2, 200 μM dNTPs, 10.0 pmol of oligonucleotides; FD1MOD and 13R (Table 4), one unit of AmpliTaq Gold, 2 μl of DNA extract and (ddH2O) to make up the final volume. Thermal cycling was as following: Initial denaturation at 94°C for 10 min, followed by 35 cycles of 94°C for 30 s, 55°C for 30 s and 72°C for 30 s with a final extension at 72°C for 10 min. Agarose electrophoresis was performed as above. The 1382 bp product were sequenced as described below in both the forward and reverse directions using the original primers and with internal primers; 515F, 91e, 11e, 16s1RRB, 907R and 342R (Table 4).
Excess primers and dNTP's were removed from the remaining PCR product using the following enzymatic method: 2.5 μl of 10× Antarctic phosphatase buffer (New England Biolabs, NEB), 10 units of Exonuclease I, E. coli (Fermentas), 2.5 units of Antarctic phosphatase (NEB) and 1.5 μl of ddH2O were added to each sample. The PCR product plus enzyme mix were incubated at 37°C for 45 min followed by 85°C for 15 min to inactivate the enzyme. DNA sequencing was performed using the Big Dye Terminator (BDT) sequencing version 3.1 (Applied Biosystems) with the following modifications: each 20 μL reaction contained 0.5 μL of BDT mix (1/16th dilution in final volume), 3.75 μL of 5× dilution buffer, 3.2 pmol of primer, 5–10 ng of DNA and ddH2O to make up the final volume. Cycle sequencing was performed using 33 cycles of 95°C for 10 s, 50°C for 10 s and 60°C for 4 min. The cycle sequencing products were purified using the sodium acetate/alcohol precipitation method as per manufacturers' instructions (Applied Biosystems). The purified products were forwarded to the Griffith university DNA sequencing facility (GUDSF), Brisbane, Australia for capillary electrophoresis using the ABI 3130 × l instrument. The sequences were assembled and trimmed to a minimum of two contiguous sequences using the Vector NTI software (Invitrogen).
Analysis of the DNA sequences was performed using the Bionumerics (Applied maths) and MEGA version 3.1 software packages . 37 reference sequences for both gyrB and 16s rRNA were deposited on GenBank (Table 1). Unknown isolates were submitted to a nucleotide Basic Local Alignment Search Tool (BLASTn)  search available at the National Centre for Biotechnology Information (NCBI) website .
Levett PN, Morey RE, Galloway RL, Steigerwalt AG: Leptospira broomii sp. nov., isolated from humans with leptospirosis. Int J Syst Evol Microbiol. 2006, 56 (Pt 3): 671-673. 10.1099/ijs.0.63783-0.
Brenner DJ, Kaufmann AF, Sulzer KR, Steigerwalt AG, Rogers FC, Weyant RS: Further determination of DNA relatedness between serogroups and serovars in the family Leptospiraceae with a proposal for Leptospira alexanderi sp. nov. and four new Leptospira genomospecies. Int J Syst Bacteriol. 1999, 49 Pt 2: 839-858.
Perolat P, Chappel RJ, Adler B, Baranton G, Bulach DM, Billinghurst ML, Letocart M, Merien F, Serrano MS: Leptospira fainei sp. nov., isolated from pigs in Australia. Int J Syst Bacteriol. 1998, 48 Pt 3: 851-858.
Yasuda PH, Steigerwalt AG, Sulzer CR, Kaufmann AF, Rogers FC, Brenner DJ: Deoxyribonucleic acid relatedness between serogroups and serovars in the family Leptospiraceae with proposals for seven new Leptospira species. Int J Syst Bacteriol. 1987, 37: 407-415.
Johnson RC, Rogers P: Differentiation of Pathogenic and Saprophytic Leptospires with 8-Azaguanine. J Bacteriol. 1964, 88: 1618-1623.
Cho JC, Tiedje JM: Bacterial species determination from DNA-DNA hybridization by using genome fragments and DNA microarrays. Appl Environ Microbiol. 2001, 67 (8): 3677-3682. 10.1128/AEM.67.8.3677-3682.2001.
Le Roux F, Gay M, Lambert C, Nicolas JL, Gouy M, Berthe F: Phylogenetic study and identification of Vibrio splendidus-related strains based on gyrB gene sequences. Dis Aquat Organ. 2004, 58 (2-3): 143-150.
Yanez MA, Catalan V, Apraiz D, Figueras MJ, Martinez-Murcia AJ: Phylogenetic analysis of members of the genus Aeromonas based on gyrB gene sequences. Int J Syst Evol Microbiol. 2003, 53 (Pt 3): 875-883. 10.1099/ijs.0.02443-0.
Coenye T, LiPuma JJ: Use of the gyrB gene for the identification of Pandoraea species. FEMS Microbiol Lett. 2002, 208 (1): 15-19. 10.1111/j.1574-6968.2002.tb11053.x.
Coenye T, Vanlaere E, LiPuma JJ, Vandamme P: Identification of genomic groups in the genus Stenotrophomonas using gyrB RFLP analysis. FEMS Immunol Med Microbiol. 2004, 40 (3): 181-185. 10.1016/S0928-8244(03)00307-9.
Itoh Y, Kawamura Y, Kasai H, Shah MM, Nhung PH, Yamada M, Sun X, Koyana T, Hayashi M, Ohkusu K, Ezaki T: dnaJ and gyrB gene sequence relationship among species and strains of genus Streptococcus. Syst Appl Microbiol. 2006
Kasai H, Ezaki T, Harayama S: Differentiation of phylogenetically related slowly growing mycobacteria by their gyrB sequences. J Clin Microbiol. 2000, 38 (1): 301-308.
Delmas J, Breysse F, Devulder G, Flandrois JP, Chomarat M: Rapid identification of Enterobacteriaceae by sequencing DNA gyrase subunit B encoding gene. Diagn Microbiol Infect Dis. 2006
Kawabata H, Sakakibara S, Imai Y, Masuzawa T, Fujita H, Tsurumi M, Sato F, Takano A, Nogami S, Kaneda K, Watanabe H: First record of Leptospira borgpetersenii isolation in the Amami Islands, Japan. Microbiol Immunol. 2006, 50 (6): 429-434.
Yamamoto S, Harayama S: Phylogenetic analysis of Acinetobacter strains based on the nucleotide sequences of gyrB genes and on the amino acid sequences of their products. Int J Syst Bacteriol. 1996, 46 (2): 506-511.
Yamamoto S, Bouvet PJ, Harayama S: Phylogenetic structures of the genus Acinetobacter based on gyrB sequences: comparison with the grouping by DNA-DNA hybridization. Int J Syst Bacteriol. 1999, 49 Pt 1: 87-95.
Schwan TG, Raffel SJ, Schrumpf ME, Policastro PF, Rawlings JA, Lane RS, Breitschwerdt EB, Porcella SF: Phylogenetic analysis of the spirochetes Borrelia parkeri and Borrelia turicatae and the potential for tick-borne relapsing fever in Florida. J Clin Microbiol. 2005, 43 (8): 3851-3859. 10.1128/JCM.43.8.3851-3859.2005.
Maeda Y, Shinohara H, Kiba A, Ohnishi K, Furuya N, Kawamura Y, Ezaki T, Vandamme P, Tsushima S, Hikichi Y: Phylogenetic study and multiplex PCR-based detection of Burkholderia plantarii, Burkholderia glumae and Burkholderia gladioli using gyrB and rpoD sequences. Int J Syst Evol Microbiol. 2006, 56 (Pt 5): 1031-1038. 10.1099/ijs.0.64184-0.
Nascimento AL, Verjovski-Almeida S, Van Sluys MA, Monteiro-Vitorello CB, Camargo LE, Digiampietri LA, Harstkeerl RA, Ho PL, Marques MV, Oliveira MC, Setubal JC, Haake DA, Martins EA: Genome features of Leptospira interrogans serovar Copenhageni. Braz J Med Biol Res. 2004, 37 (4): 459-477. 10.1590/S0100-879X2004000400003.
Nascimento AL, Ko AI, Martins EA, Monteiro-Vitorello CB, Ho PL, Haake DA, Verjovski-Almeida S, Hartskeerl RA, Marques MV, Oliveira MC, Menck CF, Leite LC, Carrer H, Coutinho LL, Degrave WM, Dellagostin OA, El-Dorry H, Ferro ES, Ferro MI, Furlan LR, Gamberini M, Giglioti EA, Goes-Neto A, Goldman GH, Goldman MH, Harakava R, Jeronimo SM, Junqueira-de-Azevedo IL, Kimura ET, Kuramae EE, Lemos EG, Lemos MV, Marino CL, Nunes LR, de Oliveira RC, Pereira GG, Reis MS, Schriefer A, Siqueira WJ, Sommer P, Tsai SM, Simpson AJ, Ferro JA, Camargo LE, Kitajima JP, Setubal JC, Van Sluys MA: Comparative genomics of two Leptospira interrogans serovars reveals novel insights into physiology and pathogenesis. J Bacteriol. 2004, 186 (7): 2164-2172. 10.1128/JB.186.7.2164-2172.2004.
Ren SX, Fu G, Jiang XG, Zeng R, Miao YG, Xu H, Zhang YX, Xiong H, Lu G, Lu LF, Jiang HQ, Jia J, Tu YF, Jiang JX, Gu WY, Zhang YQ, Cai Z, Sheng HH, Yin HF, Zhang Y, Zhu GF, Wan M, Huang HL, Qian Z, Wang SY, Ma W, Yao ZJ, Shen Y, Qiang BQ, Xia QC, Guo XK, Danchin A, Saint Girons I, Somerville RL, Wen YM, Shi MH, Chen Z, Xu JG, Zhao GP: Unique physiological and pathogenic features of Leptospira interrogans revealed by whole-genome sequencing. Nature. 2003, 422 (6934): 888-893. 10.1038/nature01597.
Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5 (2): 150-163. 10.1093/bib/5.2.150.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410. 10.1006/jmbi.1990.9999.
Kotilainen P, Jalava J, Meurman O, Lehtonen OP, Rintala E, Seppala OP, Eerola E, Nikkari S: Diagnosis of meningococcal meningitis by broad-range bacterial PCR with cerebrospinal fluid. J Clin Microbiol. 1998, 36 (8): 2205-2209.
Relman DA, Schmidt TM, MacDermott RP, Falkow S: Identification of the uncultured bacillus of Whipple's disease. N Engl J Med. 1992, 327 (5): 293-301.
Relman DA, Loutit JS, Schmidt TM, Falkow S, Tompkins LS: The agent of bacillary angiomatosis. An approach to the identification of uncultured pathogens. N Engl J Med. 1990, 323 (23): 1573-1580.
Wilbrink B, van der Heijden IM, Schouls LM, van Embden JD, Hazes JM, Breedveld FC, Tak PP: Detection of bacterial DNA in joint samples from patients with undifferentiated arthritis and reactive arthritis, using polymerase chain reaction with universal 16S ribosomal RNA primers. Arthritis Rheum. 1998, 41 (3): 535-543. 10.1002/1529-0131(199803)41:3<535::AID-ART20>3.0.CO;2-4.
Lane DJ: 16S/23S rRNA sequencing. Nucleic Acid Techniques in Bacterial Systematics. Edited by: Stackebrandt E, Goodfellow M. 1991, Chichester , John Wiley and Sons
The authors wish to acknowledge Queensland Health for providing funding and for their on-going support of the WHO/FAO/OIE Collaborating Centre for Reference & Research on Leptospirosis.
The author(s) declare that they have no competing interests.
AS was responsible for design of the study, conducting the molecular experiments and the preparation of the manuscript. MD and MS provided laboratory support by providing culture, maintaining culture collections and contributed to the editing of the manuscript. LS approved the research study/funding, provided intellectual input and contributed to the editing of the manuscript. All authors have read and approved the final manuscript.