- Research article
- Open Access
Corynebacterium ulcerans 0102 carries the gene encoding diphtheria toxin on a prophage different from the C. diphtheriae NCTC 13129 prophage
© Sekizuka et al.; licensee BioMed Central Ltd. 2012
Received: 10 January 2012
Accepted: 14 May 2012
Published: 14 May 2012
Corynebacterium ulcerans can cause a diphtheria-like illness, especially when the bacterium is lysogenized with a tox gene-carrying bacteriophage that produces diphtheria toxin. Acquisition of toxigenicity upon phage lysogenization is a common feature of C. ulcerans and C. diphtheriae. However, because of a lack of C. ulcerans genome information, a detailed comparison of prophages has not been possible between these two clinically important and closely related bacterial species.
We determined the whole genome sequence of the toxigenic C. ulcerans 0102 isolated in Japan. The genomic sequence showed a striking similarity with that of Corynebacterium pseudotuberculosis and, to a lesser extent, with that of C. diphtheriae. The 0102 genome contained three distinct prophages. One of these, ΦCULC0102-I, was a tox-positive prophage containing genes in the same structural order as for tox-positive C. diphtheriae prophages. However, the primary structures of the individual genes involved in the phage machinery showed little homology between the two counterparts.
Taken together, these results suggest that the tox-positive prophage in this strain of C. ulcerans has a distinct origin from that of C. diphtheriae NCTC 13129.
A diphtheria-like infectious disease caused by Corynebacterium ulcerans is increasing in clinical importance in developed countries and is now regarded as “diphtheria” in Europe [1, 2]. Infection with C. ulcerans occurs in a wide range of hosts, including cats, dogs, pigs, cows, and whales [3–9]. The first clearly documented case of zoonotic transmission involved a dog, as reported by Lartigue et al. . This is in contrast to the causative agent of classical diphtheria, C. diphtheriae, whose host species is thought to be limited to humans . Nevertheless, the two species share a common feature: upon lysogenization of tox-encoding bacteriophages, they become toxigenic and are able to produce the potent diphtheria toxin [1, 10]. This toxin is known to contribute to disease progression, occasionally leading to death. It is encoded by a single gene designated tox, situated inside prophages lysogenized in the bacterial genome of C. diphtheriae. The prophages are capable of induction, by ultraviolet light or DNA-damaging agents such as mitomycin C, and yield β-, δ-, ω- and other functional bacteriophage particles . Some types of bacteriophages can infect both C. diphtheriae and C. ulcerans[13–16]. Furthermore, the C. ulcerans tox gene is also encoded in a genome region surrounded by phage attachment (att) sites conserved between the two species [7, 16]. The nucleotide sequences of C. ulcerans tox genes were published by Sing et al. They showed some diversity in the genetic sequence among C. ulcerans strains, in contrast to the highly conserved C. diphtheriae tox gene [17, 18].
In 2003, the nucleotide sequence of the whole genome of C. diphtheriae strain NCTC13129 was reported . The sequence information revealed some striking features of the bacterial genome, such as the presence of as many as 13 pathogenicity islands (PAIs) , uncommon among C. diphtheriae strains . The presence of a tox-positive prophage flanked by the att regions was confirmed and supported the findings of previous reports . Despite comparable clinical importance, the genomic sequence of toxigenic C. ulcerans has not yet been reported. In the present study, we determined the nucleotide sequence of the toxigenic C. ulcerans isolate 0102 genome, obtained in 2001 from the pharyngeal pseudomembrane of a 52-year-old woman presenting with a sore throat and fever. This was the first toxigenic C. ulcerans infection reported in Japan. This patient had been living with nearly 20 cats before the onset of illness . Details of the bacteriological characteristics of the isolate have been described elsewhere . Our analysis was especially directed towards the structure of the tox-positive prophage because of its unexpectedly novel structure.
Genome sequence and genomic information for C. ulcerans 0102
To determine the complete genome sequence of C. ulcerans 0102, obtained short reads were assembled into five contigs by de novo assembly. Each gap was filled by direct PCR and sequencing. A circular chromosome sequence of C. ulcerans 0102 represents 2,579,188 bp, with a G + C content of 53.4% (Additional file 1) and corresponds to the predicted restriction fragment profiles obtained by PFGE analysis (Additional file 2). The chromosome possesses 2,349 coding sequences, 51 tRNA genes, and 4 rrn rRNA operons.
Comparative genome analysis of three pathogenic Corynebacterium spp
The tox-positive prophage of C. ulcerans 0102
The two tox- positive prophages share the same structural features, with genes aligned in an ‘integrase - packaging - head - tail - lysis - toxin’ orientation (Figure 2). Pair-wise alignment of the prophages indicates a high similarity in the region encoding the putative integrase, the 3′-ends of CULC0102_0211 and CULC0102_0212, tox, and the attachment sites (Figure 2). The major phage machineries encoded in the internal phage region showed low similarity at the nucleotide and amino acid levels (less than 18%) between C. ulcerans 0102 and C. diphtheriae NCTC13129.
Whole-genome sequencing has revealed that the C. ulcerans 0102 genome is composed of 2,579,188 bp with a G + C content of 53.4%. These values are similar to those recently reported for C. ulcerans strains 809 (2,502,095 bp, 53.3% G + C) and BR-AD22 (2,606,374 bp, 53.4% G + C) . C. ulcerans 0102 shares many common features with the two previously reported strains, including 12 virulence factors. Strain 0102 is distinctive with respect to the features of prophages integrated in its genome. It possesses a unique tox-positive prophage, ΦCULC0102-I, in its chromosome (Figure 1 and Additional file 1). In the same position of the recently reported C. ulcerans 809 genome exists a remnant phage-related integrase (intC) gene  (Figure 2). The C. ulcerans 0102 prophage differs from the corresponding prophage in C. diphtheriae. Although the integrase and tox gene sequences of ΦCULC0102-I showed high similarity to those of the corynephage encoding tox in C. diphtheriae NCTC 13129, the major phage machinery genes in ΦCULC0102-I are distinct from those in other corynephages in C. diphtheriae (Figure 2). This suggests that C. ulcerans 0102 did not immediately acquire the C. diphtheriae tox-positive corynephage.
The C. diphtheriae tox gene is highly conserved among temporally and geographically diverse strains , therefore greater variation in tox genes from C. ulcerans isolates suggests that this strain might have acquired the tox gene before C. diphtheriae.
In a recent report, whole genome sequence analysis of non-toxigenic C. ulcerans 809 and BR-AD22 , the β-corynephage-like truncated integrases (CULC809_00176 and CULC22_00173) are located adjacent to the tRNAArg gene, similar to ΦCULC0102-I in C. ulcerans 0102 and C. diphtheriae. The tRNAArg gene (CULC0102_t08) appears to be a ‘hotspot’ for the acquisition of ΦCULC0102-I-like prophages by homologous integrase.
The whole genome sequences of C. ulcerans 809 and BR-AD22 contain possible virulence factors, such as corynebacterial protease (CP40), phospholipase D (Pld), neuraminidase (NanH), venom serine protease (Vsp1), trypsin-like serine protease (TspA), Rpf interacting protein (RpfI), cell wall-associated hydrolase (CwlH), and five surface-anchored proteins (SpaB–F) . The SpaA-type pilin, encoded by the spaABC srtA gene cluster, is considered to play a crucial role in adhesion of C. diphtheriae. The gene encoding the shaft protein of SpaA-type pilin (spaA) was absent in C. ulcerans 0102, a feature consistent with previous findings in C. ulcerans 809 and BR-AD2 . As SpaB and SpaC proteins, which are assumed to be present in all three C. ulcerans strains, can contribute to host-cell adhesion in the absence of SpaA , this may imply a common mechanism of cell adhesion by C. ulcerans.
The C. ulcerans 809 strain was isolated from a patient with a rapid fatal pulmonary infection. The 809 strain-unique virulence factor (shiga toxin-like ribosome-binding protein, Rbp) is located adjacent to the truncated integrase (CULC809_00176) and corresponds to the integrase of ΦCULC0102-I. It appears that virulence factors have been acquired as a cassette gene in the ΦCULC0102-I-like prophage. It is intriguing to note that the 0102 strain does not carry the 809 strain-unique virulence factors (Rbp and the additional venom serine protease, Vsp2), but instead carries the tox gene on ΦCULC0102-I, which resulted in a diphtheria-like illness in a 52-year-old woman.
Isolates of C. ulcerans are generally obtained from a diverse range of animals, including humans. Isolation of a human pathogen C. diphtheriae from animals has been reported previously, although it is rare . The tox gene might be frequently transmitted through common prophages with the aid of the highly homologous regions among Corynebacterium spp., including C. diphtheriae and C. ulcerans isolated from animal sources.
Toxigenic C. ulcerans is an emerging pathogen that can be transmitted from animals to humans . In the host organism, as well as in C. diphtheriae, the tox gene  is encoded by prophages. Through genome sequencing, we have identified a novel structure in a tox-positive C. ulcerans prophage with no significant sequence homology to those in C. diphtheriae. This suggests distinct origins of the prophages and thus may also explain the difference in the primary structures of their tox genes. The tox-positive bacteriophages may increase the dissemination risk of toxigenic C. ulcerans isolates, therefore, C. ulcerans isolates from both human and animal sources should be investigated further to determine the level of variation.
This research was not carried out on humans. No experimental research on animals was carried out.
Preparation of genomic DNA
Genomic DNA was isolated by conventional methods, using phenol extraction and ethanol precipitation from heat-killed bacterial cells propagated in brain-heart infusion liquid medium.
Short-read DNA sequencing using an Illumina Genome Analyzer IIx
DNA libraries of the ~600 bp insert length of C. ulcerans 0102 were prepared using a genomic DNA Sample Prep Kit (Illumina, San Diego, CA, USA). DNA clusters were generated on a slide using a Cluster Generation Kit (ver. 4) on an Illumina Cluster Station (Illumina), according to the manufacturer’s instructions. Sequencing runs for 80-mer short reads were performed using an Illumina Genome Analyzer IIx (GA IIx) and TruSeq SBS kit v5. Fluorescent images were analyzed using the Illumina base-calling pipeline RTA2.6/SCS2.8 to obtain FASTQ-formatted sequence data.
De novo assembly of short DNA reads and gap-closing
The 80-mer reads were assembled (parameters k64, n51, c32.1373) using ABySS-pe v1.2.0 . Predicted gaps were amplified with a specific PCR primer pair, followed by Sanger DNA sequencing using a BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA).
Validation of the complete genome sequence using short-read mapping and pulsed-field gel electrophoresis (PFGE)
To validate the genome sequence, 40–mer short reads were re-aligned with the sequence using Maq software (ver. 0.7.1) and the easyrun Perl-command . Read alignment was inspected using the MapView graphical alignment viewer . PFGE analysis was performed to validate the predicted restriction fragment profiles from the complete genome sequence, according to De Zoysa et al. . Bacterial cells were lysed with lysozyme and protease , embedded in plugs, digested with the restriction endonuclease Sfi I (New England Biolabs, Ipswitch, MA, USA) and electrophoresed in a CHEF DRII apparatus (Bio-Rad, Hercules, CA, USA) at 11°C with a pulse time of 5–20 s for the first 20 h and 1–5 s for the following 18 h.
Annotation and pair-wise alignment analysis
Gene prediction from the complete sequence was performed using the NCBI Prokaryotic Genomes Automatic Annotation Pipeline (PGAAP; http://www.ncbi.nlm.nih.gov/genomes/static/pipeline.html). Several of the suggested errors were revised manually. Pseudogenes that were identified by PGAAP were checked using the read-mapping correction described above. Genomic information, such as nucleic acid variations and circular representation, was analyzed using IMC-GE software (Insilicobiology, Yokohama, Japan). A BLASTN homology search  was performed for the whole chromosome sequences of C. pseudotuberculosis FRC41 (accession no. NC_014329), C. ulcerans 0102, and C. diphtheriae NCTC 13129 (accession no. NC_002935). Aligned images of the homologous regions were visualized with the ACT program .
Phylogenetic analyses of all nucleotide sequences were conducted using the neighbor-joining method with 1,000-times bootstrapping in ClustalW2 . FigTree ver. 1.3.1 (http://tree.bio.ed.ac.uk/software/figtree/) software was used to display the generated tree.
Nucleotide sequence accession numbers
The complete chromosome sequence for the C. ulcerans 0102 strain has been deposited in the DNA Data Bank of Japan (DDBJ; accession no. AP012284).
The authors are grateful to Akio Hatanaka, Atsuhiro Tsunoda and Kenji Ooe for the 0102 clinical isolate. This work was supported by grants for Research on Emerging and Re-emerging Infectious Diseases (H23 Shinko-Ippan-007 and H22-Shinko-Ippan-010), from the Ministry of Health, Labour and Welfare, Japan.
- Bonnet JM, Begg NT: Control of diphtheria: guidance for consultants in communicable disease control. Commun Dis Public Health. 1999, 2: 242-249.PubMedGoogle Scholar
- European Centre for Disease Prevention and Control: Diphtheria. Surveillance Report: Annual epidemiological report on communicable diseases in Europe 2010. 2010, 133-135.Google Scholar
- Dias AASO, Silva FC, Pereira GA, Souza MC, Camello TCF, Damasceno JALD, Pacheco LGC, Miyoshi A, Azevedo VA, Hirata R, et al.: Corynebacterium ulcerans isolated from an asymptomatic dog kept in an animal shelter in the metropolitan area of Rio de Janeiro, Brazil. Vector Borne Zoonotic Dis. 2010, 10: 743-748. 10.1089/vbz.2009.0132.PubMedView ArticleGoogle Scholar
- Katsukawa C, Kawahara R, Inoue K, Ishii A, Yamagishi H, Kida K, Nishino S, Nagahama S, Komiya T, Iwaki M, Takahashi M: Toxigenic Corynebacterium ulcerans Isolated from the domestic dog for the first time in Japan. Jpn J Infect Dis. 2009, 62: 171-172.PubMedGoogle Scholar
- Lartigue M-F, Monnet X, Le Flèche A, Grimont PAD, Benet J-J, Durrbach A, Fabre M, Nordmann P: Corynebacterium ulcerans in an immunocompromised patient with diphtheria and her dog. J Clin Microbiol. 2005, 43: 999-1001. 10.1128/JCM.43.2.999-1001.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Schuhegger R, Schoerner C, Dlugaiczyk J, Lichtenfeld I, Trouillier A, Zeller-Peronnet V, Busch U, Berger A, Kugler R, Hörmansdorfer S, Sing A: Pigs as source for toxigenic Corynebacterium ulcerans. Emerg Infect Dis. 2009, 15: 1314-1315. 10.3201/eid1508.081568.PubMedPubMed CentralView ArticleGoogle Scholar
- Seto Y, Komiya T, Iwaki M, Kohda T, Mukamoto M, Takahashi M, Kozaki S: Properties of corynephage attachment site and molecular epidemiology of Corynebcterium ulcerans isolated from humans and animals in Japan. Jpn J Infect Dis. 2008, 61: 116-122.PubMedGoogle Scholar
- De Zoysa A, Hawkey PM, Engler K, George R, Mann G, Reilly W, Taylor D, Efstratiou A: Characterization of toxigenic Corynebacterium ulcerans strains isolated from humans and domestic cats in the United Kingdom. J Clin Microbiol. 2005, 43: 4377-10.1128/JCM.43.9.4377-4381.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Yoshimura Y, Yamamoto A, Komiya T: A case of axillary lymph node abscess caused by percutaneous infection of Corynebacterium ulcerans through scratch by a pus-discharging cat, June 2010 (in Japanese). Infect Agents Surveillance Rep. 2010, 31: 331-Google Scholar
- Murphy JR: Chapter 32 Corynebacterium diphtheriae. Medical Microbiology. Edited by: Baron S. 1996, University of Texas Medical Branch at Galveston, Galveston, 4Google Scholar
- Pappenheimer AM, Gill DM: Diphtheria. Recent studies have clarified the molecular mechanisms involved in its pathogenesis. Science. 1973, 182: 353-358. 10.1126/science.182.4110.353.PubMedView ArticleGoogle Scholar
- Rappuoli R, Michel JL, Murphy JR: Integration of corynebacteriophages: tox+, xtox+ and gtox+ into two attachment sites on the Corynebacterium diphtheriae chromosome. J Bacteriol. 1983, 153: 1202-1210.PubMedPubMed CentralGoogle Scholar
- Ishii-Kanei C, Uchida T, Yoneda M: Isolation of a cured strain from Corynebacterium diphtheriae PW8. Infect Immun. 1979, 25: 1081-1083.PubMedPubMed CentralGoogle Scholar
- Cianciotto NP, Groman NB: Extended host range of a β-related corynebacteriophage. FEMS Microbiol Lett. 1996, 140: 221-225.PubMedGoogle Scholar
- Oram M, Woolston JE, Jacobson AD, Holmes RK, Oram DM: Bacteriophage-based vectors for site-specific insertion of DNA in the chromosome of Corynebacteria. Gene. 2007, 391: 53-62. 10.1016/j.gene.2006.12.003.PubMedPubMed CentralView ArticleGoogle Scholar
- Cianciotto N, Rappuoli R, Groman N: Detection of homology to the beta bacteriophage integration site in a wide variety of Corynebacterium spp. J Bacrteriol. 1986, 168: 103-108.Google Scholar
- Sing A, Bierschenk S, Heesemann J: Classical diphtheria caused by Corynebacterium ulcerans in Germany: amino acid sequence differences between diphtheria toxins from Corynebacterium diphtheriae and C. ulcerans. Clin Infect Dis. 2005, 40: 325-326. 10.1086/426687.PubMedView ArticleGoogle Scholar
- Sing A, Hogardt M, Bierschenk S, Heesemann J: Detection of differences in the nucleotide and amino acid sequences of diphtheria toxin from Corynebacterium diphtheriae and Corynebacterium ulcerans causing extrapharyngeal infections. J Clin Microbiol. 2003, 41: 4848-4851. 10.1128/JCM.41.10.4848-4851.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Cerdeño-Tárraga A-M, Efstratiou A, Dover LG, Holden MTG, Pallen M, Bentley SD, Besra GS, Churcher C, James KD, De Zoysa A, et al.: The complete genome sequence and analysis of Corynebacterium diphtheriae NCTC13129. Nucl Acids Res. 2003, 31: 6516-6523. 10.1093/nar/gkg874.PubMedPubMed CentralView ArticleGoogle Scholar
- Iwaki M, Komiya T, Yamamoto A, Ishiwa A, Nagata N, Arakawa Y, Takahashi M: Genome organization and pathogenicity of Corynebacterium diphtheriae C7(−) and PW8 strains. Infect Immun. 2010, 78: 3791-3800. 10.1128/IAI.00049-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Cianciotto N, Serwold-Davis T, Groman N, Ratti G, Rappuoli R: DNA sequence homology between attB-related sites of Corynebacterium diphtheriae, Corynebacterium ulcerans, Corynebacterium glutamicum, and the attP site of gamma-corynephage. FEMS Microbiol Lett. 1990, 66: 299-301.Google Scholar
- Hatanaka A, Tsunoda A, Okamoto M, Ooe K, Nakamura A, Miyakoshi M, Komiya T, Takahashi M: Corynebacterium ulcerans diphtheria in Japan. Emerg Infect Dis. 2003, 9: 752-753.PubMedPubMed CentralView ArticleGoogle Scholar
- Komiya T, Seto Y, De Zoysa A, Iwaki M, Hatanaka A, Tsunoda A, Arakawa Y, Kozaki S, Takahashi M: Two Japanese Corynebacterium ulcerans isolates from the same hospital: ribotype, toxigenicity and serum antitoxin titre. J Med Microbiol. 2010, 59: 1497-1504. 10.1099/jmm.0.022491-0.PubMedView ArticleGoogle Scholar
- Trost E, Al-Dilaimi A, Papavasiliou P, Schneider J, Viehoever P, Burkovski A, Soares SC, Almeida SS, Dorella FA, Miyoshi A, et al.: Comparative analysis of two complete Corynebacterium ulcerans genomes and detection of candidate virulence factors. BMC Genomics. 2011, 12: 383-10.1186/1471-2164-12-383.PubMedPubMed CentralView ArticleGoogle Scholar
- Brüssow H, Canchaya C, Hardt W-D: Phages and the evolution of bacterial pathogens: from genomic rearrangements to lysogenic conversion. Microbiol Mol Biol Rev. 2004, 68: 560-602. 10.1128/MMBR.68.3.560-602.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Ochman H, Lawrence JG, Groisman EA: Lateral gene transfer and the nature of bacterial innovation. Nature. 2000, 405: 299-304. 10.1038/35012500.PubMedView ArticleGoogle Scholar
- Liu Y, Harrison PM, Kunin V, Gerstein M: Comprehensive analysis of pseudogenes in prokaryotes: widespread gene decay and failure of putative horizontally transferred genes. Genome Biol. 2004, 5: r64-10.1186/gb-2004-5-9-r64.PubMedPubMed CentralView ArticleGoogle Scholar
- Katsukawa C, Komiya T, Yamagishi H, Ishii A, Nishino S, Nagahama S, Iwaki M, Yamamoto A, Takahashi M: Prevalence of Corynebacterium ulcerans in dogs in Osaka, Japan. J Med Microbiol. 2012, 61: 266-273. 10.1099/jmm.0.034868-0.PubMedView ArticleGoogle Scholar
- Nakao H, Mazurova IK, Glushkevich T, Popovic T: Analysis of heterogeneity of Corynebacterium diphtheriae toxin gene, tox, and its regulatory element, dtxR, by direct sequencing. Res Microbiol. 1997, 148: 45-54. 10.1016/S0923-2508(97)81899-2.PubMedView ArticleGoogle Scholar
- Mandlik A, Swierczynski A, Das A, Ton-That H: Corynebacterium diphtheriae employs specific minor pilins to target human pharyngeal epithelial cells. Mol Microbiol. 2007, 64: 111-124. 10.1111/j.1365-2958.2007.05630.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Hall AJ, Cassiday PK, Bernard KA, Bolt F, Steigerwalt AG, Bixler D, Pawloski LC, Whitney AM, Iwaki M, Baldwin A, et al.: Novel Corynebacterium diphtheriae in domestic cats. Emerg Infect Dis. 2010, 16: 688-691.PubMedPubMed CentralView ArticleGoogle Scholar
- Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJM, Birol İ: ABySS: a parallel assembler for short read sequence data. Genome Res. 2009, 19: 1117-1123. 10.1101/gr.089532.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Ruan J, Durbin R: Mapping short DNA sequencing reads and calling variants using mapping quality scores. Genome Res. 2008, 18: 1851-1858. 10.1101/gr.078212.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Bao H, Guo H, Wang J, Zhou R, Lu X, Shi S: MapView: visualization of short reads alignment on a desktop computer. Bioinformatics. 2009, 25: 1554-1555. 10.1093/bioinformatics/btp255.PubMedView ArticleGoogle Scholar
- De Zoysa A, Efstratiou A, George RC, Jahkola M, Vuopio-Varkila J, Deshevoi S, Tseneva GY, Rikushin Y: Molecular epidemiology of Corynebacterium diphtheriae from northwestern Russia and surrounding countries studied by using ribotyping and pulsed-field gel electrophoresis. J Clin Microbiol. 1995, 33: 1080-1083.PubMedPubMed CentralGoogle Scholar
- Murrey BE, Singh KV, Heath JD, Sharma BR, Weinstock GM: Comparison of genomic DNAs of different enterococcal isolates using restriction endonucleases with infrequent recognition sites. J Clin Microbiol. 1990, 28: 2059-2063.Google Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView ArticleGoogle Scholar
- Carver T, Berriman M, Tivey A, Patel C, Böhme U, Barrell BG, Parkhill J, Rajandream M-A: Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. Bioinformatics. 2008, 24: 2672-2676. 10.1093/bioinformatics/btn529.PubMedPubMed CentralView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al.: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.PubMedView ArticleGoogle Scholar
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.