Drive against hotspot motifs in primates implicates the PRDM9 gene in meiotic recombination

Unexpectedly low recombination rates and presence of hotspots in termite genomes

Meiotic recombination is a fundamental evolutionary process that facilitates adaptation and the removal of deleterious genetic variation. Social Hymenoptera exhibit some of the highest recombination rates among metazoans, whereas high recombination rates have not been found among nonsocial species from this insect order. It is unknown whether elevated recombination rates are a ubiquitous feature of all social insects. In many metazoan taxa, recombination is mainly restricted to hotspots a few kilobases in length. However, little is known about the prevalence of recombination hotspots in insect genomes. Here we infer recombination rate and its fine-scale variation across the genomes of two social species from the insect order Blattodea: the termites Macrotermes bellicosus and Cryptotermes secundus We used linkage disequilibrium-based methods to infer recombination rate. We infer that recombination rates are close to 1 cM/Mb in both species, similar to the average metazoan rate. We also observe a highly punctate distribution of recombination in both termite genomes, indicative of the presence of recombination hotspots. We infer the presence of full-length PRDM9 genes in the genomes of both species, which suggests recombination hotspots in termites might be determined by PRDM9, as they are in mammals. We also find that recombination rates in genes are correlated with inferred levels of germline DNA methylation. The finding of low recombination rates in termites indicates that eusociality is not universally connected to elevated recombination rate. We speculate that the elevated recombination rates in social Hymenoptera are instead promoted by intense selection among haploid males.

A Pedigree-Based Map of Crossovers and Noncrossovers in Aye-Ayes (Daubentonia madagascariensis)

Gaining a better understanding of the rates and patterns of meiotic recombination is crucial for improving evolutionary genomic modeling, with applications ranging from demographic to selective inference. Although previous research has provided important insights into the landscape of crossovers in humans and other haplorrhines, our understanding of both the considerably more common outcome of recombination (i.e. noncrossovers) as well as the landscapes in more distantly related primates (i.e. strepsirrhines) remains limited owing to difficulties associated with both the identification of noncrossover tracts as well as species sampling. Thus, in order to elucidate recombination patterns in this understudied branch of the primate clade, we here characterize crossover and noncrossover landscapes in aye-ayes utilizing whole-genome sequencing data from six three-generation pedigrees and three two-generation multi-sibling families, and in so doing provide novel insights into this important evolutionary process shaping genomic diversity in one of the world's most critically endangered primate species.

Comparative Phylogenetics Reveal Clade-specific Drivers of Recombination Rate Evolution Across Vertebrates

Meiotic recombination is an integral cellular process, required for the production of viable gametes. Recombination rate is a fundamental genomic parameter, modulating genomic responses to selection. Our increasingly detailed understanding of its molecular underpinnings raises the prospect that we can gain insight into trait divergence by examining the molecular evolution of recombination genes from a pathway perspective, as in mammals, where protein-coding changes in later stages of the recombination pathway are connected to divergence in intra-clade recombination rate. Here, we leverage increased availability of avian and teleost genomes to reconstruct the evolution of the recombination pathway across two additional vertebrate clades: birds, which have higher and more variable rates of recombination and similar divergence times to mammals, and teleost fish, which have much deeper divergence times. Rates of molecular evolution of recombination genes are highly correlated between vertebrate clades and significantly elevated compared to control panels, suggesting that they experience similar selective pressures. Avian recombination genes are significantly more likely to exhibit signatures of positive selection than other clades, unrestricted to later stages of the pathway. Signatures of positive selection in genes linked to recombination rate variation in mammalian populations and those with signatures of positive selection across the avian phylogeny are highly correlated. In contrast, teleost fish recombination genes have significantly less evidence of positive selection despite high intra-clade recombination rate variability. Gaining clade-specific understanding of patterns of variation in recombination genes can elucidate drivers of recombination rate and thus, factors influencing genetic diversity, selection efficacy, and species divergence.

Mixed Outcomes in Recombination Rates After Domestication: Revisiting Theory and Data

The process of domestication has altered many phenotypes. Selection on these phenotypes has long been hypothesised to indirectly select for increases in the genome-wide recombination rate. This hypothesis is potentially consistent with theory on the evolution of the recombination rate, but empirical support has been unclear. We review relevant theory, lab-based experiments, and data comparing recombination rates in wild progenitors and their domesticated counterparts. We utilise population sequencing data and a deep learning method to infer genome-wide recombination rates for new comparisons of chicken/red junglefowl, sheep/mouflon, and goat/bezoar. We find evidence of increased recombination in domesticated goats compared to bezoars but more mixed results in chicken and generally decreased recombination in domesticated sheep compared to mouflon. Our results add to a growing body of literature in plants and animals that finds no consistent evidence of an increase in genome-wide recombination with domestication.

Phenogenomic resources immortalized in a panel of wild-derived strains of five species of house mice

The house mouse, Mus musculus, is a widely used animal model in biomedical research, with classical laboratory strains (CLS) being the most frequently employed. However, the limited genetic variability in CLS hinders their applicability in evolutionary studies. Wild-derived strains (WDS), on the other hand, provide a suitable resource for such investigations. This study quantifies genetic and phenotypic data of 101 WDS representing 5 species, 3 subspecies, and 8 natural Y consomic strains and compares them with CLS. Genetic variability was estimated using whole mtDNA sequences, the Prdm9 gene, and copy number variation at two sex chromosome-linked genes. WDS exhibit a large natural variation with up to 2173 polymorphic sites in mitogenomes, whereas CLS display 92 sites. Moreover, while CLS have two Prdm9 alleles, WDS harbour 46 different alleles. Although CLS resemble M. m. domesticus and M. m. musculus WDS, they differ from them in 10 and 14 out of 16 phenotypic traits, respectively. The results suggest that WDS can be a useful tool in evolutionary and biomedical studies with great potential for medical applications.

Keeping it safe: control of meiotic chromosome breakage

Meiotic cells introduce numerous programmed DNA double-strand breaks (DSBs) into their genome to stimulate crossover recombination. DSB numbers must be high enough to ensure each homologous chromosome pair receives the obligate crossover required for accurate meiotic chromosome segregation. However, every DSB also increases the risk of aberrant or incomplete DNA repair, and thus genome instability. To mitigate these risks, meiotic cells have evolved an intricate network of controls that modulates the timing, levels, and genomic location of meiotic DSBs. This Review summarizes our current understanding of these controls with a particular focus on the mechanisms that prevent meiotic DSB formation at the wrong time or place, thereby guarding the genome from potentially catastrophic meiotic errors.

Exploring the genetic mechanisms driving KIR diversification

Killer cell immunoglobulin-like receptors (KIRs) are key modulators of natural killer cell activity, displaying either activating or inhibitory potential upon recognition of major histocompatibility complex (MHC) class I molecules. The genomic organization of KIR genes is complex, involving copy number variation and allelic polymorphism, which is probably due to their coevolution with highly polymorphic MHC ligands. The KIR diversity is reflected by more than 70 similar region configurations encountered in humans, generated through meiotic recombination events. Rhesus macaques happen to display even more diversity, and over 100 distinct configurations were identified in a relatively small cohort of animals. More than half of these region configurations feature hybrid KIR genes, suggesting a more pronounced mode of diversification in macaques. The molecular mechanism facilitating meiotic rearrangements in the KIR region is poorly understood. Examination of 21 rhesus macaque and 14 human KIR region configurations revealed the presence of long terminal repeats and PRDM9 binding motifs associated with recombination hotspots. The variable DNA recognition patterns of PRDM9 could potentially contribute to the differing recombination activities documented for the KIR region in humans and macaques. The diversification process of the KIR repertoire in natural killer cells is fundamentally distinct from the mechanisms generating T and B cell receptor diversity or MHC polymorphisms. This sophisticated recombination machinery preserves the functional integrity by the frequent generation of in-frame KIR genes. A diverse KIR repertoire contributes to the protection of individuals and populations against pathogen evasion and subversion.

Germline mutation rates and fine-scale recombination parameters in zebra finch

Most of our understanding of the fundamental processes of mutation and recombination stems from a handful of disparate model organisms and pedigree studies of mammals, with little known about other vertebrates. To gain a broader comparative perspective, we focused on the zebra finch (Taeniopygia castanotis), which, like other birds, differs from mammals in its karyotype (which includes many micro-chromosomes), in the mechanism by which recombination is directed to the genome, and in aspects of ontogenesis. We collected genome sequences from three generation pedigrees that provide information about 80 meioses, inferring 202 single-point de novo mutations, 1,088 crossovers, and 275 non-crossovers. On that basis, we estimated a sex-averaged mutation rate of 5.0 × 10-9 per base pair per generation, on par with mammals that have a similar generation time (~2-3 years). Also as in mammals, we found a paternal germline mutation bias at later stages of gametogenesis (of 1.7:1) but no discernible difference between sexes in early development. Examining recombination patterns, we found that the sex-averaged crossover rate on macro-chromosomes is 0.93 cM/Mb, with a pronounced enrichment of crossovers near telomeres. In contrast, non-crossover rates are more uniformly distributed. On micro-chromosomes, sex-averaged crossover rates are substantially higher (3.96 cM/Mb), in accordance with crossover homeostasis, and both crossover and non-crossover events are more uniformly distributed. At a finer scale, recombination events overlap CpG islands more often than expected by chance, as expected in the absence of PRDM9. Estimates of the degree of GC-biased gene conversion (59%), the mean non-crossover conversion tract length (~32 bp), and the non-crossover-to-crossover ratio (5.4:1) are all comparable to those reported in primates and mice. Therefore, properties of germline mutation and recombination resolutions remain similar over large phylogenetic distances.

High-recombining genomic regions affect demography inference based on ancestral recombination graphs

Multiple methods of demography inference are based on the ancestral recombination graph. This powerful approach uses observed mutations to model local genealogies changing along chromosomes by historical recombination events. However, inference of underlying genealogies is difficult in regions with high recombination rate relative to mutation rate due to the lack of mutations representing genealogies. Despite the prevalence of high-recombining genomic regions in some organisms, such as birds, its impact on demography inference based on ancestral recombination graphs has not been well studied. Here, we use population genomic simulations to investigate the impact of high-recombining regions on demography inference based on ancestral recombination graphs. We demonstrate that inference of effective population size and the time of population split events is systematically affected when high-recombining regions cover wide breadths of the chromosomes. Excluding high-recombining genomic regions can practically mitigate this impact, and population genomic inference of recombination maps is informative in defining such regions although the estimated values of local recombination rate can be biased. Finally, we confirm the relevance of our findings in empirical analysis by contrasting demography inferences applied for a bird species, the Eurasian blackcap (Sylvia atricapilla), using different parts of the genome with high and low recombination rates. Our results suggest that demography inference methods based on ancestral recombination graphs should be carried out with caution when applied in species whose genomes contain long stretches of high-recombining regions.

The KMT2 complex protein ASH2L is required for meiotic prophase progression but dispensable for mitosis in differentiated spermatogonia

ASH2L is a core component of KMT2 complexes, crucial for H3K4 trimethylation. However, its role in spermatogenesis remains elusive. Here, we demonstrate an essential role of Ash2l for meiotic prophase but dispensable for mitosis in differentiated spermatogonia. Using a germ cell-specific Ash2l knockout mouse model, we reveal that Ash2l deficiency leads to meiotic arrest and sterility in both sexes. Ash2l-deficient spermatocytes exhibit failures in chromosomal synapsis associated with persistent DMC1 foci and γH2AX, resulting in meiocyte loss due to apoptosis. Conversely, Ash2l-deficient differentiated spermatogonia show normal development. Mechanistically, Ash2l deficiency results in a global loss of H3K4me3 in promoter regions and significantly decreases expression of thousands of genes. Among these are genes involved in epigenetic silencing pathways, such as H3K9 di-methylation, DNA methylation and piRNA pathways, that are crucial for transposon repression during meiotic prophase I progression. Supporting this, we observe that Ash2l mutant spermatocytes display ectopic expression of LINE1-ORF1P. Our findings therefore reveal the previously unappreciated role of ASH2L-dependent H3K4me3 modification in spermatogenesis and provide clues to the molecular mechanisms in epigenetic disorders underlying male infertility.

KRAB zinc-finger proteins regulate endogenous retroviruses to sculpt germline transcriptomes and genome evolution

As transposable elements (TEs) coevolved with the host genome, the host genome exploited TEs as functional regulatory elements of gene expression. Here we show that a subset of KRAB domain-containing zinc-finger proteins (KZFPs), which are highly expressed in mitotically dividing spermatogonia, repress the enhancer function of endogenous retroviruses (ERVs) and that the release from KZFP-mediated repression allows activation of ERV enhancers upon entry into meiosis. This regulatory feature is observed for independently evolved KZFPs and ERVs in mice and humans, suggesting evolutionary conservation in mammals. Further, we show that KZFP-targeted ERVs are underrepresented on the sex chromosomes in meiosis, suggesting that meiotic sex chromosome inactivation (MSCI) may antagonize the coevolution of KZFPs and ERVs in mammals. Our study uncovers a mechanism by which a subset of KZFPs regulate ERVs to sculpt germline transcriptomes. We propose that epigenetic programming during the transition from mitotic spermatogonia to meiotic spermatocytes facilitates the coevolution of KZFPs and TEs on autosomes and is antagonized by MSCI.

Regulatory Plasticity of the Human Genome

Evolutionary turnover in noncoding regions has driven phenotypic divergence during past speciation events and continues to facilitate environmental adaptation through variants. We used a deep learning model to identify the substrates of regulatory turnover using genome-wide mutations mimicking three evolutionary pathways: recent history (human-chimp substitutions), modern population (human population variation), and mutational susceptibility (random mutations). We observed enhancer turnover in approximately 6% of the whole genome, with more than 80% of the novel activity arising from repurposing of enhancers between cell types. Frequency of turnover in a cell type is remarkably similar across the three pathways, despite only ∼19% overlap in the source regions. The majority of turnover loci were found to be localized within 100 kb of a gene, with the highest turnover occurring near neurodevelopmental genes including CNTNAP2, NPAS3, and AUTS2. Flanking enhancers of these genes undergo high turnover irrespective of the mutational model pathway, suggesting a high plasticity in neurocognitive evolution. Based on susceptibility to random mutations, these enhancers were identified as vulnerable by nature and feature a higher abundance of cell type-specific transcription factor binding sites. Our findings suggest that enhancer repurposing within vulnerable loci drives regulatory innovation while keeping the core regulatory networks intact.

Rapid evolution of recombination landscapes during the divergence of cichlid ecotypes in Lake Masoko

Variation of recombination rate along the genome is of crucial importance to rapid adaptation and organismal diversification. Many unknowns remain regarding how and why recombination landscapes evolve in nature. Here, we reconstruct recombination maps based on linkage disequilibrium and use subsampling and simulations to derive a new measure of recombination landscape evolution: the Population Recombination Divergence Index (PRDI). Using PRDI, we show that fine-scale recombination landscapes differ substantially between two cichlid fish ecotypes of Astatotilapia calliptera that diverged only ~2,500 generations ago. Perhaps surprisingly, recombination landscape differences are not driven by divergence in terms of allele frequency (FST) and nucleotide diversity (Δ(π)): although there is some association, we observe positive PRDI in regions where FST and Δ(π) are zero. We found a stronger association between the evolution of recombination and 47 large haplotype blocks that are polymorphic in Lake Masoko, cover 21% of the genome, and appear to include multiple inversions. Among haplotype blocks, there is a strong and clear association between the degree of recombination divergence and differences between ecotypes in heterozygosity, consistent with recombination suppression in heterozygotes. Overall, our work provides a holistic view of changes in population recombination landscapes during the early stages of speciation with gene flow.

Young KRAB-zinc finger gene clusters are highly dynamic incubators of ERV-driven genetic heterogeneity in mice

KRAB-zinc finger proteins (KZFPs) comprise the largest family of mammalian transcription factors, rapidly evolving within and between species. Most KZFPs repress endogenous retroviruses (ERVs) and other retrotransposons, with KZFP gene numbers correlating with the ERV load across species, suggesting coevolution. How new KZFPs emerge in response to ERV invasions is currently unknown. Using a combination of long-read sequencing technologies and genome assembly, we present a first detailed comparative analysis of young KZFP gene clusters in the mouse lineage, which has undergone recent KZFP gene expansion and ERV infiltration. Detailed annotation of KZFP genes in a cluster on Mus musculus Chromosome 4 revealed parallel expansion and diversification of this locus in different mouse strains (C57BL/6J, 129S1/SvImJ and CAST/EiJ) and species (Mus spretus and Mus pahari). Our data supports a model by which new ERV integrations within young KZFP gene clusters likely promoted recombination events leading to the emergence of new KZFPs that repress them. At the same time, ERVs also increased their numbers by duplication instead of retrotransposition alone, unraveling a new mechanism for ERV enrichment at these loci.

In vitro reconstitution of meiotic DNA double-strand-break formation

The Spo11 complex catalyses the formation of DNA double-strand breaks (DSBs), initiating meiotic recombination-a process that is essential for fertility and genetic diversity1,2. Although the function of Spo11 has been known for 27 years, previous efforts to reconstitute DSB formation in vitro have been unsuccessful. Here we biochemically characterize the mouse SPO11-TOP6BL protein complex, and show that this complex cleaves DNA and covalently attaches to the 5' terminus of DNA breaks in vitro. Using a point-mutation strategy, we reveal that Mg2+ is essential for the DNA-cleavage activity of this complex in vitro, as confirmed by knock-in mice carrying a point mutation in SPO11 that disrupts its binding to Mg2+, thereby abolishing DSB formation. However, the activity of the SPO11 complex is ATP-independent. We also present evidence that the mouse SPO11 complex is biochemically distinct from the ancestral topoisomerase VI. Our findings establish a mechanistic framework for understanding the first steps of meiotic recombination.

SPO11 dimers are sufficient to catalyse DNA double-strand breaks in vitro

SPO11 initiates meiotic recombination through the induction of programmed DNA double-strand breaks (DSBs)1,2, but this catalytic activity has never been reconstituted in vitro3,4. Here, using Mus musculus SPO11, we report a biochemical system that recapitulates all the hallmarks of meiotic DSB formation. We show that SPO11 catalyses break formation in the absence of any partners and remains covalently attached to the 5' broken strands. We find that target site selection by SPO11 is influenced by the sequence, bendability and topology of the DNA substrate, and provide evidence that SPO11 can reseal single-strand DNA breaks. In addition, we show that SPO11 is monomeric in solution and that cleavage requires dimerization for the reconstitution of two hybrid active sites. SPO11 and its partner TOP6BL form a 1:1 complex that catalyses DNA cleavage with an activity similar to that of SPO11 alone. However, this complex binds DNA ends with higher affinity, suggesting a potential role after cleavage. We propose a model in which additional partners of SPO11 required for DSB formation in vivo assemble biomolecular condensates that recruit SPO11-TOP6BL, enabling dimerization and cleavage. Our work establishes SPO11 dimerization as the fundamental mechanism that controls the induction of meiotic DSBs.

Mutation and recombination parameters in zebra finch are similar to those in mammals

Most of our understanding of the fundamental processes of mutation and recombination stems from a handful of disparate model organisms and pedigree studies of mammals, with little known about other vertebrates. To gain a broader comparative perspective, we focused on the zebra finch (Taeniopygia castanotis), which, like other birds, differs from mammals in its karyotype (which includes many micro-chromosomes), in the mechanism by which recombination is directed to the genome, and in aspects of ontogenesis. We collected genome sequences from three generation pedigrees that provide information about 80 meioses, inferring 202 single-point de novo mutations, 1,174 crossovers, and 275 non-crossovers. On that basis, we estimated a sex-averaged mutation rate of 5.0 × 10-9 per base pair per generation, on par with mammals that have a similar generation time (~2-3 years). Also as in mammals, we found a paternal germline mutation bias at later stages of gametogenesis (of 1.7:1) but no discernible difference between sexes in early development. Examining recombination patterns, we found that the sex-averaged crossover rate on macro-chromosomes (1.05 cM/Mb) is again similar to values observed in mammals, as is the spatial distribution of crossovers, with a pronounced enrichment near telomeres. In contrast, non-crossover rates are more uniformly distributed. On micro-chromosomes, sex-averaged crossover rates are substantially higher (4.21 cM/Mb), as expected from crossover homeostasis, and both crossover and non-crossover events are more uniformly distributed. At a finer scale, recombination events overlap CpG islands more often than expected by chance, as expected in the absence of PRDM9. Despite differences in the mechanism by which recombination events are specified and the presence of many micro-chromosomes, estimates of the degree of GC-biased gene conversion (59%), the mean non-crossover conversion tract length (~32 bp), and the non-crossover-to-crossover ratio (5.4:1) are all comparable to those reported in primates and mice. The similarity of mutation and recombination properties in zebra finch to those in mammals suggest that they are conserved by natural selection.

Variability of PRDM9 in buffaloes

The buffalo population raised in Brazil tend to show loss of genetic variability over generations, with significant estimates of inbreeding depression. Besides mating genetically distant individuals, other tools can be used to maintain/increase the genetic variability of the population, such as the use of PRDM9 genotypes. The PRDM9 gene promotes the creation of crossing-over points across the genome, with each allele promoting the creation of a different hotspot. Thus, increasing the frequency of less frequent alleles in the population, allows the emergence of new haplotypes and increases genetic variability. So, this study aimed to characterize the alleles of the PRDM9 gene circulating in the Murrah, Jaffarabadi, and Mediterranean breeds and verify their potential impact on genetic diversity management within the populations. The three alleles (B, C and D) were found in the three breeds at different frequencies, as well as the genotypic frequencies. The mating of different homozygous genotypes and genotypes carrying less frequent alleles may increase recombination rates and population variability. Four described variants and one new variant for allele D were found by sequencing. It was verified that it is possible to mate sires and dams with different PRDM9 genotypes in order to try to increase genetic variability in buffalo populations, improving the matings choices in buffalo breeding, helping to maintain production levels.

The recombination landscape of the barn owl, from families to populations

Homologous recombination is a meiotic process that generates diversity along the genome and interacts with all evolutionary forces. Despite its importance, studies of recombination landscapes are lacking due to methodological limitations and limited data. Frequently used approaches include linkage mapping based on familial data that provides sex-specific broad-scale estimates of realized recombination and inferences based on population linkage disequilibrium that reveal a more fine-scale resolution of the recombination landscape, albeit dependent on the effective population size and the selective forces acting on the population. In this study, we use a combination of these 2 methods to elucidate the recombination landscape for the Afro-European barn owl (Tyto alba). We find subtle differences in crossover placement between sexes that lead to differential effective shuffling of alleles. Linkage disequilibrium-based estimates of recombination are concordant with family-based estimates and identify large variation in recombination rates within and among linkage groups. Larger chromosomes show variation in recombination rates, while smaller chromosomes have a universally high rate that shapes the diversity landscape. We find that recombination rates are correlated with gene content, genetic diversity, and GC content. We find no conclusive differences in the recombination landscapes between populations. Overall, this comprehensive analysis enhances our understanding of recombination dynamics, genomic architecture, and sex-specific variation in the barn owl, contributing valuable insights to the broader field of avian genomics.

Inferring fine-scale mutation and recombination rate maps in aye-ayes ( Daubentonia madagascariensis )

The rate of input of new genetic mutations, and the rate at which that variation is reshuffled, are key evolutionary processes shaping genomic diversity. Importantly, these rates vary not just across populations and species, but also across individual genomes. Despite previous studies having demonstrated that failing to account for rate heterogeneity across the genome can bias the inference of both selective and neutral population genetic processes, mutation and recombination rate maps have to date only been generated for a relatively small number of organisms. Here, we infer such fine-scale maps for the aye-aye ( Daubentonia madagascariensis ) - a highly endangered strepsirrhine that represents one of the earliest splits in the primate clade, and thus stands as an important outgroup to the more commonly-studied haplorrhines - utilizing a recently released fully-annotated genome combined with high-quality population sequencing data. We compare our indirectly inferred rates to previous pedigree-based estimates, finding further evidence of relatively low mutation and recombination rates in aye-ayes compared to other primates.

Genetics of Recombination Rate Variation Within and Between Species

Recombination diversifies the genomes of offspring, influences the evolutionary dynamics of populations, and ensures that chromosomes segregate properly during meiosis. Individuals recombine at different rates but observed levels of variation in recombination rate remain mostly unexplained. Genetic dissection of differences in recombination rate within and between species provides a powerful framework for understanding how this trait evolves. In this Perspective, I amalgamate published findings from genetic studies of variation in the genome-wide number of crossovers within and between species, and I use exploratory analyses to identify preliminary patterns. The narrow-sense heritability of crossover count is consistently low, indicating limited resemblance among relatives and predicting a weak response to short-term selection. Variants associated with crossover number within populations span the range of minor allele frequency. The size of the additive effect of recombination-associated variants, along with a negative correlation between this effect and minor allele frequency, raises the prospect that mutations inducing phenotypic shifts larger than a few crossovers are deleterious, though the contributions of methodological biases to these patterns deserve investigation. Quantitative trait loci that contribute to differences between populations or species alter crossover number in both directions, a pattern inconsistent with selection toward a constant optimum for this trait. Building on this characterization of genetic variation in crossover number within and between species, I describe fruitful avenues for future research. Better integrating recombination rate into quantitative genetics will reveal the balance of evolutionary forces responsible for genetic variation in this trait that shapes inheritance.

Comparative genetic mapping and a consensus interspecific genetic map reveal strong synteny and collinearity within the Citrus genus

Introduction

Useful germplasm for citrus breeding includes all sexually compatible species of the former genera Citrus, Clymenia, Eremocitrus, Fortunella, Microcitrus, Oxanthera, and Poncirus, now merged in the single Citrus genus. An improved knowledge on the synteny/collinearity between the genome of these different species, and on their recombination landscapes, is essential to optimize interspecific breeding schemes.

Method

We have performed a large comparative genetic mapping study including several main clades of the Citrus genus. It concerns five species (C. maxima, C. medica, C. reticulata, C. trifoliata and C. glauca), two horticultural groups resulting from interspecific admixture (clementine and lemon) and two recent interspecific hybrids (C. australis x C. australasica and C. maxima x C. reticulata). The nine individual genetic maps were established from GBS data of 1,216 hybrids.

Results and discussion

The number of SNPs mapped for each parent varies from 760 for C. medica to 4,436 for the C. maxima x C. reticulata hybrid, with an average of 2,162.3 markers by map. Their comparison with C. clementina v1.0 assembly and inter-map comparisons revealed a high synteny and collinearity between the nine genetic maps. Non-Mendelian segregation was frequent and specific for each parental combination. The recombination landscape was similar for the nine mapped parents, and large genomic regions with very low recombination were identified. A consensus genetic map was successfully established. It encompasses 10,756 loci, including 7,915 gene-based markers and 2,841 non-genic SNPs. The anchoring of the consensus map on 15 published citrus chromosome-scale genome assemblies revealed a high synteny and collinearity for the most recent assemblies, whereas discrepancies were observed for some older ones. Large structural variations do not seem to have played a major role in the differentiation of the main species of the Citrus genus. The consensus genetic map is a useful tool to check the accuracy of genome assemblies, identify large structural variation and focus on analyzing potential relationships with phenotypic variations. It should also be a reference framework to integrate the positions of QTLs and useful genes identified in different analyses.

Multiple losses of aKRAB from PRDM9 coincide with a teleost-specific intron size distribution

BACKGROUND

Primary transcripts are largely comprised of intronic sequences that are excised and discarded shortly after synthesis. In vertebrates, the shape of the intron size distribution is largely constant; however, most teleost fish have a diverged log-bimodal 'teleost distribution' (TD) that is seen only in teleosts. How the TD evolved and to what extent this was affected by adaptative or non-adaptive mechanisms is unknown.

RESULTS

Here, we show that the TD has evolved independently at least six times and that its appearance is linked to the loss of the aKRAB domain from PRDM9. We determined intron size distributions and identified PRDM9 orthologues from annotated genomes in addition to scanning 1193 teleost assemblies for the aKRAB domain. We show that a diverged form of PRDM9 ( β ) is predominant in teleosts whereas the α version is absent from most species. Only a subset of PRDM9- α proteins contain aKRAB, and hence, it is present only in a small number of teleost lineages. Almost all lineages lacking aKRAB (but no species with) had TDs.

CONCLUSIONS

In mammals, PRDM9 defines the sites of meiotic recombination through a mechanism that increases structural variance and depends on aKRAB. The loss of aKRAB is likely to have shifted the locations of both recombination and structural variance hotspots. Our observations suggest that the TD evolved as a side-effect of these changes and link recombination to the evolution of intron size illustrating how genome architectures can evolve in the absence of selection.

Reconstruction of the human amylase locus reveals ancient duplications seeding modern-day variation

Previous studies suggested that the copy number of the human salivary amylase gene, AMY1, correlates with starch-rich diets. However, evolutionary analyses are hampered by the absence of accurate, sequence-resolved haplotype variation maps. We identified 30 structurally distinct haplotypes at nucleotide resolution among 98 present-day humans, revealing that the coding sequences of AMY1 copies are evolving under negative selection. Genomic analyses of these haplotypes in archaic hominins and ancient human genomes suggest that a common three-copy haplotype, dating as far back as 800,000 years ago, has seeded rapidly evolving rearrangements through recurrent nonallelic homologous recombination. Additionally, haplotypes with more than three AMY1 copies have significantly increased in frequency among European farmers over the past 4000 years, potentially as an adaptive response to increased starch digestion.

A Minimal Hybrid Sterility Genome Assembled by Chromosome Swapping Between Mouse Subspecies (Mus musculus)

Hybrid sterility is a reproductive isolation barrier between diverging taxa securing the early steps of speciation. Hybrid sterility is ubiquitous in the animal and plant kingdoms, but its genetic control is poorly understood. In our previous studies, we have uncovered the sterility of hybrids between musculus and domesticus subspecies of the house mouse, which is controlled by the Prdm9 gene, the X-linked Hstx2 locus, and subspecific heterozygosity for genetic background. To further investigate this form of genic-driven chromosomal sterility, we constructed a simplified hybrid sterility model within the genome of the domesticus subspecies by swapping domesticus autosomes with their homologous partners from the musculus subspecies. We show that the "sterility" allelic combination of Prdm9 and Hstx2 can be activated by a musculus/domesticus heterozygosity of as few as two autosomes, Chromosome 17 (Chr 17) and Chr 18 and is further enhanced when another heterosubspecific autosomal pair is present, whereas it has no effect on meiotic progression in the pure domesticus genome. In addition, we identify a new X-linked hybrid sterility locus, Hstx3, at the centromeric end of Chr X, which modulates the incompatibility between Prdm9 and Hstx2. These results further support our concept of chromosomal hybrid sterility based on evolutionarily accumulated divergence between homologous sequences. Based on these and previous results, we believe that future studies should include more information on the mutual recognition of homologous chromosomes at or before the first meiotic prophase in interspecific hybrids, as this may serve as a general reproductive isolation checkpoint in mice and other species.

The Genetic Architecture of Recombination Rates is Polygenic and Differs Between the Sexes in Wild House Sparrows (Passer domesticus)

Meiotic recombination through chromosomal crossing-over is a fundamental feature of sex and an important driver of genomic diversity. It ensures proper disjunction, allows increased selection responses, and prevents mutation accumulation; however, it is also mutagenic and can break up favorable haplotypes. This cost-benefit dynamic is likely to vary depending on mechanistic and evolutionary contexts, and indeed, recombination rates show huge variation in nature. Identifying the genetic architecture of this variation is key to understanding its causes and consequences. Here, we investigate individual recombination rate variation in wild house sparrows (Passer domesticus). We integrate genomic and pedigree data to identify autosomal crossover counts (ACCs) and intrachromosomal allelic shuffling (r¯intra) in 13,056 gametes transmitted from 2,653 individuals to their offspring. Females had 1.37 times higher ACC, and 1.55 times higher r¯intra than males. ACC and r¯intra were heritable in females and males (ACC h2 = 0.23 and 0.11; r¯intra  h2 = 0.12 and 0.14), but cross-sex additive genetic correlations were low (rA = 0.29 and 0.32 for ACC and r¯intra). Conditional bivariate analyses showed that all measures remained heritable after accounting for genetic values in the opposite sex, indicating that sex-specific ACC and r¯intra can evolve somewhat independently. Genome-wide models showed that ACC and r¯intra are polygenic and driven by many small-effect loci, many of which are likely to act in trans as global recombination modifiers. Our findings show that recombination rates of females and males can have different evolutionary potential in wild birds, providing a compelling mechanism for the evolution of sexual dimorphism in recombination.

Novel crossover and recombination hotspots massively spread across primate genomes

BACKGROUND

The recombination landscape and subsequent natural selection have vast consequences forevolution and speciation. However, most of the crossover and recombination hotspots are yet to be discovered. We previously reported the relevance of C and G trinucleotide two-repeat units (CG-TTUs) in crossovers and recombination.

METHODS

On a genome-wide scale, here we mapped all combinations of A and T trinucleotide two-repeat units (AT-TTUs) in human, consisting of AATAAT, ATAATA, ATTATT, TTATTA, TATTAT, and TAATAA. We also compared a number of the colonies formed by the AT-TTUs (distance between consecutive AT-TTUs < 500 bp) in several other primates and mouse.

RESULTS

We found that the majority of the AT-TTUs (> 96%) resided in approximately 1.4 million colonies, spread throughout the human genome. In comparison to the CG-TTU colonies, the AT-TTU colonies were significantly more abundant and larger in size. Pure units and overlapping units of the pure units were readily detectable in the same colonies, signifying that the units were the sites of unequal crossover. We discovered dynamic sharedness of several of the colonies across the primate species studied, which mainly reached maximum complexity and size in human.

CONCLUSIONS

We report novel crossover and recombination hotspots of the finest molecular resolution, massively spread and shared across the genomes of human and several other primates. With respect to crossovers and recombination, these genomes are far more dynamic than previously envisioned.

The GC-content at the 5' ends of human protein-coding genes is undergoing mutational decay

BACKGROUND

In vertebrates, most protein-coding genes have a peak of GC-content near their 5' transcriptional start site (TSS). This feature promotes both the efficient nuclear export and translation of mRNAs. Despite the importance of GC-content for RNA metabolism, its general features, origin, and maintenance remain mysterious. We investigate the evolutionary forces shaping GC-content at the transcriptional start site (TSS) of genes through both comparative genomic analysis of nucleotide substitution rates between different species and by examining human de novo mutations.

RESULTS

Our data suggests that GC-peaks at TSSs were present in the last common ancestor of amniotes, and likely that of vertebrates. We observe that in apes and rodents, where recombination is directed away from TSSs by PRDM9, GC-content at the 5' end of protein-coding gene is currently undergoing mutational decay. In canids, which lack PRDM9 and perform recombination at TSSs, GC-content at the 5' end of protein-coding is increasing. We show that these patterns extend into the 5' end of the open reading frame, thus impacting synonymous codon position choices.

CONCLUSIONS

Our results indicate that the dynamics of this GC-peak in amniotes is largely shaped by historic patterns of recombination. Since decay of GC-content towards the mutation rate equilibrium is the default state for non-functional DNA, the observed decrease in GC-content at TSSs in apes and rodents indicates that the GC-peak is not being maintained by selection on most protein-coding genes in those species.

YY1 is involved in homologous recombination inhibition at guanine quadruplex sites in human cells

Homologous recombination (HR) is a key process for repairing DNA double strand breaks and for promoting genetic diversity. However, HR occurs unevenly across the genome, and certain genomic features can influence its activity. One such feature is the presence of guanine quadruplexes (G4s), stable secondary structures widely distributed throughout the genome. These G4s play essential roles in gene transcription and genome stability regulation. Especially, elevated G4 levels in cells deficient in the Bloom syndrome helicase (BLM) significantly enhance HR at G4 sites, potentially threatening genome stability. Here, we investigated the role of G4-binding protein Yin Yang-1 (YY1) in modulating HR at G4 sites in human cells. Our results show that YY1's binding to G4 structures suppresses sister chromatid exchange after BLM knockdown, and YY1's chromatin occupancy negatively correlates with the overall HR rate observed across the genome. By limiting RAD51 homolog 1 (RAD51) access, YY1 preferentially binds to essential genomic regions, shielding them from excessive HR. Our findings unveil a novel role of YY1-G4 interaction, revealing novel insights into cellular mechanisms involved in HR regulation.

Insights into non-crossover recombination from long-read sperm sequencing

Meiotic recombination is a fundamental process that generates genetic diversity by creating new combinations of existing alleles. Although human crossovers have been studied at the pedigree, population and single-cell level, the more frequent non-crossover events that lead to gene conversion are harder to study, particularly at the individual level. Here we show that single high-fidelity long sequencing reads from sperm can capture both crossovers and non-crossovers, allowing effectively arbitrary sample sizes for analysis from one male. Using fifteen sperm samples from thirteen donors we demonstrate variation between and within donors for the rates of different types of recombination. Intriguingly, we observe a tendency for non-crossover gene conversions to occur upstream of nearby PRDM9 binding sites, whereas crossover locations have a slight downstream bias. We further provide evidence for two distinct non-crossover processes. One gives rise to the vast majority of non-crossovers with mean conversion tract length under 50bp, which we suggest is an outcome of standard PRDM9-induced meiotic recombination. In contrast ~2% of non-crossovers have much longer mean tract length, and potentially originate from the same process as complex events with more than two haplotype switches, which is not associated with PRDM9 binding sites and is also seen in somatic cells.

Understanding the Genetic Basis of Variation in Meiotic Recombination: Past, Present, and Future

Meiotic recombination is a fundamental feature of sexually reproducing species. It is often required for proper chromosome segregation and plays important role in adaptation and the maintenance of genetic diversity. The molecular mechanisms of recombination are remarkably conserved across eukaryotes, yet meiotic genes and proteins show substantial variation in their sequence and function, even between closely related species. Furthermore, the rate and distribution of recombination shows a huge diversity within and between chromosomes, individuals, sexes, populations, and species. This variation has implications for many molecular and evolutionary processes, yet how and why this diversity has evolved is not well understood. A key step in understanding trait evolution is to determine its genetic basis-that is, the number, effect sizes, and distribution of loci underpinning variation. In this perspective, I discuss past and current knowledge on the genetic basis of variation in recombination rate and distribution, explore its evolutionary implications, and present open questions for future research.

Fine-scale contemporary recombination variation and its fitness consequences in adaptively diverging stickleback fish

Despite deep evolutionary conservation, recombination rates vary greatly across the genome and among individuals, sexes and populations. Yet the impact of this variation on adaptively diverging populations is not well understood. Here we characterized fine-scale recombination landscapes in an adaptively divergent pair of marine and freshwater populations of threespine stickleback from River Tyne, Scotland. Through whole-genome sequencing of large nuclear families, we identified the genomic locations of almost 50,000 crossovers and built recombination maps for marine, freshwater and hybrid individuals at a resolution of 3.8 kb. We used these maps to quantify the factors driving variation in recombination rates. We found strong heterochiasmy between sexes but also differences in recombination rates among ecotypes. Hybrids showed evidence of significant recombination suppression in overall map length and in individual loci. Recombination rates were lower not only within individual marine-freshwater-adaptive loci, but also between loci on the same chromosome, suggesting selection on linked gene 'cassettes'. Through temporal sampling along a natural hybrid zone, we found that recombinants showed traits associated with reduced fitness. Our results support predictions that divergence in cis-acting recombination modifiers, whose functions are disrupted in hybrids, may play an important role in maintaining differences among adaptively diverging populations.

Paleolithic Gene Duplications Primed Adaptive Evolution of Human Amylase Locus Upon Agriculture

Starch digestion is a cornerstone of human nutrition. The amylase genes code for the starch-digesting amylase enzyme. Previous studies suggested that the salivary amylase (AMY1) gene copy number increased in response to agricultural diets. However, the lack of nucleotide resolution of the amylase locus hindered detailed evolutionary analyses. Here, we have resolved this locus at nucleotide resolution in 98 present-day humans and identified 30 distinct haplotypes, revealing that the coding sequences of all amylase gene copies are evolving under negative selection. The phylogenetic reconstruction suggested that haplotypes with three AMY1 gene copies, prevalent across all continents and constituting about 70% of observed haplotypes, originated before the out-of-Africa migrations of ancestral modern humans. Using thousands of unique 25 base pair sequences across the amylase locus, we showed that additional AMY1 gene copies existed in the genomes of four archaic hominin genomes, indicating that the initial duplication of this locus may have occurred as far back 800,000 years ago. We similarly analyzed 73 ancient human genomes dating from 300 - 45,000 years ago and found that the AMY1 copy number variation observed today existed long before the advent of agriculture (~10,000 years ago), predisposing this locus to adaptive increase in the frequency of higher amylase copy number with the spread of agriculture. Mechanistically, the common three-copy haplotypes seeded non-allelic homologous recombination events that appear to be occurring at one of the fastest rates seen for tandem repeats in the human genome. Our study provides a comprehensive population-level understanding of the genomic structure of the amylase locus, identifying the mechanisms and evolutionary history underlying its duplication and copy number variability in relation to the onset of agriculture.

High prevalence of PRDM9-independent recombination hotspots in placental mammals

In many mammals, recombination events are concentrated in hotspots directed by a sequence-specific DNA-binding protein named PRDM9. Intriguingly, PRDM9 has been lost several times in vertebrates, and notably among mammals, it has been pseudogenized in the ancestor of canids. In the absence of PRDM9, recombination hotspots tend to occur in promoter-like features such as CpG islands. It has thus been proposed that one role of PRDM9 could be to direct recombination away from PRDM9-independent hotspots. However, the ability of PRDM9 to direct recombination hotspots has been assessed in only a handful of species, and a clear picture of how much recombination occurs outside of PRDM9-directed hotspots in mammals is still lacking. In this study, we derived an estimator of past recombination activity based on signatures of GC-biased gene conversion in substitution patterns. We quantified recombination activity in PRDM9-independent hotspots in 52 species of boreoeutherian mammals. We observe a wide range of recombination rates at these loci: several species (such as mice, humans, some felids, or cetaceans) show a deficit of recombination, while a majority of mammals display a clear peak of recombination. Our results demonstrate that PRDM9-directed and PRDM9-independent hotspots can coexist in mammals and that their coexistence appears to be the rule rather than the exception. Additionally, we show that the location of PRDM9-independent hotspots is relatively more stable than that of PRDM9-directed hotspots, but that PRDM9-independent hotspots nevertheless evolve slowly in concert with DNA hypomethylation.

Bridging the gap between the evolutionary dynamics and the molecular mechanisms of meiosis: A model based exploration of the PRDM9 intra-genomic Red Queen

Molecular dissection of meiotic recombination in mammals, combined with population-genetic and comparative studies, have revealed a complex evolutionary dynamic characterized by short-lived recombination hotspots. Hotspots are chromosome positions containing DNA sequences where the protein PRDM9 can bind and cause crossing-over. To explain these fast evolutionary dynamic, a so-called 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 hotspot extinction (the hotspot conversion paradox), followed by positive selection favoring mutant PRDM9 alleles recognizing new sequence motifs. Although this model predicts many empirical observations, the exact causes of the positive selection acting on new PRDM9 alleles is still not well understood. In this direction, experiment on mouse hybrids have suggested that, in addition to targeting double strand breaks, PRDM9 has another role during meiosis. Specifically, PRDM9 symmetric binding (simultaneous binding at the same site on both homologues) would facilitate homology search and, as a result, the pairing of the homologues. Although discovered in hybrids, this second function of PRDM9 could also be involved in the evolutionary dynamic observed within populations. To address this point, here, we present a theoretical model of the evolutionary dynamic of meiotic recombination integrating current knowledge about the molecular function of PRDM9. Our modeling work gives important insights into the selective forces driving the turnover of recombination hotspots. Specifically, the reduced symmetrical binding of PRDM9 caused by the loss of high affinity binding sites induces a net positive selection eliciting new PRDM9 alleles recognizing new targets. The model also offers new insights about the influence of the gene dosage of PRDM9, which can paradoxically result in negative selection on new PRDM9 alleles entering the population, driving their eviction and thus reducing standing variation at this locus.

Chromosome-level assembly of the gray fox (Urocyon cinereoargenteus) confirms the basal loss of PRDM9 in Canidae

Reference genome assemblies have been created from multiple lineages within the Canidae family; however, despite its phylogenetic relevance as a basal genus within the clade, there is currently no reference genome for the gray fox (Urocyon cinereoargenteus). Here, we present a chromosome-level assembly for the gray fox (U. cinereoargenteus), which represents the most contiguous, non-domestic canid reference genome available to date, with 90% of the genome contained in just 34 scaffolds and a contig N50 and scaffold N50 of 59.4 and 72.9 Megabases, respectively. Repeat analyses identified an increased number of simple repeats relative to other canids. Based on mitochondrial DNA, our Vermont sample clusters with other gray fox samples from the northeastern United States and contains slightly lower levels of heterozygosity than gray foxes on the west coast of California. This new assembly lays the groundwork for future studies to describe past and present population dynamics, including the delineation of evolutionarily significant units of management relevance. Importantly, the phylogenetic position of Urocyon allows us to verify the loss of PRDM9 functionality in the basal canid lineage, confirming that pseudogenization occurred at least 10 million years ago.

2023 ASHG Scientific Achievement Award

This article is based on the address given by the author at the 2023 meeting of The American Society of Human Genetics (ASHG) in Washington, D.C. A video of the original address can be found at the ASHG website.

Natural variation in the zinc-finger-encoding exon of Prdm9 affects hybrid sterility phenotypes in mice

PRDM9-mediated reproductive isolation was first described in the progeny of Mus musculus musculus (MUS) PWD/Ph and Mus musculus domesticus (DOM) C57BL/6J inbred strains. These male F1 hybrids fail to complete chromosome synapsis and arrest meiosis at prophase I, due to incompatibilities between the Prdm9 gene and hybrid sterility locus Hstx2. We identified 14 alleles of Prdm9 in exon 12, encoding the DNA-binding domain of the PRDM9 protein in outcrossed wild mouse populations from Europe, Asia, and the Middle East, 8 of which are novel. The same allele was found in all mice bearing introgressed t-haplotypes encompassing Prdm9. We asked whether 7 novel Prdm9 alleles in MUS populations and the t-haplotype allele in 1 MUS and 3 DOM populations induce Prdm9-mediated reproductive isolation. The results show that only combinations of the dom2 allele of DOM origin and the MUS msc1 allele ensure complete infertility of intersubspecific hybrids in outcrossed wild populations and inbred mouse strains examined so far. The results further indicate that MUS mice may share the erasure of PRDM9msc1 binding motifs in populations with different Prdm9 alleles, which implies that erased PRDM9 binding motifs may be uncoupled from their corresponding Prdm9 alleles at the population level. Our data corroborate the model of Prdm9-mediated hybrid sterility beyond inbred strains of mice and suggest that sterility alleles of Prdm9 may be rare.

The molecular machinery of meiotic recombination

Meiotic recombination, a cornerstone of eukaryotic diversity and individual genetic identity, is essential for the creation of physical linkages between homologous chromosomes, facilitating their faithful segregation during meiosis I. This process requires that germ cells generate controlled DNA lesions within their own genome that are subsequently repaired in a specialised manner. Repair of these DNA breaks involves the modulation of existing homologous recombination repair pathways to generate crossovers between homologous chromosomes. Decades of genetic and cytological studies have identified a multitude of factors that are involved in meiotic recombination. Recent work has started to provide additional mechanistic insights into how these factors interact with one another, with DNA, and provide the molecular outcomes required for a successful meiosis. Here, we provide a review of the recent developments with a focus on protein structures and protein-protein interactions.

Patterns of recombination in snakes reveal a tug-of-war between PRDM9 and promoter-like features

In some mammals, notably humans, recombination occurs almost exclusively where the protein PRDM9 binds, whereas in vertebrates lacking an intact PRDM9, such as birds and canids, recombination rates are elevated near promoter-like features. To determine whether PRDM9 directs recombination in nonmammalian vertebrates, we focused on an exemplar species with a single, intact PRDM9 ortholog, the corn snake (Pantherophis guttatus). Analyzing historical recombination rates along the genome and crossovers in pedigrees, we found evidence that PRDM9 specifies the location of recombination events, but we also detected a separable effect of promoter-like features. These findings reveal that the uses of PRDM9 and promoter-like features need not be mutually exclusive and instead reflect a tug-of-war that is more even in some species than others.

Recombination hotspots in Neurospora crassa controlled by idiomorphic sequences and meiotic silencing

Genes regulating recombination in specific chromosomal intervals of Neurospora crassa were described in the 1960s, but the mechanism is still unknown. For each of the rec-1, rec-2, and rec-3 genes, a single copy of the putative dominant allele, for example, rec-2SL found in St Lawrence OR74 A wild type, reduces recombination in chromosomal regions specific to that gene. However, when we sequenced the recessive allele, rec-2LG (derived from the Lindegren 1A wild type), we found that a 10 kb region in rec-2SL strains was replaced by a 2.7 kb unrelated sequence, making the "alleles" idiomorphs. When we introduced sad-1, a mutant lacking the RNA-dependent RNA polymerase that silences unpaired coding regions during meiosis into crosses heterozygous rec-2SL/rec-2LG, it increased recombination, indicating that meiotic silencing of a gene promoting recombination is responsible for dominant suppression of recombination. Consistent with this, mutation of rec-2LG by Repeat-Induced Point mutation generated an allele with multiple stop codons in the predicted rec-2 gene, which does not promote recombination and is recessive to rec-2LG. Sad-1 also relieves suppression of recombination in relevant target regions, in crosses heterozygous for rec-1 alleles and in crosses heterozygous for rec-3 alleles. We conclude that for all 3 known rec genes, 1 allele appears dominant only because meiotic silencing prevents the product of the active, "recessive," allele from stimulating recombination during meiosis. In addition, the proposed amino acid sequence of REC-2 suggests that regulation of recombination in Neurospora differs from any currently known mechanism.

Secondary Contact, Introgressive Hybridization, and Genome Stabilization in Sticklebacks

Advances in genomic studies have revealed that hybridization in nature is pervasive and raised questions about the dynamics of different genetic and evolutionary factors following the initial hybridization event. While recent research has proposed that the genomic outcomes of hybridization might be predictable to some extent, many uncertainties remain. With comprehensive whole-genome sequence data, we investigated the genetic introgression between 2 divergent lineages of 9-spined sticklebacks (Pungitius pungitius) in the Baltic Sea. We found that the intensity and direction of selection on the introgressed variation has varied across different genomic elements: while functionally important regions displayed reduced rates of introgression, promoter regions showed enrichment. Despite the general trend of negative selection, we identified specific genomic regions that were enriched for introgressed variants, and within these regions, we detected footprints of selection, indicating adaptive introgression. Geographically, we found the selection against the functional changes to be strongest in the vicinity of the secondary contact zone and weaken as a function of distance from the initial contact. Altogether, the results suggest that the stabilization of introgressed variation in the genomes is a complex, multistage process involving both negative and positive selection. In spite of the predominance of negative selection against introgressed variants, we also found evidence for adaptive introgression variants likely associated with adaptation to Baltic Sea environmental conditions.

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.

Metabolic stress-induced long ncRNA transcription governs the formation of meiotic DNA breaks in the fission yeast fbp1 gene

Meiotic recombination is a pivotal process that ensures faithful chromosome segregation and contributes to the generation of genetic diversity in offspring, which is initiated by the formation of double-strand breaks (DSBs). The distribution of meiotic DSBs is not uniform and is clustered at hotspots, which can be affected by environmental conditions. Here, we show that non-coding RNA (ncRNA) transcription creates meiotic DSBs through local chromatin remodeling in the fission yeast fbp1 gene. The fbp1 gene is activated upon glucose starvation stress, in which a cascade of ncRNA-transcription in the fbp1 upstream region converts the chromatin configuration into an open structure, leading to the subsequent binding of transcription factors. We examined the distribution of meiotic DSBs around the fbp1 upstream region in the presence and absence of glucose and observed several new DSBs after chromatin conversion under glucose starvation conditions. Moreover, these DSBs disappeared when cis-elements required for ncRNA transcription were mutated. These results indicate that ncRNA transcription creates meiotic DSBs in response to stress conditions in the fbp1 upstream region. This study addressed part of a long-standing unresolved mechanism underlying meiotic recombination plasticity in response to environmental fluctuation.

Fine-Scale Map Reveals Highly Variable Recombination Rates Associated with Genomic Features in the Eurasian Blackcap

Recombination is responsible for breaking up haplotypes, influencing genetic variability, and the efficacy of selection. Bird genomes lack the protein PR domain-containing protein 9, a key determinant of recombination dynamics in most metazoans. Historical recombination maps in birds show an apparent stasis in positioning recombination events. This highly conserved recombination pattern over long timescales may constrain the evolution of recombination in birds. At the same time, extensive variation in recombination rate is observed across the genome and between different species of birds. Here, we characterize the fine-scale historical recombination map of an iconic migratory songbird, the Eurasian blackcap (Sylvia atricapilla), using a linkage disequilibrium-based approach that accounts for population demography. Our results reveal variable recombination rates among and within chromosomes, which associate positively with nucleotide diversity and GC content and negatively with chromosome size. Recombination rates increased significantly at regulatory regions but not necessarily at gene bodies. CpG islands are associated strongly with recombination rates, though their specific position and local DNA methylation patterns likely influence this relationship. The association with retrotransposons varied according to specific family and location. Our results also provide evidence of heterogeneous intrachromosomal conservation of recombination maps between the blackcap and its closest sister taxon, the garden warbler. These findings highlight the considerable variability of recombination rates at different scales and the role of specific genomic features in shaping this variation. This study opens the possibility of further investigating the impact of recombination on specific population-genomic features.

Novel Insights into the Landscape of Crossover and Noncrossover Events in Rhesus Macaques (Macaca mulatta)

Meiotic recombination landscapes differ greatly between distantly and closely related taxa, populations, individuals, sexes, and even within genomes; however, the factors driving this variation are yet to be well elucidated. Here, we directly estimate contemporary crossover rates and, for the first time, noncrossover rates in rhesus macaques (Macaca mulatta) from four three-generation pedigrees comprising 32 individuals. We further compare these results with historical, demography-aware, linkage disequilibrium-based recombination rate estimates. From paternal meioses in the pedigrees, 165 crossover events with a median resolution of 22.3 kb were observed, corresponding to a male autosomal map length of 2,357 cM-approximately 15% longer than an existing linkage map based on human microsatellite loci. In addition, 85 noncrossover events with a mean tract length of 155 bp were identified-similar to the tract lengths observed in the only other two primates in which noncrossovers have been studied to date, humans and baboons. Consistent with observations in other placental mammals with PRDM9-directed recombination, crossover (and to a lesser extent noncrossover) events in rhesus macaques clustered in intergenic regions and toward the chromosomal ends in males-a pattern in broad agreement with the historical, sex-averaged recombination rate estimates-and evidence of GC-biased gene conversion was observed at noncrossover sites.

Meiotic DNA breaks drive multifaceted mutagenesis in the human germ line

Meiotic recombination commences with hundreds of programmed DNA breaks; however, the degree to which they are accurately repaired remains poorly understood. We report that meiotic break repair is eightfold more mutagenic for single-base substitutions than was previously understood, leading to de novo mutation in one in four sperm and one in 12 eggs. Its impact on indels and structural variants is even higher, with 100- to 1300-fold increases in rates per break. We uncovered new mutational signatures and footprints relative to break sites, which implicate unexpected biochemical processes and error-prone DNA repair mechanisms, including translesion synthesis and end joining in meiotic break repair. We provide evidence that these mechanisms drive mutagenesis in human germ lines and lead to disruption of hundreds of genes genome wide.

Characterization of meiotic recombination intermediates through gene knockouts in founder hybrid mice

Mammalian meiotic recombination proceeds via repair of hundreds of programmed DNA double-strand breaks, which requires choreographed binding of RPA, DMC1, and RAD51 to single-stranded DNA substrates. High-resolution in vivo binding maps of these proteins provide insights into the underlying molecular mechanisms. When assayed in F1-hybrid mice, these maps can distinguish the broken chromosome from the chromosome used as template for repair, revealing more mechanistic detail and enabling the structure of the recombination intermediates to be inferred. By applying CRISPR-Cas9 mutagenesis directly on F1-hybrid embryos, we have extended this approach to explore the molecular detail of recombination when a key component is knocked out. As a proof of concept, we have generated hybrid biallelic knockouts of Dmc1 and built maps of meiotic binding of RAD51 and RPA in them. DMC1 is essential for meiotic recombination, and comparison of these maps with those from wild-type mice is informative about the structure and timing of critical recombination intermediates. We observe redistribution of RAD51 binding and complete abrogation of D-loop recombination intermediates at a molecular level in Dmc1 mutants. These data provide insight on the configuration of RPA in D-loop intermediates and suggest that stable strand exchange proceeds via multiple rounds of strand invasion with template switching in mouse. Our methodology provides a high-throughput approach for characterization of gene function in meiotic recombination at low animal cost.

Take a walk on the KRAB side

Canonical Krüppel-associated box (KRAB)-containing zinc finger proteins (KZFPs) act as major repressors of transposable elements (TEs) via the KRAB-mediated recruitment of the heterochromatin scaffold KRAB-associated protein (KAP)1. KZFP genes emerged some 420 million years ago in the last common ancestor of coelacanth, lungfish, and tetrapods, and dramatically expanded to give rise to lineage-specific repertoires in contemporary species paralleling their TE load and turnover. However, the KRAB domain displays sequence and function variations that reveal repeated diversions from a linear TE-KZFP trajectory. This Review summarizes current knowledge on the evolution of KZFPs and discusses how ancestral noncanonical KZFPs endowed with variant KRAB, SCAN or DUF3669 domains have been utilized to achieve KAP1-independent functions.

Variability in the PRDM9 gene in Sindhi cattle

BACKGROUND

Sindhi is a dual-purpose breed adapted to tropical environments. However, this breed has the smallest total population among indicine breeds in Brazil and the smallest effective number. In addition, the inbreeding coefficient is higher than 6.25% in ~ 60% of the population. Therefore, alternatives to increase genetic diversity are important. Within this context, the PRDM9 gene is particularly interesting since it is involved in meiotic recombination events, consequently enhancing genetic variability in the population by increasing the number of circulating haplotypes. Each allele of the gene induces recombination at a different hotspot. The larger the number of circulating alleles, the higher the recombination rate and the greater the genetic variability.

METHODS

The aim of this study was to characterize alleles of the PRDM9 gene in Sindhi cattle. The region of the zinc finger domains of the gene was amplified by PCR, genotyped, and sequenced for allele identification in 50 Sindhi animals.

RESULTS

Three alleles (A-cattle1, B-cattle14, and C-cattle19) and six genotypes (AA, BB, CC, AB, AC, and BC) were identified.

CONCLUSION

The allele variation of the PRDM9 gene in the Sindhi breed enables to guide the mating of animals with different genotypes/alleles and to promote genetic variability by recombination if there is intralocus variability.