Eukaryotic Cell, July 2008, p. 1075-1084, Vol. 7, No. 7
1535-9778/08/$08.00+0 doi:10.1128/EC.00062-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
| MINIREVIEW |

Aberdeen Fungal Group, Institute of Medical Sciences, Aberdeen, United Kingdom
|
|
|---|
The ability to discriminate between microbial strain types has been, for very many years, a topic of great importance to epidemiologists. Strain typing facilitates tracing of sources and routes of transmission of infection outbreaks. However, a paradox exists at the heart of strain typing. The question most often posed by an epidemiologist is whether two microbial isolates are the same type, but the methodology used for typing is designed to reveal differences between strain types and is usually assessed statistically for its success in differentiation. Identity is inferred by implication when two isolates cannot be distinguished; however, increasing technological finesse means that two isolates classed as identical types in one system may be shown to differ in another. One possible resolution to the paradox of proving strain identity would be to compare whole-genome sequences for two isolates, something that is already a technical possibility for selected isolates and which seems likely to become a more common event in the near future. However, because sequence differences may arise rapidly in highly variable regions of the genome, evidence for strain identity may never achieve its theoretical perfection.
MLST, while a long way from providing complete genome sequences, generates evidence for similarities and differences between isolates from sequences determined, typically, for six to eight independent chromosomal genes. The MLST approach, originally devised by Maiden and colleagues (41), offers a number of advantages over most earlier phenotypic and genotypic typing methods. Sequencing data offer a higher level of reproducibility and minimal subjectivity in analysis compared with other technologies, and the sequences can be stored in Internet databases, offering an unprecedented degree of portability and accessibility to all interested users. (MLST databases can be accessed at two main host sites, http://pubmlst.org/ and http://www.mlst.net/.) Moreover, MLST is based on the same conceptual premise as multilocus enzyme electrophoresis, allelic variation, which means that MLST provides information relevant to population genetic studies as well as to epidemiology. The field of Candida epidemiology has already been substantially explored by a range of technologies. MLST data confirm the many previous studies in the field and offer an increment of refinement in the phylogenetic detail with which epidemiology can be explored experimentally.
Strain typing as a routine clinical laboratory procedure is more familiar to bacteriologists than to mycologists, which is perhaps why, with the exception of Aspergillus fumigatus (2), MLST data only for Candida species and Cryptococcus neoformans have been made available on the central MLST databases. In fact, the MLST approach to strain delineation was used with Coccidioides immitis (31) at least as early as the year when bacteriologists first set up their Web-based MLST systems (41). MLST has been applied to epidemiological and phylogenetic research with other fungal opportunist pathogens, including Cryptococcus gattii (8, 26), Fusarium solani (74), and Batrachochytrium dendrobatidis (43), but the authors of these publications have not utilized the central MLST Web repository for their data, which makes it less easy for others to add or reanalyze the data. The opposite applies to MLST for Cryptococcus neoformans var. grubii (35), the latest fungal addition to the central MLST database (http://www.mlst.net/).
Genes chosen to be sequenced for bacterial MLST traditionally encoded "housekeeping" functions, i.e., metabolic activities that are subject to stabilizing selection, because diversifying selection may obscure relationships among isolates (40). The SNPs for most of the genes used for MLST with Candida species generate ratios (dN/dS) of nonsynonymous to synonymous amino acid changes in their encoded peptides of less than 1.0, the criterion for stabilizing selection. However, we agree with Taylor and Fisher (68) that the choice of MLST loci needs to provide as much sequence diversity as possible; our scheme for C. krusei (30) includes loci with dN/dS ratios of >1 for precisely this reason.
The choice of genes with potential value for MLST can be made easily when full genomic information is available to help locate potential SNPs. In the absence of genomic information, the approach involves some hit-and-miss testing of a set of putatively distinct isolates by sequencing genes for which information is available, sometimes with the use of degenerate primers designed against conserved regions of likely housekeeping genes. This, for example, was the approach we used to devise the gene set for C. krusei MLST (30). The eight-gene set for typing C. dubliniensis was drawn directly from the set used to type C. albicans but with redesigned primers and, for two genes, with extended lengths of DNA sequenced (42).
In practice, to type an isolate by MLST, DNA is extracted; six (for C. glabrata, C. krusei, and C. tropicalis), seven (for C. albicans), or eight (for C. dubliniensis) PCRs are set up with the appropriate primer sets; and the PCR products are purified and used for bidirectional sequencing reactions (normally done with automated capillary sequencing equipment). The sequencing chromatograms are checked visually and scrutinized for heterozygous SNPs (C. albicans, C. dubliniensis, C. krusei, and C. tropicalis). Each sequence is compared with the reference database for the species. If it matches an existing sequence, it is assigned that sequence's allele number (C. glabrata) or (diploid) genotype number (C. albicans, C. dubliniensis, C. krusei, and C. tropicalis). If it is a new sequence, details are sent to the Web database curator, who assigns the sequence a new allele/genotype number. The set of six to eight allele/genotype numbers defines a strain type, numbered once again by reference to the database or by addition to it. The end product of the process is therefore the sequence type (ST) for haploid organisms or a diploid ST (DST) for diploid organisms. If two isolates have different STs or DSTs, this is evidence that the strain types are different, although the STs/DSTs themselves provide no information on the extent of SNP differences between two sequences. As already explained, if two isolates have the same ST or DST, this is evidence that the strains are indistinguishable by MLST but not necessarily that they are identical strains. Table 1 provides an overview of the current status of MLST data for five Candida spp.
|
View this table: [in a new window] |
TABLE 1. Overview of the contributions of sequenced gene fragments to MLST for five Candida speciesa
|
Strain typing methodology, as the name implies, results in delineation of strain types: multilocus STs in the case of MLST. The term ST (40) can be applied to all organisms; however, as already indicated, the term DST is normally used for diploid species, such as C. albicans, C. dubliniensis, C. krusei, and C. tropicalis. More problematic is the appellation of the individual gene sequences that combine to create the STs or DSTs. For haploid organisms these can properly be referred to as alleles, but for diploid organisms the sequences have been named "genotypes" ever since the earliest publication on C. albicans MLST (6). This inevitably generates a complex nomenclature for the MLST operation, with terms such as "genotype" and "ST" being capable of creating ambiguous meanings when taken out of context.
Diploid species provide an advantage for MLST work, in that heterozygous SNPs add a third level of information to the strain differentiation process. More useful for phylogenetic analyses are "haplotypes," which are the exact biallelic combinations in a diploid strain. Haplotypes can be easily determined for a given gene fragment only when both alleles are identical or when they contain a single heterozygous SNP. With two or more heterozygous SNPs in a sequenced fragment, haplotype determination requires separate cloning and sequencing of each allele. However, with sufficient sequence data available, dedicated computer programs can estimate haplotype probabilities with apparently commendable accuracy.
The Candida literature has for several years employed the word "clade" to describe a cluster of genetically related isolates. In its strict phylogenetic sense, "clade" should refer to a common ancestral species and all its evolutionary progeny and is therefore a less appropriate term to define genetically related groups of strains without any robustly defined evolutionary lineage. However, the word "cluster" is already taken in the term "clonal cluster," which defines groups of strains delimited by software such as eBURST (see below), and can therefore cause ambiguity. In phylograms of analyzed MLST data, the branches at least resemble true phylogenetic clades. We will therefore join the rest of the Candida typing field in continuing informally to "take the liberty" (59) of defining major sets of genetically related isolates as clades.
A phenomenon well known for many years in C. albicans is evolutionary change represented by small alterations in the sequences of some genes. The phenomenon was first described in 1995 (36) and was called "microevolution." Although this name seems to be consistent with the phenomenon described, population biologists take "microevolution" to refer to small-scale changes in allele frequencies, and the term seems to have taken on many shades of meaning in recent years. We have therefore coined the term "microvariation" to describe sequence changes at the level of individual gene fragments (45), and "microvariation" will be used instead of "microevolution" in this review.
500 bp used for sequencing. However, the need for sequencing to be as accurate as possible remains; for MLST both DNA strands are routinely sequenced to ensure that ambiguous base readings can be called as correctly as possible. Prospective evaluations of sequencing accuracy indicated a reproducibility of 99.72% for C. albicans MLST (63) and 100% for C. dubliniensis MLST (42). One problem with sequencing a diploid genome concerns heterozygous SNPs. In many instances a result is easily interpreted as indicating heterozygosity when two peaks in a sequencing chromatogram overlap with equal heights approximately half the size of neighboring peaks. However, the majority of heterozygous SNPs appear in practice with overlapping peaks of less than perfectly equal heights. To avoid erroneous interpretation of double peaks as heterozygous SNPs when they may be no more than sequencing noise, we routinely take account of noise in neighboring positions and carefully compare the chromatograms for both strands sequenced. Our approach is illustrated in Fig. 1. For MLST with a diploid fungus, visual scrutiny of sequencing chromatograms is a highly labor-intensive aspect of the procedure. Inaccurate calls of heterozygosity at a polymorphic site undoubtedly occur and are probably the main source of irreproducibility in MLST data, even though such irreproducibility is demonstrably very low.
![]() View larger version (22K): [in a new window] |
FIG. 1. Sequencing chromatograms for bases 2 to 13 in the ACC1 fragment of C. albicans, indicating heterozygosity. For each isolate the forward and reverse chromatograms are shown. For isolates A and B, position 8 would be read as heterozygous AG, even though the A and G peaks coincide perfectly only in one of the four sets. For isolate C, the very small A peak at position 8 would be read as noise, by comparison with adjacent bases, and the result read as homozygous G for MLST purposes.
|
![]() View larger version (18K): [in a new window] |
FIG. 2. Rates of acquisition of new genotypes and new DSTs for C. albicans MLST. (a) As the number of sequenced isolates increases, the rate of discovery of new genotypes falls off fairly rapidly (pink, VPS13; green, ZWF1b; brown, SYA1; black, AAT1a; blue, ADP1; red, PMIb; yellow, ACC1). (b) The rate of acquisition of new DSTs falls off much more slowly, reflecting the high discriminatory power of the method.
|
Consider two C. albicans isolates with the same sequences for six of the seven MLST genes but with the seventh gene in one isolate containing, say, seven heterozygous SNPs while the other isolate has no heterozygous SNPs in that gene. This sort of situation is very common in the MLST database, and indeed loss of heterozygosity has been recognized for many years as a means by which C. albicans can generate microvariation and population diversity (15, 25, 48). At the level of the seven gene fragments sequenced for MLST, the difference between the two isolates is 1/7 (14.3%). At the level of nucleotide polymorphisms, the difference is 7 SNPs out of 185 total known SNPs across the seven genes sequenced (Table 1), or only 3.8%. Two isolates would need to differ in 26 of the 185 SNPs to appear as different as two isolates with 1/7 gene fragments different. We are unaware that analysis at either the locus or the SNP level is more "correct" in the context of MLST. In multilocus enzyme electrophoresis, differences would appear only at the equivalent of the allele level, whereas it can be argued that MLST allows greater finesse by consideration of numbers of SNP differences. The problem of comparing SNPs becomes further complicated by the diploid genomes of many Candida species. For clustering purposes, homozygous SNPs (e.g., CC versus TT) have been scored as twice the value of heterozygous differences (e.g., either CC or TT versus CT) (6, 44, 63). This reduces the calculated interisolate difference for the seven-heterozygous versus seven-homozygous SNP example to just 2%.
Analysis of strain differences based on whole-genotype or allele differences is relatively straightforward. The software called eBURST was written precisely for this purpose (23). It analyzes sets of genotype (allele) data and joins isolates that differ only in one of the whole set of genes sequenced in a "clonal cluster." The output from an eBURST analysis for an entire MLST data set is visualized in the form of a "snapshot" of (D)STs (Fig. 3) in which the structures of clonal clusters and the prevalence of "singletons" (isolates that do not belong to any cluster) can be reviewed. Although the default setting for an eBURST analysis is to generate clonal clusters in which (D)STs are linked to others that differ in a single allele, it is permissible, and sometimes more revealing, to use a more relaxed criterion for a cluster, such as differences in two alleles. The theory behind the eBURST software is that the putative founder of a clonal cluster is the most probable ancestor from which all other linked types have evolved. Recent research has shown, however, that while this premise appears to hold well for many organisms for which the snapshot shows few singletons, the snapshots for the three diploid Candida species, with their high prevalence of singletons (Fig. 3), are typical of populations with a high rate of recombination relative to mutation and for which eBURST does not reliably indicate ancestry (71). This finding explains why eBURST applied to a small population of multiple C. albicans isolates from the same set of patients entirely failed to correspond to likely ancestral relationships (45). For the three diploid Candida species, then, eBURST is a helpful clustering tool but otherwise does not contribute to phylogenetic analyses.
![]() View larger version (45K): [in a new window] |
FIG. 3. Computer printouts of eBURST snapshots for 1,564 isolates of C. albicans (a), 126 isolates of C. krusei (b), and 242 isolates of C. tropicalis (c). In each snapshot, lines join pairs of isolates that show sequence differences in only one of the panel of genes used for MLST. The positions of the clonal clusters and singletons and the lengths of lines are without relevance. For all three species, the snapshots indicate several clusters of related strains but also a high proportion of singletons, indicative of a high rate of recombination relative to mutation (70).
|
It may reasonably be objected that the basis used for delineation of C. albicans clades, i.e., an arbitrary p distance cutoff, has no objective reality. (The cutoff of 0.04 used for C. albicans [44, 63] was chosen because it separates clades 2 and 4, which are well discriminated by other typing methods.) We acknowledge that a small minority of isolates can shift to different clades as the data set grows and also that very small clades may be assimilated into others. However, the assignments of isolates in the largest clade groupings appear to be robust and stable and have not altered as the size of the database has increased.
Bootstrap values for cluster nodes in the UPGMA dendrogram have proved unhelpful in delineating clade boundaries because they become very low in large data sets of (mainly) closely related strains, except for the smallest, tightest subclusters. We have attempted to use more robust phylogenetic clustering approaches for MLST data, such as maximum-parsimony and minimum-likelihood trees, but these methods require very long computational times or even fail, with some software packages, when applied to large numbers of isolates. Others using these methods have applied them only to small sets or subsets of isolates, typically fewer than 120 (8, 26, 31, 35). When larger isolate sets have been analyzed, e.g., the 218 F. solani isolates analyzed by Zhang et al. (74), only three sequenced loci were tested. We too have run maximum-parsimony analyses with the seven MLST gene loci run separately with small subsets of the C. albicans panel that represent the strain diversity found experimentally (44), but that approach involves our selecting isolates by their positions in the UPGMA dendrogram, thus defeating the object of establishing clades by objective criteria.
PHASE-reconstructed C. albicans haplotypes were analyzed for possible recombination events detectable by comparison of all pairs of SNPs in each haplotype. These occurred often enough within the species that its propagation cannot be regarded as exclusively clonal (44). This finding accords with the high prevalence of singleton strains in eBURST snapshots (71). In a related analysis the data for seven genes sequenced for C. albicans MLST showed a lack of congruence suggestive of separate evolutionary histories for each gene (44). The phylogenetic data therefore clearly point to significant levels of recombination, and even possibly sexual matings, in the evolutionary history of C. albicans.
For C. dubliniensis, where so far only 50 isolates have been typed, the data indicated a purely clonal mode of reproduction (42). One notable feature of the first C. dubliniensis MLST study is the relative lack of SNP diversity in this species compared with C. albicans (42). From Table 1 it can be seen that more than 6% of the bases sequenced for C. albicans are SNPs, compared with fewer than 1% for C. dubliniensis. Even when only the first 86 isolates of C. albicans had been typed, 2.9% of the bases sequenced were SNPs. It is not possible similarly to combine and compare MLST data for other Candida species, since they are based on sequencing of different genes. However, for C. parapsilosis, the near-total absence of SNPs was interpreted as indicting a relatively recently evolved species (64); perhaps the low proportion of SNPs among the C. dubliniensis genes sequenced for MLST also points to an evolutionary origin more recent than that for C. albicans. However, on the evidence from comparisons of nuclear DNA sequences (17) and from whole-genome phylogenies (24), it seems more likely that both species arose from a common ancestor. The fact that the positions of the MLST SNPs in the same genes from the two species are mostly different (42) also argues for a common ancestor rather than sequential origin. Like C. albicans and C. dubliniensis, C. guilliermondii and C. fermentati are two very closely related species. Isolates of the two species could be easily differentiated by sequencing five genes, although minimal intraspecies variation was seen, possibly because the genes sequenced were likely to be highly conserved (32).
The principal observation common to MLST surveys of C. albicans and C. glabrata is the tendency for isolates from individual geographical areas to cocluster within different clades. This leads to notably high or low representation of particular geographical sources within some clades, which is a predictable phenomenon, since strains given time to evolve in separate geographical compartments will inevitably develop locally associated traits. However, because most of the pathogenic Candida species are carried as human commensals, movements of human populations, particularly those that have occurred during the most recent centuries, are likely to have contributed to admixture of strain populations across the planet. Thus, among 103 isolates of C. glabrata, those from Europe were significantly underrepresented in MLST group III while isolates from Japan were significantly overrepresented in MLST group IV (19).
Among C. albicans clades, clade 2 contains a very high proportion of isolates from the United Kingdom, while clade 11 is enriched with isolates from continental Europe. Isolates from Pacific Rim countries dominate clades 14 and 17 (44). Only clade 1 isolates show a fully worldwide distribution, although even among clade 1 strains, subclusters with geographical specificity can be seen. This is depicted in Fig. 4, where the isolates comprising C. albicans clade 1 are seen to contain two eBURST clonal clusters, one with only 6% of isolates from Asian countries and the other with 78% of Asian isolates. A subset analysis of the C. albicans MLST data based only on isolates from Japan, England/Wales, and North America confirmed clonal clusters with geographical specificity, emphasizing in particular the high prevalence of Japanese isolates in clades other than 1 to 4 (62). It is entirely possible that the perception of clades 1 to 4 as the "major" ones is an artifactual consequence of the high numerical predominance of strains of non-Asian origin in C. albicans typing surveys.
![]() View larger version (13K): [in a new window] |
FIG. 4. Enlarged view of the clustering of 462 isolates in C. albicans clade 1 as defined by a UPGMA dendrogram drawn for 1,433 isolates currently in our database. Two eBURST clonal clusters that constitute the majority of clade 1 isolates are indicated; the second cluster is heavily enriched with isolates from southeast Asia.
|
Other associations of isolate properties with MLST clades include examples of resistance to antifungal agents. The majority of C. albicans isolates with resistance to flucytosine have been found in Ca3/MLST clade 1 and, so far, all clade 1 isolates resistant to flucytosine have had a mutation (R101C) in FUR1, i.e., a common mechanism of resistance within the clade, while flucytosine-resistant isolates from other clades have not shown this mutation (18, 63). Similarly, a higher number of isolates resistant to amphotericin B was found in Ca3 clade SA, which corresponds to MLST clade 4 (4). In C. dubliniensis, a clade originally defined by DNA fingerprinting and confirmed by MLST was enriched with flucytosine-resistant isolates (1, 42).
Clade-associated properties extend to other features unrelated directly to the genes sequenced for MLST and therefore indicating evolutionary codivergence in a species with predominantly clonal propagation. Significant differences between clades have been found for alleles of tandem repeat sequences in the genes ALS3, ALS5, and ALS6, which encode cell wall adhesins (47, 75). Recent work in our laboratory extends this list of clade-associated tandem repeat alleles to ALS1 and ALS7 and also other genes encoding glycosylphosphatidylinositol-anchored proteins, including HYR1 and HYR3 (D. M. MacCallum et al., unpublished data).
Analysis of a European subset of C. albicans isolates (the subset was chosen to minimize geographical influence on clade differences) showed a significantly higher proportion of isolates from superficial infections and commensal carriage compared with blood isolates in clade 1 compared with the four other largest clades (44). This raises the intriguing proposition that the apparent numerical dominance and pan-global distribution of clade 1 C. albicans may be a consequence of slightly higher fitness to colonize mucosal surfaces. Schmid and colleagues have found a "general-purpose" set of isolates, defined by Ca3 typing, whose members have a greater propensity than other types to cause infections (27, 53). This set of isolates seems to correlate largely with clade 1. A small statistical advantage for clade 1 isolates as colonizers or infecting strains is compatible with their numerical dominance among the whole population.
Two small but well differentiated sets of C. albicans isolates merit special mention. The first (clade 13 [44]) comprises almost entirely isolates that were originally described as a new species, C. africana (69). The main characteristics of these isolates were their predominantly (but not exclusively) African origin and their genital, typically vaginal, site of isolation. Among the 16 isolates presently included in clade 13, all from different patients and from 10 different countries, only three DSTs are represented (14 isolates are DST 182), 14 isolates are from genitalia (one other comes from sputum and one from blood culture), and the most dissimilar isolate in the clade, which is still within a p distance of 0.02, comes from Japan (reference 44 and unpublished data). The second small strain set of interest comprises only four isolates, all of them DST 1031 and all originally identified as C. stellatoidea type 1, i.e., strains unable to assimilate sucrose, including CBS1905, the C. stellatoidea type strain (29a). It is perhaps worth noting that the "C. stellatoidea" isolates join with five singleton isolates and with 13 isolates designated clade 16 to form a set of isolates more dissimilar than any others from the "mainstream" clades of C. albicans. They have extremely diverse geographical origins.
MLST has similarly confirmed consistent intraindividual strain types, microvariation, and occasional strain replacement (5, 10, 45). C. albicans strain types were known to be transmissible between sex partners (54), and MLST has extended the range of transmission to include families (5).
Loss of genetic heterozygosity as a mechanism underpinning microvariation in C. albicans (15, 25, 52, 66, 70) is inevitably picked up also by MLST (5, 45). However, one striking feature of MLST analyses of successive isolates from the same patient is the finding of loss of heterozygosity at one or more sequenced loci, followed by regain of heterozygosity (45). While loss of heterozygosity can be readily accounted for by phenomena now well demonstrated in C. albicans, i.e., chromosomal loss and reduplication or mitotic recombination (25, 33, 56, 58, 70, 73), the regain of heterozygosity suggests that populations colonizing or infecting an individual may contain mixtures of heterozygous and homozygous variants (45). C. albicans colonies with different morphologies plated from a single sample were shown to have identical or similar Ca3 fingerprinting patterns (28). We recently typed by MLST randomly picked, separate C. albicans colonies from primary isolation plates and found that the rates of zygosity variation (and even of strain replacement) within a sample are similar to those seen between samples when single clones are typed from sequential isolations of the fungus from a single patient (Jacobsen et al., submitted). All these findings combine to suggest that populations of C. albicans even within the same site in an individual patient may be mixtures of strain types.
Finally, an MLST study of C. albicans isolates from a variety of animals has confirmed and extended earlier work, based on Ca3 fingerprinting (20), showing that strains colonizing or infecting animals are related to those found in human samples but with clear evidence of significant differences, suggestive of a separate evolution in the different hosts (Jacobsen et al., submitted).
It is possible to view the constant addition of new isolates to the now-large C. albicans database as a biological resource, but in reality the steady increase in numbers of strains typed has ultimately contributed more to the statistical robustness of conclusions drawn than to the creation of new conclusions. The databases for the other typeable Candida species are much smaller than that for C. albicans, so these will benefit from further growth. Isolates of C. albicans that are indistinguishable by MLST can be differentiated by other approaches such as microsatellite typing (Jacobsen and Odds, unpublished data), and SNP arrays for typing Candida strains based on MLST sequences have already emerged (38) in the same way as they have for other organisms, from bacteria to mammals. The present MLST scheme may need this type of further adaptation to be of clinical utility, for example, to identify putative outbreaks of C. albicans infection within an institution. In a survey of candidemia in which all C. albicans isolates were typed by MLST, indistinguishable DSTs were found for some bloodstream isolates from hospitals in different cities (46). This is unlikely to be evidence for interhospital strain transmission; it probably is just a reminder of the difficulty mentioned at the start of this review in proving that isolates are the same.
We owe a particular debt of gratitude to (in alphabetical order) Jude Bain, Amanda Davidson, Arianna Tavanti, and Julie Whyte, who also undertook all the tedious work of sequencing, and to Duncan Shaw, whose facility with genetics and statistics has contributed hugely to our analyses of MLST data. Our colleagues Marie-Elisabeth Bougnoux, Christophe d'Enfert, and Shu-Ying Li, in Paris and Taiwan, have all been true pioneers, generating masses of data and shrewd observations in the Candida MLST world. Finally, thanks are due to Neil Gow and Martin Maiden, whose encouragement prompted our laboratory's research into typing fungi by multilocus sequencing.
Published ahead of print on 2 May 2008. ![]()
|
|
|---|
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»