Skip to main content

Comparative analysis of two Neisseria gonorrhoeae genome sequences reveals evidence of mobilization of Correia Repeat Enclosed Elements and their role in regulation



The Correia Repeat Enclosed Element (CREE) of the Neisseria spp., with its inverted repeat and conserved core structure, can generate a promoter sequence at either or both ends, can bind IHF, and can bind RNase III and either be cleaved by it or protected by it. As such, the presence of this element can directly control the expression of adjacent genes. Previous work has shown differences in regulation of gene expression between neisserial strains and species due to the presence of a CREE. These interruptions perhaps remove the expression of CREE-associated genes from ancestral neisserial regulatory networks.


Analysis of the chromosomal locations of the CREE in Neisseria gonorrhoeae strain FA1090 and N. gonorrhoeae strain NCCP11945 has revealed that most of the over 120 copies of the element are conserved in location between these genome sequences. However, there are some notable exceptions, including differences in the presence and sequence of CREE 5' of copies of the opacity protein gene opa, differences in the potential to bind IHF, and differences in the potential to be cleaved by RNase III.


The presence of CREE insertions in one strain relative to the other, CREE within a prophage region, and CREE disrupting coding sequences, provide strong evidence of mobility of this element in N. gonorrhoeae. Due to the previously demonstrated role of these elements in altering transcriptional control and the findings from comparing the two gonococcal genome sequences, it is suggested that regulatory differences orchestrated by CREE contribute to the differences between strains and also between the closely related yet clinically distinct species N. gonorrhoeae, Neisseria meningitidis, and Neisseria lactamica.


The genome sequences of the Neisseria spp. contain 100 or more copies of a repetitive sequence that has demonstrated roles in gene regulation and is involved in chromosomal changes such as gene inactivation and rearrangements. The Correia Repeat (CR) of 26 bp was first described in 1986 [1], and has only been found in Neisseria spp., both pathogenic (Neisseria gonorrhoeae, the gonococcus, and Neisseria meningitidis, the meningococcus) and commensal (Neisseria lactamica) [2, 3]. Often it is found as an inverted repeat with a characteristic core [4], the Correia Repeat Enclosed Element (CREE) [5], which is most commonly 153–157 bp or 104–108 bp (Figure 1). The CREE has features of an IS element including the presence of the inverted repeat CR sequences [4] and duplication of the target sequence [6, 7]. In addition, there are vast differences in their placement in the neisserial genome sequences and some coding sequences have been disrupted by CREE [3, 5, 79]. This suggests they are, or have been, mobile genetic elements. Unlike IS elements that encode their own mobilizing transposase, the CREE have never been demonstrated to encode proteins and the mobilizing mechanisms for CR and CREE distribution in the chromosome have not been determined.

Figure 1

CREE sequence structures. Comparative illustration of the sequence structures of different lengths of CREE. Sequences in the key are based on consensus sequences of all aligned N. gonorrhoeae CREE. Three different length types of CREE are present in the gonococcal genome sequences, however the shortest (69–73 bp) have two different structure types. An IHF binding site is only possible in the longest of the CREE (153–157 bp).

Two functional σ70 promoters have been described as being generated by CREE insertion 5' of genes [10, 11] (Figure 2). The first of these, the Black promoter, is at the end of the CREE closest to the gene [10]. The Black promoter is comprised of a -35 and a partial -10 element within the CREE. Insertion in some locations completes the -10 promoter element (Figure 2). Functional Black promoters are present 5' of uvrB [10], drg [12], lst [13], and mtrCDE [14]. The second σ70 promoter, the Snyder promoter, was identified at the end of the CREE farthest from the gene. In this case the -35 element comes from the native sequence and the -10 is contained completely within the first 6 bases of the CREE (Figure 2). A functional Snyder promoter was first identified 5' of dcw gene dcaC in N. lactamica [11].

Figure 2

CREE sequence features. For illustration purposes a consensus sequence of the longer CREE is shown, here 156 bp. The ways in which the two different CREE associated promoters are generated from native sequence (blue line) and CREE sequence (red line) are shown for the Snyder promoter (left) [11] and the Black promoter (right) [10]. Both of these would drive transcription from left to right. Given favourable native sequence it might be possible to generate both promoters on the opposite strand as well, transcribing right to left. The IHF-binding site within the CREE is shown (orange), as are the inverted repeats that potentially form mRNA hairpin substrates for RNase III (green).

The endoribonuclease RNase III is involved in the processing of stable RNAs such as rRNA and some mRNA transcripts [15]. The stem-loop of the CREE generated by the inverted repeats (Figures 1 and 2) has been determined to be a binding site for RNase III when the CREE sequence is present within the mRNA [16]. The binding of RNase III regulates gene expression post-transcriptionally either through cleavage or protection from cleavage by RNase III, depending on the CR inverted repeat sequence [14, 16, 17].

The CREE may have other functions as well. For example, some of the longer 153–157 bp CREE contain an IHF binding site (Figures 1 and 2). This may be involved in end synapsis during element transposition [6] or participate in the modulation of regulation of associated genes [14]. CREE may also be hotspots for genomic recombination and rearrangement [5, 7, 18, 19]. Additionally, CREE are associated with gene loss through deletion [18, 20] and gene disruption [8, 21, 22].

Despite all of these known functions, we still do not know whether CREE are currently mobile within these species or whether their mobilization occurred in the distant past in the evolution of these bacteria and the mobilizing element has now been lost. None are present in known or suspected regions of horizontally transferred DNA [6, 9, 23, 24], which suggests that mobility has been lost. However, diversity in the distribution in the N. meningitidis sequenced chromosomes [5, 8, 9, 25] and rearrangements between the genome sequences strongly suggest that CREE are mobile in meningococci.

The state of the CREE in N. gonorrhoeae could previously only be investigated in the context of its comparison with N. meningitidis due to the availability of only one gonococcal genome sequence [GenBank: RefSeq NC_002946]. In the meningococcus over 250 CREE have been identified [5, 8, 9, 25], yet less than half this number have been reported for gonococcal strain FA1090 [5, 6]. Is this the result of a species difference in CREE copy number between the gonococcus and the meningococcus? Or, has the relatively longer laboratory propagation of strain FA1090 resulted in the deletion of genetic material including CREE sequences? If the latter is the case then comparison with a more recently acquired and less passaged gonococcal isolate should provide evidence of CREE-mediated deletion events. Likewise, variations in CREE locations between two gonococcal genome sequences might suggest that mobilization has occurred in N. gonorrhoeae since its split with the other neisserial species and that such CREE movement might still be occurring. The recent publication of the complete genome sequence of N. gonorrhoeae strain NCCP11945, an isolate from a 2002 vaginal smear [26], has opened the way to better understand this human sexually-transmitted pathogen through comparative genome analysis. Here we investigate the CREE of N. gonorrhoeae strains FA1090 and NCCP11945.

Results and discussion

Identification of CREE in N. gonorrhoeae strain FA1090 and N. gonorrhoeae strain NCCP11945

A new analysis of the N. gonorrhoeae strain FA1090 genome sequence was conducted, identifying a total of 123 CREE and 9 single CR sequences (see Additional file 1). A total of 131 CREE and 6 single CR sequences were identified in the genome sequence of N. gonorrhoeae strain NCCP11945 (see Additional file 2). The CREE are fairly evenly distributed in the genomes (Figure 3). For both genome sequences, 17 of the CREE contain a perfect IHF binding site based on the published CREE IHF binding site sequence [6], although these 17 are different between the two strains (see Additional file 3). The IHF-binding site-containing CREE are 154 to 157 bp, except for a 143 bp variant that has an 11 bp deletion in its CR in strain NCCP11945 (Figure 3 and Additional file 3). The potential to be cleaved by RNase III is dependent upon the symmetry of the CR ends of the CREE [17]. In strain FA1090, 92 have symmetrical CR ends and 100 CREE have symmetrical ends in strain NCCP11945.

Figure 3

CREE distribution in the chromosomes. Distribution of the CREE (x-axis) and their lengths (y-axis) in the genome sequences of N. gonorrhoeae strain FA1090 (A) and N. gonorrhoeae strain NCCP11945 (B). Length variants are shown with larger spots and are colour coded: green – 97 bp CREE due to a 9 bp deletion; purple – 110 bp CREE due to a 5 bp duplication; blue – 143 bp CREE due to an 11 bp deletion; red – 169 bp CREE due to a 64 bp duplication in a 105 bp CREE.

Most CREE locations are nearly identical in N. gonorrhoeae strains FA1090 and NCCP11945

Of the 123 CREE in strain FA1090, 120 are at the same locations in strain NCCP11945 (see Additional file 4). Of these, 83 (69.16%) are nearly identical with 3 or fewer nucleotide differences in sequence and no differences in length. Some of these are 5' of genes that are integral to gonococcal biology (Table 1).

Table 1 CREE 5' of important gonococcal genes

The CREE have become associated with various systems within the gonococcus including those involved in LPS biosynthesis, pilus expression, DNA repair, iron acquisition, adhesion, competence, and pathogen-specific gene regulation. It is possible that variations in the presence of CREE associated with such genes may account for differences in strain behaviour and phenotypes that differ between the gonococcus and the meningococcus. For example, in N. gonorrhoeae the MtrCDE efflux pump system, involved in antibiotic resistance [27] and in vitro survival [28], is regulated by the repressor MtrR [27] and the activator MtrA [29]. In N. meningitidis, a CREE has inserted within the regulatory region for mtrCDE and removed the expression of these genes from control by MtrR and MtrA [14]. In this case, the CREE is present between the native promoter [29] and the genes, thus the CREE sequence is part of the mRNA transcript and its cleavage by RNase III has been demonstrated [14]. This CREE also contains an IHF-binding site and its presence in the CREE was shown to have a negative impact on transcription levels [14].

CREE location differences between the gonococcal strains

Of the CREE identified in the genome sequences of N. gonorrhoeae strain FA1090 and N. gonorrhoeae strain NCCP11945, thirteen are only present in one of the strains (see Additional files 1, 2, and 4). Four of these are within larger regions of difference. Fragments of CREE are found at three other sites. The remaining six are clear insertions of the element, often at a TA target site as has been described [7]. However, two insertions appear to have occurred at CA sites (strain NCCP11945 positions 1,577,266 and 2,185,647) and indeed some CREE have terminal CA sequences rather than TA (see Additional files 1 and 2). These insertions of CREE in just one of the gonococcal strains support the mobility of the element within the gonococcus and suggest that the mobilization mechanism for the CREE could not have been lost before the speciation event that generated N. gonorrhoeae and N. meningitidis.

Also identified was one CREE in strain FA1090 that corresponded to two different CREE in strain NCCP11945 (Figure 4). The associated chromosomal rearrangement has not created fragments of CREE in either strain.

Figure 4

CREE and chromosomal rearrangement. Alignment of the N. gonorrhoeae strain FA1090 (top) and strain NCCP11945 (bottom) genome sequences using progressive Mauve [43]. Several chromosomal rearrangements and inversions of large regions of sequence can be seen represented here by coloured blocks. Those below the mid-line in strain NCCP11945 are inverted relative to the strain FA1090 sequence. The arrows indicate the left and right ends of a single CREE in strain FA1090 and their alignment with two separate CREE at a chromosomal rearrangement point relative to strain NCCP11945.

N. gonorrhoeae have less than half the CREE copies of N. meningitidis

It is clear that the gonococci differ from their meningococcal relatives in the number of CREE copies sustained in their genome sequences, with N. meningitidis genomes having over 250 [8, 9, 25]. The commensal N. lactamica, a close relative of the pathogenic neisserial species, has been reported to have three times fewer CREE than N. meningitidis [2]. A quick search suggests that there are approximately 100 CREE in the N. lactamica ST-640 genome sequence (data not presented). Therefore, the commensal and the gonococcus share similar CREE copy numbers.

Coupled with the functions of the CREE, this suggests that N. gonorrhoeae and N. lactamica retain more of the ancestral Neisseria regulatory networks than N. meningitidis, where many of these may now be under the transcriptional control of CREE. This is certainly true for the efflux pump system encoded by mtrCDE, where both the gonococcal repression [27] and activation [29] systems have been subverted through CREE insertion in the meningococcus [14]. It has been difficult to assign a species-specific gene set to N. gonorrhoeae, N. meningitidis, and N. lactamica, largely due to horizontal exchange between these naturally competent species [30]. The small number of genetic islands identified that were thought to be unique to one species are either present in both pathogenic species or are not present in all strains of the species in which they were originally found. For example, the Gonococcal Genetic Island is not present in all gonococci [31] and was found in strains of N. meningitidis [24]. The meningococcal capsule is absent from some N. meningitidis strains including some that have caused invasive disease [32, 33]. Other Islands of Horizontal Transfer found in N. meningitidis are strain-specific [9, 25]. Added to this, only six 'virulence' genes present in all pathogenic Neisseria genome sequences were absent from the non-pathogen N. lactamica [34]. Evidence increasingly supports the idea that regulation rather than gene complement differentiates the species.

The Correia Repeat of 26 bp is a fragment of CREE sequence

The single CR sequences found in the two gonococcal genome sequences were often associated with a partial core sequence (6 out of 9 and 4 out of 6). In addition, the vast majority of the CR are part of CREE (see Additional files 1 and 2). This suggests that the 26 bp repeats identified by Correia et al. [1] are fragments of CREE, having lost the remainder of the sequence through deletion events. Cases of partial CR (<26 bp) are also found in strain FA1090 that correspond to whole CREE 5' of NGK_2270 (mafA) and between NGK_2168 and NGK_2169 in strain NCCP11945 (see Additional file 3). The CREE itself should therefore be seen as the functional unit, rather than the chance occurrence of an inverted pair of CR. This is especially evident in light of the conservation of the core sequences (Figure 1).

Characteristics of the CREE

The gonococcal CREE are evenly distributed in the chromosomes (Figure 3). Most are of the 154–156 bp or 105–107 bp types, however a small number are shorter than these at 69–73 bp. Sequence feature characteristics of the different lengths of CREE illustrated in Figure 1 show a general conservation of CREE structure based on length. Although some length variants shown in Figure 3 were found, in each case these are modifications of one of the basic CREE structures through deletions and tandem repeat duplications. All CREE lengths can have the variations in the inverted repeat end sequences (Figure 1; AACAAAAA or AAATTTAAA) that have been described previously [6]. Symmetry of these inverted repeats generates the potential for RNase III cleavage [17]. The potential for an IHF-binding site is present only in the longer of the CREE (Figure 1). CREE core sequence directionality has a CCGGTACGG end and a TCAGGACAA end (Figure 1). The 71 to 73 bp CREE are the exception to this structure, having two CCGGTACGG ends, yet these still retain a directionality to their core sequence (Figure 1).

CREE 5' of genes and roles in regulation

A total of 18 potential Black promoters and 15 potential Snyder promoters were found associated with N. gonorrhoeae strain FA1090 CREE, with 7 of these CREE containing both (see Additional file 1). While few potential Black and Snyder promoters were identified, there are 76 CREE in strain FA1090 positioned 5' of a gene. This suggests that the presence of the element in this location influences the regulation or expression of these gene(s). It is likely that RNase III has a role in post-transcriptional regulation of these genes. Of these 76 CREE, 57 have symmetrical inverted repeat ends and are therefore potential substrates for RNase III cleavage if the CREE sequence is part of the transcript. Asymmetrical sites are also influenced by RNase III through binding, therefore RNase III is thought to influence the longevity of all mRNAs containing CREE sequence. Seven CREE differ between the strains in their end symmetry, with four of these located 5' of annotated CDSs (NGO0407, NGO0452, NGO1221, and NGO1347). Likewise, nine CREE 5' of CDSs contain IHF-binding sites that may influence the binding and action of other proteins, including RNA polymerase and RNase III.

CREE associations with copies of the opacity protein gene opa

The Opa proteins of the Neisseria spp. are a family of phase variably expressed, antigenically variable outer membrane proteins involved in attachment and invasion of host cells. Multiple copies of the opa gene are found in neisserial genomes and variations in their sequences can mediate different interactions with different host cell receptors [35, 36]. There are 11 opa genes in strain FA1090 [37] and in strain NCCP11945 (Table 2). Most of the strain FA1090 opa genes are unannotated and in several cases the strain NCCP11945 annotation has not identified the initiation codon due to the frame-shift generated by the phase variable CTTCT tract within the gene [38]. The majority of the opa genes are associated with a 5' CREE, which in some cases contains an IHF-binding site sequence consensus (Table 2).

Table 2 The opa-associated CREE

In all cases, there appears to be a σ70 promoter between the CREE and opa, which would mean that the CREE sequence is not part of the mRNA transcript and that it is not therefore targeted by RNase III. There are sequence differences in the IHF-binding sites between the strains, but whether the sites that differ from the published IHF-binding site sequence are still able to bind IHF is not known. The CREE 5' of NGK_1799 in strain NCCP11945 is not present in strain FA1090, where a 17 bp fragment of the end of the CREE remains 5' of NGO1513. The CREE 5' of NGK_1847 (unannotated in strain FA1090) is only 105 bp and therefore does not carry the IHF-binding site (Figure 1).

Indeed, there are no CREE at all present 5' of the opa copies NGK_0096, NGK_0847, NGK_2170, and NGK_2410 (all unannotated in strain FA1090) (Table 2). The regions upstream of these genes are otherwise similar, which would suggest that the CREE was mobile after the gene duplication events that generated 11 copies of opa, three to four times the number of opa genes found in N. meningitidis [8, 9, 25]. Indeed all of the copies of opa found in N. meningitidis strain Z2491 have a CREE 98–99 bp 5' (NMA1676, NMA1890, NMA2043), with similarly placed CREE for three out of four opa genes in meningococcal strain MC58 (NMB0442, NMB1465, NMB1636, but not NMB0926) [8, 9, 34]. When the phase variable tract is ON in opa the influences of the various upstream sequences, including the CREE, may have further regulatory effects on expression of this important adhesin, perhaps through binding of IHF.

CREE as insertional inactivators of genes

Of the 5 previously reported genes disrupted by CREE insertion [8, 21, 22], 4 are associated with CREE in these two gonococcal genome sequences (Table 3). In the case of NMA2121, this CDS is not present in strain FA1090 as it is part of a Minimal Mobile Element [22, 39, 40], which facilitates horizontal exchange of gene cassettes between genome sequences.

Table 3 Disruption of annotated coding regions by CREE

Three additional coding sequence disruptions by CREE were identified in strains FA1090 and NCCP11945 (Table 3). In two of these, the disruption is shared by both strains, while the third has no CREE in this region in strain NCCP11945. The fact that the CREE has inserted into identifiable coding regions supports the hypothesis that this element is, or was, mobile. That disruption of some CDSs is seen only in the gonococcus supports mobilization of the CREE within N. gonorrhoeae since its split from N. meningitidis.

CREE within potential regions of horizontal transfer

It has been proposed that CREE are not found in regions of horizontal transfer and that their absence can be taken as additional evidence for the horizontal origin of a region, given their otherwise even distribution in the chromosome [6, 9, 23, 24]. The CREE sequences 38 bp 5' of NGO1004, 213 bp 5' of NGO1006, and 21 bp 3' of NGO1020, present in both gonococcal strains, are within a 13 kb region between rpoD and tRNA-Arg containing 17 CDSs annotated as "putative phage associated proteins". This demonstrates that the CREE has been mobile at some time since the acquisition of this prophage sequence in the Neisseria. While it might be possible to use the absence of CREE as additional evidence indicating foreign sequence origin, as for the meningococcal IHTs [9] and the Gonococcal Genetic Island [24], the presence of CREE does not preclude the possibility of horizontal transfer of the region.


The CREE are less numerous in the gonococcal genome sequences (~120 copies) than in the meningococcal genome sequences (~250 copies). Many of the gonococcal CREE (71 of 113 in strain FA1090) are 5' of genes and in this position are likely to be involved in the regulation of genes through generation of promoter sequences, binding of IHF, and mRNA stability control by RNase III. CREE are associated with virulence- and host survival-associated determinants such as Opa, HpuA, and MafB, as well as regulators of factors important in pathogenesis. Differences in CREE insertions between the two strains and the presence of CREE within a prophage region, upstream of multiple copies of opa, and within coding sequences provide strong evidence of mobility of this element in the gonococcus and therefore since the speciation event that differentiated the gonococcus from the meningococcus. The regulatory influence of the CREE and the copy number differences between N. gonorrhoeae, N. meningitidis, and N. lactamica may contribute to the different behaviours of these pathogens through differences in their regulatory networks.


CREE identification in the genome sequences

The genome sequences of N. gonorrhoeae strain FA1090 [GenBank: RefSeq NC_002946] and N. gonorrhoeae strain NCCP11945 [GenBank: RefSeq NC_011035] were searched for CR sequences using the tools within x BASE[41]. Each genome sequence was searched using 'fuzznuc' for three different patterns based on the conserved CR sequences identified from previous studies of this element [5, 6, 8, 25]: TATAG [CT]GGATTAACAAAAATCAGGAC; TATAG [CT]GGATTAAATTTAAACCGGTAC; TATAG [CT]GGATTAACAAAAACCGGTAC; and TATAG [CT]GGATTAAATTTAAATCAGGAC. Square brackets will find either of the letters within them. In each case, up to 3 mismatches were allowed. Increasing to 4 or 5 mismatches did not identify any additional CREE. Pattern finding has been successfully used previously to identify CREE [3, 5]. Each identified location was then investigated manually to determine if the site contained a lone CR or a CREE. Each location was catalogued and the length and sequence of the CREE was determined (see Additional files 1 and 2). As a confirmation, searches with the above sequences were also conducted using BLASTN, as has been done previously [57], against the gonococcal genome sequences and also against N. lactamica ST-640 (data not presented). IHF sites were found by searching within the CREE sequences for the sequence CTTCAGCACCTTA [6] (see Additional file 3). RNase III cleavage potential was determined based on CR sequence and symmetry between the ends of the CREE.

Genome sequence comparative analysis

Using x BASE[41], WebACT[42], and progressive Mauve[43], the genome sequence data from gonococcal strains FA1090 and NCCP11945 were comparatively aligned and visualized. For each CREE location the presence of a corresponding CREE in the other genome sequence, if present, was catalogued (see Additional file 4).

CREE sequence alignment

For each site where CREE were found in both the gonococcal genome sequences, the sequences of the CREE were aligned using EMBOSS needle [44] and the calculated percent identity was recorded (see Additional file 4). All of the identified CREE sequences were aligned together using ClustalW2 [45] from EMBL-EBI.

Other nomenclature used

Full length inverted repeats

Correia repeats; Correia elements (CE); Correia; Correia full (~150 bp); Correia internal deletion (~105 bp); Neisseria miniature insertion sequences (nemis); neisserial repetitive sequences; pseudo-transposable neisserial small repetitive elements (SRE); N. gonorrhoeae inverted repeats; Correia sequences; 152-bp repetitive sequence.

Single repeats

Correia repeats; Correia repeat unit; 26-bp repeat; 26-bp neisserial repeat (NR).


  1. 1.

    Correia FF, Inouye S, Inouye M: A 26-base-pair repetitivesequence specific for Neisseria gonorrhoeae and Neisseria meningitidis genomic DNA. J Bacteriol. 1986, 167 (3): 1009-1015.

    PubMed Central  CAS  PubMed  Google Scholar 

  2. 2.

    De Gregorio E, Abrescia C, Carlomagno MS, Di Nocera PP: Asymmetrical distribution of Neisseria miniature insertion sequence DNA repeats among pathogenic and nonpathogenic Neisseria strains. Infect Immun. 2003, 71 (7): 4217-4221. 10.1128/IAI.71.7.4217-4221.2003.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  3. 3.

    Qvarnstrom Y, Swedberg G: Variations in gene organization and DNA uptake signal sequence in the folP region between commensal and pathogenic Neisseria species. BMC Microbiol. 2006, 6: 11-10.1186/1471-2180-6-11.

    PubMed Central  Article  PubMed  Google Scholar 

  4. 4.

    Correia FF, Inouye S, Inouye M: A family of small repeated elements with some transposon-like properties in the genome of Neisseria gonorrhoeae. J Biol Chem. 1988, 263 (25): 12194-12198.

    CAS  PubMed  Google Scholar 

  5. 5.

    Liu SV, Saunders NJ, Jeffries A, Rest RF: Genome analysis and strain comparison of correia repeats and correia repeat-enclosed elements in pathogenic Neisseria. J Bacteriol. 2002, 184 (22): 6163-6173. 10.1128/JB.184.22.6163-6173.2002.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  6. 6.

    Buisine N, Tang CM, Chalmers R: Transposon-like Correia elements: structure, distribution and genetic exchange between pathogenic Neisseria sp. FEBS Lett. 2002, 522 (1–3): 52-58. 10.1016/S0014-5793(02)02882-X.

    Article  CAS  PubMed  Google Scholar 

  7. 7.

    Mazzone M, De Gregorio E, Lavitola A, Pagliarulo C, Alifano P, DiNocera PP: Whole-genome organization and functional properties of miniature DNA insertion sequences conserved in pathogenic Neisseriae. Gene. 2001, 278 (1–2): 211-222. 10.1016/S0378-1119(01)00725-9.

    Article  CAS  PubMed  Google Scholar 

  8. 8.

    Parkhill J, Achtman M, James KD, Bentley SD, Churcher C, Klee SR, Morelli G, Basham D, Brown D, Chillingworth T, et al: Complete DNA sequence of a serogroup A strain of Neisseria meningitidis Z2491. Nature. 2000, 404 (6777): 502-506. 10.1038/35006655.

    Article  CAS  PubMed  Google Scholar 

  9. 9.

    Tettelin H, Saunders NJ, Heidelberg J, Jeffries AC, Nelson KE, Eisen JA, Ketchum KA, Hood DW, Peden JF, Dodson RJ, et al: Complete genome sequence of Neisseria meningitidis serogroup B strain MC58. Science. 2000, 287 (5459): 1809-1815. 10.1126/science.287.5459.1809.

    Article  CAS  PubMed  Google Scholar 

  10. 10.

    Black CG, Fyfe JA, Davies JK: A promoter associated with the neisserial repeat can be used to transcribe the uvrB gene from Neisseria gonorrhoeae. J Bacteriol. 1995, 177 (8): 1952-1958.

    PubMed Central  CAS  PubMed  Google Scholar 

  11. 11.

    Snyder LA, Shafer WM, Saunders NJ: Divergence and transcriptional analysis of the division cell wall (dcw) gene cluster in Neisseria spp. Mol Microbiol. 2003, 47 (2): 431-442. 10.1046/j.1365-2958.2003.03204.x.

    Article  CAS  PubMed  Google Scholar 

  12. 12.

    Cantalupo G, Bucci C, Salvatore P, Pagliarulo C, Roberti V, Lavitola A, Bruni CB, Alifano P: Evolution and function of the neisserial dam-replacing gene. FEBS Lett. 2001, 495 (3): 178-183. 10.1016/S0014-5793(01)02388-2.

    Article  CAS  PubMed  Google Scholar 

  13. 13.

    Packiam M, Shell DM, Liu SV, Liu YB, McGee DJ, Srivastava R, Seal S, Rest RF: Differential expression and transcriptional analysis of the alpha-2,3-sialyltransferase gene in pathogenic Neisseria spp. Infect Immun. 2006, 74 (5): 2637-2650. 10.1128/IAI.74.5.2637-2650.2006.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  14. 14.

    Rouquette-Loughlin CE, Balthazar JT, Hill SA, Shafer WM: Modulation of the mtrCDE-encoded efflux pump gene complex of Neisseria meningitidis due to a Correia element insertion sequence. Mol Microbiol. 2004, 54 (3): 731-741. 10.1111/j.1365-2958.2004.04299.x.

    Article  CAS  PubMed  Google Scholar 

  15. 15.

    Nikolaev N, Silengo L, Schlessinger D: A role for ribonuclease 3 in processing of ribosomal ribonucleic acid and messenger ribonucleic acid precursors in Escherichia coli. J Biol Chem. 1973, 248 (22): 7967-7969.

    CAS  PubMed  Google Scholar 

  16. 16.

    De Gregorio E, Abrescia C, Carlomagno MS, Di Nocera PP: The abundant class of nemis repeats provides RNA substrates for ribonuclease III in Neisseriae. Biochim Biophys Acta. 2002, 1576 (1–2): 39-44.

    Article  CAS  PubMed  Google Scholar 

  17. 17.

    De Gregorio E, Abrescia C, Carlomagno MS, Di Nocera PP: Ribonuclease III-mediated processing of specific Neisseria meningitidis mRNAs. Biochem J. 2003, 374 (Pt 3): 799-805. 10.1042/BJ20030533.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  18. 18.

    Claus H, Elias J, Meinhardt C, Frosch M, Vogel U: Deletion of the meningococcal fetA gene used for antigen sequence typing of invasive and commensal isolates from Germany: frequencies and mechanisms. J Clin Microbiol. 2007, 45 (9): 2960-2964. 10.1128/JCM.00696-07.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  19. 19.

    Davidsen T, Tonjum T: Meningococcal genome dynamics. Nat Rev Microbiol. 2006, 4 (1): 11-22. 10.1038/nrmicro1324.

    Article  CAS  PubMed  Google Scholar 

  20. 20.

    Ende van der A, Hopman CT, Dankert J: Deletion of porA by recombination between clusters of repetitive extragenic palindromic sequences in Neisseria meningitidis. Infect Immun. 1999, 67 (6): 2928-2934.

    PubMed Central  PubMed  Google Scholar 

  21. 21.

    Klee SR, Nassif X, Kusecek B, Merker P, Beretti JL, Achtman M, Tinsley CR: Molecular and biological analysis of eight genetic islands that distinguish Neisseria meningitidis from the closely related pathogen Neisseria gonorrhoeae. Infect Immun. 2000, 68 (4): 2082-2095. 10.1128/IAI.68.4.2082-2095.2000.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  22. 22.

    Snyder LA, Davies JK, Saunders NJ: Microarray genomotyping of key experimental strains of Neisseria gonorrhoeae reveals gene complement diversity and five new neisserial genes associated with Minimal Mobile Elements. BMC Genomics. 2004, 5 (1): 23-10.1186/1471-2164-5-23.

    PubMed Central  Article  PubMed  Google Scholar 

  23. 23.

    Snyder LA, Davies JK, Ryan CS, Saunders NJ: Comparativeoverview of the genomic and genetic differences between the pathogenic Neisseria strains and species. Plasmid. 2005, 54 (3): 191-218. 10.1016/j.plasmid.2005.04.005.

    Article  CAS  PubMed  Google Scholar 

  24. 24.

    Snyder LA, Jarvis SA, Saunders NJ: Complete and variant forms of the 'gonococcal genetic island' in Neisseria meningitidis. Microbiology. 2005, 151 (Pt 12): 4005-4013. 10.1099/mic.0.27925-0.

    Article  CAS  PubMed  Google Scholar 

  25. 25.

    Bentley SD, Vernikos GS, Snyder LA, Churcher C, Arrowsmith C, Chillingworth T, Cronin A, Davis PH, Holroyd NE, Jagels K, et al: Meningococcal genetic variation mechanisms viewed through comparative analysis of serogroup C strain FAM18. PLoS Genet. 2007, 3 (2): e23-10.1371/journal.pgen.0030023.

    PubMed Central  Article  PubMed  Google Scholar 

  26. 26.

    Chung GT, Yoo JS, Oh HB, Lee YS, Cha SH, Kim SJ, Yoo CK: Complete genome sequence of Neisseria gonorrhoeae NCCP11945. J Bacteriol. 2008, 190 (17): 6035-6036. 10.1128/JB.00566-08.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  27. 27.

    Hagman KE, Pan W, Spratt BG, Balthazar JT, Judd RC, Shafer WM: Resistance of Neisseria gonorrhoeae to antimicrobial hydrophobic agents is modulated by the mtrRCDE efflux system. Microbiology. 1995, 141 (Pt 3): 611-622.

    Article  CAS  PubMed  Google Scholar 

  28. 28.

    Jerse AE, Sharma ND, Simms AN, Crow ET, Snyder LA, Shafer WM: A gonococcal efflux pump system enhances bacterial survival in a female mouse model of genital tract infection. Infect Immun. 2003, 71 (10): 5576-5582. 10.1128/IAI.71.10.5576-5582.2003.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  29. 29.

    Rouquette C, Harmon JB, Shafer WM: Induction of the mtrCDE-encoded efflux pump system of Neisseria gonorrhoeae requires MtrA, an AraC-like protein. Mol Microbiol. 1999, 33 (3): 651-658. 10.1046/j.1365-2958.1999.01517.x.

    Article  CAS  PubMed  Google Scholar 

  30. 30.

    Maiden MC: Population genomics: diversity and virulence in the Neisseria. Curr Opin Microbiol. 2008

    Google Scholar 

  31. 31.

    Dillard JP, Seifert HS: A variable genetic island specific for Neisseria gonorrhoeae is involved in providing DNA for natural transformation and is found more often in disseminated infection isolates. Mol Microbiol. 2001, 41 (1): 263-277. 10.1046/j.1365-2958.2001.02520.x.

    Article  CAS  PubMed  Google Scholar 

  32. 32.

    Findlow H, Vogel U, Mueller JE, Curry A, Njanpop-Lafourcade BM, Claus H, Gray SJ, Yaro S, Traore Y, Sangare L, et al: Three cases of invasive meningococcal disease caused by a capsule null locus strain circulating among healthy carriers in Burkina Faso. J Infect Dis. 2007, 195 (7): 1071-1077. 10.1086/512084.

    Article  PubMed  Google Scholar 

  33. 33.

    Hoang LM, Thomas E, Tyler S, Pollard AJ, Stephens G, Gustafson L, McNabb A, Pocock I, Tsang R, Tan R: Rapid and fatal meningococcal disease due to a strain of Neisseria meningitidis containing the capsule null locus. Clin Infect Dis. 2005, 40 (5): e38-42. 10.1086/427875.

    Article  PubMed  Google Scholar 

  34. 34.

    Snyder LA, Saunders NJ: The majority of genes in the pathogenic Neisseria species are present in non-pathogenic Neisseria lactamica, including those designated as 'virulence genes'. BMC Genomics. 2006, 7: 128-10.1186/1471-2164-7-128.

    PubMed Central  Article  PubMed  Google Scholar 

  35. 35.

    Hauck CR, Meyer TF: 'Small' talk: Opa proteins as mediators of Neisseria-host-cell communication. Curr Opin Microbiol. 2003, 6 (1): 43-49. 10.1016/S1369-5274(03)00004-3.

    Article  CAS  PubMed  Google Scholar 

  36. 36.

    Gray-Owen SD: Neisserial Opa proteins: impact on colonization, dissemination and immunity. Scand J Infect Dis. 2003, 35 (9): 614-618. 10.1080/00365540310016042.

    Article  CAS  PubMed  Google Scholar 

  37. 37.

    Snyder LA, Butcher SA, Saunders NJ: Comparative whole-genome analyses reveal over 100 putative phase-variable genes in the pathogenic Neisseria spp. Microbiology. 2001, 147 (Pt 8): 2321-2332.

    Article  CAS  PubMed  Google Scholar 

  38. 38.

    Muralidharan K, Stern A, Meyer TF: The control mechanism of opacity protein expression in the pathogenic Neisseriae. Antonie Van Leeuwenhoek. 1987, 53 (6): 435-440. 10.1007/BF00415499.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Saunders NJ, Snyder LA: The minimal mobile element. Microbiology. 2002, 148 (Pt 12): 3756-3760.

    Article  CAS  PubMed  Google Scholar 

  40. 40.

    Snyder LA, McGowan S, Rogers M, Duro E, O'Farrell E, Saunders NJ: The repertoire of minimal mobile elements in the Neisseria species and evidence that these are involved in horizontal gene transfer in other bacteria. Mol Biol Evol. 2007, 24 (12): 2802-2815. 10.1093/molbev/msm215.

    Article  CAS  PubMed  Google Scholar 

  41. 41.

    Chaudhuri RR, Loman NJ, Snyder LA, Bailey CM, Stekel DJ, Pallen MJ: x BASE2: a comprehensive resource for comparative bacterial genomics. Nucleic Acids Res. 2008, D543-546. 36 Database

  42. 42.

    Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J: ACT: the Artemis Comparison Tool. Bioinformatics. 2005, 21 (16): 3422-3423. 10.1093/bioinformatics/bti553.

    Article  CAS  PubMed  Google Scholar 

  43. 43.

    Darling AC, Mau B, Blattner FR, Perna NT: Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14 (7): 1394-1403. 10.1101/gr.2289704.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  44. 44.

    Needleman SB, Wunsch CD: A general method applicable to the search for similarities in the amino acid sequence of two proteins. J Mol Biol. 1970, 48 (3): 443-453. 10.1016/0022-2836(70)90057-4.

    Article  CAS  PubMed  Google Scholar 

  45. 45.

    Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.

    Article  CAS  PubMed  Google Scholar 

Download references


LS was supported by a BBSRC grant (BBE111791) awarded to MP and D.J. Stekel. This research made use of x BASE: a comprehensive resource for comparative bacterial genomics [41], which is also supported by this BBSRC grant.

Author information



Corresponding author

Correspondence to Lori AS Snyder.

Additional information

Authors' contributions

LS conceived of the study, carried out the comparative genome analyses, data collection, and drafted the manuscript. JC and MP assisted in the interpretation of the data and suggested improvements to the presentation of the results and their discussion in the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional File 1: CREE in N. gonorrhoeae strain FA1090. This spreadsheet contains the locations, sequence, and features of all of the CREE identified in the genome sequence of N. gonorrhoeae strain FA1090. (XLS 60 KB)

Additional File 2: CREE in N. gonorrhoeae strain NCCP11945. This spreadsheet contains the locations, sequence, and features of all of the CREE identified in the genome sequence of N. gonorrhoeae strain NCCP11945. (XLS 46 KB)

Additional File 3: Comparison of long CREE and IHF-binding sites between N. gonorrhoeae strains FA1090 and NCCP11945. This spreadsheet contains comparative data on the long CREE in N. gonorrhoeae strain FA1090 and N. gonorrhoeae strain NCCP11945 and the presence of an IHF-binding site that corresponds to the sequence published [6]. (XLS 48 KB)

Additional File 4: Comparison of CREE locations between N. gonorrhoeae strains FA1090 and NCCP11945. This spreadsheet contains the CREE that are in common between N. gonorrhoeae strain FA1090 and N. gonorrhoeae strain NCCP11945. (XLS 76 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Snyder, L.A., Cole, J.A. & Pallen, M.J. Comparative analysis of two Neisseria gonorrhoeae genome sequences reveals evidence of mobilization of Correia Repeat Enclosed Elements and their role in regulation. BMC Genomics 10, 70 (2009).

Download citation


  • Inverted Repeat
  • Neisseria Gonorrhoeae
  • Neisseria Meningitidis
  • Strain FA1090
  • Gonorrhoeae Strain