Virulence factor rtx in Legionella pneumophila, evidence suggesting it is a modular multifunctional protein
© D'Auria et al; licensee BioMed Central Ltd. 2008
Received: 25 July 2007
Accepted: 14 January 2008
Published: 14 January 2008
The repeats in toxin (Rtx) are an important pathogenicity factor involved in host cells invasion of Legionella pneumophila and other pathogenic bacteria. Its role in escaping the host immune system and cytotoxic activity is well known. Its repeated motives and modularity make Rtx a multifunctional factor in pathogenicity.
The comparative analysis of rtx gene among 6 strains of L. pneumophila showed modularity in their structures. Among compared genomes, the N-terminal region of the protein presents highly dissimilar repeats with functionally similar domains. On the contrary, the C-terminal region is maintained with a fashionable modular configuration, which gives support to its proposed role in adhesion and pore formation. Despite the variability of rtx among the considered strains, the flanking genes are maintained in synteny and similarity.
In contrast to the extracellular bacteria Vibrio cholerae, in which the rtx gene is highly conserved and flanking genes have lost synteny and similarity, the gene region coding for the Rtx toxin in the intracellular pathogen L. pneumophila shows a rapid evolution. Changes in the rtx could play a role in pathogenicity. The interplay of the Rtx toxin with host membranes might lead to the evolution of new variants that are able to escape host cell defences.
Legionella pneumophila is a gram negative, gamma-proteobacteria organism whose natural hosts are amoebae and protozoa. This bacterium can infect humans by inhalation of aerosols [1, 2] entering alveolar macrophages causing the well-known, and often lethal, Legionnaires' disease (LD) or Legionellosis. Despite the great number of isolates of L. pneumophila, the ones belonging to serogroup 1 are responsible of about 80 to 90% of cases of Legionellosis . The first critical event during infection by L. pneumophila involves the macrophages by the action of the type IV secretion system, which prevents the fusion of the phagosome with the lysosome and its acidification [4, 5]. It has been demonstrated that these events start very early after the infection . Several mechanisms play an important role in the formation of infection vacuoles. Legionella enters the macrophages by vacuoles that are morphologically similar to macropinosomes by an unusual mechanisms called "coiling phagocytosis" [6, 7]. The vacuole is immediately surrounded by vesicles and mitochondria and moves toward the endoplasmic reticulum escaping or delaying fusion with the lysosome . At this stage, the vacuoles offer a perfect niche for bacteria to multiply safely away from the lysosome. Legionella is also able to mediate the delayed entrance of the cell in apoptosis by modulating the activity of caspase-3 and other effectors [9, 10]. In all these stages the Dot/Icm (Defective in organelle trafficking/Intracellular multiplication) system, involved in the formation of type IV secretor machinery, is the main player acting on the transfer of a series of effectors in the host cell [11–14].
One of the first events in the pathogenic cascade is pore formation, which seems to be caused by a toxin belonging to the Rtx family ("repeats-in toxin") . It was demonstrated that the rtxA gene in L. pneumophila is strictly related to pathogenicity, and its main role involves adherence to the host membranes, thus facilitating all the molecular trafficking of the bacteria during infection processes . In other bacteria like Bordetella pertussis, Escherichia coli or Actinobacillus actinomycetemcomitans, proteins belonging to the Rtx family are also described as effectors of immune cell lyses and its action is often mediated by specific host membrane receptors [17, 18]. RtxA in L. pneumophila is a large protein (around 7.000 amino acids) with several repeated structures belonging to, at least, three protein family domains. In L. pneumophila strains Lens [GenBank:CR628337] and Paris [GenBank:CR628336] a correlation could exists between the number of repeats and greater invasion and virulence properties .
In this work, we analyse the mosaic structure of an RtxA toxin from the highly virulent L. pneumophila strain 2300/99. It was isolated in Alcoy (Spain) and was retrieved in several outbreaks in 1999 and 2000, during which more than 200 patients were infected, a dozen of whom died. In all cases, transmission was due to aerosol inhalation from out-door installations . The comparative analysis of the rtx locus of closely related L. pneumophila serogroup 1 strains showed the existence of a long tandem repeated domain of variable copy number and sequence. Furthermore, we have studied similarity and gene-order conservation of genes flanking the rtx region, finding remarkably high levels of rearrangements and diversity of the rtxA gene as compared to those from flanking regions. This pattern is completely opposite to that found in several strains of Vibrio cholerae, a phylogenetically close extracellular pathogen.
rtxA structure summary. The columns report respectively: strain with accession number in parenthesis; GenBank locus tag used to identify the gene; length in aminoacids; number and types of repeat and number and kind of domains identified by PFAM search.
Domains at N-terminal region
26 (type b1), 9 (type b2)
30 (type a)
16 (type c)
6 (type d)
After searching in the PFAM database, several kinds of adhesion related domains were identified as part of repeats type a (Paris), b1 (Lens), c1 and c2 (Corby and Alcoy). No domains related to adhesion were identified in rtx gene of the Philadelphia strain, while only the c1 and c2 repeats were phylogenetically related (see Additional file 1). Therefore, the way we approached the modular structure of the repeats was by looking at the function of the domains involved in adhesion, with only one exception: the AA100 sequence which spans the region located after the repeats, so it was not possible to include it in our description.
In the four completed genomes (Paris, Lens, Philadelphia and Corby) as well as in the contigs sequenced belonging to the Alcoy and AA100 strains, a von Willebrand factor type A domain (VWA) was identified, subsequent to the regions with tandem repeats. In addition, several blocks of tandem repeats identified as HemolysinCabind domains were also found. These blocks were formed by a number of different repeats: 3+2 repeats in Paris, 1+3+2 in Lens, 3+3+2 in Corby, 3+3+2 in Alcoy, 3+3+2 in Philadelphia and 1+3+2 in AA100. The latter domain was previously described and considered responsible for the virulent activity of the RtxA protein [16, 22–24]. A summary of all these structural features of the rtx present in different Legionella strains is shown in Table 1.
Comparative genomics and phylogenetic analyses
Strains Corby and Alcoy, on the other hand, are quite similar in their rtx gene at both protein domains and nucleotide level (Figure 1). Similarly to the other Legionella genomes, flanking genes are highly conserved both in order and orientation. N-terminal repeats are phylogenetically closely related (see Additional file 1, section C1+C2).
The alignment of repeated domains of the genomes considered here is practically impossible at nucleotide level, and even at amino acid level it is very complicated. Additional file 4 shows the amino acid alignment for each kind of repeat using CLUSTLALW and corrected by eye. The domains identified by PFAM are highlighted. It is not possible to identify any phylogenetic relationships either among the repeats or looking at the sole adhesion domains.
Here we report a comparative analysis of the Rtx toxin, as well as a fine analysis of repeats, identified in this protein in strains of L. pneumophila serogroup 1 from four completed genomes (strains Lens, Paris, Philadelphia and Corby), one shotgun ongoing sequencing project (strain Alcoy 2300/99) and one contig coming from a cosmid of the Legionella strain AA100. All these strains are known to be virulent [19, 21, 26, 27].
The rtx genes analyzed in the six genomes studied, present modularity. The toxin appears to be clearly divided in two regions, the N-terminal, involved in adhesion, and the C-terminal region, involved in adhesion and pore formation in the host membranes. The repeats in the N-terminal region analysed by PFAM database searches, have different kinds of adhesion domains. In all the N-terminal type a repeats of the Paris strain, similarity was found with the "Thrombospondin type 3 repeat" (TSP_3) of the human endothelial cell. In humans, this domain has been shown to bind fibrinogen, fibronectin, laminin and type V collagen [28, 29]. In the Lens strain, repetitions of type b1 show domains similar to the "Chlamydia polymorphic membrane protein" (Chlam_PMP). This obligate intracellular human pathogen causes infection of the upper and lower respiratory tract but the role of this membrane protein is still unknown . No similarity with specific domains was found for b2 type repeat of the Lens strain. The last kind of repeat-associated domains, identified in type c1 and c2 of Alcoy and Corby strains, were similar to the "Haemagglutinin" (HIM) domain that was found in invasins and hemagglutinins, and is associated with the Hep_Hag repeats .
Two other types of domains were commonly identified in all the rtx genes: the VWA domain involved in adhesion processes [32, 33], and another tandem repeated domain related to cytotoxic activity, the haemolysin calcium-binding (HemolysinCabind) site. The latter domain commonly brings a nonapeptide that is related to the adhesion with other host cell surfaces or vacuole membranes and pores formation [23, 34].
Except for Alcoy and Corby, which are very similar, the N-terminal repeats and their modular structure are highly variable among strains. Despite these differences, the flanking genes at 5' and 3' are highly conserved in sequence and order, which suggests that rtx undergoes a particular intra-genic evolution. Various examples of large non-interspaced repeats within CDS (Coding DNA Sequences) regions were described in bacteria such as E. coli and Bacillus subtilis, where recombination events were used to explain the distribution of large repeats among genomes [26, 35]. The case of rtx gene in L. pneumophila is particularly interesting due to the high number of observed intra-genic repeats. The origin of these repeats is yet unknown, as no similar sequences (by BLAST searches) have been identified in published data. As described in gene conversion models, DNA can enter to become part of a given gene  and afterwards, concerted evolution could be responsible for promoting the different adhesion domains. In fact, although the repeats type a, b, c, and d do not show any evolutionary relationships, the presence of different adhesion domains points towards a functional evolutionary advantage.
In V. cholerae the configuration of the rtx gene and its flanking regions is extremely different from that observed in L. pneumophila. Comparative analyses of rtx gene among five strains show a high level of conservation and the genes located at 5' of rtx, are strictly conserved, both in sequence and order, whereas the synteny, and sometimes nucleotide similarity, of those located at 3' do not follow the same pattern (see Additional file 3).
In Legionella spp. it was previously thought that only strains containing an active rtxA gene were able to produce infection in humans, and that mutants with a frame-shift inactivating rtxA protein were reduced for entry into host cells and pore formation in host membranes . rtxA seems to play a relevant role in the pathogenic activity of Legionella, although it also depends on the particular type of host. The variety of repeats and its homogeneous nature at rtxA N-terminal region of virulent strains of L. pneumophila seems to be acquired by the two mechanisms involved in concerted evolution: intra-genic gene conversion and/or unequal crossing over. As previously described for similar models in other organisms [36, 37], these mechanisms can be an important source of creating antigenic variability in Legionella, affording the ability to increase the host range and escape from the host's immune defences.
L. pneumophila serogroup 1 strain (numbered as 2300/99) was isolated from sputum of a patient with Legionnaires' disease and associated to LD outbreak detected in Alcoy (Alicante, Spain) in 1999, from hereinafter strain Alcoy.
DNA treatments library construction and sequencing
DNA from the L. pneumophila Alcoy was extracted as described in Ausubel et al. at "Centro Nacional de Microbiología del Instituto de Salud Carlos III (Majadahonda, Madrid, Spain)". Briefly, DNA was sheared by sonication for the construction of two libraries for inserts ranging from 1 to 2 kb and from 2 to 10 kb. Fragments were separated by "Pulsed Field Gel Electrophoresis" applying following conditions: Voltage, 2 V/cm; initial switch time, 0.1 s; final switch time, 1 s; temperature, 14°C; run time, 22 h; angle, 120°. After cutting bands out for the corresponding sizes, DNA was recovered from agarose by electroelution and subsequently purified by phenol/chloroform purification. Ends of fragments were flushed and tailed for TA cloning with "Single dA™ Tailing Kit" (Novagen, #69282-3). After flushing and tailing, DNA was precipitated with 2 Vol of ethanol 96% and 0.1 Vol of sodium acetate 3 M, and eluted in dH2O. Next, fragments were cloned with "TOPO XL PCR Cloning Kit" (Invitrogen, #K4700-10). Plasmid DNA purification was done with a "Montage Plasmid Miniprep 96 kit" (Millipore, #LSKP09624) in a "MULTIPROBE II-Robot Liquid Handling System". Sequencing reactions were carried out by the "ABI PRISM Big Dye Terminator v3.0 Ready Reaction Cycle Sequencing Kit" (Applied Biosystems, #4336919) and solved by the "3730XL Genetic Analyzer" (Applied Biosystems).
Sequences used in this paper are: L. pneumophila strains Paris [GenBank:CR628336], Lens [GenBank:CR628337], Corby [GenBank:CP000675], Alcoy [GenBank:EU054322], Philadelphia [GenBank:AE017354], AA100 [GenBank:AF057703]. Vibrio cholerae strains are: O395 [GenBank:NC_009457], O1 biovar eltor strain N16961 [GenBank:NC_002505], NCTC 8457 [GenBank:AAWD01000018], MZO-3 [GenBank:AAUU01000029] and MO10 [GenBank:DS179636].
The shotgun sequences assembly of the ongoing sequencing project of L. pneumophila Alcoy was carried out by the "Staden Package release v1.6.0" . ORF prediction of the yet unfinished strain Alcoy was performed by the use of "Glimmer 3" . Sequence manipulation and annotation was done with "Artemis" software . Comparative analyses were carried out by "BLAST" (Basic Local Alignment Search Tool) suite  and the visualization of the results obtained by means of the "ACT" (Artemis Comparative Tool) . Repetitions were identified by the software "Tandem Repeats Finder" . Finally, domain identification was carried out using the PFAM database .
For this analysis we selected three house keeping genes (acn, aconitate hydratase; groEL, 60 kDa chaperonin; recA, DNA recombination protein) and three fast evolving genes (flaA, flagellin; proA, zinc metalloprotease; mompS, major outer membrane protein) . The nucleotide sequences of each gene were concatenated and a tree based on Maximum Likelihood using Transition-Transversion plus Gamma substitution (as suggested by ModelTest analysis ) was obtained by the MEGA3.1  software. Reliability of the monophyletic groups was tested by a bootstrap test with 500 replicates.
Alcoy strain accession numbers
The genome project of Legionella pneumophila 2300/99 Alcoy is maintained at NCBI with project ID 18743. The new Alcoy sequences used in this work are available in GenBank: rtx gene contig, accession number EU054322. The genes used for MLST analysis have the following accession numbers: EU221241 (acn), EU221242 (flaA), EU221243 (groEL), EU221244 (mompS), EU221245 (proA), EU221246 (recA).
This work has been funded by grant BMC2006-06003 from MEC to AL and by contract with Conselleria de Sanidad of Valencian Government to AM. Nuria Jiménez is recipient of a fellowship from Carlos III and Giuseppe D'Auria has a research contract from CIBERESP. Sequencing was carried out using facilities of the SCSIE from University of Valencia.
- Lu H, Clarke M: Dynamic properties of Legionella-containing phagosomes in Dictyostelium amoebae. Cell Microbiol. 2005, 7: 995-1007. 10.1111/j.1462-5822.2005.00528.x.PubMedView ArticleGoogle Scholar
- Sabria M, Alvarez J, Dominguez A, Pedrol A, Sauca G, Salleras L: A community outbreak of Legionnaires' disease: evidence of a cooling tower as the source. Clin Microbiol Infect. 2006, 12: 642-647. 10.1111/j.1469-0691.2006.01447.x.PubMedView ArticleGoogle Scholar
- Yu VL, Plouffe JF, Pastoris MC, Stout JE, Schousboe M, Widmer A: Distribution of Legionella species and serogroups isolated by culture in patients with sporadic community-acquired legionellosis: an international collaborative survey. J Infect Dis. 2002, 186: 127-128. 10.1086/341087.PubMedView ArticleGoogle Scholar
- Horwitz MA: The Legionnaires' disease bacterium (Legionella pneumophila) inhibits phagosome-lysosome fusion in human monocytes. J Exp Med. 1983, 158: 2108-2126. 10.1084/jem.158.6.2108.PubMedView ArticleGoogle Scholar
- Amer AO, Swanson MS: Autophagy is an immediate macrophage response to Legionella pneumophila. Cell Microbiol. 2005, 7: 765-778. 10.1111/j.1462-5822.2005.00509.x.PubMedPubMed CentralView ArticleGoogle Scholar
- Roy CR, Berger KH, Isberg RR: Legionella pneumophila DotA protein is required for early phagosome trafficking decisions that occur within minutes of bacterial uptake. Mol Microbiol. 1998, 28: 663-674. 10.1046/j.1365-2958.1998.00841.x.PubMedView ArticleGoogle Scholar
- Watarai M, Derre I, Kirby J, Growney JD, Dietrich WF, Isberg RR: Legionella pneumophila is internalized by a macropinocytotic uptake pathway controlled by the Dot/Icm system and the mouse Lgn1 locus. J Exp Med. 2001, 194: 1081-1096. 10.1084/jem.194.8.1081.PubMedPubMed CentralView ArticleGoogle Scholar
- Horwitz MA: Formation of a novel phagosome by the Legionnaires' disease bacterium (Legionella pneumophila) in human monocytes. J Exp Med. 1983, 158: 1319-1331. 10.1084/jem.158.4.1319.PubMedView ArticleGoogle Scholar
- Gao LY, Abu KY: Activation of caspase 3 during Legionella pneumophila-induced apoptosis. Infect Immun. 1999, 67: 4886-4894.PubMedPubMed CentralGoogle Scholar
- bu-Zant A, Santic M, Molmeret M, Jones S, Helbig J, Abu KY: Incomplete activation of macrophage apoptosis during intracellular replication of Legionella pneumophila. Infect Immun. 2005, 73: 5339-5349. 10.1128/IAI.73.9.5339-5349.2005.View ArticleGoogle Scholar
- Bardill JP, Miller JL, Vogel JP: IcmS-dependent translocation of SdeA into macrophages by the Legionella pneumophila type IV secretion system. Mol Microbiol. 2005, 56: 90-103. 10.1111/j.1365-2958.2005.04539.x.PubMedView ArticleGoogle Scholar
- bu-Zant A, Jones S, Asare R, Suttles J, Price C, Graham J: Anti-apoptotic signalling by the Dot/Icm secretion system of L. pneumophila. Cell Microbiol. 2007, 9: 246-264. 10.1111/j.1462-5822.2006.00785.x.View ArticleGoogle Scholar
- Buscher BA, Conover GM, Miller JL, Vogel SA, Meyers SN, Isberg RR: The DotL protein, a member of the TraG-coupling protein family, is essential for Viability of Legionella pneumophila strain Lp02. J Bacteriol. 2005, 187: 2927-2938. 10.1128/JB.187.9.2927-2938.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Zink SD, Pedersen L, Cianciotto NP, bu-Kwaik Y: The Dot/Icm type IV secretion system of Legionella pneumophila is essential for the induction of apoptosis in human macrophages. Infect Immun. 2002, 70: 1657-1663. 10.1128/IAI.70.3.1657-1663.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Cirillo SL, Bermudez LE, El-Etr SH, Duhamel GE, Cirillo JD: Legionella pneumophila entry gene rtxA is involved in virulence. Infect Immun. 2001, 69: 508-517. 10.1128/IAI.69.1.508-517.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Samrakandi MM, Cirillo SL, Ridenour DA, Bermudez LE, Cirillo JD: Genetic and phenotypic differences between Legionella pneumophila strains. J Clin Microbiol. 2002, 40: 1352-1362. 10.1128/JCM.40.4.1352-1362.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Lally ET, Kieba IR, Sato A, Green CL, Rosenbloom J, Korostoff J: RTX toxins recognize a beta2 integrin on the surface of human target cells. J Biol Chem. 1997, 272: 30463-30469. 10.1074/jbc.272.48.30463.PubMedView ArticleGoogle Scholar
- Martin C, Requero MA, Masin J, Konopasek I, Goni FM, Sebo P: Membrane restructuring by Bordetella pertussis adenylate cyclase toxin, a member of the RTX toxin family. J Bacteriol. 2004, 186: 3760-3765. 10.1128/JB.186.12.3760-3765.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Cazalet C, Rusniok C, Bruggemann H, Zidane N, Magnier A, Ma L: Evidence in the Legionella pneumophila genome for exploitation of host cell functions and high genome plasticity. Nat Genet. 2004, 36: 1165-1173. 10.1038/ng1447.PubMedView ArticleGoogle Scholar
- Fernandez JA, Lopez P, Orozco D, Merino J: Clinical study of an outbreak of Legionnaire's disease in Alcoy, Southeastern Spain. Eur J Clin Microbiol Infect Dis. 2002, 21: 729-735. 10.1007/s10096-002-0819-9.PubMedView ArticleGoogle Scholar
- Cirillo SL, Lum J, Cirillo JD: Identification of novel loci involved in entry by Legionella pneumophila. Microbiology. 2000, 146 (Pt 6): 1345-1359.PubMedView ArticleGoogle Scholar
- Cirillo SL, Yan L, Littman M, Samrakandi MM, Cirillo JD: Role of the Legionella pneumophila rtxA gene in amoebae. Microbiology. 2002, 148: 1667-1677.PubMedView ArticleGoogle Scholar
- Economou A, Hamilton WD, Johnston AW, Downie JA: The Rhizobium nodulation gene nodO encodes a Ca2(+)-binding protein that is exported without N-terminal cleavage and is homologous to haemolysin and related proteins. EMBO J. 1990, 9: 349-354.PubMedPubMed CentralGoogle Scholar
- Sanchez-Magraner L, Viguera AR, Garcia-Pacios M, Garcillan MP, Arrondo JL, de la CF: The calcium-binding C-terminal domain of Escherichia coli alpha-hemolysin is a major determinant in the surface-active properties of the protein. J Biol Chem. 2007, 282: 11827-11835. 10.1074/jbc.M700547200.PubMedView ArticleGoogle Scholar
- Gaia V, Fry NK, Harrison TG, Peduzzi R: Sequence-based typing of Legionella pneumophila serogroup 1 offers the potential for true portability in legionellosis outbreak investigation. J Clin Microbiol. 2003, 41: 2932-2939. 10.1128/JCM.41.7.2932-2939.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Jepras RI, Fitzgeorge RB, Baskerville A: A comparison of virulence of two strains of Legionella pneumophila based on experimental aerosol infection of guinea-pigs. J Hyg (Lond). 1985, 95: 29-38.View ArticleGoogle Scholar
- Chien M, Morozova I, Shi S, Sheng H, Chen J, Gomez SM: The genomic sequence of the accidental pathogen Legionella pneumophila. Science. 2004, 305: 1966-1968. 10.1126/science.1099776.PubMedView ArticleGoogle Scholar
- Lawler J, Hynes RO: The structure of human thrombospondin, an adhesive glycoprotein with multiple calcium-binding sites and homologies with several different proteins. J Cell Biol. 1986, 103: 1635-1648. 10.1083/jcb.103.5.1635.PubMedView ArticleGoogle Scholar
- Kvansakul M, Adams JC, Hohenester E: Structure of a thrombospondin C-terminal fragment reveals a novel calcium core in the type 3 repeats. EMBO J. 2004, 23: 1223-1233. 10.1038/sj.emboj.7600166.PubMedPubMed CentralView ArticleGoogle Scholar
- Pedersen AS, Christiansen G, Birkelund S: Differential expression of Pmp10 in cell culture infected with Chlamydia pneumoniae CWL029. FEMS Microbiol Lett. 2001, 203: 153-159. 10.1111/j.1574-6968.2001.tb10834.x.PubMedView ArticleGoogle Scholar
- Tiyawisutsri R, Holden MT, Tumapa S, Rengpipat S, Clarke SR, Foster SJ: Burkholderia Hep_Hag autotransporter (BuHA) proteins elicit a strong antibody response during experimental glanders but not human melioidosis. BMC Microbiol. 2007, 7: 19-10.1186/1471-2180-7-19.PubMedPubMed CentralView ArticleGoogle Scholar
- Ruggeri ZM, Ware J: von Willebrand factor. FASEB J. 1993, 7: 308-316.PubMedGoogle Scholar
- Colombatti A, Bonaldo P, Doliana R: Type A modules: interacting domains found in several non-fibrillar collagens and in other extracellular matrix proteins. Matrix. 1993, 13: 297-306.PubMedView ArticleGoogle Scholar
- Valeva A, Walev I, Kemmer H, Weis S, Siegel I, Boukhallouk F: Binding of Escherichia coli hemolysin and activation of the target cells is not receptor-dependent. J Biol Chem. 2005, 280: 36657-36663. 10.1074/jbc.M507690200.PubMedView ArticleGoogle Scholar
- Rocha EP, Danchin A, Viari A: Analysis of long repeats in bacterial genomes reveals alternative evolutionary mechanisms in Bacillus subtilis and other competent prokaryotes. Mol Biol Evol. 1999, 16: 1219-1230.PubMedView ArticleGoogle Scholar
- Santoyo G, Romero D: Gene conversion and concerted evolution in bacterial genomes. FEMS Microbiol Rev. 2005, 29: 169-183. 10.1016/j.femsre.2004.10.004.PubMedView ArticleGoogle Scholar
- Deitsch KW, Moxon ER, Wellems TE: Shared themes of antigenic variation and virulence in bacterial, protozoal, and fungal infections. Microbiol Mol Biol Rev. 1997, 61: 281-293.PubMedPubMed CentralGoogle Scholar
- Ausubel FM, Brent R, Kingston RE, Moore DD, Seidman JG, Smith JA: Current protocols in Molecular Biology. 1989, Greene Publishing Associated and Wiley Interscience. New YorkGoogle Scholar
- Staden R, Beal KF, Bonfield JK: The Staden package, 1998. Methods Mol Biol. 2000, 132: 115-130.PubMedGoogle Scholar
- Delcher AL, Harmon D, Kasif S, White O, Salzberg SL: Improved microbial gene identification with GLIMMER. Nucleic Acids Res. 1999, 27: 4636-4641. 10.1093/nar/27.23.4636.PubMedPubMed CentralView ArticleGoogle Scholar
- Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, Rajandream MA: Artemis: sequence visualization and annotation. Bioinformatics. 2000, 16: 944-945. 10.1093/bioinformatics/16.10.944.PubMedView ArticleGoogle 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 TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J: ACT: the Artemis Comparison Tool. Bioinformatics. 2005, 21: 3422-3423. 10.1093/bioinformatics/bti553.PubMedView ArticleGoogle Scholar
- Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27: 573-580. 10.1093/nar/27.2.573.PubMedPubMed CentralView ArticleGoogle Scholar
- Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T: Pfam: clans, web tools and services. Nucleic Acids Res. 2006, 34: D247-D251. 10.1093/nar/gkj149.PubMedPubMed CentralView ArticleGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.PubMedView ArticleGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5: 150-163. 10.1093/bib/5.2.150.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.