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

Evolutionary and structural aspects of Solanaceae RNases T2

Plant RNases T2 are involved in several physiological and developmental processes, including inorganic phosphate starvation, senescence, wounding, defense against pathogens, and the self-incompatibility system. Solanaceae RNases form three main clades, one composed exclusively of S-RNases and two that include S-like RNases. We identified several positively selected amino acids located in highly flexible regions of these molecules, mainly close to the B1 and B2 substrate-binding sites in S-like RNases and the hypervariable regions of S-RNases. These differences between S- and S-like RNases in the flexibility of amino acids in substrate-binding regions are essential to understand the RNA-binding process. For example, in the S-like RNase NT, two positively selected amino acid residues (Tyr156 and Asn134) are located at the most flexible sites on the molecular surface. RNase NT is induced in response to tobacco mosaic virus infection; these sites may thus be regions of interaction with pathogen proteins or viral RNA. Differential selective pressures acting on plant ribonucleases have increased amino acid variability and, consequently, structural differences within and among S-like RNases and S-RNases that seem to be essential for these proteins play different functions.

Pathogenic Neisseria Bind the Complement Protein CFHR5 via Outer Membrane Porins

Neisseria meningitidis and Neisseria gonorrhoeae are important human pathogens that have evolved to bind the major negative regulator of the complement system, complement factor H (CFH). However, little is known about the interaction of pathogens with CFH-related proteins (CFHRs) which are structurally similar to CFH but lack the main complement regulatory domains found in CFH. Insights into the role of CFHRs have been hampered by a lack of specific reagents. We generated a panel of CFHR-specific monoclonal antibodies and demonstrated that CFHR5 was bound by both pathogenic Neisseria spp. We showed that CFHR5 bound to PorB expressed by both pathogens in the presence of sialylated lipopolysaccharide and enhanced complement activation on the surface of N. gonorrhoeae. Our study furthered our understanding of the interactions of CFHRs with bacterial pathogens and revealed that CFHR5 bound the meningococcus and gonococcus via similar mechanisms.

Mosaic Evolution of Beta-Barrel-Porin-Encoding Genes in Escherichia coli

Bacterial porin-encoding genes are often found under positive selection. Local recombination has also been identified in a few of them to facilitate bacterial rapid adaptation, although it remains unknown whether it is a common evolutionary mechanism for the porins or outer membrane proteins in Gram-negative bacteria. In this study, we investigated the beta-barrel (β-barrel) porin-encoding genes in Escherichia coli that were reported under positive Darwinian selection. Besides fhuA that was found with ingenic local recombination previously, we identified four other genes, i.e., lamB, ompA, ompC, and ompF, all showing the similar mosaic evolution patterns. Comparative analysis of the protein sequences disclosed a list of highly variable regions in each family, which are mostly located in the convex of extracellular loops and coinciding with the binding sites of bacteriophages. For each of the porin families, mosaic recombination leads to unique combinations of the variable regions with different sequence patterns, generating diverse protein groups. Structural modeling indicated a conserved global topology among the different porins, with the extracellular surface varying a lot due to individual or combinatorial variable regions. The conservation of global tertiary structure would ensure the channel activity, while the wide diversity of variable regions may represent selection to avoid the invasion of phages, antibiotics or immune surveillance factors. Our study identified multiple bacterial porin genes with mosaic evolution. We hypothesize that this could be generalized strategy for outer membrane proteins to both maintain normal life processes and evade the attack of unfavored factors rapidly. IMPORTANCE Microevolution studies can disclose more elaborate evolutionary mechanisms of genes, appearing especially important for genes with multifaceted function such as those encoding outer membrane proteins. However, in most cases, the gene is considered as a whole unit, and the evolutionary patterns are disclosed. Here, we report that multiple bacterial porin proteins follow mosaic evolution, with local ingenic recombination combined with spontaneous mutations based on positive Darwinian selection, and conservation for most structural regions. This could represent a common mechanism for bacterial outer membrane proteins. The variable regions within each porin family showed large coincidence with the binding sites of bacteriophages, antibiotics, and immune factors and therefore would represent effective targets for the development of new antibacterial agents or vaccines.

GenomegaMap: Within-Species Genome-Wide dN/dS Estimation from over 10,000 Genomes

The dN/dS ratio provides evidence of adaptation or functional constraint in protein-coding genes by quantifying the relative excess or deficit of amino acid-replacing versus silent nucleotide variation. Inexpensive sequencing promises a better understanding of parameters, such as dN/dS, but analyzing very large data sets poses a major statistical challenge. Here, I introduce genomegaMap for estimating within-species genome-wide variation in dN/dS, and I apply it to 3,979 genes across 10,209 tuberculosis genomes to characterize the selection pressures shaping this global pathogen. GenomegaMap is a phylogeny-free method that addresses two major problems with existing approaches: 1) It is fast no matter how large the sample size and 2) it is robust to recombination, which causes phylogenetic methods to report artefactual signals of adaptation. GenomegaMap uses population genetics theory to approximate the distribution of allele frequencies under general, parent-dependent mutation models. Coalescent simulations show that substitution parameters are well estimated even when genomegaMap’s simplifying assumption of independence among sites is violated. I demonstrate the ability of genomegaMap to detect genuine signatures of selection at antimicrobial resistance-conferring substitutions in Mycobacterium tuberculosis and describe a novel signature of selection in the cold-shock DEAD-box protein A gene deaD/csdA. The genomegaMap approach helps accelerate the exploitation of big data for gaining new insights into evolution within species.

Genome-wide analyses reveal genes subject to positive selection in Toxoplasma gondii

Toxoplasma gondii is an important protozoan pathogen that infects many wild and domestic animals and causes infections in immunocompromised humans. However, there has been little investigation of the molecular evolutionary trajectories of this pathogenic protozoa using comparative genomics data. Here, we employed a comparative evolutionary genomics approach to identify genes that are under site- and lineage-specific positive selection in nine strains of T. gondii, including two closely related species, Neospora caninum and Hammondia hammondi. Based on the analyses of five coccidian core genomes, 4.5% of the 5788 core genome genes showed strong signals for positive selection in the site model. In addition, the branch-site model analyses in the nine T. gondii core genomes indicated that 2 to 20 genes underwent significant positive selection along each lineage leading to T. gondii strains. Many of the protein products encoded by the positively selected genes are secretory or surface proteins that have previously been implicated in host pathogenesis. The adaptive changes in these positively selected genes might be related to dynamic interactions between the host immune systems and might play a crucial role in the infection and pathogenic processes of T. gondii.

Heterogeneity in non-epitope loop sequence and outer membrane protein complexes alters antibody binding to the major porin protein PorB in serogroup B Neisseria meningitidis

PorB is a well-characterized outer membrane protein that is common among Neisseria species and is required for survival. A vaccine candidate, PorB induces antibody responses that are directed against six variable surface-exposed loops that differ in sequence depending on serotype. Although Neisseria meningitidis is naturally competent and porB genetic mosaicism provides evidence for strong positive selection, the sequences of PorB serotypes commonly associated with invasive disease are often conserved, calling into question the interaction of specific PorB loop sequences in immune engagement. In this report, we provide evidence that antibody binding to a PorB epitope can be altered by sequence mutations in non-epitope loops. Through the construction of hybrid PorB types and PorB molecular dynamics simulations, we demonstrate that loops both adjacent and non-adjacent to the epitope loop can enhance or diminish antibody binding, a phenotype that correlates with serum bactericidal activity. We further examine the interaction of PorB with outer membrane-associated proteins, including PorA and RmpM. Deletion of these proteins alters the composition of PorB-containing native complexes and reduces antibody binding and serum killing relative to the parental strain, suggesting that both intramolecular and intermolecular PorB interactions contribute to host adaptive immune evasion.

Tracing evolutionary relicts of positive selection on eight malaria-related immune genes in mammals

Plasmodium-induced malaria widely infects primates and other mammals. Multiple past studies have revealed that positive selection could be the main evolutionary force triggering the genetic diversity of anti-malaria resistance-associated genes in human or primates. However, researchers focused most of their attention on the infra-generic and intra-specific genome evolution rather than analyzing the complete evolutionary history of mammals. Here we extend previous research by testing the evolutionary link of natural selection on eight candidate genes associated with malaria resistance in mammals. Three of the eight genes were detected to be affected by recombination, including TNF-α, iNOS and DARC. Positive selection was detected in the rest five immunogenes multiple times in different ancestral lineages of extant species throughout the mammalian evolution. Signals of positive selection were exposed in four malaria-related immunogenes in primates: CCL2, IL-10, HO1 and CD36. However, selection signals of G6PD have only been detected in non-primate eutherians. Significantly higher evolutionary rates and more radical amino acid replacement were also detected in primate CD36, suggesting its functional divergence from other eutherians. Prevalent positive selection throughout the evolutionary trajectory of mammalian malaria-related genes supports the arms race evolutionary hypothesis of host genetic response of mammalian immunogenes to infectious pathogens.

Receptor-like kinases in plant innate immunity

Plants employ a highly effective surveillance system to detect potential pathogens, which is critical for the success of land plants in an environment surrounded by numerous microbes. Recent efforts have led to the identification of a number of immune receptors and components of immune receptor complexes. It is now clear that receptor-like kinases (RLKs) and receptor-like proteins (RLPs) are key pattern-recognition receptors (PRRs) for microbe- and plant-derived molecular patterns that are associated with pathogen invasion. RLKs and RLPs involved in immune signaling belong to large gene families in plants and have undergone lineage specific expansion. Molecular evolution and population studies on phytopathogenic molecular signatures and their receptors have provided crucial insight into the co-evolution between plants and pathogens. [Figure: see text] Jian-Min Zhou (Corresponding author).

The role of selection in shaping diversity of natural M. tuberculosis populations

Mycobacterium tuberculosis (M.tb), the cause of tuberculosis (TB), is estimated to infect a new host every second. While analyses of genetic data from natural populations of M.tb have emphasized the role of genetic drift in shaping patterns of diversity, the influence of natural selection on this successful pathogen is less well understood. We investigated the effects of natural selection on patterns of diversity in 63 globally extant genomes of M.tb and related pathogenic mycobacteria. We found evidence of strong purifying selection, with an estimated genome-wide selection coefficient equal to -9.5 × 10(-4) (95% CI -1.1 × 10(-3) to -6.8 × 10(-4)); this is several orders of magnitude higher than recent estimates for eukaryotic and prokaryotic organisms. We also identified different patterns of variation across categories of gene function. Genes involved in transport and metabolism of inorganic ions exhibited very low levels of non-synonymous polymorphism, equivalent to categories under strong purifying selection (essential and translation-associated genes). The highest levels of non-synonymous variation were seen in a group of transporter genes, likely due to either diversifying selection or local selective sweeps. In addition to selection, we identified other important influences on M.tb genetic diversity, such as a 25-fold expansion of global M.tb populations coincident with explosive growth in human populations (estimated timing 1684 C.E., 95% CI 1620-1713 C.E.). These results emphasize the parallel demographic histories of this obligate pathogen and its human host, and suggest that the dominant effect of selection on M.tb is removal of novel variants, with exceptions in an interesting group of genes involved in transportation and defense. We speculate that the hostile environment within a host imposes strict demands on M.tb physiology, and thus a substantial fitness cost for most new mutations. In this respect, obligate bacterial pathogens may differ from other host-associated microbes such as symbionts.

Increased knowledge of Francisella genus diversity highlights the benefits of optimised DNA-based assays

Background

Recent advances in sequencing technologies offer promising tools for generating large numbers of genomes, larger typing databases and improved mapping of environmental bacterial diversity. However, DNA-based methods for the detection of Francisella were developed with limited knowledge about genetic diversity. This, together with the high sequence identity between several Francisella species, means there is a high risk of false identification and detection of the highly virulent pathogen Francisella tularensis. Moreover, phylogenetic reconstructions using single or limited numbers of marker sequences often result in incorrect tree topologies and inferred evolutionary distances. The recent growth in publicly accessible whole-genome sequences now allows evaluation of published genetic markers to determine optimal combinations of markers that minimise both time and laboratory costs.

Results

In the present study, we evaluated 38 previously published DNA markers and the corresponding PCR primers against 42 genomes representing the currently known diversity of the genus Francisella. The results highlight that PCR assays for Francisella tularensis are often complicated by low specificity, resulting in a high probability of false positives. A method to select a set of one to seven markers for obtaining optimal phylogenetic resolution or diagnostic accuracy is presented.

Conclusions

Current multiple-locus sequence-typing systems and detection assays of Francisella, could be improved by redesigning some of the primers and reselecting typing markers. The use of only a few optimally selected sequence-typing markers allows construction of phylogenetic topologies with almost the same accuracy as topologies based on whole-genome sequences.

The amino acid sequence of Neisseria lactamica PorB surface-exposed loops influences Toll-like receptor 2-dependent cell activation

Toll-like receptors (TLRs) play a major role in host mucosal and systemic defense mechanisms by recognizing a diverse array of conserved pathogen-associated molecular patterns (PAMPs). TLR2, with TLR1 and TLR6, recognizes structurally diverse bacterial products such as lipidated factors (lipoproteins and peptidoglycans) and nonlipidated proteins, i.e., bacterial porins. PorB is a pan-neisserial porin expressed regardless of organisms' pathogenicity. However, commensal Neisseria lactamica organisms and purified N. lactamica PorB (published elsewhere as Nlac PorB) induce TLR2-dependent proinflammatory responses of lower magnitude than N. meningitidis organisms and N. meningitidis PorB (published elsewhere as Nme PorB). Both PorB types bind to TLR2 in vitro but with different apparent specificities. The structural and molecular details of PorB-TLR2 interaction are only beginning to be unraveled and may be due to electrostatic attraction. PorB molecules have significant strain-specific sequence variability within surface-exposed regions (loops) putatively involved in TLR2 interaction. By constructing chimeric recombinant PorB loop mutants in which surface-exposed loop residues have been switched between N. lactamica PorB and N. meningitidis PorB, we identified residues in loop 5 and loop 7 that influence TLR2-dependent cell activation using HEK cells and BEAS-2B cells. These loops are not uniquely responsible for PorB interaction with TLR2, but NF-κB and MAP kinases signaling downstream of TLR2 recognition are likely influenced by a hypothetical "TLR2-binding signature" within the sequence of PorB surface-exposed loops. Consistent with the effect of purified PorB in vitro, a chimeric N. meningitidis strain expressing N. lactamica PorB induces lower levels of interleukin 8 (IL-8) secretion than wild-type N. meningitidis, suggesting a role for PorB in induction of host cell activation by whole bacteria.

Identification of innate immunity elicitors using molecular signatures of natural selection

The innate immune system is an ancient and broad-spectrum defense system found in all eukaryotes. The detection of microbial elicitors results in the up-regulation of defense-related genes and the elicitation of inflammatory and apoptotic responses. These innate immune responses are the front-line barrier against disease because they collectively suppress the growth of the vast majority of invading microbes. Despite their critical role, we know remarkably little about the diversity of immune elicitors. To address this paucity, we reasoned that hosts are more likely to evolve recognition to "core" pathogen proteins under strong negative selection for the maintenance of essential cellular functions, whereas repeated exposure to host-defense responses will impose strong positive selective pressure for elicitor diversification to avoid host recognition. Therefore, we hypothesized that novel bacterial elicitors can be identified through these opposing forces of natural selection. We tested this hypothesis by examining the genomes of six bacterial phytopathogens and identifying 56 candidate elicitors that have an excess of positively selected residues in a background of strong negative selection. We show that these positively selected residues are atypically clustered, similar to patterns seen in the few well-characterized elicitors. We then validated selected candidate elicitors by showing that they induce Arabidopsis thaliana innate immunity in functional (virulence suppression) and cellular (callose deposition) assays. These finding provide targets for the study of host-pathogen interactions and applied research into alternative antimicrobial treatments.

Population structure in the Neisseria, and the biological significance of fuzzy species

Phenotypic and genetic variation in bacteria can take bewilderingly complex forms even within a single genus. One of the most intriguing examples of this is the genus Neisseria, which comprises both pathogens and commensals colonizing a variety of body sites and host species, and causing a range of disease. Complex relatedness among both named species and previously identified lineages of Neisseria makes it challenging to study their evolution. Using the largest publicly available collection of bacterial sequence data in combination with a population genetic analysis and experiment, we probe the contribution of inter-species recombination to neisserial population structure, and specifically whether it is more common in some strains than others. We identify hybrid groups of strains containing sequences typical of more than one species. These groups of strains, typical of a fuzzy species, appear to have experienced elevated rates of inter-species recombination estimated by population genetic analysis and further supported by transformation experiments. In particular, strains of the pathogen Neisseria meningitidis in the fuzzy species boundary appear to follow a different lifestyle, which may have considerable biological implications concerning distribution of novel resistance elements and meningococcal vaccine development. Despite the strong evidence for negligible geographical barriers to gene flow within the population, exchange of genetic material still shows directionality among named species in a non-uniform manner.

Induction of immune responses by purified outer membrane protein complexes from Neisseria meningitidis

A broad-spectrum vaccine against disease caused by serogroup B of Neisseria meningitidis is still a challenge due to antigenic variability. In the present study outer membrane protein complexes and their components were analysed using non-denaturing 2D electrophoresis and identified using LC/MS-MS and MALDI-TOF. Outer membrane protein complexes were purified from both the wild-type strain H44/76 and their knock-out mutants lacking PorA, PorB, RmpM or FetA. The immune responses elicited by the whole outer membrane vesicles (OMV) and the purified complexes were analysed for bactericidal activity, antibody surface binding, antibody-mediated C3b/iC3b deposition, membrane attack complex (MAC) deposition and induction of opsonophagocytosis, both on the homologous and several heterologous strains. The main antigenic complexes found were homomeric, formed by the 60 kDa chaperonin (MSP63) or PorB, or heteromeric, formed by different combinations of PorA, PorB and/or RmpM. The lack of some of these proteins in the OMVs from the knock-out mutants did not affect significantly the immune responses analysed except MAC, which was significantly reduced in the anti-PorA- and anti-PorB- sera, and bactericidal activity, which was absent in the anti-PorA- serum. The sera against purified native complexes showed variable activities against the homologous strain, with greatest responses observed for anti-chaperonin and anti-PorA/PorB/RmpM sera. When tested against heterologous strains, the only anti-complex serum showing consistent responses was that against the 60 kDa chaperonin. The comparison of the responses elicited by the different sera suggests an important role of conformational epitopes, present only in native complexes, in the induction of more effective responses against N. meningitidis.

Trends of the major porin gene (ompF) evolution: insight from the genus Yersinia

OmpF is one of the major general porins of Enterobacteriaceae that belongs to the first line of bacterial defense and interactions with the biotic as well as abiotic environments. Porins are surface exposed and their structures strongly reflect the history of multiple interactions with the environmental challenges. Unfortunately, little is known on diversity of porin genes of Enterobacteriaceae and the genus Yersinia especially. We analyzed the sequences of the ompF gene from 73 Yersinia strains covering 14 known species. The phylogenetic analysis placed most of the Yersinia strains in the same line assigned by 16S rDNA-gyrB tree. Very high congruence in the tree topologies was observed for Y. enterocolitica, Y. kristensenii, Y. ruckeri, indicating that intragenic recombination in these species had no effect on the ompF gene. A significant level of intra- and interspecies recombination was found for Y. aleksiciae, Y. intermedia and Y. mollaretii. Our analysis shows that the ompF gene of Yersinia has evolved with nonrandom mutational rate under purifying selection. However, several surface loops in the OmpF porin contain positively selected sites, which very likely reflect adaptive diversification Yersinia to their ecological niches. To our knowledge, this is a first investigation of diversity of the porin gene covering the whole genus of the family Enterobacteriaceae. This study demonstrates that recombination and positive selection both contribute to evolution of ompF, but the relative contribution of these evolutionary forces are different among Yersinia species.

Interplay of recombination and selection in the genomes of Chlamydia trachomatis

Background

Chlamydia trachomatis is an obligate intracellular bacterial parasite, which causes several severe and debilitating diseases in humans. This study uses comparative genomic analyses of 12 complete published C. trachomatis genomes to assess the contribution of recombination and selection in this pathogen and to understand the major evolutionary forces acting on the genome of this bacterium.

Results

The conserved core genes of C. trachomatis are a large proportion of the pan-genome: we identified 836 core genes in C. trachomatis out of a range of 874-927 total genes in each genome. The ratio of recombination events compared to mutation (ρ/θ) was 0.07 based on ancestral reconstructions using the ClonalFrame tool, but recombination had a significant effect on genetic diversification (r/m = 0.71). The distance-dependent decay of linkage disequilibrium also indicated that C. trachomatis populations behaved intermediately between sexual and clonal extremes. Fifty-five genes were identified as having a history of recombination and 92 were under positive selection based on statistical tests. Twenty-three genes showed evidence of being under both positive selection and recombination, which included genes with a known role in virulence and pathogencity (e.g., ompA, pmps, tarp). Analysis of inter-clade recombination flux indicated non-uniform currents of recombination between clades, which suggests the possibility of spatial population structure in C. trachomatis infections.

Conclusions

C. trachomatis is the archetype of a bacterial species where recombination is relatively frequent yet gene gains by horizontal gene transfer (HGT) and losses (by deletion) are rare. Gene conversion occurs at sites across the whole C. trachomatis genome but may be more often fixed in genes that are under diversifying selection. Furthermore, genome sequencing will reveal patterns of serotype specific gene exchange and selection that will generate important research questions for understanding C. trachomatis pathogenesis.

Reviewers

This article was reviewed by Dr. Jeremy Selengut, Dr. Lee S. Katz (nominated by Dr. I. King Jordan) and Dr. Arcady Mushegian.

Variation in the Neisseria meningitidis FadL-like protein: an evolutionary model for a relatively low-abundance surface antigen

The molecular diversity of a novel Neisseria meningitidis antigen, encoded by the ORF NMB0088 of MC58 (FadL-like protein), was assessed in a panel of 64 diverse meningococcal strains. The panel consisted of strains belonging to different serogroups, serotypes, serosubtypes and MLST sequence types, of different clinical sources, years and countries of isolation. Based on the sequence variability of the protein, the FadL-like protein has been divided into four variant groups in this species. Antigen variants were associated with specific serogroups and MLST clonal complexes. Maximum-likelihood analyses were used to determine the relationships among sequences and to compare the selection pressures acting on the encoded protein. Furthermore, a model of population genetics and molecular evolution was used to detect natural selection in DNA sequences using the non-synonymous : synonymous substitution (d(N) : d(S)) ratio. The meningococcal sequences were also compared with those of the related surface protein in non-pathogenic commensal Neisseria species to investigate potential horizontal gene transfer. The N. meningitidis fadL gene was subject to only weak positive selection pressure and was less diverse than meningococcal major outer-membrane proteins. The majority of the variability in fadL was due to recombination among existing alleles from the same or related species that resulted in a discrete mosaic structure in the meningococcal population. In general, the population structuring observed based on the FadL-like membrane protein indicates that it is under intermediate immune selection. However, the emergence of a new subvariant within the hyperinvasive lineages demonstrates the phenotypic adaptability of N. meningitidis, probably in response to selective pressure.

Evolution of bacterial genes: evidences of positive Darwinian selection and fixation of base substitutions in virulence genes of Helicobacter pylori

Gene diversity in Helicobacter pylori from different origins results in a phylogeographic differentiation, and this genetic variation among populations might be driven by random drift or by selective forces. However, only the selective forces would contribute to adaptation of the bacteria to the physiology and environment of its local host and to its association with gastroduodenal diseases. We studied evolutionary forces acting on variable regions of virulence genes cagA, babA and oipA, which present geographic differences among H. pylori strains from different human groups. Gene sequences in H. pylori strains from Asia, Europe and America were analysed using state of the art analytical methods like the Maximum Likelihood method. The rate and nature of polymorphisms in these virulence genes were also compared among populations using the AMOVA and McDonald-Kreitman tests. We found strong and significant positive selection acting on variable regions of cagA, babA and oipA. We found in cagA from Asian strains regions under positive selection, which localised in amino acid sites defining the Asian fingerprint for this gene and in sites with important biological activity. Different evolutionary forces are acting on the variable region of virulence genes; they partly explain the source of genetic diversity and the differences in risk for gastroduodenal diseases among different human populations.

Variation of the factor H-binding protein of Neisseria meningitidis

There is currently no comprehensive meningococcal vaccine, due to difficulties in immunizing against organisms expressing serogroup B capsules. To address this problem, subcapsular antigens, particularly the outer-membrane proteins (OMPs), are being investigated as candidate vaccine components. If immunogenic, however, such antigens are often antigenically variable, and knowledge of the extent and structuring of this diversity is an essential part of vaccine formulation. Factor H-binding protein (fHbp) is one such protein and is included in two vaccines under development. A survey of the diversity of the fHbp gene and the encoded protein in a representative sample of meningococcal isolates confirmed that variability in this protein is structured into two or three major groups, each with a substantial number of alleles that have some association with meningococcal clonal complexes and serogroups. A unified nomenclature scheme was devised to catalogue this diversity. Analysis of recombination and selection on the allele sequences demonstrated that parts of the gene are subject to positive selection, consistent with immune selection on the protein generating antigenic variation, particularly in the C-terminal region of the peptide sequence. The highest levels of selection were observed in regions corresponding to epitopes recognized by previously described bactericidal monoclonal antibodies.

Genome wide evolutionary analyses reveal serotype specific patterns of positive selection in selected Salmonella serotypes

Background

The bacterium Salmonella enterica includes a diversity of serotypes that cause disease in humans and different animal species. Some Salmonella serotypes show a broad host range, some are host restricted and exclusively associated with one particular host, and some are associated with one particular host species, but able to cause disease in other host species and are thus considered "host adapted". Five Salmonella genome sequences, representing a broad host range serotype (Typhimurium), two host restricted serotypes (Typhi [two genomes] and Paratyphi) and one host adapted serotype (Choleraesuis) were used to identify core genome genes that show evidence for recombination and positive selection.

Results

Overall, 3323 orthologous genes were identified in all 5 Salmonella genomes analyzed. Use of four different methods to assess homologous recombination identified 270 genes that showed evidence for recombination with at least one of these methods (false discovery rate [FDR] <10%). After exclusion of genes with evidence for recombination, site and branch specific models identified 41 genes as showing evidence for positive selection (FDR <20%), including a number of genes with confirmed or likely roles in virulence and ompC, a gene encoding an outer membrane protein, which has also been found to be under positive selection in other bacteria. A total of 8, 16, 7, and 5 genes showed evidence for positive selection in Choleraesuis, Typhi, Typhimurium, and Paratyphi branch analyses, respectively. Sequencing and evolutionary analyses of four genes in an additional 42 isolates representing 23 serotypes confirmed branch specific positive selection and recombination patterns.

Conclusion

Our data show that, among the four serotypes analyzed, (i) less than 10% of Salmonella genes in the core genome show evidence for homologous recombination, (ii) a number of Salmonella genes are under positive selection, including genes that appear to contribute to virulence, and (iii) branch specific positive selection contributes to the evolution of host restricted Salmonella serotypes.

Identification of new meningococcal serogroup B surface antigens through a systematic analysis of neisserial genomes

The difficulty of inducing an effective immune response against the Neisseria meningitidis serogroup B capsular polysaccharide has lead to the search for vaccines for this serogroup based on outer membrane proteins. The availability of the first meningococcal genome (MC58 strain) allowed the expansion of high-throughput methods to explore the protein profile displayed by N. meningitidis. By combining a pan-genome analysis with an extensive experimental validation to identify new potential vaccine candidates, genes coding for antigens likely to be exposed on the surface of the meningococcus were selected after a multistep comparative analysis of entire Neisseria genomes. Eleven novel putative ORF annotations were reported for serogroup B strain MC58. Furthermore, a total of 20 new predicted potential pan-neisserial vaccine candidates were produced as recombinant proteins and evaluated using immunological assays. Potential vaccine candidate coding genes were PCR-amplified from a panel of representative strains and their variability analyzed using maximum likelihood approaches for detecting positive selection. Finally, five proteins all capable of inducing a functional antibody response vs N. meningitidis strain CU385 were identified as new attractive vaccine candidates: NMB0606 a potential YajC orthologue, NMB0928 the neisserial NlpB (BamC), NMB0873 a LolB orthologue, NMB1163 a protein belonging to a curli-like assembly machinery, and NMB0938 (a neisserial specific antigen) with evidence of positive selection appreciated for NMB0928. The new set of vaccine candidates and the novel proposed functions will open a new wave of research in the search for the elusive neisserial vaccine.

Investigation on the effect of immune selection on resistance to bactericidal antibodies to group B meningococci in vitro

The induction of resistance by immune selective pressure to bactericidal antibodies from humans immunized with Novartis recombinant meningococcal group B vaccines was assessed. Serum bactericidal antibody titers against selected bacteria were within assay variability through a selection event frequency of 1 in 10(-5). No change in antigen expression was observed by Western blotting.

Meningococcal protein antigens and vaccines

The development of a comprehensive vaccine against meningococcal disease has been challenging. Recent developments in molecular genetics have provided both explanations for these challenges and possible solutions. Since genome sequence data became available there has been a marked increase in number of protein antigens that have been suggested as prospective vaccine components. This review catalogues the proposed vaccine candidates and examines the evidence for their inclusion in potential protein vaccine formulations.

Diversifying and stabilizing selection of sialidase and N-acetylneuraminate catabolism in Mycoplasma synoviae

Sialidase activity varies widely among strains and tends to correlate with strain virulence in the avian pathogen Mycoplasma synoviae. To characterize the forms of selection acting on enzymes required for sialic acid scavenging and catabolism, the ratios of nonsynonymous (K(a)) to synonymous (K(s)) mutation frequency were calculated for codons in the sialidase gene of 16 strains of M. synoviae and for its nearly identical homolog in four strains of Mycoplasma gallisepticum. The K(a)/K(s) (omega) values for the linked genes required for nutritive N-acetylneuraminate catabolism (nanA, nagC, nanE, nagA, and nagB) from nine strains of M. synoviae were also determined. To provide context, omega was determined for all corresponding genes of 26 strains of Clostridium perfringens and Streptococcus pneumoniae. Bayesian models of sequence evolution showed that only the sialidase of M. synoviae was under significant (P < 0.001) diversifying selection, while the M. synoviae genes for N-acetylneuraminate catabolism and all genes examined from M. gallisepticum, C. perfringens, and S. pneumoniae were under neutral to stabilizing selection. Diversifying selection acting on the sialidase of M. synoviae, but not on the sialidase of M. gallisepticum or the sialidases or other enzymes essential for sialic acid scavenging in other Firmicutes, is evidence that variation in specific activity of the enzyme is perpetuated by a nonnutritive function in M. synoviae that is influenced by the genomic context of the organism.

Distribution and genetic variability of three vaccine components in a panel of strains representative of the diversity of serogroup B meningococcus

With the aim of studying the molecular diversity of the antigens of a new recombinant vaccine against meningococcus serogroup B, the three genes coding for the main vaccine components GNA (Genome-derived Neisseria Antigen) 1870 (fHbp, factor H Binding Protein), GNA1994 (NadA, Neisseria adhesin A) and GNA2132 were sequenced in a panel of 85 strains collected worldwide and selected as representative of the serogroup B meningococcal diversity. No correlations were found between vaccine antigen variability and serogroup, geographic area and year of isolation. Although a relevant clustering was found with MLST clonal complexes, each showing an almost specific antigen variant repertoire, the prediction of the antigen assortment was not possible on the basis of MLST alone. Therefore, classification of meningococcus on the basis of MLST only is not sufficient to predict vaccine antigens diversity. Sequencing each gene in the different strains will be important to evaluate antigen conservation and assortment and to allow a future prediction of potential vaccine coverage.

Cloning of immunodominant membrane protein genes of phytoplasmas and their in planta expression

Phytoplasmas are plant pathogenic bacteria that cause devastating yield losses in diverse crops worldwide. Although the understanding of the pathogen biology is important in agriculture, the inability to culture phytoplasmas has hindered their full characterization. Previous studies demonstrated that immunodominant membrane proteins could be classified into three types, immunodominant membrane protein (Imp), immunodominant membrane protein A (IdpA), and antigenic membrane protein (Amp), and they are nonhomologous to each other. Here, cloning and sequencing of imp-containing genomic fragments were performed for several groups of phytoplasma including the aster yellows and rice yellow dwarf groups, for which an imp sequence has not previously been reported. Sequence comparison analysis revealed that Imps are highly variable among phytoplasmas, and clear positive selection was observed in several Imps, suggesting that Imp has important roles in host-phytoplasma interactions. As onion yellows (OY) phytoplasma was known to have Amp as the immunodominant membrane protein, the protein accumulation level of Imp in planta was measured compared with that of Amp. The resulting accumulation of Imp was calculated as approximately one-tenth that of Amp, being consistent with the immunodominant property of Amp in OY. It is suggested that an ancestral type of immunodominant membrane protein could be Imp, and subsequently the expression level of Amp or IdpA is increased in several phytoplasma groups.

Pilus operon evolution in Streptococcus pneumoniae is driven by positive selection and recombination

BACKGROUND

The evolution of bacterial organelles involved in host-pathogen interactions is subject to intense and competing selective pressures due to the need to maintain function while escaping the host immune response. To characterize the interplay of these forces in an important pathogen, we sequenced the rlrA islet, a chromosomal region encoding for a pilus-like structure involved in adherence to lung epithelial cells in vitro and in colonization in a murine model of infection, in 44 clinical isolates of Streptococcus pneumoniae.

RESULTS

We found that the rrgA and rrgB genes, encoding the main structural components of the pilus, are under the action of positive selection. In contrast, the rrgC gene, coding for a component present in low quantities in the assembled pilus, and the srtB, srtC and srtD genes, coding for three sortase enzymes essential for pilus assembly but probably not directly exposed to the host immune system, show no evidence of positive selection. We found several events of homologous recombination in the region containing these genes, identifying 4 major recombination hotspots. An analysis of the most recent recombination events shows a high level of mosaicism of the region coding for the rrgC, srtB, srtC and srtD genes.

CONCLUSIONS

In the rlrA islet, the genes coding for proteins directly exposed to the host immune response are under the action of positive selection, and exist in distinct forms in the population of circulating strains. The genes coding for proteins not directly exposed on the surface of the bacterial cell are more conserved probably due to the homogenizing effect of recombination.

Accelerated evolutionary rates in tropical and oceanic parmelioid lichens (Ascomycota)

Background

The rate of nucleotide substitutions is not constant across the Tree of Life, and departures from a molecular clock have been commonly reported. Within parmelioid lichens, the largest group of macrolichens, large discrepancies in branch lengths between clades were found in previous studies. Using an extended taxon sampling, we test for presence of significant rate discrepancies within and between these clades and test our a priori hypothesis that such rate discrepancies may be explained by shifts in moisture regime or other environmental conditions.

Results

In this paper, the first statistical evidence for accelerated evolutionary rate in lichenized ascomycetes is presented. Our results give clear evidence for a faster rate of evolution in two Hypotrachyna clades that includes species occurring in tropical and oceanic habitats in comparison with clades consisting of species occurring in semi-arid and temperate habitats. Further we explore potential links between evolutionary rates and shifts in habitat by comparing alternative Ornstein-Uhlenbeck models.

Conclusion

Although there was only weak support for a shift at the base of a second tropical clade, where the observed nucleotide substitution rate is high, overall support for a shift in environmental conditions at cladogenesis is very strong. This suggests that speciation in some lichen clades has proceeded by dispersal into a novel environment, followed by radiation within that environment. We found moderate support for a shift in moisture regime at the base of one tropical clade and a clade occurring in semi-arid regions and a shift in minimum temperature at the base of a boreal-temperate clade.

Evolutionary analysis of endopolygalacturonase-encoding genes of Botrytis cinerea

Sequence analysis of five of the six endopolygalacturonase-encoding genes (Bcpg1, Bcpg2, Bcpg3, Bcpg4, Bcpg5) from 32 strains of Botrytis cinerea showed marked gene to gene differences in the amount of among-strains diversity. Bcpg4 was almost invariable in all strains; Bcpg3 and Bcpg5 showed a moderate variability, similar to that of non-pathogenicity-associated genes examined in other studies. Conversely, Bcpg1 and Bcpg2 were highly variable and were shown to be under positive selection based on the McDonald-Kreitman test and likelihood ratio test. The evolution of the five endopolygalacturonase genes is explained by their different ecophysiological role. Diversification and balancing selection, as detected in Bcpg1 and Bcpg2, can be used by the pathogen to escape recognition by the host and delay plant reaction in the early phases of infection. The analysis of the polymorphisms and the location of the sites with high probability of being positively selected highlighted the relevance of variability of the BcPG1 and BcPG2 proteins at their C-terminal end. By contrast, the absence of variability in Bcpg4 suggests that the efficiency of the product of this gene is critical for B. cinerea growth in late phases of infection or during intraspecific competition, thus markedly affecting strain fitness.

Genome-wide analyses reveal lineage specific contributions of positive selection and recombination to the evolution of Listeria monocytogenes

Background

The genus Listeria includes two closely related pathogenic and non-pathogenic species, L. monocytogenes and L. innocua. L. monocytogenes is an opportunistic human foodborne and animal pathogen that includes two common lineages. While lineage I is more commonly found among human listeriosis cases, lineage II appears to be overrepresented among isolates from foods and environmental sources. This study used the genome sequences for one L. innocua strain and four L. monocytogenes strains representing lineages I and II, to characterize the contributions of positive selection and recombination to the evolution of the L. innocua/L. monocytogenes core genome.

Results

Among the 2267 genes in the L. monocytogenes/L. innocua core genome, 1097 genes showed evidence for recombination and 36 genes showed evidence for positive selection. Positive selection was strongly associated with recombination. Specifically, 29 of the 36 genes under positive selection also showed evidence for recombination. Recombination was more common among isolates in lineage II than lineage I; this trend was confirmed by sequencing five genes in a larger isolate set. Positive selection was more abundant in the ancestral branch of lineage II (20 genes) as compared to the ancestral branch of lineage I (9 genes). Additional genes under positive selection were identified in the branch separating the two species; for this branch, genes in the role category "Cell wall and membrane biogenesis" were significantly more likely to have evidence for positive selection. Positive selection of three genes was confirmed in a larger isolate set, which also revealed occurrence of multiple premature stop codons in one positively selected gene involved in flagellar motility (flaR).

Conclusion

While recombination and positive selection both contribute to evolution of L. monocytogenes, the relative contributions of these evolutionary forces seem to differ by L. monocytogenes lineages and appear to be more important in the evolution of lineage II, which seems to be found in a broader range of environments, as compared to the apparently more host adapted lineage I. Diversification of cell wall and membrane biogenesis and motility-related genes may play a particularly important role in the evolution of L. monocytogenes.

Variation in the Neisseria lactamica porin, and its relationship to meningococcal PorB

One potential vaccine strategy in the fight against meningococcal disease involves the exploitation of outer-membrane components of Neisseria lactamica, a commensal bacterium closely related to the meningococcus, Neisseria meningitidis. Although N. lactamica shares many surface structures with the meningococcus, little is known about the antigenic diversity of this commensal bacterium or the antigenic relationships between N. lactamica and N. meningitidis. Here, the N. lactamica porin protein (Por) was examined and compared to the related PorB antigens of N. meningitidis, to investigate potential involvement in anti-meningococcal immunity. Relationships among porin sequences were determined using distance-based methods and F(ST), and maximum-likelihood analyses were used to compare the selection pressures acting on the encoded proteins. These analyses demonstrated that the N. lactamica porin was less diverse than meningococcal PorB and although it was subject to positive selection, this was not as strong as the positive selection pressures acting on the meningococcal porin. In addition, the N. lactamica porin gene sequences and the protein sequences of the loop regions predicted to be exposed to the human immune system were dissimilar to the corresponding sequences in the meningococcus. This suggests that N. lactamica Por, contrary to previous suggestions, may have limited involvement in the development of natural immunity to meningococcal disease and might not be effective as a meningococcal vaccine component.

Phytoplasmas: bacteria that manipulate plants and insects

TAXONOMY

Superkingdom Prokaryota; Kingdom Monera; Domain Bacteria; Phylum Firmicutes (low-G+C, Gram-positive eubacteria); Class Mollicutes; Candidatus (Ca.) genus Phytoplasma.

HOST RANGE

Ca. Phytoplasma comprises approximately 30 distinct clades based on 16S rRNA gene sequence analyses of approximately 200 phytoplasmas. Phytoplasmas are mostly dependent on insect transmission for their spread and survival. The phytoplasma life cycle involves replication in insects and plants. They infect the insect but are phloem-limited in plants. Members of Ca. Phytoplasma asteris (16SrI group phytoplasmas) are found in 80 monocot and dicot plant species in most parts of the world. Experimentally, they can be transmitted by approximately 30, frequently polyphagous insect species, to 200 diverse plant species.

DISEASE SYMPTOMS

In plants, phytoplasmas induce symptoms that suggest interference with plant development. Typical symptoms include: witches' broom (clustering of branches) of developing tissues; phyllody (retrograde metamorphosis of the floral organs to the condition of leaves); virescence (green coloration of non-green flower parts); bolting (growth of elongated stalks); formation of bunchy fibrous secondary roots; reddening of leaves and stems; generalized yellowing, decline and stunting of plants; and phloem necrosis. Phytoplasmas can be pathogenic to some insect hosts, but generally do not negatively affect the fitness of their major insect vector(s). In fact, phytoplasmas can increase fecundity and survival of insect vectors, and may influence flight behaviour and plant host preference of their insect hosts.

DISEASE CONTROL

The most common practices are the spraying of various insecticides to control insect vectors, and removal of symptomatic plants. Phytoplasma-resistant cultivars are not available for the vast majority of affected crops.

Plasmodium falciparum antigenic variation. Mapping mosaic var gene sequences onto a network of shared, highly polymorphic sequence blocks

Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) is a potentially important family of immune targets, encoded by an extremely diverse gene family called var. Understanding of the genetic organization of var genes is hampered by sequence mosaicism that results from a long history of non-homologous recombination. Here we have used software designed to analyse social networks to visualize the relationships between large collections of short var sequences tags sampled from clinical parasite isolates. In this approach, two sequences are connected if they share one or more highly polymorphic sequence blocks. The results show that the majority of analysed sequences including several var-like sequences from the chimpanzee parasite Plasmodium reichenowi can be either directly or indirectly linked together in a single unbroken network. However, the network is highly structured and contains putative subgroups of recombining sequences. The major subgroup contains the previously described group A var genes, previously proposed to be genetically distinct. Another subgroup contains sequences found to be associated with rosetting, a parasite virulence phenotype. The mosaic structure of the sequences and their division into subgroups may reflect the conflicting problems of maximizing antigenic diversity and minimizing epitope sharing between variants while maintaining their host cell binding functions.

Molecular characterization of Neisseria meningitidis isolates using a resequencing DNA microarray

Neisseria meningitidis is a major cause of both meningitis and septicemia. Typically, isolates are characterized by using a combination of immunological phenotyping, using monoclonal and polyclonal antisera, and Sanger nucleotide sequencing of epitope-encoding variable regions, although these methods can be both time-consuming and limited by reagent availability. Herein, we describe and evaluate a novel microarray to define the porB and porA serotypes of N. meningitidis by the resequencing of variable regions in a single hybridization reaction. PCR products for each gene were amplified, pooled in equimolar concentrations, hybridized to the microarray, and analyzed using Affymetrix GeneChip DNA Analysis Software. Resequencing of the microarray data was then validated by comparison with sequencing data. Molecular profiles were generated for 50 isolates that were combinations of phenotypically typeable (ie, PorA and PorB) and non-typeable (PorB only) isolates. Microarray-generated profiles from isolates with a PorB phenotype were concordant with predicted profiles compared with a previously described typing scheme. In addition, 42% (8 of 19) of previously non-typeable samples were assigned a PorB type when tested using the microarray. The remaining isolates were novel types for which no typing antisera are currently available. The porA data were 97% concordant with Sanger nucleotide sequencing. These results suggest that that microarray resequencing may be a useful tool for the characterization of meningococci, particularly for those isolates that cannot be phenotyped, offering an alternative to conventional sequencing methods.

Recombination and positive selection contribute to evolution of Listeria monocytogenes inlA

The surface molecule InlA interacts with E-cadherin to promote invasion of Listeria monocytogenes into selected host cells. DNA sequencing of inlA for 40 L. monocytogenes isolates revealed 107 synonymous and 45 nonsynonymous substitutions. A frameshift mutation in a homopolymeric tract encoding part of the InlA signal peptide was identified in three lineage II isolates, which also showed reduced ability to invade human intestinal epithelial cells. Phylogenies showed clear separation of inlA sequences into lineages I and II. Thirteen inlA recombination events, predominantly involving lineage II strains as recipients (12 events), were detected and a number of amino acid residues were shown to be under positive selection. Four of the 45 non-synonymous changes were found to be under positive selection with posterior probabilities >95 %. Mapping of polymorphic and positively selected amino acid sites on the partial crystal structure for InlA showed that the internalin surface of the leucine-rich repeat (LRR) region that faces the InlA receptor E-cadherin does not include any polymorphic sites; all polymorphic and positively selected amino acids mapped to the outer face of the LRR region or to other InlA regions. The data show that (i) inlA is highly polymorphic and evolution of inlA involved a considerable number of recombination events in lineage II isolates; (ii) positive selection at specific amino acid sites appears to contribute to evolution of inlA, including fixation of recombinant events; and (iii) single-nucleotide deletions in a lineage II-specific 3' homopolymeric tract in inlA lead to complete loss of InlA or to production of truncated InlA, which conveys reduced invasiveness. In conclusion, inlA has a complex evolutionary history, which is consistent with L. monocytogenes' natural history as an environmental pathogen with broad host-range, including its adaptation to environments and hosts where different inlA alleles may provide a selective advantage or where inlA may not be required.

Haplotypic diversity in human CEACAM genes: effects on susceptibility to meningococcal disease

Adhesion between the opacity-associated adhesin (Opa) proteins of Neisseria meningitidis and human carcino-embryonic antigen cell adhesion molecule (CEACAM) proteins is an important stage in the pathogenesis of meningococcal disease, a globally important bacterial infection. Most disease is caused by a small number of meningococcal genotypes known as hyperinvasive lineages. As these are also carried asymptomatically, acquisition of them alone cannot explain why only some hosts develop meningococcal disease. Our aim was to determine whether genetic diversity in CEACAM is associated with susceptibility to meningococcal disease. Frequency distributions of alleles, genotypes and haplotypes were compared in four CEACAM genes in 384 case samples and 190 controls. Linkage disequilibrium among polymorphic sites, haplotype structures and relationships were also analysed. A number of polymorphisms were observed in CEACAM genes but the diversity of CEACAM1, to which most Opa proteins bind, was lower, and a small number of high-frequency haplotypes were detected. Dose-dependent associations of three CEACAM haplotypes with meningococcal disease were observed, with the effect of carrying these haplotypes amplified in homozygous individuals. Two haplotypes were protective while one haplotype in CEACAM6 was associated with a twofold increase in disease susceptibility. These data imply that human CEACAM may be one determinant of human susceptibility to meningococcal disease.

Helicobacter pylori evolution: lineage- specific adaptations in homologs of eukaryotic Sel1-like genes

Geographic partitioning is postulated to foster divergence of Helicobacter pylori populations as an adaptive response to local differences in predominant host physiology. H. pylori's ability to establish persistent infection despite host inflammatory responses likely involves active management of host defenses using bacterial proteins that may themselves be targets for adaptive evolution. Sequenced H. pylori genomes encode a family of eight or nine secreted proteins containing repeat motifs that are characteristic of the eukaryotic Sel1 regulatory protein, whereas the related Campylobacter and Wolinella genomes each contain only one or two such "Sel1-like repeat" (SLR) genes ("slr genes"). Signatures of positive selection (ratio of nonsynonymous to synonymous mutations, dN/dS = omega > 1) were evident in the evolutionary history of H. pylori slr gene family expansion. Sequence analysis of six of these slr genes (hp0160, hp0211, hp0235, hp0519, hp0628, and hp1117) from representative East Asian, European, and African H. pylori strains revealed that all but hp0628 had undergone positive selection, with different amino acids often selected in different regions. Most striking was a divergence of Japanese and Korean alleles of hp0519, with Japanese alleles having undergone particularly strong positive selection (omegaJ > 25), whereas alleles of other genes from these populations were intermingled. Homology-based structural modeling localized most residues under positive selection to SLR protein surfaces. Rapid evolution of certain slr genes in specific H. pylori lineages suggests a model of adaptive change driven by selection for fine-tuning of host responses, and facilitated by geographic isolation. Characterization of such local adaptations should help elucidate how H. pylori manages persistent infection, and potentially lead to interventions tailored to diverse human populations.

Detecting the form of selection in the outer membrane protein C of Enterobacter aerogenes strains and Salmonella species

The types of selective pressure operating on the outer membrane protein C (ompC) of Enterobacter aerogenes strains, the causative agent for nosocomial infections, and Salmonella sp., the hazardous pathogen are investigated using the maximum likelihood-based codon substitution models. Although the rate of amino acid replacement to the silent substitution (omega) across the entire codon sites of ompC of E. aerogenes (omega=0.3194) and Salmonella sp. (omega=0.2047) indicate that the gene is subjected to purifying selection (i.e. omega<1), approximately 3.7% of ompC codon sites in E. aerogenes (omega=21.52) are under the influence of positive Darwinian selection (i.e. omega>1). Such contrast in the intensity of selective pressures in both pathogens could be associated with the differential response to the adverse environmental changes. In E. aerogenes, majority of the positively selected sites are located in the hypervariable cell-surface-exposed domains whereas the trans-membrane domains are functionally highly constrained.

Positively selected codons in immune-exposed loops of the vaccine candidate OMP-P1 of Haemophilus influenzae

The high levels of variation in surface epitopes can be considered as an evolutionary hallmark of immune selection. New computational tools enable analysis of this variation by identifying codons that exhibit high rates of amino acid changes relative to the synonymous substitution rate. In the outer membrane protein P1 of Haemophilus influenzae, a vaccine candidate for nontypeable strains, we identified four codons with this attribute in domains that did not correspond to known or assumed B- and T-cell epitopes of OMP-P1. These codons flank hypervariable domains and do not appear to be false positives as judged from parsimony and maximum likelihood analyses. Some closely spaced positively selected codons have been previously considered part of a transmembrane domain, which would render this region unsuited for inclusion in a vaccine. Secondary structure analysis, three-dimensional structural database searches, and homology modeling using FadL of E. coli as a structural homologue, however, revealed that all positively selected codons are located in or near extracellular looping domains. The spacing and level of diversity of these positively selected and exposed codons in OMP-P1 suggest that vaccine targets based on these and conserved flanking residues may provide broad coverage in H. influenzae.

Molecular typing of meningococci: recommendations for target choice and nomenclature

The diversity and dynamics of Neisseria meningitidis populations generate a requirement for high resolution, comprehensive, and portable typing schemes for meningococcal disease surveillance. Molecular approaches, specifically DNA amplification and sequencing, are the methods of choice for various reasons, including: their generic nature and portability, comprehensive coverage, and ready implementation to culture negative clinical specimens. The following target genes are recommended: (1) the variable regions of the antigen-encoding genes porA and fetA and, if additional resolution is required, the porB gene for rapid investigation of disease outbreaks and investigating the distribution of antigenic variants; (2) the seven multilocus sequence typing loci-these data are essential for the most effective national, and international management of meningococcal disease, as well as being invaluable in studies of meningococcal population biology and evolution. These targets have been employed extensively in reference laboratories throughout the world and validated protocols have been published. It is further recommended that a modified nomenclature be adopted of the form: serogroup: PorA type: FetA type: sequence type (clonal complex), thus: B: P1.19,15: F5-1: ST-33 (cc32).

Sequencing of the porB gene: a step toward a true characterization of Neisseria meningitidis

Variations in class 2/3 (PorB) proteins form the basis for meningococcal serotyping. Antibodies against these proteins are bactericidal, making serotyping results useful not only for epidemiological surveillance of meningococcal disease but also for identifying potential vaccine components. A total of 20 to 60% of meningococcal B and C isolates from any given population are nontypeable (NT) using a panel of monoclonal antibodies. To analyze the mechanisms responsible for the nonserotypeability characteristic in Neisseria meningitidis, we (i) established the nucleotide sequences of porB gene in 146 meningococcal strains (95 not recognized by the serotyping panel), (ii) identified 18 new allelic variants of the porB gene, (iii) correlated allelic variants with serotypes, (iv) suggest the nontypeability characteristic in those 95 NT strains, and (v) reject the possibility of variation in the levels of PorB expression.

Proteomic analysis of a meningococcal outer membrane vesicle vaccine prepared from the group B strain NZ98/254

In the absence of a suitable carbohydrate-based vaccine, outer membrane vesicle (OMV) vaccines have been used to disrupt outbreaks of serogroup B meningococcal disease for more than 20 years. Proteomic technology provides physical methods with the potential to assess the composition and consistency of these complex vaccines. 2-DE, combined with MS, were used to generate a proteome map of an OMV vaccine, developed to disrupt a long-running outbreak of group B disease in New Zealand. Seventy four spots from the protein map were identified including the outer membrane protein (OMP) antigens: PorA, PorB, RmpM and OpcA. Protein identification indicates that, in addition to OMPs, OMV vaccines contain periplasmic, membrane-associated and cytoplasmic proteins. 2-D-DIGE technology highlighted differences between preclinical development batches of vaccines from two different manufacturers.

Who ate whom? Adaptive Helicobacter genomic changes that accompanied a host jump from early humans to large felines

Helicobacter pylori infection of humans is so old that its population genetic structure reflects that of ancient human migrations. A closely related species, Helicobacter acinonychis, is specific for large felines, including cheetahs, lions, and tigers, whereas hosts more closely related to humans harbor more distantly related Helicobacter species. This observation suggests a jump between host species. But who ate whom and when did it happen? In order to resolve this question, we determined the genomic sequence of H. acinonychis strain Sheeba and compared it to genomes from H. pylori. The conserved core genes between the genomes are so similar that the host jump probably occurred within the last 200,000 (range 50,000-400,000) years. However, the Sheeba genome also possesses unique features that indicate the direction of the host jump, namely from early humans to cats. Sheeba possesses an unusually large number of highly fragmented genes, many encoding outer membrane proteins, which may have been destroyed in order to bypass deleterious responses from the feline host immune system. In addition, the few Sheeba-specific genes that were found include a cluster of genes encoding sialylation of the bacterial cell surface carbohydrates, which were imported by horizontal genetic exchange and might also help to evade host immune defenses. These results provide a genomic basis for elucidating molecular events that allow bacteria to adapt to novel animal hosts.

Who ate whom? Adaptive Helicobacter genomic changes that accompanied a host jump from early humans to large felines

Helicobacter pylori infection of humans is so old that its population genetic structure reflects that of ancient human migrations. A closely related species, Helicobacter acinonychis, is specific for large felines, including cheetahs, lions, and tigers, whereas hosts more closely related to humans harbor more distantly related Helicobacter species. This observation suggests a jump between host species. But who ate whom and when did it happen? In order to resolve this question, we determined the genomic sequence of H. acinonychis strain Sheeba and compared it to genomes from H. pylori. The conserved core genes between the genomes are so similar that the host jump probably occurred within the last 200,000 (range 50,000-400,000) years. However, the Sheeba genome also possesses unique features that indicate the direction of the host jump, namely from early humans to cats. Sheeba possesses an unusually large number of highly fragmented genes, many encoding outer membrane proteins, which may have been destroyed in order to bypass deleterious responses from the feline host immune system. In addition, the few Sheeba-specific genes that were found include a cluster of genes encoding sialylation of the bacterial cell surface carbohydrates, which were imported by horizontal genetic exchange and might also help to evade host immune defenses. These results provide a genomic basis for elucidating molecular events that allow bacteria to adapt to novel animal hosts.

Positive selection acting on a surface membrane protein of the plant-pathogenic phytoplasmas

Phytoplasmas are plant-pathogenic bacteria that cause numerous diseases. This study shows a strong positive selection on the phytoplasma antigenic membrane protein (Amp). The ratio of nonsynonymous to synonymous substitutions was >1 with all the methods we tested. The clear positive selections imply an important biological role for Amp in host-bacterium interactions.

Rapid Adaptive Evolution of the Tumor Suppressor Gene Pten in an Insect Lineage

The Pten gene was initially identified in humans as a tumor suppressor. It has since been shown to play important roles in the control of cell size, cell motility, apoptosis, and organ size, and it has also been implicated in aging. Pten is highly conserved among organisms as diverse as nematodes, insects, and vertebrates. In contrast, a phylogenetic analysis by maximum likelihood of a 133-amino acid region showed an average nonsynonymous-to-synonymous rate ratio of 10.4 for Pten in the lineage leading to parasitoid wasps of the Nasonia genus, indicating very strong positive selection. A previous study identified Pten as a potential QTL candidate gene for differences in male wing size in Nasonia. Most of the amino acid replacements that occurred in the Nasonia lineage cluster in a small region of the protein surface, suggesting that they might be involved in an interaction between Pten and another protein. The phenotypic changes due to Pten are not yet known, although it is not associated with known differences in male wing size. Introgression of Pten from one species to another does affect longevity, but a causal relationship is not established.

Positive selection of the Hrp pilin HrpE of the plant pathogen Xanthomonas

The plant-pathogenic bacterium Xanthomonas campestris pv. vesicatoria possesses a type III secretion (TTS) system which is encoded by the 23-kb hrp (hypersensitive response and pathogenicity) gene cluster. The TTS system is necessary for pathogenicity in susceptible hosts and induction of the hypersensitive response in resistant plants. At the cell surface, the TTS system is associated with an extracellular filamentous structure, the Hrp pilus, which serves as a conduit for the transfer of bacterial proteins into the plant cell cytosol. The major pilus component, the HrpE pilin, is unique to xanthomonads. Previous work showed that HrpE contains two regions: a hypervariable surface-exposed domain, including the N-terminal secretion signal, and a C-terminal polymerization domain. In this study, the evolutionary rate of the hrpE gene was analyzed. Twenty-one alleles were cloned, sequenced, and compared with five known hrpE alleles. The ratio of synonymous (K(s)) and nonsynonymous (K(a)) substitution rates shows that parts of the HrpE N terminus are subjected to positive selection and the C terminus is subjected to purifying selection. The trade-off between positive and purifying selection at the very-N terminus allowed us to ascertain the amphipathic alpha-helical nature of the TTS signal. This is the first report of a surface structure from a plant-pathogenic bacterium that evolved under the constraint of positive selection and hints to the evolutionary adaptation of this extracellular appendage to avoid recognition by the plant defense surveillance system.

Estimating diversifying selection and functional constraint in the presence of recombination

Models of molecular evolution that incorporate the ratio of nonsynonymous to synonymous polymorphism (dN/dS ratio) as a parameter can be used to identify sites that are under diversifying selection or functional constraint in a sample of gene sequences. However, when there has been recombination in the evolutionary history of the sequences, reconstructing a single phylogenetic tree is not appropriate, and inference based on a single tree can give misleading results. In the presence of high levels of recombination, the identification of sites experiencing diversifying selection can suffer from a false-positive rate as high as 90%. We present a model that uses a population genetics approximation to the coalescent with recombination and use reversible-jump MCMC to perform Bayesian inference on both the dN/dS ratio and the recombination rate, allowing each to vary along the sequence. We demonstrate that the method has the power to detect variation in the dN/dS ratio and the recombination rate and does not suffer from a high false-positive rate. We use the method to analyze the porB gene of Neisseria meningitidis and verify the inferences using prior sensitivity analysis and model criticism techniques.