This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Berent, L. M.
Right arrow Articles by Messick, J. B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Berent, L. M.
Right arrow Articles by Messick, J. B.

 Previous Article  |  Next Article 

Infection and Immunity, June 2003, p. 3657-3662, Vol. 71, No. 6
0019-9567/03/$08.00+0     DOI: 10.1128/IAI.71.6.3657-3662.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.

Physical Map and Genome Sequencing Survey of Mycoplasma haemofelis (Haemobartonella felis)

Linda M. Berent{dagger} and Joanne B. Messick*

Department of Veterinary Pathobiology, College of Veterinary Medicine, University of Illinois, Urbana, Illinois 61802

Received 11 October 2002/ Returned for modification 7 January 2003/ Accepted 11 March 2003


arrow
ABSTRACT
 
Mycoplasma haemofelis is an uncultivable red-cell pathogen of cats. Isolated M. haemofelis DNA was used to create a bacterial artificial chromosome library and physical map. Random sequencing of this material revealed 75 genes that had not been previously reported for M. haemofelis or any other hemotrophic mycoplasma.


arrow
TEXT
 
Mycoplasma haemofelis is a member of a newly defined group of mycoplasmas that parasitize the red blood cells of animals and humans. M. haemofelis, which was formerly known as Haemobartonella felis strain OH (18), was chosen for use in these molecular studies because it is the highly pathogenic causative agent of the syndrome historically reported as feline infectious anemia (8, 11). This syndrome is often associated with feline retroviral infections and may serve as a model for hemotrophic mycoplasmal infection in immunosuppressed human patients (13; M. I. Duarte, M. S. Oliveira, M. A. Shikanai-Yasuda, O. N. Mariano, C. F. Takakura, C. Pagliari, and C. E. Corbett, Letter, J. Infect. Dis. 165:976-977, 1992). This paper describes the genome size of M. haemofelis, the creation of the first genome encyclopedia and physical map of an uncultured mycoplasma, and the collection of M. haemofelis genome survey sequences (GSSs) for the purpose of gene discovery.

To acquire large quantities of M. haemofelis DNA for the following experiments, a feline leukemia virus- and feline immunodeficiency virus-seronegative cat was splenectomized and infected with M. haemofelis. When more than 60% of the red cells contained organisms, blood was aseptically drawn from the jugular vein into a syringe containing 1 ml of anticoagulant-citrate-dextrose solution per 5 ml of blood. The organisms were released from the red cells and embedded in agarose plugs, and the DNA was extracted by using previously described protocols (1, 19). Related organisms Mycoplasma genitalium (ATCC 49895) and Mycoplasma haemosuis (15), as well as a commercial DNA marker, were used as size controls in subsequent experiments.

Gamma radiation was used to linearize the bacterial DNA prior to pulsed-field gel electrophoresis (PFGE) to determine the full genome length of M. haemofelis (20, 31). A 1.0% 0.5x Tris-borate-EDTA gel was run for 24 h at 14°C and 6 V/cm, with 60- to 120-s switch times and a field angle of 120°. This PFGE was repeated four times with M. haemofelis plugs from separate blood collections. To confirm the identify of the bands seen, the Southern blotted DNA was probed by using a 393-bp digoxigenin-labeled fragment of the M. haemofelis 16S rRNA gene which spanned hypervariable regions 1 to 3 (1, 17). Membranes were prehybridized in 10 to 20 ml of PerfectHyb Plus (Sigma-Aldrich Corp., St. Louis, Mo.) at 68°C for at least 5 min. Labeled probe was boiled for 10 min and added to the prehybridization solution, and the probe was allowed to hybridize overnight at 68°C. Stringency washes and signal detection with CDP-star were performed according to the manufacturer's instructions (Boehringer Mannheim Biochemicals [Roche Molecular Biochemicals], Indianapolis, Ind.).

The average size of the linearized M. haemofelis chromosome was 1,199 kb (standard deviation, 13.5 kb). This falls well within the published range (580 to 1,400 kb) for the Mycoplasma pneumoniae subgroup of mycoplasmas (14, 20). When membranes were probed with a labeled fragment of the 16S rRNA gene of M. haemofelis, a strong signal was seen with the M. haemofelis DNA and weaker signals were observed for M. haemosuis and M. genitalium. This weak heterologous binding of our probe to the other mycoplasmas was expected due to the sequence similarities between the hypervariable regions of the 16S rRNA genes.

The creation of the bacterial artificial chromosome (BAC) library was done by using the Escherichia coli host strain DH10ß and the vector pBeloBAC11 (28). Full-length M. haemofelis genomic DNA was partially cut with HindIII (24), size-selected by using PFGE, and isolated from the gel by using standard techniques (26). A molar ratio of vector to insert of 5:1 was used in the ligation experiments. An aliquot of the ligation mix was added to freshly prepared DH10ß electrocompetent cells, and the cells were transformed in a 0.2-cm cuvette at 200 {Omega}, 25 µF, and 12.5 V/cm.

A probe made from full-length M. haemofelis genomic DNA was used to screen colony hybridization blots for clones containing M. haemofelis DNA and eliminate clones possibly containing E. coli or cat genomic DNA. One hundred and fifty colonies that exhibited strong hybridization signals were then used for mapping experiments. The average size of these clones was 42.6 kb, and there was a 99.7% probability that the entire genome was represented (6).

The BAC clones containing M. haemofelis DNA were cut with restriction endonucleases NotI, NruI (New England BioLabs, Beverly, Mass.), and SalI (Promega, Madison, Wis.) in either single or double digestions according to the manufacturers' instructions. The enzymes were chosen such that the pBeloBAC11 backbone would be excised and would serve as an internal control for complete digestion. These digestions were run on a 1.0% gel in 0.5x Tris-borate-EDTA buffer at 14°C for 15 h with switch times of 5 to 15 s at 6 V/cm and at an angle of 120o. The DNA was then Southern blotted onto a positively charged nylon membrane by using standard techniques (26).

To define a starting point for the restriction map, the M. haemofelis 16S rRNA gene probe was used to screen dot blots of BAC clones and identify a contig of six clones containing this gene. The probe should hybridize just before the NruI site at nucleotide 506 within the 16S rRNA gene (1, 25). Based on the restriction patterns of the clones and the placement of the 16S rRNA gene, it was determined that these six clones represented a single contig and that only one copy of this gene was present in M. haemofelis. This NruI site was designated as point zero for the purposes of this mapping project. Probes were then made from restriction fragments flanking the 16S rRNA gene and used to identify adjacent BAC clones. The order of the various restriction fragments from each BAC clone was determined by examining the restriction enzyme fingerprints and by hybridizing single and double restriction digests with probes made from the fragments in question. After repeated cycles of hybridizations of the dot blots, pulsed-field gel blots, and new probe selections, a complete genomic circle was closed. A total of 24 BAC clones were identified which represented the entire genome of M. haemofelis with minimal overlap (the genomic encyclopedia). The restriction map of M. haemofelis contains 28 NruI sites, 21 SalI sites, and 9 NotI sites.

A selection of BAC clones containing known regions on the M. haemofelis physical map were subcloned into the vector pGEM-3Zf(+) (Promega) and partially sequenced. A total of 624 GSSs were successfully generated. An average length of 450 bases per GSS was obtained, and this yielded 283,351 bp of genomic sequence. After elimination of duplicate sequences and combination of contigs, 430 unique GSSs were deposited in GenBank. The average guanine-plus-cytosine (G+C) percentage was 38.5% and showed little variation among the subclones. This G+C percentage is less than that found in M. pneumoniae (41%) (12) and greater than that seen in M. genitalium, Mycoplasma pulmonis, and Ureaplasma urealyticum (26 to 32%) (3, 9, 10). It also confirms that there is no contamination by E. coli DNA, which has a G+C percentage of 48 to 52% (2), or mammalian DNA, which has a G+C percentage of 52 to 54% (16).

Nine different contigs were identified that spanned overlapping areas on the genomic encyclopedia map. These contigs confirmed the overlap between clones 6.115 and 6.105, 6.67 and 6.25, 6.25 and 6.31, and 6.31 and 6.88 (Fig. 1). In addition, no contigs were detected between clones that did not overlap in the physical map. The 16S rRNA gene was found in two subclones of BAC 6.115 (as expected), and the other known M. haemofelis gene (rnpA) was found on the opposite side of the physical map in BAC clone 6.30. Finally, no GSSs were detected which represent feline or E. coli DNA present in the nonredundant database.



View larger version (33K):
[in this window]
[in a new window]
  		FIG. 1. M. haemofelis genomic encyclopedia and map. White bars indicate BAC clones that were not subcloned and gray bars indicate BAC clones used for subcloning and GSS generation. The genes discovered within each BAC clone are listed in random order by the abbreviations listed in Table 2. The scale is in kilobases.

The GSSs were compared to known sequences in the nonredundant database by using the BLASTn, BLASTx, and tBLASTx programs (National Center for Biotechnology Information; www.ncbi.nlm.nih.gov) with the mycoplasmal translation code. Among the GSSs obtained, 165 (26.4%) had a significant BLAST score ([E] of <=10-5) (27) with at least one entry in GenBank (Table 1). The majority of the most significant BLAST scores were for genes and proteins in the Mycoplasma and Ureaplasma genera. The remaining GSSs with significant BLAST scores were split evenly between members of the Bacillus-Clostridium group and other miscellaneous bacteria. Four hundred fifty-nine (73.6%) of the M. haemofelis GSSs failed to yield a significant database match ([E] of >=10-5). These results indicate that the M. haemofelis genome likely encodes a large number of unique proteins. Unlike most other mycoplasmas, which tend to colonize mucosal surfaces, the host-adapted survival of M. haemofelis is achieved through surface parasitism of the red blood cell, and it is therefore not surprising that this bacterium contains many unique genes.


View this table:
[in this window]
[in a new window]
  		TABLE 1. Database match categories of GSSs from M. haemofelis

The 165 GSSs with significant BLAST scores were further scrutinized to eliminate GSSs of the same gene. After deleting redundant GSSs, 77 out of 165 remained as nonredundant GSSs, representing 75 distinct genes in M. haemofelis that had not been previously reported. Two GSSs represent the 16S rRNA and RNase P RNA genes of M. haemofelis, both of which had been previously described (GenBank accession numbers U95297 and AF407210). Sequences were further grouped according to the cluster of orthologous groups (COG) and the presumptive gene function (29). The largest number of genes found were related to translation, ribosomal structure, and biogenesis (25%). Other groups with significant numbers of genes included those related to DNA replication (14%), carbohydrate metabolism (9%), transcription (6%), and nucleotide transport (6%) (Table 2).


View this table:
[in this window]
[in a new window]
  		TABLE 2. Putative genes identified in the GSSs of M. haemofelis, the BAC clone in which the gene is located, the affiliated COG, the functional category, and the GenBank accession number(s)a

Several putative virulence factors were detected in this analysis which merit further study. These include membrane lipoproteins, possibly related to antigenic variation systems (5), and other well-described virulence factors such as VacB and MgpA-like proteins (4, 7).

The superoxide dismutase (SOD) gene appears to be unique to M. haemofelis among the mycoplasmas. There are studies describing SOD activity in some mycoplasmas (21), but there are currently no available reports of a mycoplasma SOD gene sequence available in GenBank or from the fully sequenced mycoplasmas (3, 9, 10, 12). The best BLAST scores found when comparing the putative M. haemofelis SOD gene to the nonredundant database were with bacteria from the Bacillus-Clostridium group. This suggests that the SOD gene was present in a common ancestor and either that other mycoplasmas may have lost the gene or that a product of a different gene has assumed the role of SOD production.

Of particular interest are two genes found in M. haemofelis that are involved in purine biosynthesis and encode proteins similar to IMP (inosine-5'-monophosphate) dehydrogenase and GMP (guanosine 3',5'-monophosphate) synthase. There are no reports describing the enzymatic activity of these proteins in Mollicutes spp., and thus it was assumed that all Mollicutes spp. require a source of preformed guanine (23, 30). The mere presence of these genes in M. haemofelis does not definitively prove that the enzymes have the predicted activity (22). If M. haemofelis could be grown, the culture could be examined for the enzyme activity, but with M. haemofelis we are limited to using an expression vector to study this protein in the future.

This genomic encyclopedia and collection of survey sequences provides the first physical map and genomic survey of an uncultured hemotrophic mycoplasma. It is our hope that this information will ultimately lead to the development of more effective treatments and detection of M. haemofelis by characterizing its virulence and biosynthetic pathways.

Nucleotide sequence accession numbers. The 430 unique GSSs were deposited in GenBank under the accession numbers BH792926 through BH793355.


arrow
FOOTNOTES
 
* Corresponding author. Mailing address: Department of Veterinary Pathobiology, College of Veterinary Medicine, 2001 S. Lincoln Ave., University of Illinois, Urbana, IL 61802. Phone: (217) 333-5372. Fax: (217) 333-5342. E-mail: jmessick{at}cvm.uiuc.edu. Back

Editor: J. T. Barbieri

{dagger} Present address: University of Missouri, College of Veterinary Medicine, Columbia, MO 65211. Back


arrow
REFERENCES
 
    1
  1. Berent, L. M., J. B. Messick, and S. K. Cooper. 1998. Detection of Haemobartonella felis in cats with experimentally induced acute and chronic infections, using a polymerase chain reaction assay. Am. J. Vet. Res. 59:1215-1220.[Medline]
    2. 2
  2. Blattner, F. R., G. Plunkett III, C. A. Bloch, N. T. Perna, V. Burland, M. Riley, J. Collado-Vides, J. D. Glasner, C. K. Rode, G. F. Mayhew, J. Gregor, N. W. Davis, H. A. Kirkpatrick, M. A. Goeden, D. J. Rose, B. Mau, and Y. Shao. 1997. The complete genome sequence of Escherichia coli K-12. Science 277:1453-1474.[Abstract/Free Full Text]
    3. 3
  3. Chambaud, I., R. Heilig, S. Ferris, V. Barbe, D. Samson, F. Galisson, I. Moszer, K. Dybvig, H. Wroblewski, A. Viari, E. P. Rocha, and A. Blanchard. 2001. The complete genome sequence of the murine respiratory pathogen Mycoplasma pulmonis. Nucleic Acids Res. 29:2145-2153.[Abstract/Free Full Text]
    4. 4
  4. Cheng, Z. F., Y. Zuo, Z. Li, K. E. Rudd, and M. P. Deutscher. 1998. The vacB gene required for virulence in Shigella flexneri and Escherichia coli encodes the exoribonuclease RNase R. J. Biol. Chem. 273:14077-14080.[Abstract/Free Full Text]
    5. 5
  5. Citti, C., and K. S. Wise. 1995. Mycoplasma hyorhinis vlp gene transcription: critical role in phase variation and expression of surface lipoproteins. Mol. Microbiol. 18:649-660.[CrossRef][Medline]
    6. 6
  6. Clarke, L., and J. Carbon. 1976. A colony bank containing synthetic Col El hybrid plasmids representative of the entire E. coli genome. Cell 9:91-99.[CrossRef][Medline]
    7. 7
  7. Dallo, S. F., A. Chavoyo, C. J. Su, and J. B. Baseman. 1989. DNA and protein sequence homologies between the adhesins of Mycoplasma genitalium and Mycoplasma pneumoniae. Infect. Immun. 57:1059-1065.[Abstract/Free Full Text]
    8. 8
  8. Foley, J. E., S. Harrus, A. Poland, B. Chomel, and N. C. Pedersen. 1998. Molecular, clinical, and pathologic comparison of two distinct strains of Haemobartonella felis in domestic cats. Am. J. Vet. Res. 59:1581-1588.[Medline]
    9. 9
  9. Fraser, C. M., J. D. Gocayne, O. White, M. D. Adams, R. A. Clayton, R. D. Fleischmann, C. J. Bult, A. R. Kerlavage, G. Sutton, J. M. Kelley, J. L. Fritschman, J. F. Weidman, K. V. Small, M. Sandusky, D. Nguyen, T. R. Utterback, D. M. Saudek, C. A. Phillips, J. M. Merrick, J.-F. Tomb, B. A. Dougherty, K. F. Bott, P.-C. Hu, T. S. Lucier, S. N. Peterson, H. O. Smith, C. A. I. Hutchinson, and J. C. Venter. 1995. The minimal gene complement of Mycoplasma genitalium. Science 270:397-403.[Abstract/Free Full Text]
    10. 10
  10. Glass, J. I., E. J. Lefkowitz, J. S. Glass, C. R. Heiner, E. Y. Chen, and G. H. Cassell. 2000. The complete sequence of the mucosal pathogen Ureaplasma urealyticum. Nature 407:757-762.[CrossRef][Medline]
    11. 11
  11. Harvey, J. W., and G. JM. 1977. Experimental feline haemobartonellosis. J. Am. Anim. Hosp. Assoc. 13:28-38.
    12. 12
  12. Himmelreich, R., H. Hilbert, H. Plagens, E. Pirkl, B.-C. Li, and R. Herrmann. 1996. Complete sequence analysis of the genome of the bacterium Mycoplasma pneumoniae. Nucleic Acids Res. 24:4420-4449.[Abstract/Free Full Text]
    13. 13
  13. Kallick, C. A., S. Levin, K. T. Reddi, and W. L. Landau. 1972. Systemic lupus erythematosus associated with Haemobartonella-like organisms. Nat. New Biol. 236:145-146.[Medline]
    14. 14
  14. Maniloff, J. 1992. Mycoplasmas: molecular biology and pathogenesis. American Society for Microbiology, Washington, D.C.
    15. 15
  15. Messick, J. B., G. Smith, L. Berent, and S. Cooper. 2000. Genome size of Eperythrozoon suis and hybridization with 16S rRNA gene. Can. J. Microbiol. 46:1082-1086.[CrossRef][Medline]
    16. 16
  16. Mouchiroud, D., and G. Bernardi. 1993. Compositional properties of coding sequences and mammalian phylogeny. J. Mol. Evol. 37:109-116.[CrossRef][Medline]
    17. 17
  17. Neefs, J. M., Y. Van de Peer, L. Hendriks, and R. De Wachter. 1990. Compilation of small ribosomal subunit RNA sequences. Nucleic Acids Res. 18(Suppl.):2237-2317.
    18. 18
  18. Neimark, H., K. E. Johansson, Y. Rikihisa, and J. G. Tully. 2001. Proposal to transfer some members of the genera Haemobartonella and Eperythrozoon to the genus Mycoplasma with descriptions of ‘Candidatus Mycoplasma haemofelis’, ‘Candidatus Mycoplasma haemomuris’, ‘Candidatus Mycoplasma haemosuis’ and ‘Candidatus Mycoplasma wenyonii’. Int. J. Syst. Evol. Microbiol. 51:891-899.[Abstract]
    19. 19
  19. Neimark, H., D. Mitchelmore, and R. H. Leach. 1998. An approach to characterizing uncultivated prokaryotes: the Grey Lung agent and proposal of a Candidatus taxon for the organism, ‘Candidatus Mycoplasma ravipulmonis’. Int. J. Syst. Bacteriol. 48:389-394.[Abstract/Free Full Text]
    20. 20
  20. Neimark, H. C., and C. S. Lange. 1990. Pulse-field electrophoresis indicates full-length Mycoplasma chromosomes range widely in size. Nucleic Acids Res. 18:5443-5448.[Abstract/Free Full Text]
    21. 21
  21. O'Brien, S. J., J. M. Simonson, M. W. Grabowski, and M. F. Barile. 1981. Analysis of multiple isoenzyme expression among twenty-two species of Mycoplasma and Acholeplasma. J. Bacteriol. 146:222-232.[Abstract/Free Full Text]
    22. 22
  22. Pollack, J. D. 2001. Ureaplasma urealyticum: an opportunity for combinatorial genomics. Trends Microbiol. 9:169-175.[CrossRef][Medline]
    23. 23
  23. Pollack, J. D., M. V. Williams, and R. N. McElhaney. 1997. The comparative metabolism of the mollicutes (Mycoplasmas): the utility for taxonomic classification and the relationship of putative gene annotation and phylogeny to enzymatic function in the smallest free-living cells. Crit. Rev. Microbiol. 23:269-354.[Medline]
    24. 24
  24. Riethman, H., B. W. Birren, and A. Gnirke. 1997. Preparation, manipulation, and mapping of HMW DNA, p. 83-248. In B. W. Birren, E. D. Green, S. Klapholz, R. M. Myers, and J. Roskams (ed.), Genome analysis—a laboratory manual, vol. 1. Cold Spring Harbor Laboratory Press, Plainview, N.Y.
    25. 25
  25. Rikihisa, Y., M. Kawahara, B. Wen, G. Kociba, P. Fuerst, F. Kawamori, C. Suto, S. Shibata, and M. Futohashi. 1997. Western immunoblot analysis of Haemobartonella muris and comparison of 16S rRNA gene sequences of H. muris, H. felis, and Eperythrozoon suis. J. Clin. Microbiol. 35:823-829.[Abstract]
    26. 26
  26. Sambrook, J., and D. W. Russell. 2001. Molecular cloning: a laboratory manual, 3rd ed. Cold Spring Harbor Press, Cold Spring Harbor, N.Y.
    27. 27
  27. Sanchez, D. O., R. O. Zandomeni, S. Cravero, R. E. Verdun, E. Pierrou, P. Faccio, G. Diaz, S. Lanzavecchia, F. Aguero, A. C. Frasch, S. G. Andersson, O. L. Rossetti, O. Grau, and R. A. Ugalde. 2001. Gene discovery through genomic sequencing of Brucella abortus. Infect. Immun. 69:865-868.[Abstract/Free Full Text]
    28. 28
  28. Shizuya, H., B. Birren, U. J. Kim, V. Mancino, T. Slepak, Y. Tachiiri, and M. Simon. 1992. Cloning and stable maintenance of 300-kilobase-pair fragments of human DNA in Escherichia coli using an F-factor-based vector. Proc. Natl. Acad. Sci. USA 89:8794-8797.[Abstract/Free Full Text]
    29. 29
  29. Tatusov, R. L., E. V. Koonin, and D. J. Lipman. 1997. A genomic perspective on protein families. Science 278:631-637.[Abstract/Free Full Text]
    30. 30
  30. Tryon, V. V., and J. D. Pollack. 1985. Distinctions in Mollicutes purine metabolism: pyrophosphate-dependent nucleoside kinase and dependence on guanylate salvage. Int. J. Syst. Bacteriol. 35:497-501.[Abstract/Free Full Text]
    31. 31
  31. van der Bliek, A. M., C. R. Lincke, and P. Borst. 1988. Circular DNA of 3T6R50 double minute chromosomes. Nucleic Acids Res. 16:4841-4851.[Abstract/Free Full Text]

Infection and Immunity, June 2003, p. 3657-3662, Vol. 71, No. 6
0019-9567/03/$08.00+0     DOI: 10.1128/IAI.71.6.3657-3662.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Berent, L. M.
Right arrow Articles by Messick, J. B.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Berent, L. M.
Right arrow Articles by Messick, J. B.

Copyright © 2009 by the American Society for Microbiology.   For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»