What are the genomic drivers of the rapid evolution of PRDM9?

The role of recombination dynamics in shaping signatures of direct and indirect selection across the Ficedula flycatcher genome†

Recombination is a central evolutionary process that reshuffles combinations of alleles along chromosomes, and consequently is expected to influence the efficacy of direct selection via Hill-Robertson interference. Additionally, the indirect effects of selection on neutral genetic diversity are expected to show a negative relationship with recombination rate, as background selection and genetic hitchhiking are stronger when recombination rate is low. However, owing to the limited availability of recombination rate estimates across divergent species, the impact of evolutionary changes in recombination rate on genomic signatures of selection remains largely unexplored. To address this question, we estimate recombination rate in two Ficedula flycatcher species, the taiga flycatcher (Ficedula albicilla) and collared flycatcher (Ficedula albicollis). We show that recombination rate is strongly correlated with signatures of indirect selection, and that evolutionary changes in recombination rate between species have observable impacts on this relationship. Conversely, signatures of direct selection on coding sequences show little to no relationship with recombination rate, even when restricted to genes where recombination rate is conserved between species. Thus, using measures of indirect and direct selection that bridge micro- and macro-evolutionary timescales, we demonstrate that the role of recombination rate and its dynamics varies for different signatures of selection.

Meiosis: Dances Between Homologs

The raison d'être of meiosis is shuffling of genetic information via Mendelian segregation and, within individual chromosomes, by DNA crossing-over. These outcomes are enabled by a complex cellular program in which interactions between homologous chromosomes play a central role. We first provide a background regarding the basic principles of this program. We then summarize the current understanding of the DNA events of recombination and of three processes that involve whole chromosomes: homolog pairing, crossover interference, and chiasma maturation. All of these processes are implemented by direct physical interaction of recombination complexes with underlying chromosome structures. Finally, we present convergent lines of evidence that the meiotic program may have evolved by coupling of this interaction to late-stage mitotic chromosome morphogenesis.

The Recombination Hotspot Paradox: Co-evolution between PRDM9 and its target sites

Recombination often concentrates in small regions called recombination hotspots where recombination is much higher than the genome's average. In many vertebrates, including humans, gene PRDM9 specifies which DNA motifs will be the target for breaks that initiate recombination, ultimately determining the location of recombination hotspots. Because the sequence that breaks (allowing recombination) is converted into the sequence that does not break (preventing recombination), the latter sequence is over-transmitted to future generations and recombination hotspots are self-destructive. Given their self-destructive nature, recombination hotspots should eventually become extinct in genomes where they are found. While empirical evidence shows that individual hotspots do become inactive over time (die), hotspots are abundant in many vertebrates: a contradiction called the Recombination Hotspot Paradox. What saves recombination hotspots from their foretold extinction? Here we formulate a co-evolutionary model of the interaction among sequence-specific gene conversion, fertility selection, and recurrent mutation. We find that allelic frequencies oscillate leading to stable limit cycles. From a biological perspective this means that when fertility selection is weaker than gene conversion, it cannot stop individual hotspots from dying but can save them from extinction by driving their re-activation (resuscitation). In our model, mutation balances death and resuscitation of hotspots, thus maintaining their number over evolutionary time. Interestingly, we find that multiple alleles result in oscillations that are chaotic and multiple targets in oscillations that are asynchronous between targets thus helping to maintain the average genomic recombination probability constant. Furthermore, we find that the level of expression of PRDM9 should control for the fraction of targets that are hotspots and the overall temperature of the genome. Therefore, our co-evolutionary model improves our understanding of how hotspots may be replaced, thus contributing to solve the Recombination Hotspot Paradox. From a more applied perspective our work provides testable predictions regarding the relation between mutation probability and fertility selection with life expectancy of hotspots.

Incompatibility and Interchangeability in Molecular Evolution

There is remarkable variation in the rate at which genetic incompatibilities in molecular interactions accumulate. In some cases, minor changes-even single-nucleotide substitutions-create major incompatibilities when hybridization forces new variants to function in a novel genetic background from an isolated population. In other cases, genes or even entire functional pathways can be horizontally transferred between anciently divergent evolutionary lineages that span the tree of life with little evidence of incompatibilities. In this review, we explore whether there are general principles that can explain why certain genes are prone to incompatibilities while others maintain interchangeability. We summarize evidence pointing to four genetic features that may contribute to greater resistance to functional replacement: (1) function in multisubunit enzyme complexes and protein-protein interactions, (2) sensitivity to changes in gene dosage, (3) rapid rate of sequence evolution, and (4) overall importance to cell viability, which creates sensitivity to small perturbations in molecular function. We discuss the relative levels of support for these different hypotheses and lay out future directions that may help explain the striking contrasts in patterns of incompatibility and interchangeability throughout the history of molecular evolution.

Fertility of hybrids of dromedary and Bactrian camels: A possible role of conserved architecture of zinc finger domain of recombination regulator PRDM9

Recombination regulator, PRDM9, has been regarded as the most rapidly evolving gene in the genomes of many metazoans, in addition to being acknowledged as the sole speciation gene in vertebrates. It has become the focus of many scientific investigations because of exceptional numerical and sequence variability in its zinc finger (ZF) domain within and across species that contributes to reproductive isolation between species. This study is the maiden attempt to explore the architecture of PRDM9 ZF domain in two Camelid species (Camelus dromedarius and Camelus bactrianus). Sequence analysis revealed highly conserved domain architecture with presence of 3 and 4 ZFs in dromedary and Bactrian camels, respectively. Typical evolutionary features of PRDM9 ZF domain i.e. concerted evolution and positive selection were invariably absent in both the one-humped dromedary and the two-humped Bactrian camels. Fertility of hybrids of dromedary and Bactrian camels, despite being taxonomically distinct species can be attributed to the lack of sequence variability in PRDM9 in these species. Phylogenetic analysis underpinned clear demarcation of camels from other livestock species. The results of the present study defy what has been learnt so far about PRDM9 and add to the enigma surrounding the most intriguing gene in the genome.

PRDM9 losses in vertebrates are coupled to those of paralogs ZCWPW1 and ZCWPW2

We take a phylogenetic approach to search for molecular partners of PRDM9, a key meiotic recombination gene, by leveraging the fact that the complete PRDM9 gene has been lost at least 13 times independently in vertebrates. We identify two genes, ZCWPW1 and its paralog ZCWPW2, whose presence or absence across vertebrates is coupled to that of PRDM9. ZCWPW1 was recently shown to be recruited to sites of PRDM9 binding and to aid in the repair of double strand breaks. ZCWPW2 is likely recruited to sites of PRDM9 binding as well; its tight coevolution with PRDM9 across vertebrates suggests that it too plays an important role in mammals and beyond, either in double strand break formation or repair.

In most mammals and likely throughout vertebrates, the gene PRDM9 specifies the locations of meiotic double strand breaks; in mice and humans at least, it also aids in their repair. For both roles, many of the molecular partners remain unknown. Here, we take a phylogenetic approach to identify genes that may be interacting with PRDM9 by leveraging the fact that PRDM9 arose before the origin of vertebrates but was lost many times, either partially or entirely—and with it, its role in recombination. As a first step, we characterize PRDM9 domain composition across 446 vertebrate species, inferring at least 13 independent losses. We then use the interdigitation of PRDM9 orthologs across vertebrates to test whether it coevolved with any of 241 candidate genes coexpressed with PRDM9 in mice or associated with recombination phenotypes in mammals. Accounting for the phylogenetic relationship among a subsample of 189 species, we find two genes whose presence and absence is unexpectedly coincident with that of PRDM9: ZCWPW1, which was recently shown to facilitate double strand break repair, and its paralog ZCWPW2, as well as, more tentatively, TEX15 and FBXO47. ZCWPW2 is expected to be recruited to sites of PRDM9 binding; its tight coevolution with PRDM9 across vertebrates suggests that it is a key interactor within mammals and beyond, with a role either in recruiting the recombination machinery or in double strand break repair.

Snake Recombination Landscapes Are Concentrated in Functional Regions despite PRDM9

Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions-a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions-a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages.

The Role of Gene Conversion between Transposable Elements in Rewiring Regulatory Networks

Nature has found many ways to utilize transposable elements (TEs) throughout evolution. Many molecular and cellular processes depend on DNA-binding proteins recognizing hundreds or thousands of similar DNA motifs dispersed throughout the genome that are often provided by TEs. It has been suggested that TEs play an important role in the evolution of such systems, in particular, the rewiring of gene regulatory networks. One mechanism that can further enhance the rewiring of regulatory networks is nonallelic gene conversion between copies of TEs. Here, we will first review evidence for nonallelic gene conversion in TEs. Then, we will illustrate the benefits nonallelic gene conversion provides in rewiring regulatory networks. For instance, nonallelic gene conversion between TE copies offers an alternative mechanism to spread beneficial mutations that improve the network, it allows multiple mutations to be combined and transferred together, and it allows natural selection to work efficiently in spreading beneficial mutations and removing disadvantageous mutations. Future studies examining the role of nonallelic gene conversion in the evolution of TEs should help us to better understand how TEs have contributed to evolution.

The expression of histone methyltransferases and distribution of selected histone methylations in testes of yak and cattle-yak hybrid

Interspecies hybridization exists widely in nature and plays an important role in animal evolution and adaptation. It is commonly recognized that male offspring of interspecies hybrid are often sterile, which presents a crucial way of reproductive isolation. Currently, the mechanisms underlying interspecies hybrid male sterility are not well understood. Cattle-yak, progeny of yak (Bos grunniens) and cattle (Bos taurus) cross, is a unique animal model for investigating hybrid male sterility. Because histone modifications are vital for spermatogenesis, herein, we examined expressions of histone methyltransferases (HMTs) and distributions of histone methylations in the yak and cattle-yak testis. Histological examination of seminiferous tubules revealed that gonocytes and spermatocytes were established normally, however, spermatogenesis was arrested at the meiosis phase began at 10 months after birth in the hybrids. SUV420H1 was the only HMT examined showing a significant enrichment in cattle-yak testes at 3 months. Relative expressions of MLL5, SETDB1 and SUV420H1 were increased while SETDB2 and EZH2 were decreased in cattle-yak testes at 10 months. Relative concentrations of MLL5 and SUV420H1 were again increased while EHMT2 and PRDM9 expressions were decreased at 24 months. Immunofluorescent detection of selected histone methylations in cross-sections of testicular tissues or meiotic chromosomes demonstrated that depletion of H3K4me3 and significant enrichment of H3K27me3 and H4K20me3 were observed in Sertoli cells of cattle-yak. Levels and localizations of H3K4me3, H3K9me1, H3K9me3 and H4K20me3 were strikingly different in meiotic chromosomes of cattle-yak spermatocytes. These results highlighted the potential roles of histone methylations in spermatogenic failure and hybrid male sterility.

A rapidly evolved domain, the SCML2 DNA-binding repeats, contributes to chromatin binding of mouse SCML2†

Genes involved in sexual reproduction diverge rapidly as a result of reproductive fitness. Here, we identify a novel protein domain in the germline-specific Polycomb protein SCML2 that is required for the establishment of unique gene expression programs after the mitosis-to-meiosis transition in spermatogenesis. We term this novel domain, which is comprised of rapidly evolved, DNA-binding repeat units of 28 amino acids, the SCML2 DNA-binding (SDB) repeats. These repeats are acquired in a specific subgroup of the rodent lineage, having been subjected to positive selection in the course of evolution. Mouse SCML2 has two DNA-binding domains: one is the SDB repeats and the other is an RNA-binding region, which is conserved in human SCML2. For the recruitment of SCML2 to target loci, the SDB repeats cooperate with the other functional domains of SCML2 to bind chromatin. The cooperative action of these domains enables SCML2 to sense DNA hypomethylation in an in vivo chromatin environment, thereby enabling SCML2 to bind to hypomethylated chromatin. We propose that the rapid evolution of SCML2 is due to reproductive adaptation, which has promoted species-specific gene expression programs in spermatogenesis.

Where to Cross Over? Defining Crossover Sites in Plants

It is believed that recombination in meiosis serves to reshuffle genetic material from both parents to increase genetic variation in the progeny. At the same time, the number of crossovers is usually kept at a very low level. As a consequence, many organisms need to make the best possible use from the one or two crossovers that occur per chromosome in meiosis. From this perspective, the decision of where to allocate rare crossover events becomes an important issue, especially in self-pollinating plant species, which experience limited variation due to inbreeding. However, the freedom in crossover allocation is significantly limited by other, genetic and non-genetic factors, including chromatin structure. Here we summarize recent progress in our understanding of those processes with a special emphasis on plant genomes. First, we focus on factors which influence the distribution of recombination initiation sites and discuss their effects at both, the single hotspot level and at the chromosome scale. We also briefly explain the aspects of hotspot evolution and their regulation. Next, we analyze how recombination initiation sites translate into the development of crossovers and their location. Moreover, we provide an overview of the sequence polymorphism impact on crossover formation and chromosomal distribution.

Aberrant PRDM9 expression impacts the pan-cancer genomic landscape

The binding of PRDM9 to chromatin is a key step in the induction of DNA double-strand breaks associated with meiotic recombination hotspots; it is normally expressed solely in germ cells. We interrogated 1879 cancer samples in 39 different cancer types and found that PRDM9 is unexpectedly expressed in 20% of these tumors even after stringent gene homology correction. The expression levels of PRDM9 in tumors are significantly higher than those found in healthy neighboring tissues and in healthy nongerm tissue databases. Recurrently mutated regions located within 5 Mb of the PRDM9 loci, as well as differentially expressed genes in meiotic pathways, correlate with PRDM9 expression. In samples with aberrant PRDM9 expression, structural variant breakpoints frequently neighbor the DNA motif recognized by PRDM9, and there is an enrichment of structural variants at sites of known meiotic PRDM9 activity. This study is the first to provide evidence of an association between aberrant expression of the meiosis-specific gene PRDM9 with genomic instability in cancer.

PRDM9, a driver of the genetic map

During meiosis, maternal and paternal chromosomes undergo exchanges by homologous recombination. This is essential for fertility and contributes to genome evolution. In many eukaryotes, sites of meiotic recombination, also called hotspots, are regions of accessible chromatin, but in many vertebrates, their location follows a distinct pattern and is specified by PR domain-containing protein 9 (PRDM9). The specification of meiotic recombination hotspots is achieved by the different activities of PRDM9: DNA binding, histone methyltransferase, and interaction with other proteins. Remarkably, PRDM9 activity leads to the erosion of its own binding sites and the rapid evolution of its DNA-binding domain. PRDM9 may also contribute to reproductive isolation, as it is involved in hybrid sterility potentially due to a reduction of its activity in specific heterozygous contexts.

Connecting theory and data to understand recombination rate evolution

Meiotic recombination is necessary for successful gametogenesis in most sexually reproducing organisms and is a fundamental genomic parameter, influencing the efficacy of selection and the fate of new mutations. The molecular and evolutionary functions of recombination should impose strong selective constraints on the range of recombination rates. Yet, variation in recombination rate is observed on a variety of genomic and evolutionary scales. In the past decade, empirical studies have described variation in recombination rate within genomes, between individuals, between sexes, between populations and between species. At the same time, theoretical work has provided an increasingly detailed picture of the evolutionary advantages to recombination. Perhaps surprisingly, the causes of natural variation in recombination rate remain poorly understood. We argue that empirical and theoretical approaches to understand the evolution of recombination have proceeded largely independently of each other. Most models that address the evolution of recombination rate were created to explain the evolutionary advantage of recombination rather than quantitative differences in rate among individuals. Conversely, most empirical studies aim to describe variation in recombination rate, rather than to test evolutionary hypotheses. In this Perspective, we argue that efforts to integrate the rich bodies of empirical and theoretical work on recombination rate are crucial to moving this field forward. We provide new directions for the development of theory and the production of data that will jointly close this gap.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

The Red Queen model of recombination hot-spot evolution: a theoretical investigation

In humans and many other species, recombination events cluster in narrow and short-lived hot spots distributed across the genome, whose location is determined by the Zn-finger protein PRDM9. To explain these fast evolutionary dynamics, an intra-genomic Red Queen model has been proposed, based on the interplay between two antagonistic forces: biased gene conversion, mediated by double-strand breaks, resulting in hot-spot extinction, followed by positive selection favouring new PRDM9 alleles recognizing new sequence motifs. Thus far, however, this Red Queen model has not been formalized as a quantitative population-genetic model, fully accounting for the intricate interplay between biased gene conversion, mutation, selection, demography and genetic diversity at the PRDM9 locus. Here, we explore the population genetics of the Red Queen model of recombination. A Wright–Fisher simulator was implemented, allowing exploration of the behaviour of the model (mean equilibrium recombination rate, diversity at the PRDM9 locus or turnover rate) as a function of the parameters (effective population size, mutation and erosion rates). In a second step, analytical results based on self-consistent mean-field approximations were derived, reproducing the scaling relations observed in the simulations. Empirical fit of the model to current data from the mouse suggests both a high mutation rate at PRDM9 and strong biased gene conversion on its targets.

This article is part of the themed issue ‘Evolutionary causes and consequences of recombination rate variation in sexual organisms’.

Interrogating the Functions of PRDM9 Domains in Meiosis

Homologous recombination is required for proper segregation of homologous chromosomes during meiosis. It occurs predominantly at recombination hotspots that are defined by the DNA binding specificity of the PRDM9 protein. PRDM9 contains three conserved domains typically involved in regulation of transcription; yet, the role of PRDM9 in gene expression control is not clear. Here, we analyze the germline transcriptome of Prdm9^−/− male mice in comparison to Prdm9^+/+ males and find no apparent differences in the mRNA and miRNA profiles. We further explore the role of PRDM9 in meiosis by analyzing the effect of the KRAB, SSXRD, and post-SET zinc finger deletions in a cell culture expression system and the KRAB domain deletion in mice. We found that although the post-SET zinc finger and the KRAB domains are not essential for the methyltransferase activity of PRDM9 in cell culture, the KRAB domain mutant mice show only residual PRDM9 methyltransferase activity and undergo meiotic arrest. In aggregate, our data indicate that domains typically involved in regulation of gene expression do not serve that role in PRDM9, but are likely involved in setting the proper chromatin environment for initiation and completion of homologous recombination.

Construction of PRDM9 allele-specific recombination maps in cattle using large-scale pedigree analysis and genome-wide single sperm genomics

PRDM9 contributes to hybrid sterility and species evolution. However, its role is to be confirmed in cattle, a major domesticated livestock species. We previously found an association near PRDM9 with cattle recombination features, but the causative variants are still unknown. Using millions of genotyped cattle with pedigree information, we characterized five PRDM9 alleles and generated allele-specific recombination maps. By examining allele-specific recombination patterns, we observed the impact of PRDM9 on global distribution of recombination, especially in the two ends of chromosomes. We also showed strong associations between recombination hotspot regions and functional mutations within PRDM9 zinc finger domain. More importantly, we found one allele of PRDM9 to be very different from others in both protein composition and recombination landscape, indicating the causative role of this allele on the association between PRDM9 and cattle recombination. When comparing recombination maps from sperm and pedigree data, we observed similar genome-wide recombination patterns, validating the quality of pedigree-based results. Collectively, these evidence supported PRDM9 alleles as causal variants for the reported association with cattle recombination. Our study comprehensively surveyed the bovine PRDM9 alleles, generated allele-specific recombination maps, and expanded our understanding of the role of PRDM9 on genome distribution of recombination.

Birth and death of a protein

The ways in which recombination sites are determined during meiosis are becoming clearer following a phylogenomic analysis for 225 different species.

The Red Queen and King in finite populations

In antagonistic symbioses, such as host-parasite interactions, one population's success is the other's loss. In mutualistic symbioses, such as division of labor, both parties can gain, but they might have different preferences over the possible mutualistic arrangements. The rates of evolution of the two populations in a symbiosis are important determinants of which population will be more successful: Faster evolution is thought to be favored in antagonistic symbioses (the "Red Queen effect"), but disfavored in certain mutualistic symbioses (the "Red King effect"). However, it remains unclear which biological parameters drive these effects. Here, we analyze the effects of the various determinants of evolutionary rate: generation time, mutation rate, population size, and the intensity of natural selection. Our main results hold for the case where mutation is infrequent. Slower evolution causes a long-term advantage in an important class of mutualistic interactions. Surprisingly, less intense selection is the strongest driver of this Red King effect, whereas relative mutation rates and generation times have little effect. In antagonistic interactions, faster evolution by any means is beneficial. Our results provide insight into the demographic evolution of symbionts.

Ruminant-specific multiple duplication events of PRDM9 before speciation

Background

Understanding the genetic and evolutionary mechanisms of speciation genes in sexually reproducing organisms would provide important insights into mammalian reproduction and fitness. PRDM9, a widely known speciation gene, has recently gained attention for its important role in meiotic recombination and hybrid incompatibility. Despite the fact that PRDM9 is a key regulator of recombination and plays a dominant role in hybrid incompatibility, little is known about the underlying genetic and evolutionary mechanisms that generated multiple copies of PRDM9 in many metazoan lineages.

Results

The present study reports (1) evidence of ruminant-specific multiple gene duplication events, which likely have had occurred after the ancestral ruminant population diverged from its most recent common ancestor and before the ruminant speciation events, (2) presence of three copies of PRDM9, one copy (lineages I) in chromosome 1 (chr1) and two copies (lineages II & III) in chromosome X (chrX), thus indicating the possibility of ancient inter- and intra-chromosomal unequal crossing over and gene conversion events, (3) while lineages I and II are characterized by the presence of variable tandemly repeated C2H2 zinc finger (ZF) arrays, lineage III lost these arrays, and (4) C2H2 ZFs of lineages I and II, particularly the amino acid residues located at positions −1, 3, and 6 have evolved under strong positive selection.

Conclusions

Our results demonstrated two gene duplication events of PRDM9 in ruminants: an inter-chromosomal duplication that occurred between chr1 and chrX, and an intra-chromosomal X-linked duplication, which resulted in two additional copies of PRDM9 in ruminants. The observation of such duplication between chrX and chr1 is rare and may possibly have happened due to unequal crossing-over millions of years ago when sex chromosomes were independently derived from a pair of ancestral autosomes. Two copies (lineages I & II) are characterized by the presence of variable sized tandem-repeated C2H2 ZFs and evolved under strong positive selection and concerted evolution, supporting the notion of well-established Red Queen hypothesis. Collectively, gene duplication, concerted evolution, and positive selection are the likely driving forces for the expansion of ruminant PRDM9 sub-family.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0892-4) contains supplementary material, which is available to authorized users.

Prophase I: Preparing Chromosomes for Segregation in the Developing Oocyte

Formation of an oocyte involves a specialized cell division termed meiosis. In meiotic prophase I (the initial stage of meiosis), chromosomes undergo elaborate events to ensure the proper segregation of their chromosomes into gametes. These events include processes leading to the formation of a crossover that, along with sister chromatid cohesion, forms the physical link between homologous chromosomes. Crossovers are formed as an outcome of recombination. This process initiates with programmed double-strand breaks that are repaired through the use of homologous chromosomes as a repair template. The accurate repair to form crossovers takes place in the context of the synaptonemal complex, a protein complex that links homologous chromosomes in meiotic prophase I. To allow proper execution of meiotic prophase I events, signaling processes connect different steps in recombination and synapsis. The events occurring in meiotic prophase I are a prerequisite for proper chromosome segregation in the meiotic divisions. When these processes go awry, chromosomes missegregate. These meiotic errors are thought to increase with aging and may contribute to the increase in aneuploidy observed in advanced maternal age female oocytes.

The pioneering role of PRDM9 indel mutations in tarsier evolution

PRDM9 is currently the sole speciation gene found in vertebrates causing hybrid sterility probably due to incompatible alleles. Its role in defining the double strand break loci during the meiotic prophase I is crucial for proper chromosome segregation. Therefore, the rapid turnover of the loci determining zinc finger array seems to be causative for incompatibilities. We here investigated the zinc finger domain-containing exon of PRDM9 in 23 tarsiers. Tarsiers, the most basal extant haplorhine primates, exhibit two frameshifting indels at the 5'-end of the array. The first mutation event interrupts the reading frame and function while the second compensates both. The fixation of this allele variant in tarsiers led to hypothesize that de- and reactivation of the zinc finger domain drove the speciation in early haplorhine or tarsiiform primates. Moreover, the high allelic diversity within Tarsius points to multiple effects of genetic drift reflecting their phylogeographic history since the Miocene.

Evidence of positive selection and concerted evolution in the rapidly evolving PRDM9 zinc finger domain in goats and sheep

Meiotic recombination contributes to augmentation of genetic diversity, exclusion of deleterious alleles and proper segregation of chromatids. PRDM9 has been identified as the gene responsible for specifying the location of recombination hotspots during meiosis and is also the only known vertebrate gene associated with reproductive isolation between species. PRDM9 encodes a protein with a highly variable zinc finger (ZF) domain that varies between as well as within species. In the present study, the ZF domain of PRDM9 on chromosome 1 was characterized for the first time in 15 goat breeds and 25 sheep breeds of India. A remarkable variation in the number and sequence of ZF domains was observed. The number of ZF repeats in the ZF array varied from eight to 12 yielding five homozygous and 10 heterozygous genotypes. The number of different ZF domains was 84 and 52 producing 36 and 26 unique alleles in goats and sheep respectively. The posterior mean of dN/dS or omega values were calculated using the codeml tool of pamlx to identify amino acids that are evolving positively in goats and sheep, as positions -1, +3 and +6 in the ZF domain have been reported to experience strong positive selection across different lineages. Our study identified sites -5, -1, +3, +4 and +6 to be experiencing positive selection. Small ruminant zinc fingers were also found to be evolving under concerted evolution. Our results demonstrate the existence of a vast diversity of PRDM9 in goats and sheep, which is in concert with reports in many metazoans.

The Meiotic Recombination Activator PRDM9 Trimethylates Both H3K36 and H3K4 at Recombination Hotspots In Vivo

In many mammals, including humans and mice, the zinc finger histone methyltransferase PRDM9 performs the first step in meiotic recombination by specifying the locations of hotspots, the sites of genetic recombination. PRDM9 binds to DNA at hotspots through its zinc finger domain and activates recombination by trimethylating histone H3K4 on adjacent nucleosomes through its PR/SET domain. Recently, the isolated PR/SET domain of PRDM9 was shown capable of also trimethylating H3K36 in vitro, raising the question of whether this reaction occurs in vivo during meiosis, and if so, what its function might be. Here, we show that full-length PRDM9 does trimethylate H3K36 in vivo in mouse spermatocytes. Levels of H3K4me3 and H3K36me3 are highly correlated at hotspots, but mutually exclusive elsewhere. In vitro, we find that although PRDM9 trimethylates H3K36 much more slowly than it does H3K4, PRDM9 is capable of placing both marks on the same histone molecules. In accord with these results, we also show that PRDM9 can trimethylate both K4 and K36 on the same nucleosomes in vivo, but the ratio of K4me3/K36me3 is much higher for the pair of nucleosomes adjacent to the PRDM9 binding site compared to the next pair further away. Importantly, H3K4me3/H3K36me3-double-positive nucleosomes occur only in regions of recombination: hotspots and the pseudoautosomal (PAR) region of the sex chromosomes. These double-positive nucleosomes are dramatically reduced when PRDM9 is absent, showing that this signature is PRDM9-dependent at hotspots; the residual double-positive nucleosomes most likely come from the PRDM9-independent PAR. These results, together with the fact that PRDM9 is the only known mammalian histone methyltransferase with both H3K4 and H3K36 trimethylation activity, suggest that trimethylation of H3K36 plays an important role in the recombination process. Given the known requirement of H3K36me3 for double strand break repair by homologous recombination in somatic cells, we suggest that it may play the same role in meiosis.

Genetic recombination is the meiotic process by which novel combinations of alleles are passed on to the next generation. This process accelerates evolution by creating genetic diversity, and is also essential for successful meiosis. In mammals, the enzyme PRDM9 initiates recombination and determines the subset of sites within the genome—called recombination hotspots—that recombine. PRDM9 does this by binding DNA at hotspots and placing epigenetic marks. Previously, PRDM9 was only known to place the H3K4me3 mark at hotspots in living cells. Here, we show that PRDM9 places the H3K36me3 mark at hotspots, and that H3K4me3 and H3K36me3 coincide only at regions of recombination in germ cells. We prove that this coincidence is driven almost entirely by PRDM9; there is dramatically less coincidence between H3K4me3 and H3K36me3 when PRDM9 is absent. These results reveal a new enzymatic function for PRDM9, and a new epigenetic signature at hotspots that may restrict recombination to these sites. Since aberrant recombination can cause aneuploidy resulting in fetal loss, and abnormal genome rearrangements that underlie many congenital syndromes and some cancers, a thorough understanding of this fundamental process has potentially far-reaching implications.

Zinc Finger Domain of the PRDM9 Gene on Chromosome 1 Exhibits High Diversity in Ruminants but Its Paralog PRDM7 Contains Multiple Disruptive Mutations

PRDM9 is the sole hybrid sterility gene identified so far in vertebrates. PRDM9 gene encodes a protein with an immensely variable zinc-finger (ZF) domain that determines the site of meiotic recombination hotspots genome-wide. In this study, the terminal ZF domain of PRDM9 on bovine chromosome 1 and its paralog on chromosome 22 were characterized in 225 samples from five ruminant species (cattle, yak, mithun, sheep and goat). We found extraordinary variation in the number of PRDM9 zinc fingers (6 to 12). We sequenced PRDM9 ZF encoding region from 15 individuals (carrying the same ZF number in both copies) and found 43 different ZF domain sequences. Ruminant zinc fingers of PRDM9 were found to be diversifying under positive selection and concerted evolution, specifically at positions involved in defining their DNA-binding specificity, consistent with the reports from other vertebrates such as mice, humans, equids and chimpanzees. ZF-encoding regions of the PRDM7, a paralog of PRDM9 on bovine chromosome 22 and on unknown chromosomes in other studied species were found to contain 84 base repeat units as in PRDM9, but there were multiple disruptive mutations after the first repeat unit. The diversity of the ZFs suggests that PRDM9 may activate recombination hotspots that are largely unique to each ruminant species.

High-Resolution Mapping of Crossover and Non-crossover Recombination Events by Whole-Genome Re-sequencing of an Avian Pedigree

Recombination is an engine of genetic diversity and therefore constitutes a key process in evolutionary biology and genetics. While the outcome of crossover recombination can readily be detected as shuffled alleles by following the inheritance of markers in pedigreed families, the more precise location of both crossover and non-crossover recombination events has been difficult to pinpoint. As a consequence, we lack a detailed portrait of the recombination landscape for most organisms and knowledge on how this landscape impacts on sequence evolution at a local scale. To localize recombination events with high resolution in an avian system, we performed whole-genome re-sequencing at high coverage of a complete three-generation collared flycatcher pedigree. We identified 325 crossovers at a median resolution of 1.4 kb, with 86% of the events localized to <10 kb intervals. Observed crossover rates were in excellent agreement with data from linkage mapping, were 52% higher in male (3.56 cM/Mb) than in female meiosis (2.28 cM/Mb), and increased towards chromosome ends in male but not female meiosis. Crossover events were non-randomly distributed in the genome with several distinct hot-spots and a concentration to genic regions, with the highest density in promoters and CpG islands. We further identified 267 non-crossovers, whose location was significantly associated with crossover locations. We detected a significant transmission bias (0.18) in favour of 'strong' (G, C) over 'weak' (A, T) alleles at non-crossover events, providing direct evidence for the process of GC-biased gene conversion in an avian system. The approach taken in this study should be applicable to any species and would thereby help to provide a more comprehensive portray of the recombination landscape across organism groups.

PR Domain-containing Protein 7 (PRDM7) Is a Histone 3 Lysine 4 Trimethyltransferase

PR domain-containing protein 7 (PRDM7) is a primate-specific histone methyltransferase that is the result of a recent gene duplication of PRDM9. The two proteins are highly homologous, especially in the catalytic PR/SET domain, where they differ by only three amino acid residues. Here we report that PRDM7 is an efficient methyltransferase that selectively catalyzes the trimethylation of H3 lysine 4 (H3K4) both in vitro and in cells. Through selective mutagenesis we have dissected the functional roles of each of the three divergent residues between the PR domains of PRDM7 and PRDM9. These studies indicate that after a single serine to tyrosine mutation at residue 357 (S357Y), PRDM7 regains the substrate specificities and catalytic activities similar to its evolutionary predecessor, including the ability to efficiently methylate H3K36.

Widespread Positive Selection Drives Differentiation of Centromeric Proteins in the Drosophila melanogaster subgroup

Rapid evolution of centromeric satellite repeats is thought to cause compensatory amino acid evolution in interacting centromere-associated kinetochore proteins. Cid, a protein that mediates kinetochore/centromere interactions, displays particularly high amino acid turnover. Rapid evolution of both Cid and centromeric satellite repeats led us to hypothesize that the apparent compensatory evolution may extend to interacting partners in the Condensin I complex (i.e., SMC2, SMC4, Cap-H, Cap-D2, and Cap-G) and HP1s. Missense mutations in these proteins often result in improper centromere formation and aberrant chromosome segregation, thus selection for maintained function and coevolution among proteins of the complex is likely strong. Here, we report evidence of rapid evolution and recurrent positive selection in seven centromere-associated proteins in species of the Drosophila melanogaster subgroup, and further postulate that positive selection on these proteins could be a result of centromere drive and compensatory changes, with kinetochore proteins competing for optimal spindle attachment.

Recombining without Hotspots: A Comprehensive Evolutionary Portrait of Recombination in Two Closely Related Species of Drosophila

Meiotic recombination rate varies across the genome within and between individuals, populations, and species in virtually all taxa studied. In almost every species, this variation takes the form of discrete recombination hotspots, determined in some mammals by a protein called PRDM9. Hotspots and their determinants have a profound effect on the genomic landscape, and share certain features that extend across the tree of life. Drosophila, in contrast, are anomalous in their absence of hotspots, PRDM9, and other species-specific differences in the determination of recombination. To better understand the evolution of meiosis and general patterns of recombination across diverse taxa, we present a truly comprehensive portrait of recombination across time, combining recently published cross-based contemporary recombination estimates from each of two sister species with newly obtained linkage-disequilibrium-based historic estimates of recombination from both of these species. Using Drosophila pseudoobscura and Drosophila miranda as a model system, we compare recombination rate between species at multiple scales, and we suggest that Drosophila replicate the pattern seen in human-chimpanzee in which recombination rate is conserved at broad scales. We also find evidence of a species-wide recombination modifier(s), resulting in both a present and historic genome-wide elevation of recombination rates in D. miranda, and identify broad scale effects on recombination from the presence of an inversion. Finally, we reveal an unprecedented view of the distribution of recombination in D. pseudoobscura, illustrating patterns of linked selection and where recombination is taking place. Overall, by combining these estimation approaches, we highlight key similarities and differences in recombination between Drosophila and other organisms.

Germline methylation patterns determine the distribution of recombination events in the dog genome

The positive-regulatory domain containing nine gene, PRDM9, which strongly associates with the location of recombination events in several vertebrates, is inferred to be inactive in the dog genome. Here, we address several questions regarding the control of recombination and its influence on genome evolution in dogs. First, we address whether the association between CpG islands (CGIs) and recombination hotspots is generated by lack of methylation, GC-biased gene conversion (gBGC), or both. Using a genome-wide dog single nucleotide polymorphism data set and comparisons of the dog genome with related species, we show that recombination-associated CGIs have low CpG mutation rates, and that CpG mutation rate is negatively correlated with recombination rate genome wide, indicating that nonmethylation attracts the recombination machinery. We next use a neighbor-dependent model of nucleotide substitution to disentangle the effects of CpG mutability and gBGC and analyze the effects that loss of PRDM9 has on these rates. We infer that methylation patterns have been stable during canid genome evolution, but that dog CGIs have experienced a drastic increase in substitution rate due to gBGC, consistent with increased levels of recombination in these regions. We also show that gBGC is likely to have generated many new CGIs in the dog genome, but these mostly occur away from genes, whereas the number of CGIs in gene promoter regions has not increased greatly in recent evolutionary history. Recombination has a major impact on the distribution of CGIs that are detected in the dog genome due to the interaction between methylation and gBGC. The results indicate that germline methylation patterns are the main determinant of recombination rates in the absence of PRDM9.

Strong artificial selection in domestic mammals did not result in an increased recombination rate

Recombination rates vary in intensity and location at the species, individual, sex and chromosome levels. Despite the fundamental biological importance of this process, the selective forces that operate to shape recombination rate and patterns are unclear. Domestication offers a unique opportunity to study the interplay between recombination and selection. In domesticates, intense selection for particular traits is imposed on small populations over many generations, resulting in organisms that differ, sometimes dramatically, in morphology and physiology from their wild ancestor. Although earlier studies suggested increased recombination rate in domesticates, a formal comparison of recombination rates between domestic mammals and their wild congeners was missing. In order to determine broad-scale recombination rate, we used immunolabeling detection of MLH1 foci as crossover markers in spermatocytes in three pairs of closely related wild and domestic species (dog and wolf, goat and ibex, and sheep and mouflon). In the three pairs, and contrary to previous suggestions, our data show that contemporary recombination rate is higher in the wild species. Subsequently, we inferred recombination breakpoints in sequence data for 16 genomic regions in dogs and wolves, each containing a locus associated with a dog phenotype potentially under selection during domestication. No difference in the number and distribution of recombination breakpoints was found between dogs and wolves. We conclude that our data indicate that strong directional selection did not result in changes in recombination in domestic mammals, and that both upper and lower bounds for crossover rates may be tightly regulated.

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.

In eukaryotic genomes, recombination plays a central role by ensuring the proper segregation of chromosomes during meiosis and increasing genetic diversity at the population scale. Recombination events are not uniformly distributed along chromosomes, but cluster in narrow regions called hotspots. The absence of overlap between human and chimpanzee hotspots indicates that the location of these hotspots evolves rapidly. However, the reasons for this rapid dynamic are still unknown. To gain insight into the processes driving the evolution of recombination hotspots we analyzed the recent history of human hotspots, using the genome of a closely related archaic hominid, Denisovan. We searched for genomic signatures of past recombination activity and compared them to present-day patterns of recombination in humans. Our results show that human hotspots are younger than previously thought and that they are not conserved in Denisovans. Moreover, we confirm that hotspots are subject to a self-destruction process, due to biased gene conversion. We quantified this process, and showed that its intensity is strong enough to cause the fast turnover of human hotspots.

Unraveling recombination rate evolution using ancestral recombination maps

Recombination maps of ancestral species can be constructed from comparative analyses of genomes from closely related species, exemplified by a recently published map of the human-chimpanzee ancestor. Such maps resolve differences in recombination rate between species into changes along individual branches in the speciation tree, and allow identification of associated changes in the genomic sequences. We describe how coalescent hidden Markov models are able to call individual recombination events in ancestral species through inference of incomplete lineage sorting along a genomic alignment. In the great apes, speciation events are sufficiently close in time that a map can be inferred for the ancestral species at each internal branch - allowing evolution of recombination rate to be tracked over evolutionary time scales from speciation event to speciation event. We see this approach as a way of characterizing the evolution of recombination rate and the genomic properties that influence it.

A high-density linkage map enables a second-generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high-density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best-order map contained 4215 markers, with a total distance of 3132 cm and a mean genetic distance between markers of 0.12 cm. Facilitated by the array being designed to include markers from most scaffolds, we obtained a second-generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super-scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cm/Mb and correlated closely with chromosome size, from 2 cm/Mb for chromosomes >100 Mb to >10 cm/Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cm was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.

Initiation of meiotic recombination: how and where? Conservation and specificities among eukaryotes

Meiotic recombination is essential for fertility in most sexually reproducing species. This process also creates new combinations of alleles and has important consequences for genome evolution. Meiotic recombination is initiated by the formation of DNA double-strand breaks (DSBs), which are repaired by homologous recombination. DSBs are catalyzed by the evolutionarily conserved SPO11 protein, assisted by several other factors. Some of them are absolutely required, whereas others are needed only for full levels of DSB formation and may participate in the regulation of DSB timing and frequency as well as the coordination between DSB formation and repair. The sites where DSBs occur are not randomly distributed in the genome, and remarkably distinct strategies have emerged to control their localization in different species. Here, I review the recent advances in the components required for DSB formation and localization in the various model organisms in which these studies have been performed.

Transmission distortion affecting human noncrossover but not crossover recombination: a hidden source of meiotic drive

Meiotic recombination ensures the correct segregation of homologous chromosomes during gamete formation and contributes to DNA diversity through both large-scale reciprocal crossovers and very localised gene conversion events, also known as noncrossovers. Considerable progress has been made in understanding factors such as PRDM9 and SNP variants that influence the initiation of recombination at human hotspots but very little is known about factors acting downstream. To address this, we simultaneously analysed both types of recombinant molecule in sperm DNA at six highly active hotspots, and looked for disparity in the transmission of allelic variants indicative of any cis-acting influences. At two of the hotspots we identified a novel form of biased transmission that was exclusive to the noncrossover class of recombinant, and which presumably arises through differences between crossovers and noncrossovers in heteroduplex formation and biased mismatch repair. This form of biased gene conversion is not predicted to influence hotspot activity as previously noted for SNPs that affect recombination initiation, but does constitute a powerful and previously undetected source of recombination-driven meiotic drive that by extrapolation may affect thousands of recombination hotspots throughout the human genome. Intriguingly, at both of the hotspots described here, this drive favours strong (G/C) over weak (A/T) base pairs as might be predicted from the well-established correlations between high GC content and recombination activity in mammalian genomes.

The B-subdomain of the Xenopus laevis XFIN KRAB-AB domain is responsible for its weaker transcriptional repressor activity compared to human ZNF10/Kox1

The Krüppel-associated box (KRAB) domain interacts with the nuclear hub protein TRIM28 to initiate or mediate chromatin-dependent processes like transcriptional repression, imprinting or suppression of endogenous retroviruses. The prototype KRAB domain initially identified in ZNF10/KOX1 encompasses two subdomains A and B that are found in hundreds of zinc finger transcription factors studied in human and murine genomes. Here we demonstrate for the first time transcriptional repressor activity of an amphibian KRAB domain. After sequence correction, the updated KRAB-AB domain of zinc finger protein XFIN from the frog Xenopus laevis was found to confer transcriptional repression in reporter assays in Xenopus laevis A6 kidney cells as well as in human HeLa, but not in the minnow Pimephales promelas fish cell line EPC. Binding of the XFIN KRAB-AB domain to human TRIM28 was demonstrated in a classical co-immunoprecipitation approach and visualized in a single-cell compartmentalization assay. XFIN-AB displayed reduced potency in repression as well as lower strength of interaction with TRIM28 compared to ZNF10 KRAB-AB. KRAB-B subdomain swapping between the two KRAB domains indicated that it was mainly the KRAB-B subdomain of XFIN that was responsible for its lower capacity in repression and binding to human TRIM28. In EPC fish cells, ZNF10 and XFIN KRAB repressor activity could be partially restored to low levels by adding exogenous human TRIM28. In contrast to XFIN, we did not find any transcriptional repression activity for the KRAB-like domain of human PRDM9 in HeLa cells. PRDM9 is thought to harbor an evolutionary older domain related to KRAB whose homologs even occur in invertebrates. Our results support the notion that functional bona fide KRAB domains which confer transcriptional repression and interact with TRIM28 most likely co-evolved together with TRIM28 at the beginning of tetrapode evolution.

Diversity of Prdm9 zinc finger array in wild mice unravels new facets of the evolutionary turnover of this coding minisatellite

In humans and mice, meiotic recombination events cluster into narrow hotspots whose genomic positions are defined by the PRDM9 protein via its DNA binding domain constituted of an array of zinc fingers (ZnFs). High polymorphism and rapid divergence of the Prdm9 gene ZnF domain appear to involve positive selection at DNA-recognition amino-acid positions, but the nature of the underlying evolutionary pressures remains a puzzle. Here we explore the variability of the Prdm9 ZnF array in wild mice, and uncovered a high allelic diversity of both ZnF copy number and identity with the caracterization of 113 alleles. We analyze features of the diversity of ZnF identity which is mostly due to non-synonymous changes at codons -1, 3 and 6 of each ZnF, corresponding to amino-acids involved in DNA binding. Using methods adapted to the minisatellite structure of the ZnF array, we infer a phylogenetic tree of these alleles. We find the sister species Mus spicilegus and M. macedonicus as well as the three house mouse (Mus musculus) subspecies to be polyphyletic. However some sublineages have expanded independently in Mus musculus musculus and M. m. domesticus, the latter further showing phylogeographic substructure. Compared to random genomic regions and non-coding minisatellites, none of these patterns appears exceptional. In silico prediction of DNA binding sites for each allele, overlap of their alignments to the genome and relative coverage of the different families of interspersed repeated elements suggest a large diversity between PRDM9 variants with a potential for highly divergent distributions of recombination events in the genome with little correlation to evolutionary distance. By compiling PRDM9 ZnF protein sequences in Primates, Muridae and Equids, we find different diversity patterns among the three amino-acids most critical for the DNA-recognition function, suggesting different diversification timescales.

Meiotic recombination in mammals: localization and regulation

During meiosis, a programmed induction of DNA double-strand breaks (DSBs) leads to the exchange of genetic material between homologous chromosomes. These exchanges increase genome diversity and are essential for proper chromosome segregation at the first meiotic division. Recent findings have highlighted an unexpected molecular control of the distribution of meiotic DSBs in mammals by a rapidly evolving gene, PR domain-containing 9 (PRDM9), and genome-wide analyses have facilitated the characterization of meiotic DSB sites at unprecedented resolution. In addition, the identification of new players in DSB repair processes has allowed the delineation of recombination pathways that have two major outcomes, crossovers and non-crossovers, which have distinct mechanistic roles and consequences for genome evolution.

Meiotic recombination strongly influences GC-content evolution in short regions in the mouse genome

Meiotic recombination is known to influence GC-content evolution in large regions of mammalian genomes by favoring the fixation of G and C alleles and increasing the rate of A/T to G/C substitutions. This process is known as GC-biased gene conversion (gBGC). Until recently, genome-wide measures of fine-scale recombination activity were unavailable in mice. Additionally, comparative studies focusing on mouse were limited as the closest organism with its genome fully sequenced was rat. Here, we make use of the recent mapping of double strand breaks (DSBs), the first step of meiotic recombination, in the mouse genome and of the sequencing of mouse closely related subspecies to analyze the fine-scale evolutionary signature of meiotic recombination on GC-content evolution in recombination hotspots, short regions that undergo extreme rates of recombination. We measure substitution rates around DSB hotspots and observe that gBGC is affecting a very short region (≈ 1 kbp) in length around these hotspots. Furthermore, we can infer that the locations of hotspots evolved rapidly during mouse evolution.

The role of gene conversion in preserving rearrangement hotspots in the human genome

Hotspots of non-allelic homologous recombination (NAHR) have a crucial role in creating genetic diversity and are also associated with dozens of genomic disorders. Recent studies suggest that many human NAHR hotspots have been preserved throughout the evolution of primates. NAHR hotspots are likely to remain active as long as the segmental duplications (SDs) promoting NAHR retain sufficient similarity. Here, we propose an evolutionary model of SDs that incorporates the effect of gene conversion and compare it with a null model that assumes SDs evolve independently without gene conversion. The gene conversion model predicts a much longer lifespan of NAHR hotspots compared with the null model. We show that the literature on copy number variants (CNVs) and genomic disorders, and also the results of additional analysis of CNVs, are all more consistent with the gene conversion model.

Evolution of recombination and genome structure in eusocial insects

Eusocial Hymenoptera, such as the European honey bee, Apis mellifera, have the highest recombination rates of multicellular animals.(1) Recently, we showed(2) that a side-effect of recombination in the honey bee, GC biased gene conversion (bGC), helps maintain the unusual bimodal GC-content distribution of the bee genome by increasing GC-content in high recombination areas while low recombination areas are losing GC-content because of biased AT mutations and low rates of bGC. Although the very high recombination rate of A. mellifera makes GC-content evolution easier to study, the pattern is consistent with results found in many other species including mammals and yeast.(3) Also consistent across phyla is the association of higher genetic diversity and divergence with high GC and high recombination areas.(4) (,) (5) Finally, we showed that genes overexpressed in the brains of workers cluster in GC-rich genomic areas with the highest rates of recombination and molecular evolution.(2) In this Addendum we present a conceptual model of how eusociality and high recombination rates may co-evolve.

Characterization of Prdm9 in equids and sterility in mules

Prdm9 (Meisetz) is the first speciation gene discovered in vertebrates conferring reproductive isolation. This locus encodes a meiosis-specific histone H3 methyltransferase that specifies meiotic recombination hotspots during gametogenesis. Allelic differences in Prdm9, characterized for a variable number of zinc finger (ZF) domains, have been associated with hybrid sterility in male house mice via spermatogenic failure at the pachytene stage. The mule, a classic example of hybrid sterility in mammals also exhibits a similar spermatogenesis breakdown, making Prdm9 an interesting candidate to evaluate in equine hybrids. In this study, we characterized the Prdm9 gene in all species of equids by analyzing sequence variation of the ZF domains and estimating positive selection. We also evaluated the role of Prdm9 in hybrid sterility by assessing allelic differences of ZF domains in equine hybrids. We found remarkable variation in the sequence and number of ZF domains among equid species, ranging from five domains in the Tibetan kiang and Asiatic wild ass, to 14 in the Grevy’s zebra. Positive selection was detected in all species at amino acid sites known to be associated with DNA-binding specificity of ZF domains in mice and humans. Equine hybrids, in particular a quartet pedigree composed of a fertile mule showed a mosaic of sequences and number of ZF domains suggesting that Prdm9 variation does not seem by itself to contribute to equine hybrid sterility.

Programmed induction of DNA double strand breaks during meiosis: setting up communication between DNA and the chromosome structure

During the first meiotic prophase, hundreds of DNA double strand breaks (DSBs) are deliberately self-inflicted along chromosomes in order to promote homologous recombination between homologs. These DSBs, catalyzed by the evolutionary conserved Spo11 protein, are highly regulated. Recent studies in yeast and mammals have identified key components involved in meiotic DSB formation. In mammals, the DNA binding specificity of PRDM9 determines where DSB occur, whereas in yeast, Spo11 acts in regions which one important feature is chromatin accessibility. However, DSB formation requires additional proteins located on chromosome axes, and the Saccharomyces cerevisiae protein, Spp1 has been recently identified to make the link between axes and DSB sites. These recent findings open exciting routes to understanding how the requirement to regulate DSBs along and between homologs is achieved.

A prominent role of KRAB-ZNF transcription factors in mammalian speciation?

The mechanisms of speciation have been one of the most debated topics in evolutionary biology. Among all reproductive barriers, postzygotic reproductive isolation is perhaps the one that has attracted the most attention from geneticists. Despite remarkable advances in the identification of loci involved in Drosophila speciation, little is known about the genes, functions, and biochemical interactions of the molecules underlying hybrid sterility and inviability in mammals. Here, we discuss the main evolutionary and molecular features that make transcription factors (TFs), especially the family of zinc finger proteins with a Krüppel-associated box domain (KRAB-ZNF), strong candidates to play an important role in postzygotic reproductive isolation. Motivated by the recent identification of the gene encoding PR domain zinc finger protein 9 (Prdm9; a KRAB-ZNF gene) as the first hybrid sterility gene identified in mammals, we further propose integrative approaches to study KRAB-ZNF genes with the main goal of characterizing the molecular pathways and interactions involved in hybrid incompatibilities.

Zinc finger binding motifs do not explain recombination rate variation within or between species of Drosophila

In humans and mice, the Cys(2)His(2) 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(2)His(2) zinc fingers to predict nucleotide binding motifs for all Cys(2)His(2) 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.

Epigenetic remodeling of meiotic crossover frequency in Arabidopsis thaliana DNA methyltransferase mutants

Meiosis is a specialized eukaryotic cell division that generates haploid gametes required for sexual reproduction. During meiosis, homologous chromosomes pair and undergo reciprocal genetic exchange, termed crossover (CO). Meiotic CO frequency varies along the physical length of chromosomes and is determined by hierarchical mechanisms, including epigenetic organization, for example methylation of the DNA and histones. Here we investigate the role of DNA methylation in determining patterns of CO frequency along Arabidopsis thaliana chromosomes. In A. thaliana the pericentromeric regions are repetitive, densely DNA methylated, and suppressed for both RNA polymerase-II transcription and CO frequency. DNA hypomethylated methyltransferase1 (met1) mutants show transcriptional reactivation of repetitive sequences in the pericentromeres, which we demonstrate is coupled to extensive remodeling of CO frequency. We observe elevated centromere-proximal COs in met1, coincident with pericentromeric decreases and distal increases. Importantly, total numbers of CO events are similar between wild type and met1, suggesting a role for interference and homeostasis in CO remodeling. To understand recombination distributions at a finer scale we generated CO frequency maps close to the telomere of chromosome 3 in wild type and demonstrate an elevated recombination topology in met1. Using a pollen-typing strategy we have identified an intergenic nucleosome-free CO hotspot 3a, and we demonstrate that it undergoes increased recombination activity in met1. We hypothesize that modulation of 3a activity is caused by CO remodeling driven by elevated centromeric COs. These data demonstrate how regional epigenetic organization can pattern recombination frequency along eukaryotic chromosomes.

Genetic variants in REC8, RNF212, and PRDM9 influence male recombination in cattle

We use >250,000 cross-over events identified in >10,000 bovine sperm cells to perform an extensive characterization of meiotic recombination in male cattle. We map Quantitative Trait Loci (QTL) influencing genome-wide recombination rate, genome-wide hotspot usage, and locus-specific recombination rate. We fine-map three QTL and present strong evidence that genetic variants in REC8 and RNF212 influence genome-wide recombination rate, while genetic variants in PRDM9 influence genome-wide hotspot usage.

Homologous recombination is an essential cellular process that determines proper chromosome segregation during meiosis, affects fertility, and influences evolvability. Nevertheless, the components of the recombination apparatus remain incompletely characterized in mammals. One approach to identify such components is to identify the genes that underlie inherited variation in recombination phenotypes. In addition to providing mechanistic insights, this would allow the study of the evolutionary forces that shape the recombination process. In this paper, we take advantage of genotypes for 50,000 genome-wide SNP markers to measure four recombination phenotypes (genome-wide recombination rate, genome-wide hotspot usage, locus-specific recombination rate, genome-wide cross-over interference) for >750 bulls on the basis of >250,000 cross-overs detected in sperm cells transmitted to >10,000 sons. We quantify the heritability and scan the genome for Quantitative Trait Loci (QTL) influencing each one of these recombination phenotypes. We perform a detailed genetic analysis of three such QTL, thereby providing evidence that genetic variants in REC8 and RNF212 influence genome-wide recombination rate, while genetic variants in an X-linked PRDM9 paralogue influence genome-wide hotspot usage.