Extensive sex differences at the initiation of genetic recombination

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.

Principles of chromosome organization for meiotic recombination

In meiotic cells, chromosomes are organized as chromatin loop arrays anchored to a protein axis. This organization is essential to regulate meiotic recombination, from DNA double-strand break (DSB) formation to their repair. In mammals, it is unknown how chromatin loops are organized along the genome and how proteins participating in DSB formation are tethered to the chromosome axes. Here, we identify three categories of axis-associated genomic sites: PRDM9 binding sites, where DSBs form; binding sites of the insulator protein CTCF; and H3K4me3-enriched sites. We demonstrate that PRDM9 promotes the recruitment of MEI4 and IHO1, two proteins essential for DSB formation. In turn, IHO1 anchors DSB sites to the axis components HORMAD1 and SYCP3. We discovered that IHO1, HORMAD1, and SYCP3 are associated at the DSB ends during DSB repair. Our results highlight how interactions of proteins with specific genomic elements shape the meiotic chromosome organization for recombination.

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.

Mapping crossover events of mouse meiotic recombination by restriction fragment ligation-based Refresh-seq

Single-cell whole-genome sequencing methods have undergone great improvements over the past decade. However, allele dropout, which means the inability to detect both alleles simultaneously in an individual diploid cell, largely restricts the application of these methods particularly for medical applications. Here, we develop a new single-cell whole-genome sequencing method based on third-generation sequencing (TGS) platform named Refresh-seq (restriction fragment ligation-based genome amplification and TGS). It is based on restriction endonuclease cutting and ligation strategy in which two alleles in an individual cell can be cut into equal fragments and tend to be amplified simultaneously. As a new single-cell long-read genome sequencing method, Refresh-seq features much lower allele dropout rate compared with SMOOTH-seq. Furthermore, we apply Refresh-seq to 688 sperm cells and 272 female haploid cells (secondary polar bodies and parthenogenetic oocytes) from F1 hybrid mice. We acquire high-resolution genetic map of mouse meiosis recombination at low sequencing depth and reveal the sexual dimorphism in meiotic crossovers. We also phase the structure variations (deletions and insertions) in sperm cells and female haploid cells with high precision. Refresh-seq shows great performance in screening aneuploid sperm cells and oocytes due to the low allele dropout rate and has great potential for medical applications such as preimplantation genetic diagnosis.

3D genome remodeling and homologous pairing during meiotic prophase of mouse oogenesis and spermatogenesis

During meiosis, the chromatin and transcriptome undergo prominent switches. Although recent studies have explored the genome reorganization during spermatogenesis, the chromatin remodeling in oogenesis and characteristics of homologous pairing remain largely elusive. We comprehensively compared chromatin structures and transcriptomes at successive substages of meiotic prophase in both female and male mice using low-input high-through chromosome conformation capture (Hi-C) and RNA sequencing (RNA-seq). Compartments and topologically associating domains (TADs) gradually disappeared and slowly recovered in both sexes. We found that homologs adopted different sex-conserved pairing strategies prior to and after the leptotene-to-zygotene transition, changing from long interspersed nuclear element (LINE)-enriched compartments B to short interspersed nuclear element (SINE)-enriched compartments A. We complemented marker genes and predicted the sex-specific meiotic sterile genes for each substage. This study provides valuable insights into the similarities and distinctions between sexes in chromosome architecture, homologous pairing, and transcriptome during meiotic prophase of both oogenesis and spermatogenesis.

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.

Genomic and demographic processes differentially influence genetic variation across the human X chromosome

Many forces influence genetic variation across the genome including mutation, recombination, selection, and demography. Increased mutation and recombination both lead to increases in genetic diversity in a region-specific manner, while complex demographic patterns shape patterns of diversity on a more global scale. While these processes act across the entire genome, the X chromosome is particularly interesting because it contains several distinct regions that are subject to different combinations and strengths of these forces: the pseudoautosomal regions (PARs) and the X-transposed region (XTR). The X chromosome thus can serve as a unique model for studying how genetic and demographic forces act in different contexts to shape patterns of observed variation. We therefore sought to explore diversity, divergence, and linkage disequilibrium in each region of the X chromosome using genomic data from 26 human populations. Across populations, we find that both diversity and substitution rate are consistently elevated in PAR1 and the XTR compared to the rest of the X chromosome. In contrast, linkage disequilibrium is lowest in PAR1, consistent with the high recombination rate in this region, and highest in the region of the X chromosome that does not recombine in males. However, linkage disequilibrium in the XTR is intermediate between PAR1 and the autosomes, and much lower than the non-recombining X. Finally, in addition to these global patterns, we also observed variation in ratios of X versus autosomal diversity consistent with population-specific evolutionary history as well. While our results were generally consistent with previous work, two unexpected observations emerged. First, our results suggest that the XTR does not behave like the rest of the recombining X and may need to be evaluated separately in future studies. Second, the different regions of the X chromosome appear to exhibit unique patterns of linked selection across different human populations. Together, our results highlight profound regional differences across the X chromosome, simultaneously making it an ideal system for exploring the action of evolutionary forces as well as necessitating its careful consideration and treatment in genomic analyses.

Comprehensive dissection of meiotic DNA double-strand breaks and crossovers in cucumber

Meiotic recombination drives genetic diversity and crop genome optimization. In plant breeding, parents with favorable traits are crossed to create elite varieties. Different hybridizations produce diverse types of segment reshuffling between homologous chromosomes. However, little is known about the factors that cause hybrid-specific changes in crossovers (COs). Here, we constructed 2 F2 populations from crosses between a semiwild and 2 domesticated cucumber (Cucumis sativus) accessions and examined CO events. COs mainly occurred around genes and differed unevenly along chromosomes between the 2 hybrids. Fine-scale CO distributions were suppressed in regions of heterozygous structural variations (SVs) and were accelerated by high sequence polymorphism. C. sativus RADiation sensitive 51A (CsRAD51A) binding, histone H3 lysine 4 trimethylation (H3K4me3) modification, chromatin accessibility, and hypomethylation were positively associated with global CO landscapes and in local DNA double-strand break (DSB) hotspots and genes. The frequency and suppression of COs could be roughly predicted based on multiomic information. Differences in CO events between hybrids could be partially traced to distinct genetic and epigenetic features and were significantly associated with specific DSB hotspots and heterozygous SVs. Our findings identify the genomic and epigenetic features that contribute to CO formation and hybrid-specific divergence in cucumber and provide theoretical support for selecting parental combinations and manipulating recombination events at target genomic regions during plant breeding.

Sexual dimorphic regulation of recombination by the synaptonemal complex in C. elegans

In sexually reproducing organisms, germ cells faithfully transmit the genome to the next generation by forming haploid gametes, such as eggs and sperm. Although most meiotic proteins are conserved between eggs and sperm, many aspects of meiosis are sexually dimorphic, including the regulation of recombination. The synaptonemal complex (SC), a large ladder-like structure that forms between homologous chromosomes, is essential for regulating meiotic chromosome organization and promoting recombination. To assess whether sex-specific differences in the SC underpin sexually dimorphic aspects of meiosis, we examined Caenorhabditis elegans SC central region proteins (known as SYP proteins) in oogenesis and spermatogenesis and uncovered sex-specific roles for the SYPs in regulating meiotic recombination. We find that SC composition, specifically SYP-2, SYP-3, SYP-5, and SYP-6, is regulated by sex-specific mechanisms throughout meiotic prophase I. During pachytene, both oocytes and spermatocytes differentially regulate the stability of SYP-2 and SYP-3 within an assembled SC. Further, we uncover that the relative amount of SYP-2 and SYP-3 within the SC is independently regulated in both a sex-specific and a recombination-dependent manner. Specifically, we find that SYP-2 regulates the early steps of recombination in both sexes, while SYP-3 controls the timing and positioning of crossover recombination events across the genomic landscape in only oocytes. Finally, we find that SYP-2 and SYP-3 dosage can influence the composition of the other SYPs in the SC via sex-specific mechanisms during pachytene. Taken together, we demonstrate dosage-dependent regulation of individual SC components with sex-specific functions in recombination. These sexual dimorphic features of the SC provide insights into how spermatogenesis and oogenesis adapted similar chromosome structures to differentially regulate and execute recombination.

Modeling Recombination Rate as a Quantitative Trait Reveals New Insight into Selection in Humans

Meiotic recombination is both a fundamental biological process required for proper chromosomal segregation during meiosis and an important genomic parameter that shapes major features of the genomic landscape. However, despite the central importance of this phenotype, we lack a clear understanding of the selective pressures that shape its variation in natural populations, including humans. While there is strong evidence of fitness costs of low rates of recombination, the possible fitness costs of high rates of recombination are less defined. To determine whether a single lower fitness bound can explain the variation in recombination rates observed in human populations, we simulated the evolution of recombination rates as a sexually dimorphic quantitative trait. Under each scenario, we statistically compared the resulting trait distribution with the observed distribution of recombination rates from a published study of the Icelandic population. To capture the genetic architecture of recombination rates in humans, we modeled it as a moderately complex trait with modest heritability. For our fitness function, we implemented a hyperbolic tangent curve with several flexible parameters to capture a wide range of existing hypotheses. We found that costs of low rates of recombination alone are likely insufficient to explain the current variation in recombination rates in both males and females, supporting the existence of fitness costs of high rates of recombination in humans. With simulations using both upper and lower fitness boundaries, we describe a parameter space for the costs of high recombination rates that produces results consistent with empirical observations.

DEtail-seq is an ultra-efficient and convenient method for meiotic DNA break profiling in multiple organisms

Programmed DNA double-strand break (DSB) formation is a crucial step in meiotic recombination, yet techniques for high-efficiency and precise mapping of the 3' ends of DSBs are still in their infancy. Here, we report a novel technique, named DNA End tailing and sequencing (DEtail-seq), which can directly and ultra-efficiently characterize the 3' ends of meiotic DSBs with near single-nucleotide resolution in a variety of species, including yeast, mouse, and human. We find that the 3' ends of meiotic DSBs are stable without significant resection in budding yeast. Meiotic DSBs are strongly enriched in de novo H3K4me3 peaks in the mouse genome at leptotene stage. We also profile meiotic DSBs in human and find DSB hotspots are enriched near the common fragile sites during human meiosis, especially at CCCTC-binding factor (CTCF)-associated enhancers. Therefore, DEtail-seq provides a powerful method to detect DSB ends in various species, and our results provide new insights into the distribution and regulation of meiotic DSB hotspots.

Divergent patterns of meiotic double strand breaks and synapsis initiation dynamics suggest an evolutionary shift in the meiosis program between American and Australian marsupials

In eutherian mammals, hundreds of programmed DNA double-strand breaks (DSBs) are generated at the onset of meiosis. The DNA damage response is then triggered. Although the dynamics of this response is well studied in eutherian mammals, recent findings have revealed different patterns of DNA damage signaling and repair in marsupial mammals. To better characterize these differences, here we analyzed synapsis and the chromosomal distribution of meiotic DSBs markers in three different marsupial species (Thylamys elegans, Dromiciops gliorides, and Macropus eugenii) that represent South American and Australian Orders. Our results revealed inter-specific differences in the chromosomal distribution of DNA damage and repair proteins, which were associated with differing synapsis patterns. In the American species T. elegans and D. gliroides, chromosomal ends were conspicuously polarized in a bouquet configuration and synapsis progressed exclusively from the telomeres towards interstitial regions. This was accompanied by sparse H2AX phosphorylation, mainly accumulating at chromosomal ends. Accordingly, RAD51 and RPA were mainly localized at chromosomal ends throughout prophase I in both American marsupials, likely resulting in reduced recombination rates at interstitial positions. In sharp contrast, synapsis initiated at both interstitial and distal chromosomal regions in the Australian representative M. eugenii, the bouquet polarization was incomplete and ephemeral, γH2AX had a broad nuclear distribution, and RAD51 and RPA foci displayed an even chromosomal distribution. Given the basal evolutionary position of T. elegans, it is likely that the meiotic features reported in this species represent an ancestral pattern in marsupials and that a shift in the meiotic program occurred after the split of D. gliroides and the Australian marsupial clade. Our results open intriguing questions about the regulation and homeostasis of meiotic DSBs in marsupials. The low recombination rates observed at the interstitial chromosomal regions in American marsupials can result in the formation of large linkage groups, thus having an impact in the evolution of their genomes.

Epitope tag-specific differences in the detection of COSA-1 marked crossover sites in C. elegans spermatocytes

Nascent crossover sites in C. elegans meiocytes can be cytologically detected using epitope-tagged versions of the pro-crossover protein COSA-1. In spermatocytes, differences exist between cytologically-detected and genetically-detected double crossover rates. Here, we examine nascent crossovers using both GFP- and OLLAS-tagged COSA-1. Similar to previous work, we find that most late pachytene spermatocytes display 5 COSA-1 foci, indicating one crossover per autosome bivalent. However, we detected more nuclei with >5 COSA-1 foci using OLLAS::COSA-1, reflecting some bivalents having 2 COSA-1 foci. These results demonstrate tag-specific differences in the detection of COSA-1 marked nascent crossovers in spermatocytes.

Rapid Evolution of the Fine-scale Recombination Landscape in Wild House Mouse (Mus musculus) Populations

Meiotic recombination is an important evolutionary force and an essential meiotic process. In many species, recombination events concentrate into hotspots defined by the site-specific binding of PRMD9. Rapid evolution of Prdm9's zinc finger DNA-binding array leads to remarkably abrupt shifts in the genomic distribution of hotspots between species, but the question of how Prdm9 allelic variation shapes the landscape of recombination between populations remains less well understood. Wild house mice (Mus musculus) harbor exceptional Prdm9 diversity, with >150 alleles identified to date, and pose a particularly powerful system for addressing this open question. We employed a coalescent-based approach to construct broad- and fine-scale sex-averaged recombination maps from contemporary patterns of linkage disequilibrium in nine geographically isolated wild house mouse populations, including multiple populations from each of three subspecies. Comparing maps between wild mouse populations and subspecies reveals several themes. First, we report weak fine- and broad-scale recombination map conservation across subspecies and populations, with genetic divergence offering no clear prediction for recombination map divergence. Second, most hotspots are unique to one population, an outcome consistent with minimal sharing of Prdm9 alleles between surveyed populations. Finally, by contrasting aggregate hotspot activity on the X versus autosomes, we uncover evidence for population-specific differences in the degree and direction of sex dimorphism for recombination. Overall, our findings illuminate the variability of both the broad- and fine-scale recombination landscape in M. musculus and underscore the functional impact of Prdm9 allelic variation in wild mouse populations.

TOPOVIBL-REC114 interaction regulates meiotic DNA double-strand breaks

Meiosis requires the formation of programmed DNA double strand breaks (DSBs), essential for fertility and for generating genetic diversity. DSBs are induced by the catalytic activity of the TOPOVIL complex formed by SPO11 and TOPOVIBL. To ensure genomic integrity, DNA cleavage activity is tightly regulated, and several accessory factors (REC114, MEI4, IHO1, and MEI1) are needed for DSB formation in mice. How and when these proteins act is not understood. Here, we show that REC114 is a direct partner of TOPOVIBL, and identify their conserved interacting domains by structural analysis. We then analyse the role of this interaction by monitoring meiotic DSBs in female and male mice carrying point mutations in TOPOVIBL that decrease or disrupt its binding to REC114. In these mutants, DSB activity is strongly reduced genome-wide in oocytes, and only in sub-telomeric regions in spermatocytes. In addition, in mutant spermatocytes, DSB activity is delayed in autosomes. These results suggest that REC114 is a key member of the TOPOVIL catalytic complex, and that the REC114/TOPOVIBL interaction ensures the efficiency and timing of DSB activity.

High-resolution mapping reveals hotspots and sex-biased recombination in Populus trichocarpa

Fine-scale meiotic recombination is fundamental to the outcome of natural and artificial selection. Here, dense genetic mapping and haplotype reconstruction were used to estimate recombination for a full factorial Populus trichocarpa cross of 7 males and 7 females. Genomes of the resulting 49 full-sib families (N = 829 offspring) were resequenced, and high-fidelity biallelic SNP/INDELs and pedigree information were used to ascertain allelic phase and impute progeny genotypes to recover gametic haplotypes. The 14 parental genetic maps contained 1,820 SNP/INDELs on average that covered 376.7 Mb of physical length across 19 chromosomes. Comparison of parental and progeny haplotypes allowed fine-scale demarcation of cross-over regions, where 38,846 cross-over events in 1,658 gametes were observed. Cross-over events were positively associated with gene density and negatively associated with GC content and long-terminal repeats. One of the most striking findings was higher rates of cross-overs in males in 8 out of 19 chromosomes. Regions with elevated male cross-over rates had lower gene density and GC content than windows showing no sex bias. High-resolution analysis identified 67 candidate cross-over hotspots spread throughout the genome. DNA sequence motifs enriched in these regions showed striking similarity to those of maize, Arabidopsis, and wheat. These findings, and recombination estimates, will be useful for ongoing efforts to accelerate domestication of this and other biomass feedstocks, as well as future studies investigating broader questions related to evolutionary history, perennial development, phenology, wood formation, vegetative propagation, and dioecy that cannot be studied using annual plant model systems.

Unwinding during stressful times: Mechanisms of helicases in meiotic recombination

Successful meiosis I requires that homologous chromosomes be correctly linked before they are segregated. In most organisms this physical linkage is achieved through the generation of crossovers between the homologs. Meiotic recombination co-opts and modifies the canonical homologous recombination pathway to successfully generate crossovers One of the central components of this pathway are a number of conserved DNA helicases. Helicases couple nucleic acid binding to nucleotide hydrolysis and use this activity to modify DNA or protein-DNA substrates. During meiosis I it is necessary for the cell to modulate the canonical DNA repair pathways in order to facilitate the generation of interhomolog crossovers. Many of these meiotic modulations take place in pathways involving DNA helicases, or with a meiosis specific helicase. This short review explores what is currently understood about these helicases, their interaction partners, and the role of regulatory modifications during meiosis I. We focus in particular on the molecular structure and mechanisms of these helicases.

Analysis of the Genomic Sequence of ABO Allele Using Next-Generation Sequencing Method

Background

Although many molecular diagnostic methods have been used for ABO genotyping, there are few reports on the full-length genomic sequence analysis of the ABO gene. Recently, next-generation sequencing (NGS) has been shown to provide fast and high-throughput results and is widely used in the clinical laboratory. Here, we established an NGS method for analyzing the sequence of the start codon to the stop codon in the ABO gene.

Study Design and Methods

Two pairs of primers covering the partial 5’-untranslated region (UTR) to 3’-UTR of the ABO gene were designed. The sequences covering from the start codon to the stop codon of the ABO gene were amplified using these primers, and an NGS method based on the overlap amplicon was developed. A total of 110 individuals, including 88 blood donors with normal phenotypes and 22 ABO subtypes, were recruited and analyzed. All these specimens were first detected by serological tests and then determined by polymerase chain reaction sequence-based typing (PCR-SBT) and NGS. The sequences, including all the intron regions for the specimens, were analyzed by bioinformatics software.

Results

Among the 88 blood donors with a normal phenotype, 48 homozygous individuals, 39 heterozygous individuals, and one individual with a novel O allele were found according to the results of the PCR-SBT method. Some single-nucleotide variants (SNV) in intronic regions were found to be specific for different ABO alleles from 48 homozygous individuals using the NGS method. Sequences in the coding region of all specimens using the NGS method were the same as those of the PCR-SBT method. Three intronic SNVs were found to be associated with the ABO subtypes, including one novel intronic SNV (c.28+5956T>A). Moreover, six specimens were found to exhibit DNA recombination.

Conclusion

An NGS method was established to analyze the sequence from the start codon to the stop codon of the ABO gene. Two novel ABO alleles were identified, and DNA recombination was found to exist in the ABO alleles.

Adaptive Control of the Meiotic Recombination Landscape by DNA Site-dependent Hotspots With Implications for Evolution

Meiosis is an essential component of the sexual life cycle in eukaryotes. The independent assortment of chromosomes in meiosis increases genetic diversity at the level of whole chromosomes and meiotic recombination increases genetic diversity within chromosomes. The resulting variability fuels evolution. Interestingly, global mapping of recombination in diverse taxa revealed dramatic changes in its frequency distribution between closely related species, subspecies, and even isolated populations of the same species. New insight into mechanisms for these evolutionarily rapid changes has come from analyses of environmentally induced plasticity of recombination in fission yeast. Many different DNA sites, and where identified their binding/activator proteins, control the positioning of recombination at hotspots. Each different class of hotspots functions as an independently controlled rheostat that modulates rates of recombination over a broad dynamic range in response to changing conditions. Together, this independent modulation can rapidly and dramatically alter the global frequency distribution of recombination. This process likely contributes substantially to (i.e., can largely explain) evolutionarily rapid, Prdm9-independent changes in the recombination landscape. Moreover, the precise control mechanisms allow cells to dynamically favor or disfavor newly arising combinations of linked alleles in response to changing extracellular and intracellular conditions, which has striking implications for the impacts of meiotic recombination on evolution.

Deep learning identifies and quantifies recombination hotspot determinants

MOTIVATION

Recombination is one of the essential genetic processes for sexually reproducing organisms, which can happen more frequently in some regions, called recombination hotspots. Although several factors, such as PRDM9 binding motifs, are known to be related to the hotspots, their contributions to the recombination hotspots have not been quantified, and other determinants are yet to be elucidated. Here, we propose a computational method, RHSNet, based on deep learning and signal processing, to identify and quantify the hotspot determinants in a purely data-driven manner, utilizing datasets from various studies, populations, sexes and species.

RESULTS

RHSNet can significantly outperform other sequence-based methods on multiple datasets across different species, sexes and studies. In addition to being able to identify hotspot regions and the well-known determinants accurately, more importantly, RHSNet can quantify the determinants that contribute significantly to the recombination hotspot formation in the relation between PRDM9 binding motif, histone modification and GC content. Further cross-sex, cross-population and cross-species studies suggest that the proposed method has the generalization power and potential to identify and quantify the evolutionary determinant motifs.

AVAILABILITY AND IMPLEMENTATION

https://github.com/frankchen121212/RHSNet.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

3D chromatin remodelling in the germ line modulates genome evolutionary plasticity

Chromosome folding has profound impacts on gene regulation, whose evolutionary consequences are far from being understood. Here we explore the relationship between 3D chromatin remodelling in mouse germ cells and evolutionary changes in genome structure. Using a comprehensive integrative computational analysis, we (i) reconstruct seven ancestral rodent genomes analysing whole-genome sequences of 14 species representatives of the major phylogroups, (ii) detect lineage-specific chromosome rearrangements and (iii) identify the dynamics of the structural and epigenetic properties of evolutionary breakpoint regions (EBRs) throughout mouse spermatogenesis. Our results show that EBRs are devoid of programmed meiotic DNA double-strand breaks (DSBs) and meiotic cohesins in primary spermatocytes, but are associated in post-meiotic cells with sites of DNA damage and functional long-range interaction regions that recapitulate ancestral chromosomal configurations. Overall, we propose a model that integrates evolutionary genome reshuffling with DNA damage response mechanisms and the dynamic spatial genome organisation of germ cells.

The role of genome folding in the heritability and evolvability of structural variations is not well understood. Here the authors investigate the impact of the three-dimensional genome topology of germ cells in the formation and transmission of gross structural genomic changes detected from comparing whole-genome sequences of 14 rodent species.

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

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

Differential effects of two catalytic mutations on full-length PRDM9 and its isolated PR/SET domain reveal a case of pseudomodularity

PRDM9 is a DNA-binding histone methyltransferase that designates and activates recombination hotspots in mammals by locally trimethylating lysines 4 and 36 of histone H3. In mice, we recently reported two independently produced point mutations at the same residue, Glu360Pro (Prdm9EP) and Glu360Lys (Prdm9EK), which severely reduce its H3K4 and H3K36 methyltransferase activities in vivo. Prdm9EP is slightly less hypomorphic than Prdm9EK, but both mutations reduce both the number and amplitude of PRDM9-dependent H3K4me3 and H3K36me3 peaks in spermatocytes. While both mutations cause infertility with complete meiotic arrest in males, Prdm9EP, but not Prdm9EK, is compatible with some female fertility. When we tested the effects of these mutations in vitro, both Prdm9EP and Prdm9EK abolished H3K4 and H3K36 methyltransferase activity in full-length PRDM9. However, in the isolated PRDM9 PR/SET domain, these mutations selectively compromised H3K36 methyltransferase activity, while leaving H3K4 methyltransferase activity intact. The difference in these effects on the PR/SET domain vs the full-length protein shows that PRDM9 is not an intrinsically modular enzyme; its catalytic domain is influenced by its tertiary structure and possibly by its interactions with DNA and other proteins in vivo. These two informative mutations illuminate the enzymatic chemistry of PRDM9, and potentially of PR/SET domains in general, reveal the minimal threshold of PRDM9-dependent catalytic activity for female fertility, and potentially have some practical utility for genetic mapping and genomics.

Constraint and divergence in the evolution of male and female recombination rates in fishes

Recombination is a fundamental feature of sexual reproduction across eukaryotes, yet recombination rates are highly variable both within and between species. In particular, sex differences in recombination rate between males and females (heterochiasmy) are more often the rule than the exception, but despite the prevalence of heterochiasmy the ultimate causes of global patterns of heterochiasmy remain unclear. Here, we assemble a comprehensive dataset of sex-specific recombination rate estimates for 61 fish species, and combine this with information on sex determination, fertilization mode, and sexual dimorphism to test competing theories for the causes and evolution of heterochiasmy. We find that sex differences in recombination rate are evolutionary labile, with frequent shifts in the direction and magnitude of heterochiasmy. This rapid turnover does not appear to be driven by simple neutral processes and is inconsistent with nonadaptive explanations for heterochiasmy, including biological sex differences in meiosis. Although patterns of heterochiasmy across the phylogeny indicate a potential role for adaptive processes, we are unable to directly link variation in heterochiasmy with proxies for sexual selection or sexual conflict across species, indicating that these effects-if present-are either subtle or complex. Finally, we show evidence for correlated rates of recombination rate evolution between males and females, indicating the potential for genetic constraints and sexual conflict over the recombination landscape.

Transcriptional progression during meiotic prophase I reveals sex-specific features and X chromosome dynamics in human fetal female germline

During gametogenesis in mammals, meiosis ensures the production of haploid gametes. The timing and length of meiosis to produce female and male gametes differ considerably. In contrast to males, meiotic prophase I in females initiates during development. Hence, the knowledge regarding progression through meiotic prophase I is mainly focused on human male spermatogenesis and female oocyte maturation during adulthood. Therefore, it remains unclear how the different stages of meiotic prophase I between human oogenesis and spermatogenesis compare. Analysis of single-cell transcriptomics data from human fetal germ cells (FGC) allowed us to identify the molecular signatures of female meiotic prophase I stages leptotene, zygotene, pachytene and diplotene. We have compared those between male and female germ cells in similar stages of meiotic prophase I and revealed conserved and specific features between sexes. We identified not only key players involved in the process of meiosis, but also highlighted the molecular components that could be responsible for changes in cellular morphology that occur during this developmental period, when the female FGC acquire their typical (sex-specific) oocyte shape as well as sex-differences in the regulation of DNA methylation. Analysis of X-linked expression between sexes during meiotic prophase I suggested a transient X-linked enrichment during female pachytene, that contrasts with the meiotic sex chromosome inactivation in males. Our study of the events that take place during meiotic prophase I provide a better understanding not only of female meiosis during development, but also highlights biomarkers that can be used to study infertility and offers insights in germline sex dimorphism in humans.

A pathogenic DMC1 frameshift mutation causes nonobstructive azoospermia but not primary ovarian insufficiency in humans

Nonobstructive azoospermia (NOA) and diminished ovarian reserve (DOR) are two disorders that can lead to infertility in males and females. Genetic factors have been identified to contribute to NOA and DOR. However, the same genetic factor that can cause both NOA and DOR remains largely unknown. To explore the candidate pathogenic gene that causes both NOA and DOR, we conducted whole-exome sequencing (WES) in a non-consanguineous family with two daughters with DOR and a son with NOA. We detected one pathogenic frameshift variant (NM_007068:c.28delG, p. Glu10Asnfs*31) following a recessive inheritance mode in a meiosis gene DMC1 (DNA meiotic recombinase 1). Clinical analysis showed reduced antral follicle number in both daughters with DOR, but metaphase II oocytes could be retrieved from one of them. For the son with NOA, no spermatozoa were found after microsurgical testicular sperm extraction. A further homozygous Dmc1 knockout mice study demonstrated total failure of follicle development and spermatogenesis. These results revealed a discrepancy of DMC1 action between mice and humans. In humans, DMC1 is required for spermatogenesis but is dispensable for oogenesis, although the loss of function of this gene may lead to DOR. To our knowledge, this is the first report on the homozygous frameshift mutation as causative for both NOA and DOR and demonstrating that DMC1 is dispensable in human oogenesis.

Applications of Single-Cell DNA Sequencing

Over the past decade, genomic analyses of single cells-the fundamental units of life-have become possible. Single-cell DNA sequencing has shed light on biological questions that were previously inaccessible across diverse fields of research, including somatic mutagenesis, organismal development, genome function, and microbiology. Single-cell DNA sequencing also promises significant future biomedical and clinical impact, spanning oncology, fertility, and beyond. While single-cell approaches that profile RNA and protein have greatly expanded our understanding of cellular diversity, many fundamental questions in biology and important biomedical applications require analysis of the DNA of single cells. Here, we review the applications and biological questions for which single-cell DNA sequencing is uniquely suited or required. We include a discussion of the fields that will be impacted by single-cell DNA sequencing as the technology continues to advance.

Chromosome-wide characterization of meiotic noncrossovers (gene conversions) in mouse hybrids

During meiosis, the recombination-initiating DNA double-strand breaks (DSBs) are repaired by crossovers or noncrossovers (gene conversions). While crossovers are easily detectable, noncrossover identification is hampered by the small size of their converted tracts and the necessity of sequence polymorphism. We report identification and characterization of a mouse chromosome-wide set of noncrossovers by next-generation sequencing of 10 mouse intersubspecific chromosome substitution strains. Based on 94 identified noncrossovers, we determined the mean length of a conversion tract to be 32 bp. The spatial chromosome-wide distribution of noncrossovers and crossovers significantly differed, although both sets overlapped the known hotspots of PRDM9-directed histone methylation and DNA DSBs, thus supporting their origin in the standard DSB repair pathway. A significant deficit of noncrossovers descending from asymmetric DSBs proved their proposed adverse effect on meiotic recombination and pointed to sister chromatids as an alternative template for their repair. The finding has implications for the molecular mechanism of hybrid sterility in mice from crosses between closely related Mus musculus musculus and Mus musculus domesticus subspecies.

Meiotic recombination mirrors patterns of germline replication in mice and humans

Genetic recombination generates novel trait combinations, and understanding how recombination is distributed across the genome is key to modern genetics. The PRDM9 protein defines recombination hotspots; however, megabase-scale recombination patterning is independent of PRDM9. The single round of DNA replication, which precedes recombination in meiosis, may establish these patterns; therefore, we devised an approach to study meiotic replication that includes robust and sensitive mapping of replication origins. We find that meiotic DNA replication is distinct; reduced origin firing slows replication in meiosis, and a distinctive replication pattern in human males underlies the subtelomeric increase in recombination. We detected a robust correlation between replication and both contemporary and historical recombination and found that replication origin density coupled with chromosome size determines the recombination potential of individual chromosomes. Our findings and methods have implications for understanding the mechanisms underlying DNA replication, genetic recombination, and the landscape of mammalian germline variation.

Runaway GC Evolution in Gerbil Genomes

Recombination increases the local GC-content in genomic regions through GC-biased gene conversion (gBGC). The recent discovery of a large genomic region with extreme GC-content in the fat sand rat Psammomys obesus provides a model to study the effects of gBGC on chromosome evolution. Here, we compare the GC-content and GC-to-AT substitution patterns across protein-coding genes of four gerbil species and two murine rodents (mouse and rat). We find that the known high-GC region is present in all the gerbils, and is characterized by high substitution rates for all mutational categories (AT-to-GC, GC-to-AT, and GC-conservative) both at synonymous and nonsynonymous sites. A higher AT-to-GC than GC-to-AT rate is consistent with the high GC-content. Additionally, we find more than 300 genes outside the known region with outlying values of AT-to-GC synonymous substitution rates in gerbils. Of these, over 30% are organized into at least 17 large clusters observable at the megabase-scale. The unusual GC-skewed substitution pattern suggests the evolution of genomic regions with very high recombination rates in the gerbil lineage, which can lead to a runaway increase in GC-content. Our results imply that rapid evolution of GC-content is possible in mammals, with gerbil species providing a powerful model to study the mechanisms of gBGC.

Sycp1 Is Not Required for Subtelomeric DNA Double-Strand Breaks but Is Required for Homologous Alignment in Zebrafish Spermatocytes

In meiotic prophase I, homologous chromosomes are bound together by the synaptonemal complex, in which two axial elements are connected by transverse filaments and central element proteins. In human and zebrafish spermatocytes, homologous recombination and assembly of the synaptonemal complex initiate predominantly near telomeres. In mice, synapsis is not required for meiotic double-strand breaks (DSBs) and homolog alignment but is required for DSB repair; however, the interplay of these meiotic events in the context of peritelomeric bias remains unclear. In this study, we identified a premature stop mutation in the zebrafish gene encoding the transverse filament protein Sycp1. In sycp1 mutant zebrafish spermatocytes, axial elements were formed and paired at chromosome ends between homologs during early to mid-zygonema. However, they did not synapse, and their associations were mostly lost in late zygotene- or pachytene-like stages. In sycp1 mutant spermatocytes, γH2AX signals were observed, and Dmc1/Rad51 and RPA signals appeared predominantly near telomeres, resembling wild-type phenotypes. We observed persistent localization of Hormad1 along the axis in sycp1 mutant spermatocytes, while the majority of Iho1 signals appeared and disappeared with kinetics similar to those in wild-type spermatocytes. Notably, persistent Iho1 foci were observed in spo11 mutant spermatocytes, suggesting that Iho1 dissociation from axes occurs in a DSB-dependent manner. Our results demonstrated that Sycp1 is not required for peritelomeric DSB formation but is necessary for complete pairing of homologs in zebrafish meiosis.

Two telomeric ends of acrocentric chromosome play distinct roles in homologous chromosome synapsis in the fetal mouse oocyte

In mammalian oocytes, proper chromosome segregation at the first meiotic division is dictated by the presence and site of homologous chromosome recombination, which takes place in fetal life. Our current understanding of how homologous chromosomes find each other and initiate synapsis, which is prerequisite for homologous recombination, is limited. It is known that chromosome telomeres are anchored into the nuclear envelope (NE) at the early meiotic prophase I (MPI) and move along NE to facilitate homologous chromosome search and pairing. However, the mouse (Mus musculus) carries all acrocentric chromosomes with one telomeric end close to the centromere (subcentromeric telomere; C-telomere) and the other far away from the centromere (distal telomere; D-telomere), and how C- and D-telomeres participate in chromosome pairing and synapsis during the MPI progression is not well understood. Here, we found in the mouse oocyte that C- and D-telomeres transiently clustered in one area, but D-telomeres soon separated together from C-telomeres and then dispersed to preferentially initiate synapsis, while C-telomeres remained in clusters and synapsed at the last. In the Spo11 null oocyte, which is deficient in SPO11-dependent DSBs formation and homologous synapsis, the pattern of C- and D-telomere clustering and resolution was not affected, but synapsis was more frequently initiated at C-telomeres. These results suggest that SPO11 suppresses the early synapsis between C-telomeres in clusters.

Exo1 recruits Cdc5 polo kinase to MutLγ to ensure efficient meiotic crossover formation

Crossovers generated during the repair of programmed meiotic double-strand breaks must be tightly regulated to promote accurate homolog segregation without deleterious outcomes, such as aneuploidy. The Mlh1-Mlh3 (MutLγ) endonuclease complex is critical for crossover resolution, which involves mechanistically unclear interplay between MutLγ and Exo1 and polo kinase Cdc5. Using budding yeast to gain temporal and genetic traction on crossover regulation, we find that MutLγ constitutively interacts with Exo1. Upon commitment to crossover repair, MutLγ-Exo1 associate with recombination intermediates, followed by direct Cdc5 recruitment that triggers MutLγ crossover activity. We propose that Exo1 serves as a central coordinator in this molecular interplay, providing a defined order of interaction that prevents deleterious, premature activation of crossovers. MutLγ associates at a lower frequency near centromeres, indicating that spatial regulation across chromosomal regions reduces risky crossover events. Our data elucidate the temporal and spatial control surrounding a constitutive, potentially harmful, nuclease. We also reveal a critical, noncatalytic role for Exo1, through noncanonical interaction with polo kinase. These mechanisms regulating meiotic crossovers may be conserved across species.

PRDM9 activity depends on HELLS and promotes local 5-hydroxymethylcytosine enrichment

Meiotic recombination starts with the formation of DNA double-strand breaks (DSBs) at specific genomic locations that correspond to PRDM9-binding sites. The molecular steps occurring from PRDM9 binding to DSB formation are unknown. Using proteomic approaches to find PRDM9 partners, we identified HELLS, a member of the SNF2-like family of chromatin remodelers. Upon functional analyses during mouse male meiosis, we demonstrated that HELLS is required for PRDM9 binding and DSB activity at PRDM9 sites. However, HELLS is not required for DSB activity at PRDM9-independent sites. HELLS is also essential for 5-hydroxymethylcytosine (5hmC) enrichment at PRDM9 sites. Analyses of 5hmC in mice deficient for SPO11, which catalyzes DSB formation, and in PRDM9 methyltransferase deficient mice reveal that 5hmC is triggered at DSB-prone sites upon PRDM9 binding and histone modification, but independent of DSB activity. These findings highlight the complex regulation of the chromatin and epigenetic environments at PRDM9-specified hotspots.

Sexual dimorphism in the meiotic requirement for PRDM9: A mammalian evolutionary safeguard

In many mammals, genomic sites for recombination are determined by the histone methyltransferase PRMD9. Some mouse strains lacking PRDM9 are infertile, but instances of fertility or semifertility in the absence of PRDM9 have been reported in mice, canines, and a human female. Such findings raise the question of how the loss of PRDM9 is circumvented to maintain fertility. We show that genetic background and sex-specific modifiers can obviate the requirement for PRDM9 in mice. Specifically, the meiotic DNA damage checkpoint protein CHK2 acts as a modifier allowing female-specific fertility in the absence of PRDM9. We also report that, in the absence of PRDM9, a PRDM9-independent recombination system is compatible with female meiosis and fertility, suggesting sex-specific regulation of meiotic recombination, a finding with implications for speciation.

Meiotic Double-Strand Break Processing and Crossover Patterning Are Regulated in a Sex-Specific Manner by BRCA1-BARD1 in Caenorhabditis elegans

Meiosis is regulated in a sex-specific manner to produce two distinct gametes, sperm and oocytes, for sexual reproduction. To determine how meiotic recombination is regulated in spermatogenesis, we analyzed the meiotic phenotypes of mutants in the tumor suppressor E3 ubiquitin ligase BRC-1-BRD-1 complex in Caenorhabditis elegans male meiosis. Unlike in mammals, this complex is not required for meiotic sex chromosome inactivation, the process whereby hemizygous sex chromosomes are transcriptionally silenced. Interestingly, brc-1 and brd-1 mutants show meiotic recombination phenotypes that are largely opposing to those previously reported for female meiosis. Fewer meiotic recombination intermediates marked by the recombinase RAD-51 were observed in brc-1 and brd-1 mutants, and the reduction in RAD-51 foci could be suppressed by mutation of nonhomologous-end-joining proteins. Analysis of GFP::RPA-1 revealed fewer foci in the brc-1brd-1 mutant and concentration of BRC-1-BRD-1 to sites of meiotic recombination was dependent on DNA end resection, suggesting that the complex regulates the processing of meiotic double-strand breaks to promote repair by homologous recombination. Further, BRC-1-BRD-1 is important to promote progeny viability when male meiosis is perturbed by mutations that block the pairing and synapsis of different chromosome pairs, although the complex is not required to stabilize the RAD-51 filament as in female meiosis under the same conditions. Analyses of crossover designation and formation revealed that BRC-1-BRD-1 inhibits supernumerary COs when meiosis is perturbed. Together, our findings suggest that BRC-1-BRD-1 regulates different aspects of meiotic recombination in male and female meiosis.

Mechanisms for the Approach/Avoidance Decision Applied to Autism

As a neurodevelopmental disorder with serious lifelong consequences, autism has received considerable attention from neuroscientists and geneticists. We present a hypothesis of mechanisms plausibly affected during brain development in autism, based on neural pathways that are associated with social behavior and connect the prefrontal cortex (PFC) to the basal ganglia (BG). We consider failure of social approach in autism as a special case of imbalance in the fundamental dichotomy between behavioral approach and avoidance. Differential combinations of genes mutated, differences in the timing of their impact during development, and graded degrees of hormonal influences may help explain the heterogeneity in symptomatology in autism and predominance in boys.

The Configuration of RPA, RAD51, and DMC1 Binding in Meiosis Reveals the Nature of Critical Recombination Intermediates

Meiotic recombination proceeds via binding of RPA, RAD51, and DMC1 to single-stranded DNA (ssDNA) substrates created after formation of programmed DNA double-strand breaks. Here we report high-resolution in vivo maps of RPA and RAD51 in meiosis, mapping their binding locations and lifespans to individual homologous chromosomes using a genetically engineered hybrid mouse. Together with high-resolution microscopy and DMC1 binding maps, we show that DMC1 and RAD51 have distinct spatial localization on ssDNA: DMC1 binds near the break site, and RAD51 binds away from it. We characterize inter-homolog recombination intermediates bound by RPA in vivo, with properties expected for the critical displacement loop (D-loop) intermediates. These data support the hypothesis that DMC1, not RAD51, performs strand exchange in mammalian meiosis. RPA-bound D-loops can be resolved as crossovers or non-crossovers, but crossover-destined D-loops may have longer lifespans. D-loops resemble crossover gene conversions in size, but their extent is similar in both repair pathways.

Ensuring meiotic DNA break formation in the mouse pseudoautosomal region

Sex chromosomes in males of most eutherian species share only a diminutive homologous segment, the pseudoautosomal region (PAR), wherein double-strand break (DSB) formation, pairing, and crossing over must occur for correct meiotic segregation^1,2. How cells ensure PAR recombination is unknown. Here we delineate an unexpected dynamic ultrastructure of the PAR and identify controlling cis- and trans-acting factors that make this the hottest area of DSB formation in the male mouse genome. Before break formation, multiple DSB-promoting factors hyper-accumulate in the PAR, its chromosome axes elongate, and the sister chromatids separate. These phenomena are linked to heterochromatic mo-2 minisatellite arrays and require MEI4 and ANKRD31 proteins but not axis components REC8 or HORMAD1. We propose that the repetitive PAR sequence confers unique chromatin and higher order structures crucial for recombination. Chromosome synapsis triggers collapse of the elongated PAR structure and, remarkably, oocytes can be reprogrammed to display spermatocyte-like PAR DSB levels simply by delaying or preventing synapsis. Thus, sexually dimorphic behavior of the PAR rests in part on kinetic differences between the sexes for a race between maturation of PAR structure, DSB formation, and completion of pairing and synapsis. Our findings establish a mechanistic paradigm of sex chromosome recombination during meiosis.

foxl3, a sexual switch in germ cells, initiates two independent molecular pathways for commitment to oogenesis in medaka

Germ cells have the ability to differentiate into eggs and sperm and must determine their sexual fate. In vertebrates, the mechanism of commitment to oogenesis following the sexual fate decision in germ cells remains unknown. Forkhead-box protein L3 (foxl3) is a switch gene involved in the germline sexual fate decision in the teleost fish medaka (Oryzias latipes). Here, we show that foxl3 organizes two independent pathways of oogenesis regulated by REC8 meiotic recombination protein a (rec8a), a cohesin component, and F-box protein (FBP) 47 (fbxo47), a subunit of E3 ubiquitin ligase. In mutants of either gene, germ cells failed to undergo oogenesis but developed normally into sperm in testes. Disruption of rec8a resulted in arrest at a meiotic pachytenelike stage specifically in females, revealing a sexual difference in meiotic progression. Analyses of fbxo47 mutants showed that this gene regulates transcription factors that facilitate folliculogenesis: LIM homeobox 8 (lhx8b), factor in the germline α (figla), and newborn ovary homeobox (nobox). Interestingly, we found that the fbxo47 pathway ensures that germ cells do not deviate from an oogenic pathway until they reach diplotene stage. The mutant phenotypes together with the timing of their expression imply that germline feminization is established during early meiotic prophase I.

Dual histone methyl reader ZCWPW1 facilitates repair of meiotic double strand breaks in male mice

Meiotic crossovers result from homology-directed repair of DNA double-strand breaks (DSBs). Unlike yeast and plants, where DSBs are generated near gene promoters, in many vertebrates DSBs are enriched at hotspots determined by the DNA binding activity of the rapidly evolving zinc finger array of PRDM9 (PR domain zinc finger protein 9). PRDM9 subsequently catalyzes tri-methylation of lysine 4 and lysine 36 of Histone H3 in nearby nucleosomes. Here, we identify the dual histone methylation reader ZCWPW1, which is tightly co-expressed during spermatogenesis with Prdm9, as an essential meiotic recombination factor required for efficient repair of PRDM9-dependent DSBs and for pairing of homologous chromosomes in male mice. In sum, our results indicate that the evolution of a dual histone methylation writer/reader (PRDM9/ZCWPW1) system in vertebrates remodeled genetic recombination hotspot selection from an ancestral static pattern near genes towards a flexible pattern controlled by the rapidly evolving DNA binding activity of PRDM9.

DNA methylation in bull spermatozoa: evolutionary impacts, interindividual variability, and contribution to the embryo

The DNA methylome of spermatozoa results from a unique epigenetic reprogramming crucial for chromatin compaction and the protection of the paternal genetic heritage. Although bull semen is widely used for artificial insemination (AI), little is known about the sperm epigenome in cattle. The purpose of this review is to synthetize recent work on the bull sperm methylome in light of the knowledge accumulated in humans and model species. We will address sperm-specific DNA methylation features and their potential evolutionary impacts, with particular emphasis on hypomethylated regions and repetitive elements. We will review recent examples of interindividual variability and intra-individual plasticity of the bull sperm methylome as related to fertility and age, respectively. Finally, we will address paternal methylome reprogramming after fertilization, as well as the mechanisms potentially involved in epigenetic inheritance, and provide some examples of disturbances that alter the dynamics of reprogramming in cattle. Because the selection of AI bulls is closely based on their genotypes, we will also discuss the complex interplay between sequence polymorphism and DNA methylation, which represents both a difficulty in addressing the role of DNA methylation in shaping phenotypes and an opportunity to better understand genome plasticity.

ATM and PRDM9 regulate SPO11-bound recombination intermediates during meiosis

Meiotic recombination is initiated by SPO11-induced double-strand breaks (DSBs). In most mammals, the methyltransferase PRDM9 guides SPO11 targeting, and the ATM kinase controls meiotic DSB numbers. Following MRE11 nuclease removal of SPO11, the DSB is resected and loaded with DMC1 filaments for homolog invasion. Here, we demonstrate the direct detection of meiotic DSBs and resection using END-seq on mouse spermatocytes with low sample input. We find that DMC1 limits both minimum and maximum resection lengths, whereas 53BP1, BRCA1 and EXO1 play surprisingly minimal roles. Through enzymatic modifications to END-seq, we identify a SPO11-bound meiotic recombination intermediate (SPO11-RI) present at all hotspots. We propose that SPO11-RI forms because chromatin-bound PRDM9 asymmetrically blocks MRE11 from releasing SPO11. In Atm-/- spermatocytes, trapped SPO11 cleavage complexes accumulate due to defective MRE11 initiation of resection. Thus, in addition to governing SPO11 breakage, ATM and PRDM9 are critical local regulators of mammalian SPO11 processing.

Refined spatial temporal epigenomic profiling reveals intrinsic connection between PRDM9-mediated H3K4me3 and the fate of double-stranded breaks

Meiotic recombination is initiated by the formation of double-strand breaks (DSBs), which are repaired as either crossovers (COs) or noncrossovers (NCOs). In most mammals, PRDM9-mediated H3K4me3 controls the nonrandom distribution of DSBs; however, both the timing and mechanism of DSB fate control remain largely undetermined. Here, we generated comprehensive epigenomic profiles of synchronized mouse spermatogenic cells during meiotic prophase I, revealing spatiotemporal and functional relationships between epigenetic factors and meiotic recombination. We find that PRDM9-mediated H3K4me3 at DSB hotspots, coinciding with H3K27ac and H3K36me3, is intimately connected with the fate of the DSB. Our data suggest that the fate decision is likely made at the time of DSB formation: earlier formed DSBs occupy more open chromatins and are much more competent to proceed to a CO fate. Our work highlights an intrinsic connection between PRDM9-mediated H3K4me3 and the fate decision of DSBs, and provides new insight into the control of CO homeostasis.

Sex Differences in the Recombination Landscape

Sex differences in overall recombination rates are well known, but little theoretical or empirical attention has been given to how and why sexes differ in their recombination landscapes: the patterns of recombination along chromosomes. In the first scientific review of this phenomenon, we find that recombination is biased toward telomeres in males and more uniformly distributed in females in most vertebrates and many other eukaryotes. Notable exceptions to this pattern exist, however. Fine-scale recombination patterns also frequently differ between males and females. The molecular mechanisms responsible for sex differences remain unclear, but chromatin landscapes play a role. Why these sex differences evolve also is unclear. Hypotheses suggest that they may result from sexually antagonistic selection acting on coding genes and their regulatory elements, meiotic drive in females, selection during the haploid phase of the life cycle, selection against aneuploidy, or mechanistic constraints. No single hypothesis, however, can adequately explain the evolution of sex differences in all cases. Sex-specific recombination landscapes have important consequences for population differentiation and sex chromosome evolution.

Genomic Structure of Hstx2 Modifier of Prdm9-Dependent Hybrid Male Sterility in Mice

F1 hybrids between mouse inbred strains PWD and C57BL/6 represent the most thoroughly genetically defined model of hybrid sterility in vertebrates. Hybrid male sterility can be fully reconstituted from three components of this model, the Prdm9 gene, intersubspecific homeology of Mus musculus musculus and Mus musculus domesticus autosomes, and the X-linked Hstx2 locus. Hstx2 modulates the extent of Prdm9-dependent meiotic arrest and harbors two additional factors responsible for intersubspecific introgression-induced oligospermia (Hstx1) and meiotic recombination rate (Meir1). To facilitate positional cloning and to overcome the recombination suppression within the 4.3 Mb encompassing the Hstx2 locus, we designed Hstx2-CRISPR and SPO11/Cas9 transgenes aimed to induce DNA double-strand breaks specifically within the Hstx2 locus. The resulting recombinant reduced the Hstx2 locus to 2.70 Mb (chromosome X: 66.51-69.21 Mb). The newly defined Hstx2 locus still operates as the major X-linked factor of the F1 hybrid sterility, and controls meiotic chromosome synapsis and meiotic recombination rate. Despite extensive further crosses, the 2.70 Mb Hstx2 interval behaved as a recombination cold spot with reduced PRDM9-mediated H3K4me3 hotspots and absence of DMC1-defined DNA double-strand-break hotspots. To search for structural anomalies as a possible cause of recombination suppression, we used optical mapping and observed high incidence of subspecies-specific structural variants along the X chromosome, with a striking copy number polymorphism of the microRNA Mir465 cluster. This observation together with the absence of a strong sterility phenotype in Fmr1 neighbor (Fmr1nb) null mutants support the role of microRNA as a likely candidate for Hstx2.

A high-resolution map of non-crossover events reveals impacts of genetic diversity on mammalian meiotic recombination

During meiotic recombination, homologue-templated repair of programmed DNA double-strand breaks (DSBs) produces relatively few crossovers and many difficult-to-detect non-crossovers. By intercrossing two diverged mouse subspecies over five generations and deep-sequencing 119 offspring, we detect thousands of crossover and non-crossover events genome-wide with unprecedented power and spatial resolution. We find that both crossovers and non-crossovers are strongly depleted at DSB hotspots where the DSB-positioning protein PRDM9 fails to bind to the unbroken homologous chromosome, revealing that PRDM9 also functions to promote homologue-templated repair. Our results show that complex non-crossovers are much rarer in mice than humans, consistent with complex events arising from accumulated non-programmed DNA damage. Unexpectedly, we also find that GC-biased gene conversion is restricted to non-crossover tracts containing only one mismatch. These results demonstrate that local genetic diversity profoundly alters meiotic repair pathway decisions via at least two distinct mechanisms, impacting genome evolution and Prdm9-related hybrid infertility.

During meiotic recombination, genetic information is transferred or exchanged between parental chromosome copies. Using a large hybrid mouse pedigree, the authors generated high-resolution maps of these transfer/exchange events and discovered new properties governing their processing and resolution.

The histone modification reader ZCWPW1 is required for meiosis prophase I in male but not in female mice

Meiosis is a specialized type of cell division that creates haploid germ cells and ensures their genetic diversity through homologous recombination. We show that the H3K4me3 reader ZCWPW1 is specifically required for meiosis prophase I progression in male but not in female germ cells in mice. Loss of Zcwpw1 in male mice caused a complete failure of synapsis, resulting in meiotic arrest at the zygotene to pachytene stage, accompanied by incomplete DNA double-strand break repair and lack of crossover formation, leading to male infertility. In oocytes, deletion of Zcwpw1 only somewhat slowed down meiosis prophase I progression; Zcwpw1^-/- oocytes were able to complete meiosis, and Zcwpw1^-/- female mice had normal fertility until mid-adulthood. We conclude that the H3K4me3 reader ZCWPW1 is indispensable for meiosis synapsis in males but is dispensable for females. Our results suggest that ZCWPW1 may represent a previously unknown, sex-dependent epigenetic regulator of germ cell meiosis in mammals.

Heterogeneous transposable elements as silencers, enhancers and targets of meiotic recombination

During meiosis, DNA double-strand breaks are initiated by the topoisomerase-like enzyme SPO11 and are repaired by inter-sister chromatid and inter-homologue DNA repair pathways. Genome-wide maps of initiating DNA double-strand breaks and inter-homologue repair events are now available for a number of mammalian, fungal and plant species. In mammals, PRDM9 specifies the location of meiotic recombination initiation via recognition of specific DNA sequence motifs by its C2H2 zinc finger array. In fungi and plants, meiotic recombination appears to be initiated less discriminately in accessible chromatin, including at gene promoters. Generally, meiotic crossover is suppressed in highly repetitive genomic regions that are made up of transposable elements (TEs), to prevent deleterious non-allelic homologous recombination events. However, recent and older studies have revealed intriguing relationships between meiotic recombination initiation and repair, and transposable elements. For instance, gene conversion events have been detected in maize centromeric retroelements, mouse MULE-MuDR DNA transposons undergo substantial meiotic recombination initiation, Arabidopsis Helitron TEs are among the hottest of recombination initiation hotspots, and human TE sequences can modify the crossover rate at adjacent PRDM9 motifs in cis. Here, we summarize the relationship between meiotic recombination and TEs, discuss recent insights from highly divergent eukaryotes and highlight outstanding questions in the field.