| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
,
Division of Human Biology, Fred Hutchinson Cancer Research Center, Seattle, Washington,1 Program in Molecular and Cellular Biology, University of Washington, Seattle, Washington2
Received 26 July 2006/ Returned for modification 11 October 2006/ Accepted 2 November 2006
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
While mechanisms of virulence with respect to pathogen-host interactions are poorly understood, several important mechanisms of survival required to initiate the infection have been characterized (42). H. pylori absolutely requires a family of genes involved in the production of urease, an enzyme critical for neutralizing the pH around the organism during exposure to the acidic lumen of the stomach. H. pylori also requires a large set of motility genes related to the successful production and operation of its flagellum. Without motility, H. pylori cannot penetrate the mucous layer which protects the stomach epithelium from the acid it produces. While the bacteria can persist in the mucous layer, they also attach tightly to gastric epithelial cells via a number of adhesins. Particularly important are two outer membrane proteins, BabA and SabA, that bind oligosaccharide antigens present on cellular subpopulations of the gastric epithelium (43). Successful colonization thus requires intimate interaction with the epithelium itself.
Several virulence factors are thought to be important once contact with the host cell epithelium is established. VacA, a secreted protein, has a vast array of activities linked to it, including membrane insertion, anion-conducting channel activity, alteration of transepithelial resistance, inhibition of antigen processing, and induction of apoptosis (50). Presumably, the loss of membrane integrity in the host cell increases nutrient availability at the host cell surface. The cag pathogenicity island (PAI) is composed of 27 genes. The cag PAI contains homologues to type IV secretion system (T4SS) components in other gram-negative organisms such as Agrobacterium tumefaciens and Legionella pneumophila (plant and human pathogens, respectively). The presence of the cag PAI correlates with more-serious disease outcomes, implying that when functioning, it is important in pathogenesis. The CagA protein (also encoded in the island) is translocated to the host cell cytoplasm by the T4SS (7). Once CagA is translocated into epithelial cells, it remains associated with the host membrane and becomes tyrosine phosphorylated on C-terminal repeat motifs, Glu-Pro-Ile-Tyr-Ala (EPIYA motifs), by proteins of the Src family of tyrosine kinases. Both phosphorylated and unphosphorylated CagA proteins have been reported to induce cell signaling pathways resulting in altered spreading, migration, and adhesion of epithelial cells. The Ras/mitogen-activated kinase kinase/extracellular signal-regulated kinase (MEK/ERK) and Src homology 2 (SH2) domain containing tyrosine phosphatase (SHP-2) pathways are some of the pathways reported to be activated by cagA-positive strains, explaining observations of increased cell proliferation during H. pylori infection, a hallmark for a precancerous state.
Several screens for colonization or virulence factors of H. pylori have been described using libraries of insertional mutants made either by the transposition of H. pylori DNA clones propagated in Escherichia coli (shuttle mutagenesis) (26, 31, 44, 45) or by the integration of plasmids containing random small pieces of H. pylori DNA (5, 6). Both these techniques take advantage of H. pylori bacteria's ability to integrate homologous DNA by recombination. So far, genetic screens have focused mostly on in vitro phenotypes thought to be important for infection of the human stomach, such as motility, acid survival, urease activity, adherence to gastric epithelial cells, and the ability to take up exogenous DNA (reviewed in reference 23). Recently, the first in vivo screen using a gerbil model of infection was reported (31). That study queried 960 mutants, giving 252 candidate mutants corresponding to 47 genes (due to redundancy of the original 960-mutant pool). Both the previously characterized colonization genes (such as those that encode motility factors) and the novel components (such as collagenase) were identified. None of the screens to date has been saturating, and none has produced either novel T4SS-dependent genes or genes required for the regulation of virulence gene expression.
We recently developed a genome-saturating library of transposon mutants and a method to track pools of these transposon mutants (MATT), using a whole-genome microarray (54). We have initiated a screen for H. pylori colonization factors by using these tools in a mouse infection model (32). Here we present the results of our initial screen in two different mouse-adapted strain backgrounds. Our goal was to identify genes contributing to bacterial colonization and/or persistence in the stomach.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Molecular biology. DNA manipulation (restriction digests, PCR, agarose gel electrophoresis, and Southern blotting) was performed according to standard procedures (3). Genomic DNA was prepared from H. pylori by using a Wizard genomic DNA preparation kit (Promega). A list of primers used for PCR and sequencing is given in Table S1 in the supplemental material.
Mouse-adapted strains. The mouse-adapted strains used in this study were NSH79 and NSH57. NSH79 is a mouse-adapted H. pylori strain obtained after multiple passages through mice of a mixture of strains, including strain SS1 described by Lee and coworkers (32). We distinguish the NSH79 strain from strain SS1 because microarray-based comparative genomic DNA hybridization revealed that the NSH79 clone has numerous genetic changes (gene acquisition and loss) relative to the SS1 strain (N. R. Salama, L. J. Thompson, A. Lee, and S. Falkow, unpublished observations) (Fig. 1). NSH57 is a mouse-adapted derivative of strain G27 (15). NSH57 was generated by inoculating approximately 5 x 108 G27 bacteria, as described below, into four FVB mice. At 3 weeks, mice were sacrificed, and the output from one mouse was pooled (passage 1), amplified by growth in liquid culture, and used to inoculate four new FVB mice (passage 2). For the third passage, 30 mice were infected. These mice were sacrificed at 3 weeks, and 26 of the 30 inoculated mice were found to have detectable infection, with bacterial loads ranging from 3.7 to 5.3 log CFU/g. Bacteria from three mice giving bacterial loads of >5.0 log CFU/g were separately inoculated into C57BL/6 mice. At 3 weeks, the mice were sacrificed. All mice were infected with bacterial loads ranging from 2.6 to 5.4 log CFU/g. Single-colony isolates were obtained from mice with >5.0 log CFU/g, and one isolate, designated NSH57, was used for further study.
|
Mouse infections. Female C57BL/6 mice, 24 to 28 days old, were obtained from Charles River Laboratories and certified free of endogenous Helicobacter sp. infection. For infection of pools of mutants, frozen stocks of the pools of 48 clones were plated on solid medium to give at least 1,000 individual colonies. After 2 (NSH57) or 3 (NSH79) days of growth, the bacteria were harvested directly from the plates and resuspended in BB10 medium to give approximately 5 x 109 bacteria per ml. Cultures were examined by microscopy to ensure a spiral shape and motility indicative of mid- to late logarithmic growth phase. Four mice were inoculated by oral gavage with 0.5 to 1 ml. After dosing, an aliquot of the inoculum was plated on HB-CHL plates to generate the input pool. The mice were housed in sterilized microisolator cages with irradiated PMI 5053 rodent chow, autoclaved corncob bedding, and acidified, reverse-osmosis purified water provided ad libitum. All studies were done under practices and procedures of Animal Biosafety Level 2. The facility is fully accredited by the Association for Assessment and Accreditation of Laboratory Animal Care, International, and all activities were approved by the Institutional Animal Care and Use Committee. Two mice at 1 week and two mice at 1 month were euthanized by inhalation of CO2. The glandular stomach was removed and cut along the greater curvature; any remaining food was removed, and tissue samples were placed in 2 ml of BB10 and homogenized with a Powergen 125 homogenizer (Fisher Scientific). The homogenate was cleared by centrifugation at 500 rpm in a clinical centrifuge (model 5810R; Eppendorf), and 200 µl was plated to generate the output pool.
For competition experiments, each indicated null mutant strain and the parental wild-type strain were grown from frozen stock in liquid culture to mid- to late logarithmic growth phase. The wild-type and mutant bacteria were combined to give approximately 2.5 x 108 bacteria of each in 5-ml and 1-ml amounts and inoculated by oral gavage in each of four mice. After inoculation, a portion of the inoculum was plated on HB-BAC plates and HB-CHL plates to enumerate the actual number of wild-type and mutant bacteria present in the inocula, respectively. After 1 week, the mice were euthanized and the stomach removed as described above. The stomach was then cut in half longitudinally, and one half of the stomach was placed in 0.5 ml BB10 for homogenization with disposable pellet pestles (Fisher Scientific). Dilutions of the homogenate were plated on HB-BAC to enumerate total bacteria and on HB-CHL plates to enumerate mutant bacteria.
For the 50% infective dose (ID50) experiments, the indicated strains were grown in liquid culture to mid- to late logarithmic phase and concentrated to give approximately 5 x 108 bacteria per ml. Serial 10-fold dilutions were prepared and inoculated into each of five mice. After 7 to 10 days, mice were euthanized and processed as described above for competition experiments with plating on HB-BAC plates.
Microarray-based comparative genomic hybridization (CGH array). Microarray design and hybridization conditions have been described previously (28, 52). Each isolate was analyzed on at least two microarrays, generating four potential data points for each gene. Data points were excluded due to low signal or slide abnormalities, and only those genes with two measurements were analyzed. Data were normalized using the default-computed normalization of the Stanford Microarray Database (4), and the mean of the log 2(red/green ratio) was computed. The cutoff for the absence of a gene was defined as a log 2(red/green ratio) of 1.0 based on test hybridizations (52). Data were simplified into a binary score, analyzed with CLUSTER software (http://rana.lbl.gov/EisenSoftware.htm), and displayed with TREEVIEW (19). The complete data set is available in Table S3 in the supplemental material, and raw data are available at http://genome-www5.stanford.edu.
Microarray analysis of mutant pools.
An adaptation of MATT (54) was used to estimate the location of all 48 mutants in each pool (outlined in Fig. 2). Briefly, genomic DNA was prepared from the input and output pools of each mouse. Sequences flanking one side of the transposons in each pool were amplified in a two-step PCR and labeled for microarray hybridization as previously described (54). The transposon primer N2 was redesigned and named N3 (5'-CTTTAATACGACGGGCAATTTGCACTTCAG-3'), based on hybridization optimization of known mutant pools (data not shown). Additionally, only primer CEKG2C was used as the random anchored primer for amplification. For all experiments (each mouse represents a single experiment, with an input pool and an output pool), the input pool DNA was labeled with Cy-3 (green [G]) and the output pool was labeled with Cy-5 (red [R]). Therefore, the red/green (R/G) ratio for each gene provides a measure of fitness for any given gene mutated in the input pool: a yellow spot (R/G, 1) present in both the input and the output represents no phenotype; a green spot (R/G, <1) present in the input but absent in the output represents a candidate mutant. For the set of four infections, the input pool DNAs were identical. Microarray scanning and analysis were performed on a GenePix 4000B scanner (Axon) using GenePix Pro 5.0 software (Axon). Spots were filtered for slide abnormalities (flag, 0) and uniformity of spot intensity (regression correlation of the spot,
0.6). Because each pool contained only 48 insertions, most of the 1,660 gene spots on the array were expected to give a background signal. Previous experiments using MATT to map transposon insertions had revealed that not all insertions are detected and that in some cases multiple genes adjacent to the transposon insertion get amplified (54). Thus, we expected roughly 100 spots in each experiment array to give a strong signal, but the actual number would vary from experiment to experiment. To select spots for analysis, we first computed the average background-subtracted signal intensity in the Cy-3 channel (input pool) of all the spots except the brightest 100 spots. We then chose for further analysis spots whose intensity was greater than 4 standard deviations from this mean. The average number of spots selected was 150 (range, 89 to 164). The log2 R/G was collected for each of these gene spots. Two data sets were collected for each experiment (mouse infection): one where the amplification of flanking DNA was performed from the left side of the transposon (S primers) and a second where the amplification of flanking DNA was performed from the right arm of the transposon (N primers). The data from all arrays were then collated using a relational database (Access software; Microsoft). To define an insertion site, we required that its gene spot be labeled from both sides of the transposon insertion, and thus there were data for the same gene spot in both data sets. These genes were assigned a probability of insertion of 1. To account for transposon insertions where the insertion lay near the end of a gene, we also defined insertions where, after the gene spots were arranged in chromosomal order, there was signal for a gene spot in the reaction from one side of the transposon and signal for the nearest neighboring gene spot in the reaction from the other side of the transposon. In this case, we gave each of the two neighboring genes a probability of insertion of 0.5. We then queried the color of the insertion marking the gene spots. If the spot was yellow [present in both the input and the output, log2(R/G)
2], we assigned a positive value; and if it was green [present in the input, absent in the output, log2(R/G)
2], we assigned a negative value. Thus, insertions with no apparent phenotype (present in both the input and the output) received a value of 1 or 0.5, and insertions with a mutant phenotype (absent in the output) received a value of 1 or 0.5. The data for all the pools tested in each strain background were analyzed separately. For each strain background, the data were summed across all the pools, and genes with a value of less than 1.5 (a value of 0.5 in each of three mice) were considered candidate mutants.
|
Generation of isogenic knockouts. Genes whose transposon locations were verified by PCR and sequencing were considered for further analysis. Null alleles were constructed using a vector-free allelic replacement strategy, as described previously, to generate alleles where a chloramphenicol antibiotic resistance cassette replaced most of the coding sequencing of the gene while preserving the start and stop codons (12, 54). The resistance cassette contains its own promoter but lacks a transcriptional terminator and, in all cases, was cloned in the same direction of transcription as the native gene. Briefly, N-terminal and C-terminal fragments of each gene were amplified using the upstream and downstream primers combined with gene internal primers such that both fragments were 250 to 400 bp (see Table S1 in the supplemental material). The gene internal primers had sequences complementary to the full lengths of primers C1 and C2, used to amplify the Campylobacter coli Cat gene, appended to the 5' end. PCR products from each of three individual PCRs (the N terminus, the C terminus, and the Cat gene) were purified using DNA Clean & Concentrator-5 (Zymo Research) or by agarose gel electrophoresis followed by QIAEX II gel extraction (QIAGEN) if there were multiple bands. A final PCR was performed with approximately 100 ng of each of the three PCR products as template and the upstream and downstream primers to generate the knockout cassette. This final PCR product was verified by agarose gel electrophoresis, and 10 µl was directly used for the natural transformation (66) of each strain background (NSH79 and NSH57) and selected on HB-CHL media. Four to eight clones were evaluated by PCR to confirm replacement of the wild-type allele with the null allele, urease, and motility phenotypes. A single clone was used for infection experiments.
AGS coculture experiments. The human gastric adenocarcinoma cell line AGS was maintained in 10% CO2 in Dulbecco's minimal essential medium supplemented with 10% fetal bovine serum (FBS). Cells were seeded at a density of 1 x 105 in 24-well plates. H. pylori strains were grown at 37°C overnight in BB10 to an optical density at 600 nm (OD600) of 0.4 to 0.8, harvested, and resuspended in DB medium (81% Dulbecco's minimum essential medium, 9% Brucella broth, 10% FBS) at 1 x 106 to 2 x 106 bacteria/ml, and 1 ml was used to inoculate each well. At each time point, the supernatant was harvested, centrifuged, and frozen for interleukin 8 (IL-8) analysis (Biotrak enzyme-linked immunosorbent assay system; Amersham Biosciences). To detect phosphorylated and total CagA, the remaining cells and bacteria in each well were washed three times with 1 mM Na3VO4 in phosphate-buffered saline (Invitrogen) and lysed with 100 µl 2x sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) sample buffer (0.25 M Tris-Cl [pH 6.8], 4% glycerol, 4% SDS, 0.001% bromophenol blue, 2% 2-mercaptoethanol). Samples were resolved by SDS-PAGE and transferred to polyvinylidene difluoride membranes. To visualize phosphorylated CagA, membranes were probed with anti-phosphotyrosine-PY20 (BD Transduction Laboratories) or 4G10 (Upstate), followed by 1:10,000 goat anti-mouse horseradish peroxidase (HRP) (Amersham Biosciences). Immunoreactive proteins were visualized using ECL Plus Western blotting detection reagent (Amersham Biosciences). To measure CagA, the same membranes were probed with 1:10,000 anti-CagA polyclonal antibody pAs (2), followed by 1:10,000 goat anti-rabbit HRP (Amersham Biosciences), and visualized using ECL Plus.
Motility assays. Individual clones and the wild-type parent strain were inoculated into plates containing 28 g/liter Brucella broth, 5% FBS, and 0.4% agar. After 4 days of growth, the plates were evaluated for halo production.
Urease assays. Bacteria were resuspended from freshly growing plates (24 to 48 h) in 100 µl urease broth (Difco) in 96-well plates in duplicate. The wells were covered with tape and incubated at 37°C in 14% CO2 for 24 h. Wells where the pH-sensitive indicator dye changed from yellow to pink were scored urease positive.
Growth rate determination. Null mutant or isogenic wild-type bacteria were inoculated from 2 to 3 days of growth on plates into 2 ml BB10 medium for overnight cultures. The resulting cultures were checked for spiral morphology and motility and then diluted to an OD600 of 0.1. Aliquots were removed every 3 to 6 h to determine the OD600. The resulting growth curves were compared graphically using Excel (Microsoft).
Microarray data accession numbers. Raw data for each array were submitted to the Stanford Microarray Database for normalization and warehousing (4). The raw data are available at http://smd.stanford.edu/cgi-bin/publication/viewPublication.pl?pub_no=564.
| RESULTS |
|---|
|
|
|---|
As described in Materials and Methods, NSH57 and NSH79 were both derived from multiple mouse and laboratory passages of previously characterized H. pylori isolates (G27 [15] and SS1 [32], respectively). We reported previously that the adaptation of human clinical isolates of H. pylori to animal hosts can result in gene deletion, particularly in the cag PAI (27, 30, 60). Therefore, we performed comparative genomic DNA hybridizations to the whole-genome cDNA microarrays to detect whole-gene-level alterations relative to those of the starting strains (see Table S3 in the supplemental material). No changes in genomic content were detected with NSH57 relative to the genome of the parental G27. Strains NSH79 and SS1, however, differed in the presence and absence of 114 genes, indicating that they were quite genetically distinct and that NSH79 likely is unrelated to SS1. Hierarchical clustering of the strains based on the presence and absence of genes revealed that NSH79 was as equally distinct from SS1 and its pre-mouse clinical isolate (60) as it was from G27 and NSH57 (Fig. 1A). We did not find any major alterations in the cag PAI (see Table S3 in the supplemental material) among any of these strains.
Even in the absence of gene deletion, several studies have shown that many mouse-adapted strains, including SS1, display attenuation of the proinflammatory activity of the cag PAI T4SS (17, 49, 63). To check the activity of the T4SS in our two mouse-adapted strains, we assayed both IL-8 release and host cell translocation-dependent tyrosine phosphorylation of CagA after coculture of bacteria with AGS gastric epithelial tissue culture cells (Fig. 1). Both G27 and the mouse-adapted variant of this strain, NSH57, showed robust induction of IL-8 and translocation of CagA into host cells (Fig. 1B and D and see Fig. S1 in the supplemental material). In contrast, strain NSH79 showed less induction of IL-8, and we were unable to detect phosphorylated CagA in spite of robust CagA expression, even at 24 h (Fig. 1B, C, and D). Thus, NSH79 appears to have a partially attenuated cag PAI T4SS and/or a phosphorylation-defective cagA, while the NSH57 cag PAI T4SS is fully functional.
Determination of pool sizes. In order to screen as many mutants as possible using the lowest number of mice, we wanted to infect mice with pools of mutants and monitor the fate of these mutant clones using MATT. The number of mutants that can be screened simultaneously depends on the number of clones that can independently establish infection in the stomach mucosa at one time (41). To address this issue, we performed pilot infections with 24, 48, and 96 random clones from the NSH79 mutant library for 2 weeks. To do this, each pool of transposon mutants was used to infect three mice. Genomic DNA was prepared from input bacteria and output bacteria from each of the three mice. This DNA was labeled for MATT analysis either from one mouse individually or by combining the bacterial outputs from all three infected mice. As summarized in Table 1, MATT was able to detect an increasing number of clones as the input number of clones increased from 24 to 48. However, less than the expected number of clones was detected when a mixture of 96 independent clones was used for the infection. Additionally, there were six clones that disappeared from the output pool when one mouse was analyzed individually, but these clones were present when the outputs from the three mice were pooled. This indicates that somewhere between 48 and 96 clones, the criteria of independent action of clones fail to be satisfied. Thus, we have shown that MATT can be used to track clones during mouse infection and that 48 clones can be screened simultaneously in a single mouse.
|
Each mouse was analyzed separately, using a modification of MATT (54) as outlined in Fig. 2. This method uses semirandom PCR to amplify the DNA next to transposon insertions. Pooled genomic DNA from the input and output bacteria from each mouse was amplified and labeled with different fluorescent dyes separately and then cohybridized on a single microarray. Transposon insertions were assigned if detected by the same gene spot in separate hybridization experiments initiated from the left and right sides of the transposon. Genes were assigned a probability of insertion value of 1 (1 for a mutant phenotype, see Materials and Methods) if the insertion mapped to the middle of the gene. If the insertion lay at the end or between genes, each of the two neighboring genes was assigned a value of 0.5 (or 0.5). For each strain background, data were summed across all the pools and genes with a value less than 1.5 (a value of 0.5 in three mice) were given further consideration. The summed data for each strain background are given in Table S4 in the supplemental material.
Using MATT, we were able to track the behavior of 1,378 of the 2,400 initial transposon insertions which mapped to 758 genes. By requiring observation in the microarray hybridizations from both sides of the transposon, we excluded data for many insertions. We felt that this loss of sensitivity was balanced by the elimination of false-positive results. Some genes had insertions observed in multiple pools. Of the 758 genes measured in our screen, only 250 were measured in both strains. This confirms that neither screen was saturating. As expected, most transposon mutant clones were present in both the input and the output pools, but 223 genes (29%) showed a colonization phenotype (see Table S4 in the supplemental material). Of the candidate mutant genes, 98 were queried in both strain backgrounds. For 23 genes, the data from both strain backgrounds indicated a phenotype in colonization after transposon insertion-mediated gene inactivation. However, 75 clones were eliminated during infection in one strain background but not in the other. These discrepancies could arise from the locations of transposon insertion being different for the two screens, where in one case it inactivates a gene while in another it does not. Alternatively, there may be functionally redundant genes present in one strain but not in the other.
Candidate colonization genes fall into expected and unexpected function categories. The genes we identified as candidate colonization factors in our screen fell into a variety of function categories, some of which were expected based on previous work with H. pylori or other pathogenic organisms, including motility, lipopolysaccharide or exopolysaccharide biosynthesis and modification, urease production, iron homeostasis, and stress response (Table 2). The largest number of genes (81 genes) does not have informative homologies based on DNA sequence alignment and is annotated as hypothetical. Only 21 of these genes are conserved in other sequenced organisms or have recognizable domains (e.g., GTP binding). The second most frequent category (19) was motility and chemotaxis, including a regulator of flagellar biosynthetic gene expression (HP0703, flgR), genes involved in flagellar biosynthesis (12), three methyl-accepting chemotaxis proteins, and three genes involved in transducing signals to the flagellar motor. A number of genes (14) can be grouped together because they are involved in biosynthesis, including three genes involved in nucleotide metabolism, six genes involved in amino acid biosynthesis, three genes involved in biosynthesis of cofactors, and two genes involved in lipid or fatty acid metabolism. Similarly, we identified 11 genes annotated as transporters of various nutrients and a putative efflux pump. Identification of these last two function classes is consistent with a host environment that is limited in the variety of nutrients essential to H. pylori growth compared to the rich media used to propagate the mutant pools in vitro.
|
We identified six genes from the cag PAI in our screen. In general, these genes showed discordant data in the two strain backgrounds, and as described below, we confirmed this strain-specific phenotype for some of these genes. Two other genes have annotations that implicate them in virulence: a paralogue of the vacA cytotoxin (HP0609) and an orphan homologue of the Agrobacterium tumafaciens type IV secretion apparatus subunit virB4 (jhp0917). We identified eight genes that fell into the categories of protein synthesis and modification and included three genes involved in lipoprotein biogenesis, a putative sialoglycoprotease, and three genes predicted to affect protein translation. Seven genes are annotated to be involved in DNA metabolism, including DNA repair proteins and putative recombination enzymes. Seven of the transposon insertions showing an attenuated phenotype mapped to one of several transposase genes of the IS605 insertion sequences which are present in multiple copies in the genome. The mechanism by which these insertions affect colonization is likely indirect. Six genes are annotated to be involved in iron homeostasis, including the ferric uptake regulation protein (fur, HP1027), the nonheme iron-containing ferritin (pfr, HP0653), and four genes with homology to outer membrane iron transporters. Four genes are annotated to be involved in stress response, including a superoxide dismutase (HP0389), a chaperone (HP0210), a peroxyredoxin (HP0136), and a periplasmic acid phosphatase (HP1285). In addition to the flgR and fur genes, we identified two additional candidate transcriptional regulators in our screen: the two-component response regulator HP1043 and tenA (HP1287), a gene shown to induce expression of extracellular enzymes in Bacillus subtilis (47). Finally, we detected a single gene that was annotated to be involved in cell division (fic, HP1159).
Validation of two pools. We initially focused on a complete analysis of the candidate mutants from two mutant pools from the screen of the NSH79 strain background (NSH79-2 and NSH79-3). We detected transposon insertions in 44 different genes from these two pools, three of which were detected in both pools. Of these 44, 12 genes were suspected to be involved in colonization because they were not detected in the output pool (Table 3). Consistent with the results from the entire screen, the majority (seven genes) had no annotation, based on nucleotide sequence homology, and one is annotated as a conserved hypothetical integral membrane protein. One gene was annotated as a nonfunctional type II restriction endonuclease, presumably because it lies adjacent to a putative type II N6 adenine-specific methyl transferase (HP0369), though the sequence of HP0368 has degenerated to the point that no conserved endonuclease domain can be identified.
|
When the data from these pools were compared to the larger data set, some of these genes appeared in multiple pools. Additionally, some were detected in both screens. One gene, HP0503, showed a colonization phenotype in both strain backgrounds, while HP0022, HP0368, HP0529 (cag9), and HP1028 were not candidate mutants in the NSH57 strain background. This discordance could result from strain-specific differences in genetic redundancy or phenotypic buffering or from differences in the site of transposon insertion in the different clones. Transposon insertions at the very 5' or the 3' 10 to 15% portion of the gene have been suggested not to severely affect gene function (29).
To evaluate the robustness of the array data, we first identified the exact location of each transposon in the region predicted by the arrays. Gene-specific primers were designed for the 5' and 3' ends of each gene and were used in combination with the transposon-specific primers to amplify genomic fragments flanking the insertions. These PCR products were then sequenced. The position of the transposon insertion for each candidate mutant is given in Table 3. We made independent null alleles of each of these genes by replacing most of the coding sequence with a nonpolar antibiotic resistance cassette. Null alleles were made in both the NSH79 and NSH57 mouse-adapted strain backgrounds to determine if the phenotype was strain specific or more universal. We tested these mutants in a 1:1 competition experiment with wild-type bacteria. After 1 week of infection, we harvested the stomach and determined the competitive index (CFU mutant bacteria/CFU wild-type bacteria). This competitive index number was corrected by dividing the actual input ratio enumerated from plating the inocula after infection. As seen in Fig. 3, all the gene loci tested gave a phenotype in the NSH79 strain background, which was the strain background of the pools being analyzed. Many of these gene loci also gave a phenotype in the second strain background (NSH57). The one gene locus, HP0503, which was predicted to give a phenotype in both strains did so, as did HP0137, for which there were no data from the NSH57 screen. When the four genes that gave discordant data were examined, HP0022 and HP0529 showed a strain-specific phenotype as expected. HP0368 and HP1028, however, gave a phenotype in both strain backgrounds. We examined the location of the transposon for the clones detected in the NSH57 screen and found an insertion in gene HP0368 at nucleotide 40 (out of 399) of the coding sequence, and the insertion in gene HP1028 was at nucleotide 150 (out of 450). For the NSH57 screen, the insertion in HP1028 appeared attenuated in some but not in all mice and thus was not considered a candidate mutant. In addition to HP0022 and HP0529, four additional genes showed a strain-specific phenotype.
|
|
Two of the genes we retested are part of the cag PAI. HP0529 encodes a homologue of the Agrobacterium tumefaciens virB8 gene, which is a core component of the T4SS, and HP0538 encodes CagN, a protein that localizes to the inner membrane but has no known effects on the induction of IL-8 or CagA translocation into host cells (8, 20). Mutation of these genes displayed a phenotype only in the NSH79 strain background. While the NSH79 cag PAI T4SS appears partially attenuated (Fig. 1B), induction of host cell IL-8 secretion was further diminished by the mutation of the virB8 homologue HP0529 and not by the mutation of cagN (HP0538). This indicates that NSH79 induced IL-8 secretion by virtue of its partially functional cag PAI T4SS and not other bacterial components (Fig. 1D). We were unable to assess the contributions of HP0529 and HP0538 to CagA translocation in NSH79 because the wild-type parent showed no detectable CagA phosphorylation in our assay system, though it remains possible that some CagA does become translocated into host cells during mouse infection. These results suggest that a partially functional cag PAI T4SS provides a competitive advantage in the NSH79 strain background but not in the NSH57 strain background.
| DISCUSSION |
|---|
|
|
|---|
While screening of very-high-complexity libraries in single animals has been possible for systemic infection models such as interperitoneal mouse inoculation with Salmonella enterica and Mycobacterium tuberculosis (13, 55), colonization of mucosal surfaces often results in survival bottlenecks such that clones no longer show independent action (37, 40). These bottlenecks presumably result from innate defenses of mucosal barriers. In our study, we empirically determined an optimal pool size of 48 clones. This allowed us to screen 1,200 transposon clones from two libraries in different strain backgrounds for a total of 2,400 clones. We used whole-genome microarrays to monitor the behavior of individual clones in the pools and to identify those that could not survive. Our method of detecting the transposon insertion site on the microarray relied on a semirandom PCR amplification step of the DNA sequences flanking the transposon, using an anchored random primer expected to anneal approximately once in every 1,000 bases in the H. pylori genome (54). In practice, not all transposon-flanking sequences are amplified efficiently. We further required a microarray signal in the amplifications from both sides of the transposon to eliminate technical artifacts. This resulted in the detection of only 1,378 clones out of 2,400, or 57%. These insertions mapped to 758 genes, of which 223 had a predicted colonization defect. A recent signature-tagged mutagenesis-based screen of H. pylori stomach colonization using a gerbil infection model queried 960 transposon mutants and identified 252 candidate mutants corresponding to 47 genes (31). These two in vivo screens yielded similar percentages of candidate mutants (29% in our screen versus 26% in the previous screen). This suggests that a high percentage of the H. pylori genome is devoted to the ability to colonize and persist in the stomach.
Comparison to other studies. When comparing our data to the above-mentioned screen for colonization of the gerbil stomach, we collected data for 22 of the 47 in vivo essential genes described previously (31). We detected a colonization defect for 14 of these genes, including four urease genes, five genes involved in motility or chemotaxis, alpha-ketoglutarate permease (HP1091), UDP-glucose 4-epimerase (HP0360), a gene with homology to an exopolysaccharide biosynthesis protein (HP0366), and two hypothetical protein-encoding genes (HP0486 and HP1525). None of these fourteen mutants were retested in the previous work, so our data represent the first independent analysis of many of these genes' roles in colonization. The remaining eight genes for which there were data in our screen did not make our cutoff point as candidate mutants, either because they were not measured in all the mice of a particular pool or because there were conflicting data from different pools. Since we do not know the precise site of the transposon insertion, it may be that the insertions we measured were not inactivating. Alternatively, since the previous work was done using a different bacterial strain background and a different host, these genes may have strain-specific or host-specific effects.
In addition to the above-mentioned screen, several groups have assessed the roles of individual genes in colonization, and many of these genes were identified by our screen. Examples of genes previously shown to be compromised for colonization in the mouse model, when mutated and identified as colonization defective in our screen, include the urease genes (16, 22, 62), the nickel transporter nikA gene, the fur gene (10), the bacterial ferritin pfr gene (64), the ß-1-4 galactosyl transferase gene (18), several chemotaxis genes (21, 59), and the flagellar sheath-localized adhesin hpaA gene (11). A hallmark of H. pylori infection is infiltration of neutrophils that produce reactive oxygen and nitrogen species, exposing both the host and the bacteria to DNA and protein damage. Thus, genes involved in protection from oxidative stress or repair of DNA damage might be expected to play an essential role in vivo. We found candidate mutants linked to several stress-associated genes, including superoxide dismutase (sodB) and a peroxyredoxin-thioredoxin-like protein-encoding gene (HP0136), which previous work had shown to be essential (56, 65), and a gene not previously tested, that encoding a heat shock protein (HSP90) homologue (HP0210). Catalase (katA) (25) and
-glutamyltranspeptidase (14, 24), an enzyme involved in glutathione metabolism, showed subtle or variable effects on colonization and persistence in our screen and in the published literature. Similarly, two genes involved in DNA repair, ruvC (34) and the endonuclease III gene (46, 51), gave variable colonization phenotypes at early time points in other studies and gave inconsistent results in our screen. However, transposon mutants mapping to two DNA glycosylases (HP1347 and HP0602) were attenuated in our screen, suggesting that the ability to repair DNA damage is likely very important for persistence in vivo.
Another group of genes that has been investigated in some detail in animal models is the cag PAI which encodes a T4SS that translocates CagA into host cells and induces proinflammatory cytokine production. This pathogenicity island is variably present in human clinical isolates and has been associated with more-severe pathological outcomes. In humans it is thought not to be required for establishing infection since it is not present in all strains. As with human infection, some mouse-colonizing strains do not contain the PAI (60). Murine adaptation of cag+ human clinical isolates often results in attenuation of PAI activity (49), and this is the case for NSH79, one of the strain backgrounds used in our experiments. While many insertions within the PAI did not give a colonization phenotype, we confirmed by construction of independent null alleles that disruption of the HP0529, an Agrobacterium tumefaciens virB8 homologue, and the HP0538 (cagN) genes conferred a competitive defect in the NSH79 background. Others have observed a competitive defect for PAI mutants at early, but not later, time points in the mouse model (35). The mechanism for this attenuation is not clear. There may be an energy cost to producing a nonfunctional PAI, although we could not measure any growth attenuation in vitro. Further analysis of the kinetics of infection and the precise localization of the bacteria may shed light on this puzzle.
One advantage of the microarray method for tracking transposon mutant behavior over most signature-tagged mutagenesis strategies is that it provides information for genes that have been tested and have not given a phenotype. For example, another group evaluated the importance of the Entner-Doudoroff pathway in vivo and found no defect for the mutation of edd (6-phosphogluconate dehydratase) (67). We found similar results for both edd and eda, a gene encoding another enzyme in this pathway. Similarly the amidases amiE and amiF previously were suggested to not be important for colonization when deleted in the SS1 strain background (9), similar to our data for the NSH79 background. In the NSH57 background, however, we observed a phenotype for both genes. The importance of these enzymes has been suggested by the profound role ammonia production plays both in acid protection and as a nitrogen source for H. pylori. Indeed, early testing of urease mutants in vivo showed a requirement for this enzyme even when the stomach pH was neutralized, suggesting that ammonia plays multiple roles during infection (16). AmiE is a highly expressed protein in the H. pylori cell (33, 61), and both amiE and amiF are regulated by pH (10, 38) and iron availability (39). The strain-specific requirement for the amidases may result from genetic variation in the various pathways that govern ammonia production and consumption. This hypothesis will require further testing.
New colonization gene candidates. While our colonization screen was not saturating, we queried more than half of the genome, allowing a global view of the H. pylori cellular processes important for establishing and maintaining infection in the stomach. As detailed above, many of the functional categories identified support results from similar analyses of mucosal surface colonization by other bacterial pathogens. Indeed, specific examples from each of these categories have been investigated on an individual gene basis in H. pylori, and our results confirm and extend these studies. Perhaps the more significant impact of this work is the identification of candidate colonization genes that may be specific to H. pylori infection of the stomach. The largest class of genes that showed attenuated colonization potential when interrupted by transposon insertion consisted of hypothetical genes that may be uniquely important to H. pylori biology. We validated the requirement for eight of these gene loci by making independent null alleles, assaying competition with wild-type bacteria, and in some cases, determining the dose for 50% infection. At this point, the mechanism by which these genes contribute to the infectious process is not clear, and it remains possible that disruption of the targeted genes affects the expression of neighboring genes that mediate the phenotype. In a few cases, we could measure strain-specific effects on urease activity or motility, which are processes known to be required for colonization. One of the genes (HP1028) that gave a strong phenotype in both strain backgrounds we tested was suggested as a putative virulence factor from a bioinformatics-based systems biology analysis of genes linked to those of the cag PAI (57). While that study did not include functional studies, the authors suggest that this putative secreted protein may have toxin-like properties based on protein folding prediction software.
Many of the genes in our overall screen as well as those we independently validated showed strain-specific colonization phenotypes. H. pylori strains differ substantially in both their gene complements and sequence polymorphisms. Thus, it is perhaps not surprising that the genetic requirements for colonization would differ in different strain backgrounds. Furthermore, recent studies have shown that strains differ in their anatomic distribution within the stomach even in the murine model (1), which could contribute to a differential requirement for some genes. It is tempting to focus on those genes required in a multiple strain background as being the most important for colonization. Strain-specific genes, however, may regulate biological pathways that are equally important but show differential phenotypic buffering by the actions of unique genes or alleles in different strains. More-detailed analysis of the genes identified in this screen should yield a deeper appreciation of the pathogenic potential of various H. pylori strains and illuminate additional targets for antimicrobial strategies.
| ACKNOWLEDGMENTS |
|---|
This work was supported by Public Health Service grant AI054423 from the National Institute of Allergy and Infectious Diseases and by a grant from the Pew Charitable Trusts Program in the Biomedical Sciences.
| FOOTNOTES |
|---|
Published ahead of print on 13 November 2006. ![]()
Supplemental material for this article may be found at http://iai.asm.org/. ![]()
| REFERENCES |
|---|
|
|
|---|
| 1. | Akada, J. K., K. Ogura, D. Dailidiene, G. Dailide, J. M. Cheverud, and D. E. Berg. 2003. Helicobacter pylori tissue tropism: mouse-colonizing strains can target different gastric niches. Microbiology 149:1901-1909. |
| 2. | Amieva, M. R., N. R. Salama, L. S. Tompkins, and S. Falkow. 2002. Helicobacter pylori enter and survive within multivesicular vacuoles of epithelial cells. Cell. Microbiol. 4:677-690.[CrossRef][Medline] |
| 3. | Ausebel, F., R. Brent, R. Kingston, D. Moore, J. Seidman, J. Smith, and K. Stuhl (ed.). 1997. Short protocols in molecular biology, 3rd ed. John Wiley & Sons, New York, NY. |
| 4. | Ball, C. A., I. A. Awad, J. Demeter, J. Gollub, J. M. Hebert, T. Hernandez-Boussard, H. Jin, J. C. Matese, M. Nitzberg, F. Wymore, Z. K. Zachariah, P. O. Brown, and G. Sherlock. 2005. The Stanford Microarray Database accommodates additional microarray platforms and data formats. Nucleic Acids Res. 33:D580-D582. |
| 5. | Bijlsma, J. J., A. L. M. Lie, I. C. Nootenboom, C. M. Vandenbroucke-Grauls, and J. G. Kusters. 2000. Identification of loci essential for the growth of Helicobacter pylori under acidic conditions. J. Infect. Dis. 182:1566-1569.[CrossRef][Medline] |
| 6. | Bijlsma, J. J. E., C. M. J. E. Vandenbroucke-Grauls, S. H. Phadnis, and J. G. Kusters. 1999. Identification of virulence genes of Helicobacter pylori by random insertion mutagenesis. Infect. Immun. 67:2433-2440. |
| 7. | Bourzac, K. M., and K. Guillemin. 2005. Helicobacter pylori-host cell interactions mediated by type IV secretion. Cell. Microbiol. 7:911-919.[CrossRef][Medline] |
| 8. | Bourzac, K. M., L. A. Satkamp, and K. Guillemin. 2006. The Helicobacter pylori cag pathogenicity island protein CagN is a bacterial membrane-associated protein that is processed at its C terminus. Infect. Immun. 74:2537-2543. |
| 9. | Bury-Moné, S., S. Skouloubris, C. Dauga, J.-M. Thiberge, D. Dailidiene, D. E. Berg, A. Labigne, and H. De Reuse. 2003. Presence of active aliphatic amidases in Helicobacter species able to colonize the stomach. Infect. Immun. 71:5613-5622. |
| 10. | Bury-Mone, S., J. M. Thiberge, M. Contreras, A. Maitournam, A. Labigne, and H. De Reuse. 2004. Responsiveness to acidity via metal ion regulators mediates virulence in the gastric pathogen Helicobacter pylori. Mol. Microbiol. 53:623-638.[CrossRef][Medline] |
| 11. | Carlsohn, E., J. Nyström, I. Bolin, C. L. Nilsson, and A.-M. Svennerholm. 2006. HpaA is essential for Helicobacter pylori colonization in mice. Infect. Immun. 74:920-926. |
| 12. | Chalker, A. F., H. W. Minehart, N. J. Hughes, K. K. Koretke, M. A. Lonetto, K. K. Brinkman, P. V. Warren, A. Lupas, M. J. Stanhope, J. R. Brown, and P. S. Hoffman. 2001. Systematic identification of selective essential genes in Helicobacter pylori by genome prioritization and allelic replacement mutagenesis. J. Bacteriol. 183:1259-1268. |
| 13. | Chan, K., C. C. Kim, and S. Falkow. 2005. Microarray-based detection of Salmonella enterica serovar Typhimurium transposon mutants that cannot survive in macrophages and mice. Infect. Immun. 73:5438-5449. |
| 14. | Chevalier, C., J. M. Thiberge, R. L. Ferrero, and A. Labigne. 1999. Essential role of Helicobacter pylori gamma-glutamyltranspeptidase for the colonization of the gastric mucosa of mice. Mol. Microbiol. 31:1359-1372.[CrossRef][Medline] |
| 15. | Covacci, A., S. Censini, M. Bugnoli, R. Petracca, D. Burroni, G. Macchia, A. Massone, E. Papini, Z. Xiang, N. Figura, et al. 1993. Molecular characterization of the 128-kDa immunodominant antigen of Helicobacter pylori associated with cytotoxicity and duodenal ulcer. Proc. Natl. Acad. Sci. USA 90:5791-5795. |
| 16. | Eaton, K. A., C. L. Brooks, D. R. Morgan, and S. Krakowka. 1991. Essential role of urease in pathogenesis of gastritis induced by Helicobacter pylori in gnotobiotic piglets. Infect. Immun. 59:2470-2475. |
| 17. | Eaton, K. A., D. Kersulyte, M. Mefford, S. J. Danon, S. Krakowka, and D. E. Berg. 2001. Role of Helicobacter pylori cag region genes in colonization and gastritis in two animal models. Infect. Immun. 69:2902-2908. |
| 18. | Eaton, K. A., S. M. Logan, P. E. Baker, R. A. Peterson, M. A. Monteiro, and E. Altman. 2004. Helicobacter pylori with a truncated lipopolysaccharide O chain fails to induce gastritis in SCID mice injected with splenocytes from wild-type C57BL/6J mice. Infect. Immun. 72:3925-3931. |
| 19. | Eisen, M. B., P. T. Spellman, P. O. Brown, and D. Botstein. 1998. Cluster analysis and display of genome-wide expression patterns. Proc. Natl. Acad. Sci. USA 95:14863-14868. |
| 20. | Fischer, W., J. Puls, R. Buhrdorf, B. Gebert, S. Odenbreit, and R. Haas. 2001. Systematic mutagenesis of the Helicobacter pylori cag pathogenicity island: essential genes for CagA translocation in host cells and induction of interleukin-8. Mol. Microbiol. 42:1337-1348.[CrossRef][Medline] |
| 21. | Foynes, S., N. Dorrell, S. J. Ward, R. A. Stabler, A. A. McColm, A. N. Rycroft, and B. W. Wren. 2000. Helicobacter pylori possesses two CheY response regulators and a histidine kinase sensor, CheA, which are essential for chemotaxis and colonization of the gastric mucosa. Infect. Immun. 68:2016-2023. |
| 22. | Ghiara, P., M. Marchetti, M. J. Blaser, M. K. Tummuru, T. L. Cover, E. D. Segal, L. S. Tompkins, and R. Rappuoli. 1995. Role of the Helicobacter pylori virulence factors vacuolating cytotoxin, CagA, and urease in a mouse model of disease. Infect. Immun. 63:4154-4160.[Abstract] |
| 23. | Guillemin, K. J., and N. R. Salama. 2002. Helicobacter pylori functional genomics, p. 291-320. In B. Wren and N. Dorrell (ed.), Functional microbial genomics, vol. 33. Academic Press, London, United Kingdom.[CrossRef] |
| 24. | Guo, B. P., and J. J. Mekalanos. 2002. Rapid genetic analysis of Helicobacter pylori gastric mucosal colonization in suckling mice. Proc. Natl. Acad. Sci. USA 99:8354-8359. |
| 25. | Harris, A. G., J. E. Wilson, S. J. Danon, M. F. Dixon, K. Donegan, and S. L. Hazell. 2003. Catalase (KatA) and KatA-associated protein (KapA) are essential to persistent colonization in the Helicobacter pylori SS1 mouse model. Microbiology 149:665-672. |
| 26. | Hofreuter, D., S. Odenbreit, G. Henke, and R. Haas. 1998. Natural competence for DNA transformation in Helicobacter pylori: identification and genetic characterization of the comB locus. Mol. Microbiol. 28:1027-1038.[CrossRef][Medline] |
| 27. | Israel, D. A., N. Salama, C. N. Arnold, S. F. Moss, T. Ando, H. P. Wirth, K. T. Tham, M. Camorlinga, M. J. Blaser, S. Falkow, and R. M. Peek, Jr. 2001. Helicobacter pylori strain-specific differences in genetic content, identified by microarray, influence host inflammatory responses. J. Clin. Investig. 107:611-620.[CrossRef][Medline] |
| 28. | Israel, D. A., N. Salama, U. Krishna, U. M. Rieger, J. C. Atherton, S. Falkow, and R. M. Peek, Jr. 2001. Helicobacter pylori genetic diversity within the gastric niche of a single human host. Proc. Natl. Acad. Sci. USA 98:14625-14630. |
| 29. | Jacobs, M. A., A. Alwood, I. Thaipisuttikul, D. Spencer, E. Haugen, S. Ernst, O. Will, R. Kaul, C. Raymond, R. Levy, L. Chun-Rong, D. Guenthner, D. Bovee, M. V. Olson, and C. Manoil. 2003. Comprehensive transposon mutant library of Pseudomonas aeruginosa. Proc. Natl. Acad. Sci. USA 100:14339-14344. |
| 30. | Joyce, E. A., K. Chan, N. R. Salama, and S. Falkow. 2002. Redefining bacterial populations: a post-genomic reformation. Nat. Rev. Genet. 3:462-473.[Medline] |
| 31. | Kavermann, H., B. P. Burns, K. Angermuller, S. Odenbreit, W. Fischer, K. Melchers, and R. Haas. 2003. Identification and characterization of Helicobacter pylori genes essential for gastric colonization. J. Exp. Med. 197:813-822. |
| 32. | Lee, A., J. O'Rourke, M. C. De Ungria, B. Robertson, G. Daskalopoulos, and M. F. Dixon. 1997. A standardized mouse model of Helicobacter pylori infection: introducing the Sydney strain. Gastroenterology 112:1386-1397.[CrossRef][Medline] |
| 33. | Lock, R. A., S. J. Cordwell, G. W. Coombs, B. J. Walsh, and G. M. Forbes. 2001. Proteome analysis of Helicobacter pylori: major proteins of type strain NCTC 11637. Pathology 33:365-374.[Medline] |
| 34. | Loughlin, M. F., F. M. Barnard, D. Jenkins, G. J. Sharples, and P. J. Jenks. 2003. Helicobacter pylori mutants defective in RuvC Holliday junction resolvase display reduced macrophage survival and spontaneous clearance from the murine gastric mucosa. Infect. Immun. 71:2022-2031. |
| 35. | < |