Previous Article | Next Article ![]()
Journal of Clinical Microbiology, November 2006, p. 4125-4135, Vol. 44, No. 11
0095-1137/06/$08.00+0 doi:10.1128/JCM.01231-06
U.S. Department of Agriculture, Agricultural Research Service, Western Regional Research Center, Produce Safety and Microbiology Unit, 800 Buchanan Street, Albany, California 94710
Received 15 June 2006/ Returned for modification 12 August 2006/ Accepted 18 August 2006
|
|
|---|
|
|
|---|
The complete genomic sequence of the C. jejuni strain NCTC 11168 established an origin for understanding C. jejuni genetic differences that would facilitate identification of those determinants that might contribute to the C. jejuni pathogenesis. Using C. jejuni strain NCTC 11168 microarrays, several groups have indexed the complete gene contents of several C. jejuni strains in relation to strain NCTC 11168 (3, 11, 12, 20, 23). These gene indexing studies identified regions of variability between strain NCTC 11168 and other C. jejuni strains, such as the lipooligosaccharide biosynthesis (LOS), capsular biosynthesis (CAP), flagellar modification (FM), and DNA restriction/modification (R/M) loci. Furthermore, DNA sequencing of these strain-specific loci identified variability in both the DNA sequence of common genes and gene composition at each locus (6, 10, 14, 18). Indeed, the complete genomic sequence of a second strain of C. jejuni, RM1221, illustrated the breadth of intraspecies genome diversity (4). The genome of RM1221 is syntenic with the genome of C. jejuni NCTC 11168 but is disrupted by four genomic islands and smaller gene clusters (4). The four genomic islands in strain RM1221 are referred to as Campylobacter jejuni-integrated elements (CJIEs). The first genomic island, CJIE1, located upstream of argC (CJE0275), is a Campylobacter Mu-like phage (also termed CMLP1) encoding several proteins with similarity to bacteriophage Mu and other Mu-like prophage proteins (15). CJIE2 and CJIE4 have several genes predicted to encode phage-related endonucleases, methylases, or repressors and are integrated into the 3' end of arginyl- and methionyl-tRNA genes, respectively. Finally, CJIE3 (integrated into the 3' end of an arginyl-tRNA) may be an integrated plasmid. This is based on the observation that 73% (45/62) of the CJIE3 predicted proteins show sequence similarity to those encoded on the Campylobacter coli RM2228 megaplasmid (4) or other Campylobacter plasmids (2, 17). Also of note, 23% (14/62) of the putative CJIE3 products are similar to proteins found within the 71-kb pathogenicity island of Helicobacter hepaticus (HHGI1) (4). These findings suggest that these unique integrated elements in strain RM1221 may contribute to an additional level of diversity in C. jejuni.
In this study, a comparative genomic analysis was conducted from a collection of 67 C. jejuni and 12 C. coli strains that were obtained from various geographical locations and clinical and veterinary sources. By using a PCR-based assay and DNA microarrays, we demonstrated that the four CJIEs identified in strain RM1221 were present also in other Campylobacter strains that were examined in the present study. The gene indexing analysis performed by using DNA microarrays revealed that several of the coding sequences within these four CJIEs in 26 C. jejuni strains were absent or highly divergent compared to strain RM1221, demonstrating an even greater degree of variability within these integrated elements. Furthermore, we observed that, in several C. jejuni strains possessing the Mu-like CJIE1, CJIE1 was located at sites distinct from its locations within strain RM1221. Together, these data provide greater insights into the degree of genomic diversity and suggest additional genomic regions to be utilized in differentiating C. jejuni strains.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Bacterial strains used in this study and PCR analysis of CJIE contentaa
|
To determine whether the integrated elements detected in other Campylobacter strains were in the same chromosomal location relative to those present in RM1221, genomic DNA from these strains was amplified with primer sets flanking CJIE1 and CJIE3. The flanking primer sets are CJ432 (5'-ATG GGG TAC GAG TGA TGA GTA TAA-3') and CJ504 (5'-TCT TTT CAA ATT CAA AAC TAA AGC-3') for CJIE1 and CJ1230-2 (5'-AAG ATA GCC CTT ATG GTA CAC TTT AT-3') and CJ1299-2 (5'-TTA GGC GTG ATT GCT TAT CTT ATT-3') for CJIE3. The expected amplification product sizes for NCTC 11168 (negative control) are 581 bp for CJIE1 and 681 bp for CJIE3. The genomic DNA of strains containing elements at these positions would not be expected to amplify under standard PCR conditions.
PCR reagents were supplied by Epicenter (Madison, WI). PCR oligonucleotides were purchased from Operon (Alameda, CA). Each amplification mixture contained 50 ng genomic DNA, 1x MasterAmp Taq PCR buffer, 1x MasterAmp Taq enhancer, 2.5 mM MgCl2, 2.5 mM concentrations of each deoxynucleoside triphosphate (dNTP), 50 pmol of each primer, and 1 U MasterAmp Taq DNA polymerase. PCR was performed using a Tetrad thermal cycler (Bio-Rad, Hercules, CA) with the following amplification parameters: 30 cycles of 1 min at 94°C, 2 min at 50°C, and 3 min at 72°C and a final extension at 72°C for 5 min. All PCR samples were analyzed by gel electrophoresis, and positive samples were identified based on the presence of bands of the appropriate size.
Inverse PCR. Genomic DNA (1 µg) was digested to completion with HindIII, and 50 ng of digested genomic DNA was then self-ligated overnight. The primers used for each reaction were mu_inv1a (5'-TTA GCA AAA ATC CCA AAT AAT AG-3') and mu_inv1b (5'-AAT TAT TTA GGG ATA GTT ACA TG-3'). Inverse PCRs were performed on an ABI 9700 thermocycler with the following settings: 1 min at 94°C, 2 min at 50°C, 3 min at 72°C (30 cycles). Each amplification mixture contained 2.5 ng genomic DNA, 1x MasterAmp PCR buffer (Epicenter), 1x MasterAmp PCR enhancer (Epicenter), 2.5 mM MgCl2, 250 µM concentrations of each dNTP, 50 pmol each primer, and 1 U polymerase (Epicenter). Flanking sequences were obtained by sequencing each amplicon using the mu_inv1b primer. The insertion point for each CMLP1 phage was determined by BLASTN analysis against C. jejuni (strains RM1221 and NCTC 11168) sequences.
Construction of the C. jejuni DNA microarray. DNA fragments of individual open reading frames (ORFs) were amplified using the Sigma-Genosys (The Woodlands, TX) C. jejuni ORFmer set specific for strain NCTC 11168 coding sequences and primers from Operon Technologies (Alameda, CA) and designed with ArrayDesigner 2.0 (Premier Biosoft, Palo Alto, CA) specific for unique strain RM1221 coding sequences. Each PCR mixture (total reaction volume, 100 µl) consisted of 1x MasterAmp Taq PCR buffer, 1x MasterAmp Taq Enhancer, 2.5 mM MgCl2, 200 µM concentrations of each dNTP, forward and reverse primers at 0.2 µM each, 0.5 U of MasterAmp Taq DNA polymerase (Epicenter), and approximately 50 ng of genomic DNA (either NCTC 11168 or RM1221). Thermal cycling was performed using a Tetrad thermal cycler (Bio-Rad, Hercules, CA) with the following amplification parameters: 30 cycles of 25 s at 94°C, 25 s at 52°C, and 2 min at 72°C and a final extension at 72°C for 5 min. PCR products were analyzed by gel electrophoresis in a 1% (wt/vol) agarose gel (containing 0.5 µg of ethidium bromide ml1) in 1x Tris-acetate-EDTA buffer. DNA bands were examined under UV illumination. We successfully amplified a total of 1,530 and 227 PCR products from strains NCTC 11168 and RM1221, respectively. These PCR products were purified on a QIAGEN 8000 robot using a Qiaquick 96-well Biorobot kit (QIAGEN, Valencia, CA), dried, and resuspended to an average concentration of 0.1 to 0.2 µg µl1 in 20 µl of 50% dimethyl sulfoxide containing 0.3x saline-sodium citrate (SSC) (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate). All of the PCR probes were then spotted in duplicate on GAPSII slides (Corning) using an OmniGrid Accent (GeneMachines, Ann Arbor, MI), producing a final array that contained a total of 3,514 features.
Preparation and fluorescent labeling of genomic DNA. Genomic DNA from C. jejuni was isolated as described previously (18) or purified using the QIAGEN Dneasy kit (QIAGEN) according to the manufacturer's specifications. For each microarray hybridization reaction mixture, genomic DNAs from the reference strains (an equal amount of DNA from C. jejuni strain NCTC 11168 and C. jejuni strain RM1221) and a test strain were fluorescently labeled with indodicarbocyanine (Cy5)-dUTP and indocarbocyanine (Cy3)-dUTP, respectively. An aliquot (2 µg) of DNA was mixed with 5 µl 10x NEBlot labeling buffer containing random sequence octamer oligonucleotides (NEB, Beverly, MA.) and water to a final volume of 41 µl. This mixture was heated to 95°C for 5 min and then cooled for 5 min on at 4°C. After this time, the remainder of the labeling reaction components were added: 5 µl of 10x dNTP labeling mix (1.2 mM [each] dATP, dGTP, and dCTP; 0.5 mM dTTP in 10 mM Tris, pH 8.0; 1 mM EDTA), 3 µl of Cy3 dUTP or Cy5 dUTP (GE Biosciences, Piscataway, NJ), and 1 µl of Klenow fragment. The labeling reaction mixtures were incubated overnight at 37°C. Labeled DNA was purified from unincorporated label using Qiaquick PCR cleanup kits and dried by vacuum.
Microarray hybridization. Labeled reference and test DNAs were combined in 45 µl Pronto! cDNA hybridization solution (Corning, Corning, NY) and heated to 95°C for 5 min. Then, 15 µl of the hybridization mixture was put onto a microarray slide and sealed with a coverslip. The microarray slide was placed in a hybridization chamber (Corning) and incubated at 42°C for 18 h. Following hybridization, the slides were washed twice in 2x SSC, 0.1% sodium dodecyl sulfate at 42°C for 10 min, followed by twice in 1x SSC at room temperature for 10 min, and finally twice in 0.2x SSC at room temperature for 5 min. The microarray slides were dried by centrifugation at 300 x g for 10 min before scanning. At least two hybridization reactions were performed for each test strain.
Microarray data analysis. Microarrays were scanned and analyzed as previously described by Anjum et al. (1) with modifications. DNA microarrays were scanned using an Axon GenePix 4000B microarray laser scanner (Axon Instruments, Inc., Union City, CA). Features and the local background intensities were detected and quantified with GenePix 4.0 software (Axon Instruments, Inc.). Poor features were excluded from further analysis if they contained abnormalities or were within regions of high fluorescent background. The data were filtered so that spots with a reference signal lower than the background plus 2 standard deviations of the background were discarded. Signal intensities were corrected by subtracting the local background, and then the Cy5/Cy3 ratios were calculated. To compensate for unequal dye incorporation, data normalization was performed as described previously (1). Briefly, the median natural log (ln) (Cy5/Cy3) for each block was brought to 0 (one block being defined as the group of spots printed by the same pin) with the following equation: ln(Ti) = ln(Cy5i/Cy3i) c, where T is the centered ratio, i is the gene index, Cy5 and Cy3 are the Cy5 and Cy3 intensities, respectively, and c is the 50th percentile of all Cy5/Cy3 ratios (1). It should be noted that dye-swapping was performed for 6 strains, and no effects due to differences in dye incorporation on downstream analysis were observed. The comparative genomic hybridization (CGH) analysis of either strain NCTC 11168 or RM1221 versus the reference (NCTC 11168 and RM1221 combined) defined ratio cutoffs for the present, divergent, and absent genes. Based on these data, we defined the status of a gene as present when the Cy3/Cy5 (test/reference) intensity ratio was >0.6, as divergent when the Cy3/Cy5 intensity ratio was between 0.6 and 0.3, and absent when the Cy3/Cy5 intensity ratio was <0.3. These values are similar to intensity ratios utilized in other microarray studies (1, 5, 23). The NCTC 11168 and RM1221 strain-specific spots hybridized to only half of the reference DNA (the Cy5-labeled mixture of NCTC 11168 and RM1221 DNA), increasing the Cy3/Cy5 ratio by twofold. The ratios for these spots were therefore divided by 2 before determining the status of the gene. The presence, divergence, and absence status for all genes was converted into trinary scores (present, 2; divergent, 1; absent, 0). The trinary gene scores for each replicate for all strains were analyzed further with GeneSpring microarray analysis software version 7.3 (Agilent Technologies, Redwood City, CA) and subjected to average-linkage hierarchical clustering with the standard correlation and bootstrapping.
|
|
|---|
Validation of the custom C. jejuni microarrays. A multistrain C. jejuni microarray was constructed based on sequence data from the genome strains NCTC 11168 and RM1221 (4, 16). The DNA microarray comprised 1,530 genes from NCTC 11168 and 227 genes from RM1221 (see Materials and Methods). To validate this DNA microarray in specifically distinguishing the genome strains, a reference DNA mixture containing equal amounts of genomic DNA from both NCTC 11168 and RM1221 was labeled with Cy5 and was then competitively hybridized with either Cy3-labeled NCTC 11168 or Cy3-labeled RM1221 genomic DNA. This CGH analysis of either strain NCTC 11168 or RM1221 versus both strains NCTC 11168 and RM1221 combined (reference) defined the ratio cutoffs for the present, divergent, and absent genes. Based on these data, we defined a gene's status as present when the Cy3/Cy5 (test/reference) intensity ratio was >0.6, as divergent when the intensity ratio was between 0.6 and 0.3, and absent when the intensity ratio was <0.3.
Figure 1 depicts the gene status for all RM1221 genes based on the normalized Cy3/Cy5 fluorescence intensity ratio using strain NCTC 11168 as the test strain (Cy3 labeled). These results demonstrated that a gene present in both strains produced a signal intensity ratio of around 1, as visualized in Fig. 1. These results also confirmed that the genomic DNA of either strain RM1221 or NCTC 11168 hybridized equally well to features on the DNA microarray representing genes common to both strains. In contrast, features on the DNA microarray targeting the regions and genes known to be distinct between strains NCTC 11168 and RM1221 were identified by a low-intensity ratio of <0.3, as visualized in Fig. 1. In previous microarray studies (3, 20, 23), most absent or divergent genes clustered in several genomic regions in the NCTC 11168 genome. Specifically, Taboada et al. (23) defined 16 intraspecies hypervariable genomic regions described in Table 2, including the LOS, CAP, and DNA R/M genes. Interestingly, comparing the genomes of both sequenced strains demonstrated that strain RM1221 possessed genes distinct from strain NCTC 11168 in 11 of these 16 variable regions and was highly divergent from strain NCTC 11168 in 13 of the 16 variable regions. An additional 17th hypervariable region (Cj0258-Cj0263) identified in previous microarray reports (11, 20) and by sequence comparison (4) was also found in our study to be distinct between strains NCTC 11168 and RM1221. As shown in Fig. 1, the results obtained in our DNA microarray analysis correlated well with sequence analysis (4) and verified 14 regions in strain RM1221 that were distinct from strain NCTC 11168, such as the CAP, LOS, and FM regions. Additionally, our DNA microarrays verified other divergent regions containing 20 or more genes with fluorescence intensity ratios of <0.3 that were determined to be absent in strain NCTC 11168. These divergent regions, also shown in Fig. 1, corresponded to the four CJIEs, described previously (4), which are present in strain RM1221.
![]() View larger version (27K): [in a new window] |
FIG. 1. Genomic index of RM1221 genes in strain NCTC 11168. C. jejuni gene status was determined by DNA microarray analysis. The gene status is color coded as follows: blue, present; yellow, variable; red, absent; gray, no data. For cutoffs of absence and presence predictions, refer to Materials and Methods. The genes are represented on the circular RM1221 chromosome. Prominent hypervariable regions are indicated (Table 2), including the four CJIEs.
|
|
View this table: [in a new window] |
TABLE 2. Intraspecies hypervariable regions in 35 C. jejuni strains
|
![]() View larger version (45K): [in a new window] |
FIG. 2. Detailed genomic index of NCTC 11168 (HS:2) and RM1221 (HS:53) capsular regions in 35 strains of C. jejuni. The numbers at the bottom of the figure correspond to the strains in Table 1. The gene status is color coded as follows: blue, present; yellow, variable; red, absent; gray, no data. For cutoffs of absence and presence predictions, refer to Materials and Methods.
|
![]() View larger version (55K): [in a new window] |
FIG. 3. Detailed genomic index of the four CJIEs of RM1221 in 35 strains of Campylobacter jejuni. The numbers at the bottom of the figure correspond to the strains in Table 1. The gene status is color coded as follows: blue, present; yellow, variable; red, absent; gray, no data. For cutoffs of absence and presence predictions, refer to Materials and Methods. An average linkage hierarchical clustering of the C. jejuni strains was compiled in GeneSpring version 7.3 from the CGH data for each element with the standard coefficient correlation and bootstrapping. (A) CJIE1; (B) CJIE2; (C) CJIE3; (D) CJIE4. Two sets of genes, CJE0590 to CJE0598 (CJIE2) and CJE1444 to CJE1451 (CJIE4), show a greater than 95% sequence identity between these elements. Due to this identity, we observed that the CJIE2-PCR positive strains, RM1048, RM1852, and RM2240, cross hybridized with the gene region of CJIE4. Similarly, CJIE4 PCR-positive strains, RM1046, RM1170, RM1413, RM1449, RM2229, and RM2232, cross hybridized with the gene region of CJIE2 (data not shown). We also observed this cross hybridization for CJE0550 and CJE1421, genes encoding putative DNA methyltransferases. Thus, all of these genes were removed from the cluster analysis and are not represented in the figure.
|
The plasmid-like CJIE3 was represented by the most PCR-positive strains in our sample set, yet this element appears to be the most diverse region, forming four clusters, (Fig. 3C). As with CJIE1, CJIE2, and CJIE4, the divergent CJIE3 elements in the other strains appeared to have a mosaic arrangement of genes. Particular CJIE3 gene modules were found in both PCR-positive and PCR-negative strains. All 14 PCR-positive strains possessed at least one divergent gene in the CJIE3 element compared to strain RM1221, with two strains, RM2240 and RM1852, forming a subcluster with RM1221 (distance score = 0.053). Four PCR-positive strains, RM1163, RM1508, RM1413, and RM1503, possessed less than 22% (7/33) of the CJIE3 genes present on the array; these CJIE3 genes were mostly contiguous. By CGH, we were able to identify six PCR-negative strains (RM1170, RM1244, RM1409, RM1246, RM2769, and RM1248), each harboring three or more genes from CJIE3. Of particular note was isolate RM1248, which subclustered with seven CJIE3 PCR-positive strains (RM1047, RM1050, RM1188, RM1449, RM1850, RM2229, and RM2232) (distance score = 0.512) and shared 3 gene modules. The CGH analysis for this strain indicated that it possessed gene CJE1094, the gene used in the CJIE3 PCR. The negative result for this strain in our PCR-based assay was probably due to screen sequence differences at the primer binding site(s). Recently, strain RM1047 (TGH 9011) was shown to possess genes that encoded products similar to those encoded by the pathogenicity island of H. hepaticus, genes that are absent from strain NCTC 11168 (21). In fact, several of the TGH 9011 genes have DNA sequence homology to the CJIE3 genes of strain RM1221 (Table 3) and were found to be positive based on the CGH analysis. The other five PCR-negative strains (RM1170, RM1244, RM1246, RM1409, and RM2769) possessed at least two genes with sequence similarity to Campylobacter plasmids (2, 17) (Table 3). Further analysis would be required to determine if CJIE3 genes in these other strains are chromosomal, as they are in C. jejuni strain RM1221, or plasmid borne.
|
View this table: [in a new window] |
TABLE 3. CJIE3 genes with similarity to Campylobacter plasmids or C. jejuni strain TGH 9011 contigs
|
The other three elements, CJIE2, CJIE3, and CJIE4, were discovered to be adjacent to tRNA genes. This is commonly observed in other bacterial pathogens, including the Salmonella pathogenicity island 1 and the locus of enterocyte effacement in Escherichia coli (9). The chromosomal location of CJIE3 was determined for several C. jejuni strains that were positive for this element. We designed a PCR primer set, CJ1230-2/CJ1299-2, to the regions flanking CJIE3 in strain RM1221. Strains NCTC 11168 and 81116 both lack CJIE3; therefore, a PCR amplicon was produced with this primer set. As expected, no fragment was amplified by PCR when strain RM1221 was examined due to the presence of the approximately 40-kb CJIE3. Of 23 C. jejuni strains examined for the chromosomal location of CJIE3, only 7 yielded an amplicon by PCR, indicating that CJIE3 in each of these 7 strains was inserted at a site different from strain RM1221 or that the CJIE3 genes in these strains were plasmid borne. In contrast, 16 C. jejuni strains possessed the element at the same position as strain RM1221, since the DNA of these strains did not amplify with the primer set, CJ1230-2/CJ1299-2. The fact that many of the CJIE3s in the tested strains are in the same chromosomal location as strain RM1221 may be due to the nature of the insertion. If CJIE3 requires insertion into a tRNA site, as in RM1221, then only a limited number of potential integration sites will be available. Determination of the locations of CJIE2 and CJIE4 in different C. jejuni strains will require additional analysis.
Concluding remarks. While several CGH studies have demonstrated regions of high variability between strains of C. jejuni and strain NCTC 11168 (3, 11, 12, 20, 23), the CGH results reported here, using C. jejuni microarrays containing genes from both strain NCTC 11168 and strain RM1221, enabled us to examine additional levels of diversity in C. jejuni. Our analysis demonstrated the occurrence of four CJIEs, RM1221-like genomic elements, in other strains of C. jejuni. Although a PCR screen also permitted the identification of C. jejuni strains possessing CJIEs, results obtained from CGH analysis by DNA microarrays revealed gene diversity between PCR-positive strains within the corresponding CJIEs and also identified several CJIE genes in PCR-negative strains. The presence of bacteriophage-related genes within the CJIEs suggests that various bacteriophage that lysogenize C. jejuni contribute to an additional genomic diversity in C. jejuni. The modular pattern of CJIE genes in other strains of C. jejuni supports the model for bacteriophage genome dynamics, as originally proposed by Hendrix et al. (7), in which a large common bacteriophage gene pool results in the creation of mosaic bacteriophage genomes through horizontal exchange. Indeed, the contribution of bacteriophages and their mosaic structure to the diversification of the E. coli and Salmonella enterica chromosomes has also been observed by CGH (22, 24). In the present study, we also demonstrated that in all of the C. jejuni strains examined, CJIE1, the Mu-like element, was located in different chromosomal locations, possibly due to random insertion events, providing another mechanism to increase the genomic diversity of C. jejuni.
Additionally, the global genomic nature of our CGH analysis allowed us to corroborate the clustered nature of highly divergent genes, as observed in previous studies (3, 11, 12, 20, 23). Most highly divergent genes that were identified in our study occurred either within the four CJIEs or within the 17 intraspecies hypervariable regions. Furthermore, with distinct genes from the hypervariable regions of both strains NCTC 11168 and RM1221, the CGH results demonstrated that many of the C. jejuni strains in this study were divergent not only from strain NCTC 11168 but also from RM1221. Future work will further explore if these genomic regions could be developed as informative epidemiological tools in distinguishing various strains of C. jejuni.
This work was supported by the United States Department of Agriculture, Agricultural Research Service CRIS project 5325-42000-041.
Published ahead of print on 30 August 2006. ![]()
|
|
|---|
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»