Distribution of Recombination Hotspots in the Human Genome – A Comparison of Computer Simulations with Real Data

Recombination is the main cause of genetic diversity. Thus, errors in this process can lead to chromosomal abnormalities. Recombination events are confined to narrow chromosome regions called hotspots in which characteristic DNA motifs are found. Genomic analyses have shown that both recombination hotspots and DNA motifs are distributed unevenly along human chromosomes and are much more frequent in the subtelomeric regions of chromosomes than in their central parts. Clusters of motifs roughly follow the distribution of recombination hotspots whereas single motifs show a negative correlation with the hotspot distribution. To model the phenomena related to recombination, we carried out computer Monte Carlo simulations of genome evolution. Computer simulations generated uneven distribution of hotspots with their domination in the subtelomeric regions of chromosomes. They also revealed that purifying selection eliminating defective alleles is strong enough to cause such hotspot distribution. After sufficiently long time of simulations, the structure of chromosomes reached a dynamic equilibrium, in which number and global distribution of both hotspots and defective alleles remained statistically unchanged, while their precise positions were shifted. This resembles the dynamic structure of human and chimpanzee genomes, where hotspots change their exact locations but the global distributions of recombination events are very similar.     

Patterns of recombination activity on mouse chromosome 11 revealed by high resolution mapping

The success of high resolution genetic mapping of disease predisposition and quantitative trait loci in humans and experimental animals depends on the positions of key crossover events around the gene of interest. In mammals, the majority of recombination occurs at highly delimited 1-2 kb long sites known as recombination hotspots, whose locations and activities are distributed unevenly along the chromosomes and are tightly regulated in a sex specific manner. The factors determining the location of hotspots started to emerge with the finding of PRDM9 as a major hotspot regulator in mammals, however, additional factors modulating hotspot activity and sex specificity are yet to be defined. To address this limitation, we have collected and mapped the locations of 4829 crossover events occurring on mouse chromosome 11 in 5858 meioses of male and female reciprocal F1 hybrids of C57BL/6J and CAST/EiJ mice. This chromosome was chosen for its medium size and high gene density and provided a comparison with our previous analysis of recombination on the longest mouse chromosome 1. Crossovers were mapped to an average resolution of 127 kb, and thirteen hotspots were mapped to <8 kb. Most crossovers occurred in a small number of the most active hotspots. Females had higher recombination rate than males as a consequence of differences in crossover interference and regional variation of sex specific rates along the chromosome. Comparison with chromosome 1 showed that recombination events tend to be positioned in similar fashion along the centromere-telomere axis but independently of the local gene density. It appears that mammalian recombination is regulated on at least three levels, chromosome-wide, regional, and at individual hotspots, and these regulation levels are influenced by sex and genetic background but not by gene content.     

Assembly of eukaryotic algal chromosomes in yeast

Background

Synthetic genomic approaches offer unique opportunities to use powerful yeast and Escherichia coli genetic systems to assemble and modify chromosome-sized molecules before returning the modified DNA to the target host. For example, the entire 1 Mb Mycoplasma mycoides chromosome can be stably maintained and manipulated in yeast before being transplanted back into recipient cells. We have previously demonstrated that cloning in yeast of large (>?~?150 kb), high G?+?C (55%) prokaryotic DNA fragments was improved by addition of yeast replication origins every ~100 kb. Conversely, low G?+?C DNA is stable (up to at least 1.8 Mb) without adding supplemental yeast origins. It has not been previously tested whether addition of yeast replication origins similarly improves the yeast-based cloning of large (> 150 kb) eukaryotic DNA with moderate G?+?C content. The model diatom Phaeodactylum tricornutum has an average G?+?C content of 48% and a 27.4 Mb genome sequence that has been assembled into chromosome-sized scaffolds making it an ideal test case for assembly and maintenance of eukaryotic chromosomes in yeast.

Results

We present a modified chromosome assembly technique in which eukaryotic chromosomes as large as ~500 kb can be assembled from cloned ~100 kb fragments. We used this technique to clone fragments spanning P. tricornutum chromosomes 25 and 26 and to assemble these fragments into single, chromosome-sized molecules. We found that addition of yeast replication origins improved the cloning, assembly, and maintenance of the large chromosomes in yeast. Furthermore, purification of the fragments to be assembled by electroelution greatly increased assembly efficiency.

Conclusions

Entire eukaryotic chromosomes can be successfully cloned, maintained, and manipulated in yeast. These results highlight the improvement in assembly and maintenance afforded by including yeast replication origins in eukaryotic DNA with moderate G?+?C content (48%). They also highlight the increased efficiency of assembly that can be achieved by purifying fragments before assembly.

     

Zinc Finger Binding Motifs Do Not Explain Recombination Rate Variation within or between Species of Drosophila

In humans and mice, the Cys₂His₂ zinc finger protein PRDM9 binds to a DNA sequence motif enriched in hotspots of recombination, possibly modifying nucleosomes, and recruiting recombination machinery to initiate Double Strand Breaks (DSBs). However, since its discovery, some researchers have suggested that the recombinational effect of PRDM9 is lineage or species specific. To test for a conserved role of PRDM9-like proteins across taxa, we use the Drosophila pseudoobscura species group in an attempt to identify recombination associated zinc finger proteins and motifs. We leveraged the conserved amino acid motifs in Cys₂His₂ zinc fingers to predict nucleotide binding motifs for all Cys₂His₂ zinc finger proteins in Drosophila pseudoobscura and identified associations with empirical measures of recombination rate. Additionally, we utilized recombination maps from D. pseudoobscura and D. miranda to explore whether changes in the binding motifs between species can account for changes in the recombination landscape, analogous to the effect observed in PRDM9 among human populations. We identified a handful of potential recombination-associated sequence motifs, but the associations are generally tenuous and their biological relevance remains uncertain. Furthermore, we found no evidence that changes in zinc finger DNA binding explains variation in recombination rate between species. We therefore conclude that there is no protein with a DNA sequence specific human-PRDM9-like function in Drosophila. We suggest these findings could be explained by the existence of a different recombination initiation system in Drosophila.     

Extent of the activity domain and possible roles of FtsK in the Escherichia coli chromosome terminus

Escherichia coli FtsK protein couples cell division and chromosome segregation. It is a component of the septum essential for cell division. It also acts during chromosome dimer resolution by XerCD-specific recombination at the dif site, with two distinct activities: DNA translocation oriented by skewed sequence elements and direct activation of Xer recombination. Dimer resolution requires that the skewed elements polarize in opposite directions 30-50 kb on either side of dif. This constitutes the DIF domain, approximately coincident with the region where replication terminates. The observation that the ftsK1 mutation increases recombination near dif was exploited to determine whether the chromosome region on which FtsK acts is limited to the DIF domain. A monitoring of recombination activity at multiple loci in a 350 kb region to the left of dif revealed (i) zones of differing activities unconnected to dimer resolution and (ii) a constant 10-fold increase of recombination in the 250 kb region adjacent to dif in the ftsK1 mutant. The latter effect allows definition of an FTSK domain whose total size is at least fourfold that of the DIF domain. Additional analyses revealed that FtsK activity responds to polarization in the whole FTSK domain and that displacement of the region where replication terminates preserves differences between recombination zones. Our interpretation is that translocation by FtsK occurs mostly on DNA belonging to a specifically organized domain of the chromosome, when physical links between either dimeric or still intercatenated chromosomes force this DNA to run across the septum at division.     

The recombinational anatomy of a mouse chromosome

Among mammals, genetic recombination occurs at highly delimited sites known as recombination hotspots. They are typically 1-2 kb long and vary as much as a 1,000-fold or more in recombination activity. Although much is known about the molecular details of the recombination process itself, the factors determining the location and relative activity of hotspots are poorly understood. To further our understanding, we have collected and mapped the locations of 5,472 crossover events along mouse Chromosome 1 arising in 6,028 meioses of male and female reciprocal F1 hybrids of C57BL/6J and CAST/EiJ mice. Crossovers were mapped to a minimum resolution of 225 kb, and those in the telomere-proximal 24.7 Mb were further mapped to resolve individual hotspots. Recombination rates were evolutionarily conserved on a regional scale, but not at the local level. There was a clear negative-exponential relationship between the relative activity and abundance of hotspot activity classes, such that a small number of the most active hotspots account for the majority of recombination. Females had 1.2x higher overall recombination than males did, although the sex ratio showed considerable regional variation. Locally, entirely sex-specific hotspots were rare. The initiation of recombination at the most active hotspot was regulated independently on the two parental chromatids, and analysis of reciprocal crosses indicated that parental imprinting has subtle effects on recombination rates. It appears that the regulation of mammalian recombination is a complex, dynamic process involving multiple factors reflecting species, sex, individual variation within species, and the properties of individual hotspots.     

Whole‐genome patterns of linkage disequilibrium across flycatcher populations clarify the causes and consequences of fine‐scale recombination rate variation in birds

Recombination rate is heterogeneous across the genome of various species and so are genetic diversity and differentiation as a consequence of linked selection. However, we still lack a clear picture of the underlying mechanisms for regulating recombination. Here we estimated fine‐scale population recombination rate based on the patterns of linkage disequilibrium across the genomes of multiple populations of two closely related flycatcher species (Ficedula albicollis and F. hypoleuca). This revealed an overall conservation of the recombination landscape between these species at the scale of 200 kb, but we also identified differences in the local rate of recombination despite their recent divergence (<1 million years). Genetic diversity and differentiation were associated with recombination rate in a lineage‐specific manner, indicating differences in the extent of linked selection between species. We detected 400–3,085 recombination hotspots per population. Location of hotspots was conserved between species, but the intensity of hotspot activity varied between species. Recombination hotspots were primarily associated with CpG islands (CGIs), regardless of whether CGIs were at promoter regions or away from genes. Recombination hotspots were also associated with specific transposable elements (TEs), but this association appears indirect due to shared preferences of the transposition machinery and the recombination machinery for accessible open chromatin regions. Our results suggest that CGIs are a major determinant of the localization of recombination hotspots, and we propose that both the distribution of TEs and fine‐scale variation in recombination rate may be associated with the evolution of the epigenetic landscape.     

Mapping candidate hotspots of meiotic recombination in segments of human DNA cloned in the yeast<Emphasis Type="Italic"> Saccharomyces cerevisiae</Emphasis>

The hotspots of meiotic recombination in the human genome can be localized by genetic techniques. The resolution of these techniques is in the range of kilobases and depends on the density of the physical markers identifying allelic variants of the chromosomal loci. We thought it would be interesting to localize these sites with higher resolution. Assuming that some human chromosomal sites conserve their propensity for recombination when cloned in yeast, we localized the hotspots of recombination in several yeast artificial chromosomes (YACs) carrying human DNA. A number of potential recombination hotspots could be identified in the clones studied. Among them there are two classes of sites that are particularly recombination prone also in human meiotic cells: sites associated with CpG islands and sites located in the vicinity of long minisatellite sequences.Communicated by G. P. Georgiev     

How Hot Are Drosophila Hotspots? Examining Recombination Rate Variation and Associations with Nucleotide Diversity,Divergence, and Maternal Age in Drosophila pseudoobscura

Fine scale meiotic recombination maps have uncovered a large amount of variation in crossover rate across the genomes of many species, and such variation in mammalian and yeast genomes is concentrated to <5kb regions of highly elevated recombination rates (10–100x the background rate) called “hotspots.” Drosophila exhibit substantial recombination rate heterogeneity across their genome, but evidence for these highly-localized hotspots is lacking. We assayed recombination across a 40Kb region of Drosophila pseudoobscura chromosome 2, with one 20kb interval assayed every 5Kb and the adjacent 20kb interval bisected into 10kb pieces. We found that recombination events across the 40kb stretch were relatively evenly distributed across each of the 5kb and 10kb intervals, rather than concentrated in a single 5kb region. This, in combination with other recent work, indicates that the recombination landscape of Drosophila may differ from the punctate recombination pattern observed in many mammals and yeast. Additionally, we found no correlation of average pairwise nucleotide diversity and divergence with recombination rate across the 20kb intervals, nor any effect of maternal age in weeks on recombination rate in our sample.     

The Red Queen Model of Recombination Hotspots Evolution in the Light of Archaic and Modern Human Genomes

Recombination is an essential process in eukaryotes, which increases diversity by disrupting genetic linkage between loci and ensures the proper segregation of chromosomes during meiosis. In the human genome, recombination events are clustered in hotspots, whose location is determined by the PRDM9 protein. There is evidence that the location of hotspots evolves rapidly, as a consequence of changes in PRDM9 DNA-binding domain. However, the reasons for these changes and the rate at which they occur are not known. In this study, we investigated the evolution of human hotspot loci and of PRDM9 target motifs, both in modern and archaic human lineages (Denisovan) to quantify the dynamic of hotspot turnover during the recent period of human evolution. We show that present-day human hotspots are young: they have been active only during the last 10% of the time since the divergence from chimpanzee, starting to be operating shortly before the split between Denisovans and modern humans. Surprisingly, however, our analyses indicate that Denisovan recombination hotspots did not overlap with modern human ones, despite sharing similar PRDM9 target motifs. We further show that high-affinity PRDM9 target motifs are subject to a strong self-destructive drive, known as biased gene conversion (BGC), which should lead to the loss of the majority of them in the next 3 MYR. This depletion of PRDM9 genomic targets is expected to decrease fitness, and thereby to favor new PRDM9 alleles binding different motifs. Our refined estimates of the age and life expectancy of human hotspots provide empirical evidence in support of the Red Queen hypothesis of recombination hotspots evolution.     

A multi-step strategy for BAC recombineering of large DNA fragments

Recombineering techniques have been developed to modify bacterial artificial chromosomes (BACs) via bacterial homologous recombination systems, simplifying the molecular manipulations of large DNA constructs. However, precise modifications of a DNA fragment larger than 2-3 kb by recombineering remain a difficult task, due to technical limitations in PCR amplification and purification of large DNA fragments. Here, we describe a new recombineering strategy for the replacement of large DNA fragments using the commonly utilized phage/Red recombination host system. This approach involved the introduction of rare restriction enzyme sites and positive selection markers into the ends of a large DNA fragment, followed by its release from the donor BAC construct and integration into an acceptor BAC. We have successfully employed this method to precisely swap a number of large DNA fragments ranging from 6 to 40 kb between two BAC constructs. Our results demonstrated that this new strategy was highly effective in the manipulations of large genomic DNA fragments and therefore should advance the conventional BAC recombineering technology to the next level.     

Small, repetitive DNAs contribute significantly to the expanded mitochondrial genome of cucumber

Closely related cucurbit species possess eightfold differences in the sizes of their mitochondrial genomes. We cloned mitochondrial DNA (mtDNA) fragments showing strong hybridization signals to cucumber mtDNA and little or no signal to watermelon mtDNA. The cucumber mtDNA clones carried short (30-53 bp), repetitive DNA motifs that were often degenerate, overlapping, and showed no homology to any sequences currently in the databases. On the basis of dot-blot hybridizations, seven repetitive DNA motifs accounted for >13% (194 kb) of the cucumber mitochondrial genome, equaling >50% of the size of the Arabidopsis mitochondrial genome. Sequence analysis of 136 kb of cucumber mtDNA revealed only 11.2% with significant homology to previously characterized mitochondrial sequences, 2.4% to chloroplast DNA, and 15% to the seven repetitive DNA motifs. The remaining 71.4% of the sequence was unique to the cucumber mitochondrial genome. There was <4% sequence colinearity surrounding the watermelon and cucumber atp9 coding regions, and the much smaller watermelon mitochondrial genome possessed no significant amounts of cucumber repetitive DNAs. Our results demonstrate that the expanded cucumber mitochondrial genome is in part due to extensive duplication of short repetitive sequences, possibly by recombination and/or replication slippage.     

Active and Inactive Transplacement of the M26 Recombination Hotspot in Schizosaccharomyces Pombe

The ade6-M26 mutation of the fission yeast Schizosaccharomyces pombe creates a meiotic recombination hotspot that elevates ade6 intragenic recombination ~10-15-fold. A heptanucleotide sequence including the M26 point mutation is required but not sufficient for hotspot activity. We studied the effects of plasmid and chromosomal context on M26 hotspot activity. The M26 hotspot was inactive on a multicopy plasmid containing M26 embedded within 3.0 or 5.9 kb of ade6 DNA. Random S. pombe genomic fragments totaling ~7 Mb did not activate the M26 hotspot on a plasmid. M26 hotspot activity was maintained when 3.0-, 4.4-, and 5.9-kb ade6-M26 DNA fragments, with various amounts of non-S. pombe plasmid DNA, were integrated at the ura4 chromosomal locus, but only in certain configurations relative to the ura4 gene and the cointegrated plasmid DNA. Several integrations created new M26-independent recombination hotspots. In all cases the non-ade6 DNA was located >1 kb from the M26 site, and in some cases >2 kb. Because the chromosomal context effect was transmitted over large distances, and did not appear to be mediated by a single discrete DNA sequence element, we infer that the local chromatin structure has a pronounced effect on M26 hotspot activity.     

Double-strand breaks on artificial chromosomes in yeast

Yeast artificial chromosomes composed primarily of bacteriophage λ DNA exhibit very low levels of meiotic crossing over compared with similarly sized intervals of natural yeast DNA. When these recombinationally quiet chromosomes were augmented with a 12.5 kb insert of sequences from yeast chromosome VIII, genetic studies demonstrated that the artificial chromosomes had acquired recombination properties characteristic of this region of chromosome VIII. On authentic yeast chromosomes, most meiotic recombination events are initiated at sites where the DNA is cleaved to create a double-strand break (DSB). This report describes physical analyses that were carried out to examine the relationship between DSB sites and the recombination behavior of the artificial chromosomes. The results show that DSBs are rare on these artificial chromosomes, except for the 12.5 kb insert. Mapping of the DSB sites shows that their positions correlate with the previously determined positions of DSB sites on chromosome VIII. Deletion of two characterized chromosome VIII DSB sites from the 12.5 kb insert on the artificial chromosome resulted in the loss of the predicted DSB fragments and a reduction in crossing over between artificial chromosomes. Received: 15 May 1998; in revised form: 26 September 1999 / Accepted: 18 November 1999     

Properties of natural double-strand-break sites at a recombination hotspot in Saccharomyces cerevisiae

This study addresses three questions about the properties of recombination hotspots in Saccharomyces cerevisiae: How much DNA is required for double-strand-break (DSB) site recognition? Do naturally occurring DSB sites compete with each other in meiotic recombination? What role does the sequence located at the sites of DSBs play? In S. cerevisiae, the HIS2 meiotic recombination hotspot displays a high level of gene conversion, a 3''-to-5'' conversion gradient, and two DSB sites located approximately 550 bp apart. Previous studies of hotspots, including HIS2, suggest that global chromosome structure plays a significant role in recombination activity, raising the question of how much DNA is sufficient for hotspot activity. We find that 11.5 kbp of the HIS2 region is sufficient to partially restore gene conversion and both DSBs when moved to another yeast chromosome. Using a variety of different constructs, studies of hotspots have indicated that DSB sites compete with one another for DSB formation. The two naturally occurring DSBs at HIS2 afforded us the opportunity to examine whether or not competition occurs between these native DSB sites. Small deletions of DNA at each DSB site affect only that site; analyses of these deletions show no competition occurring in cis or in trans, indicating that DSB formation at each site at HIS2 is independent. These small deletions significantly affect the frequency of DSB formation at the sites, indicating that the DNA sequence located at a DSB site can play an important role in recombination initiation.     

DNA sequence-mediated, evolutionarily rapid redistribution of meiotic recombination hotspots

Hotspots regulate the position and frequency of Spo11 (Rec12)-initiated meiotic recombination, but paradoxically they are suicidal and are somehow resurrected elsewhere in the genome. After the DNA sequence-dependent activation of hotspots was discovered in fission yeast, nearly two decades elapsed before the key realizations that (A) DNA site-dependent regulation is broadly conserved and (B) individual eukaryotes have multiple different DNA sequence motifs that activate hotspots. From our perspective, such findings provide a conceptually straightforward solution to the hotspot paradox and can explain other, seemingly complex features of meiotic recombination. We describe how a small number of single-base-pair substitutions can generate hotspots de novo and dramatically alter their distribution in the genome. This model also shows how equilibrium rate kinetics could maintain the presence of hotspots over evolutionary timescales, without strong selective pressures invoked previously, and explains why hotspots localize preferentially to intergenic regions and introns. The model is robust enough to account for all hotspots of humans and chimpanzees repositioned since their divergence from the latest common ancestor.     

Structural features of 3C6/C7 interband chromatin organization in Drosophila melanogaster polytene chromosomes

Mechanisms of chromatin decompaction in interbands of Drosophila polytene chromosomes have been studied. Using the example of interband 3C6/C7 of the X chromosome, we investigated the ability of different DNA segments to form an interband in a new genetic environment. We applied site-specific FLP recombination between two transposons with FRT-sites that allows introducing the DNA fragments from the interband 3C6/C7 into pICon(dv) transposon located in cytologically well-characterized 84F region of chromosome 3 followed by electron microscopic analysis of changes in the region caused by insertion of the DNA fragments into the transposon. It was shown that the insertion of a 276-bp DNA fragment from the 3C6/C7 region into the pICon(dv) transposon leads to the formation of a new interband between two thin bands represented by the transposon material. This DNA fragment is the known minimal sequence that is necessary and sufficient for interband generation. In addition, the sequence containing three copies repeated in tandem of 0.9 kb DNA from the interband 3C6/C7, including the 276-bp fragment, were integrated in the transposon. The presence of introduced DNA fragments did not change the morphology of the resulting interband. It was shown that the sites of DNase I hypersensitivity were saved in transposon sequences introduced into a new genetic environment. The data obtained allow analysis to be started of specific factors (proteins, DNA motifs, etc.) that determine the formation of decompacted chromatin in a certain interband region and chromomeric organization of interphase chromosomes in Drosophila as a whole.     

Organelle DNA accumulation in the recently evolved papaya sex chromosomes

Sex chromosomes are a pair of specialized chromosomes containing a sex determination region that is suppressed for recombination. Without recombination, Y chromosomes are thought to accumulate repetitive DNA sequences which contribute to their degeneration. A pair of primitive sex chromosomes controls sex type in papaya with male and hermaphrodite determined by the slightly different male-specific region of the Y (MSY) and hermaphrodite-specific region of Y^h (HSY) chromosomes, respectively. Here, we show that the papaya HSY and MSY in the absence of recombination have accumulated nearly 12 times the amount of chloroplast-derived DNA than the corresponding region of the X chromosome and 4 times the papaya genome-wide average. Furthermore, a chloroplast genome fragment containing the rsp15 gene has been amplified 23 times in the HSY, evidence of retrotransposon-mediated duplication. Surprisingly, mitochondria-derived sequences are less abundant in the X and HSY compared to the whole genome. Shared organelle integrations are sparse between X and HSY, with only 11 % of chloroplast and 12 % of mitochondria fragments conserved, respectively, suggesting that the accelerated accumulation of organelle DNA occurred after the HSY was suppressed for recombination. Most of the organelle-derived sequences have divergence times of <7 MYA, reinforcing this notion. The accumulated chloroplast DNA is evidence of the slow degeneration of the HSY.