Previous Article | Next Article ![]()
Journal of Clinical Microbiology, March 2006, p. 783-789, Vol. 44, No. 3
0095-1137/06/$08.00+0 doi:10.1128/JCM.44.3.783-789.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Seattle Biomedical Research Institute, Seattle, Washington,1 McGill University Health Center, Montreal, Canada,2 U.S. Environmental Protection Agency, National Exposure Research Laboratory, Cincinnati, Ohio,3 National Institute of Public Health and the Environment, Bilthoven, The Netherlands,4 Departmento de Microbiologia, Immunologia e Parasitologia, Universidade Federal de São Paulo, Escola Paulista de Medicina, UNIFESP-EPM, São Paulo, Brazil5
Received 2 November 2005/ Returned for modification 14 December 2005/ Accepted 28 December 2005
|
|
|---|
|
|
|---|
To assess the clinical incidence of infections involving strain 104, we examined a sample of clinical isolates of M. avium from Southern California and elsewhere for genotypic identity with the genome sequence strain. Although these isolates were initially classified as M. avium subsp. avium, new data suggest that true M. avium subsp. avium strains infect a narrow spectrum of avian hosts and are described genotypically by a characteristic restriction fragment length polymorphism with IS1245, by the presence of IS901, and by a limited temperature range in which they can be grown (15). Mijs and colleagues recommended the subspecies name "hominissuis" for the more genetically polymorphic M. avium strains that grow in a wider temperature range of 24 to 45°C and classically infect pigs and humans. The M. avium isolates used in this study are phenotypically consistent with M. avium subsp. hominissuis (15). However, because most of the isolates were not assessed to determine whether they have the genotypic characteristics of M. avium subsp. hominissuis, they are referred to in this study by the umbrella term "M. avium."
M. avium isolates are known for their genotypic polymorphism (1-3, 8, 12, 19, 20, 24). High-resolution pathogen genotyping approaches, such as pulsed-field gel electrophoresis and restriction fragment length polymorphism (RFLP) with IS1245, have rarely linked MAC isolates outside of geographically restricted samples. Because these methods are time consuming, labor intensive, and not portable (i.e., difficult to standardize between laboratories), few molecular epidemiologic studies of MAC have been conducted that encompass broader geographic samples. Recently, a commercial repetitive sequence-based PCR (rep-PCR) system was shown to deliver resolving power equivalent to that of RFLP with IS1245 (6). This method is rapid and reproducible enough for analysis of large samples, but its cost can be restrictive. To stratify strains of interest for high-resolution typing, preliminary analysis with a moderate resolution typing method can be cost effective. For example, a two-tier approach was used for many years in molecular epidemiological analysis of Mycobacterium tuberculosis. Spoligotyping was used to presumptively link isolates prior to high-resolution typing by the more labor-intensive RFLP method (4, 10, 26). To search a large strain sample for genotypic identity to strain 104, a high-throughput tool was needed that could be applied like spoligotyping to presumptively link isolates prior to high-resolution typing.
Sequencing and annotation of multiple mycobacterial genomes has helped to identify large sequence polymorphisms (LSPs), defined as insertions or deletions comprising at least one ORF. Varying patterns of LSP occurrence offer an alternative method of typing mycobacteria in a quick, PCR-based method (7, 9). LSPs are thought to be relatively stable and unique events, and they have been used to identify genetic linkages within and between Mycobacterium species (11, 16, 25). As a simple, inexpensive method to interrogate strains for stable polymorphisms, LSP typing by PCR can streamline the molecular epidemiological analysis of a large sample, especially when the sample set is screened for isolates that are genotypically related to a reference strain. It can be expected that isolates that share the reference strain's LSP type are more likely than other isolates to share the reference strain's high-resolution DNA fingerprint (e.g., IS1245 RFLP or rep-PCR patterns). We used this strategy to screen a large clinical isolate sample set for epidemiologically unlinked isolates that share genotypic identity with the genome sequence strain M. avium strain 104.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. M. avium subsp. avium clinical isolates (n = 208)
|
LSP typing PCR strategy.
The 4 LSPs examined in this study were designated LSP2, LSP7, HSD, and LSPP5. This nomenclature was consistent with prior publications (14, 22, 23). Each consisted of a genomic region of
2.5 kb that was present in some M. avium strains and absent in others. Criteria for their selection are described in Results. PCR primers used to interrogate genomic DNA for the presence or absence of each region are summarized in Table 2. Genomic DNA from each isolate was interrogated with 4 multiplex PCRs, each of which tested for one LSP.
|
View this table: [in a new window] |
TABLE 2. PCR primers used for LSP typing
|
The presence of the HSD region was detected in a similar manner. When the HSD region was present, an amplicon was primed by the LFP and LIP, as in LSP2 and LSP7. When the HSD region was absent, an amplicon was primed by the LFP and an opposing primer that spanned the deletion junction (bridging primer).
Preliminary analysis of LSPP5 found multiple polymorphisms in both border regions. Therefore, a set of 4 primers was used to detect the presence or absence of this LSP. When the region was present, two internal primers (LSPP5 LIP and right internal primer) amplified an internal segment. When it was absent, flanking primers (left and right) amplified across the junction as in LSP2 and LSP7.
Amplicon sizes expected from the multiplex PCRs testing for the presence and absence of each region were as follows: LSP2, 954 bp (present) and 834 bp (absent); LSP7, 181 bp (present) and 329 bp (absent); HSD, 480 bp (present) and 343 bp (absent); LSPP5, 794 bp (present) and 728 bp (absent).
PCR conditions. Except for HSD, the PCR was performed by combining 5 µl of 5x GC-rich PCR buffer, 2 µl of 10 mM nucleotides, 1.5 µl of 10 mM primer mix, 2.5 µl of GC-rich resolution solution (Roche Diagnostics, Basel, Switzerland), 1 µl of 10 µg/ml template, 1 U of Taq enzyme, and sterile water to make 25 µl per reaction mixture. The HSD PCR mixture required the addition of 0.6 µl of 2.5 mM MgCl2, which lowered the GC-rich resolution solution to 2 µl per reaction. PCR was performed with a PTC-200 Peltier thermocycler (MJ Research, Boston, MA). Both LSP7 and LSPP5 used an initial denaturation step of 95°C for 3 min, followed by 36 cycles of 1 min at 58°C, followed by 1 min at 72°C. PCR conditions for HSD began with an initial denaturation step at 95°C for 3 min followed by 36 cycles of denaturation at 95°C for 45 s, annealing at 55°C for 1 min, and elongation at 72°C for 1 min. PCR conditions for LSP2 included an initial denaturation step at 94°C for 3 min, followed by 36 cycles of 45 s each at 94°C, 57°C, and 72°C. Amplicons were analyzed by 1% agarose gel electrophoresis with ethidium bromide staining in comparison to a 1-kb ladder (Sigma-Aldrich Corp., St. Louis, MO). Sizing was confirmed by the ladder as well as by correspondence to well-characterized positive- or negative-control strains (M. avium strains 104 and HMC02, and M. avium subsp. paratuberculosis strain K10).
LSP typing nomenclature. An amplification result corresponding precisely in size to the presence of a region was denoted with the number "2." A result corresponding precisely in size to the absence of a region was denoted with the number "1." If no amplification or anomalous bands were observed, the locus was assigned an "X" (indeterminate). Thus, each isolate was assigned a 4-digit LSP type, with results for each locus ordered as follows: HSD, LSP2, LSP7, LSPP5. For example, a strain lacking all 4 regions was denoted "t1111," whereas a strain with genomic material at LSP2 and LSP7, but not at the other two LSPs, was denoted "t1221."
High-resolution strain typing by rep-PCR. The DiversiLab System (Bacterial Barcodes, Athens, GA) was used for high-resolution strain typing, following the manufacturer's protocols. rep-PCR takes advantage of repetitive elements interspersed throughout bacterial genomes. When amplified by PCR, these elements generate highly discriminative genomic fingerprints. The DiversiLab Mycobacterium kit has been shown to deliver resolving power equivalent to that of restriction fragment length polymorphism targeting IS1245 in MAC (5). rep-PCR required 2 µl of 25 to 50 ng/µl genomic DNA. Dilutions were performed with sterile water.
|
|
|---|
To avoid bias toward regions that are present in strain 104, two additional regions were identified that are absent in strain 104 relative to other strains. One such region, HSD, was characterized in a previous genetic analysis of a Seattle clinical isolate, HMC02 (14). It is
2,631 bp long and codes for a cytochrome P450-like gene and two type I restriction-modification genes, hsdR and hsdM, which are disrupted by insertion sequences in some clones. In strain 104, there is a deletion of this region that starts within the P450-like gene and extends beyond hsdR and hsdM to a right terminus that has not yet been mapped. Preliminary analysis showed the mapped terminus within the P450-like gene to be well conserved with almost no polymorphism between strains that carry the deletion.
The fourth LSP, LSPP5, was identified for this study by in silico analysis. It is a 14,471-bp region that was found by CROSS_MATCH analysis to be present in M. avium subsp. paratuberculosis strain K10 but absent in strain 104. It carries 15 ORFs, including multiple transposase genes, indicative of insertion elements. As preliminary PCR analysis revealed the bordering regions to be variable between M. avium strains, PCR primers were designed to hybridize to the conserved internal and external regions within and outside the variable border regions.
The 2 LSPs present in strain 104, LSP2 and LSP7, mapped at positions 5,240,295 to 5,263,079 and 3,210,765 to 3,242,533, respectively, in the most recent draft of the M. avium 104 genome sequence (www.tigr.org). The genomic region that amplified between the flanking primers in LSPP5 corresponded to positions 2,622,423 to 2,623,150 in the 104 genome. The HSD deletion junction mapped to base pair 1,537,166.
Frequency of precise LSP types. Strain 104 plus 207 additional human clinical isolates of M. avium were evaluated by PCR for the 4 LSPs. Each isolate was assigned an "LSP type" that consisted of 4 binary scores for the presence (coded as "2") or absence (coded as "1") of the 4 loci. The scores were listed in the order HSD, LSP2, LSP7, and LSPP5. Of the 16 possible LSP types, 14 were found in the sample. Of the 208 isolates, 174 (84%) yielded amplicons at all four sites of the precise sizes expected for "present" or "absent" results (Fig. 1). The other 34 strains yielded indeterminate results at one or more loci. Most indeterminate results consisted of failure to amplify at one or more loci or the occurrence of multiple bands at one or more loci. Indeterminate results occurred relatively frequently at LSP2 (7.2%) and LSPP5 (6.3%) and less frequently at LSP7 (4.8%) and HSD (3.8%). The largest geographical sample sets had the greatest variety of LSP types. The Southern California and The Netherlands samples exhibited 10 and 8, respectively, of the 16 possible LSP types.
![]() View larger version (11K): [in a new window] |
FIG. 1. Distribution of LSP types. The number of isolates with each LSP type is shown. Blanks indicate that no isolates were found with a given LSP type at a given site. The t1221 LSP types are highlighted by the box.
|
Identification of t1221 isolates. Strain 104 exhibited the t1221 LSP type. In total, 19 t1221 isolates were identified, 12 (63%) of which came from Southern California, a region that provided 75 (36%) of the 208 isolates examined. Of the remaining 7 t1221 isolates, 5 came from Seattle, and 1 (each) came from The Netherlands and Brazil (Fig. 1). The United States t1221 isolates originated from 4 of the 6 Southern California sites and 1 of the 2 Seattle sites. In addition to the 19 t1221 strains, a single indeterminate strain from Southern California, W214 (t1X21), had an LSP type identical to strain 104 at 3 of the 4 loci but failed to produce an amplicon for LSP2.
High-resolution linkages within LSP types.
To identify t1221 isolates that were genotypically identical to strain 104, the 18 additional t1221 strains plus the t1X21 indeterminate strain W214 were typed by rep-PCR, a high-resolution method (5). Nine t1221 isolates exhibited
92% similarity to strain 104 (Fig. 2), a level of similarity that was previously established as the threshold for strain identity within MAC (5). Strain W214 also fell into this cluster. The remaining eight t1221 isolates, including isolates from Brazil and The Netherlands, showed <92% similarity on rep-PCR.
![]() View larger version (43K): [in a new window] |
FIG. 2. rep-PCR analysis of t1221 isolates and other clinical isolates. The dendrogram and gel-like images were generated by the DiversiLab software, which also assigned numbers (1 though 24) to the samples in the set. The strain identifications, LSP types, and dates and locations of collection were added by us, as was the vertical dotted line that marks the 92% similarity cutoff. The Southern California sample collection dates spanned a decade, with the 100 and 500 series collected in the mid-1980s. The 104-like cluster included isolates numbered 5 through 15.
|
To assess the occurrence of the 104-like rep-PCR type in strains that did not share the t1221 LSP type, rep-PCR patterns of 65 non-t1221 strains, representing 11 of the other 13 LSP types that occurred within the sample, were compared to that of strain 104. These strains had been typed in separate studies (40 strains) or were chosen at random (25 strains) for this comparison. None of these strains exhibited
92% similarity to strain 104 (Fig. 3).
![]() View larger version (41K): [in a new window] |
FIG. 3. rep-PCR analysis of non-t1221 isolates compared to strain 104. The dendrogram and gel-like images were generated by the DiversiLab software. The strain identifications, LSP types, and 92% similarity cutoff line were added by us, as for Fig. 2. Strain 104 is isolate number 14 on the figure.
|
|
|
|---|
The 17-year time span was an unexpected finding. We are not aware of prior reports of this nature, and there are several possible models to explain it. A single environmental source of infection may exist, and the 11 patients may have acquired their infections from that single source at different time points over 17 years. Alternatively, the 104-like strain may have stably colonized multiple environmental sites in the region, resulting in infection of epidemiologically unlinked patients. A third scenario is that all 11 patients were exposed to a single point source in the early 1980s and became stably but asymptomatically colonized. Subsequently, over the course of the ensuing 17 years, they may have become immunocompromised and developed MAC disease at differing times and locales. These divergent models highlight an unanswered question in MAC disease. Namely, is the pathogen acquired early in life and asymptomatically carried until host immunity fails, or is MAC acquisition transient and cleared unless a preexisting susceptibility exists? To answer these questions and to better define models of MAC acquisition, large numbers of clinical and environmental isolates will have to be analyzed from multiple sites. The PCR-based LSP typing would greatly facilitate such a study.
LSP typing successfully narrowed the search for the 104-like genotype, as indicated by the fact that 104-like isolates were common among t1221 and t1X21 strains but were not found in a sample of 65 strains with other LSP types. However, it was also noted that the non-t1221 sample included strains that appeared similar or identical to each other by rep-PCR but had differing LSP types (Fig. 3). This observation indicates that LSP types can, in some instances, segregate rapidly and independently of other genotypic parameters.
In contrast to a previous LSP-based study conducted on M. tuberculosis (7), which targeted LSPs that were unique to a single Beijing family outbreak strain, the LSPs used in this study of M. avium strains were broadly heterogeneous among isolates from 3 continents, as evidenced by the variety in LSP types observed. Therefore, the specific LSP tests used in this study can be applied to any molecular epidemiologic analysis of clinical M. avium.
In theory, both the resolving power and the throughput of LSP typing can be further improved. Analysis of additional LSPs will improve the resolution of the test, and probes for these LSPs can be incorporated into filter-based, high-throughput formats similar to that described for "deligotyping" of M. tuberculosis (9). Alternatively, these four LSP loci could be combined with the recently described M. avium variable number of tandem repeats-mycobacterial interspersed repetitive unit loci (21) for a multimodal PCR-based genotyping method. An advantage of LSP typing is that its binary readout is portable, meaning that results can be expressed in spreadsheet form and compared between laboratories. As illustrated by the identification of 104-like clones, LSP typing can be used as a high-throughput, moderate-resolution tool to expedite searches for shared genotypic identity.
Genome sequences were made available by www.tigr.org. We are indebted to Clark Inderlied for providing some of the Southern California isolates and sharing information about their derivation. The United States Environmental Protection Agency, through its Office of Research and Development, collaborated in the research described here. It has been subjected to Agency review and approved for publication.
|
|
|---|
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»