Population genetics of Neisseria gonorrhoeae in a high-prevalence community using a hypervariable outer membrane porB and 13 slowly evolving housekeeping genes

Baltimore, Md., is an urban community with a high prevalence of Neisseria gonorrhoeae. Due to partially protective immune responses, introduction of new strains from other host populations, and exposure of N. gonorrhoeae to antibiotics, the phenotypic and genotypic characteristics of the circulating strains can fluctuate over time. Understanding the overall genetic diversity and population structure of N. gonorrhoeae is essential for informing public health interventions to eliminate this pathogen. We studied gonococci population genetics in Baltimore by analyzing a hypervariable and strongly selected outer membrane porB gene and 13 slowly evolving and presumably neutral housekeeping genes (abcZ, adk, aroE, fumC, gdh, glnA, gnd, pdhC, pgm, pilA, ppk, pyrD, and serC) in 204 isolates collected in 1991, 1996, and 2001 from male and female patients of two public sexually transmitted diseases clinics. Genetic diversity (), recombination (C), growth (g), population structure, and adaptive selection under codon-substitution and amino acid property models were estimated and compared between these two gene classes. Estimates of the F(ST) fixation index and the chi(2) test of sequence absolute frequencies revealed significant temporal substructuring for both gene types. Baltimore's N. gonorrhoeae populations have increased since 1991 as indicated by consistent positive values of g. Female patients showed similar or lower levels of and C than male patients. Within the MLST housekeeping genes, levels of and C ranged from 0.001-0.013 and 0.000-0.018, respectively. Overall recombination seems to be the dominant force driving evolution in these populations. All loci showed amino acid sites and physicochemical properties under adaptive (or positive-destabilizing) selection, rejecting the generally assumed hypothesis of stabilizing selection for these MLST genes. Within the porB gene, protein I B showed higher and C values than protein I A. Directional positive selection possibly mediated by the immune system operates to a significant extent in the protein I sequences, as indicated by the distribution of the positively selected sites in the surface-exposed loops. Thirteen amino acid physicochemical properties seem to drive protein evolution of the PI porins in N. gonorrhoeae.     

Sequence analysis and relationships between meningococcal class 3 serotype proteins and other porins from pathogenic and non-pathogenic Neisseria species

The presence of highly conserved regions within previously determined porin gene sequences from Neisseria meningitidis and Neisseria gonorrhoeae permitted the construction of oligonucleotide primers for PCR amplification of other neisserial porin genes. Although two separate porin genes (porA and porB) are present in N. meningitidis only a single fragment, corresponding to porB, could be amplified from this species. The amplified porB genes from four different meningococcal serotypes, which express the class 3 outer membrane protein, were sequenced. Amplified fragments corresponding to porin genes from N. lactamica and N. sicca were also sequenced. In common with the known neisserial porins, models of the organisation of the predicted proteins indicated trans-membrane structures with eight surface exposed loops. In the meningococcal class 3 proteins the main regions of sequence variation, which must be responsible for serotype specificity, were located on loops 5 and 7. A phylogenetic analysis of the family of porins from the Neisseria confirmed the close relationship of the meningococcal class 3 protein with the gonococcal PIA protein, while the gonococcal PIB protein was shown to be closely related to the N. lactamica porin. The close relationship seen between porins of the pathogenic and non-pathogenic Neisseriae identified no obvious virulence-associated regions in the proteins, but did suggest that the current nomenclature for neisserial porin genes may need reviewing.     

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.     

Porin-mediated antibiotic resistance in Neisseria gonorrhoeae: ion, solute, and antibiotic permeation through PIB proteins with penB mutations

Neisseria gonorrhoeae has two porins, PIA and PIB, whose genes (porA and porB, respectively) are alleles of a single por locus. We recently demonstrated that penB mutations at positions 120 and 121 in PIB, which are presumed to reside in loop 3 that forms the pore constriction zone, confer intermediate-level resistance to penicillin and tetracycline (M. Olesky, M. Hobbs, and R. A. Nicholas, Antimicrob. Agents Chemother. 46:2811-2820, 2002). In the present study, we investigated the electrophysiological properties as well as solute and antibiotic permeation rates of recombinant PIB proteins containing penB mutations (G120K, G120D/A121D, G120P/A121P, and G120R/A121H). In planar lipid bilayers, the predominant conducting state of each porin variant was 30 to 40% of the wild type, even though the anion selectivity and maximum channel conductance of each PIB variant was similar to that of the wild type. Liposome-swelling experiments revealed no significant differences in the permeation of sugars or beta-lactam antibiotics through the wild type or PIB variants. Although these results are seemingly contradictory with the ability of these variants to increase antibiotic resistance, they are consistent with MIC data showing that these porin mutations confer resistance only in strains containing an mtrR mutation, which increases expression of the MtrC-MtrD-MtrE efflux pump. Moreover, both the mtrR and penB mutations were required to decrease in vivo permeation rates below those observed in the parental strain containing either mtrR or porin mutations alone. Thus, these data demonstrate a novel mechanism of porin-mediated resistance in which mutations in PIB have no affect on antibiotic permeation alone but instead act synergistically with the MtrC-MtrD-MtrE efflux pump in the development of antibiotic resistance in gonococci.     

Genetic variation versus recombination rate in a structured population of mice

The correlation between genetic variation and recombination rate was investigated in a structured mouse population. Nucleotide sequence data from 19 autosomal DNA loci from eight inbred strains of mouse (Mus musculus) sampled from three major subspecies were analyzed. The recombination rate was estimated from the comparison of genetic and physical map distances between markers flanking a 10-cM region of each locus. The strains were categorized into four groups (subpopulations) based on geography. By partitioning the genetic diversity into within-group and among-group variation, we detected a positive correlation between the recombination rate and nucleotide diversity within groups. The level of nucleotide differentiation among groups (G(ST)) showed a negative correlation with the rate of recombination. There was no significant correlation between recombination rate and nucleotide diversity when data from different subpopulations were pooled. No correlation was detected between recombination rate and nucleotide divergence of M. musculus and M. spicilegus. These patterns deviate from the strict neutral expectation under the constant nucleotide substitution rate, and they are likely to have been formed either by a hitchhiking effect of positively selected mutants or by background selection of deleterious mutants occurring in a subdivided population. Our series of comparisons show that because a real population always has some structure, incorporation of its information is important in detecting non-neutral evolution.     

The coalescent process in models with selection, recombination and geographic subdivision

A population genetic model with a single locus at which balancing selection acts and many linked loci at which neutral mutations can occur is analysed using the coalescent approach. The model incorporates geographic subdivision with migration, as well as mutation, recombination, and genetic drift of neutral variation. It is found that geographic subdivision can affect genetic variation even with high rates of migration, providing that selection is strong enough to maintain different allele frequencies at the selected locus. Published sequence data from the alcohol dehydrogenase locus of Drosophila melanogaster are found to fit the proposed model slightly better than a similar model without subdivision.     

J-Western forms of Helicobacter pylori cagA constitute a distinct phylogenetic group with a widespread geographic distribution

Chronic infection with Helicobacter pylori strains expressing the bacterial oncoprotein CagA confers an increased risk of gastric cancer. While much is known about the ancestry and molecular evolution of Western, East Asian, and Amerindian cagA sequences, relatively little is understood about a fourth group, known as "J-Western," which has been detected mainly in strains from Okinawa, Japan. We show here that J-Western cagA sequences have a more widespread global distribution than previously recognized, occur in strains with multiple different ancestral origins (based on multilocus sequence typing [MLST] analysis), and did not arise recently. As shown by comparisons of Western and J-Western forms of CagA, there are 45 fixed or nearly fixed amino acid differences, and J-Western forms contain a unique 4-amino-acid insertion. The mean nucleotide diversity of synonymous sites (π(s)) is slightly lower in the J-Western group than in the Western and East Asian groups (0.066, 0.086, and 0.083, respectively), which suggests that the three groups have comparable, but not equivalent, effective population sizes. The reduced π(s) of the J-Western group is attributable to ancestral recombination events within the 5' region of cagA. Population genetic analyses suggest that within the cagA region encoding EPIYA motifs, the East Asian group underwent a marked reduction in effective population size compared to the Western and J-Western groups, in association with positive selection. Finally, we show that J-Western cagA sequences are found mainly in strains producing m2 forms of the secreted VacA toxin and propose that these functionally interacting proteins coevolved to optimize the gastric colonization capacity of H. pylori.     

Genetic Diversity and Mosaicism at the por Locus of Neisseria gonorrhoeae

The por genes of the predominant serovars of Neisseria gonorrhoeae circulating in a high-frequency transmitter core group located in Nairobi, Kenya, were examined for nucleotide sequence polymorphism. The level of por gene diversity did not differ significantly between core group-derived gonococcal strains and gonococcal strains originating elsewhere. However, por mosaicism appeared to be more frequent among core group-derived strains, suggesting that recombination of different por sequences may be a important strategy by which N. gonorrhoeae generates por gene diversity within core group populations. Despite extensive sequence variability, por expressed by gonococcal isolates of different geographic origin exhibited conserved patterns of nucleotide change, suggesting that diversity among por alleles may also be finite.     

Molecular evolution of key genes for type II secretion in Legionella pneumophila

Given the role of type II protein secretion system (T2S) in the ecology and pathogenesis of Legionella pneumophila, it is possible that this system is a target for adaptive evolution. The population genetic structure of L.pneumophila was inferred from the partial sequences of rpoB and from the complete sequence of three T2S structural components (lspD, lspE and pilD) and from two T2S effectors critical for intracellular infection of protozoa (proA and srnA) of 37 strains isolated from natural and man-made environments and disease-related from worldwide sources. A phylogenetic analysis was obtained for the concatenated alignment and for each individual locus. Seven main groups were identified containing the same L.pneumophila strains, suggesting an ancient divergence for each cluster and indicating that the operating selective pressures have equally affected the evolution of the five genes. Although linkage disequilibrium analysis indicate a clonal nature for population structure in this sample, our results indicate that recombination is a common phenomenon among T2S-related genes on this species, as 24 of the 37 L.pneumophila isolates contained at least one locus in which recombination was identified. Furthermore, neutral selection acting on the analysed T2S-related genes emerged as a clear result, namely on T2S effectors, ProA and SrnA, indicating that they are probably implicated in conserved virulence mechanisms through legionellae hosts.     

Recombination and selection at Brassica self-incompatibility loci

In Brassica species, self-incompatibility is controlled genetically by haplotypes involving two known genes, SLG and SRK, and possibly an as yet unknown gene controlling pollen incompatibility types. Alleles at the incompatibility loci are maintained by frequency-dependent selection, and diversity at SLG and SRK appears to be very ancient, with high diversity at silent and replacement sites, particularly in certain "hypervariable" portions of the genes. It is important to test whether recombination occurs in these genes before inferences about function of different parts of the genes can be made from patterns of diversity within their sequences. In addition, it has been suggested that, to maintain the relationship between alleles within a given S-haplotype, recombination is suppressed in the S-locus region. The high diversity makes many population genetic measures of recombination inapplicable. We have analyzed linkage disequilibrium within the SLG gene of two Brassica species, using published coding sequences. The results suggest that intragenic recombination has occurred in the evolutionary history of these alleles. This is supported by patterns of synonymous nucleotide diversity within both the SLG and SRK genes, and between domains of the SRK gene. Finally, clusters of linkage disequilibrium within the SLG gene suggest that hypervariable regions are under balancing selection, and are not merely regions of relaxed selective constraint.     

Evolutionary dynamics of Ralstonia solanacearum

We investigated the genetic diversity, extent of recombination, natural selection, and population divergence of Ralstonia solanacearum samples obtained from sources worldwide. This plant pathogen causes bacterial wilt in many crops and constitutes a serious threat to agricultural production due to its very wide host range and aggressiveness. Five housekeeping genes, dispersed around the chromosome, and three virulence-related genes, located on the megaplasmid, were sequenced from 58 strains belonging to the four major phylogenetic clusters (phylotypes). Whereas genetic variation is high and consistent for all housekeeping loci studied, virulence-related gene sequences are more diverse. Phylogenetic and statistical analyses suggest that this organism is a highly diverse bacterial species containing four major, deeply separated evolutionary lineages (phylotypes I to IV) and a weaker subdivision of phylotype II into two subgroups. Analysis of molecular variations showed that the geographic isolation and spatial distance have been the significant determinants of genetic variation between phylotypes. R. solanacearum displays high clonality for housekeeping genes in all phylotypes (except phylotype III) and significant levels of recombination for the virulence-related egl and hrpB genes, which are limited mainly to phylotype strains III and IV. Finally, genes essential for species survival are under purifying selection, and those directly involved in pathogenesis might be under diversifying selection.     

The potential protective effect of immunization with outer-membrane protein I from Neisseria gonorrhoeae

Immunization of rabbits with outer membranes (OM) of Neisseria gonorrhoeae produced antibodies directed against outer-membrane proteins PI and PIII. The antibodies directed against PIII reacted equally well on Western blots with all strains tested, but antibodies directed against PI reacted only with the homologous strain. When purified PIB was used for immunization the immune response was quite different: the sera obtained reacted with both homologous and heterologous PIB types and also reacted with strains expressing PIA. Western blotting of peptides produced by sequential cleavage of PIB revealed that the antigenic determinants recognized by anti-OM sera were predominantly located in the central surface-exposed region of PIB, as is the epitope recognized by the protective anti-PIB monoclonal antibody SM24. In contrast antibodies produced by immunization with purified PI reacted with antigenic determinants in the N-terminal portion of PIB. Nevertheless these determinants are accessible to immune attack on the native protein since the anti-PI sera were opsonic and were strongly bactericidal for both PIA- and PIB-expressing strains.     

Linking phylogenetics with population genetics to reconstruct the geographic origin of a species

Reconstructing ancestral geographic origins is critical for understanding the long-term evolution of a species. Bayesian methods have been proposed to test biogeographic hypotheses while accommodating uncertainty in phylogenetic reconstruction. However, the problem that certain taxa may have a disproportionate influence on conclusions has not been addressed. Here, we infer the geographic origin of Drosophila simulans using 2,014 bp of the period locus from 63 lines collected from 18 countries. We also analyze two previously published datasets, alcohol dehydrogenase related and NADH:ubiquinone reductase 75 kDa subunit precursor. Phylogenetic inferences of all three loci support Madagascar as the geographic origin of D. simulans. Our phylogenetic conclusions are robust to taxon resampling and to the potentially confounding effects of recombination. To test our phylogenetically derived hypothesis we develop a randomization test of the population genetics prediction that sequences from the geographic origin should contain more genetic polymorphism than those from derived populations. We find that the Madagascar population has elevated genetic polymorphism relative to non-Madagascar sequences. These data are corroborated by mitochondrial DNA sequence data.     

Evolutionary Dynamics of Ralstonia solanacearum

We investigated the genetic diversity, extent of recombination, natural selection, and population divergence of Ralstonia solanacearum samples obtained from sources worldwide. This plant pathogen causes bacterial wilt in many crops and constitutes a serious threat to agricultural production due to its very wide host range and aggressiveness. Five housekeeping genes, dispersed around the chromosome, and three virulence-related genes, located on the megaplasmid, were sequenced from 58 strains belonging to the four major phylogenetic clusters (phylotypes). Whereas genetic variation is high and consistent for all housekeeping loci studied, virulence-related gene sequences are more diverse. Phylogenetic and statistical analyses suggest that this organism is a highly diverse bacterial species containing four major, deeply separated evolutionary lineages (phylotypes I to IV) and a weaker subdivision of phylotype II into two subgroups. Analysis of molecular variations showed that the geographic isolation and spatial distance have been the significant determinants of genetic variation between phylotypes. R. solanacearum displays high clonality for housekeeping genes in all phylotypes (except phylotype III) and significant levels of recombination for the virulence-related egl and hrpB genes, which are limited mainly to phylotype strains III and IV. Finally, genes essential for species survival are under purifying selection, and those directly involved in pathogenesis might be under diversifying selection.     

Influence of conserved and hypervariable genetic markers on genotyping circulating strains of Neisseria gonorrhoeae

Presently there is no vaccine against Neisseria gonorrhoeae and therefore accurate information on gonococcal transmission plays a crucial role for interventions designed to limit the spread of infections caused by this microorganism. We evaluated the impact of two different categories of genetic markers, (i) concatenated sequences of 10 housekeeping genes and (ii) hypervariable porB DNA sequences, on the genetic relatedness and subsequently on genotyping analysis of this human pathogen. Eighty gonococcal isolates from Canada, China, the US, Argentina, Venezuela and Chile, collected over different times, were analyzed. Our results show that the choice of genetic marker had a profound effect on the interpretation of genotyping results associated with N. gonorrhoeae. The concatenated sequences of the housekeeping genes preserved the genetic relatedness of closely related isolates, enabling detection of the predominant strains circulating within a community (Saskatchewan, Canada) over an extended period of time. In contrast, a genetic marker based on antigen gene, porB, may lead to a failure to detect these predominant circulating strains. Based on the analysis of the DNA sequences of the 10 housekeeping genes, we identified two major clonal complexes, CC33 and CC22, which comprised STs from China, and Argentina as well as two STs from Canada. Several minor clonal complexes were observed among isolates from Saskatchewan. eBURST analysis suggested that the majority of the tested gonococcal isolates from Saskatchewan, Canada were endemic, with only a couple of genotypes introduced.     

High intraspecific recombination rate in a native population of Candidatus pelagibacter ubique (SAR11)

Recombination is an important process in microbial evolution. Rates of recombination with extracellular DNA matter because models of microbial population structure are profoundly influenced by the degree to which recombination is occurring within the population. Low rates of recombination may be sufficient to ensure the lateral propagation of genes that have a high selective advantage without disrupting the clonal pattern of inheritance for other genes. High rates of recombination potentially can obscure clonal patterns, leading to linkage equilibrium, and give microbial populations a population genetic structure more akin to sexually interbreeding eukaryotic populations. We examined eight loci from nine strains of candidatus Pelagibacter ubique (SAR11), isolated from a single 2L niskin sample of natural seawater, for evidence of genetic recombination between strains. The Shimodaira-Hasegawa test revealed significant phylogenetic incongruence in seven of the genes, indicating that frequent recombination obscures phylogenetic signals from the linear inheritance of genes in this population. Statistical evidence for intragenic recombination was found for six loci. An informative sites matrix showed extensive evidence for a widespread breakdown of linkage disequilibrium. Although the mechanisms of genetic transfer in native SAR11 populations are unknown, we measured recombination rates, rho, that are much higher than point mutation rates, theta, as a source of genetic diversity in this clade. The eukaryotic model of species sharing a common pool of alleles is more apt for this SAR11 population than a strictly clonal model of inheritance in which allelic diversity is controlled by periodic selection.     

Genetic structure of marine Borrelia garinii and population admixture with the terrestrial cycle of Lyme borreliosis

Despite the importance of population structure for the epidemiology of pathogenic bacteria, the spatial and ecological heterogeneity of these populations is often poorly characterized. Here, we investigated the genetic diversity and population structure of the Lyme borreliosis (LB) spirochaete Borrelia garinii in its marine cycle involving colonial seabirds and different host races of the seabird tick Ixodes uriae. Multilocus sequence analyses (MLSA) on eight chromosomal and two plasmid loci (ospA and ospC) indicate that B. garinii circulating in the marine system is highly diverse. Microevolution in marine B. garinii seems to be mainly clonal, but recombination and selection do occur. Sequence types were not evenly distributed among geographic regions, with substantial population subdivision between Atlantic and Pacific Ocean basins. However, no geographic structuring was evident within regions. Results of selection analyses and phylogenetic discordance between chromosomal and plasmid loci indicate adaptive evolution is likely occurring in this system, but no pattern of host or vector-associated divergence was found. Recombination analyses showed evidence for population admixture between terrestrial and marine strains, suggesting that LB spirochaetes are exchanged between these enzootic cycles. Importantly, our results highlight the need to explicitly consider the marine system for a complete understanding of the evolutionary ecology and global epidemiology of Lyme borreliosis.     

Multilocus patterns of polymorphism and selection across the X chromosome of Caenorhabditis remanei

Natural selection and neutral processes such as demography, mutation, and gene conversion all contribute to patterns of polymorphism within genomes. Identifying the relative importance of these varied components in evolution provides the principal challenge for population genetics. To address this issue in the nematode Caenorhabditis remanei, I sampled nucleotide polymorphism at 40 loci across the X chromosome. The site-frequency spectrum for these loci provides no evidence for population size change, and one locus presents a candidate for linkage to a target of balancing selection. Selection for codon usage bias leads to the non-neutrality of synonymous sites, and despite its weak magnitude of effect (N(e)s approximately 0.1), is responsible for profound patterns of diversity and divergence in the C. remanei genome. Although gene conversion is evident for many loci, biased gene conversion is not identified as a significant evolutionary process in this sample. No consistent association is observed between synonymous-site diversity and linkage-disequilibrium-based estimators of the population recombination parameter, despite theoretical predictions about background selection or widespread genetic hitchhiking, but genetic map-based estimates of recombination are needed to rigorously test for a diversity-recombination relationship. Coalescent simulations also illustrate how a spurious correlation between diversity and linkage-disequilibrium-based estimators of recombination can occur, due in part to the presence of unbiased gene conversion. These results illustrate the influence that subtle natural selection can exert on polymorphism and divergence, in the form of codon usage bias, and demonstrate the potential of C. remanei for detecting natural selection from genomic scans of polymorphism.     

ospC diversity in Borrelia burgdorferi: different hosts are different niches

The outer surface protein C (ospC) locus of the Lyme disease bacterium, Borrelia burgdorferi, is at least an order of magnitude more variable than other genes in the species. This variation is classified into 22 ospC major groups, 15 of which are found in the northeastern United States. The frequency distributions of ospC within populations suggest that this locus is under balancing selection. In multiple-niche polymorphism, a type of balancing selection, diversity within a population can be maintained when the environment is heterogeneous and no one genotype has the highest fitness in all environments. Genetically different individuals within vertebrate species and different vertebrate species constitute diverse environments for B. burgdorferi. We examined four important host species of B. burgdorferi and found that the strains that infected each species had different sets of ospC major groups. We found no variation among conspecific hosts in the ospC major groups of their infecting strains. These results suggest multiple niches create balancing selection at the ospC locus.     

Hitchhiking effect of a beneficial mutation spreading in a subdivided population

A central problem in population genetics is to detect and analyze positive natural selection by which beneficial mutations are driven to fixation. The hitchhiking effect of a rapidly spreading beneficial mutation, which results in local removal of standing genetic variation, allows such an analysis using DNA sequence polymorphism. However, the current mathematical theory that predicts the pattern of genetic hitchhiking relies on the assumption that a beneficial mutation increases to a high frequency in a single random-mating population, which is certainly violated in reality. Individuals in natural populations are distributed over a geographic space. The spread of a beneficial allele can be delayed by limited migration of individuals over the space and its hitchhiking effect can also be affected. To study this effect of geographic structure on genetic hitchhiking, we analyze a simple model of directional selection in a subdivided population. In contrast to previous studies on hitchhiking in subdivided populations, we mainly investigate the range of sufficiently high migration rates that would homogenize genetic variation at neutral loci. We provide a heuristic mathematical analysis that describes how the genealogical structure at a neutral locus linked to the locus under selection is expected to change in a population divided into two demes. Our results indicate that the overall strength of genetic hitchhiking--the degree to which expected heterozygosity decreases--is diminished by population subdivision, mainly because opportunity for the breakdown of hitchhiking by recombination increases as the spread of the beneficial mutation across demes is delayed when migration rate is much smaller than the strength of selection. Furthermore, the amount of genetic variation after a selective sweep is expected to be unequal over demes: a greater reduction in expected heterozygosity occurs in the subpopulation from which the beneficial mutation originates than in its neighboring subpopulations. This raises a possibility of detecting a "hidden" geographic structure of population by carefully analyzing the pattern of a selective sweep.