A combination of cis and trans control can solve the hotspot conversion paradox

There is growing evidence that in a variety of organisms the majority of meiotic recombination events occur at a relatively small fraction of loci, known as recombination hotspots. If hotspot activity results from the DNA sequence at or near the hotspot itself (in cis), these hotspots are expected to be rapidly lost due to biased gene conversion, unless there is strong selection in favor of the hotspot itself. This phenomenon makes it very difficult to maintain existing hotspots and even more difficult for new hotspots to evolve; it has therefore come to be known as the "hotspot conversion paradox." I develop an analytical framework for exploring the evolution of recombination hotspots under the forces of selection, mutation, and conversion. I derive the general conditions under which cis- and trans-controlled hotspots can be maintained, as well as those under which new hotspots controlled by both a cis and a trans locus can invade a population. I show that the conditions for maintenance of and invasion by trans- or cis-plus-trans-controlled hotspots are broader than for those controlled entirely in cis. Finally, I show that a combination of cis and trans control may allow for long-lived polymorphisms in hotspot activity, the patterns of which may explain some recently observed features of recombination hotspots.     

The Red Queen theory of recombination hotspots

Recombination hotspots are small chromosomal regions, where meiotic crossover events happen with high frequency. Recombination is initiated by a double‐strand break (DSB) that requires the intervention of the molecular repair mechanism. The DSB repair mechanism may result in the exchange of homologous chromosomes (crossover) and the conversion of the allelic sequence that breaks into the one that does not break (biased gene conversion). Biased gene conversion results in a transmission advantage for the allele that does not break, thus preventing recombination and rendering recombination hotspots transient. How is it possible that recombination hotspots persist over evolutionary time (maintaining the average chromosomal crossover rate) when they are self‐destructive? This fundamental question is known as the recombination hotspot paradox and has attracted much attention in recent years. Yet, that attention has not translated into a fully satisfactory answer. No existing model adequately explains all aspects of the recombination hotspot paradox. Here, we formulate an intragenomic conflict model resulting in Red Queen dynamics that fully accounts for all empirical observations regarding the molecular mechanisms of recombination hotspots, the nonrandom targeting of the recombination machinery to hotspots and the evolutionary dynamics of hotspot turnover.     

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.     

PRDM9 Drives Evolutionary Erosion of Hotspots in Mus musculus through Haplotype-Specific Initiation of Meiotic Recombination

Meiotic recombination generates new genetic variation and assures the proper segregation of chromosomes in gametes. PRDM9, a zinc finger protein with histone methyltransferase activity, initiates meiotic recombination by binding DNA at recombination hotspots and directing the position of DNA double-strand breaks (DSB). The DSB repair mechanism suggests that hotspots should eventually self-destruct, yet genome-wide recombination levels remain constant, a conundrum known as the hotspot paradox. To test if PRDM9 drives this evolutionary erosion, we measured activity of the Prdm9 ^Cst allele in two Mus musculus subspecies, M.m. castaneus, in which Prdm9^Cst arose, and M.m. domesticus, into which Prdm9^Cst was introduced experimentally. Comparing these two strains, we find that haplotype differences at hotspots lead to qualitative and quantitative changes in PRDM9 binding and activity. Using Mus spretus as an outlier, we found most variants affecting PRDM9^Cst binding arose and were fixed in M.m. castaneus, suppressing hotspot activity. Furthermore, M.m. castaneus×M.m. domesticus F1 hybrids exhibit novel hotspots, with large haplotype biases in both PRDM9 binding and chromatin modification. These novel hotspots represent sites of historic evolutionary erosion that become activated in hybrids due to crosstalk between one parent''s Prdm9 allele and the opposite parent''s chromosome. Together these data support a model where haplotype-specific PRDM9 binding directs biased gene conversion at hotspots, ultimately leading to hotspot erosion.     

Multimer Formation Explains Allelic Suppression of PRDM9 Recombination Hotspots

Genetic recombination during meiosis functions to increase genetic diversity, promotes elimination of deleterious alleles, and helps assure proper segregation of chromatids. Mammalian recombination events are concentrated at specialized sites, termed hotspots, whose locations are determined by PRDM9, a zinc finger DNA-binding histone methyltransferase. Prdm9 is highly polymorphic with most alleles activating their own set of hotspots. In populations exhibiting high frequencies of heterozygosity, questions remain about the influences different alleles have in heterozygous individuals where the two variant forms of PRDM9 typically do not activate equivalent populations of hotspots. We now find that, in addition to activating its own hotspots, the presence of one Prdm9 allele can modify the activity of hotspots activated by the other allele. PRDM9 function is also dosage sensitive; Prdm9 ^+/- heterozygous null mice have reduced numbers and less active hotspots and increased numbers of aberrant germ cells. In mice carrying two Prdm9 alleles, there is allelic competition; the stronger Prdm9 allele can partially or entirely suppress chromatin modification and recombination at hotspots of the weaker allele. In cell cultures, PRDM9 protein variants form functional heteromeric complexes which can bind hotspots sequences. When a heteromeric complex binds at a hotspot of one PRDM9 variant, the other PRDM9 variant, which would otherwise not bind, can still methylate hotspot nucleosomes. We propose that in heterozygous individuals the underlying molecular mechanism of allelic suppression results from formation of PRDM9 heteromers, where the DNA binding activity of one protein variant dominantly directs recombination initiation towards its own hotspots, effectively titrating down recombination by the other protein variant. In natural populations with many heterozygous individuals, allelic competition will influence the recombination landscape.     

The M26 hotspot of Schizosaccharomyces pombe stimulates meiotic ectopic recombination and chromosomal rearrangements.

Homologous recombination is increased during meiosis between DNA sequences at the same chromosomal position (allelic recombination) and at different chromosomal positions (ectopic recombination). Recombination hotspots are important elements in controlling meiotic allelic recombination. We have used artificially dispersed copies of the ade6 gene in Schizosaccharomyces pombe to study hotspot activity in meiotic ectopic recombination. Ectopic recombination was reduced 10-1000-fold relative to allelic recombination, and was similar to the low frequency of ectopic recombination between naturally repeated sequences in S. pombe. The M26 hotspot was active in ectopic recombination in some, but not all, integration sites, with the same pattern of activity and inactivity in ectopic and allelic recombination. Crossing over in ectopic recombination, resulting in chromosomal rearrangements, was associated with 35-60% of recombination events and was stimulated 12-fold by M26. These results suggest overlap in the mechanisms of ectopic and allelic recombination and indicate that hotspots can stimulate chromosomal rearrangements.     

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.     

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.     

Live hot, die young: transmission distortion in recombination hotspots

There is strong evidence that hotspots of meiotic recombination in humans are transient features of the genome. For example, hotspot locations are not shared between human and chimpanzee. Biased gene conversion in favor of alleles that locally disrupt hotspots is a possible explanation of the short lifespan of hotspots. We investigate the implications of such a bias on human hotspots and their evolution. Our results demonstrate that gene conversion bias is a sufficiently strong force to produce the observed lack of sharing of intense hotspots between species, although sharing may be much more common for weaker hotspots. We investigate models of how hotspots arise, and find that only models in which hotspot alleles do not initially experience drive are consistent with observations of rather hot hotspots in the human genome. Mutations acting against drive cannot successfully introduce such hotspots into the population, even if there is direct selection for higher recombination rates, such as to ensure correct segregation during meiosis. We explore the impact of hotspot alleles on patterns of haplotype variation, and show that such alleles mask their presence in population genetic data, making them difficult to detect.     

Recombination hotspots and single-stranded DNA binding proteins couple DNA translocation to DNA unwinding by the AddAB helicase-nuclease

AddAB is a helicase-nuclease that processes double-stranded DNA breaks for repair by homologous recombination. This process is modulated by Chi recombination hotspots: specific DNA sequences that attenuate the nuclease activity of the translocating AddAB complex to promote downstream recombination. Using a combination of kinetic and imaging techniques, we show that AddAB translocation is not coupled to DNA unwinding in the absence of single-stranded DNA binding proteins because nascent single-stranded DNA immediately re-anneals behind the moving enzyme. However, recognition of recombination hotspot sequences during translocation activates unwinding by coupling these activities, thereby ensuring the downstream formation of single-stranded DNA that is required for RecA-mediated recombinational repair. In addition to their implications for the mechanism of double-stranded DNA break repair, these observations may affect our implementation and interpretation of helicase assays and our understanding of helicase mechanisms in general.     

Predicting patterns of long‐term adaptation and extinction with population genetics

Population genetics struggles to model extinction; standard models track the relative rather than absolute fitness of genotypes, while the exceptions describe only the short‐term transition from imminent doom to evolutionary rescue. But extinction can result from failure to adapt not only to catastrophes, but also to a backlog of environmental challenges. We model long‐term adaptation to long series of small challenges, where fitter populations reach higher population sizes. The population's long‐term fitness dynamic is well approximated by a simple stochastic Markov chain model. Long‐term persistence occurs when the rate of adaptation exceeds the rate of environmental deterioration for some genotypes. Long‐term persistence times are consistent with typical fossil species persistence times of several million years. Immediately preceding extinction, fitness declines rapidly, appearing as though a catastrophe disrupted a stably established population, even though gradual evolutionary processes are responsible. New populations go through an establishment phase where, despite being demographically viable, their extinction risk is elevated. Should the population survive long enough, extinction risk later becomes constant over time.     

Molecular mechanisms for environmentally induced and evolutionarily rapid redistribution (plasticity) of meiotic recombination

It has long been known (circa 1917) that environmental conditions, as well as speciation, can affect dramatically the frequency distribution of Spo11/Rec12-dependent meiotic recombination. Here, by analyzing DNA sequence-dependent meiotic recombination hotspots in the fission yeast Schizosaccharomyces pombe, we reveal a molecular basis for these phenomena. The impacts of changing environmental conditions (temperature, nutrients, and osmolarity) on local rates of recombination are mediated directly by DNA site-dependent hotspots (M26, CCAAT, and Oligo-C). This control is exerted through environmental condition-responsive signal transduction networks (involving Atf1, Pcr1, Php2, Php3, Php5, and Rst2). Strikingly, individual hotspots modulate rates of recombination over a very broad dynamic range in response to changing conditions. They can range from being quiescent to being highly proficient at promoting activity of the basal recombination machinery (Spo11/Rec12 complex). Moreover, each different class of hotspot functions as an independently controlled rheostat; a condition that increases the activity of one class can decrease the activity of another class. Together, the independent modulation of recombination rates by each different class of DNA site-dependent hotspots (of which there are many) provides a molecular mechanism for highly dynamic, large-scale changes in the global frequency distribution of meiotic recombination. Because hotspot-activating DNA sites discovered in fission yeast are conserved functionally in other species, this process can also explain the previously enigmatic, Prdm9-independent, evolutionarily rapid changes in hotspot usage between closely related species, subspecies, and isolated populations of the same species.     

A coalescent model of recombination hotspots

Recent experimental findings suggest that the assumption of a homogeneous recombination rate along the human genome is too naive. These findings point to block-structured recombination rates; certain regions (called hotspots) are more prone than other regions to recombination. In this report a coalescent model incorporating hotspot or block-structured recombination is developed and investigated analytically as well as by simulation. Our main results can be summarized as follows: (1) The expected number of recombination events is much lower in a model with pure hotspot recombination than in a model with pure homogeneous recombination, (2) hotspots give rise to large variation in recombination rates along the genome as well as in the number of historical recombination events, and (3) the size of a (nonrecombining) block in the hotspot model is likely to be overestimated grossly when estimated from SNP data. The results are discussed with reference to the current debate about block-structured recombination and, in addition, the results are compared to genome-wide variation in recombination rates. A number of new analytical results about the model are derived.     

Two types of recombination hotspots in bacteriophage T4: one requires DNA damage and a replication origin and the other does not

Recombination hotspots have previously been discovered in bacteriophage T4 by two different approaches, marker rescue recombination from heavily damaged phage genomes and recombination during co-infection by two undamaged phage genomes. The phage replication origin ori(34) is located in a region that has a hotspot in both assays. To determine the relationship between the origin and the two kinds of hotspots, we generated phage carrying point mutations that should inactivate ori(34) but not affect the gene 34 reading frame (within which ori(34) is located). The mutations eliminated the function of the origin, as judged by both autonomous replication of plasmids during T4 infection and two-dimensional gel analysis of phage genomic replication intermediates. As expected from past studies, the ori(34) mutations also eliminated the hotspot for marker rescue recombination from UV-irradiated genomes. However, the origin mutations had no effect on the recombination hotspot that is observed with co-infecting undamaged phage genomes, demonstrating that some DNA sequence other than the origin is responsible for inflated recombination between undamaged genomes. The hotspots for marker rescue recombination may result from a replication fork restart process that acts upon origin-initiated replication forks that become blocked at nearby DNA damage. The two-dimensional gel analysis also revealed phage T4 replication intermediates not previously detected by this method, including origin theta forms.     

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.     

A test of the CoHR motif associated with meiotic double-strand breaks in Saccharomyces cerevisiae

Meiotic recombination is not random along chromosomes; rather, there are preferred regions for initiation called hotspots. Although the general properties of meiotic hotspots are known, the requirements at the DNA sequence level for the determination of hotspot activity are still unclear. The sequence of six known hotspots in Saccharomyces cerevisiae was compared to identify a common homology region (CoHR). They reported that the locations of CoHR sequences correspond to mapped double-strand break (DSB) sites along three chromosomes (I, III, VI). We report here that a deletion of CoHR at HIS2, a hotspot used to identify the motif, has no significant effect on recombination. In the absence of CoHR, DSB formation occurs at a high frequency and at the same sequences as in wild-type strains. In cases where the deletion of sequences containing the CoHR motif has been shown to reduce recombination, we propose that it may be a reflection of the location of the deletion, rather than the loss of CoHR, per se.