Previous Article | Next Article ![]()
Journal of Bacteriology, September 2005, p. 6386-6395, Vol. 187, No. 18
0021-9193/05/$08.00+0 doi:10.1128/JB.187.18.6386-6395.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
McGill University Health Centre, Montreal, Quebec H3G 1A4, Canada,1 Veterinary Laboratories Agency (Weybridge), Woodham Lane, New Haw, Addlestone, Surrey KT15 3NB, United Kingdom,2 Grupo de Genetica de Micobacterias, Departamento de Microbiologia Medicina Preventiva y Salud Publica, Universidad de Zaragoza, Spain,3 MRC Centre for Molecular and Cellular Biology, Department of Medical Biochemistry, Faculty of Health Sciences, University of Stellenbosch, P.O. Box 19063, Tygerberg 7505, South Africa,4 Mycobacteria Reference Unit, Diagnostic Laboratory of Infectious Diseases and Perinatal Screening, National Institute of Public Health and the Environment, Bilthoven, The Netherlands,5 Australian Reference Laboratory for Bovine Tuberculosis, Department of Agriculture, South Perth 6151, Australia6
Received 25 April 2005/ Accepted 27 June 2005
|
|
|---|
|
|
|---|
With the availability of genomic data from sequenced strains of the MTC (7, 11, 16, 18) and tools of comparative genomics (5, 45), it is now recognized that large sequence polymorphisms (LSPs) deleted from M. tuberculosis serve as accurate markers for diagnostic testing (23, 42) and, with the exception of mycobacteriophage DNA, as evolutionary markers (5, 8, 22, 24, 35, 36, 39, 50). As such, the presence or absence of these LSPs has been used to reveal associations between strains and host populations (22, 34, 37). By extension, genomic information about any LSPs specific to M. bovis and "M. bovis-like" organisms could guide in assessing the prevalence of the organisms through space and host type. A practical utility of this approach would be to provide a secure genomic definition for M. bovis and to help determine whether M. bovis is a single discrete genomic entity or instead a genomic continuum comprising strains that vary according to a number of properties, including (but not limited to) virulence and host preference.
We therefore selected M. bovis and M. bovis-like organisms isolated from various mammalian host types of diverse geographic origins, which present an array of genetic profiles and laboratory properties. The DNA of these organisms was then interrogated via Affymetrix GeneChip, followed by confirmatory PCR and sequencing. The results obtained reveal the previously unknown evolution of the MTC, provide markers for diagnostic testing and molecular epidemiologic assessment, and point to both geographic and host-specific forms of genomic variability.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Characteristics of bacterial isolates used in this study
|
GeneChip analysis. The M. tuberculosis Affymetrix GeneChip is an oligonucleotide array representing sequences from the M. tuberculosis genome (45). When applied towards MTC genomic interrogation, the GeneChip is capable of identifying LSPs of hundreds of nucleotides or larger with respect to M. tuberculosis. Eight micrograms of MTC DNA was extracted, fragmented, biotin labeled, and hybridized to the GeneChip as previously described (45). Where necessary, extracted DNA was amplified using a GenomiPhi DNA amplification kit (Amersham Biosciences) to ensure that an appropriate quantity of DNA could be advanced to GeneChip and/or downstream PCR analysis. Affymetrix Microarray Suite (MAS) was used to analyze the comparative genomic data and ultimately suggest candidate deleted regions (34, 36, 37, 39).
PCR amplification and sequencing across deletions. Candidate deleted regions inferred by GeneChip analysis were pursued as described in previous publications (34, 36, 37, 39). In brief, primers were designed to amplify regions in M. tuberculosis H37Rv flanking the putative deletion. The test strain(s) presenting the candidate deleted region, together with referent strains of M. tuberculosis and/or M. bovis BCG (included as DNA controls), was concurrently subjected to PCR amplification. PCR amplicons were run on 2% agarose gel, where products different from the expected base pair size of M. tuberculosis H37Rv could be visualized. All products obtained by amplification across a deleted region were subsequently sequenced by dideoxy terminal sequencing at the McGill University and Genome Quebec Innovation Center. Sequence results were compared by BLAST analysis to M. tuberculosis H37Rv with TubercuList (http://genolist.pasteur.fr/TubercuList/) and to M. bovis AF2122/97 with BoviList (http://genolist.pasteur.fr/BoviList/) to confirm whether the amplified MTC DNA harbored a specific deletion event. Annotation of open reading frames (ORFs) affected by deletion events was determined by using TubercuList. Nomenclature of deleted regions disclosed in this study, hereafter referred to as "new deletions," follows previously established guidelines (7, 32, 34) and primarily denotes the MTC lineage from which strains harboring that region of difference (RD) can be specifically isolated. Nomenclature for new deletions further reflects an already published RD should it overlap the genetic region affected by a newfound RD, deleted genes for which a function has already been ascribed, the host's geographic source, and/or is numbered according to its order of discovery.
Analysis of deletions. Primers designed to amplify new deletions are provided as supplemental material (Table S4). A genomic deletion, herein defined as a region of the ancestral genome that has been deleted only from a subgroup of clonally related organisms, is only considered a unidirectional deletion event if confirmed by PCR and sequence analysis (22, 35). For isolates applied to GeneChip, regions were called present if denoted as such via MAS. Conversely, new deletions among interrogated isolates were sequence confirmed to determine whether the exact same genomic event was being identified. The distributions of deletions already described as absent from M. bovis (5, 18, 20) were based on MAS analysis and were not pursued by further analysis. Deletion events were afterward assigned to the previously characterized bacterial isolates. Regions of difference specific to M. tuberculosis (8), namely, TbD1, RvD1, RvD2, RvD3, RvD4, and RvD5, are not available for GeneChip analysis. To confirm the presence of these regions across isolates applied to GeneChip, PCR primers designed to amplify genes internal to these deletion events were applied (Table S4).
|
|
|---|
A total of 29 new deletions were confirmed, affecting 152,480 bp and 164 ORFs (Table 2). The distribution of these deletions across isolates is presented in Table 3, providing distinct genomic groupings observed to correlate with spoligotypes that resembled each other (Table 1), and is independently discussed below for isolates classified as M. pinnipedii, oryx bacillus, M. caprae, and M. bovis sensu stricto (Fig. 1).
|
View this table: [in a new window] |
TABLE 2. Confirmed new deletions among isolates, determined by GeneChip analysisa
|
|
View this table: [in a new window] |
TABLE 3. Distribution of large sequence polymorphisms among isolatesa
|
![]() View larger version (16K): [in a new window] |
FIG. 1. Revised deletion-based phylogeny of the MTC, based on deleted regions demonstrated through genomic analysis. The vertical axis presents the stepwise accumulation of unidirectional deletion events (RDs and N [new]-RDs) characterized among members of the MTC. Clustered along each horizontal axis are organisms for which one or more genomic deletions specific to this evolutionary branch have been revealed. Individual deletions revealed herein for M. pinnipedii, oryx bacillus, M. caprae, and multiple genomovars of M. bovis [(a) through (d)], and information about the ORFs they affect, are described in the supplemental material (Table S5). Beyond RD7 to RD10, M. microti, M. africanum (RD7 to RD10 deleted), M. pinnipedii, the dassie bacillus, and the oryx bacillus do not present a unidirectional deletion event shared with M. caprae and/or M. bovis; as a result, their respective order within this MTC phylogenetic model based on genomic deletions is arbitrary.
|
Because M. pinnipedii a975 revealed a unique genomic profile compared to those of other MTC members, deletions observed as specific to it, namely, RD2seal and RDpin, were tested for across a heterogeneous panel of isolates also characterized as M. pinnipedii (Table S1). With all isolates tested lacking both regions, these analyses revealed a common deletion profile for M. pinnipedii, irrespective of geography or pinniped host type. In all, the genomic content of M. pinnipedii isolates suggests a unique phylogenetic position within the MTC (13), genomically distinct from M. bovis (Fig. 1).
Deletions specific to oryx bacillus. In addition to deletion events shared with other MTC members (Table 3), we detected in both oryx bacillus isolates (kremer 24 and kremer 69) five new deletions not observed in other MTC members (Table 2). Of note, all five new deletions from these isolates (which have 18 and 19 copies of IS6110 [Table 1]) present an IS6110 sequence at their deletion junction, suggesting their mediation of these deletion events (15, 33). Three deletion events specific to M. tuberculosis, RvD2, RvD3, and RvD4 (8), each overlap independent deletion events unique to oryx bacillus, RDoryx_1, RDoryx_wag22, and RD5oryx, respectively.
Deletions specific to oryx bacillus could not be tested across a heterogeneous panel of isolates also characterized as oryx bacillus, because these isolates are rarely studied. Potentially, other MTCs isolated from an oryx host could reveal a common deletion profile for oryx bacillus irrespective of geography or oryx host type, though this remains to be determined. Until then, and in agreement with previous analysis genetically differentiating it from M. bovis (8, 27), including the observation that these two isolates present the only spoligotype pattern to have retained spacers 40 to 43, the oryx bacillus appears to comprise an exclusive phylogenetic position within the MTC (Fig. 1).
Deletions specific to M. caprae.
M. caprae M57, isolated from a caprine host in Spain, was missing seven deleted regions not previously documented, here named RDcap_Spain1 through RDcap_Spain7 (Table 2). Though the region flanking RvD1 appeared intact from GeneChip analysis, primers internal to RvD1 failed to amplify in M57. Seeing that primers for this region amplified as present in all other MTC isolates tested here and elsewhere for MTC other than M. tuberculosis (8), the absence of this internal sequence as it relates to a deletion event specific to isolates of M. caprae requires further investigation. Though M57 contains seven IS6110 elements, these deletions do not appear to be generally mediated by insertion elements but involve the rearrangement of other repeat sequences in the genome, namely, PE (Pro-Glu motif) and PPE (Pro-Pro-Glu motif) genes. Functionally noteworthy among these deletions is RDcap_Spain2, which affects the genes gca (Rv0112), gmhA (Rv0113), and gmhB (Rv0114) involved in the biosynthesis of nucleotide-activated glycero-manno-heptose precursors of lipopolysaccharides (51). Also notable is RDcap_Spain7, which deletes 11 genes (Rv3884c to Rv3894c), including the mycosin mycP2 (9), and the ESAT-6-like proteins encoded by esxC and esxD (19). Genes within this region overlap two independent deletions in M. bovis described below [RDpan and RDbovis(a)_
pan], suggesting it as a hotspot of genomic loss. M. caprae 60312, isolated from a human with TB who originally came from Cambodia, is not missing any of these seven deletions but instead lacks three genomic regions unique to it, here named RDcap_Asia1 to RDcap_Asia3 (Table 2). The specificities of these deletions from M. caprae are confirmed by the observation that they are not observed as deleted in any M. pinnipedii, oryx bacillus, or M. bovis isolates tested, but their distribution in other settings requires a greater panel of isolates. Importantly, both M. caprae M57 and M. caprae 60312 have retained the RD4 region as intact (Table 3), corroborating their phylogenetic distinction from M. bovis (1, 8, 35) (Fig. 1).
Because M. caprae M57 possessed a new set of deletions potentially specific to caprine isolates from Spain, these deletions were tested for across a diverse panel of isolates previously characterized as M. caprae from different farms across Spain to assess their geographic prevalence (Table S2). A common deletion profile for Spanish isolates of M. caprae was discovered, in that all isolates were missing RDcap_Spain2, RDcap_Spain3, RDcap_Spain5, RDcap_Spain6, and RDcap_Spain7, suggesting regional dominance of a clonal genotype.
Deletions specific to M. bovis. In deleting RD7, RD8, RD9, RD10, RD5, RD12, RD13, N-RD25, and RD4 (Table 3), the remaining 23 isolates demonstrated the anticipated deletion profile of M. bovis (8, 35). Additionally, 13 new deletions were identified among all M. bovis isolates interrogated via GeneChip (Table 2). Assessing the distribution of these new deletions revealed four distinct genomovars of M. bovis, here referred to as M. bovis (a) through (d), characterized by the presence/absence of N-RD17 (45), a 2,405-bp deletion flanking the PAN promoter region used to identify mycobacteria (44), hereafter referred to as RDpan, and RD17 (21) (Fig. 1).
M. bovis (a).
Unlike all other isolates having deleted RD4, only M. bovis a143 and a163 have retained N-RD17 and RD17 (Table 3). Of note, both of these isolates have deleted RDbovis(a)_
pan, an independent genetic deletion event with a different deletion location from what has been described for RDpan by Rauzier and colleagues among isolates of M. bovis (44). Furthermore, only a143 and a163 share RDbovis(a)_kdp, a deletion truncating three genes (Rv1028A to Rv1030) annotated as components of the high-affinity ATP-driven potassium transport (or KDP) system. Together, these data suggest that a143 and a163, both isolated from captive bovine hosts inTanzania, represent members of an ancestral M. bovis lineage (Fig. 1).
M. bovis (b). Though no isolates belonging to this genomovar were directly interrogated for analysis here, M. bovis (b) represents the initial lineage lacking N-RD17. The existence of this genomovar is supported by the observation that M. bovis BCG vaccines previously interrogated via GeneChip analysis (36) share this same deletion profile, pointing to their origin from M. bovis presenting this genotype (Fig. 1).
M. bovis (c). The deletion of RDpan was observed for a number of M. bovis isolates having one to five copies of IS6110 and infecting multiple animal species (Table 1). Though all retained RD17, isolates belonging to this intermediate genomovar of M. bovis variably presented other new deletions (Table 3).
Three isolates from different hosts in Canada with the same spoligotype pattern shared the genomic deletion of RDbovis(c)_virS affecting Rv3082c to Rv3085. This deletion includes a gene annotated as a virulence regulating transcriptional regulator (virS) that has also been independently deleted from the attenuated dassie bacillus (34). A different deletion, RDbovis(c)_Kruger, was restricted to a subset of isolates from Kruger National Park in South Africa (659, 1462, 1595F, and 734). As was observed with the Canadian deletion, the regional success of a particular clone having deleted RDbovis(c)_Kruger can be inferred from the observation that only isolates from Kruger National Park presented this deletion but that different hosts were implicated (Table 3). RDbovis(c)_Kruger was the largest deletion revealed in this study, affecting 23,877 bp and 20 ORFs, including two genes annotated as transcriptional regulators (Rv1358 and Rv1359) as well as the anti-anti sigma factor rsfA (Rv1365c) involved in regulating SigF (4). The fact that other strains from South Africa (9969-1 and 1136) have not deleted RDbovis(c)_Kruger but instead have deleted RDpan and RD17 indicates at least two genomovars of M. bovis in this region.
M. bovis (d). This most derivative group of M. bovis was observed for 11 isolates in which RD17 was deleted (Table 3 and Fig. 1). The previously described deletion of RD17 affects treY (Rv1563c), truncating one of three known pathways for the synthesis of trehalose (14, 38). GeneChip analysis revealed seven new deletions unique to isolates missing RD17 (Tables 2 and 3). Of these, one deletion, RDbovis(d)_sigK, disrupts sigK (Rv0445c), an extracytoplasmic sigma factor recently implicated in the control of MPB70/83 production by BCG strains (10). The disruption of sigK function suggests a possible explanation for the occasional lab result where isolates of M. bovis present as negative for MPB70/83 production (29). Another deletion, RD12HUP, was only observed as deleted from a South African M. bovis isolate (1136), affecting 4,675 bp and Rv3109 to Rv3115. RD12HUP removes a cluster of genes annotated as being involved in molybdenum cofactor biosynthesis, just upstream of the RD12 region (Rv3117 to Rv3121) already missing from all isolates of M. caprae and M. bovis.
Specific to isolates of farmed East Asian water buffalo from Australia (a100 and a181), the deletion of RDbovis(d)_buff1 and RDbovis(d)_buff2 potentially indicates a strain of M. bovis that preferentially infects these hosts. Of functional interest from other studies, RDbovis(d)_buff1 deletes 1,052 bp affecting fadD30 (Rv0404), a gene annotated as involved in lipid degradation and found to be essential for in vitro growth in two independent studies (28, 46).
GeneChip analysis of GB M. bovis strains from badgers was performed on isolates from three GB-dominant spoligotype groups (48). Badger isolate 59/0173/02 from Gloucestershire showed the RDbovis(d)_0173 deletion, a lesion that removes genes involved in lipid and intermediary metabolism (Rv1928c to Rv1936) and includes one gene predicted to encode a transcriptional regulatory protein of the AraC family (Rv1931c). The East Sussex badger isolate (61/1160/97) showed a deletion, RDbovis(d)_1160_1, that removes genes encoding a predicted regulator of the ArsR family (Rv0576) and others encoding proteins of unknown function (Rv0576 to Rv0585c). Badger isolate 61/1160/97 also presented a second deletion of almost 1.8 kb that encompassed Rv3403c to Rv3406 and removes a TetR family regulator and a putative dioxygenase [RDbovis(d)_1160_2]. To determine if deletions discovered by array analysis were absent from other M. bovis strains, we then checked a panel of 154 isolates of diverse spoligotypes from cattle, badger, and deer hosts (Table S3). Irrespective of the host type from which these strains were isolated, it was found that new deletions were lost only from the subset of strains sharing the same spoligotype pattern as those applied to array analysis, thereby offering no evidence for a host-specific deletion profile. These data, in agreement with previous genetic analysis of the population structure of M. bovis in GB by spoligotyping and variable number of tandem repeats (48), provide little evidence for the appearance of host-adapted clones in GB.
|
|
|---|
Genomic analyses have consistently described M. bovis as the most derivative MTC member (8, 35), though only in this report have LSPs been sought to study M. bovis itself. From this, a phylogenetic stratification of M. bovis genotypes (Fig. 1) offers a structure upon which to base future M. bovis study. Though it is unfortunate that isolates from the genomovar of M. bovis (b) were not analyzed here, the fact that BCG strains belong to this genomovar indicates that genetic variability between BCG strains and sequenced strains of either M. tuberculosis or M. bovis may result from polymorphisms specific to this genetic lineage and not necessarily mutations incurred during in vitro passage of BCG. The most derived genomovar of M. bovis characterized here was called M. bovis (d) and includes the sequenced strain of M. bovis AF2122/97 (18), which now de facto serves as the referent profile for M. bovis. Two practical consequences of these findings pertain to genetic and phenotypic studies of M. bovis. First, in genomic comparisons of M. tuberculosis, M. bovis, and BCG, it is noteworthy that M. bovis AF2122/97 harbors its own specific polymorphisms; therefore, the ideal referent for the genetic analysis of BCG-specific evolution should be a virulent isolate of M. bovis (b). Second, it merits reconsideration whether phenotypic properties generally ascribed to M. bovis, such as preferential growth with supplementation of pyruvate (25), constitutive production of MPB70/83 (29), and virulence in rabbits (30), are general properties of all genomovars of the organism.
A considerable number of genes disrupted through these deletions have predicted functions of potential biological impact, but their relevance for phenotypes such as virulence will clearly require targeted study. As a means of exploring the impact of these deletions, the functional classification of genes that have been deleted can be compared to their proportional distribution in the H37Rv genome (Table S6 in the supplemental material). In contrast to a null model where genetic loss is random across functional classifications, several differences are observed. As would be expected from a highly repetitive sequence, genes classified as "insertion sequences and phages" and "PE/PPE genes" are often implicated in deletion events. In contrast, genes classified as "information pathways" and "intermediary metabolism and respiration" have been relatively spared from deletion. Though deletions from M. tuberculosis sensu stricto (50) also revealed the deletion of mobile genetic elements to be more frequent than expected by chance, a high deletion rate of genes involved in intermediary metabolism and respiration is contrary to what is observed from study here. Thought to confer a selective advantage only during certain stages of infection or transmission ("intermittently essential"), the deletion of genes involved in intermediary metabolism and respiration may not present a viable strategy in the long term for derived MTC. To further explore the functional consequences of deletions, these same genes can be examined in comparison to published lists of genes essential in vitro (46) and in vivo (47) by transposon site hybridization studies (Table S5 in the supplemental material). Of genes essential in vitro, eight were implicated in the deletions characterized here (fadD30, gca, Rv1232c, Rv2813, moaA1, moaC1, moaD1, and Rv3113, where moaC1, moaD1, and Rv3113 are involved in two independent deletions), a somewhat surprising result, as the organisms typed had been isolated by culture in microbiology labs. In addition, five genes deemed essential for in vivo growth were lost in the newly described deletions (Rv1371, Rv1930c, Rv1931c, Rv2812, and Rv3114, where Rv3114 is involved in two independent deletions). Together with the loss of genes Rv1514c, Rv1563c, Rv1974, and cobL (Rv2072c) in deletions previously documented (RD4, RD17, RD7, and RD9, respectively), this indicates that a total of nine in vivo essential genes have been deleted from isolates manifestly capable of causing disease. While these data might be used to challenge the validity of the mouse model for determining genes that impact in vivo pathogenesis, one should consider that although five genes within the RD1 region deleted from BCG vaccines (Rv3871 to Rv3879c) are also considered essential in vivo (47), it is established that natural RD1 deletion mutants cause disease in voles (7) and dassies (34) yet are of attenuated virulence in a number of hosts, including humans (43). Together, these data suggest that the evolution of MTC organisms across different hosts has involved the deletion of genes elsewhere considered essential.
From the genomic comparison of M. bovis AF2122/97 to M. tuberculosis H37Rv, it was noted that the genomic evolution of MTC members appears to parallel the massive gene decay documented in the leprosy bacillus, M. leprae (12, 18). Consistent with this, examination of the ORFs in the newfound deletions reveals that the majority (98%) are absent or no longer functional in the genome of M. leprae (Table S5). The amount of genome loss between M. bovis isolates can also be contrasted with published observations on M. tuberculosis. A similar GeneChip study of 100 M. tuberculosis clones recently quantified the amount of DNA deleted from isolates with respect to the H37Rv genome and estimated the difference in gene content between two M. tuberculosis strains that are not closely related as between 40 and 100 genes (50). From this, it was concluded that the observed magnitude of genomic variability caused by LSPs within M. tuberculosis sensu stricto is similar to that found between M. tuberculosis H37Rv and M. bovis AF2122/97. LSP genomic variability among the M. bovis isolates interrogated here (restricted to samples having deleted RD4) can be quantified in a similar manner by using the sequenced strain M. bovis AF2122/97 as the referent. On average, one region affecting seven ORFs (or 8,483 bp) was deleted per isolate, with the maximum number of deleted genes being 23. Knowing that five genes are missing from regions in AF2122/97 but not in other M. bovis isolates (N-RD17, RDpan, and RD17), we can approximate that the difference in gene content between two M. bovis strains that are not closely related is on the order of 12 genes and is unlikely to exceed 28 genes. Incorporating the overall genomic variability observed across M. bovis strains, these estimates should be irrespective of the M. bovis referent genome chosen for analysis. These estimates place the diversity in gene content among M. bovis sensu stricto as less than that for M. tuberculosis sensu stricto, consistent with M. bovis being a more recent clone having evolved for less time than has M. tuberculosis (8, 35).
The approach to MTC speciation has traditionally been based upon phenotypic distinction, itself viewed as indicative of phylogenetic separation. Because of documented phenotypic convergence across the MTC, which limits the ability of in vitro tests to predict MTC genotype (37), reliance is increasingly being placed on molecular divergence from known species for diagnostic work (23, 42). With this change, standard phenotypic tests may serve as a means of describing certain physiological properties of samples, while PCR-based testing for genomic deletions such as those described here can become a simple and straightforward means of rapidly classifying organisms. With this comes the ability to address the generally held belief that, unlike the principally human-specific M. tuberculosis, M. bovis universally infects and is able to sustain disease across a broad range of animal species. Our data, from a limited number of isolates, suggests this may be partly true. In Canada and South Africa, isolates sharing the same genomic profile caused disease in a variety of hosts, pointing to geographically successful clones that are not host restricted at least in terms of establishing infection, although the capacity of M. bovis to persist in these alternative host types has been inconsistent and requires further investigation. In contrast to these cross-species events, our data confirm the genomic distinction of MTC organisms from pinniped and caprine hosts and suggest other distinct clades that may be host specific, notably, isolates from oryx and the East Asian water buffalo. Further study using these genomic markers and others across a broader sample of clinical isolates will illuminate a greater appreciation of the diversity of M. bovis organisms and, in turn, enable a more focused understanding of the challenges faced for TB control.
This work was supported by Canadian Institutes of Health Research (CIHR) grant MOP 36054 and the Department of Environment, Food, and Rural Affairs (United Kingdom). S.M. is the recipient of a studentship award from the Fonds de la Recherche en Santé du Québec. M.A.B. is a New Investigator of the CIHR.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»