Recombination hotspots as a point process

The variation of the recombination rate along chromosomal DNA is one of the important determinants of the patterns of linkage disequilibrium. A number of inferential methods have been developed which estimate the recombination rate and its variation from population genetic data. The majority of these methods are based on modelling the genealogical process underlying a sample of DNA sequences and thus explicitly include a model of the demographic process. Here we propose a different inferential procedure based on a previously introduced framework where recombination is modelled as a point process along a DNA sequence. The approach infers regions containing putative hotspots based on the inferred minimum number of recombination events; it thus depends only indirectly on the underlying population demography. A Poisson point process model with local rates is then used to infer patterns of recombination rate estimation in a fully Bayesian framework. We illustrate this new approach by applying it to several population genetic datasets, including a region with an experimentally confirmed recombination hotspot.     

The Contribution of Alu Elements to Mutagenic DNA Double-Strand Break Repair

Alu elements make up the largest family of human mobile elements, numbering 1.1 million copies and comprising 11% of the human genome. As a consequence of evolution and genetic drift, Alu elements of various sequence divergence exist throughout the human genome. Alu/Alu recombination has been shown to cause approximately 0.5% of new human genetic diseases and contribute to extensive genomic structural variation. To begin understanding the molecular mechanisms leading to these rearrangements in mammalian cells, we constructed Alu/Alu recombination reporter cell lines containing Alu elements ranging in sequence divergence from 0%-30% that allow detection of both Alu/Alu recombination and large non-homologous end joining (NHEJ) deletions that range from 1.0 to 1.9 kb in size. Introduction of as little as 0.7% sequence divergence between Alu elements resulted in a significant reduction in recombination, which indicates even small degrees of sequence divergence reduce the efficiency of homology-directed DNA double-strand break (DSB) repair. Further reduction in recombination was observed in a sequence divergence-dependent manner for diverged Alu/Alu recombination constructs with up to 10% sequence divergence. With greater levels of sequence divergence (15%-30%), we observed a significant increase in DSB repair due to a shift from Alu/Alu recombination to variable-length NHEJ which removes sequence between the two Alu elements. This increase in NHEJ deletions depends on the presence of Alu sequence homeology (similar but not identical sequences). Analysis of recombination products revealed that Alu/Alu recombination junctions occur more frequently in the first 100 bp of the Alu element within our reporter assay, just as they do in genomic Alu/Alu recombination events. This is the first extensive study characterizing the influence of Alu element sequence divergence on DNA repair, which will inform predictions regarding the effect of Alu element sequence divergence on both the rate and nature of DNA repair events.     

Investigating the mechanism of chromosomal deletion: characterization of 39 deletion breakpoints in introns 47 and 48 of the human dystrophin gene

The region of the dystrophin gene containing introns 45-50 is characterized by a high rate of recombination events that give rise to large deletions causing dystrophinopathy. The nucleotide sequence of this intronic region has recently been released in GenBank. With the aim of further understanding the mechanism favoring the occurrence of these deletions, we have characterized the distribution of introns 47 and 48 deletion endpoints in 39 dystrophinopathy patients. In 14 of these patients we were able to sequence the break junction. On these sequences we were able to identify several intronic motifs that could predispose to DNA double-strand breaks. Our results, combined with other literature data, show that unequal homologous recombination is a very poorly represented event in the dystrophin gene, whereas junction features are suggestive of a model of recombination in which DNA double-strand breaks are incorrectly repaired by a nonhomologous end-joining mechanism. The correlation among recombination rate, deletion frequency, and percentage of repetitive elements is discussed.     

Nonrandom localization of recombination events in human alpha satellite repeat unit variants: implications for higher-order structural characteristics within centromeric heterochromatin.

Tandemly repeated DNA families appear to undergo concerted evolution, such that repeat units within a species have a higher degree of sequence similarity than repeat units from even closely related species. While intraspecies homogenization of repeat units can be explained satisfactorily by repeated rounds of genetic exchange processes such as unequal crossing over and/or gene conversion, the parameters controlling these processes remain largely unknown. Alpha satellite DNA is a noncoding tandemly repeated DNA family found at the centromeres of all human and primate chromosomes. We have used sequence analysis to investigate the molecular basis of 13 variant alpha satellite repeat units, allowing comparison of multiple independent recombination events in closely related DNA sequences. The distribution of these events within the 171-bp monomer is nonrandom and clusters in a distinct 20- to 25-bp region, suggesting possible effects of primary sequence and/or chromatin structure. The position of these recombination events may be associated with the location within the higher-order repeat unit of the binding site for the centromere-specific protein CENP-B. These studies have implications for the molecular nature of genetic recombination, mechanisms of concerted evolution, and higher-order structure of centromeric heterochromatin.     

Bounds on the minimum number of recombination events in a sample history

Recombination is an important evolutionary factor in many organisms, including humans, and understanding its effects is an important task facing geneticists. Detecting past recombination events is thus important; this article introduces statistics that give a lower bound on the number of recombination events in the history of a sample, on the basis of the patterns of variation in the sample DNA. Such lower bounds are appropriate, since many recombination events in the history are typically undetectable, so the true number of historical recombinations is unobtainable. The statistics can be calculated quickly by computer and improve upon the earlier bound of Hudson and Kaplan 1985. A method is developed to combine bounds on local regions in the data to produce more powerful improved bounds. The method is flexible to different models of recombination occurrence. The approach gives recombination event bounds between all pairs of sites, to help identify regions with more detectable recombinations, and these bounds can be viewed graphically. Under coalescent simulations, there is a substantial improvement over the earlier method (of up to a factor of 2) in the expected number of recombination events detected by one of the new minima, across a wide range of parameter values. The method is applied to data from a region within the lipoprotein lipase gene and the amount of detected recombination is substantially increased. Further, there is strong clustering of detected recombination events in an area near the center of the region. A program implementing these statistics, which was used for this article, is available from http://www.stats.ox.ac.uk/mathgen/programs.html.     

Sequence requirements for DNA rearrangements induced by the terminal repeat of herpes simplex virus type 1 KOS DNA.

We investigated the sequence requirements for the site-specific DNA cleavages and recombinational genome isomerization events driven by the terminal repeat or a sequence of herpes simplex virus type 1 KOS DNA by inserting a series of mutated a sequences into the thymidine kinase locus in the intact viral genome. Our results indicate that sequences located at both extremities of the a sequence contribute to these events. Deletions entering from the Ub side of the a sequence progressively reduced the frequency of DNA rearrangements, and further deletion of the internal DR2 repeat array had an additional inhibitory effect. This deletion series allowed us to map the pac1 site-specific DNA cleavage signal specifying the S-terminal cleavage to a sequence that is conserved among herpesvirus genomes. Constructs lacking this signal were unable to directly specify the S-terminal cleavage event but retained a reduced ability to give rise to S termini following recombination with intact a sequences. Deletions entering from the Uc side demonstrated that the copy of direct repeat 1 located adjacent to the Uc region plays an important role in the DNA rearrangements induced by the a sequence: mutants lacking this sequence displayed a reduced frequency of novel terminal and recombinational inversion fragments, and further deletions of the Uc region had a relatively minor additional effect. By using a construct in which site-specific cleavage was directed to heterologous DNA sequences, we found that the recombination events leading to genome segment inversion did not occur at the sites of DNA cleavage used by the cleavage-packaging machinery. This observation, coupled with the finding that completely nonoverlapping portions of the a sequence retained detectable recombinational activity, suggests that inter-a recombination does not occur by cleavage-ligation at a single specific site in herpes simplex virus type 1 strain KOS. The mutational sensitivity of the extremities of the a sequence leads us to hypothesize that the site-specific DNA breaks induced by the cleavage-packaging system stimulate the initiation of recombination.     

'Illegitimate' recombination events in polyoma-transformed rat cells

In the LPT line of polyoma (Py)-transformed rat cells, amplification of the integrated viral DNA and of cell nucleotide sequences flanking the viral integration site, can be induced either spontaneously or by treatment with carcinogens. We show here that the amplified DNA includes interspersed viral and cellular sequences generated by 'illegitimate' recombination events. Genomic libraries have been prepared in phage lambda vectors from LPT cells treated with the inducing agent mitomycin C and from untreated LPT cells. Four phages, including viral-cell DNA recombinants, have been isolated from these libraries. Sequencing through the recombination sites revealed the following characteristics: (i) The crossover points map at four different positions in the viral DNA and at four different positions in the flanking cell DNA. (ii) There are very short homologous sequences of 1, 2, or 4 bp, at the recombination sites. (iii) Aside from the exchanges between the viral and the cellular DNA, no further rearrangements occurred around the new viral-cellular DNA junctions. (iv) Next to the recombination sites, there are blocks of homopurine-homopyrimidine sequences, which may assume a structure that differs from the Watson-Crick double helix. (v) Clustered homologous sequence blocks of up to 10 bp are present less than 200 bp away from the recombination sites. These homologies are not in register. Based on these results, we propose a model that may account for these recombination events and, more generally, for recombination events that occur during gene amplification in mammalian cells.     

Inference of homologous recombination in bacteria using whole-genome sequences

Bacteria and archaea reproduce clonally, but sporadically import DNA into their chromosomes from other organisms. In many of these events, the imported DNA replaces an homologous segment in the recipient genome. Here we present a new method to reconstruct the history of recombination events that affected a given sample of bacterial genomes. We introduce a mathematical model that represents both the donor and the recipient of each DNA import as an ancestor of the genomes in the sample. The model represents a simplification of the previously described coalescent with gene conversion. We implement a Monte Carlo Markov chain algorithm to perform inference under this model from sequence data alignments and show that inference is feasible for whole-genome alignments through parallelization. Using simulated data, we demonstrate accurate and reliable identification of individual recombination events and global recombination rate parameters. We applied our approach to an alignment of 13 whole genomes from the Bacillus cereus group. We find, as expected from laboratory experiments, that the recombination rate is higher between closely related organisms and also that the genome contains several broad regions of elevated levels of recombination. Application of the method to the genomic data sets that are becoming available should reveal the evolutionary history and private lives of populations of bacteria and archaea. The methods described in this article have been implemented in a computer software package, ClonalOrigin, which is freely available from http://code.google.com/p/clonalorigin/.     

Development of a highly efficient gene targeting system for Fusarium graminearum using the disruption of a polyketide synthase gene as a visible marker

We cloned a polyketide synthase gene (pks12) from Fusarium graminearum, a devastating fungal pathogen of cereals. Transformation-mediated gene disruption led to an easily detectable albino phenotype of the disruptants. We used the disruption of the pks12 gene as a visible marker for transformation-mediated homologous recombination and optimized the transformation procedure to achieve a high rate of homologous recombination. In combination with the published genomic sequence data and the generation of expressed sequence tags (ESTs) for F. graminearum, this is a useful tool to investigate this important plant pathogen on a molecular level. Optimized transformation of F. graminearum resulted in at least 93% homologous recombination events when the homologous genomic DNA fragment in the vector had a size of approximately 800bp and was linearized in the middle. Using a genomic sequence of approximately 500bp in the transformation vector, 70% of the transformants still exhibited homologous recombination. On the contrary, no more than 10% homologous recombination events were observed when less than 400bp DNA fragments were used. We co-transformed F. graminearum with two different vectors. One vector harboured a DNA insert homologous to the pks12 gene, while the other vector consisted of the same vector backbone carrying the selection marker specific for F. graminearum. About 70% of the transformants had a disrupted pks12 gene, and all of these showed an integration of the second vector into the pks disruption vector. Therefore, the time-consuming construction of a single transformation vector can be avoided; furthermore, it is now easily feasible to express a gene construct at a defined and mutated genomic site.     

Extensive homologous recombination between introduced and native regulatory plastid DNA elements in transplastomic plants

Homologous recombination within plastids directs plastid genome transformation for foreign gene expression and study of plastid gene function. Though transgenes are generally efficiently targeted to their desired insertion site, unintended homologous recombination events have been observed during plastid transformation. To understand the nature and abundance of these recombination events, we analyzed transplastomic tobacco lines derived from three different plastid transformation vectors utilizing two different loci for foreign gene insertion. Two unintended recombinant plastid DNA species were formed from each regulatory plastid DNA element included in the transformation vector. Some of these recombinant DNA species accumulated to as much as 10–60% of the amount of the desired integrated transgenic sequence in T0 plants. Some of the recombinant DNA species undergo further, “secondary” recombination events, resulting in an even greater number of recombinant plastid DNA species. The abundance of novel recombinant DNA species was higher in T0 plants than in T1 progeny, indicating that the ancillary recombination events described here may have the greatest impact during selection and regeneration of transformants. A line of transplastomic tobacco was identified containing an antibiotic resistance gene unlinked from the intended transgene insertion as a result of an unintended recombination event, indicating that the homologous recombination events described here may hinder efficient recovery of plastid transformants containing the desired transgene. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Homologous and nonhomologous recombination in monkey cells.

Though recombinational events are important for the proper functioning of most cells, little is known about the frequency and mechanisms of recombination in mammalian cells. We have used simian virus 40 (SV40)-pBR322 hybrid plasmids constructed in vitro as substrates to detect and quantitate intramolecular homologous and nonhomologous recombination events in cultured monkey cells. Excision of wild-type or defective SV40 DNAs by recombination from these plasmids was scored by the viral plaque assay, in either the absence or the presence of DNA from a temperature-sensitive helper virus. Several independent products of homologous and nonhomologous recombination have been isolated and characterized at the DNA sequence level. We find that neither DNA replication of the recombination substrate nor SV40 large T antigen is essential for either homologous or nonhomologous recombination involving viral or pBR322 sequences.     

The role of the mismatch repair machinery in regulating mitotic and meiotic recombination between diverged sequences in yeast

Nonidentical recombination substrates recombine less efficiently than do identical substrates in yeast, and much of this inhibition can be attributed to action of the mismatch repair (MMR) machinery. In this study an intron-based inverted repeat assay system has been used to directly compare the rates of mitotic and meiotic recombination between pairs of 350-bp substrates varying from 82% to 100% in sequence identity. The recombination rate data indicate that sequence divergence impacts mitotic and meiotic recombination similarly, although subtle differences are evident. In addition to assessing recombination rates as a function of sequence divergence, the endpoints of mitotic and meiotic recombination events involving 94%-identical substrates were determined by DNA sequencing. The endpoint analysis indicates that the extent of meiotic heteroduplex DNA formed in a MMR-defective strain is 65% longer than that formed in a wild-type strain. These data are consistent with a model in which the MMR machinery interferes with the formation and/or extension of heteroduplex intermediates during recombination.     

Inverted duplication-transposition event in mammalian cells at an illegitimate recombination join.

Illegitimate recombination events in mammalian cells often contain extraneous nucleotides or filler DNA at the recombinant joins. The polyomavirus-transformed cell line 7axB has previously been found to contain 37 base pairs (bp) of filler DNA at one virus-host join of the single insert of integrated viral DNA (A. Hayday, H. E. Ruley, and M. Fried, J. Virol. 44:67-77, 1982). By using a synthetic oligomer of these 37 bp as a probe, we demonstrated that this filler DNA is an inverted duplication of a single-copy rat sequence found 650 bp upstream from this virus-host join. The other virus-host join appears to be the result of a simple illegitimate recombination event between viral and host sequences. This is the first identification of filler DNA as a transposed copy of a chromosomal sequence. The relevance of the recombination events studied to cellular rearrangements and viral integration is discussed.     

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.     

重复序列及异染色质化在植物性染色体重组抑制中的作用

重组抑制是植物性染色体由常染色体进化而来的前提条件,性别决定位点区域发生的重组抑制使早期的性染色体发生了退化和分化。研究表明,重组抑制的产生和染色体上一系列行为的发生有着密切的关系,如重复序列的累积、异染色质化及DNA的甲基化。转座因子和卫星DNA等重复序列的累积使早期植物性染色体形态和分子结构发生了分化,同时还导致性染色体的异染色质化,抑制了性染色体间的重组的发生。文章综述了这一领域的进展,并对DNA甲基化在植物性染色体重组抑制形成过程中可能的作用进行了简要分析。     

Products and implied mechanism of H chain switch recombination

The Ig H chain switch is a DNA recombination event. The recombination occurs between two or more switch regions, areas of tandem sequence duplication that lie upstream of the corresponding H chain C region genes. We have determined the DNA sequence at four recombination sites in three molecularly cloned, rearranged switch regions. All eight donor and recipient recombination sites are at the common pentamers GGGGT, GAGCT, and GGTGG. One of the switch recombination events is an inversion of S gamma 3 sequences. Another of the recombinational events is an internal S gamma 1 deletion, which may be switch enzyme mediated. These results, together with other switch recombination site sequences, suggest that switch recombination is mediated by cutting enzymes with modest specificity and religation enzymes with no specificity.     

Transfer of chimeric plasmids among Salmonella typhimurium strains by P22 transduction

Salmonella typhimurium bacteriophage P22 transduced plasmids having P22 sequences inserted in the vector pBR322 with high frequency. Analysis of the structure of the transducing particle DNA and the transduced plasmids indicates that this plasmid transduction involves two homologous recombination events. In the donor cell, a single recombination between the phage and the homologous sequences on the plasmid inserted the plasmid into the phage chromosome, which was then packaged by headfuls into P22 particles. The transducing particle DNA contained duplications of the region of homology flanking the integrated plasmid vector sequences and lacked some phage genes. When these defective phage genomes containing the inserted plasmid infected a recipient cell, recombination between the duplicated regions regenerated the plasmid. A useful consequence of this sequence of events was that genetic markers in the region of homology were readily transferred from phage to plasmid. Plasmid transduction required homology between the phage and the plasmid, but did not depend on the presence of any specific P22 sequence in the plasmid. When the infecting P22 carried a DNA sequence homologous to the ampicillin resistance region of pBR322, the vector plasmid having no P22 insert could be transduced. P22-mediated transduction is a useful way to transfer chimeric plasmids, since most S. typhimurium strains are poorly transformed by plasmid DNA.     

A Novel Computational Method Identifies Intra- and Inter-Species Recombination Events in Staphylococcus aureus and Streptococcus pneumoniae

Advances in high-throughput DNA sequencing technologies have determined an explosion in the number of sequenced bacterial genomes. Comparative sequence analysis frequently reveals evidences of homologous recombination occurring with different mechanisms and rates in different species, but the large-scale use of computational methods to identify recombination events is hampered by their high computational costs. Here, we propose a new method to identify recombination events in large datasets of whole genome sequences. Using a filtering procedure of the gene conservation profiles of a test genome against a panel of strains, this algorithm identifies sets of contiguous genes acquired by homologous recombination. The locations of the recombination breakpoints are determined using a statistical test that is able to account for the differences in the natural rate of evolution between different genes. The algorithm was tested on a dataset of 75 genomes of Staphylococcus aureus and 50 genomes comprising different streptococcal species, and was able to detect intra-species recombination events in S. aureus and in Streptococcus pneumoniae. Furthermore, we found evidences of an inter-species exchange of genetic material between S. pneumoniae and Streptococcus mitis, a closely related commensal species that colonizes the same ecological niche. The method has been implemented in an R package, Reco, which is freely available from supplementary material, and provides a rapid screening tool to investigate recombination on a genome-wide scale from sequence data.