Dynamics of linkage disequilibrium in bacterial genomes undergoing transformation and/or conjugation

Bacteria may undergo recombinational exchange either by conjugation followed by crossing over, or by transformation of small segments of DNA into the cell followed by incorporation into the chromosome by gene conversion. These two forms of recombination may have very different consequences on the patterns of linkage disequilibrium seen within bacterial genomes. In this paper deterministic recursions are obtained for three linked loci in populations having these two forms of recombination. Both neutral genetic variation and the case of one selected gene are considered. It is shown that the two forms of exchange have identical consequences on two-locus linkage disequilibria, but that three-locus disequilibria can have different behaviors. Hitchhiking also has different consequences on the pattern of disequilibrium seen between linked neutral genes in the region of the selected locus. Inference of the relative importance of these two modes of recombination from static samples of DNA sequences will hinge on the relationship between linkage map distance and disequilibria.     

Corn and humans: recombination and linkage disequilibrium in two genomes of similar size

Two species with genomes of almost identical size, maize and human, have different evolutionary histories, and as a result their genomes differ greatly in their content of retroelements, average size of the genes and amount of genetic diversity. However, there are also significant similarities: they both have undergone bottlenecks during their recent history and seem to have non-uniform distribution of recombination events. The human genome has been shown to contain large linkage blocks characterized by a limited number of haplotypes. A similar linkage block structure is likely to exist in maize. Although highly diverse maize populations show rapid decline of linkage disequilibrium, as in humans, it is possible to define populations with strong linkage disequilibrium, suitable for whole-genome scan association mapping. The genetic diversity and lack of sequence homology found in maize influences recombinational properties and local linkage disequilibrium levels but also challenges our understanding of the relationship between the genome sequence and species definition.     

The relative contributions of recombination and mutation to the divergence of clones of Neisseria meningitidis.

Multilocus sequence typing (MLST) is a recently developed nucleotide sequence-based method for the definitive assignment of isolates within bacterial populations to specific clones. MLST uses the same principles as multilocus enzyme electrophoresis and provides data that can be used to investigate aspects of the population genetics and evolution of bacterial species. We used an MLST data set consisting of the sequences of approximately 450-bp fragments from seven housekeeping loci from a large strain collection of Neisseria meningitidis to estimate the relative impact of recombination compared with point mutation in the diversification of N. meningitidis clonal complexes. 126 meningococcal isolates were assigned to 10 clonal complexes, 9 of which contained minor clonal variants. The allelic variation within each complex was classified as a recombinational exchange or a putative point mutation through a comparison of the sequences of each variant allele with that of the allele typically found in the clonal complex. The nine clonal complexes contained a total of 23 allelic variants, and analysis of the sequences of these variant alleles revealed that a single nucleotide site in a meningococcal housekeeping gene is at least 80-fold more likely to change as a result of recombination than as a result of mutation. This value is estimated to be 10-50-fold for Escherichia coli and approximately 50-fold for Streptococcus pneumoniae.     

Expected Linkage Disequilibrium for a Neutral Locus Linked to a Chromosomal Arrangement

The expected value of the squared linkage disequilibrium is derived for a neutral locus associated with a chromosomal arrangement that is maintained in the population by strong balancing selection. For a given value of recombination, the expected squared linkage disequilibrium is shown to decrease as the intensity of selection maintaining the arrangement increases. The transient behavior of the expected square linkage disequilibrium is also derived. This theory applies to loci that are closely linked to inversions in Drosophila species and to loci closely linked to the differential segments of the translocation complexes in ring-forming species of Oenothera. In both cases the strong linkage disequilibria that have been observed in natural populations can be explained by random drift.     

Ecological separation and genetic isolation of Neisseria gonorrhoeae and Neisseria meningitidis

BACKGROUND: Classifying bacteria into species is problematic. Most microbiologists consider species to be groups of isolates that share some arbitrary degree of relatedness of biochemical or molecular (such as DNA sequence) features and that, ideally, are clearly delineated from all other groups of isolates. The main problem in applying to bacteria a biological concept of species based on the ability or inability of their genes to recombine, is that recombination appears to be rare in bacteria in nature, as indicated by the strong linkage disequilibrium between alleles found in most bacterial populations. However, there are some naturally transformable bacteria in which assortative recombination appears to be so frequent that alleles are in, or close to, linkage equilibrium. For these recombining populations a biological concept of species might be applicable. RESULTS: Populations of Neisseria gonorrhoeae and Neisseria meningitidis from Spain were analysed by multilocus enzyme electrophoresis. The data indicate that assortative recombination occurs frequently within populations, but not between populations. Similarly, the sequences of two house-keeping genes show no evidence of intragenic recombination between N. gonorrhoeae and N. meningitidis. CONCLUSIONS: N. gonorrhoeae and N. meningitidis represent extremely closely related 'sexual' populations that appear to be genetically isolated in nature, and thus conform to the biological concept of species. The extreme uniformity of N. gonorrhoeae house-keeping genes suggests that this species may have arisen recently as a clone of N. meningitidis that could colonize the genital tract. Ecological isolation - of populations that can colonize the genital tract from those that can colonize the nasopharynx - may have been an important component in speciation, leading to a lower frequency of recombination between species than within species.     

Variation in recombination rate across the genome: evidence and implications

Recent data from humans and other species provide convincing evidence of variation in recombination rate in different genomic regions. Comparison of physical and genetic maps reveals variation on a scale of megabases, with substantial differences between sexes. Recombination is often suppressed near centromeres and elevated near telomeres, but neither of these observations is true for all chromosomes. In humans, patterns of linkage disequilibrium and experimental measures of recombination from sperm-typing reveal dramatic hotspots of recombination on a scale of kilobases. Genome-wide variation in the amount of crossing-over may be due to variation in the density of hotspots, the intensity of hotspots, or both. Theoretical models of selection and linkage predict that genetic variation will be reduced in regions of low recombination, and this prediction is supported by data from several species. Heterogeneity in rates of crossing-over provides both an opportunity and a challenge for identifying disease genes: as associations occur in blocks, genomic regions containing disease loci may be identified with relatively few markers, yet identifying the causal mutations is unlikely to be achieved through associations alone.     

Evidence for rare recombination at the gametophytic self-incompatibility locus

The gametophytic self-incompatibility locus has been thought to be a nonrecombining genomic region. Inferences have been made, however, about the functional importance of different parts of the S-locus, based on differences in the levels of variability along the gene, and this is valid only if recombination occurs. It is thus important to test whether recombination occurs within and near the S-locus. Several recent attempts to test this have reached conflicting conclusions. In this study, we examine a large data set on sequence variation at the S-locus in several species with gametophytic self-incompatibility systems, in the Solanaceae, Rosaceae and Scrophulariaceae. We use the longest sequences available to test for recombination based on linkage disequilibrium between polymorphic sites in the S-locus. The relationship between linkage disequilibrium and physical distance between the sites suggests rare intragenic exchange in the evolutionary history of four species of Solanaceae and two species of Rosaceae.     

Recombination, Diversity and Allele Sharing of Infectivity Proteins Between Bartonella Species from Rodents

The alpha-Proteobacterium Bartonella is a common parasite of voles and mice, giving rise to short-lived (4 weeks to 2 months) infections. Here, we report high sequence diversity in genes of the VirB/VirD type IV secretion system (T4SS), amongst Bartonella from natural rodent populations in NE Poland. The VirB5 protein is predicted to consist of three conserved alpha helices separated by loops of variable length which include numerous indels. The C-terminal domain includes repeat stretches of KEK residues, reflecting underlying homopolymeric stretches of adenine residues. A total of 16 variants of VirB5, associated with host identity, but not bacterial taxon, were identified from 22 Bartonella isolates. One was clearly a recombinant from two others, another included an insertion of two KEK repeats. The virB5 gene appears to evolve via both mutation and recombination, as well as slippage mediated insertion/deletion events. The recombinational units are thought to be relatively short, as there was no evidence of linkage disequilibrium between virB5 and the bepA locus only 5.5 kb distant. The diversity of virB5 is assumed to be related to immunological role of this protein in Bartonella infections; diversity of virB5 may assist persistence of Bartonella in the rodent population, despite the relatively short (3-4 weeks) duration of individual infections. It is clear from the distribution of virB5 and bepA alleles that recombination within and between clades is widespread, and frequently crosses the boundaries of conventionally recognised Bartonella species.     

The beta -globin recombinational hotspot reduces the effects of strong selection around HbC, a recently arisen mutation providing resistance to malaria

Recombination is expected to reduce the effect of selection on the extent of linkage disequilibrium (LD), but the impact that recombinational hotspots have on sites linked to selected mutations has not been investigated. We empirically determine chromosomal linkage phase for 5.2 kb spanning the beta -globin gene and hotspot. We estimate that the HbC mutation, which is positively selected because of malaria, originated <5,000 years ago and that selection coefficients are 0.04-0.09. Despite strong selection and the recent origin of the HbC allele, recombination (crossing-over or gene conversion) is observed within 1 kb 5' of the selected site on more than one-third of the HbC chromosomes sampled. The rapid decay in LD upstream of the HbC allele demonstrates the large effect the ss-globin hotspot has in mitigating the effects of positive selection on linked variation.     

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.     

Recombinational exchanges at the capsular polysaccharide biosynthetic locus lead to frequent serotype changes among natural isolates of Streptococcus pneumoniae

Serotype 19F variants of the major Spanish multiresistant serotype 23F clone of Streptococcus pneumoniae have been proposed to have arisen by recombinational exchanges at the capsular biosynthetic locus. Members of the Spanish multiresistant serotype 23F clone and the serotype 19F variants were confirmed to be essentially identical in overall genotype, as they were indistinguishable by REP-PCR, and had identical sequences at three polymorphic housekeeping genes. Eight serotype 19F variants were studied and all had large recombinational replacements at the capsular biosynthetic locus. In all cases, one of the recombinational cross-over points appeared to be upstream of dexB which flanks one end of the capsular locus, and in six of the variants the other cross-over point was downstream of aliA , which flanks the other end of the locus. In two strains a recombinational cross-over point between the introduced serotype 19F capsular region and that of the Spanish serotype 23F clone could be clearly identified, within cpsN in one strain and within cpsM in the other. The differences in the recombinational junctions and sequence polymorphisms within the introduced capsular genes, suggested that the eight serotype 19F variants emerged on at least four separate occasions. Changes in capsular type by recombination may therefore be relatively frequent in pneumococci and this has implications for the long-term efficacy of conjugate pneumococcal vaccines that will protect against only a limited number of serotypes.     

Evolving Concepts of Bacterial Species

The same evolutionary forces that cause diversification in sexual eukaryotes are expected to cause diversification in bacteria. However, in bacteria, the wider variety of mechanisms for gene exchange (or lack thereof) increases the range of expected diversity patterns compared to those of sexual organisms. Two parallel concepts for bacterial speciation have developed, based on ecological divergence or barriers to recombination in turn. Recent evidence from DNA sequence data shows that both processes can generate independently evolving groups that are equivalent to sexual species and that represent separate arenas within which recombination (when it occurs), selection and drift occur. It remains unclear, however, how often different processes act in concert to generate simple units of diversity, or whether a more complex model of diversity is required, specifying hierarchical levels at which different cohesive processes operate. We advocate an integrative approach that evaluates the effects of multiple evolutionary forces on diversity patterns. There is also great potential for laboratory studies of bacterial evolution that test evolutionary mechanisms inferred from population genetic analyses of multi-locus and genome sequence data.     

Analysis of genetic structure in soil populations of Rhizobium leguminosarum recovered from the USA and the UK

Although many studies have shown that animal-associated bacterial species exhibit linkage disequilibrium at chromosomal loci, recent studies indicate that both animal-associated and soil-borne bacterial species can display a nonclonal genetic structure in which alleles at chromosomal loci are in linkage equilibrium. To examine the situation in soil-borne species further, we compared genetic structure in two soil populations of Rhizobium leguminosarum bv. trifolii and two populations of R. leguminosarum bv. viciae from two sites in Oregon, with genetic structure in R. leguminosarum bv. viciae populations recovered from peas grown at a site in Washington, USA, and at a site in Norfolk, UK. A total of 234 chromosomal types (ET) were identified among 682 strains analysed for allelic variation at 13 enzyme-encoding chromosomal loci by multilocus enzyme electrophoresis (MLEE). Chi-square tests for heterogeneity of allele frequencies showed that the populations were not genetically uniform. A comparison of the genetic diversity within combined and individual populations confirmed that the Washington population was the primary cause of genetic differentiation between the populations. Each individual population exhibited linkage disequilibrium, with the magnitude of the disequilibrium being greatest in the Washington population and least in the UK population of R. leguminosarum bv. viciae. Linkage disequilibrium in the UK population was created between two clusters of 9 and 23 ETs, which, individually, were in linkage equilibrium. Strong linkage disequilibrium between the two major clusters of 8 and 12 ETs in the Washington population was caused by the low genetic diversity of the ETs within each cluster relative to the inter-cluster genetic distance. Because neither the magnitude of genetic diversity nor of linkage disequilibrium increased as hierarchical combinations of the six local populations were analysed, we conclude that the populations have not been isolated from each other for sufficient time, nor have they been exposed to enough selective pressure to develop unique multilocus genetic structure.     

Relative roles of mutation and recombination in generating allelic polymorphism at an MHC class II locus in Peromyscus maniculatus

The MHC class II loci encoding cell surface antigens exhibit extremely high allelic polymorphism. There is considerable uncertainty in the literature over the relative roles of recombination and de novo mutation in generating this diversity. We studied class II sequence diversity and allelic polymorphism in two populations of Peromyscus maniculatus, which are among the most widespread and abundant mammals of North America. We find that intragenic recombination (or gene conversion) has been the predominant mode for the generation of allelic polymorphism in this species, with the amount of population recombination per base pair exceeding mutation by at least an order of magnitude during the history of the sample. Despite this, patchwork motifs of sites with high linkage disequilibrium are observed. This does not appear to be consistent with the much larger amount of recombination versus mutation in the history of the sample, unless the recombination rate is highly non-uniform over the sequence or selection maintains certain sites in linkage disequilibrium. We conclude that selection is most likely to be responsible for preserving sequence motifs in the presence of abundant recombination.     

The influence of linkage and inbreeding on patterns of nucleotide sequence diversity at duplicate alcohol dehydrogenase loci in wild barley (Hordeum vulgare ssp. spontaneum)

Patterns of nucleotide sequence diversity are analyzed for three duplicate alcohol dehydrogenase loci (adh1-adh3) within a species-wide sample of 25 accessions of wild barley (Hordeum vulgare ssp. spontaneum). The adh1 and adh2 loci are tightly linked (recombination fraction <0.01) while the adh3 locus is inherited independently. Wild barley is predominantly self-fertilizing (approximately 98%), and as a consequence, effective recombination is restricted by the extreme reduction in heterozygosity. Large reductions in effective recombination, in turn, widen the conditions for linkage to influence nucleotide sequence diversity through the action of selective sweeps or background selection. These considerations would appear to predict (1) homogeneity in patterns of nucleotide sequence diversity, especially between closely linked loci, and (2) extensive linkage disequilibrium relative to random-mating species. In contrast to these expectations, the wild barley data reveal heterogeneity in patterns of nucleotide sequence diversity and levels of linkage disequilibrium that are indistinguishable from those observed at adh1 in maize, an outbreeding grass species.     

Linkage disequilibrium under skewed offspring distribution among individuals in a population

Correlations in coalescence times between two loci are derived under selectively neutral population models in which the offspring of an individual can number on the order of the population size. The correlations depend on the rates of recombination and random drift and are shown to be functions of the parameters controlling the size and frequency of these large reproduction events. Since a prediction of linkage disequilibrium can be written in terms of correlations in coalescence times, it follows that the prediction of linkage disequilibrium is a function not only of the rate of recombination but also of the reproduction parameters. Low linkage disequilibrium is predicted if the offspring of a single individual frequently replace almost the entire population. However, high linkage disequilibrium can be predicted if the offspring of a single individual replace an intermediate fraction of the population. In some cases the model reproduces the standard Wright-Fisher predictions. Contrary to common intuition, high linkage disequilibrium can be predicted despite frequent recombination, and low linkage disequilibrium under infrequent recombination. Simulations support the analytical results but show that the variance of linkage disequilibrium is very large.     

Heterogeneity in the Frequency and Characteristics of Homologous Recombination in Pneumococcal Evolution

The bacterium Streptococcus pneumoniae (pneumococcus) is one of the most important human bacterial pathogens, and a leading cause of morbidity and mortality worldwide. The pneumococcus is also known for undergoing extensive homologous recombination via transformation with exogenous DNA. It has been shown that recombination has a major impact on the evolution of the pathogen, including acquisition of antibiotic resistance and serotype-switching. Nevertheless, the mechanism and the rates of recombination in an epidemiological context remain poorly understood. Here, we proposed several mathematical models to describe the rate and size of recombination in the evolutionary history of two very distinct pneumococcal lineages, PMEN1 and CC180. We found that, in both lineages, the process of homologous recombination was best described by a heterogeneous model of recombination with single, short, frequent replacements, which we call micro-recombinations, and rarer, multi-fragment, saltational replacements, which we call macro-recombinations. Macro-recombination was associated with major phenotypic changes, including serotype-switching events, and thus was a major driver of the diversification of the pathogen. We critically evaluate biological and epidemiological processes that could give rise to the micro-recombination and macro-recombination processes.     

Evolution of recombination due to random drift

In finite populations subject to selection, genetic drift generates negative linkage disequilibrium, on average, even if selection acts independently (i.e., multiplicatively) upon all loci. Negative disequilibrium reduces the variance in fitness and hence, by Fisher's (1930) fundamental theorem, slows the rate of increase in mean fitness. Modifiers that increase recombination eliminate the negative disequilibria that impede selection and consequently increase in frequency by "hitchhiking." Thus, stochastic fluctuations in linkage disequilibrium in finite populations favor the evolution of increased rates of recombination, even in the absence of epistatic interactions among loci and even when disequilibrium is initially absent. The method developed within this article allows us to quantify the strength of selection acting on a modifier allele that increases recombination in a finite population. The analysis indicates that stochastically generated linkage disequilibria do select for increased recombination, a result that is confirmed by Monte Carlo simulations. Selection for a modifier that increases recombination is highest when linkage among loci is tight, when beneficial alleles rise from low to high frequency, and when the population size is small.     

Estimating recombinational parameters in Streptococcus pneumoniae from multilocus sequence typing data

Multilocus sequence typing (MLST) is a highly discriminatory molecular typing method that defines isolates of bacterial pathogens using the sequences of approximately 450-bp internal fragments of seven housekeeping genes. This technique has been applied to 575 isolates of Streptococcus pneumoniae and identifies a number of discrete clonal complexes. These clonal complexes are typically represented by a single group of isolates sharing identical alleles at all seven loci, plus single-locus variants that differ from this group at only one out of the seven loci. As MLST is highly discriminatory, the members of each clonal complex can be assumed to have a recent common ancestor, and the molecular events that give rise to the single-locus variants can be used to estimate the relative contributions of recombination and mutation to clonal divergence. By comparing the sequences of the variant alleles within each clonal complex with the allele typically found within that clonal complex, we estimate that recombination has generated new alleles at a frequency approximately 10-fold higher than mutation, and that a single nucleotide site is approximately 50 times more likely to change through recombination than mutation. We also demonstrate how to estimate the average length of recombinational replacements from MLST data.     

Haplotype fine mapping by evolutionary trees

To refine the location of a disease gene within the bounds provided by linkage analysis, many scientists use the pattern of linkage disequilibrium between the disease allele and alleles at nearby markers. We describe a method that seeks to refine location by analysis of "disease" and "normal" haplotypes, thereby using multivariate information about linkage disequilibrium. Under the assumption that the disease mutation occurs in a specific gap between adjacent markers, the method first combines parsimony and likelihood to build an evolutionary tree of disease haplotypes, with each node (haplotype) separated, by a single mutational or recombinational step, from its parent. If required, latent nodes (unobserved haplotypes) are incorporated to complete the tree. Once the tree is built, its likelihood is computed from probabilities of mutation and recombination. When each gap between adjacent markers is evaluated in this fashion and these results are combined with prior information, they yield a posterior probability distribution to guide the search for the disease mutation. We show, by evolutionary simulations, that an implementation of these methods, called "FineMap," yields substantial refinement and excellent coverage for the true location of the disease mutation. Moreover, by analysis of hereditary hemochromatosis haplotypes, we show that FineMap can be robust to genetic heterogeneity.