The Use of High-Throughput DNA Sequencing in the Investigation of Antigenic Variation: Application to Neisseria Species

Antigenic variation occurs in a broad range of species. This process resembles gene conversion in that variant DNA is unidirectionally transferred from partial gene copies (or silent loci) into an expression locus. Previous studies of antigenic variation have involved the amplification and sequencing of individual genes from hundreds of colonies. Using the pilE gene from Neisseria gonorrhoeae we have demonstrated that it is possible to use PCR amplification, followed by high-throughput DNA sequencing and a novel assembly process, to detect individual antigenic variation events. The ability to detect these events was much greater than has previously been possible. In N. gonorrhoeae most silent loci contain multiple partial gene copies. Here we show that there is a bias towards using the copy at the 3′ end of the silent loci (copy 1) as the donor sequence. The pilE gene of N. gonorrhoeae and some strains of Neisseria meningitidis encode class I pilin, but strains of N. meningitidis from clonal complexes 8 and 11 encode a class II pilin. We have confirmed that the class II pili of meningococcal strain FAM18 (clonal complex 11) are non-variable, and this is also true for the class II pili of strain NMB from clonal complex 8. In addition when a gene encoding class I pilin was moved into the meningococcal strain NMB background there was no evidence of antigenic variation. Finally we investigated several members of the opa gene family of N. gonorrhoeae, where it has been suggested that limited variation occurs. Variation was detected in the opaK gene that is located close to pilE, but not at the opaJ gene located elsewhere on the genome. The approach described here promises to dramatically improve studies of the extent and nature of antigenic variation systems in a variety of species.     

Frequency and Rate of Pilin Antigenic Variation of Neisseria meningitidis

The rates of pilin antigenic variation (Av) of two strains of Neisseria meningitidis were determined using an unbiased DNA sequencing assay. Strain MC58 underwent pilin Av at a rate similar to that of N. gonorrhoeae strain MS11 but lower than that of N. gonorrhoeae strain FA1090. Pilin Av was undetectable in strain FAM18.Neisseria meningitidis is a Gram-negative diplococcus that colonizes the nasopharynx of approximately 5 to 10% of the population and is usually nonpathogenic but can occasionally enter the bloodstream to cause septicemia and can eventually spread to the meninges, causing meningitis (15). Approximately 500,000 cases of meningococcal meningitis occur every year, with nearly 10% resulting in fatality (2).Type IV pili (TFP) are long filamentous structures protruding from the bacterial surface and are required for adherence of N. meningitidis to host cells (7). As with the TFP of the closely related pathogen Neisseria gonorrhoeae, the pili are able to undergo antigenic variation (Av). In N. gonorrhoeae, pilin Av occurs as a result of recombination between one of the multiple silent pilS copies and the expressed pilin gene (pilE). The pilS copies share significant regions of homology with pilE yet lack a promoter or ribosome-binding site and the initial 5′ coding sequence. Pilin Av relies on RecA and the RecF-like recombination pathway to catalyze gene conversion, resulting in an altered pilE sequence, carrying part of the pilS donor, and the original unaltered pilS sequence (8, 9).While the frequency of pilin Av has been measured in N. gonorrhoeae (5, 10, 12), this process has never been quantified in N. meningitidis. Two sequenced strains were picked to measure pilin Av: serogroup B strain MC58 (sequence type [ST-32] complex), isolated from an invasive infection (14), and serogroup C strain FAM18 (ST-11 complex), which was isolated from a patient with septicemia (1). In both strains, the native recA gene was replaced with the very highly conserved N. gonorrhoeae recA6 construct, which allows regulation of expression with IPTG (isopropyl-β-d-thiogalactopyranoside) (17). recA6 strains are RecA⁺ when grown with IPTG but are RecA⁻ when grown without IPTG (17). These phenotypes were confirmed by measuring the UV sensitivities and DNA transformation competence levels of both strains with or without IPTG, and both strains were shown to be piliated by transmission electron microscopy (data not shown). Bacteria were grown at 37°C with 5% CO₂ on gonococcal medium base (GCB; Difco) plus Kellogg supplements I and II (11).The pilin Av sequencing assay was performed as described previously (5, 10, 12) with slight modifications. Briefly, FAM18 and MC58 were grown on solid GCB with 1 mM IPTG, allowing for the expression of RecA, for 22 h and 12.5 h, respectively, which was estimated to produce 20 generations. For FAM18, little or no pilin Av was expected since the G-quartet-forming sequence required for pilin Av is degenerate in this strain (3). Therefore, two random progenitor colonies were picked from IPTG-enriched medium and passaged on GCB without IPTG. Between 91 and 94 colonies were isolated from each FAM18 progenitor, and the sequence of the pilE gene was determined. For MC58, seven random progenitor colonies were picked and passaged on GCB without IPTG. Between 28 and 47 progeny colonies arising from each of the seven progenitors were isolated, and the pilE gene sequence was determined. In both MC58 and FAM18, the progeny colonies were passaged on GCB two times to ensure colony clonality. A single colony from each sample was isolated, and the pilE gene was PCR amplified as described previously (13).The primers used for amplification of MC58 pilE were McPilRBS (5′-GCATTTCCTTTCCAATTAGGAG) and MC58SP3A (5′-TTCCGTACGGATAGCTTCGTC). The primers used for amplification of FAM18 pilE were FAMFOR-2 (5′-ATTACGGGTTTACGTTTGCGG) and FAMREV-2 (5′-ACGCACCTACGCCTCACCCTAC). The DNA sequence was determined for each sample (SeqWright, Houston, TX, and the Genomics Core at Northwestern University) and analyzed (MacVector; Symantec Corp.). Colonies that showed pilE sequence changes were reanalyzed to confirm the Av event.The pilin Av frequency was determined for the progeny of each progenitor by dividing the total number of detected pilin Av events by the number of progeny of each set, resulting in two values for FAM18 and seven values for MC58 (Table ​(Table1).1). The pilin Av rate was determined by dividing the pilin Av frequency by the number of generations for each sample grown in the presence of IPTG, as determined by a colony assay at the time of harvest. After growth for the same time period, the total numbers of generations for MC58 and FAM18 grown under RecA induction were approximately 19 and 23, respectively.

TABLE 1.

Frequencies and rates of pilin Av in MC58 and FAM18^a

Molecular models accounting for the gene conversion reactions mediating gonococcal pilin antigenic variation

The pilus antigenic variation (Av) system of Neisseria gonorrhoeae is one of several high-frequency variation systems that utilize gene conversion to switch between numerous forms of an antigen on the cell surface. We have tested three predictions of the first models that explain the movement of DNA during pilin Av: (i) Av requires two recombinations at short regions of identity, (ii) circular intermediates exist that carry pilE/pilS hybrid loci and (iii) these pilE/pilS hybrid loci target the pilS sequences to a recipient pilE gene. We confirm that normal pilin Av utilizes recombination at very short regions of DNA sequence identity and that these recombination events can occur independent of homologous recombination functions. We have isolated covalently closed circular DNA molecules carrying hybrid pilin loci, but propose that an alternative hybrid molecule is the intermediate of pilin Av. Our most striking finding is that transformation of isolated pilE/pilS hybrid loci targets the pilS sequences of the hybrid to a recipient pilE at frequencies much higher than normal recombination frequencies. These results show that the different steps of a model that explains pilin Av can be separately tested to support the validity of these novel models that account for the high-frequency gene conversions that mediate pilin Av.     

The frequency and rate of pilin antigenic variation in Neisseria gonorrhoeae

The pilin antigenic variation (Av) system of Neisseria gonorrhoeae (Gc) mediates unidirectional DNA recombination from silent gene copies into the pilin expression locus. A DNA sequencing assay was developed to accurately measure pilin Av in a population of Gc strain FA1090 arising from a defined pilin progenitor under non-selective culture conditions. This assay employs a piliated parental Gc variant with a recA allele whose promoter is replaced by lac-regulatory elements, allowing for controlled induction of pilin Av. From this assay, the frequency of pilin Av was measured as 0.13 recombination events per cell, with a corresponding rate of pilin Av of 4x10(-3) events per cell per generation. Most pilin variants retained the parental piliation phenotype, providing the first comprehensive analysis of piliated variants arising from a piliated progenitor. Sequence analysis of pilin variants revealed that a subset of possible recombination events predominated, which differed between piliated and non-piliated progeny. Pilin Av exhibits the highest reported frequency of any pathogenic gene conversion system and can account for the extensive pilin variation detected during human infection.     

Sequence Features Contributing to Chromosomal Rearrangements in Neisseria gonorrhoeae

Through whole genome sequence alignments, breakpoints in chromosomal synteny can be identified and the sequence features associated with these determined. Alignments of the genome sequences of Neisseria gonorrhoeae strain FA1090, N. gonorrhoeae strain NCCP11945, and N. gonorrhoeae strain TCDC-NG08107 reveal chromosomal rearrangements that have occurred. Based on these alignments and dot plot pair-wise comparisons, the overall chromosomal arrangement of strain NCCP11945 and TCDC-NG08107 are very similar, with no large inversions or translocations. The insertion of the Gonococcal Genetic Island in strain NCCP11945 is the most prominent distinguishing feature differentiating these strains. When strain NCCP11945 is compared to strain FA1090, however, 14 breakpoints in chromosomal synteny are identified between these gonococcal strains. The majority of these, 11 of 14, are associated with a prophage, IS elements, or IS-like repeat enclosed elements which appear to have played a role in the rearrangements observed. Additional rearrangements of small regions of the genome are associated with pilin genes. Evidence presented here suggests that the rearrangements of blocks of sequence are mediated by activation of prophage and associated IS elements and reintegration elsewhere in the genome or by homologous recombination between IS-like elements that have generated inversions.     

The pilE gene product of Pseudomonas aeruginosa, required for pilus biogenesis, shares amino acid sequence identity with the N-termini of type 4 prepilin proteins

A new locus required for type 4 pilus biogenesis by Pseudomonas aeruginosa has been identified. A pilE mutant, designated MJ-6, was broadly resistant to pili-specific phages and unable to translocate across solid surfaces by the pilus-dependent mechanism of twitching motility (Twt⁻). Immunoblot analysis demonstrated that MJ-6 was devoid of pili (Pil⁻) but was unaffected in the production of unassembled pilin pools. Genetic studies aimed at localizing the pilE mutation on the P. aeruginosa PAO chromosome demonstrated a strong co-linkage between MJ-6 phage resistance and the proB marker located at 71 min. Cloning of the pilE gene was facilitated by the isolation and identification of a proB⁺-containing plasmid from a PAO1 cosmid library. Upon introduction of the PA01 proB⁺ cosmid clone into MJ-6, sensitivity to pili-specific phage, twitching motility and pilus production were restored. The nucleotide sequence of a 1 kb Eco RV-Clal fragment containing the pilE region revealed a single complete open reading frame with characteristic P. aeruginosa codon bias. PilE, a protein with a molecular weight of 15278, showed significant sequence identity to the pilin precursors of P. aeruginosa and to other type 4 prepilin proteins. The region of highest homology was localized to the N-terminal 40 amino acid residues. The putative PilE N-terminus contained a seven-residue basic leader sequence followed by a consensus cleavage site for prepilin pep-tidase and a largely hydrophobic region which contained tyrosine residues (Tyr-24 and Tyr-27) previously implicated in maintaining pilin subunit-subunit interactions. The requirement of PilE in pilus biogenesis was confirmed by demonstrating that chromosomal pilE insertion mutants were pilus- and twitching-motility deficient.     

Analysis of the Piv recombinase-related gene family of Neisseria gonorrhoeae

Neisseria gonorrhoeae (the gonococcus) is an obligate human pathogen and the causative agent of the disease gonorrhea. The gonococcal pilus undergoes antigenic variation through high-frequency recombination events between unexpressed pilS silent copies and the pilin expression locus pilE. The machinery involved in pilin antigenic variation identified to date is composed primarily of genes involved in homologous recombination. However, a number of characteristics of antigenic variation suggest that one or more recombinases, in addition to the homologous recombination machinery, may be involved in mediating sequence changes at pilE. Previous work has identified several genes in the gonococcus with significant identity to the pilin inversion gene (piv) from Moraxella species and transposases of the IS110 family of insertion elements. These genes were candidates for a recombinase system involved in pilin antigenic variation. We have named these genes irg for invertase-related gene family. In this work, we characterize these genes and demonstrate that the irg genes do not complement for Moraxella lacunata Piv invertase or IS492 MooV transposase activities. Moreover, by inactivation of all eight gene copies and overexpression of one gene copy, we conclusively show that these recombinases are not involved in gonococcal pilin variation, DNA transformation, or DNA repair. We propose that the irg genes encode transposases for two different IS110-related elements given the names ISNgo2 and ISNgo3. ISNgo2 is located at multiple loci on the chromosome of N. gonorrhoeae, and ISNgo3 is found in single and duplicate copies in the N. gonorrhoeae and Neisseria meningitidis genomes, respectively.     

Role for the RecBCD recombination pathway for pilE gene variation in repair-proficient Neisseria gonorrhoeae

The role of the RecBCD recombination pathway in PilE antigenic variation in Neisseria gonorrhoeae is contentious and appears to be strain dependent. In this study, N. gonorrhoeae strain MS11 recB mutants were assessed for recombination/repair. MS11 recB mutants were found to be highly susceptible to DNA treatments that caused double-chain breaks and were severely impaired for growth; recB growth suppressor mutants arose at high frequencies. When the recombination/repair capacity of strain MS11 was compared to that of strains FA1090 and P9, innate differences were observed between the strains, with FA1090 and P9 rec⁺ bacteria presenting pronounced recombination/repair defects. Consequently, MS11 recB mutants present a more robust phenotype than the other strains that were tested. In addition, MS11 recB mutants are also shown to be defective for pilE/pilS recombination. Moreover, pilE/pilS recombination is shown to proceed with gonococci that carry inverted pilE loci. Consequently, a novel RecBCD-mediated double-chain-break repair model for PilE antigenic variation is proposed.     

Analysis in Neisseria meningitidis and other Neisseria species of genes homologous to the FKBP immunophilin family

The immunophilin family of FK506-binding proteins (FKBPs), involved in eukaryotic protein folding and cell regulation, have recently been found to have prokaryotic homologues. Genes with sequences homologous to those encoding human FKBPs were examined in Neisseria species. An FKBP DNA sequence was present, as shown by the polymerase chain reaction and Southern blotting experiments, in the chromosome of Neisseria meningitidis (14 strains) and in all 11 different commensal Neisseria spp. studied, but was not found in Neisseria gonorrhoeae (11 strains tested) or in Moraxella catarrhalis. The nucleotide and predicted protein sequences of the FKBP-encoding domain from five of the meningococcal strains were highly conserved (e.g. ≥97% homologous). The meningococcal nucleotide sequence was ≥93% homologous and the consensus meningococcal protein sequence was ≥97% homologous to FKBP sequences found in seven different commensal Neisseria spp. The meningococcal nucleotide and predicted protein sequences were ≥59% homologous to the conserved C-terminus of the human FKBP gene family. The FKBP nucleotide sequence was present as a single copy in the chromosome of commensal Neisseria spp. and in most strains of N. meningitidis. The FKBP gene was linked to the silent pilin locus, pilS, in class II-piliated meningococcal strains. In meningococcal strains expressing class I pili, the FKBP gene was linked to one of several pilS loci but not the pilE locus present in these strains. FKBP genes found in commensal Neisseria spp. were not linked to known pilin loci.