Estimating the Relative Roles of Recombination and Point Mutation in the Generation of Single Locus Variants in Campylobacter jejuni and Campylobacter coli

Single locus variants (SLVs) are bacterial sequence types that differ at only one of the seven canonical multilocus sequence typing (MLST) loci. Estimating the relative roles of recombination and point mutation in the generation of new alleles that lead to SLVs is helpful in understanding how organisms evolve. The relative rates of recombination and mutation for Campylobacter jejuni and Campylobacter coli were estimated at seven different housekeeping loci from publically available MLST data. The probability of recombination generating a new allele that leads to an SLV is estimated to be roughly seven times more than that of mutation for C. jejuni, but for C. coli recombination and mutation were estimated to have a similar contribution to the generation of SLVs. The majority of nucleotide differences (98?% for C. jejuni and 85?% for C. coli) between strains that make up an SLV are attributable to recombination. These estimates are much larger than estimates of the relative rate of recombination to mutation calculated from more distantly related isolates using MLST data. One explanation for this is that purifying selection plays an important role in the evolution of Campylobacter. A simulation study was performed to test the performance of our method under a range of biologically realistic parameters. We found that our method performed well when the recombination tract length was longer than 3?kb. For situations in which recombination may occur with shorter tract lengths, our estimates are likely to be an underestimate of the ratio of recombination to mutation, and of the importance of recombination for creating diversity in closely related isolates. A parametric bootstrap method was applied to calculate the uncertainty of these estimates.     

The impact of homologous recombination on the generation of diversity in bacteria

The imprint of demographic and selective processes on bacterial population structure needs to be evaluated as deviation from the expectations of an appropriate null neutral model. We explore the impact of varying the population mutation and recombination rates theta and rho on ideal populations, using a recently developed model of neutral drift at multiple loci. This model may be fitted to experimental data to provide estimates of these parameters, and we do so for seven bacterial species (Neisseria meningitidis, Streptococcus pneumoniae, Streptococcus pyogenes, Staphylococcus aureus, Helicobacter pylori, Burkholderia pseudomallei and Bacillus cereus), illustrating that bacterial species vary extensively in these fundamental parameters. Historically, the influence of recombination has often been estimated through its influence on the Index of Association I(A). We show that this may be relatively insensitive to changes in either mutation or recombination rates. It is known that biased sampling can lead to artificially high estimates of I(A). We therefore provide a method of precisely separating the effects of such bias and true linkage between alleles. We also demonstrate that by fitting the neutral model to experimental data, more informative and precise estimates of the relative roles of recombination and mutation may be obtained.     

A genome-wide departure from the standard neutral model in natural populations of Drosophila

We analyze nucleotide polymorphism data for a large number of loci in areas of normal to high recombination in Drosophila melanogaster and D. simulans (24 and 16 loci, respectively). We find a genome-wide, systematic departure from the neutral expectation for a panmictic population at equilibrium in natural populations of both species. The distribution of sequence-based estimates of 2Nc across loci is inconsistent with the assumptions of the standard neutral theory, given the observed levels of nucleotide diversity and accepted values for recombination and mutation rates. Under these assumptions, most estimates of 2Nc are severalfold too low; in other words, both species exhibit greater intralocus linkage disequilibrium than expected. Variation in recombination or mutation rates is not sufficient to account for the excess of linkage disequilibrium. While an equilibrium island model does not seem to account for the data, more complicated forms of population structure may. A proper test of alternative demographic models will require loci to be sampled in a more consistent fashion.     

Use of fluorescent protein to analyse recombination at three loci in Neurospora crassa

We have inserted a histone H1-GFP fusion gene adjacent to three loci on different chromosomes of Neurospora crassa and made mating pairs in which a wild type version of GFP is crossed to one with a mutation in the 5' end of GFP. The loci are his-3, am and his-5, chosen because recombination mechanisms appear to differ between his-3 and am, and because crossing over adjacent to his-5, like his-3, is regulated by rec-2. At his-3, the frequencies of crossing over between GFP and the centromere and of conversion of 5'GFP to GFP(+) are comparable to those obtained by classical recombination assays, as is the effect of rec-2 on these frequencies, suggesting that our system does not alter the process of recombination. At each locus we have obtained sufficient data, on both gene conversion and crossing over, to be able to assess the effect of deletion of any gene involved in recombination. In addition, crosses between a GFP(+) strain and one with normal sequence at all three loci have been used to measure the interval to the centromere and to show that GFP experiences gene conversion with this system. Since any gene expressed in meiosis is silenced in Neurospora if hemizygous, any of our GFP(+) strains can be used as a quick screen to determine if a gene deleted by the Neurospora Genome Project is involved in crossing over or gene conversion.     

Strong positive selection and recombination drive the antigenic variation of the PilE protein of the human pathogen Neisseria meningitidis

The PilE protein is the major component of the Neisseria meningitidis pilus, which is encoded by the pilE/pilS locus that includes an expressed gene and eight homologous silent fragments. The silent gene fragments have been shown to recombine through gene conversion with the expressed gene and thereby provide a means by which novel antigenic variants of the PilE protein can be generated. We have analyzed the evolutionary rate of the pilE gene using the nucleotide sequence of two complete pilE/pilS loci. The very high rate of evolution displayed by the PilE protein appears driven by both recombination and positive selection. Within the semivariable region of the pilE and pilS genes, recombination appears to occur within multiple small sequence blocks that lie between conserved sequence elements. Within the hypervariable region, positive selection was identified from comparison of the silent and expressed genes. The unusual gene conversion mechanism that operates at the pilE/pilS locus is a strategy employed by N. meningitidis to enhance mutation of certain regions of the PilE protein. The silent copies of the gene effectively allow "parallelized" evolution of pilE, thus enabling the encoded protein to rapidly explore a large area of sequence space in an effort to find novel antigenic variants.     

Recombination, gene conversion, and identity-by-descent at three loci

We investigate the probabilities of identity-by-descent at three loci in order to find a signature which differentiates between the two types of crossing over events: recombination and gene conversion. We use a Markov chain to model coalescence, recombination, gene conversion and mutation in a sample of size two. Using numerical analysis, we calculate the total probability of identity-by-descent at the three loci, and partition these probabilities based on a partial ordering of coalescent events at the three loci. We use these results to compute the probabilities of four different patterns of conditional identity and non-identity at the three loci under recombination and gene conversion. Although recombination and gene conversion do make different predictions, the differences are not likely to be useful in distinguishing between them using three locus patterns between pairs of DNA sequences. This implies that measures of genetic identity in larger samples will be needed to distinguish between gene conversion and recombination.     

Possible Role of Natural Selection in the Formation of Tandem-Repetitive Noncoding DNA

A simulation model of sequence-dependent amplification, unequal crossing over and mutation is analyzed. This model predicts the spontaneous formation of tandem-repetitive patterns of noncoding DNA from arbitrary sequences for a wide range of parameter values. Natural selection is found to play an essential role in this self-organizing process. Natural selection which is modeled as a mechanism for controlling the length of a nucleotide string but not the sequence itself favors the formation of tandem-repetitive structures. Two measures of sequence heterogeneity, inter-repeat variability and repeat length, are analyzed in detail. For fixed mutation rate, both inter-repeat variability and repeat length are found to increase with decreasing rates of (unequal) crossing over. The results are compared with data on micro-, mini- and satellite DNAs. The properties of minisatellites and satellite DNAs resemble the simulated structures very closely. This suggests that unequal crossing over is a dominant long-range ordering force which keeps these arrays homogeneous even in regions of very low recombination rates, such as at satellite DNA loci. Our analysis also indicates that in regions of low rates of (unequal) crossing over, inter-repeat variability is maintained at a low level at the expense of much larger repeat units (multimeric repeats), which are characteristic of satellite DNA. In contrast, the microsatellite data do not fit the proposed model well, suggesting that unequal crossing over does not act on these very short tandem arrays.     

Recombinational and mutagenic repair of psoralen interstrand cross-links in Saccharomyces cerevisiae

Psoralen photoreacts with DNA to form interstrand cross-links, which can be repaired by both nonmutagenic nucleotide excision repair and recombinational repair pathways and by mutagenic pathways. In the yeast Saccharomyces cerevisiae, psoralen cross-links are processed by nucleotide excision repair to form double-strand breaks (DSBs). In yeast, DSBs are repaired primarily by homologous recombination, predicting that cross-link and DSB repair should induce similar recombination end points. We compared psoralen cross-link, psoralen monoadduct, and DSB repair using plasmid substrates with site-specific lesions and measured the patterns of gene conversion, crossing over, and targeted mutation. Psoralen cross-links induced both recombination and mutations, whereas DSBs induced only recombination, and monoadducts were neither recombinogenic nor mutagenic. Although the cross-link- and DSB-induced patterns of plasmid integration and gene conversion were similar in most respects, they showed opposite asymmetries in their unidirectional conversion tracts: primarily upstream from the damage site for cross-links but downstream for DSBs. Cross-links induced targeted mutations in 5% of the repaired plasmids; all were base substitutions, primarily T --> C transitions. The major pathway of psoralen cross-link repair in yeast is error-free and involves the formation of DSB intermediates followed by homologous recombination. A fraction of the cross-links enter an error-prone pathway, resulting in mutations at the damage site.     

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.     

Extrachromosomal and chromosomal gene conversion in mammalian cells.

We constructed substrates to study gene conversion in mammalian cells specifically without the complication of reciprocal recombination events. These substrates contain both an insertion mutation of the neomycin resistance gene (neoX) and an internal, homologous fragment of the neo gene (neo-526), such that gene conversion from neo-526 to neoX restores a functional neo gene. Although two reciprocal recombination events can also produce an intact neo gene, these double recombination events occur much less frequently that gene conversion in mammalian cells, We used our substrates to characterize extrachromosomal gene conversion in recombination-deficient bacteria and in monkey COS cells. Chromosomal recombination was also studied after stable integration of these substrates into the genome of mouse 3T6 cells. All extrachromosomal and chromosomal recombination events analyzed in mammalian cells resulted from gene conversion. Chromosomal gene conversion events occurred at frequencies of about 10(-6) per cell generation and restored a functional neo gene without overall effects on sequence organization.     

Gene conversion during vector insertion in embryonic stem cells.

Recombination of an insertion vector into its chromosomal homologue is a conservative event in that both the chromosomal and the vector sequences are preserved. However, gene conversion may accompany homologous recombination of an insertion vector. To examine gene conversion in more detail we have determined the targeting frequencies and the structure of the recombinant alleles generated with a series of vectors which target the hprt gene in embryonic stem cells. We demonstrate that gene conversion of the introduced mutation does not significantly limit homologous recombination and that gene conversion occurs without a sequence specific bias for five different mutations. The frequency of the loss of a vector mutation and the gain of a chromosomal sequence is inversely proportional to the distance between the vector mutation and the double-strand break. The loss of a chromosomal sequence and the gain of a vector mutation occurs at a low frequency.     

Genome rearrangement in top3 mutants of Saccharomyces cerevisiae requires a functional RAD1 excision repair gene.

Saccharomyces cerevisiae cells that are mutated at TOP3, a gene that encodes a protein homologous to bacterial type I topoisomerases, have a variety of defects, including reduced growth rate, altered gene expression, blocked sporulation, and elevated rates of mitotic recombination at several loci. The rate of ectopic recombination between two unlinked, homologous loci, SAM1 and SAM2, is sixfold higher in cells containing a top3 null mutation than in wild-type cells. Mutations in either of the two other known topoisomerase genes in S. cerevisiae, TOP1 and TOP2, do not affect the rate of recombination between the SAM genes. The top3 mutation also changes the distribution of recombination events between the SAM genes, leading to the appearance of novel deletion-insertion events in which conversion tracts extend beyond the coding sequence, replacing the DNA flanking the 3' end of one SAM gene with nonhomologous DNA flanking the 3' end of the other. The effects of the top3 null mutation on recombination are dependent on the presence of an intact RAD1 excision repair gene, because both the rate of SAM ectopic gene conversion and the conversion tract length were reduced in rad1 top3 mutant cells compared with top3 mutants. These results suggest that a RAD1-dependent function is involved in the processing of damaged DNA that results from the loss of Top3 activity, targeting such DNA for repair by recombination.     

Molecular models accounting for the gene conversion reactions mediating gonococcal pilin antigenic variation

The pilus antigenic variation (Av) system of Neisseria gonorrhoeae is one of several high-frequency variation systems that utilize gene conversion to switch between numerous forms of an antigen on the cell surface. We have tested three predictions of the first models that explain the movement of DNA during pilin Av: (i) Av requires two recombinations at short regions of identity, (ii) circular intermediates exist that carry pilE/pilS hybrid loci and (iii) these pilE/pilS hybrid loci target the pilS sequences to a recipient pilE gene. We confirm that normal pilin Av utilizes recombination at very short regions of DNA sequence identity and that these recombination events can occur independent of homologous recombination functions. We have isolated covalently closed circular DNA molecules carrying hybrid pilin loci, but propose that an alternative hybrid molecule is the intermediate of pilin Av. Our most striking finding is that transformation of isolated pilE/pilS hybrid loci targets the pilS sequences of the hybrid to a recipient pilE at frequencies much higher than normal recombination frequencies. These results show that the different steps of a model that explains pilin Av can be separately tested to support the validity of these novel models that account for the high-frequency gene conversions that mediate pilin Av.     

Double-Stranded Gap Repair of DNA by Gene Conversion in Escherichia Coli

We demonstrated repair of a double-stranded DNA gap through gene conversion by a homologous DNA sequence in Escherichia coli. We made a double-stranded gap in one of the two regions of homology in an inverted orientation on a plasmid DNA molecule and introduced it into an E. coli strain which has the RecE system of recombination (genotype; sbcA23 recB21 recC22). We detected repair products by genetic selection. The repair products were those expected by the double-strand-gap repair model. Gene conversion was frequently accompanied by crossing over of the flanking sequences as in eukaryotes. This double-strand gap repair mechanism can explain plasmid recombination in the absence of an artificial double-stranded break reported in a companion study by Yamamoto et al.     

A nuclear mutation defective in mitochondrial recombination in yeast.

Homologous recombination (crossing over and gene conversion) is generally essential for heritage and DNA repair, and occasionally causes DNA aberrations, in nuclei of eukaryotes. However, little is known about the roles of homologous recombination in the inheritance and stability of mitochondrial DNA which is continuously damaged by reactive oxygen species, by-products of respiration. Here, we report the first example of a nuclear recessive mutation which suggests an essential role for homologous recombination in the stable inheritance of mitochondrial DNA. For the detection of this class of mutants, we devised a novel procedure, 'mitochondrial crossing in haploid', which has enabled us to examine many mutant clones. Using this procedure, we examined mutants of Saccharomyces cerevisiae that showed an elevated UV induction of respiration-deficient mutations. We obtained a mutant that was defective in both the omega-intron homing and Endo.SceI-induced homologous gene conversion. We found that the mutant cells are temperature sensitive in the maintenance of mitochondrial DNA. A tetrad analysis indicated that elevated UV induction of respiration-deficient mutations, recombination deficiency and temperature sensitivity are all caused by a single nuclear mutation (mhr1) on chromosome XII. The pleiotropic characteristics of the mutant suggest an essential role for the MHR1 gene in DNA repair, recombination and the maintenance of DNA in mitochondria.     

Testing models of selection and demography in Drosophila simulans

We analyze patterns of nucleotide variability at 15 X-linked loci and 14 autosomal loci from a North American population of Drosophila simulans. We show that there is significantly more linkage disequilibrium on the X chromosome than on chromosome arm 3R and much more linkage disequilibrium on both chromosomes than expected from estimates of recombination rates, mutation rates, and levels of diversity. To explore what types of evolutionary models might explain this observation, we examine a model of recurrent, nonoverlapping selective sweeps and a model of a recent drastic bottleneck (e.g., founder event) in the demographic history of North American populations of D. simulans. The simple sweep model is not consistent with the observed patterns of linkage disequilibrium nor with the observed frequencies of segregating mutations. Under a restricted range of parameter values, a simple bottleneck model is consistent with multiple facets of the data. While our results do not exclude some influence of selection on X vs. autosome variability levels, they suggest that demography alone may account for patterns of linkage disequilibrium and the frequency spectrum of segregating mutations in this population of D. simulans.     

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.     

Different pathways of homologous recombination are used for the repair of double-strand breaks within tandemly arranged sequences in the plant genome

Different DNA repair pathways that use homologous sequences in close proximity to genomic double-strand breaks (DSBs) result in either an internal deletion or a gene conversion. We determined the efficiency of these pathways in somatic plant cells of transgenic Arabidopsis lines by monitoring the restoration of the beta-glucuronidase (GUS) marker gene. The transgenes contain a recognition site for the restriction endonuclease I-SceI either between direct GUS repeats to detect deletion formation (DGU.US), or within the GUS gene to detect gene conversion using a nearby donor sequence in direct or inverted orientation (DU.GUS and IU.GUS). Without expression of I-SceI, the frequency of homologous recombination (HR) was low and similar for all three constructs. By crossing the different lines with an I-SceI expressing line, DSB repair was induced, and resulted in one to two orders of magnitude higher recombination frequency. The frequencies obtained with the DGU.US construct were about five times higher than those obtained with DU.GUS and IU.GUS, irrespective of the orientation of the donor sequence. Our results indicate that recombination associated with deletions is the most efficient pathway of homologous DSB repair in plants. However, DSB-induced gene conversion seems to be frequent enough to play a significant role in the evolution of tandemly arranged gene families like resistance genes.     

Nucleotide sequence of the chi recombinational hot spot chi +D in bacteriophage lambda.

Chi sites in bacteriophage lambda stimulate recombination promoted by the RecBC pathway of Escherichia coli. Mutations which create these sites occur at four widely separated loci in lambda. We report here the nucleotide sequence surrounding the site of one of these loci, chi D, located near the S gene. The mutations creating the active Chi site, designated chi +D, are transversions from CG to AT. This mutation, like the chi +B and chi +C mutations previously analyzed, leads to a nucleotide sequence common to all three active chi sites.     

Xrs2, a DNA Repair Gene of Saccharomyces Cerevisiae, Is Needed for Meiotic Recombination

The XRS2 gene of Saccharomyces cerevisiae has been previously identified as a DNA repair gene. In this communication, we show that XRS2 also encodes an essential meiotic function. Spore inviability of xrs2 strains is rescued by a spo13 mutation, but meiotic recombination (both gene conversion and crossing over) is highly depressed in spo13 xrs2 diploids. The xrs2 mutation suppresses spore inviability of a spo13 rad52 strain suggesting that XRS2 acts prior to RAD52 in the meiotic recombination pathway. In agreement with the genetic data, meiosis-specific double-strand breaks at the ARG4 meiotic recombination hotspot are not detected in xrs2 strains. Despite its effects on meiotic recombination, the xrs2 mutation does not prevent mitotic recombination events, including homologous integration of linear DNA, mating-type switching and radiation-induced gene conversion. Moreover, xrs2 strains display a mitotic hyper-rec phenotype. Haploid xrs2 cells fail to carry out G2-repair of gamma-induced lesions, whereas xrs2 diploids are able to perform some diploid-specific repair of these lesions. Meiotic and mitotic phenotypes of xrs2 cells are very similar to those of rad50 cells suggesting that XRS2 is involved in homologous recombination in a way analogous to that of RAD50.