A coalescence approach to gene conversion

In this paper we develop a coalescent model with intralocus gene conversion. Such models are of increasing importance in the analysis of intralocus variability and linkage disequilibrium. We derive the distribution of the waiting time until a gene conversion event occurs in a sample in terms of the distribution of the length of the transferred segment, zeta. We do not assume any specific form of the distribution of zeta. Further, given that a gene conversion event occurs we find the distribution of (sigma, tau), the end points of the transferred segment and derive results on correlations between local trees in positions chi(1) and chi(2). Among other results we show that the correlation between the branch lengths of two local trees in the coalescent with gene conversion (and no recombination) decreases toward a nonzero constant when the distance between chi(1) and chi(2) increases. Finally, we show that a model including both recombination and gene conversion might account for the lack of intralocus associations found in, e.g., Drosophila melanogaster.     

The coalescent with gene conversion

In this article we develop a coalescent model with intralocus gene conversion. The distribution of the tract length is geometric in concordance with results published in the literature. We derive a simulation scheme and deduce a number of analytical results for this coalescent with gene conversion. We compare patterns of variability in samples simulated according to the coalescent with recombination with similar patterns simulated according to the coalescent with gene conversion alone. Further, an expression for the expected number of topology shifts in a sample of present-day sequences caused by gene conversion events is derived.     

Assessing the pathogenic potential of the V(D)J recombinase by interlocus immunoglobulin light-chain gene rearrangement.

Chromosomal translocations involving antigen receptor genes and oncogenes have been observed in several forms of lymphoid malignancy. Observations of their lymphocyte-restricted occurrence and a molecular analysis of some translocation breakpoints have suggested that some of these rearrangements are generated by V(D)J recombinase activity. However, a direct correlation between this activity and the generation of such rearrangements has never been established. In addition, because these aberrant rearrangements are usually detected only after a tumor has been formed, the frequency with which the recombinase machinery generates translocations has never been assessed directly. To approach these issues, immunoglobulin light-chain gene rearrangements were induced in pre-B cells transformed by temperature-sensitive mutants of Abelson murine leukemia virus and PCR was used to identify interlocus recombinants. Vlambda Jkappa and Vkappa Jlambda rearrangements as well as signal joints resulting from the recombination of Vlambda and Jkappa coding elements were recovered and were found to be similar in structure to conventional intrachromosomal joints. Because these products were detected only when the cells were undergoing active intralocus rearrangement, they provide direct evidence that translocations can be generated by the V(D)J recombinase machinery. Dilution analyses revealed that interlocus rearrangements occur about 1,000 times less frequently than conventional intralocus rearrangements. Considering the large numbers of lymphocytes generated throughout life, aberrant rearrangements generated by the V(D)J recombinase may be relatively common.     

Distinguishing recombination and intragenic gene conversion by linkage disequilibrium patterns

Deterministic theory suggests that reciprocal recombination and intragenic, interallelic conversion have different effects on the linkage disequilibrium between a pair of genetic markers. Under a model of reciprocal recombination, the decay rate of linkage disequilibrium depends on the distance between the two markers, while under conversion the decay rate is independent of this distance, provided that conversion tracts are short. A population genetic three-locus model provides a function Q of two-locus linkage disequilibria. Viewed as a random variable, Q is the basis for a test of the relative impact of conversion and recombination. This test requires haplotype frequency data of a sufficiently variable three-locus system. One of the few examples currently available is data from the Human Leukocyte Antigen (HLA) class I genes of three Amerindian populations. We find that conversion may have played a dominant role in shaping haplotype patterns over short stretches of DNA, whereas reciprocal recombination may have played a greater role over longer stretches of DNA. However, in order to draw firm conclusions more independent data are necessary.     

Evidence against reciprocal recombination as the basis for tuf gene conversion in Salmonella enterica serovar Typhimurium

The duplicate tuf genes on the Salmonella enterica serovar Typhimurium chromosome co-evolve by a RecA-, RecB-dependent gene conversion mechanism. Gene conversion is defined as a non-reciprocal transfer of genetic information. However, in a replicating bacterial chromosome there is a possibility that a reciprocal genetic exchange between different tuf genes sitting on sister chromosomes could result in "apparent" gene conversion. We asked whether the major mechanism of tuf gene conversion was classical or apparent. We devised a genetic selection that allowed us to isolate and examine both expected products from a reciprocal recombination event between the tuf genes. Using this selection we tested within individual cultures for a correlation in the frequency of jackpots as expected if recombination were reciprocal. We found no correlation, either in the frequency of each type of recombinant product, or in the DNA sequences of the products resulting from each recombination event. We conclude that the evidence argues in favor of a non-reciprocal gene conversion mechanism as the basis for tuf gene co-evolution.     

An Attempt to Measure the Patchwork Pattern Observed Among Alleles at Major Histocompatibility Complex Loci

By means of simulations and DNA sequence analyses, standardized identity excess (a measure of linkage disequilibrium) between segregating nucleotide sites was studied as an effort to quantify the patchwork pattern among alleles of the major histocompatibility complex loci. It was found that the pattern under selective neutrality, and/or no intralocus recombination does not fit the observed pattern based on DNA sequences. However, the intensity and type of selection and the rate of recombination are difficult to estimate by comparing simulation results with the observed pattern. Received: 10 December 1999 / Accepted: 2 March 2000     

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.     

Rapid postradiation recovery of yeasts: the relationship to gamma-induced reciprocal mitotic recombination and gene conversion

gamma-induced reciprocal mitotic recombination and gene conversion have been studied under conditions inhibiting "rapid" postirradiation recovery of diploid yeast Saccharomyces cerevisiae. It turned out that, if the first postirradiation cell division occurs at higher KCl concentrations ("rapid" postirradiation recovery is inhibited), the frequency of mitotic reciprocal recombination within the gene ADE2-centromere region decreases. Keeping of irradiated cells (in the G1 phase of the cell division cycle) in water at 28 degrees C prior to plating on the selective agar containing 1.5 M KCl leads to smaller frequency of gene conversion lys2-25/lys2-22----Lys+, as compared with that for the cells immediately plated on the selective agar. Correlation has been found between the coefficient of gene conversion frequency decrease, due to postirradiation keeping in water, and "rapid" recovery efficiency. Interpretation of the data is based on the hypothesis that recombination repair of DNA double-strand breaks induced by ionizing radiation is responsible for "rapid" postirradiation recovery.     

Immunoglobulin gene ''remnant'' DNA--implications for antibody gene recombination.

Many immunoglobulin (Ig)-producing cells retain the DNA that separates Ig variable (V) and constant (C) region genes in the germline. This "remnant" DNA must be moved during the recombination process that joins V and C genes via a joining (J) segment. We have analyzed remnant DNAs in several Ig-producing cell lines. The nucleotide sequences of kappa (kappa) light chain remnant DNAs indicate close relationships to V-J joining. We find fused V kappa and J kappa recognition sequences in five remnant DNAs, suggesting reciprocal relationships to the fused V kappa and J kappa segments produced by V-J joining. However, of sixteen plasmacytoma remnant DNAs analyzed, all involve only recombination with J kappa l. Thus, in most cell lines, remnant DNAs are not directly reciprocal to recombined kappa-genes. On the other hand, our analyses of some myelomas do indicate indirect relationships between remnant DNAs and kappa-genes. Our results suggest that multiple steps of DNA recombination occur during Ig-gene rearrangement. Because remnant DNA joining sites do not exhibit the flexibility that has been observed in Ig-gene V-J joining, our findings also suggest that the joining mechanism may involve endonuclease, exonuclease and ligase activities.     

The non-reciprocality of organelle gene recombination in Chlamydomonas reinhardtii and Saccharomyces cerevisiae: some new observations and a restatement of some old problems.

Summary Organelle recombinant genotype frequencies, derived from analysis of individual mitotic zygote clones of Chlamydomonas reinhardtii and Saccharomyces cerevisiae, were subjected to two types of statistical tests in an attempt to detect the occurrence of reciprocal recombination: (i) calculation of correlation coefficients for the frequencies of two recombinant genotypes (reciprocal or non-reciprocal pairs) within individual zygote clones, and (ii) application of the chi-square test for independence to the frequencies of zygotes yielding one or the other, neither, or both of a given recombinant pair.Applying test (i), the strongest correlations are found for non-reciprocal rather than reciprocal pairs. When the data are analyzed by method (ii), some reciprocal as well as non-reciprocal pairs appear to be produced concurrently in zygote clones. However, such deviations from independence are greatest for non-reciprocal pairs. These tests yield comparable results for yeast mitochondrial and Chlamydomonas chloroplast gene recombination, and provide no convincing evidence for reciprocal genetic exchange. Explanations for the observed lack of reciprocality are discussed with reference both to our present understanding of the molecular events responsible for genetic recombination, and to the problems which may be unique to the analysis of organelle gene recombination.     

T-cell receptor gene rearrangement in primary tumors: effect of genetic background and inducing agent

The status of T-cell receptor beta and gamma genes has been assessed in a series of primary tumors induced by a chemical carcinogen or by gamma-irradiation using two inbred strains of mice. It appears that these well-characterized regimens of carcinogenesis yield T-cell tumors showing gene rearrangements consistent with a clonal origin of the tumors. Individual rearranged bands seem to represent orthodox, intralocus recombination events. A variety of rearrangement phenotypes are observed, most strikingly for the gamma genes, and differences in the degree of T-cell receptor gene rearrangements observed can be categorized according to the inducing agent and to the genetic background of the mice, with the implication that premalignant thymocytes have been captured in different stages of T-cell development. Additionally, primary tumors were shown to express significant levels of mature beta gene mRNA.     

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.     

Reciprocal control of catalysis by the tyrosine recombinases XerC and XerD: an enzymatic switch in site-specific recombination

In Xer site-specific recombination, sequential DNA strand exchange reactions are catalyzed by a heterotetrameric complex composed of two recombinases, XerC and XerD. It is demonstrated that XerC and XerD catalytic activity is controlled by an interaction involving the C-terminal end of each protein (the donor region) and an internal region close to the active site (the acceptor region). Mutations in these regions reciprocally alter the relative activity of XerC and XerD, with their combination producing synergistic effects on catalysis. The data support a model in which C-terminal intersubunit interactions contribute to coupled protein-DNA conformational changes that lead to sequential activation and reciprocal inhibition of pairs of active sites in the recombinase tetramer during recombination.     

Recombination hot spot in the human beta-globin gene cluster: meiotic recombination of human DNA fragments in Saccharomyces cerevisiae.

We describe a novel system for the analysis of sequence-specific meiotic recombination in Saccharomyces cerevisiae. A comparison of three adjacent restriction fragments from the human beta-globin locus revealed that one of them, previously hypothesized to contain a relative hot spot for genetic recombination, engages in reciprocal exchange during yeast meiosis significantly more frequently than either of the other two fragments. Removal of the longest of four potential Z-DNA-forming regions from this fragment does not affect the high frequency of genetic recombination.     

Preaching about the converted: how meiotic gene conversion influences genomic diversity

Meiotic crossover (CO) recombination involves a reciprocal exchange between homologous chromosomes. COs are often associated with gene conversion at the exchange site where genetic information is unidirectionally transferred from one chromosome to the other. COs and independent assortment of homologous chromosomes contribute significantly to the promotion of genomic diversity. What has not been appreciated is the contribution of another product of meiotic recombination, noncrossovers (NCOs), which result in gene conversion without exchange of flanking markers. Here, we review our comprehensive analysis of recombination at a highly polymorphic mouse hotspot. We found that NCOs make up ～90% of recombination events. Preferential recombination initiation on one chromosome allowed us to estimate the contribution of CO and NCO gene conversion to transmission distortion, a deviation from Mendelian inheritance in the population. While NCO gene conversion tracts are shorter, and thus have a more punctate effect, their higher frequency translates into an approximately two-fold greater contribution than COs to gene conversion-based allelic shuffling and transmission distortion. We discuss the potential impact of mammalian NCO characteristics on evolution and genomic diversity.     

Jaatha: a fast composite-likelihood approach to estimate demographic parameters

While information about a species’ demography is interesting in its own right, it is an absolute necessity for certain types of population genetic analyses. The most widely used methods to infer a species’ demographic history do not take intralocus recombination or recent divergence into account, and some methods take several weeks to converge. Here, we present Jaatha, a new composite‐likelihood method that does incorporate recent divergence and is also applicable when intralocus recombination rates are high. This new method estimates four demographic parameters. The accuracy of Jaatha is comparable to that of other currently available methods, although it is superior under certain conditions, especially when divergence is very recent. As a proof of concept, we apply this new method to estimate demographic parameters for two closely related wild tomato species, Solanum chilense and S. peruvianum. Our results indicate that these species likely diverged 1.44·N generations ago, where N is the effective population size of S. chilense, and that some introgression between these species continued after the divergence process initiated. Furthermore, S. peruvianum likely experienced a population expansion following speciation.     

Direct evidence for intrastrand DNA inversion of kappa immunoglobulin gene segments in two murine plasmacytomas.

The products of kappa immunoglobulin gene recombination have been characterized in two murine plasmacytomas to examine the relationship between V-J products and reciprocal elements. By cloning, sequencing, hybridization, and application of the polymerase chain reaction, we have established the direct relationship of the kappa recombination products in these cells. The results provide stronger support for the intrastrand mechanism of kappa gene recombination as well as demonstrating a role for secondary, corrective recombinations.     

T cell receptor beta gene sequences in the circular DNA of thymocyte nuclei: direct evidence for intramolecular DNA deletion in V-D-J joining

We have identified circular DNA containing T cell receptor (TCR) beta gene sequences in mouse thymocytes, thereby providing direct evidence for the intramolecular DNA deletion model of V-D-J joining in TCR beta genes. Two types of excision products of V-D-J joining have been identified. Type I, a circular reciprocal recombinant of normal V-D or D-J joining, contains a 7mer-7mer head-to-head structure expected from an excised product of normal V-D or D-J joining. Type II contains a D beta 2-J beta 1 structure on the circular DNA; the recombination event producing this molecule occurs between an upstream J and a downstream D segment, probably leaving the reciprocal 7mer-7mer structure on the chromosome. Some type I molecules seem to represent excision products of secondary joining after formation of the first D-J or V-D-J structure. The recombination mechanism that generates the circular DNA is discussed.     

Targeted gene knockout reveals a role in meiotic recombination for ZHP-3, a Zip3-related protein in Caenorhabditis elegans

The meiotically expressed Zip3 protein is found conserved from Saccharomyces cerevisiae to humans. In baker's yeast, Zip3p has been implicated in synaptonemal complex (SC) formation, while little is known about the protein's function in multicellular organisms. We report here the successful targeted gene disruption of zhp-3 (K02B12.8), the ZIP3 homolog in the nematode Caenorhabditis elegans. Homozygous zhp-3 knockout worms show normal homologue pairing and SC formation. Also, the timing of appearance and the nuclear localization of the recombination protein Rad-51 seem normal in these animals, suggesting proper initiation of meiotic recombination by DNA double-strand breaks. However, the occurrence of univalents during diplotene indicates that C. elegans ZHP-3 protein is essential for reciprocal recombination between homologous chromosomes and thus chiasma formation. In the absence of ZHP-3, reciprocal recombination is abolished and double-strand breaks seem to be repaired via alternative pathways, leading to achiasmatic chromosomes and the occurrence of univalents during meiosis I. Green fluorescent protein-tagged C. elegans ZHP-3 forms lines between synapsed chromosomes and requires the SC for its proper localization.     

A genetic linkage map of rye (Secale cereale L.)

 A genetic linkage map of rye composed of 91 loci (88 RFLP, two morphological and one isozyme markers) has been developed using two reciprocal crosses. The RFLP loci covering all seven chromosomes were detected by a selection of rye, wheat, barley and oat cDNA and genomic DNA probes. The level of polymorphism was dependent on the source of the clones, with a ranking of rye>wheat>barley>oat. Distorted segregations were detected in linkage groups of chromosomes 1R, 4R, 5R and 7R. When the recombination of the two reciprocal crosses was compared, no systematic increase or decrease in one or the other direction was observed suggesting that a combination of populations of reciprocal crosses is possible. Received: 5 August 1997/Accepted: 2 September 1997