Phylogenetic Evidence for Frequent Positive Selection and Recombination in the Meningococcal Surface Antigen PorB

Rachel Urwin*, Edward C. Holmes{dagger}, Andrew J. Fox{ddagger}, Jeremy P. Derrick§ and Martin C. J. Maiden*

*The Peter Medawar Building for Pathogen Research and Department of Zoology, University of Oxford;
{dagger}Department of Zoology, University of Oxford;
{ddagger}Meningococcus Reference Unit, Public Health Laboratory, Withington Hospital, Manchester;
§Department of Biomolecular Sciences, University of Manchester Institute of Science and Technology


    Abstract
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
Previous estimates of rates of synonymous (dS) and nonsynonymous (dN) substitution among Neisseria meningitidis gene sequences suggested that the surface loops of the variable outer membrane protein PorB were under only weak selection pressure from the host immune response. These findings were consistent with studies indicating that PorB variants were not always protective in immunological and microbiological assays and questioned the suitability of this protein as a vaccine component. PorB, which is expressed at high levels on the surface of the meningococcus, has been implicated in mechanisms of pathogenesis and has also been used as a typing target in epidemiological investigations. In this work, using more precise estimates of selection pressures and recombination rates, we have shown that some residues in the surface loops of PorB are under very strong positive selection, as great as that observed in human immunodeficiency virus-1 surface glycoproteins, whereas amino acids within the loops and the membrane-spanning regions of the protein are under purifying selection, presumably because of structural constraints. Congruence tests showed that recombination occurred at a rate that was not sufficient to erase all phylogenetic similarity and did not greatly bias selection analysis. Homology models of PorB structure indicated that many strongly selected sites encoded residues that were predicted to be exposed to host immune responses, implying that this protein is under strong immune selection and requires further examination as a potential vaccine candidate. These data show that phylogenetic inference can be used to complement immunological and biochemical data in the choice of vaccine candidates.


    Introduction
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
The generation of antigenic diversity has evolved as a strategy for evading immune attack in a wide range of pathogenic and commensal organisms (Deitsch, Moxon, and Wellems 1997Citation ). It is effective against both natural and artificially induced immunity and represents a major obstacle to the development of vaccines against pathogens as diverse as Plasmodium falciparum, human immunodeficiency virus (HIV), and Neisseria meningitidis, the meningococcus. In an era when many new vaccine candidates are being identified by genomic techniques, sequence data of the antigen genes obtained from population samples of pathogens can be analyzed by phylogenetic and biochemical modeling techniques to provide a picture of the evolutionary processes acting on these sequences and hence a preliminary evaluation of their vaccine potential.

The meningococcus is an appropriate model system to evaluate this approach because it is a pathogen of global significance which is genetically and antigenically diverse and for which no comprehensive vaccine exists (Pollard and Frasch 2001Citation ). Further, large genetically defined isolate collections have been assembled and models of the population biology of this organism are available (Caugant et al. 1987Citation ; Maiden et al. 1998Citation ). Amongst the candidate vaccine components proposed are the variable outer membrane proteins, the trimeric porins, which act as pores for the passage of solutes into the cell (Tommassen et al. 1990Citation ). These molecules are targeted by the host immune response and have been used in meningococcal typing schemes (Bjune et al. 1991aCitation , 1991bCitation ; Sierra et al. 1991Citation ; van der Ley and Poolman 1992Citation ; van der Ley, van der Biezen, and Poolman 1995Citation ). Unlike most other Neisseria species, the meningococcus expresses two porins, PorA and PorB. Expression of PorA is regulated at transcription and exhibits three levels depending on the length of the polyguanidine stretch in the promoter region of the porA gene (van der Ende et al. 1995Citation ), whereas there is no evidence to suggest that PorB proteins are subject to phase variation. Mutant meningococcal strains that lack PorB do not grow well (Tommassen et al. 1990Citation ), suggesting that PorB has a function essential for growth. PorB is also capable of translocating vectorially into the membranes of mammalian cells (Blake and Gotschlich 1987Citation ), and of binding ATP and GTP, which down regulates pore size and alters voltage dependence and ion selectivity (Rudel et al. 1996Citation ). These functional and structural characteristics are thought to influence the early stages of neutrophil activation and therefore implicate the PorB protein in meningococcal pathogenesis (Rudel et al. 1996Citation ).

A PorB topology model has been constructed on the basis of nucleotide sequence data (Maiden et al. 1991Citation ; van der Ley et al. 1991Citation ) and, more recently, the structural similarity between the Neisseria porins and the Escherichia coli porins OmpF and PhoE has been exploited to generate a three-dimensional homology model for Neisseria porins (Derrick et al. 1999Citation ). These models predicted eight surface exposed "loops" interspersed with highly conserved outer membrane-spanning sequences that formed a "ß-barrel" (Kleffel et al. 1985Citation ). The antigenically variable epitopes targeted in the host immune response (Saukkonen et al. 1989Citation ) were proposed to reside in the surface-exposed loops (McGuinness et al. 1990Citation ; Maiden et al. 1991Citation ). Serological and molecular characterization of the meningococcal porins have been used for epidemiological analyses of meningococcal carriage and disease (Frasch, Zollinger, and Poolman 1985Citation ; Poolman et al. 1986Citation ; Maiden et al. 1991Citation ), although the variability of these proteins means that they are not always reliable epidemiological markers (Achtman 1995Citation ; Urwin et al. 1998aCitation , 1998bCitation ).

The meningococcus possesses one of two PorB protein classes, PorB2 or PorB3, which are encoded by alternate allele classes present at the porB locus. Phylogenetic analyses show that the porB3 gene is most closely related to one of the gonococcal porin genes, porB1a (Smith, Maynard Smith, and Spratt 1995Citation ) and that these gene sequences form a clade together with gonococcal porB1b, N. lactamica por, and N. polysaccharea por gene sequences (Derrick et al. 1999Citation ). The meningococcal porB2 gene shares sequence similarity both with members of this clade (specifically, the sequence encoding the putative ATP and GTP binding site) and with the porin genes of most human commensal and animal Neisseria species, suggesting that porB2 may have arisen because of interspecies recombination (Derrick et al. 1999Citation ). Indeed, the incongruence between phylogenetic trees drawn for individual loop-encoding regions of porB3 genes (Bash et al. 1995Citation ) suggests that inter- and intraspecies recombination is a mechanism that increases genetic variation among porB alleles.

Previous comparisons of the rates and distribution of synonymous (dS) and nonsynonymous (dN) substitutions among meningococcal porin genes have concluded that, unlike the gonococcal PIA and PIB porins which were under positive selection in the surface loop regions, meningococcal porB genes were subject to only weak positive selection and purifying selection (Smith, Maynard Smith, and Spratt 1995Citation ). This observation supports the hypothesis that PorB is a less important vaccine constituent than the PorA protein. This conclusion was also drawn from studies of the bactericidal activity of monoclonal (Saukkonen et al. 1987Citation ) and polyclonal (van der Ley and Poolman 1992Citation ) antibodies in mice and immunological results from human vaccine trials (Rosenqvist et al. 1995Citation ; Perkins et al. 1998Citation ). In addition, some PorB3 proteins expressed on the surfaces of live meningococci have been reported to be poorly accessible for antibody binding (Michaelsen et al. 2001Citation ). Consequently, PorB has been deliberately excluded from some vaccine formulations (van der Ley, van der Biezen, and Poolman 1995Citation ). But the previous selection analysis was conducted on very small data sets, and dS and dN were estimated as mean values across whole or partial gene sequences, making it possible that strongly selected sites were missed in this broadscale comparison (Smith, Maynard Smith, and Spratt 1995Citation ). More recently, phylogenetic analyses carried out on a larger set of gonococcal porin gene sequences concluded that there were differences in the evolution of PIA and PIB homology groups, with positive selection driving evolution of the PIA proteins and both positive and purifying selection acting on PIB protein sequences (Posada et al. 2000Citation ).

In this work we undertook rigorous maximum likelihood analyses of selection pressures acting on a large set of meningococcal PorB sequences. A likelihood-based approach was also used to determine with more accuracy the extent of recombination in porB2 and porB3. Our study reveals that both genes are subject to exceptionally high rates of positive selection, as well as frequent recombination, which has important implications for the use of PorB as a potential vaccine candidate.


    Materials and Methods
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
Meningococcal porB Gene Sequences
A total of 324 porB gene sequences was examined: 121 from GenBank, including previously published sequences (Murakami, Gotschlich, and Seiff 1989Citation ; Wolff and Stern 1991Citation ; Feavers et al. 1992Citation ; Ward, Lambden, and Heckels 1992Citation ; Bash et al. 1995Citation ; Sacchi et al. 1998Citation ; Urwin et al. 1998aCitation , 1998bCitation ), and two from complete meningococcal genome sequences (Parkhill et al. 2000Citation ; Tettelin et al. 2000Citation ); a further 203 porB sequences from meningococcal isolates were determined de novo, including 107 isolates from globally representative strain collections, (Maiden et al. 1998Citation ), 12 PorB serotyping reference strains (Feavers et al. 1992Citation ), 50 isolates from healthy carriers, and 44 disease-causing meningococci from England and Wales.

For de novo sequencing, the propagation of isolates, preparation of DNA, porB gene amplification, and nucleotide sequence determination were as described previously (Feavers et al. 1999Citation ; Urwin et al. 1998aCitation ). The sequences were assembled with the Staden sequence analysis package (Staden 1996Citation ) and all sequences aligned manually in the Seqlab alignment program (Genetics Computer Group, Madison, Wis.) (Devereux, Haeberli, and Smithies 1984Citation ). The sequence alignment was trimmed at 5'- and 3'-ends so that all sequences began at the 5'-end with the thirteenth codon of the sequence encoding the mature protein (GAA in all sequences) and ended seven codons from the 3'-end of this sequence (ATG in porB2; GGT in porB3) because this corresponded to the length of the shortest sequences in the data set. Calculation of the number of nucleotide differences between pairs of alleles was determined using MEGA version 1.01 (Kumar, Tamura, and Nei 1994Citation ). The porB allele sequences and alignments can be viewed at http://neisseria.org/typing/porb. There were 125 unique porB allele sequences, 46 of which were porB2 sequences (named porB2-1 to porB2-46, according to a previously defined nomenclature [Feavers and Maiden 1998Citation ]) and 79, which were porB3 sequences (porB3-1 to porB3-79). Previously unpublished sequences have been submitted to GenBank, accession numbers AF520356AF520416.

Analysis of Selection Pressures
A maximum likelihood (ML) approach was used to examine selection pressures acting on the meningococcal porB genes. Here, dN and dS were examined codon-by-codon, using different models of codon substitution that differed in how dN/dS ratios (parameter {omega}) varied along sequences, as well as incorporating information about the phylogenetic relationships of the sequences in question so that comparisons are independent (Yang et al. 2000Citation ). Model M0 estimated a single {omega} parameter for all sites, whereas the M1 model divided codons into conserved sites (p0), with {omega}0 fixed at 0 and neutral sites (p1) with {omega}1 set to 1. The M2 model could account for positive selection through a third category of sites (p2) with {omega}2 estimated from the data. M3 provided a more sensitive test by estimating, from the data, {omega} values for three classes of site all of which could be >1. The M7 and M8 models both used a discrete beta distribution (with 10 categories and described by parameters p and q) to model {omega} ratios among sites, although M8, unlike M7, considered an extra class of sites for which {omega} could be >1. Nested models could be compared using a likelihood ratio test (LRT) in which twice the difference in log likelihood between models was compared with the value obtained under a {chi}2 distribution (degrees of freedom equal to the difference in the number of parameters between models). Finally, Bayesian methods were used to determine the probability that a particular codon site fell into the positively selected class. All these analyses used the CODEML program from the PAML package (Yang 1997Citation ).

Phylogenetic trees for the two data sets were constructed using the maximum likelihood method available in the PAUP* package (Swofford 1998Citation ). The HKY85 model of nucleotide substitution was used, with values for both the transition-transversion ratio and the shape parameter ({alpha}) of a gamma distribution of rate variation among sites (with eight categories) estimated during tree reconstruction (trees and parameter values available on request).

Recombination Analysis
The extent of recombination in the porB sequence data was analyzed by assessing the degree of phylogenetic congruence. Because of the large numbers of sequences available, this analysis was performed on a set of 35 randomly sampled porB2 and porB3 alleles. The porB2 and porB3 alignments were first split into two equal-sized fragments. ML phylogenetic trees were then estimated for both halves of the alignments using the procedures described above. To establish whether the trees constructed on each half of the alignment were significantly different in topology, as might be expected given frequent recombination, the difference in log likelihood ({delta}) between the ML tree for the first half of the gene and the ML tree topology for the second half of the gene fitted to the data from the first half, but with branch lengths reoptimized, were compared. The significance of the likelihood differences was assessed using two randomization tests. The first test used Monte Carlo simulation on 100 replicate data sets simulated under the ML model parameters using the program Seq-Gen (Rambaut and Grassly 1997Citation ). Maximum likelihood trees were then constructed for each of these simulated data sets, using the procedures described above, and their likelihoods compared with those of ML topology on each data set. If the {delta} values for the real data fell within this null distribution of {delta} values, then the trees constructed for each half of the gene were not significantly different in topology. In the second test, 200 random trees were created using PAUP*. The likelihoods of these trees were then estimated on the data from the first half of the porB2 and porB3 alignments, again with the reoptimization of branch lengths, and the {delta} values between these random trees and the two ML trees were then compared. If the {delta} values for the two ML trees fall within the 99th percentile of this null distribution then we may say that they are no more similar than two random trees inferred from these data (Holmes, Urwin, and Maiden 1999Citation ).

Structural Models
Structural models for the PorB2 and PorB3 proteins were generated using the software package Modeller (Sali et al. 1995Citation ), using the crystal structure of the porin Omp32 from Comamonas acidovorans as a template (Zeth et al. 2000Citation ; PDB accession 1E54). The sequence alignment was based on that given by Zeth et al. for PorB, with minor modifications.


    Results
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
Nucleotide and Amino Acid Sequence Diversity
A total of 14% (147 of 1,053, table 1 ) of the nucleotide sites of the 46 unique porB2 allele sequences were polymorphic, with percentage divergences (p distances) in pairwise comparisons ranging from 0.1% to 7.2% (mean 4.2%). For the 79 unique porB3 sequences, 15.2% (136 of 894, table 1 ) of the nucleotide sites were polymorphic, with p distances between 0.1% and 8.2% (mean 3.9%). Sequence divergence between the porB2 and porB3 allele classes ranged from 30.9% to 32.3%. The amino acid sequences of the putative surface exposed loop regions of the PorB2 and PorB3 proteins, determined according to the structural model for Neisseria porins, were longer in the PorB2 proteins. For the PorB2 proteins, putative loops V and VI were the most variable, with 25 distinct amino acid sequences (25.2% mean divergence) and 19 sequences (34.1% mean divergence), respectively. Putative loop VII was the most variable among PorB3 sequences, with 15 distinct amino acid sequences (35.9% mean divergence), although putative loops I, V, and VI were also variable, with 10 or more amino acid sequences (mean divergence >27%) in each loop region. The structural regions that interspersed the variable loops were largely conserved.


View this table:
[in this window]
[in a new window]
 
Table 1 Nucleotide and Amino Acid Sequence Diversity Among Meningococcal porB Genes

 
Analysis of Recombination Frequency
The maximum likelihood analysis of tree congruence provided evidence for recombination within both the porB2 and porB3 genes, although not at a rate sufficiently high to completely erase all phylogenetic similarity. Specifically, although the ML tree for the first half of the porB2 gene had a much higher likelihood than that of the ML tree for the second half of the gene fitted to the first half data ({delta} = 254.751), and this difference was far greater than expected by chance as determined using Monte Carlo simulation (max {delta} = 25.034), it was less than that seen in random tree topologies (99th percentile of {delta} = 327.565). Similar results were seen in porB3, where the difference between the two ML trees, {delta} = 350.571, also fell outside of the null distribution produced by Monte Carlo simulation (max {delta} = 32.876) but was less than that of random tree topologies (99th percentile of {delta} = 395.806).

Analysis of Selection Pressures Acting on the porB Genes of N. meningitidis
The maximum likelihood analysis of the selection pressures acting on the porB2 and porB3 alleles provided strong evidence for positive selection (tables 2 and 3 ). For both porB2 and porB3, the best-supported M3 model estimated {omega} (dN/dS) parameters >>1, indicative of strong positive selection. Similarly, the M8 model, which could incorporate positive selection, was significantly favored over the M7 model which did not, and estimated {omega} >> 1 for both porB2 and porB3. The strength of the inferred selection pressures acting on both porB2 and porB3 was also striking, as was the similarity in selection pressures between these allele classes. For the M3 model, two classes of positively selected sites were apparent. In the case of porB2 ~4.5% of sites fell into a relatively weakly positively selected class, where {omega}0 = 4.163, whereas ~1.1% of sites are seemingly subject to very strong positive selection with {omega}2 = 18.553. The remaining 94% of sites were highly conserved ({omega}1 = 0.067). A similar distribution of sites was apparent in porB3. Here, ~4.6% of sites had an {omega} value of 3.229, whereas 0.7% of sites were subject to much stronger positive selection pressures with {omega}2 = 13.923. In the case of porB3, 95% of sites were subject to strong selective constraints ({omega}1 = 0.033).


View this table:
[in this window]
[in a new window]
 
Table 2 Parameter Estimates for the Maximum Likelihood Analysis of Selection Pressures Acting on the porB2 Gene of Neisseria meningitidis

 

View this table:
[in this window]
[in a new window]
 
Table 3 Parameter Estimates for the Maximum Likelihood Analysis of Selection Pressures Acting on the porB3 Gene of Neisseria meningitidis

 
Bayesian methods were used to identify the sites with the highest probability of falling into the positively selected class under the M3 model (fig. 1 ). When these selected sites were mapped onto the three-dimensional structural models (fig. 2 ), it was clear that the majority fell within regions predicted to encode surface-exposed loop regions. Among PorB2 proteins there was evidence of weak or strong positive selection at multiple sites in all surface loops except putative loop II. Four residues in putative loop IV were shown to be under strong selection, six of the 23 amino acids in putative loop V were positively selected and all six selected sites in putative loop VI were under strong selection for amino acid change. For loops IV and V, there were nonselected residues at the center of the loop region, flanked by residues that were much more variable in sequence. Five selected sites were also identified within the structurally constrained membrane-spanning regions of PorB2. At position 39, the replacement of leucine with phenylalanine was always observed with a change at position 42 from methionine to isoleucine or valine. Among PorB3 sequences, most positively selected sites were located in putative surface loops I, V, VI, and VII. Seven sites in loop I were under positive selection compared with only one site in loop I of PorB2. No selected sites were identified in putative loop III, whereas three weakly selected sites were identified within the largely conserved membrane-spanning regions of PorB3.



View larger version (26K):
[in this window]
[in a new window]
 
Fig. 1.—Positive selection among PorB2 and PorB3 sequences. Translated (a) porB2-5 and (b) porB3-2 sequences are used as examples of PorB2 and PorB3 amino acid sequences. The locations of the putative surface-exposed loop sequences (I–VIII) are indicated. Amino acid residues that are subject to positive selection (probability >0.95) are shown in bold type. Sites under strong selection ({omega} = 18.553 and 13.923 for PorB2 and PorB3, respectively) are marked below the sequence with black blocks, and sites under weak selection ({omega} = 4.163 and 3.229 for PorB2 and PorB3, respectively) are denoted with white blocks

 


View larger version (42K):
[in this window]
[in a new window]
 
Fig. 2.—Ribbon diagrams of models for PorB2 (top) and PorB3 (bottom) with superposition of residues subject to positive selection. Residues under strong selection are shown in red and residues under weak selection are shown in yellow. The diagrams were produced using Molscript (Kraulis 1991)

 
Finally, because the maximum likelihood approach of Yang et al. (2000)Citation explicitly assumes a phylogenetic tree when estimating selection pressures, we attempted to assess whether high rates of recombination could have produced a false-positive signal for positive selection. We therefore repeated the analysis assuming that sequences were linked by a "star" phylogeny, where the lineages diverge simultaneously from a single root node, thereby removing the affect of phylogenetic history. For both PorB2 and PorB3, this analysis again produced significant evidence for positive selection (P < 0.000 for M3 vs. M2 and M8 vs. M7) and for extremely high {omega} values at the selected sites—maximum values of 75.052 and 11.305, respectively, under the M3 model, with generally the same sites falling into the positively selected class as in the original analysis (full results available from authors on request).


    Discussion
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
Here phylogenetic and structural modeling techniques have been used to reveal evolutionary processes acting on an established, but controversial, vaccine candidate, the meningococcal PorB protein. Maximum likelihood analysis of selection pressures provided powerful evidence for adaptive evolution of the porB genes of N. meningitidis. In both porB2 and porB3, codons were identified that had been subject to very strong selection pressure (dN/dS of 18.533 and 16.221, respectively). These values were some of the highest seen in any gene studied to date, and similar to the values observed in studies of the intrahost evolution of the HIV-1 env gene, often seen as a paradigm of immune driven positive selection (Nielson and Yang 1998Citation ; Zanotto et al. 1999Citation ). The fact that positive selection was still observed at these same sites after removing the effect of phylogenetic history by assuming a star phylogeny indicated that the analysis was unlikely to have been greatly biased by recombination. Mapping of these sites onto refined models of the protein structure of PorB confirmed that the majority of selected sites were located in regions of the protein predicted to be exposed to the host immune system. These results contradict those obtained previously (Smith, Maynard Smith, and Spratt 1995Citation ), which provided no evidence for positive selection in the putative loop regions of meningococcal PorB proteins, and in doing so highlight the need to consider each codon separately within amino acid alignments.

In addition, the analysis also revealed that although the PorB2 and PorB3 proteins were subject to very similar selection pressures they exhibited different distributions of positively selected sites. These findings were supported by flow cytometry analysis of live meningococci which demonstrated that some PorB3 variants were not easily accessible for antibody binding (Michaelsen et al. 2001Citation ), possibly because of shielding or due to PorB3 extracellular loops being shorter than those of PorB2, so that fewer residues were subject to intense immune selection. In contrast, the gonococcal porins PIA and PIB showed no significant difference in the distribution of selected sites despite reported differences in selection intensities between the two homology groups (Posada et al. 2000Citation ). Posada et al. therefore concluded that epidemiological differences between gonococci expressing PIA or PIB proteins were responsible for differences in selection rather than structural differences in the proteins. In meningococci, however, it is possible that the differences in the lengths of the surface loops of PorB2 and PorB3 porins are sufficient to affect the conformation and structure of epitopes presented by these proteins and that this will determine which sites are exposed to selective pressure from the host responses. Furthermore, although there may be some epidemiological differences between meningococci expressing PorB2 and PorB3 proteins, there is no evidence from these data that this has led to concomitant selective differences, although additional studies are required to test this hypothesis further.

The conformational effects of particular amino acid substitutions on PorB structure remain difficult to determine. The homology models used were useful for examining the approximate disposition of residues in space and could provide insights not readily apparent from a sequence alignment. The crystal structure of Omp32 from C. acidovorans was used as the basis for the homology models of PorB2 and PorB3; this is a closer homolog to the Neisseria porins than the E. coli porin crystal structures that were used previously (Derrick et al. 1999Citation ). But the size of most of the external loop regions precluded an accurate estimation of their conformations by standard homology modeling techniques, and the location of selected residues shown in external loop regions in figure 2 were approximate.

Figures 1 and 2 illustrate that the strongly and weakly selected sites are not distributed evenly across the loop regions. The presence of conserved residues at the apices of the PorB2 variable loops I, IV, and V suggested that these amino acids were not exposed to the host immune response because of protein folding, or perhaps they fulfill an important role retaining the surface loop structure. For example, the L2 loop in the OmpF protein of E. coli contributes to the stabilization of the porin trimer (Phale et al. 1998Citation ) and residues in other loop regions could play analogous roles in stabilizing the protein, using loop-loop interactions. Indeed, although the portions of polypeptide chain joining the ends of the ß-barrel strands are frequently referred to as loops, they are likely to contain regions of regular secondary structure, as is seen in other porin structures (Koebnik, Locher, and Van Gelder 2000Citation ), and this would place constraints on the sequence variation within these regions in PorB2 and PorB3 proteins. These observations are borne out by the difficulty in mimicking PorB epitopes with linear peptides (Zapata et al. 1992Citation ) and by sequence variation identified among serologically similar PorB antigens (Urwin et al. 1998aCitation , 1998bCitation ).

Although most residues within the transmembrane ß-strands of PorB2 were highly conserved, five positively selected sites were located among these structural regions. According to the structural model, two of these sites (positions 39 and 43 in fig. 1a ) were located close to one another on opposite sides of a ß-turn and were identified as positively selected with >99% probability (although falling into the weakly rather than the strongly positively selected class). A leucine residue at position 39 was invariably accompanied by a methionine residue at position 43, whereas the phenylalanine residue at position 39 was accompanied by the smaller hydrophobic residues valine or isoleucine at position 43. Within the limitations of the homology model, the side chains of the two residues may be in contact, suggesting that these are compensatory mutations. Elsewhere, the reasons for strong positive selection within the ß-barrel region were more difficult to discern. The PorB2 protein contained a number of mutations within the L3 loop, which folds back into the center of the barrel creating a constricted channel for the passage of low molecular weight solutes (Zeth et al. 2000Citation ). Residues within the L3 loop were unlikely to be subject to immune selection but are probably involved in ion selectivity: the presence of a number of mutations which resulted in changes in charge within this region of the protein are consistent with this idea (Bauer et al. 1989Citation ; Benz et al. 1989Citation ; Saint et al. 1996Citation ; Schirmer and Phale 1999Citation ).

Despite the extensive evidence for recombinational reassortment in PorB proteins, (Vázquez et al. 1995Citation ; Cooke et al. 1998Citation ; Derrick et al. 1999Citation ), this has not been so frequent as to remove all phylogenetic signals. In contrast, congruence tests between different housekeeping genes from N. meningitidis, subject to strong purifying selection, revealed that most gene trees were no more similar to one another than expected by chance (Holmes, Urwin, and Maiden 1999Citation ). Overall, the diversity of meningococcal protein antigens explored here exhibits features which phylogenetic and structural models can exploit in the elucidation of the function and vaccine potential of these molecules. The insights obtained can complement data from both genomic studies and experimental studies of immunogenicity in man and animals to provide more complete information on the interactions of pathogen proteins with the hosts and its immune system.


    Acknowledgements
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 
This work was supported by The Wellcome Trust (047072), The Royal Society, The Lister Institute, and the United Kingdom Public Health Laboratory Service. We thank Dr. Z. Yang for useful advice regarding CODEML.


    Footnotes
 
Peer Bork, Reviewing Editor

Keywords: Neisseria meningitidis porB evolution recombination selection Back

Address for correspondence and reprints: Martin C. J. Maiden, Department of Zoology, University of Oxford, South Parks Road, Oxford OX1 3PS, U.K. martin.maiden{at}zoo.ox.ac.uk Back


    References
 TOP
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 Acknowledgements
 References
 

    Achtman M., 1995 Global epidemiology of meningococcal disease Pp 159–175 in K. A. V. Cartwright, ed. Meningococcal disease. John Wiley and Sons, Chichester, U.K

    Bash M. C., K. B. Lesiak, S. D. Banks, C. E. Frasch, 1995 Analysis of Neisseria meningitidis class 3 outer membrane protein gene variable regions and type identification using genetic techniques Infect. Immun 63:1484-1490[Abstract]

    Bauer K., M. Struyve, D. Bosch, R. Benz, J. Tommassen, 1989 One single lysine residue is responsible for the special interaction between polyphosphate and the outer-membrane porin PhoE of Escherichia coli J. Biol. Chem 264:16393-16398[Abstract/Free Full Text]

    Benz R., A. Schmid, P. Van der Ley, J. Tommassen, 1989 Molecular-basis of porin selectivity—membrane experiments with OmpC-PhoE and OmpF-PhoE hybrid proteins of Escherichia coli K-12 Biochim. Biophys. Acta 981:8-14[ISI][Medline]

    Bjune G., J. K. Gronnesby, E. A. Høiby, O. Closs, H. Nokleby, 1991a. Results of an efficacy trial with an outer membrane vesicle vaccine against systemic serogroup B meningococcal disease in Norway NIPH Ann 14:125-130[Medline]

    Bjune G., E. A. Høiby, J. K. Grønnesby, et al. (16 co-authors) 1991b. Effect of outer membrane vesicle vaccine against group B meningococcal disease in Norway Lancet 338:1093-1096.[ISI][Medline]

    Blake M. S., E. C. Gotschlich, 1987 Functional and immunologic properties of pathogeneic Neisserial surface proteins Pp 377–400 in M. Inouye, ed. Bacterial outer membranes as model systems. John Wiley and Sons, Chichester, U.K

    Caugant D. A., L. F. Mocca, C. E. Frasch, L. O. Frøholm, W. D. Zollinger, R. K. Selander, 1987 Genetic structure of Neisseria meningitidis populations in relation to serogroup, serotype, and outer membrane protein pattern J. Bacteriol 169:2781-2792[ISI][Medline]

    Cooke S. J., K. Jolley, C. A. Ison, H. Young, J. E. Heckels, 1998 Naturally occurring isolates of Neisseria gonorrhoeae, which display anomalous serovar properties, express PIA/PIB hybrid porins, deletions in PIB or novel PIA molecules FEMS Microbiol. Lett 162:75-82[ISI][Medline]

    Deitsch K. W., E. R. Moxon, T. E. Wellems, 1997 Shared themes of antigenic variation and virulence in bacterial, protozoal, and fungal infections Microbiol. Mol. Biol. Rev 61:281-293[Abstract]

    Derrick J. P., R. Urwin, J. Suker, I. M. Feavers, M. C. J. Maiden, 1999 Structural and evolutionary inference from molecular variation in Neisseria porins Infect. Immun 67:2406-2413[Abstract/Free Full Text]

    Devereux J. P., P. Haeberli, O. Smithies, 1984 A comprehensive set of sequence analysis programs for the VAX Nucleic Acids Res 12:387-395[Abstract]

    Feavers I. M., S. J. Gray, R. Urwin, J. E. Russell, J. A. Bygraves, E. B. Kaczmarski, M. C. J. Maiden, 1999 Multilocus sequence typing and antigen gene sequencing in the investigation of a meningococcal disease outbreak J. Clin. Microbiol 37:3883-3887[Abstract/Free Full Text]

    Feavers I. M., M. C. J. Maiden, 1998 A gonococcal porA pseudogene: implications for understanding the evolution and pathogenicity of Neisseria gonorrhoeae Mol. Microbiol 30:647-656[ISI][Medline]

    Feavers I. M., J. Suker, A. J. McKenna, A. B. Heath, M. C. J. Maiden, 1992 Molecular analysis of the serotyping antigens of Neisseria meningitidis Infect. Immun 60:3620-3629[Abstract]

    Frasch C. E., W. D. Zollinger, J. T. Poolman, 1985 Serotype antigens of Neisseria meningitidis and a proposed scheme for designation of serotypes Rev. Infect. Dis 7:504-510[ISI][Medline]

    Holmes E. C., R. Urwin, M. C. J. Maiden, 1999 The influence of recombination on the population structure and evolution of the human pathogen Neisseria meningitidis Mol. Biol. Evol 16:741-749[Abstract]

    Kleffel B., R. M. Garavito, W. Baumeister, J. P. Rosenbusch, 1985 Secondary structure of a channel-forming protein: porin from E. coli outer membranes EMBO J 4:1589-1592[Abstract]

    Koebnik R., K. P. Locher, P. Van Gelder, 2000 Structure and function of bacterial outer membrane proteins: barrels in a nutshell Mol. Microbiol 37:239-253[ISI][Medline]

    Kraulis P. J., 1991 MOLSCRIPT: a program to produce both detailed and schematic plots of protein structures J. Appl. Cryst 24:946-950

    Kumar S., K. Tamura, M. Nei, 1994 MEGA: molecular evolutionary genetics analysis software for microcomputers Comput. Appl. Biosci 10:189-191[Abstract]

    Maiden M. C. J., J. A. Bygraves, E. Feil, et al. (12 co-authors) 1998 Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms Proc. Natl. Acad. Sci. USA 95:3140-3145.[Abstract/Free Full Text]

    Maiden M. C. J., J. Suker, A. J. McKenna, J. A. Bygraves, I. M. Feavers, 1991 Comparison of the class 1 outer membrane proteins of eight serological reference strains of Neisseria meningitidis Mol. Microbiol 5:727-736[ISI][Medline]

    McGuinness B., A. K. Barlow, I. N. Clarke, J. E. Farley, A. Anilionis, J. T. Poolman, J. E. Heckels, 1990 Deduced amino acid sequences of class 1 protein (PorA) from three strains of Neisseria meningitidis. Synthetic peptides define the epitopes responsible for serosubtype specificity J. Exp. Med 171:1871-1882[Abstract]

    Michaelsen T. E., A. Aase, J. Kolberg, E. Wedege, E. Rosenqvist, 2001 PorB3 outer membrane protein on Neisseria meningitidis is poorly accessible for antibody binding on live bacteria Vaccine 19:1526-1533[ISI][Medline]

    Murakami K., E. C. Gotschlich, M. E. Seiff, 1989 Cloning and characterisation of the structural gene for the class 2 protein of Neisseria meningitidis Infect. Immun 57:2318-2323[ISI][Medline]

    Nielson R., Z. Yang, 1998 Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene Genetics 148:929-936[Abstract/Free Full Text]

    Parkhill J., M. Achtman, K. D. James, et al. (27 co-authors) 2000 Complete DNA sequence of a serogroup A strain of Neisseria meningitidis Z2491 Nature 404:502-506.[ISI][Medline]

    Perkins B. A., K. Jonsdottir, H. Briem, et al. (18 co-authors) 1998 Immunogenicity of two efficacious outer membrane protein-based serogroup B meningococcal vaccines among young adults in Iceland J. Infect. Dis 177:683-691.[ISI][Medline]

    Phale P. S., A. Philippsen, T. Kiefhaber, R. Koebnik, V. P. Phale, T. Schirmer, J. P. Rosenbusch, 1998 Stability of trimeric OmpF porin: the contributions of the latching loop L2 Biochemistry 37:15663-15670[ISI][Medline]

    Pollard A. J., C. Frasch, 2001 Development of natural immunity to Neisseria meningitidis Vaccine 19:1327-1346[ISI][Medline]

    Poolman J. T., I. Lind, K. Jonsdottir, L. O. Frøholm, D. M. Jones, H. C. Zanen, 1986 Meningococcal serotypes and serogroup B disease in north-west Europe Lancet ii:555-557

    Posada D., K. A. Crandall, M. Nguyen, J. C. Demma, R. P. Viscidi, 2000 Population genetics of the porB gene of Neisseria gonorrhoeae: different dynamics in different homology groups Mol. Biol. Evol 17:423-436[Abstract/Free Full Text]

    Rambaut A., N. C. Grassly, 1997 Seq-Gen: an application for the Monte Carlo simulation of DNA sequence evolution along phylogenetic trees Comput. Appl. Biosci 13:235-238[Abstract]

    Rosenqvist E., E. Arne Høiby, E. Wedege, K. Bryn, J. Kolberg, A. Klem, E. Rønnild, G. Bjune, H. Nøkleby, 1995 Human antibody responses to meningococcal outer membrane antigens after three doses of the Norwegian group B meningococcal vaccine Infect. Immun 63:4642-4652[Abstract]

    Rudel T., A. Schmid, R. Benz, H. A. Kolb, F. Lang, T. F. Meyer, 1996 Modulation of Neisseria porin (PorB) by cytosolic ATP/GTP of target cells: parallels between pathogen accommodation and mitochondrial endosymbiosis Cell 85:391-402[ISI][Medline]

    Sacchi C. T., A. P. S. Lemos, A. M. Whitney, C. A. Solari, M. E. Brandt, C. E. A. Melles, C. E. Frasch, L. W. Mayer, 1998 Correlation between serological and sequencing analyses of the PorB outer membrane protein in the Neisseria meningitidis serotyping system Clin. Diagn. Lab. Immunol 5:348-354[Abstract]

    Saint N., K. L. Lou, C. Widmer, M. Luckey, T. Schirmer, J. P. Rosenbusch, 1996 Structural and functional characterization of OmpF porin mutants selected for larger pore size 2. Functional characterization J. Biol. Chem 271:20676-20680[Abstract/Free Full Text]

    Sali A., L. Potterton, F. Yuan, H. van Vlijmen, M. Karplus, 1995 Evaluation of comparative protein modeling by MODELLER Proteins 23:318-326[ISI][Medline]

    Saukkonen K., H. Abdillahi, J. T. Poolman, M. Leinonen, 1987 Protective efficacy of monoclonal antibodies to class 1 and class 3 outer membrane proteins of Neisseria meningitidis B: 15: P1.16 in infant rat infection model: new prospects for vaccine development Microb. Pathog 3:261-267[ISI][Medline]

    Saukkonen K., M. Leinonen, H. Abdillahi, J. T. Poolman, 1989 Comparative evaluation of potential components for group B meningococcal vaccine by passive protection in the infant rat and in vitro bactericidal assay Vaccine 7:325-328[ISI][Medline]

    Schirmer T., P. S. Phale, 1999 Brownian dynamics simulation of ion flow through porin channels J. Mol. Biol 294:1159-1167[ISI][Medline]

    Sierra G. V. G., H. C. Campa, N. M. Varcacel, I. L. Garcia, P. L. Izquierdo, P. F. Sotolongo, G. V. Casanueva, C. O. Rico, C. R. Rodriguez, M. H. Terry, 1991 Vaccine against group B Neisseria meningitidis: protection trial and mass vaccination results in Cuba National Institute for Public Health Annals 14:195-207.

    Smith N. H., J. Maynard Smith, B. G. Spratt, 1995 Sequence evolution of the porB gene of Neisseria gonorrhoeae and Neisseria meningitidis: evidence of positive Darwinian selection Mol. Biol. Evol 12:363-370[Abstract]

    Staden R., 1996 The Staden sequence analysis package Mol. Biotechnol 5:233-241[ISI][Medline]

    Swofford D., 1998 PAUP*: phylogenetic analysis using parsimony and other methods Sinauer, Sunderland, Mass

    Tettelin H., N. J. Saunders, J. Heidelberg, et al. (41 co-authors) 2000 Complete genome sequence of Neisseria meningitidis serogroup B strain MC58 Science 287:1809-1815.[Abstract/Free Full Text]

    Tommassen J., P. Vermeij, M. Struyvé, R. Benz, J. T. Poolman, 1990 Isolation of Neisseria meningitidis mutants deficient in class 1 (PorA) and class 3 (PorB) outer membrane proteins Infect. Immun 58:1355-1359[ISI][Medline]

    Urwin R., I. M. Feavers, D. M. Jones, M. C. J. Maiden, A. J. Fox, 1998a. Molecular variation of meningococcal serotype 4 antigen genes Epidemiol. Infect 121:95-101[ISI][Medline]

    Urwin R., A. J. Fox, M. Musilek, P. Kriz, M. C. J. Maiden, 1998b. Heterogeneity of the PorB protein in serotype 22 Neisseria meningitidis J. Clin. Microbiol 36:3680-3682[Abstract/Free Full Text]

    van der Ende A., C. T. P. Hopman, S. Zaat, B. B. Oude Essink, B. Berkhout, J. Dankert, 1995 Variable expression of class 1 outer membrane protein in Neisseria meningitidis is caused by variation in the spacing between the -10 and -35 regions of the promoter J. Bacteriol 177:2475-2480[Abstract]

    van der Ley P., J. E. Heckels, M. Virji, P. Hoogerhout, J. T. Poolman, 1991 Topology of outer membrane porins in pathogenic Neisseria species Infect. Immun 59:2963-2971[ISI][Medline]

    van der Ley P., J. T. Poolman, 1992 Construction of a multivalent meningococcal vaccine strain based on the class 1 outer membrane protein Infect. Immun 60:3156-3161[Abstract]

    van der Ley P., J. van der Biezen, J. T. Poolman, 1995 Construction of Neisseria meningitidis strains carrying multiple chromosomal copies of the porA gene for use in the production of a multivalent outer membrane vesicle vaccine Vaccine 13:401-407[ISI][Medline]

    Vázquez J. A., S. Berron, M. O'Rourke, G. Carpenter, E. Feil, N. H. Smith, B. G. Spratt, 1995 Interspecies recombination in nature: a meningococcus that has acquired a gonococcal PIB porin Mol. Microbiol 15:1001-1007[ISI][Medline]

    Ward M. J., P. R. Lambden, J. E. Heckels, 1992 Sequence analysis and relationships between meningococcal class 3 serotype proteins and other porins of pathogenic and non-pathogenic Neisseria species FEMS Microbiol. Lett 94:283-290[ISI]

    Wolff K., A. Stern, 1991 The class 3 outer membrane protein (PorB) of Neisseria meningitidis: gene sequence and homology to the gonococcal porin PIA FEMS Microbiol. Lett 83:179-185[ISI]

    Yang Z., 1997 PAML: a program package for phylogenetic analysis by maximum likelihood Comput. Appl. Biol. Sci 13:555-556[Medline]

    Yang Z., R. Nielsen, N. Goldman, A. M. K. Pedersen, 2000 Codon-substitution models for heterogenous selection pressure at amino acid sites Genetics 155:431-449[Abstract/Free Full Text]

    Zanotto P. M. de, E. G. Kallas, R. F. de Souza, E. C. Holmes, 1999 Genealogical evidence for positive selection in the nef gene of HIV-1 Genetics 153:1077-1089.[Abstract/Free Full Text]

    Zapata G. A., W. F. Vann, Y. Rubinstein, C. E. Frasch, 1992 Identification of variable region differences in Neisseria meningitidis class 3 protein sequences among five group B serotypes Mol. Microbiol 6:3493-3499[ISI][Medline]

    Zeth K., K. Diederichs, W. Welte, H. Engelhardt, 2000 Crystal structure of Omp32, the anion-selective porin from Comamonas acidovorans, in complex with a periplasmic peptide at 2.1 Å resolution Structure 8:981-992[ISI][Medline]

Accepted for publication May 6, 2002.