Previous Article | Next Article ![]()
Journal of Bacteriology, February 2005, p. 1276-1286, Vol. 187, No. 4
0021-9193/05/$08.00+0 doi:10.1128/JB.187.4.1276-1286.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Brian LeCuyer,1
Anne G. Lenich,2,
Matthew P. Lazio,1,
Donna Perkins-Balding,2
H. Steven Seifert,1 and
Anna C. Karls2,3*
Department of Microbiology and Immunology, Feinberg School of Medicine, Northwestern University, Chicago, Illinois,1 Department of Microbiology and Immunology, Emory University School of Medicine, Atlanta,2 Department of Microbiology, University of Georgia, Athens, Georgia3
Received 10 June 2004/ Accepted 4 November 2004
|
|
|---|
|
|
|---|
Gonococcal pilus variation is mediated mainly by high-frequency, nonreciprocal, homologous recombination reactions between silent copies of pilin information and the pilin expression locus pilE. Nineteen individual copies of silent pilin information are separated into six loci (pilS) throughout the chromosome of strain FA1090 (23). Recombination between pilE and pilS silent copies occurs at short regions of identity and results in alterations in the coding sequence of the pilE gene (9, 10, 26, 36). The multiple donor silent copies, in conjunction with the frequent regions where crossovers between the recombining genes can occur, contribute to the large repertoire of possible expressed pilin monomers and hence antigenically distinct pili that a single strain of the gonococcus (Gc) can produce. Most of the factors known to be required for antigenic variation at pilE are part of the cell's homologous recombination machinery. This includes recA (13), recX (34), recO, recQ (18), and recJ (32). All of these gene products also play a role in DNA repair. In addition to genes involved in homologous recombination, a potential link to cell division has been demonstrated with the reported requirement of the cell division-associated gene rdgC in pilin variation (17). All homologous recombination processes have been shown to increase under iron-limited conditions in N. gonorrhoeae (31), underscoring the possibility that a global regulatory system controlling recombination exists in this organism.
A number of characteristics of the recombination reactions responsible for pilin variation cannot be explained by homologous recombination reactions alone. First, the regions of identity between recombining segments have been shown to be as short as 7 bp (11). Second, the nonreciprocal nature of the recombination reaction, which is reflected in the exchange of DNA sequences between pilE and pilS copies resulting only in gene conversion at pilE, is not consistent with a mechanism mediated solely by generalized recombination functions (11). Finally, the 3' region of all pilin loci contains sequences reported to be similar to defined site-specific recombinase binding sites (43). Taken together, these observations support the hypothesis that one or more recombinases act in cooperation with the homologous recombination machinery to mediate gonococcal pilin antigenic variation.
Moraxella lacunata and Moraxella bovis are ocular pathogens of humans and cattle, respectively. These members of the family Neisseriaceae elaborate type IV pili which are phase or antigenically variable. Previous work has demonstrated that antigenic and phase variation of the M. bovis and M. lacunata pili is mediated exclusively by a site-specific recombinase known as Piv, for pilin inversion protein (15, 16). The PivMB and PivML proteins belong to a novel family of DNA recombinases known as the Piv/MooV family; the family is unusual because the related recombinases include both site-specific recombinases and DNA transposases (40). Many of the recombinases of this family mediate transposition of insertion sequence (IS) elements from the IS110 family, including MooV, the transposase encoded by IS492. IS110-related elements are atypical insertion sequences because they lack terminal inverted repeat sequences and because they do not appear to generate target site duplications (4). Structure-function studies and molecular modeling with the defining members of this recombinase family, Piv and MooV, have suggested that all members of this family use a similar mechanism for recombination (40; C. Carpenter, D. Perkins-Balding, and A. C. Karls, unpublished data). A gene predicted to encode a protein with weak identity (18%) and similarity (34%) to Piv was identified in N. gonorrhoeae in a screen for proteins capable of inverting Moraxella pilin DNA in Escherichia coli. This gene was named gcr for gonococcal recombinase, and a loss-of-function mutant revealed that this gene product is dispensable for pilin antigenic variation (24). More recently, Southern blot analysis has revealed the presence of multiple genes encoding proteins with significant identity to Piv and MooV in the chromosome of numerous strains of N. gonorrhoeae but not the commensal bacteria Neisseria subflava and Neisseria lactamica (3). The gonococcus is not the only member of Neisseria species with Piv/MooV homologues, as three potential members of the Piv/MooV family, PivNM-1A, PivNM-1B, and PivNM-2, have been identified in the human pathogen Neisseria meningitidis strain MC58 genome sequence (38). Like N. gonorrhoeae, N. meningitidis is capable of pilin antigenic variation (39).
In this work, we identify eight Piv/MooV homologues encoded in the genome of N. gonorrhoeae strain FA1090 and confirm that they are found only within the pathogenic Neisseria. We named these genes invertase-related genes 1 through 8 (irg1-8) based on their similarity to the Piv proteins. We demonstrate that at least some of these genes are expressed; however, a mutant strain inactivated for all eight genes and a strain that overexpresses one irg gene still exhibit wild-type levels of pilin antigenic variation, DNA transformation, and DNA repair. While biochemical assays for complementation of Piv invertase and MooV transposase reactions did not prove the recombinase activities of Irg proteins, comparative analyses of the DNA flanking the coding sequences support the hypothesis that the Piv homologues are transposases for two new members of the IS110 family of insertion sequences.
|
|
|---|
[F endA1 hsdR17 supE44 thi-1 recA1 gyrA96 relA1
(lacZYA-argF) U1691
80 dlacZ
M15] and HMS174(DE3) [recA1 hsdR Rifr gal
(DE3)] were grown in Luria-Bertani broth or agar at 37°C and used to propagate plasmids. Solid medium contained 15 g of agar per liter. FA1090 and VD300 (a derivative of strain MS11) are gonococcal laboratory strains. The remaining gonococcal strains used in this study were collected directly from patients presenting to an urban sexually transmitted disease clinic with gonorrhea. Chromosomal DNA preparations from N. gonorrhoeae laboratory strains FA19, RD5, UU1, and F62; clinical isolates 1084, 1918, 1349, 1402, 1384; N. meningitidis strains NMB and GA0929 (clinical isolate, group Y); and commensal strains Neisseria cineria and N. lactamica were graciously provided to A.C.K. by William Shafer. Gonococcal strains were grown on Gc medium base (Difco) plus Kellogg supplements (GCB) [22.2 mM glucose, 0.68 mM glutamine, 0.45 mM cocarboxylase, 1.23 mM Fe(NO3)3 (all from Sigma)] (12) at 37°C in 5% CO2. FA1090 recA6 strains were created by transformation of FA1090 with plasmid DNA carrying an isopropyl-ß-D-thiogalactopyranoside (IPTG)-regulatable gonococcal recA allele, recA6, which allows control of recA and subsequent control of recA-dependent recombination and repair processes, including pilin antigenic variation (27). An IPTG (Diagnostic Chemicals Limited) concentration of 1 mM in the medium allows for maximal induction of recA transcription and restoration of transformation competence to near-wild-type levels (27). Antibiotics were added at the following concentrations for E. coli: chloramphenicol, 30 mg/liter; erythromycin, 250 mg/liter; kanamycin, 40 mg/liter; ampicillin, 100 mg/liter; and tetracycline, 12.5 mg/liter. For N. gonorrhoeae, antibiotic concentrations were as follows: chloramphenicol, 2 mg/liter, erythromycin, 1 mg/liter, kanamycin, 50 or 250 mg/liter, and nalidixic acid, 2 mg/liter. DNA and RNA manipulations and analyses. Standard procedures were performed as described previously by Sambrook et al. (25). Plasmid DNA was isolated from strains of E. coli by using QIAGEN (Valencia, Calif.) plasmid kits. Enzymes were used according to manufacturers' directions (Promega Corp. and New England Biolabs [NEB]). E. coli strains were transformed by electroporation using the Gene Pulser II electroporation system (Bio-Rad Laboratories) according to the manufacturer's specifications. For Southern blot analysis, DNA was transferred to a Magnagraph nylon membrane according to the manufacturer's specifications (Micro Separations Inc.). Sequencing reactions were performed by using the Big-Dye Terminator Cycle Sequencing kit (Perkin-Elmer Corp.), and sequencing products were separated with an ABI Model 377 automated DNA sequencer. DNA sequence analysis was performed by using Lasergene software (DNASTAR, Inc.) and VectorNTI software (Informax, Inc.). PCR fragments were purified by using a Qiaquick PCR purification kit (QIAGEN). The BLAST program (1) was used to search the National Center for Biotechnology Information (NCBI) nonredundant database. pilE sequences were determined by amplifying pilE from the chromosome with primers PILRBS (5'-GGCTTTCCCCTTTCAATTAGGAG-3') and SP3A (5'-CCGGAACGGACGACCCCG-3') (29) using Taq polymerase and sequencing the resulting PCR product with CONST-F2 (5'-TACCAAGACTACACCGCCCG-3') (29). RNA isolation was performed by using the RNeasy Mini protocol for isolation of total RNA from bacteria (QIAGEN). DNase treatment was performed by using OnColumn DNase digestion with the RNase-Free DNase set (QIAGEN). A260 and A280 levels were taken with a spectrophotometer to assess purity. RNA was run in a 0.8% agarose gel to assess purity, and 28S rRNA bands were quantified with ImageQuant software (Molecular Dynamics, Piscataway, N.J.). RNA folding analysis was performed by using a derivative of the Vienna Folding package from the DNASTAR package. cDNA synthesis was performed as described previously by Serkin and Seifert (30). The primer used for the synthesis of cDNA was LCirgrev (5'-GACCGGCAATACTGCGC-3'). Real-time PCR was performed with a LightCycler instrument (Roche Diagnostics, Indianapolis, Ind.) using a fluorescence resonance energy transfer (FRET) detection method (45, 46). FRET requires the use of two probes: the first (donor) probe is labeled with a fluorophore at the 3' end, and the second (acceptor) probe is labeled with a fluorophore at the 5' end. When both probes properly anneal to target DNA so that they are in close proximity, excitation of the donor results in emission, which in turn excites the acceptor. In this manner, the ratio of acceptor fluorescence to donor fluorescence is measured. Hybridization probes were created in accordance with previously published guidelines (7). The distance between the donor and acceptor probes was 1 bp. The probes used were LCHybirg1 (5'-CGTAGTGAACCCGCTGAAAATAAGCAAGTATGC-3') and LCHybirg2 (5'-GAAAGCAGGTTCAAGCGAACCAAAACAG-3') for the irg genes and LC HybrecA 1 (5'-CGGGCGCATCGTCGAAATCTTCGG-3') and LC HybrecA 2 (5'-CCCGAATCCTCCGGCAAAACCACA-3') for recA. All donor probes were labeled with fluorescein isothiocyanate at the 3' end, while all acceptor probes were labeled at the 5' end with LightCycler-Red 640-N-hydroxysuccinimide ester. Acceptor probes also had phosphate groups attached to the 3' end to prevent probe extension (IT Biochem). Primer and probe melting temperatures were calculated by using the freeware program TM Utility, version 1.3 (Idaho Technology [http://www.idahotech.com/downloads_up/Tmutility_form.htm]). Mixtures in hybridization experiments contained 2 µl of LightCycler DNA Master Hybridization Probes mix (Roche Diagnostics), 3.2 µl of MgCl2 (for a final concentration of 5 mM), 2 µl of each primer (0.5 mM), 2 µl of each hybridization probe (0.2 mM), 2 µl of template, and PCR-grade sterile water for a final volume of 20 µl. Fluorescence was acquired at the end of every extension step during the PCR. Primers used for PCR were LCirg (5'-CGTAGTGAACCCGCTGAAAATAAGCAAGTATGC-3') and LCirgrev (5'-GACCGGCAATACTGCGC-3').
Identification, cloning, and inactivation of the irg copies. The Gonococcal Sequencing Project database (23) was used to identify the individual irg genes in the chromosome of N. gonorrhoeae strain FA1090. BLASTN analysis of the complete gonococcal genome was performed by using the pivNG sequence (3). Sequences greater than 54 bp with E values of 1056 or less were concluded to be Piv homologues and hence irg copies. This BLASTN analysis did not produce any hits longer than 100 bp besides the eight irg genes.
The following mutant irg genes were created and introduced into FA1090 recA6 in the order presented to create a strain with all eight irg genes knocked out. irg5, irg6, and irg8 were PCR amplified by using either PivDir (5'-TCGGGATTTACGCCGATTTG-3') with PivPrimer (5'-ATAAGGGTAATCCCTG-3') for irg6 or PivFor (5'-CGCCGAAAGGAACGTGTATGCT-3') with PivRev2 (5'-CGGATTTCAGACGCGGCAAAGCA-3') for irg5 and irg8 and cloned into the pBlunt vector (Invitrogen). Chloramphenicol and kanamycin resistance cassetts were then inserted into SnaBI and NdeI sites internal to irg5 and irg6 genes, respectively. A kanamycin resistance cassette was inserted into irg8 at the NdeI site. These vectors were then transformed into FA1090 recA6, and the appropriate antibiotic resistance was selected to inactivate the individual irg copy by insertion. irg7 was inactivated by using transposon shuttle mutagenesis (28), resulting in insertional inactivation of irg7 with an erythromycin minitransposon. irg1 and irg4 were inactivated through the creation of in-frame deletions. To create an in-frame deletion in irg1, a 4,590-bp segment of DNA containing irg1 was amplified from the chromosome of N. gonorrhoeae by using primers PivNG3B-For (5'-ACCGCAGCAAACCCTCTCCCTAAC-3') and PivNG3B-Rev (5'-CGAATACGACATCCGCCCAGAAGA-3'). This fragment was cloned into pBlunt (Invitrogen), and a 357-bp fragment internal to irg1 was removed by double digestion with AgeI and NdeI. The resulting construct was transformed into N. gonorrhoeae, and successful deletions were screened with primer PivNG3BScreen (5'-ATCAGTTTTCGGACCGGTATGCG-3'), which spans the newly formed junction. To create an in-frame deletion in irg4, primer pairs were created to generate products that when ligated together would result in the removal of 568 nucleotides internal to irg4. The primers used to amplify the 3' portion of irg4 including downstream sequence were PivNG4Right-For (5'-GAATCCGGGACAAGCGTAAGGGG-3') and PivNG4Right-Rev (5'-CCTGAAAGTCAGTGCCGATGCCG-3'), and the primers used to amplify the 5' portion of irg4 including upstream sequence were PivNG4Left-For (5'-TGAGTCAATACCTGTCGCAACGCC-3') and PivNGLeft-Rev (5'-TCCCAAGCTTTTCAACCGGTCCG-3'). The two products ligated together were cloned into pCR-BluntII (Invitrogen) and transformed into N. gonorrhoeae, and successful deletions were screened with primer PivNG4-Screen (5'-GGTTGAAAAGCTTGGGAGAATCCG-3'). irg3 was inactivated by deleting part of the gene and inserting a spectinomycin resistance cassette, pHP45
(22). A 2,278-bp segment of DNA which included irg3 was amplified from the FA1090 chromosome by using primers Irg 3 Upstream (5'-TCG TAA TGG GCG GAA AGT CC-3') and Irg 3 Downstream (5'-CGC TAA TGG GCT TCA GAC GG-3'). The PCR product was cloned into pCR-Blunt II (Invitrogen). A 571-bp fragment internal to irg3 was dropped out by double digestion with AgeI and AflII (NEB), and the vector was treated with T4 DNA polymerase (NEB). pHP45
was digested with SmaI (NEB), and the fragment containing
Spc was ligated into the blunted pCR-BluntII-irg3 vector. The resulting construct, pCR-Blunt irg3::Spc was then sequenced. FA1090 recA6 irg145678 was transformed with this plasmid, and transformants were selected on Spec45 plates. To create an in-frame deletion in irg2, primer sets were created to generate products that when ligated together would result in the removal of 764 bases internal to irg2, leaving 199 bp of the gene with an intervening 8-bp PacI site. The primers used to amplify the 5' end of irg2, including upstream sequence, were Irg2UP Uptake, which contains the Gc uptake sequence (5'-TTG CCG TCT GAA GCA ATG AGG GCG GTA CAG G-3') and Irg2 Pac I, which contains a PacI site (5'-GAT TAA TTA AAT CTA AAC CTT TTG AAT CGT TG-3'). This 835-bp PCR fragment was cloned into pCR-Blunt II (Invitrogen). For the 3' end of irg2 including downstream sequence, a 1,384-bp fragment was amplified using primers Irg 2 Pac I-1599, containing a PacI site (5'-TTA ATT AAA GCA TTG ATG CGT AAA CTC G-3'), and Irg 2 Downstream (5'-GCA CAT GCT ATT CAA ATC AAG G-3'). This PCR product was also cloned into pCR-Blunt II, and then each plasmid was digested with PacI and PstI (NEB). The fragment containing the 3' portion of the gene was gel isolated and subcloned into the linear plasmid containing the 5' portion of irg2. The resulting plasmid was sequenced and transformed into FA1090 recA6 irg1345678. The colonies were screened for the deleted irg2 gene by using a primer that spanned the newly formed junction containing the PacI site, IRG2INFRAME (5'-AGG TTT AGA TTT AAT TAA A-3'). Successful combination of all eight mutations in one strain was confirmed by Southern blot and PCR analysis (data not shown).
To create a strain to be used in overexpression analysis, irg5 was amplified and cloned into the Neisserial Insertional Complementation System vector pGCC6 (32), which allows for IPTG-inducible expression of irg5. This construct was then transformed into FA1090 recA6 at a site in the gonococcal chromosome irrelevant to pilin antigenic variation to create FA1090 recA6 NICS6::irg5.
In vivo inversion and transposition assays. Expression plasmids were created for two of the irg genes. irg2 (pAG712) and irg7 (pAG711) were amplified with Pfu polymerase (Promega) by using primer pairs 5'-CAGTTCAGCTAGCCGTAACGCCGTAGGATTGG-3' and 5'-CTCGTCTCGAGCGCCGATTTGTAACGCGATGG-3' as well as 5'-CAGTTCACATATGAATATAATCGGGCCGGACATC-3' and 5'-CTCGTCTCGAGTTGATTCAATCGGTGTCTTTCC-3' introduced into the NheI/XhoI and NdeI/XhoI sites, respectively, of pET21a (Novagen) and expressed from the T7 promoter as His6-tagged proteins in E. coli HMS174(DE3). The pivML (pAG1300) and mooV (pAG921) expression vectors, containing the recombinase genes in pET21a between the NdeI/XhoI sites, have previously been shown to exhibit in vivo inversion and precise excision activities, respectively (21, 41). Recombinase expression was induced with 0.2 mM IPTG, and samples were taken at 2 h postinduction from both induced and uninduced cultures for Western blot analysis as described previously by Tobiason et al. (40). Anti-His6 rabbit antibody (Santa Cruz Biotechnology) was used for the primary antibody.
The in vivo inversion assay was performed essentially as described previously by Tobiason et al. (40). The inversion substrate pAG862 and each recombinase expression vector were introduced into HMS174(DE3), and transformants were grown in Luria broth for 12 h without IPTG or for 5 and 12 h after induction at an optical density at 600 nm of 0.6 with 0.005 mM IPTG. The lower level of induction, compared to that of the cultures used for Western blot analysis, was necessary due to toxicity of high-level expression of the recombinases. Isolated plasmid DNA from these cultures was restricted with KpnI and SalI and electrophoresed on a 1% agarose gel. Southern blot analysis of the inversion-diagnostic restriction fragments utilized probe prepared from pMxL1 (16), which contains the full invertible segment from M. lacunata. The 32P-labeled probe was prepared by PCR as described previously by Perkins-Balding et al. (21). Three primers were used to generate probes for both orientations of the invertible segment: 5'-CGGAATTCGCTTGCAATACTGGCTCTTTAC-3', 5'-TATGCCCGCTGTCATAGTCAGCT-3',and 5'-GAGGTGATTTCGCATAATATA-3'.
The in vivo transposition assay was performed as described previously (21), using pAG921, pAG711, and pAG712 to complement the transposition substrate for precise excision and circle junction formation. The transposition substrate pAG713 was constructed by first replacing the BamHI and XhoI segments of pACYC177 with the IS492-containing BamHI/XhoI fragment from pAG949 (21). Sequence between the MfeI sites of mooV was then replaced with the chloramphenicol acetyltransferase gene (cat) from pACYC184 (PCR amplified with Pfu DNA polymerase and the primers 5'-CCTACAATTGACGGAAGATCACTTCGCAGAA-3' and 5'-GCGCAATTGAGGGCACCAATAACTGCCTT-3') to generate pAG713. The induction of recombinase expression was done with 0.005 mM IPTG; induction of higher levels of His6-tagged MooV inhibited the excision of IS492 (21).
Kinetic pilus-dependent colony variation assay. The kinetic colony phase variation assay (E.V. Sechman and H. S. Seifert, unpublished data) was performed as follows. N. gonorrhoeae strains were streaked for isolated colonies on GCB plates containing 1 mM IPTG and incubated at 37°C at 5% CO2. Isolated colonies were screened at 18, 20, 22, 24, and 26 h for the presence of sectors of variation which are indicative of pilus variation. Scores ranging from 0 to 4 were assigned to colonies based on the number of sectors of variation present per colony. Colonies that contained more than four sectors of variation were given a score of 4. Analysis was performed on at least 20 colonies per screen per experiment, and the assay was performed at least three times per strain.
UV sensitivity assay. The UV sensitivity of Gc strains was assayed as previously described (14). Gonococci were incubated for 24 h on GCB plates (in the presence of IPTG for recA6 and Neisserial Insertional Complementation System derivatives), collected with a Dacron swab, and suspended in GCB liquid. Serial dilutions were made and exposed to UV radiation at 0 to 80 J/m2 in a Stratagene UV Stratalinker 1800. Relative resistance was expressed as percent survival compared to the number of CFU in the unirradiated culture.
DNA transformation competence assay. Gonococcal strains were grown on GCB medium overnight and transferred by a Dacron swab into 37°C GCB liquid plus MgCl2 containing uncut pSY6 plasmid DNA (33), which carries a point mutation in the gonococcal gyrB gene that confers nalidixic acid resistance when recombined into the gyrB locus. After a 15-min incubation, unincorporated DNA was degraded with RQ1 DNase I (Promega) for 15 min, and the cells were diluted into GCB, incubated for 4 h, and plated for nalidixic acid-resistant and total CFU.
|
|
|---|
![]() View larger version (24K): [in a new window] |
FIG. 1. Chromosomal location of irg copies in N. gonorrhoeae strain FA1090. A schematic diagram of the N. gonorrhoeae chromosome is shown. pil loci are labeled, and irg copies are denoted by black arrows that point in the predicted direction of transcription. The chromosomal position of the start codon for each irg copy is given for the linear sequence of the N. gonorrhoeae strain FA1090 genome, starting at dnaA (GenBank accession no. AE004969 [23]).
|
![]() View larger version (89K): [in a new window] |
FIG. 2. Alignment of Irg proteins with Piv and homologues. Conserved residues are boxed in gray. Asterisks are above residues shown to be required for PivML-mediated inversion of Moraxella DNA in E. coli (40). PivNM-1A, PivNM-1B, and PivNM-2 are Piv homologues found in N. meningitidis. MooV is the transposase for IS492 in P. atlantica (21) and IS621 Tnp is the transposase for IS621 in E. coli (5).
|
To determine whether strains of N. gonorrhoeae, in addition to FA1090, contain multiple copies of the irg genes, Southern blot analysis was performed on the laboratory strains FA19, UU1, RD5, and F62 as well as five isolates from patients presenting with disseminated gonococcal infections (Fig. 3). The probe, which was generated by PCR from the full irg2 gene, should hybridize to irg1, irg2, irg3, irg4, irg5, irg6, and irg8 but should not hybridize to irg7 or to the divergent N. meningitidis Piv homologues (Fig. 3, lanes 10 and 11). The probe specificity was confirmed by the lack of hybridization with piv sequence on M. lacunata chromosomal DNA or pAG702 (41) (Fig. 3, lanes 14 and 15), while strong hybridization was seen with pAG852, in which irg5 replaces piv on the pAG702 vector (lane 16). This Southern blot analysis revealed the same number of irg gene copies on similarly migrating fragments in both laboratory isolates (Fig. 3, lanes 1 to 4) and minimally passaged clinical isolates (lanes 5 to 9). The only exceptions were laboratory strain UU1 (lane 2) and one of the clinical isolates (lane 8), which each showed one fragment with a different mobility. These shifts in restriction fragment mobility most probably reflect restriction fragment polymorphism but could also represent DNA rearrangements. An additional nine epidemiologically unrelated clinical isolates of N. gonorrhoeae were analyzed by Southern blot and found to exhibit an identical distribution of the irg copies to one another (data not shown). Therefore, the number and chromosomal positions of the irg genes are extremely well conserved within this species. The commensal Neisseria species N. cineria and N. lactamica (Fig. 3, lanes 12 and 13) do not have any DNA sequences that hybridize with the irg probe. Therefore, the irg genes are found only in the pathogenic Neisseria species.
![]() View larger version (78K): [in a new window] |
FIG. 3. Distribution of irg copies in various isolates of N. gonorrhoeae. Southern analysis was performed with chromosomal DNA preparations that were digested with ClaI, electrophoresed on a 0.8% agarose gel, transferred to a nylon membrane, and probed with denatured 32P-labeled PCR product corresponding to the irg2 gene. N. gonorrhoeae strains are as follows: lane 1, FA19; lane 2, UU1; lane 3, RD5; lane 4, F62; lane 5, 1084; lane 6, 1918; lane 7, 1349; lane 8, 1402; lane 9, 1384. N. meningitidis strains are as follows: lane 10, GA0929; lane 11, NMB. Neisseria commensals are in lanes 12 (N. cineria) and 13 (N. lactamica). Negative control lanes are lanes 14 (M. lacunata) and 15 (pivML clone pAG702). The positive control, irg5 clone pAG852, is in lane 16.
|
![]() View larger version (44K): [in a new window] |
FIG. 4. Assays for complementation of Piv invertase and MooV transposase activities by Irg2 and Irg7 in E. coli. (A) Western blot analysis of His6-tagged MooV (MooV-H6), PivML (Piv-H6), Irg2 (Irg2-H6), and Irg7 (Irg7-H6) expression. Whole-cell lysates from induced (I) and uninduced (U) cultures were assayed on the Western blot as described in Materials and Methods. The proteins range in size from 318 to 322 amino acid residues without the His6 tag; the arrow marks the approximate location of the tagged Irg proteins. The molecular weights (in thousands) of protein standards (S) are indicated. (B) His6-tagged Irg7, Irg2, and the invertase PivML were assayed for invertase activity in HMS174(DE3) by using the Piv inversion substrate (pAG862), which contains the M. lacunata pilin invertible segment. Expression of the recombinases was induced and inversion products were isolated as described in Materials and Methods. Southern blot analysis of the pAG862 KpnI/SalI restriction fragments which are diagnostic for inversion of the pilin segment is shown. Arrows indicate the restriction fragments that are unique to the plasmid substrate that has undergone Piv-mediated inversion of the pilin segment, as illustrated. Tc, tetracycline. (C) His6-tagged MooV, Irg2, Irg7, and PivML were assayed for transposase activity in HMS174(DE3) by using the plasmid substrate pAG713, containing IS492 mooV::cat that is defective for excision due to the substitution of cat for most of the mooV transposase gene. Excision of IS492 mooV::cat to yield a circular element is detected in a PCR-based assay using primers (arrows) directed outward from the ends of the element which amplify the circlejunction. The 405-bp PCR product is indicated on the agarose gel by the arrow. Lane M contains pBR322/MspI (sizes in base pairs are shown). Ap, ampicillin; Kn, kanamycin.
|
mooV::cat to yield the circular form of the element, but neither of the Irg proteins was capable of complementing for MooV-mediated excision of the IS element (Fig. 4C). The irg genes are not involved in pilin variation or general recombination. To determine whether the Irg proteins may play a role in a recombination or repair process of the gonococcus, all eight irg genes were inactivated in FA1090 recA6 by using a mixture of antimicrobial resistance cassettes, minitransposons, and in-frame deletions to create strain FA1090 recA6 irg1-8, as confirmed by PCR and Southern blot analyses (data not shown). In addition, to determine whether overexpression of an irg gene alters recombination or repair processes, irg5 was placed under Plac control and introduced into an intergenic region on the chromosome of strain FA1090 recA6, creating FA1090 recA6 NICS6::irg5. The irg mutant and overexpressor strains were generated in a recA6 strain, in which recA expression is controlled by LacI (27), so that initiation of antigenic variation could be regulated. Because different pilE sequences can lead to changes in pilin-dependent phenotypes, including antigenic variation, the pilE genes of these various constructs were sequenced. All genes were found to be identical, ensuring that any observed differences were not caused by different starting pilE sequences.
Reverse transcription (RT)-PCR analysis of FA1090 recA6 total RNA showed that at least one of these irg genes produces a detectable transcript (data not shown). A quantitative RT-PCR assay for irg mRNA levels demonstrated that the overexpressing strain (FA1090 recA6 NICS6::irg5) showed an approximately sevenfold increase in total mRNA compared to that of the parental strain (FA1090 recA6) (Fig. 5). Therefore, the parental strain FA1090 recA6 irg1-8 and FA1090 recA6 NICS6::irg5 were used to determine whether elimination or overexpression of irg affected pilin antigenic variation.
![]() View larger version (29K): [in a new window] |
FIG. 5. Quantitative real-time PCR analysis of irg transcript levels in the overexpressed strain. Transcript levels were measured by using FRET probes in a LightCycler (Roche). Black bars are levels of irg transcripts, and gray bars are levels of recA transcripts included as a control. Total RNA was isolated from FA1090 recA6, in which recA gene expression is inducible by IPTG, and FA1090 recA6 NICS6::irg5, which overexpresses irg5 upon induction with IPTG.
|
![]() View larger version (16K): [in a new window] |
FIG. 6. Inactivation or overexpression of irg genes does not affect pilin antigenic variation, DNA transformation, or DNA repair. (A) Kinetic variation assay of N. gonorrhoeae measuring the number of sectors of a colony that undergo pilin-dependent colony variation over time. Error bars represent the standard errors of the means of at least four experiments analyzing 20 individual colonies per experiment. (B) DNA transformation competence measuring the introduction of a point mutation in gyrB that confers NAL resistance. Error bars represent the standard deviations from the means for at least four replicate experiments. (C) Relative survival after irradiation with UV light. Error bars represent the standard deviations from the means for four experiments.
|
Analyses of the sequences flanking the Irg coding sequences. As members of the Piv/MooV family, the Irg proteins are likely to function as site-specific recombinases or as transposases for IS elements belonging to the IS110 family. Because our results indicate that the Irg proteins are not generalized recombinases, we investigated the possibility that the irg genes and their flanking sequences define IS elements. An unusual feature of most IS110-related elements is the absence of both target site duplication and terminal inverted repeats (4), making it difficult to define the ends of possible elements. Alignments of chromosomal sequences upstream and downstream of the irg genes, performed with the MegAlign tools of LaserGene Navigator, revealed that four (irg2, irg3, irg5, and irg6) of the eight irg genes are adjacent to long repeated sequences, encoding putative filamentous phage proteins (Fig. 7A). Preference for insertion next to repeated sequences or site-specific insertion has been reported for other members of the IS110 family (2, 21). The alignment of the flanking sequences for the irg genes reveals 14 bp of identity immediately upstream of the predicted start codons for irg2, irg3, irg4, irg5, irg6 (irg2-6), and irg8 and 133 bp of near-identity (1-bp difference in the irg8 flanking sequence) in the sequence downstream of the stop codons for irg1-6 and irg8 (Fig. 7B). The Irg proteins encoded by irg2-6 and irg8 are >99.3% identical in their primary amino acid sequence and, based on the homology in the flanking DNA sequence, are probably the transposases for an IS element belonging to the IS110 family. The name attributed to this putative element is ISNgo2 (http://www-is.biotoul.fr/IS.html). The termini of ISNgo2 cannot be definitively assigned until movement of the element has been characterized. However, because the putative ISNgo2 element containing irg2 is inserted in the repeated prophage sequence in the opposite orientation of the irg3-, irg5-, and irg6-associated ISNgo2 copies, the conserved phage sequence helps to delimit the left- and right-end sequences of the element (Fig. 7B and C). As illustrated in Fig. 7C, the 5'-ATA of the conserved left-end sequence appears to be part of the prophage target sequence and not part of ISNgo2. This definition of the left end makes ISNgo2 1,107 bp long with very short terminal inverted repeats of two base pairs (5'-AG). Possible target site duplications are not evident for five of the six copies, indicating that like other IS110-related elements, insertion of the element does not generate a duplication of target sequence. Examination of the conserved left-end sequence for a potential ribosome binding site for irg identified the sequence AGGG, ending 7 bp before the irg start codon. This sequence can base pair with the 3' end of N. gonorrhoeae 16S RNA (8) but includes one GU base pair (Fig. 7B). An additional AU base pair is contributed by the flanking sequence.
![]() View larger version (55K): [in a new window] |
FIG. 7. Sequences flanking the irg copies and definition of ISNgo2 and ISNgo3. (A) The relative order and orientation of the irg genes are shown with flanking genes; gaps are indicated by //. The irg1 and irg8 genes are inserted into sequences predicted to encode a conserved hypothetical protein and a nitroreductase, respectively (dark gray arrows). Putative prophage genes zot, rstA, and tspB (light gray arrows) are described in Discussion; this repeated sequence is nearly identical over 6,874 bp upstream of irg3 and irg5, 4,836 bp downstream of irg2, and 3,319 bp upstream of irg6. Potential prophage sequence associated with irg7 differs from the sequence associated with other irg genes but still encodes a filamentous phage structural protein with a Zot domain (zot2). Parentheses indicate truncated genes or sequences. hrpA is a helicase. (B) Thirty-six base pairs of sequence 5' to the start codon and 59 bp of sequence (with a 100-bp gap of identical sequence) 3' to the stop codon of irg1-6 and irg8 are grouped together and aligned. Thirty-seven base pairs of sequence upstream from the start codon and 59 bp of sequence (with a 119-bp gap in which 3 bp of irg7-associated sequence differs from those of pivNM-1A and pivNM-1B) downstream of the stop codon of irg7, pivNM-1A, and pivNM-1B are aligned. The completely conserved base pairs are highlighted in gray; the proposed left and right ends for ISNgo2 and ISNgo3 are boxed with a black line. Interaction of the 3' end of the N. gonorrhoeae 16S RNA with the proposed ribosome binding sites for the irg genes of ISNgo2 and ISNgo3 is shown below the aligned sequences. dRS3 repeat sequence is underlined, and the flanking prophage sequences are italicized. (C) The putative IS elements are shown with flanking sequences to illustrate their insertion relative to the repeated prophage sequences (italicized). The irg genes of ISNgo2 are black boxes, the left ends are white boxes, and the right ends are horizontally hatched boxes. The ISNgo3 elements, with the irg and pivNM genes as gray boxes and the left and right ends as white and diagonally hatched boxes, respectively, are inserted in the same orientation relative to prophage sequence (see above).
|
BLASTN searches of bacterial sequences in the NCBI database using each of the irg genes and their flanking sequences as query sequences revealed pivNM-1A and pivNM-1B and their surrounding sequences, including 1,874 bp of upstream and 155 bp of downstream sequence, as nearly identical to irg7. There were no other DNA sequences that showed significant homology to the other irg genes and their flanking sequences at the DNA level. As described above, while all the Irg proteins possess approximately the same level of homology with the Piv/MooV family proteins (
20 to 30% identity), Irg2-6 and Irg8 all share >99.3% amino acid identity with each other but only 40% identity with Irg7 (Fig. 2). Also, the sequence flanking irg7 is not homologous to that of irg2-6 and irg8. Alignment of irg7, pivNM-1A, and pivNM-1B and their flanking sequences predicts a putative 1,124-bp IS110-related element, designated ISNgo3 (http://www-is.biotoul.fr/IS.html) (Fig. 7B). The left terminus of this element is uncertain but is currently defined by homology with the left end of the ISNgo2 element. There are no terminal inverted repeats and no apparent target site duplications.
Comparison of the sequences flanking the putative elements reveals that both ISNgo2 and ISNgo3 may prefer repeated elements as target sites or have a preferred target sequence. Four of the six putative ISNgo2 elements (containing irg genes irg2, irg3, irg5, and irg6) are inserted at exactly the same site in the repeated prophage sequence (Fig. 7). One of these elements (irg5) is inserted between the prophage and type dRS3 repeat sequences (42) (Fig. 7B). The irg8-associated element is inserted into a homologue of the putative nitroreductase genes from N. meningitidis strains MC58 and Z2491 (GenBank accession numbers NC_003112 and NC_003116). Based on the N. meningitidis sequence for the uninterrupted nitroreductase gene, the ISNgo2 element is immediately adjacent to a palindromic 34-bp sequence, 5'-TATATTGCGCCCCTTCTTAGGGACGCAATATATA-3', that is present in the nitroreductase gene found in N. gonorrhoeae. Interestingly, the irg4-associated ISNgo2 copy is inserted in the same orientation relative to a sequence that is identical to 27 bp of this 34-bp palindromic sequence which is located within a truncated copy of the IS1016 transposase from N. meningitidis MC58 (Fig. 7A). The putative ISNgo3 elements are found adjacent to identical repeated prophage sequences in both N. gonorrhoeae and N. meningitidis (Fig. 7B and C).
|
|
|---|
At least one of the irg copies is actively transcribed in N. gonorrhoeae, as demonstrated by RT-PCR experiments using a primer that would detect any of the irg2-6 or irg8 gene transcripts. Possible recombinase activities for irg2 and the divergent irg7 gene products were examined in vivo (in E. coli) using plasmid-based systems for Piv-mediated site-specific inversion of the M. lacunata pilin segment and for MooV-mediated excision of IS492. Neither of the two homologues could complement Piv or MooV activities in these systems. While there are a number of possible reasons why complementation was not seen, it is most likely due to differences in DNA binding or cleavage site specificity rather than catalytic activity; Piv and MooV cannot complement for each other, as well (41). There is very little similarity between the Piv-binding sites within the inversion region of M. lacunata and the flanking sequences of the irg genes (data not shown), and the sites for MooV binding and DNA cleavage on IS492 have not yet been determined. To determine if the irg genes code for recombinases involved in antigenic variation at pilE, the eight irg copies were inactivated by using a combination of molecular techniques. To further this analysis, a strain was created which overexpressed irg5. Neither overexpression of irg5 nor inactivation of the eight irg copies affects the frequency of antigenic variation at pilE. We conclude that the hypothesis that the irg genes are recombinases involved in pilin antigenic variation is incorrect. Additionally, the irg genes do not code for recombinases involved in DNA transformation or DNA repair.
Although an involvement in site-specific recombination reactions elsewhere in the gonococcus cannot be ruled out, it is most likely that the irg genes encode transposases for IS110-related insertion sequences in the chromosome of N. gonorrhoeae. In support of this hypothesis, all the irg gene products exhibit significant homology to the transposases of the IS110 family of insertion elements. In addition, a high degree of amino acid sequence identity (>99%) is seen among the irg gene products proposed to act as the transposase for ISNgo2 and among the gene products of irg7, pivNM-1A, and pivNM1-B, which are the predicted transposases for ISNgo3. As would be predicted for an IS element, the sequences flanking the proposed transposases for ISNgo2 are identical in each copy, with the exception of only a 1-bp difference found in the right-end sequence of the irg8-associated copy. Nearly the same level of nucleotide sequence identity is seen for the predicted left and right ends of the three ISNgo3 copies. There are no conserved features that define the ends or target sites of IS110-related elements; in general, there are limited or no terminal inverted repeats, and the element does not create target site duplications. Consistent with belonging to the IS110 family, the putative ISNgo2 element has 2-bp terminal inverted repeats and ISNgo3 has no terminal inverted repeats; in addition, neither element appears to generate target site duplications. While there is no common nucleotide sequence for the termini of all the IS110-related elements, there are subgroups to the family in which the IS elements share nucleotide motifs at or near the termini that may be involved in transposase recognition and cleavage (20). We found some nucleotide sequence identity at or near the ends of IS492 with the proposed termini of the ISNgo2 element. For example, 5'-AGGNNAATC-3' is at the left terminus of ISNgo2 and 19 bp from the left end of IS492. Also, 5'-TAGTNTCT-3' is at the right terminus of ISNgo2 and the juncture of the right end of IS492 and the 5-bp target sequence. The significance of homology between the ends of IS492 and ISNgo2 will not be evident until we have further characterized the features and functionality of the ISNgo2 element.
We cannot conclude that ISNgo2 and ISNgo3 are active IS elements, since we have not observed movement. Carrick et al. (3) reported that the mobility of irg-hybridizing restriction fragments changes during the passage of strain MS11. We have attempted to repeat this observation with both MS11 and FA1090. After 2 to 3 weeks of daily serial passage, neither strain showed a change in the mobilities of irg-hybridizing fragments (data not shown). However, the association of the irg genes with repeated sequences is a feature of the N. meningitidis IS elements IS1016, IS1106, and IS1655, which are preferentially inserted in or near repeat DNA elements (19). In addition, it has recently been reported that the IS110-related elements IS621 and IS1594 target repetitive extragenic palindromic-like repeated sequences in E. coli and Anabaena sp. strain PCC7120, respectively (5).
ISNgo2 and ISNgo3 are the only new IS elements assigned to the N. gonorrhoeae genome since the discovery of the Correia elements (6). Interestingly, four of the six ISNgo2 copies and the one copy of ISNgo3 are inserted adjacent to degenerate prophage sequences that are found repeated on the N. gonorrhoeae and N. meningitidis chromosomes (23). Based on BLASTX analysis, these ISNgo-associated prophage sequences are predicted to encode a filamentous phage structural protein that has a Zot domain (zonular occludens toxin [zot]), a putative surface antigen (tspB), and a filamentous phage replication protein (rstA). Zot was characterized from the filamentous phage CTX
of pathogenic Vibrio cholerae and has the unusual activity of reversibly altering intestinal epithelial tight junctions, allowing macromolecules to pass through mucosal barriers (8, 44). The tspB-like gene exhibits approximately 80% identity with tspB from N. meningitidis which is reported in NCBI (accession number NP_284510) and EMBL (accession number O87783) databases to encode a potent T-cell- and B-cell-stimulating Neisseria-specific antigen. The role, if any, of these prophage-encoded proteins in the pathogenesis of N. gonorrhoeae is unknown.
We are certain from the inactivation and overexpression data that the irg genes have no role in pilin variation, thereby disproving the hypothesis that this gene family has a direct role in this immune evasion mechanism. However, the presence of ISNgo2 and ISNgo3 in the pathogenic Neisseria, but not the commensal Neisseria, and their association with potential prophage-encoded virulence factors suggest that they could be associated with virulence properties. The observation that N. meningitidis has copies of the ISNgo3 element, but not ISNgo2, is suggestive of the evolutionary history, although it is not predictive. Although several scenarios could explain the acquisition of these elements in the pathogenic Neisseria, our preferred model is that ISNgo3 was acquired as part of the bacteriophage genome by a progenitor of the pathogenic Neisseria after it diverged from the commensal species and that ISNgo2 was acquired separately by the progenitor of N. gonorrhoeae after it diverged from N. meningitidis. The ISNgo2 element inserted into preferred target sites associated with the now repeated, degenerate prophage sequence. The possible selective force for duplication of prophage sequences could have been increased expression of the Zot, TspB, or RstA protein, if these phage-associated proteins play a role in gonococcal pathogenesis. Regardless of the evolutionary pathway to the contemporary organisms, these sequences could be used as diagnostic markers for the pathogenic Neisseria.
We acknowledge the Gonococcal Genome Sequencing Project, B. A. Roe, L. Song, S. P. Lin, X. Yuan, S. Clifton, T. Ducey, L. Lewis,and D. W. Dyer at the University of Oklahoma (supported by USPHS/NIH grant number AI38399), particularly for assistance in assigning the irg open reading frames. We thank William Shafer for supplying the neisserial DNA. We appreciate the assistance of Deborah Tobiason, Amanda Ozin, and Jonathan Kim in cloning the irg genes and analyses of irg recombinase activity. We thank Kimberly Kline, Eric Sechman, and Brian Higgins for critical reading of the manuscript.
Present address: Molecular Genetics and Cell Biology, Committee on Microbiology, University of Chicago, Chicago, IL 60637. ![]()
Present address: Quest Diagnostics, San Francisco, CA 94107. ![]()
Present address: Carver College of Medicine, University of Iowa, Iowa City, IA 52242. ![]()
|
|
|---|
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»