Methods for estimating gene frequencies and detecting selection in bacterial populations

Maximum antigen diversification in a lyme bacterial population and evolutionary strategies to overcome pathogen diversity

Natural populations of pathogens and their hosts are engaged in an arms race in which the pathogens diversify to escape host immunity while the hosts evolve novel immunity. This co-evolutionary process poses a fundamental challenge to the development of broadly effective vaccines and diagnostics against a diversifying pathogen. Based on surveys of natural allele frequencies and experimental immunization of mice, we show high antigenic specificities of natural variants of the outer surface protein C (OspC), a dominant antigen of a Lyme Disease-causing bacterium (Borrelia burgdorferi). To overcome the challenge of OspC antigenic diversity to clinical development of preventive measures, we implemented a number of evolution-informed strategies to broaden OspC antigenic reactivity. In particular, the centroid algorithm-a genetic algorithm to generate sequences that minimize amino-acid differences with natural variants-generated synthetic OspC analogs with the greatest promise as diagnostic and vaccine candidates against diverse Lyme pathogen strains co-existing in the Northeast United States. Mechanistically, we propose a model of maximum antigen diversification (MAD) mediated by amino-acid variations distributed across the hypervariable regions on the OspC molecule. Under the MAD hypothesis, evolutionary centroids display broad cross-reactivity by occupying the central void in the antigenic space excavated by diversifying natural variants. In contrast to vaccine designs based on concatenated epitopes, the evolutionary algorithms generate analogs of natural antigens and are automated. The novel centroid algorithm and the evolutionary antigen designs based on consensus and ancestral sequences have broad implications for combating diversifying pathogens driven by pathogen-host co-evolution.

Differing lifestyles of Staphylococcus epidermidis as revealed through Bayesian clustering of multilocus sequence types

Staphylococcus epidermidis is part of the normal bacterial flora of human skin and a leading cause of infections associated with indwelling medical devices. Previous phylogenetic analyses of subgenomic data have been unable to distinguish between S. epidermidis strains with nosocomial or commensal lifestyles, despite the identification of specific phenotypes and accessory genes that may contribute to such lifestyles. To attempt to better define the population structure of this species, the international S. epidermidis multilocus sequence typing database was analyzed with the Bayesian clustering programs STRUCTURE and BAPS. A total of six genetic clusters (GCs) were identified. A local population of S. epidermidis from clinical specimens was classified according to these six GCs, and further characterized for antibiotic susceptibilities, biofilm, and various genetic markers. GC5 was abundant and significantly enriched for isolates that were resistant to four classes of antibiotics, high biofilm production, and positive for the virulence markers icaA, IS256, and sesD/bhp, indicating its potential clinical relevance. In contrast, GC2 was rare and contained the only isolates positive for the putative commensal marker, fdh. GC1 and GC6 were abundant but not significantly associated with any of the examined characteristics, except for sesF/aap and GC6. GC3 was rare and identified as a potential genetic sink that received, but did not donate, core genetic material from other GCs. In conclusion, population genetics analyses were essential for identifying clusters of strains that may differ in their adaptation to nosocomial or commensal lifestyles. These results provide a new, population genetics framework for studying S. epidermidis.

Pervasive recombination and sympatric genome diversification driven by frequency-dependent selection in Borrelia burgdorferi, the Lyme disease bacterium

How genomic diversity within bacterial populations originates and is maintained in the presence of frequent recombination is a central problem in understanding bacterial evolution. Natural populations of Borrelia burgdorferi, the bacterial agent of Lyme disease, consist of diverse genomic groups co-infecting single individual vertebrate hosts and tick vectors. To understand mechanisms of sympatric genome differentiation in B. burgdorferi, we sequenced and compared 23 genomes representing major genomic groups in North America and Europe. Linkage analysis of >13,500 single-nucleotide polymorphisms revealed pervasive horizontal DNA exchanges. Although three times more frequent than point mutation, recombination is localized and weakly affects genome-wide linkage disequilibrium. We show by computer simulations that, while enhancing population fitness, recombination constrains neutral and adaptive divergence among sympatric genomes through periodic selective sweeps. In contrast, simulations of frequency-dependent selection with recombination produced the observed pattern of a large number of sympatric genomic groups associated with major sequence variations at the selected locus. We conclude that negative frequency-dependent selection targeting a small number of surface-antigen loci (ospC in particular) sufficiently explains the maintenance of sympatric genome diversity in B. burgdorferi without adaptive divergence. We suggest that pervasive recombination makes it less likely for local B. burgdorferi genomic groups to achieve host specialization. B. burgdorferi genomic groups in the northeastern United States are thus best viewed as constituting a single bacterial species, whose generalist nature is a key to its rapid spread and human virulence.

Geographical distribution of selected and putatively neutral SNPs in Southeast Asian malaria parasites

Loci targeted by directional selection are expected to show elevated geographical population structure relative to neutral loci, and a flurry of recent papers have used this rationale to search for genome regions involved in adaptation. Studies of functional mutations that are known to be under selection are particularly useful for assessing the utility of this approach. Antimalarial drug treatment regimes vary considerably between countries in Southeast Asia selecting for local adaptation at parasite loci underlying resistance. We compared the population structure revealed by 10 nonsynonymous mutations (nonsynonymous single-nucleotide polymorphisms [nsSNPs]) in four loci that are known to be involved in antimalarial drug resistance, with patterns revealed by 10 synonymous mutations (synonymous single-nucleotide polymorphisms [sSNPs]) in housekeeping genes or genes of unknown function in 755 Plasmodium falciparum infections collected from 13 populations in six Southeast Asian countries. Allele frequencies at known nsSNPs underlying resistance varied markedly between locations (F(ST) = 0.18-0.66), with the highest frequencies on the Thailand-Burma border and the lowest frequencies in neighboring Lao PDR. In contrast, we found weak but significant geographic structure (F(ST) = 0-0.14) for 8 of 10 sSNPs. Importantly, all 10 nsSNPs showed significantly higher F(ST) (P < 8 x 10(-5)) than simulated neutral expectations based on observed F(ST) values in the putatively neutral sSNPs. This result was unaffected by the methods used to estimate allele frequencies or the number of populations used in the simulations. Given that dense single-nucleotide polymorphism (SNP) maps and rapid SNP assay methods are now available for P. falciparum, comparing genetic differentiation across the genome may provide a valuable aid to identifying parasite loci underlying local adaptation to drug treatment regimes or other selective forces. However, the high proportion of polymorphic sites that appear to be under balancing selection (or linked to selected sites) in the P. falciparum genome violates the central assumption that selected sites are rare, which complicates identification of outlier loci, and suggests that caution is needed when using this approach.

Lyme borreliosis agents and the genetics and sex of their vector, Ixodes ricinus

The tick Ixodes ricinus is responsible for the transmission and maintenance of a wide variety of pathogenic organisms in the Northern Hemisphere, among which Lyme disease represents a major threat to humans. Despite numerous studies, the epidemiology of the different bacterial species responsible for this disease remains unclear. Recent evidence for a sex-biased genetic structure of its European vector leads us to analyse the consequences of this pattern on Borrelia transmission. Here we show that male and female ticks are not equivalently infected by Borrelia burgdorferi, that Borrelia afzelii affects tick migration capabilities, especially for the most vagile sex (i.e., male) and that Lyme borreliosis agents are consequently vectorised in a much more complex way than usually thought. Such results change the epidemiological perception of Lyme borreliosis and suggest new co-evolutionary pathways between the ticks and the borrelia.

Allele Frequency-Based Analyses Robustly Map Sequence Sites Under Balancing Selection in a Malaria Vaccine Candidate Antigen

The Plasmodium falciparum apical membrane antigen 1 (AMA1) is a leading candidate for a malaria vaccine. Here, within-population analyses of alleles from 50 Thai P. falciparum isolates yield significant evidence for balancing selection on polymorphisms within the disulfide-bonded domains I and III of the surface accessible ectodomain of AMA1, a result very similar to that seen previously in a Nigerian population. Studying the frequency of nucleotide polymorphisms in both populations shows that the between-population component of variance (FST) is significantly lower in domains I and III compared to the intervening domain II and compared to 11 unlinked microsatellite loci. A nucleotide site-by-site analysis shows that sites with exceptionally high or low FST values cluster significantly into serial runs, with four runs of low values in domain I and one in domain III. These runs may map the sequences that are consistently under the strongest balancing selection from naturally acquired immune responses.

Measuring immune selection

Immune responses that kill pathogens or reduce their reproductive rate are generally important in protecting hosts from infection and disease. Pathogens that escape the full impact of such responses will survive, and any heritable genetic basis of this evasion will be selected. Due to the memory component of vertebrate immune responses, pathogens with rare alleles of a target antigen can have an advantage over those with common alleles, leading to the maintenance of a polymorphism. At the genetic level, there ought to be detectable signatures of balancing selection in the genes encoding these antigens. Here, methods for identifying these selective signatures are reviewed. Their practical utility for identifying which antigens are targets of protective immune responses is discussed.

Geographic uniformity of the Lyme disease spirochete (Borrelia burgdorferi) and its shared history with tick vector (Ixodes scapularis) in the Northeastern United States

Over 80% of reported cases of Lyme disease in the United States occur in coastal regions of northeastern and mid-Atlantic states. The genetic structure of the Lyme disease spirochete (Borrelia burgdorferi) and its main tick vector (Ixodes scapularis) was studied concurrently and comparatively by sampling natural populations of I. scapularis ticks along the East Coast from 1996 to 1998. Borrelia is genetically highly diverse at the outer surface protein ospC. Since Borrelia is highly clonal, the ospC alleles can be used to define clones. A newly designed reverse line blotting (RLB) assay shows that up to 10 Borrelia clones can infect a single tick. The clone frequencies in Borrelia populations are the same across the Northeast. On the other hand, I. scapularis populations show strong regional divergence (among northeastern, mid-Atlantic, and southern states) as well as local differentiation. The high genetic diversity within Borrelia populations and the disparity in the genetic structure between Borrelia and its tick vector are likely consequences of strong balancing selection on local Borrelia clones. Demographically, both Borrelia and I. scapularis populations in the Northeast show the characteristics of a species that has recently expanded from a population bottleneck. Major geological and ecological events, such as the last glacial maximum (18,000 years ago) and the modern-day expansion of tick habitats, are likely causes of the observed "founder effects" for the two organisms in the Northeast. We therefore conclude that the genetic structure of B. burgdorferi has been intimately shaped by the natural history of its main vector, the northern lineage of I. scapularis ticks.