ABSTRACT
The genetic structure and evolution of a novel exchangeable meningococcal genomic island was defined for the important human pathogen Neisseria meningitidis. In 125 meningococcal strains tested, one of three unrelated nucleotide sequences, designatedexl (exchangeable locus), was found between a gene required for heme utilization, hemO, andcol, encoding a putative Escherichia colicollagenase homologue. The 5′ boundary of each exlcassette was the stop codon of hemO, whereas the 3′ boundary was delineated by a 33-bp repeat containing neisserial uptake sequences located downstream of col. One of the three alternative exl cassettes contained the meningococcal hemoglobin receptor gene, hmbR (exl3). In other meningococcal strains, hmbR was absent from the genome and was replaced by either a nucleotide sequence containing a novel open reading frame, exl2, or a cassette containing exl3. The proteins encoded byexl2 and exl3 had no significant amino acid homology to HmbR but contained six motifs that are also present in the lipoprotein components of the lactoferrin (LbpB), transferrin (TbpB), and hemoglobin-haptoglobin (HpuA) uptake systems. To determine the evolutionary relationships among meningococci carryinghmbR, exl2, or exl3, isolates representing 92 electrophoretic types were examined.hmbR was found throughout the population structure ofN. meningitidis (genetic distance, >0.425), whereasexl2 and exl3 were found in clonal groups at genetic distances of <0.2. The commensal neisserial species were identified as reservoirs for all of the exl cassettes found in meningococci. The structure of these cassettes and their correlation with clonal groups emphasize the extensive gene pool and frequent horizontal DNA transfer events that contribute to the evolution and virulence of N. meningitidis.
The distinct syndrome of epidemic meningitis and meningococcemia caused by Neisseria meningitidis was first reported in 1805 (6). An interesting aspect of this naturally transformable human pathogen is the rapid evolution of the meningococcal genome through high-frequency recombination events and uptake of DNA (10). Because of this high rate of microevolution, the emergence of new clonal groups responsible for outbreaks of meningococcal disease can be readily detected through genetic typing (2). These approaches have identified the emergence of new clonal groups since the 1960s, the most recent being lineage III, originally found in The Netherlands in 1980 but now detected worldwide (2).
Conceptually, transformation and recombination of foreign DNA into the neisserial genome could also facilitate the incorporation of a new gene(s) into the meningococcal chromosome. DNA fragments containing uptake sequences are preferentially incorporated into the genomes of naturally transformable species such as Haemophilus influenzae(5), N. meningitidis, and N. gonorrhoeae (16). Uptake sequence-mediated DNA uptake has been proposed to be the mechanism by which sodC fromHaemophilus spp. has been incorporated into the meningococcal genome (13). Transformation and recombination are also proposed to be the mechanisms through which the DNA island containing the capsule biosynthesis genes was originally inserted between the conserved genes tex and galEin the genome of N. meningitidis (22; D. S. Stephens et al., unpublished data). Further evidence suggests that this process facilitates the exchange of different capsule biosynthesis loci, resulting in capsule switching (35). Apart from the capsule loci, other meningococcus-specific genetic cassettes found so far include those encoding glutathione peroxidase (19), rotamase (30), and the RTX (repeats in toxin) toxin-like Frp proteins (37). Subtractive hybridization of meningococcal and gonococcal genomes has identified at least eight genetic islands in meningococci (11). However, the extent of the gene pool, the rate of dissemination, the mechanisms of acquisition of new genes, and the association of new genes with virulence in N. meningitidis remain poorly understood.
A comparison of the meningococcal genomes of serogroup A strain Z2491 and serogroup B strain MC58 (21, 36) has revealed that other genetic islands exist within N. meningitidis. The genes encoding the hemoglobin (Hb)-haptoglobin (Hp) receptor,hpuAB, are located on a genetic island that is present in serogroup A strain Z2491 but absent in serogroup B strain MC58. A second phase-variable meningococcal Hb receptor, encoded byhmbR, however, is present in both strains, suggesting that this receptor is necessary for iron uptake from Hb in strain MC58. Conversely, Richardson and Stojiljkovic (28) have found that many clinical meningococcal isolates have lost hmbR but have retained hpuAB. In this report, we have investigated the genetic mechanism that results in the complete loss of thehmbR locus from clinical meningococcal isolates. We have identified an exchangeable DNA island (exl) consisting of multiple cassettes that are capable of being exchanged forhmbR in the meningococcal genome. The commensal neisserial species appear to be the reservoirs for these exl cassettes.
MATERIALS AND METHODS
Bacterial strains and culture conditions.Meningococcal strains were cultured under aerobic conditions with 3.5% CO2 at 37°C on GC agar (Difco) supplemented with 0.4% glucose, 0.01% glutamine, 0.2 mg of cocarboxylase/liter, and 5 mg of Fe(NO3)3/liter. The meningococcal strains used in this study were obtained from several geographic locations over a period of years and span the genetic diversity of N. meningitidis. Invasive isolates collected as part of a prospective population-based surveillance for meningococcal disease in metropolitan Atlanta, Ga., from 1988 to the present (25, 31), serogroup A meningococcal isolates representing each clonal group (I to VIII) (kindly provided by Mark Achtman), and isolates obtained from the collection of the Centers for Disease Control and Prevention, Atlanta, Ga., were examined in this study. Nasopharyngeal isolates were provided from a study of asymptomatic carriers in the metropolitan Atlanta area (Stephens et al., unpublished). Commensal Neisseria spp. have been described previously (18).
Amplification of the exl cassettes and flanking regions.The region surrounding the exl cassettes was PCR amplified from chromosomal preparations of 10 meningococcal isolates: NMB (32), 44/76 (38), 6940 (41), M981 (41), FAM18 (18), M1205 (35), M1843 (35), F8229 (23), GA0929 (Atlanta Surveillance Group, VA Medical Center, Atlanta, Ga.), and 6083 (18). Gonococcal strain FA1090 was used as a positive control. Region I (Fig.1) was PCR amplified from these strains using the primer pair pqi3A (5′-CTCGATATACGTCATACCGTAAGC-3′) and pqi11 (5′-GCGAACACGACTTTGTAGAATGC-3′).
Genetic arrangement of the hmbR locus inN. gonorrhoeae strain FA1090 (accession number AF319529 ) and N. meningitidis strains NMB (accession numbersAF319530 and AF319531 ), 44/76 (accession number AF319527 ), and FAM18 (accession number AF319528 ). Each ORF is represented by an arrow that indicates the direction of transcription. The thin black vertical line indicates the position of a 23-bp sequence which is duplicated inhemO of strain NMB. The thick grey vertical lines indicate the positions of the 33-bp repeats of orfU. Black bars represent a 112-bp conserved nucleotide sequence immediately bordering the 3′ end of hmbR. White bars indicate the locations of neisserial repeat sequences (4). Hatched bars indicate the positions of the 1,280-bp sequence containingorfYX.
Southern blot analysis.Chromosomal DNA from meningococcal strains was prepared using the method of Nath (20). Double-stranded DNA probes were labeled with digoxigenin using a Genius 2 labeling kit (Boehringer Mannheim Biochemicals) according to the manufacturer's instructions. The internal hmbR probe from N. meningitidis strain NMB was amplified with hmbR2 (5′-GAAGCTGCAACCGAAACCACACC-3′) and hmbR3 (5′-GTGGTGCTGTAGCTATACAAACCG-3′), the exl2 probe from N. meningitidis strain GA0929 was amplified with tbp1 (5′-GTCTAGCTCTTTCAGCTTGTAGTTCG-3′) and col12 (5′-GCTTTCGTCAAATCAAGGATTCAAACG-3′), and the exl3probe from N. meningitidis strain 6083 was amplified with Acol23 (5′-GAACATACTGAGCCTGCGATGAATGC-3′) and Acol28 (5′-CATCCCCATACCTGAAGGATCTTTGG-3′). Chromosomal DNA was digested with either Bsp106 or StyI overnight and electrophoresed in a 0.7% agarose gel. Southern transfer onto a nylon membrane was performed as described by Maniatis et al. (15). Hybridization and development of the Southern blots were carried out according to the instructions for the Genius system (Boehringer).
DNA sequencing.Automated sequencing was performed at the Emory DNA Sequencing Core Facility (VA Medical Center, Atlanta, Ga.). The sequence was analyzed using the GAP sequence comparison program of the Genetics Computer Group software package (version 7.3.1 for UNIX) and the Lasergene suite of programs (DNASTAR). All the sequences described in this manuscript have been submitted to GenBank (see below).
Multilocus enzyme electrophoresis.Multilocus enzyme electrophoresis was carried out by established methods (39). Briefly, constitutive cytosolic enzymes were extracted and run in 11% starch gels. The band position for each of 24 cellular enzymes was established by specific staining for each enzyme. Relative enzyme mobilities (alleles) were numbered in order of increasing anodal migration, and each unique set of alleles was defined as an electrophoretic type (ET). Strains representing ET0717, ET0715, ET0712, ET700, and ET708 were included as markers for serogroup A subgroups I, II, III, IV-1, V, VI, VII, and VIII (courtesy of Tanya Popovic, Centers for Disease Control and Prevention) but were not available for exl cassette typing.
Carriage of exl cassettes in the meningococcal population.Colony PCR was performed on each isolate using an anchor primer, pqi10 (5′-GCATTCTACAAAGTCGTGTTCGC-3′), Fpqi10 (5′-GTGAAACCTTCGGTTTGCCGGAAGG-3), or Ypqi10 (5′-GCATTTTACAAAGTCGTGTTGCG-3′), paired with a primer specific for each open reading frame (ORF) located in the exlcassette. ORF-specific primers were hmbR2 and hmbR3 forhmbR; col14 (5′-GTCAAGTTGCGCGAACGGAAATGC-3′) forexl2; Acol26 (5′-CCTGTTTTAGGTTGGTTATCGGCAGC-3′) for exl3; and A1col2 (5′-CATATCACGGTTGTTTTCAGGCTGC-3′) for exl3A. PCR amplification of an internal section of the conserved flanking gene,col, using primers col1 (5′-CGATTCGCTCGAAATCATCCACC-3′) and col2 (5′-GGTGGATGATTTCGAGCGAATCG-3′) was included as a positive control for template quality in every PCR run.
Nucleotide sequence accession numbers.The nucleotide sequences submitted to GenBank (accession numbers in parentheses) are as follows: region I from N. gonorrhoeae strain FA1090 (AF319529 ), region I from N. meningitidis strain NMB (AF319530 ), region II and III in N. meningitidis strain FAM18 (AF319528 ), region II and III from N. meningitidisstrain NMB (AF319531 ), region II and III from N. meningitidis strain 44/76 (AF319527 ); exl2 fromN. meningitidis strain GA0929 (AF319532 ), exl3from N. meningitidis strain 6083 (AF319533 ),exl2C from N. cinerea (AF319536 ),exl3L from N. lactamicaAF319537 ),exl3A from N. meningitidis strain Z5005 (AF319534 ), and exl3 from N. meningitidis strain Z5035 (AF319535 ).
RESULTS
Arrangement of the hmbR (exl1) locus in N. meningitidisWe first studied the organization of the hmbR locus because of the reported linkage of hmbR to the lipopolysaccharide heptosyltransferase encoded by rfaC (33). Using chromosome walking, direct PCR, and sequencing, we found thatrfaC was not linked to hmbR in N. meningitidis or N. gonorrhoeae (Fig. 1). DNA sequencing of the locus located upstream of hmbR an ORF which was subsequently termed hemO, as it encodes the heme oxygenase necessary for the release of iron from imported heme (40). Upstream of hemO lies a homologue ofEscherichia coli pqiA which encodes a paraquat-inducible stress protein (12). The ORF downstream ofhmbR was termed col and was predicted to encode a protein with 63% identity and 75% similarity to a putative collagenase from E. coli.
Region I contains pqiA and hemO and was shown by PCR to be identical in size in eight meningococcal isolates (see Materials and Methods). Region I in gonococcal strain FA1090 produced a larger PCR product which, upon sequencing, was shown to contain insertion element IS1016 and an associated 55-bp deletion located downstream of hemO (accession number AF319529 ). Nevertheless, a comparison of the hemO genes from gonococcal strain FA1090 and serogroup B meningococcal strain NMB (accession number AF319530 ) showed 88.5% nucleotide sequence identity between these two strains, indicating a high degree of conservation. Similarly, region II, which contains col (Fig. 1), was highly conserved at the nucleotide level in meningococcal strains FAM18 (accession number AF319528 ) and NMB (97% identity over 146 bp).
Region III was polymorphic (Fig. 1) and contained four elements that were found in different combinations in the eight meningococcal isolates and the one gonocococcal strain examined. Downstream ofcol in strain FA1090, a 524-bp novel ORF, orfU, flanked by two 33-bp inverted repeats (AATGC CGTCTGAAAACGGTTTCAGACGGCATT) containing neisserial uptake sequences (shown in bold) was found (B. A. Roe, S. Clifton, and D. W. Dyer, Gonococcal Genome Sequencing Project [http://www.genome.ou.edu ]). The organization oforfU resembles that of many insertion sequence elements. However, in meningococcal serogroup B strain NMB, orfU was replaced with a 1,280-bp sequence that contained a 418-bp ORF,orfX, and a 248-bp ORF, orfY (accession numberAF319531 ). OrfX shares 50% identity at the amino acid level with a hypothetical protein encoded by an ORF found in multiple copies in the genomes of Synechocystis spp., whereas OrfY shares 34% amino acid identity with a conserved hypothetical protein fromStreptomyces coelicolor (26). Region III of serogroup B meningococcal strain 44/76 (accession number AF319527 ) and serogroup C strain FAM18 (accession number AF319528 ) possessed hybrid sequences containing orfYX bounded by a neisserial repeat element and the 33-bp inverted repeat and either a partial or a complete orfU gene (Fig. 1). Southern analysis of these strains indicated that both orfYX and orfU were exclusively associated with this genomic location (data not shown). The various combinations and arrangements of the elements in region III suggest that it is a site for recombination. Invariably, region III was delineated by the stop codon of hemO and the 33-bp repeat found downstream of col. Region III will be subsequently referred to as the exl (exchangeable locus) cassette.
Detection and organization of exl2 andexl3 loci in N. meningitidisSouthern analysis of a panel of 10 meningococcal isolates using anhmbR probe revealed that serogroup Y strain GA0929 and serogroup W-135 strain 6083 did not contain hmbR (Fig.2A). The region betweenhemO and col from both isolates was amplified and then sequenced. Region I from strains GA0929 (accession number AF319532 ) and 6083 (accession number AF319533 ) were both 93% identical to the corresponding region in strain NMB (data not shown). Likewise, col from region II (146 bp) was 100% identical in GA0929 and 6083 and 94% identical to the same region in strain NMB. The 1,904-bp exl cassette in strain GA0929, however, diverged from the hmbR loci immediately after the stop codon of hemO (Fig.3) and contained an 809-bp ORF designatedexl2. Southern analysis of the panel of 10 meningococcal isolates using a specific probe containing an internal section ofexl2 showed that none of thehmbR-containing strains carried this ORF elsewhere in their chromosomes (Fig. 2).
(A) Carriage of exl loci in meningococcal genomes is mutually exclusive. Southern blots of 10 meningococcal strains hybridized with hmbR, exl2, andexl3 probes are shown. Lanes 1 through 10 contain genomic DNAs from strains NMB, 44/76, 6940, M981, FAM18, M1205, M1843, F8229, GA0929, and 6083, respectively. Lane L contains labeled size markers, and lane P contains unlabeled probe DNA as a positive control. (B) Genetic arrangement of the exl cassettes inN. meningitidis strainsGA0929 (accession numberAF319532 ) and 6083 (accession number AF319533 ). See the legend to Fig.1 for details.
The 5′ boundary of the exl cassettes corresponds to the stop codon of hemO. Shown is the alignment of the exl cassettes from the 3′ end ofhemO to the start codon of hmbR fromN. meningitidis strain NMB, exl2 fromN. meningitidis strain GA0929, and exl3from N. meningitidis strain 6083. The stop codon of eachhemO sequence and the start codon of eachexl ORF are marked by grey boxes. Nucleotide residues that do not match the consensus sequence are enclosed in black boxes. Gaps are represented by dashes. The position of the 23-bp sequence inhemO (nucleotides 15 to 37) that is duplicated in strain NMB is indicated by arrows.
In strain 6083, region III consisted of a second, novel 3,426-bp sequence containing a 1,368-bp ORF, exl3 (Fig. 2B). Theexl3 ORF contained a homopolymeric cytosine tract located in the 5′ end of the gene and therefore is potentially phase variable. Meningococcal isolates containing hmbR or exl2did not contain exl3 elsewhere in their genomes (Fig. 2A). Although exl2 and exl3 did not have significant nucleotide sequence identity, the translated amino acid sequences from these ORFs were 55% identical (data not shown). Both predicted Exl2 and Exl3 proteins had conserved signal peptides containing a signal peptidase II cleavage site, suggesting that both could be membrane-bound lipoproteins.
A search for amino acid sequence similarity using PSI-BLAST (1) showed that both Exl2 and Exl3 had 47% identity and 61% similarity to a putative lipoprotein (GNA2132) encoded by the ORF NMB2132 in the meningococcal genome of strain MC58 (36). GNA2132, Exl2, and Exl3 all have low sequence identity to the lipoprotein components of the meningococcal transferrin (TbpB), lactoferrin (LbpB), and Hb-Hp (HpuA) receptors (Fig.4). Mazarin et al. (17) have previously shown that the TbpB proteins, which are generally about 700 amino acids long, consist of two domains of approximately 350 amino acid residues each. Both the carboxy- and amino-terminal domains contain the same patterns of conserved motifs, which may be involved in binding the bilobed transferrin substrate (27). Alignments of the carboxy- and amino-terminal sequences of TbpB proteins from various species with the Exl2 and Exl3 sequences revealed that six identifiable motifs, corresponding to those described by Mazarin et al. (17), were found in Exl2 and Exl3 (Fig. 4). Our alignment also reveals that, like TbpB, LbpB (741 amino acids) contains two sets of the six motifs (Fig. 4). Because they are much smaller proteins, Exl2 (269 amino acids), Exl3 (455 amino acids), and GNA2132 (427 amino acids) contain only one set of the six motifs. In this respect, these proteins are similar to HpuA (341 amino acids), a protein that is involved in Hb-Hp utilization.
Alignment of Exl2 and Exl3 with iron binding lipoproteins. The carboxy termini of Exl2, Exl3, and GNA2132 are 45% identical. These three proteins have low identity (below 25%) with the transferrin binding proteins (TbpB) of N. meningitidis(NmTbpB), Moraxella catarrhalis (McTbpB),Actinobacillus pleuropneumoniae (ApTfbB), andHaemophilus influenzae (HiTbpB); the meningococcal lactoferrin receptor (NmLbpB); and the Hb-Hp receptor (NmHpuA). Six motifs were identified in the carboxy termini (A) and amino termini (B) of the NmTbpB and NmLbpB proteins. Exl2, Exl3, GNA2132, and NmHpuA contained only one copy of the six motifs. The consensus sequence for each motif is shown above the alignments. Uppercase letters indicate residues that are conserved in the carboxy- and amino-terminal motifs, whereas lowercase letters indicate residues specific for the motif in each domain. The consensus sequence for the amino-terminal motifs (B) was determined using all nine proteins (data not shown).
Commensal Neisseria species containhmbR, exl2, and exl3Stojiljkovic et al. (34) previously showed that the commensal Neisseria spp., N. polysaccharea and N. perflava, containhmbR. To determine whether exl2 andexl3 were also present in commensalNeisseria spp., Southern hybridization usingexl probes was performed on a representative sample of commensal strains. Representative isolates of N. polysaccharea, N. mucosa, and N. macaca contained a locus which hybridized to thehmbR probe, whereas the exl3 probe hybridized to N. lactamica and N. flavescens. N. cinerea hybridized to the exl2 probe. Three isolates representing the commensal species N. subflava,N. elongata, and N. sicca did not bind to any of the probes. As with meningococci, hmbR,exl2, and exl3 were mutually exclusive, with only one copy of each ORF present in each genome. The linkage ofpqiA to col in N. lactamica, N. flavescens, N. cinerea, and N. polysaccharea was confirmed by PCR. The organization of this region in the N. subflavaand N. sicca isolates that we examined could not be confirmed, even though pqiA and col were present in their genomes (data not shown).
A comparison of the nucleotide sequence of exl2 fromN. meningitidis strain GA0929 with that of its homologue inN. cinerea (accession number AF319536 ) indicated that meningococcal exl2 had been shortened by an internal deletion that removed 331 codons (or one-third of the ORF). The translational frame of meningococcal exl2 was maintained but contained a stop codon at position 203. These features suggest that meningococcal exl2 was a deletion version of the ORF designated exl2C in N. cinerea (data not shown).
The nucleotide sequences of the exl3-like ORF from N. lactamica (exl3L; accession number AF319537 ) andexl3 from N. meningitidis strain 6083 were 63% identical (data not shown). Complete analysis of an exlcassette from the serogroup A meningococcal strain Z5005, which hybridized weakly to the exl3 probe, revealed the presence of a third variant of exl3, termed exl3A(accession number AF319534 ). In contrast to the size variations observed in exl2 variants, all of the exl3variants were of approximately similar sizes. For clarity, we assigned new gene names where the nucleotide sequence of an ORF was less than 90% identical to the closest homologue.
Distributions of hmbR, exl2, andexl3 in the N. meningitidispopulation.A survey of 125 meningococcal isolates representing 92 ETs (genetic distance, >0.425) was conducted to determine the distributions of hmbR, exl2, and exl3or exl3A in the meningococcal population. Multiple isolates with the same ET contained the same exl cassette (data not shown). When categorized by serogroup, isolates of serogroups B, C, and A contained similar proportions of hmbR and exl3or exl3A (approximately 80 and 20%, respectively). In contrast, 55% of serogroup Y strains (11 of 20) examined carriedexl2, 35% (7 of 20) carried exl3 orexl3A, and 10% (2 of 20) contained hmbR. Nongroupable isolates collected from asymptomatic carriers containedhmbR, exl2, and exl3A in a ratio of 6:3:1. From this survey, no correlation of exl type with site of isolation could be found (e.g., blood, cerebrospinal fluid, or nasopharynx) (Fig. 5).
ET dendrogam showing the distribution of differentexl cassettes in the meningococcal population represented by 97 ETs. The carriage of hmbR (stippled boxes), exl2 (black boxes), and exl3 orexl3A (hatched boxes) in 92 ET isolates is shown to the left of the dendrogram. The ET designations appear on they axis in the following manner: serogroup/ET type, where serogroup is represented by the single-letter code (A, B, C, Y, W [W-135], Z, and NG [nongroupable]) and ET type is represented by the four-digit numerical code. All of the representative strains from each ET type are clinical isolates, except for those labeled with an asterisk, indicating isolation from the nasopharynx of an asymptomatic carrier. ET complexes or clonal groups are boxed in grey and are labeled accordingly. The ET clonal groups ET5 and ET17-ET24 and subgroups I, II, III, IV-1, V, VI, VII, and VIII have been described elsewhere (2). The clonal complexes ET501-ET508, A, and B are described in this study.
Analysis of the population by ET clonal group revealed that mostexl2-containing serogroup Y and nongroupable isolates were members of the ET501-ET508 clonal complex (Fig. 5). Only two ETs (ET809 and ET517) containing exl2 were outside the ET501-ET508 complex. Conversely, ET501-ET508 clonal group isolates that did not express the serogroup Y capsule did not contain exl2. These data suggest a genetic linkage between serogroup Y capsule expression and carriage of exl2 in isolates from the ET501-ET508 clonal cluster.
The exl3 and exl3A ORFs were widely distributed (genetic distance, >0.4) across the meningococcal population (Fig. 5). In many instances, isolates containing exl3 andexl3A, however, were found in clonal clusters (genetic distance, <0.2). For example, representative isolates of serogroup A subgroups I, VI, V, and VII contained exl3 andexl3A, and a previously undefined complex, termed B, consisting of five ETs representing serogroup Y isolates, containedexl3A. Complex A consisted of three ETs (ET514, ET384, and ET916) and contained isolates carrying exl3 andexl3A and expressing serogroup Y, B, and W-135 capsules. Thus, the distribution of exl3 and exl3A in the meningococcal population, unlike the distribution of exl2, was not linked to a single capsular serogroup.
Acquisition, dissemination, and microevolution ofexl2 and exl3.To understand the dissemination of the exl cassettes, we analyzed the nucleotide sequences of 10 meningococcal isolates (genetic distance, >0.4) that contained exl3 alleles (allele defined as >90% identity at the nucleotide sequence level). A 290-bp region spanning the last 56 bp of hemO, the intergenic space, and the first 124 bp of the exl3 allele was amplified and sequenced. We hypothesized that if the insertion of theexl3-containing cassettes were the result of a specific integration into this region, then polymorphisms within thehemO nucleotide sequence would be randomly distributed and not linked to polymorphisms in the exl3 sequences. However, in the 10 meningococcal isolates and an N. lactamicaisolate, the sequence polymorphisms in hemO and the intergenic space were linked to the sequence polymorphisms in theexl3, exl3A, and exl3L alleles (Fig.6A). These results indicated that homologous recombination within the conserved flanking gene, hemO, is the mechanism by which theexl cassettes have spread in meningococci.
(A) Polymorphisms in the 3′ end ofhemO and the intergenic space between this position and the exl3 start codon are characteristic of theexl3 alleles in N. meningitidis. Alignment of the nucleotide sequences from position 663 ofhemO (position 1 in this figure) to position 288 ofexl3 in 10 meningococcal isolates and N. lactamica is shown. The numbering of the polymorphic sites within the sequences is shown in vertical format above the sequence for N. lactamica exl3L (accession number AF319537 ). In this scheme, numbers 29 to 31 correspond to the stop codon of hemOand positions 164 to 166 correspond to the exl3 start codon. Nucleotide positions that are identical to those in N. lactamica exl3L are shown as dots, and dashes indicate spaces. Meningococcal isolates (ETs): B431 (375), 6083 (916), B472 (384), M0022 (504), M1978 (525), M1762 (509), Z5005 (723), M2436 (546), M2505 (548), and M2743 (800). (B) Polymorphisms in the first 650 bp ofexl3L from N. meningitidis strain B431 (ET375) and N. lactamica exl3L. Identical nucleotides (dots) and spaces (dashes) are marked according to the reference sequence of exl3 from N. meningitidisstrain 6083 (ET916). The polymorphic sites are numbered from positions 1 to 3 corresponding to the ATG start codon of exl3. Asterisks mark polymorphisms not shared between theexl3L cassettes in N. meningitidis strain B431 and N. lactamica. (C) Alignment of the completeexl3 genes from N. meningitidis strain 6083 (ET916) and N. meningitidis strain Z5035 (ET747). This alignment revealed polymorphisms in the nucleotide sequence that corresponds to the C-terminal domain of the translated protein. In this numbering scheme, positions 1 to 3 correspond to the ATG start codon ofexl3 (see above for details).
To assess microevolution within the exl cassettes, we examined both exl2 and exl3 ORFs for polymorphisms. The exl2 ORF nucleotide sequences of seven isolates, five from the ET501-ET508 clonal complex and one isolate each from unrelated ET809 and ET517 (Fig. 5), were identical. Thus, the acquisition and spread of exl2 in meningococci appears to be recent, since different alleles have not yet been produced through point mutation or homologous recombination. In contrast, the nucleotide sequences of the exl3 and exl3A ORFs showed evidence of microevolution (data not shown). The first 650 bp ofexl3 in strains representing ET916, ET504, ET747, and ET384 were identical, except for a single nucleotide at position 221 in ET384 (data not shown). ET916 and ET384 are closely related (complex A; genetic distance, 0.1) (Fig. 5) and thus would be expected to contain similar exl3 alleles. ET504 and ET747 are not closely related to ET916 or ET384 (genetic distances, 0.25 and 0.34, respectively) and are likely to have acquired the nearly identicalexl3 alleles through horizontal exchange. Examination of theexl3 ORFs from strains representing ET747 (accession numberAF319535 ) and ET916 (accession number AF319533 ) revealed that there were 69 base-pair changes, including transitions, transversions, and small deletions, over a 342-bp section of the C-terminal domain (Fig.6C). Therefore, the meningococcal population contains multiple alleles of exl3, indicating that this region has diverged by microevolution from an ancestral homologue. Similar conclusions were made regarding exl3A from isolates representing ET548, ET800, ET525, ET546, ET509, and ET723 (data not shown).
exl3L in N. meningitidis strain B431 (ET375) contained many of the polymorphisms characteristic of exl3Lin N. lactamica (Fig. 6B), in both the conserved regions (positions 1 to 297 and 450 to 650) and the central region (positions 297 to 450) of the gene. In all, 61 base-pair changes (out of 650 bp) distinguished exl3L from N. meningitidis strain B431 from that of N. lactamica. This result suggests that acquisition of the exl cassettes by meningococci has occurred multiple times.
DISCUSSION
The ability of N. meningitidis and N. gonorrhoeae to establish infections is in part due to their ability to scavenge essential nutrients, such as iron, from human tissues. Indeed, it has been recently shown that the ability to utilize transferrin is essential for gonococci to establish an infection in the human challenge model (3). Unlike many other pathogens that use siderophores to shuttle iron, meningococci and gonococci express at least four outer membrane protein receptors that are specific for transferrin (tbpAB), lactoferrin (lbpAB), Hb (hmbR), and Hb-Hp complexes (hpuAB) (7). Once heme has been transported into the cytoplasm, iron is released from the heme by a heme oxygenase-like enzyme, encoded by hemO, which is found upstream of hmbR (40).
However, during the investigation of the genetic region containinghemO, we found that many meningococcal isolates did not contain hmbR. In these isolates, hmbR was absent from the chromosome and had been replaced by two distinct DNA islands, which contain ORFs exl2 and exl3. The predicted Exl2 and Exl3 proteins had significant sequence identity with GNA2132, a highly antigenic TbpB-like protein that was recently identified as a potential vaccine candidate (24). In a broader context, Exl2, Exl3, and GNA2132 contain six motifs that are also present in TbpB, LbpB, and HpuA (Fig. 4). TbpB contains two sets of these motifs, one each in the N-terminal and C-terminal domains, which are hypothesized to be involved in binding to the bilobed substrate human transferrin (27). The function of these motifs, however, in LbpB and HpuA, which differ from TbpB in their substrate specificities, is currently unknown. The family of proteins consisting of Exl2, Exl3, and GNA2132 may represent iron binding proteins with an as-yet-undefined substrate specificity.
A comparison of the organization of the exl cassettes in meningococci and commensal Neisseria spp. suggests that each cassette arose from a specific insertion event. This concept is supported by the finding that each exl cassette is delineated by clearly defined 5′ (hemO stop codon) and 3′ (uptake sequence downstream of col) boundaries. In this respect, each cassette resembles an integron, but no other features, such as an integration site (att) or a gene encoding an integrase, could be identified (9). Since orfYXand orfU have some of the features of insertion sequence elements and are located at the 3′ boundary of each exlcassette, they may have been involved in the initial acquisition of theexl cassettes (8). Correlation of the nucleotide polymorphisms outside the 5′ boundary of the exlcassette with base-pair changes found in the cassette, however, indicates that once the cassette was inserted, it was disseminated by homologous recombination within the conserved flanking neisserial genes.
The prevalence of hmbR within the meningococcal population and in N. gonorrhoeae suggests that hmbR may have been present in a common ancestor. Although hmbR has a high degree of primary sequence conservation, alleles of this gene have been previously identified (34), indicating thathmbR has undergone microevolution within meningococci. Examination of the exl cassettes containing exl3and exl3A within the meningococcal population also revealed that these genes have undergone microevolution, and evidence of transfer of identical exl3 and exl3A alleles to distantly related strains indicates that horizontal transfer has occurred in the meningococcal population. The finding thatexl3 and exl3A are mostly confined to a few clonal groups suggests, however, that these alleles have only recently been acquired by meningococci and are slowly being dispersed by homologous recombination outside these organisms.
The detection of only one allele of the exl2-containing cassette suggests that it has only recently been acquired by meningococci. Further, the exl2-containing DNA cassette is largely restricted to strains within the ET501-ET508 clonal complex. Interestingly, the ET501-ET508 clonal complex has a relatively short history, as it has emerged over the last decade as a major causative agent of serogroup Y disease in the United States (29). The presence of the exl2-containing DNA cassette in meningococci offers the unique opportunity to track the rates of microevolution and dissemination of a novel gene within the meningococcal population.
In this report, we have defined an exchangeable DNA island that has been inserted into a single site in the meningococcal chromosome and that consists of at least three distinct cassettes. One of these cassettes encodes the Hb receptor, HmbR, while the other two cassettes encode unique proteins, Exl3 or Exl3A and Exl2. It is currently unclear why so many different strategies have been used by N. meningitidis to change the patterns of expression of the Hb receptors HpuAB and HmbR. Apart from phase-variable expression, the genes of both receptors can be completely replaced or removed from the genome. It is possible that complete removal of HmbR has been developed as a way of replacing a relatively antigenically stable outer membrane protein (34) in the meningococcal population. For many other outer membrane proteins, such as TbpB (14), a high level of antigenic variation is necessary for maintenance in the meningococcal population.
ACKNOWLEDGMENTS
This work was supported by NIH grant AI-33517–08 from the National Institute of Allergy and Infectious Diseases of the National Institutes of Health.
We thank John K. Davies and Harry Sakellaris for critical reading of the manuscript. We are also grateful to Lane Pucko for assistance in preparing the manuscript.
Notes
Editor: T. R. Kozel
FOOTNOTES
- Received 6 September 2000.
- Returned for modification 5 October 2000.
- Accepted 27 November 2000.
- Copyright © 2001 American Society for Microbiology
