Evolution of recombination in eutherian mammals: insights into mechanisms that affect recombination rates and crossover interference

Toward telomere-to-telomere cat genomes for precision medicine and conservation biology

Genomic data from species of the cat family Felidae promise to stimulate veterinary and human medical advances, and clarify the coherence of genome organization. We describe how interspecies hybrids have been instrumental in the genetic analysis of cats, from the first genetic maps to propelling cat genomes toward the T2T standard set by the human genome project. Genotype-to-phenotype mapping in cat models has revealed dozens of health-related genetic variants, the molecular basis for mammalian pigmentation and patterning, and species-specific adaptations. Improved genomic surveillance of natural and captive populations across the cat family tree will increase our understanding of the genetic architecture of traits, population dynamics, and guide a future of genome-enabled biodiversity conservation.

Modeling the evolution of recombination plasticity: A prospective review

Meiotic recombination is one of the main sources of genetic variation, a fundamental factor in the evolutionary adaptation of sexual eukaryotes. Yet, the role of variation in recombination rate and other recombination features remains underexplored. In this review, we focus on the sensitivity of recombination rates to different extrinsic and intrinsic factors. We briefly present the empirical evidence for recombination plasticity in response to environmental perturbations and/or poor genetic background and discuss theoretical models developed to explain how such plasticity could have evolved and how it can affect important population characteristics. We highlight a gap between the evidence, which comes mostly from experiments with diploids, and theory, which typically assumes haploid selection. Finally, we formulate open questions whose solving would help to outline conditions favoring recombination plasticity. This will contribute to answering the long-standing question of why sexual recombination exists despite its costs, since plastic recombination may be evolutionary advantageous even in selection regimes rejecting any non-zero constant recombination.

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.

Synaptonemal & CO analyzer: A tool for synaptonemal complex and crossover analysis in immunofluorescence images

During the formation of ova and sperm, homologous chromosomes get physically attached through the synaptonemal complex and exchange DNA at crossover sites by a process known as meiotic recombination. Chromosomes that do not recombine or have anomalous crossover distributions often separate poorly during the subsequent cell division and end up in abnormal numbers in ova or sperm, which can lead to miscarriage or developmental defects. Crossover numbers and distribution along the synaptonemal complex can be visualized by immunofluorescent microscopy. However, manual analysis of large numbers of cells is very time-consuming and a major bottleneck for recombination studies. Some image analysis tools have been created to overcome this situation, but they are not readily available, do not provide synaptonemal complex data, or do not tackle common experimental difficulties, such as overlapping chromosomes. To overcome these limitations, we have created and validated an open-source ImageJ macro routine that facilitates and speeds up the crossover and synaptonemal complex analyses in mouse chromosome spreads, as well as in other vertebrate species. It is free, easy to use and fulfills the recommendations for enhancing rigor and reproducibility in biomedical studies.

Garcias M, Kratochvíl L, Valenzuela N, Georges A, Waters PD, Ruiz-Herrera A. Meiotic chromosome dynamics and double strand break formation in reptiles

During meiotic prophase I, tightly regulated processes take place, from pairing and synapsis of homologous chromosomes to recombination, which are essential for the generation of genetically variable haploid gametes. These processes have canonical meiotic features conserved across different phylogenetic groups. However, the dynamics of meiotic prophase I in non-mammalian vertebrates are poorly known. Here, we compare four species from Sauropsida to understand the regulation of meiotic prophase I in reptiles: the Australian central bearded dragon (Pogona vitticeps), two geckos (Paroedura picta and Coleonyx variegatus) and the painted turtle (Chrysemys picta). We first performed a histological characterization of the spermatogenesis process in both the bearded dragon and the painted turtle. We then analyzed prophase I dynamics, including chromosome pairing, synapsis and the formation of double strand breaks (DSBs). We show that meiosis progression is highly conserved in reptiles with telomeres clustering forming the bouquet, which we propose promotes homologous pairing and synapsis, along with facilitating the early pairing of micro-chromosomes during prophase I (i.e., early zygotene). Moreover, we detected low levels of meiotic DSB formation in all taxa. Our results provide new insights into reptile meiosis.

Background selection under evolving recombination rates

Background selection (BGS), the effect that purifying selection exerts on sites linked to deleterious alleles, is expected to be ubiquitous across eukaryotic genomes. The effects of BGS reflect the interplay of the rates and fitness effects of deleterious mutations with recombination. A fundamental assumption of BGS models is that recombination rates are invariant over time. However, in some lineages, recombination rates evolve rapidly, violating this central assumption. Here, we investigate how recombination rate evolution affects genetic variation under BGS. We show that recombination rate evolution modifies the effects of BGS in a manner similar to a localized change in the effective population size, potentially leading to underestimation or overestimation of the genome-wide effects of selection. Furthermore, we find evidence that recombination rate evolution in the ancestors of modern house mice may have impacted inferences of the genome-wide effects of selection in that species.

Fragile, unfaithful and persistent Ys-on how meiosis can shape sex chromosome evolution

Sex-linked inheritance is a stark exception to Mendel's Laws of Heredity. Here we discuss how the evolution of heteromorphic sex chromosomes (mainly the Y) has been shaped by the intricacies of the meiotic programme. We propose that persistence of Y chromosomes in distantly related mammalian phylogroups can be explained in the context of pseudoautosomal region (PAR) size, meiotic pairing strategies, and the presence of Y-borne executioner genes that regulate meiotic sex chromosome inactivation. We hypothesise that variation in PAR size can be an important driver for the evolution of recombination frequencies genome wide, imposing constraints on Y fate. If small PAR size compromises XY segregation during male meiosis, the stress of producing aneuploid gametes could drive function away from the Y (i.e., a fragile Y). The Y chromosome can avoid fragility either by acquiring an achiasmatic meiotic XY pairing strategy to reduce aneuploid gamete production, or gain meiotic executioner protection (a persistent Y). Persistent Ys will then be under strong pressure to maintain high recombination rates in the PAR (and subsequently genome wide), as improper segregation has fatal consequences for germ cells. In the event that executioner protection is lost, the Y chromosome can be maintained in the population by either PAR rejuvenation (extension by addition of autosome material) or gaining achiasmatic meiotic pairing, the alternative is Y loss. Under this dynamic cyclic evolutionary scenario, understanding the meiotic programme in vertebrate and invertebrate species will be crucial to further understand the plasticity of the rise and fall of heteromorphic sex chromosomes.

Strategies for meiotic sex chromosome dynamics and telomeric elongation in Marsupials

During meiotic prophase I, homologous chromosomes pair, synapse and recombine in a tightly regulated process that ensures the generation of genetically variable haploid gametes. Although the mechanisms underlying meiotic cell division have been well studied in model species, our understanding of the dynamics of meiotic prophase I in non-traditional model mammals remains in its infancy. Here, we reveal key meiotic features in previously uncharacterised marsupial species (the tammar wallaby and the fat-tailed dunnart), plus the fat-tailed mouse opossum, with a focus on sex chromosome pairing strategies, recombination and meiotic telomere homeostasis. We uncovered differences between phylogroups with important functional and evolutionary implications. First, sex chromosomes, which lack a pseudo-autosomal region in marsupials, had species specific pairing and silencing strategies, with implications for sex chromosome evolution. Second, we detected two waves of γH2AX accumulation during prophase I. The first wave was accompanied by low γH2AX levels on autosomes, which correlated with the low recombination rates that distinguish marsupials from eutherian mammals. In the second wave, γH2AX was restricted to sex chromosomes in all three species, which correlated with transcription from the X in tammar wallaby. This suggests non-canonical functions of γH2AX on meiotic sex chromosomes. Finally, we uncover evidence for telomere elongation in primary spermatocytes of the fat-tailed dunnart, a unique strategy within mammals. Our results provide new insights into meiotic progression and telomere homeostasis in marsupials, highlighting the importance of capturing the diversity of meiotic strategies within mammals.

The generation of haploid gametes is a hallmark of sexual reproduction. And this is accomplished by a complex, albeit tightly regulated, reductional cell division called meiosis. Although meiosis has been extensively studied in eutherian mammal model species, our understanding of the mechanisms regulating chromosome synapsis, recombination and segregation during meiosis progression is still incomplete especially in non-eutherian mammals. To fill this gap and capture the diversity of meiotic strategies among mammals, we study previously uncharacterised representative marsupial species, an evolutionary assemblage that last shared a common ancestry more than 80 million years ago. We uncover novel, hence non-canonical, strategies for sex chromosome pairing, DNA repair, recombination and transcription. Most importantly, we reveal the uniqueness of marsupial meiosis, which includes the unprecedented detection of alternative mechanism (ALT) for the paternal control of telomere length during prophase I. Our findings suggest that ALT (previously only associated to cancer cells) could play a role in telomere homeostasis in mammalian germ cells.

Negative heterosis for meiotic recombination rate in spermatocytes of the domestic chicken Gallus gallus

Benef its and costs of meiotic recombination are a matter of discussion. Because recombination breaks allele combinations already tested by natural selection and generates new ones of unpredictable f itness, a high recombination rate is generally benef icial for the populations living in a f luctuating or a rapidly changing environment and costly in a stable environment. Besides genetic benef its and costs, there are cytological effects of recombination, both positive and negative. Recombination is necessary for chromosome synapsis and segregation. However, it involves a massive generation of double-strand DNA breaks, erroneous repair of which may lead to germ cell death or various mutations and chromosome rearrangements. Thus, the benef its of recombination (generation of new allele combinations) would prevail over its costs (occurrence of deleterious mutations) as long as the population remains suff iciently heterogeneous. Using immunolocalization of MLH1, a mismatch repair protein, at the synaptonemal complexes, we examined the number and distribution of recombination nodules in spermatocytes of two chicken breeds with high (Pervomai) and low (Russian Crested) recombination rates and their F1 hybrids and backcrosses. We detected negative heterosis for recombination rate in the F1 hybrids. Backcrosses to the Pervomai breed were rather homogenous and showed an intermediate recombination rate. The differences in overall recombination rate between the breeds, hybrids and backcrosses were mainly determined by the differences in the crossing over number in the seven largest macrochromosomes. The decrease in recombination rate in F1 is probably determined by diff iculties in homology matching between the DNA sequences of genetically divergent breeds. The suppression of recombination in the hybrids may impede gene f low between parapatric populations and therefore accelerate their genetic divergence.

Sex-specific variation in the genome-wide recombination rate

In most species that reproduce sexually, successful gametogenesis requires recombination during meiosis. The number and placement of crossovers (COs) vary among individuals, with females and males often presenting the most striking contrasts. Despite the recognition that the sexes recombine at different rates (heterochiasmy), existing data fail to answer the question of whether patterns of genetic variation in recombination rate are similar in the two sexes. To fill this gap, we measured the genome-wide recombination rate in both sexes from a panel of wild-derived inbred strains from multiple subspecies of house mice (Mus musculus) and from a few additional species of Mus. To directly compare recombination rates in females and males from the same genetic backgrounds, we applied established methods based on immunolocalization of recombination proteins to inbred strains. Our results reveal discordant patterns of genetic variation in the two sexes. Whereas male genome-wide recombination rates vary substantially among strains, female recombination rates measured in the same strains are more static. The direction of heterochiasmy varies within two subspecies, Mus musculus molossinus and Mus musculus musculus. The direction of sex differences in the length of the synaptonemal complex and CO positions is consistent across strains and does not track sex differences in genome-wide recombination rate. In males, contrasts between strains with high recombination rate and strains with low recombination rate suggest more recombination is associated with stronger CO interference and more double-strand breaks. The sex-specific patterns of genetic variation we report underscore the importance of incorporating sex differences into recombination research.

Genetic variation in recombination rate in the pig

Background

Meiotic recombination results in the exchange of genetic material between homologous chromosomes. Recombination rate varies between different parts of the genome, between individuals, and is influenced by genetics. In this paper, we assessed the genetic variation in recombination rate along the genome and between individuals in the pig using multilocus iterative peeling on 150,000 individuals across nine genotyped pedigrees. We used these data to estimate the heritability of recombination and perform a genome-wide association study of recombination in the pig.

Results

Our results confirmed known features of the recombination landscape of the pig genome, including differences in genetic length of chromosomes and marked sex differences. The recombination landscape was repeatable between lines, but at the same time, there were differences in average autosome-wide recombination rate between lines. The heritability of autosome-wide recombination rate was low but not zero (on average 0.07 for females and 0.05 for males). We found six genomic regions that are associated with recombination rate, among which five harbour known candidate genes involved in recombination: RNF212, SHOC1, SYCP2, MSH4 and HFM1.

Conclusions

Our results on the variation in recombination rate in the pig genome agree with those reported for other vertebrates, with a low but nonzero heritability, and the identification of a major quantitative trait locus for recombination rate that is homologous to that detected in several other species. This work also highlights the utility of using large-scale livestock data to understand biological processes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00643-0.

The impact of chromosomal fusions on 3D genome folding and recombination in the germ line

The spatial folding of chromosomes inside the nucleus has regulatory effects on gene expression, yet the impact of genome reshuffling on this organization remains unclear. Here, we take advantage of chromosome conformation capture in combination with single-nucleotide polymorphism (SNP) genotyping and analysis of crossover events to study how the higher-order chromatin organization and recombination landscapes are affected by chromosomal fusions in the mammalian germ line. We demonstrate that chromosomal fusions alter the nuclear architecture during meiosis, including an increased rate of heterologous interactions in primary spermatocytes, and alterations in both chromosome synapsis and axis length. These disturbances in topology were associated with changes in genomic landscapes of recombination, resulting in detectable genomic footprints. Overall, we show that chromosomal fusions impact the dynamic genome topology of germ cells in two ways: (i) altering chromosomal nuclear occupancy and synapsis, and (ii) reshaping landscapes of recombination.

Ultracontinuous Single Haplotype Genome Assemblies for the Domestic Cat (Felis catus) and Asian Leopard Cat (Prionailurus bengalensis)

In addition to including one of the most popular companion animals, species from the cat family Felidae serve as a powerful system for genetic analysis of inherited and infectious disease, as well as for the study of phenotypic evolution and speciation. Previous diploid-based genome assemblies for the domestic cat have served as the primary reference for genomic studies within the cat family. However, these versions suffered from poor resolution of complex and highly repetitive regions, with substantial amounts of unplaced sequence that is polymorphic or copy number variable. We sequenced the genome of a female F1 Bengal hybrid cat, the offspring of a domestic cat (Felis catus) x Asian leopard cat (Prionailurus bengalensis) cross, with PacBio long sequence reads and used Illumina sequence reads from the parents to phase >99.9% of the reads into the 2 species' haplotypes. De novo assembly of the phased reads produced highly continuous haploid genome assemblies for the domestic cat and Asian leopard cat, with contig N50 statistics exceeding 83 Mb for both genomes. Whole-genome alignments reveal the Felis and Prionailurus genomes are colinear, and the cytogenetic differences between the homologous F1 and E4 chromosomes represent a case of centromere repositioning in the absence of a chromosomal inversion. Both assemblies offer significant improvements over the previous domestic cat reference genome, with a 100% increase in contiguity and the capture of the vast majority of chromosome arms in 1 or 2 large contigs. We further demonstrated that comparably accurate F1 haplotype phasing can be achieved with members of the same species when one or both parents of the trio are not available. These novel genome resources will empower studies of feline precision medicine, adaptation, and speciation.

Meiotic Nuclear Architecture in Distinct Mole Vole Hybrids with Robertsonian Translocations: Chromosome Chains, Stretched Centromeres, and Distorted Recombination

Genome functioning in hybrids faces inconsistency. This mismatch is manifested clearly in meiosis during chromosome synapsis and recombination. Species with chromosomal variability can be a model for exploring genomic battles with high visibility due to the use of advanced immunocytochemical methods. We studied synaptonemal complexes (SC) and prophase I processes in 44-chromosome intraspecific (Ellobius tancrei × E. tancrei) and interspecific (Ellobius talpinus × E. tancrei) hybrid mole voles heterozygous for 10 Robertsonian translocations. The same pachytene failures were found for both types of hybrids. In the intraspecific hybrid, the chains were visible in the pachytene stage, then 10 closed SC trivalents formed in the late pachytene and diplotene stage. In the interspecific hybrid, as a rule, SC trivalents composed the SC chains and rarely could form closed configurations. Metacentrics involved with SC trivalents had stretched centromeres in interspecific hybrids. Linkage between neighboring SC trivalents was maintained by stretched centromeric regions of acrocentrics. This centromeric plasticity in structure and dynamics of SC trivalents was found for the first time. We assume that stretched centromeres were a marker of altered nuclear architecture in heterozygotes due to differences in the ancestral chromosomal territories of the parental species. Restructuring of the intranuclear organization and meiotic disturbances can contribute to the sterility of interspecific hybrids, and lead to the reproductive isolation of studied species.

MLH1 focus mapping in the guinea fowl (Numida meleagris) give insights into the crossover landscapes in birds

Crossover rates and localization are not homogeneous throughout the genomes. Along the chromosomes of almost all species, domains with high crossover rates alternate with domains where crossover rates are significantly lower than the genome-wide average. The distribution of crossovers along chromosomes constitutes the recombination landscape of a given species and can be analyzed at broadscale using immunostaining of the MLH1 protein, a component of mature recombination nodules found on synaptonemal complexes during pachytene. We scored the MLH1 foci in oocytes of the chicken and the guinea fowl and compared their frequencies in the largest bivalents. The average autosomal number of foci is 62 in the chicken and 44 in the guinea fowl. The lower number in the guinea fowl responds to the occurrence of fewer crossovers in the six largest bivalents, where most MLH1 foci occur within one-fifth of the chromosome length with high polarization towards opposite ends. The skewed distribution of foci in the guinea fowl contrast with the more uniform distribution of numerous foci in the chicken, especially in the four largest bivalents. The crossover distribution observed in the guinea fowl is unusual among Galloanserae and also differs from other, more distantly related birds. We discussed the current evidence showing that the shift towards crossover localization, as observed in the guinea fowl, was not a unique event but also occurred at different moments of bird evolution. A comparative analysis of genome-wide average recombination rates in birds shows variations within narrower limits compared to mammals and the absence of a phylogenetic trend.

Crossover Interference: Shedding Light on the Evolution of Recombination

Through recombination, genes are freed to evolve more independently of one another, unleashing genetic variance hidden in the linkage disequilibrium that accumulates through selection combined with drift. Yet crossover numbers are evolutionarily constrained, with at least one and not many more than one crossover per bivalent in most taxa. Crossover interference, whereby a crossover reduces the probability of a neighboring crossover, contributes to this homogeneity. The mechanisms by which interference is achieved and crossovers are regulated are a major current subject of inquiry, facilitated by novel methods to visualize crossovers and to pinpoint recombination events. Here, we review patterns of crossover interference and the models built to describe this process. We then discuss the selective forces that have likely shaped interference and the regulation of crossover numbers.

Molecular evolution of the meiotic recombination pathway in mammals

Meiotic recombination shapes evolution and helps to ensure proper chromosome segregation in most species that reproduce sexually. Recombination itself evolves, with species showing considerable divergence in the rate of crossing-over. However, the genetic basis of this divergence is poorly understood. Recombination events are produced via a complicated, but increasingly well-described, cellular pathway. We apply a phylogenetic comparative approach to a carefully selected panel of genes involved in the processes leading to crossovers-spanning double-strand break formation, strand invasion, the crossover/non-crossover decision, and resolution-to reconstruct the evolution of the recombination pathway in eutherian mammals and identify components of the pathway likely to contribute to divergence between species. Eleven recombination genes, predominantly involved in the stabilization of homologous pairing and the crossover/non-crossover decision, show evidence of rapid evolution and positive selection across mammals. We highlight TEX11 and associated genes involved in the synaptonemal complex and the early stages of the crossover/non-crossover decision as candidates for the evolution of recombination rate. Evolutionary comparisons to MLH1 count, a surrogate for the number of crossovers, reveal a positive correlation between genome-wide recombination rate and the rate of evolution at TEX11 across the mammalian phylogeny. Our results illustrate the power of viewing the evolution of recombination from a pathway perspective.

CENP-A binding domains and recombination patterns in horse spermatocytes

Centromeres exert an inhibitory effect on meiotic recombination, but the possible contribution of satellite DNA to this "centromere effect" is under debate. In the horse, satellite DNA is present at all centromeres with the exception of the one from chromosome 11. This organization of centromeres allowed us to investigate the role of satellite DNA on recombination suppression in horse spermatocytes at the stage of pachytene. To this aim we analysed the distribution of the MLH1 protein, marker of recombination foci, relative to CENP-A, marker of centromeric function. We demonstrated that the satellite-less centromere of chromosome 11 causes crossover suppression, similarly to satellite-based centromeres. These results suggest that the centromere effect does not depend on satellite DNA. During this analysis, we observed a peculiar phenomenon: while, as expected, the centromere of the majority of meiotic bivalent chromosomes was labelled with a single immunofluorescence centromeric signal, double-spotted or extended signals were also detected. Their number varied from 0 to 7 in different cells. This observation can be explained by positional variation of the centromeric domain on the two homologs and/or misalignment of pericentromeric satellite DNA arrays during homolog pairing confirming the great plasticity of equine centromeres.

Conservation of the genome-wide recombination rate in white-footed mice

Despite being linked to the fundamental processes of chromosome segregation and offspring diversification, meiotic recombination rates vary within and between species. Recent years have seen progress in quantifying recombination rate evolution across multiple temporal and genomic scales. Nevertheless, the level of variation in recombination rate within wild populations-a key determinant of evolution in this trait-remains poorly documented on the genomic scale. To address this notable gap, we used immunofluorescent cytology to quantify genome-wide recombination rates in males from a wild population of the white-footed mouse, Peromyscus leucopus. For comparison, we measured recombination rates in a second population of male P. leucopus raised in the laboratory and in male deer mice from the subspecies Peromyscus maniculatus bairdii. Although we found differences between individuals in the genome-wide recombination rate, levels of variation were low-within populations, between populations, and between species. Quantification of synaptonemal complex length and crossover positions along chromosome 1 using a novel automated approach also revealed conservation in broad-scale crossover patterning, including strong crossover interference. We propose stabilizing selection targeting recombination or correlated processes as the explanation for these patterns.

Interbreed variation in meiotic recombination rate and distribution in the domestic chicken Gallus gallus

The efficiency of natural and artificial selection is critically dependent on the recombination rate. However, interbreed and individual variation in recombination rate in poultry remains unknown. Conventional methods of analysis of recombination such as genetic linkage analysis, sperm genotyping and chiasma count at lampbrush chromosomes are expensive and time-consuming. In this study, we analyzed the number and distribution of recombination nodules in spermatocytes of the roosters of six chicken breeds using immunolocalization of key proteins involved in chromosome pairing and recombination. We revealed significant effects of breed (R 2 = 0.17 ; p < 0.001 ) and individual (R 2 = 0.28 ; p < 0.001 ) on variation in the number of recombination nodules. Both interbreed and individual variations in recombination rate were almost entirely determined by variation in recombination density on macrochromosomes, because almost all microchromosomes in each breed had one recombination nodule. Despite interbreed differences in the density of recombination nodules, the patterns of their distribution along homologous chromosomes were similar. The breeds examined in this study showed a correspondence between the age of the breed and its recombination rate. Those with high recombination rates (Pervomai, Russian White and Brahma) are relatively young breeds created by crossing several local breeds. The breeds displaying low recombination rate are ancient local breeds: Cochin (Indo-China), Brown Leghorn (Tuscany, Italy) and Russian Crested (the European part of Russia).

Desiccation-induced changes in recombination rate and crossover interference in Drosophila melanogaster: evidence for fitness-dependent plasticity

Meiotic recombination is evolutionarily ambiguous, as being associated with both benefits and costs to its bearers, with the resultant dependent on a variety of conditions. While existing theoretical models explain the emergence and maintenance of recombination, some of its essential features remain underexplored. Here we focus on one such feature, recombination plasticity, and test whether recombination response to stress is fitness-dependent. We compare desiccation stress effects on recombination rate and crossover interference in chromosome 3 between desiccation-sensitive and desiccation-tolerant Drosophila lines. We show that relative to desiccation-tolerant genotypes, desiccation-sensitive genotypes exhibit a significant segment-specific increase in single- and double-crossover frequencies across the pericentromeric region of chromosome 3. Significant changes (relaxation) in crossover interference were found for the interval pairs flanking the centromere and extending to the left arm of the chromosome. These results indicate that desiccation is a recombinogenic factor and that desiccation-induced changes in both recombination rate and crossover interference are fitness-dependent, with a tendency of less fitted individuals to produce more variable progeny. Such dependence may play an important role in the regulation of genetic variation in populations experiencing environmental challenges.

Per-Nucleus Crossover Covariation and Implications for Evolution

Crossing over is a nearly universal feature of sexual reproduction. Here, analysis of crossover numbers on a per-chromosome and per-nucleus basis reveals a fundamental, evolutionarily conserved feature of meiosis: within individual nuclei, crossover frequencies covary across different chromosomes. This effect results from per-nucleus covariation of chromosome axis lengths. Crossovers can promote evolutionary adaptation. However, the benefit of creating favorable new allelic combinations must outweigh the cost of disrupting existing favorable combinations. Covariation concomitantly increases the frequencies of gametes with especially high, or especially low, numbers of crossovers, and thus might concomitantly enhance the benefits of crossing over while reducing its costs. A four-locus population genetic model suggests that such an effect can pertain in situations where the environment fluctuates: hyper-crossover gametes are advantageous when the environment changes while hypo-crossover gametes are advantageous in periods of environmental stasis. These findings reveal a new feature of the basic meiotic program and suggest a possible adaptive advantage.

A rigorous measure of genome-wide genetic shuffling that takes into account crossover positions and Mendel's second law

Comparative studies in evolutionary genetics rely critically on evaluation of the total amount of genetic shuffling that occurs during gamete production. Such studies have been hampered by the absence of a direct measure of this quantity. Existing measures consider crossing-over by simply counting the average number of crossovers per meiosis. This is qualitatively inadequate, because the positions of crossovers along a chromosome are also critical: a crossover toward the middle of a chromosome causes more shuffling than a crossover toward the tip. Moreover, traditional measures fail to consider shuffling from independent assortment of homologous chromosomes (Mendel's second law). Here, we present a rigorous measure of genome-wide shuffling that does not suffer from these limitations. We define the parameter [Formula: see text] as the probability that the alleles at two randomly chosen loci are shuffled during gamete production. This measure can be decomposed into separate contributions from crossover number and position and from independent assortment. Intrinsic implications of this metric include the fact that [Formula: see text] is larger when crossovers are more evenly spaced, which suggests a selective advantage of crossover interference. Utilization of [Formula: see text] is enabled by powerful emergent methods for determining crossover positions either cytologically or by DNA sequencing. Application of our analysis to such data from human male and female reveals that (i) [Formula: see text] in humans is close to its maximum possible value of 1/2 and that (ii) this high level of shuffling is due almost entirely to independent assortment, the contribution of which is ∼30 times greater than that of crossovers.

Meiotic Recombination: Genetics' Good Old Scalpel

In the era of genome engineering, a new study returns to classical genetics to decipher genotype-phenotype relationships in unprecedented throughput and with unprecedented accuracy. Capitalizing on natural variation in yeast strains and frequent meiotic recombination, She and Jarosz (2018) dissect and map to nucleotide resolution, simple and complex determinants of diverse phenotypic traits.

Connecting theory and data to understand recombination rate evolution

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

What drives the evolution of condition-dependent recombination in diploids? Some insights from simulation modelling

While the evolutionary advantages of non-zero recombination rates have prompted diverse theoretical explanations, the evolution of essential recombination features remains underexplored. We focused on one such feature, the condition dependence of recombination, viewed as the variation in within-generation sensitivity of recombination to external (environment) and/or internal (genotype) conditions. Limited empirical evidence for its existence comes mainly from diploids, whereas theoretical models show that it only easily evolves in haploids. The evolution of condition-dependent recombination can be explained by its advantage for the selected system (indirect effect), or by benefits to modifier alleles, ensuring this strategy regardless of effects on the selected system (direct effect). We considered infinite panmictic populations of diploids exposed to a cyclical two-state environment. Each organism had three selected loci. Examining allele dynamics at a fourth, selectively neutral recombination modifier locus, we frequently observed that a modifier allele conferring condition-dependent recombination between the selected loci displaced the allele conferring the optimal constant recombination rate. Our simulations also confirm the results of theoretical studies showing that condition-dependent recombination cannot evolve in diploids on the basis of direct fitness-dependent effects alone. Therefore, the evolution of condition-dependent recombination in diploids can be driven by indirect effects alone, i.e. by modifier effects on the selected system.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Reorganization of the Y Chromosomes Enhances Divergence in Israeli Mole Rats Nannospalax ehrenbergi (Spalacidae, Rodentia): Comparative Analysis of Meiotic and Mitotic Chromosomes

The Y chromosome in mammals is variable, even in closely related species. Middle East blind mole rats Nannospalax ehrenbergi demonstrate autosomal variability, which probably leads to speciation. Here, we compare the mitotic and meiotic chromosomes of mole rats. For the first time, we studied the behavior of their sex chromosomes in the meiotic prophase I using electron microscopy and immunocytochemical analysis. Unexpectedly, the sex chromosomes of the 52- and 60-chromosome forms of mole rats showed different synaptic and recombination patterns due to distinct locations of the centromeres on the Y chromosomes. The absence of recombination in the 60-chromosome form, the asymmetric synapsis, and the short-term disturbance in the synaptic co-orientation of the telomeric regions of the X and Y chromosomes were revealed as specific features of mole rat sex bivalents. We suggest several scenarios of Y chromosome alteration in connection with species differentiation in mole rats.

Chromosome Synapsis and Recombination in Male-Sterile and Female-Fertile Interspecies Hybrids of the Dwarf Hamsters (Phodopus, Cricetidae)

Hybrid sterility is an important step in the speciation process. Hybrids between dwarf hamsters Phodopus sungorus and P.campbelli provide a good model for studies in cytological and genetic mechanisms of hybrid sterility. Previous studies in hybrids detected multiple abnormalities of spermatogenesis and a high frequency of dissociation between the X and Y chromosomes at the meiotic prophase. In this study, we found that the autosomes of the hybrid males and females underwent paring and recombination as normally as their parental forms did. The male hybrids showed a significantly higher frequency of asynapsis and recombination failure between the heterochromatic arms of the X and Y chromosomes than the males of the parental species. Female hybrids as well as the females of the parental species demonstrated a high incidence of centromere misalignment at the XX bivalent and partial asynapsis of the ends of its heterochromatic arms. In all three karyotypes, recombination was completely suppressed in the heterochromatic arm of the X chromosome, where the pseudoautosomal region is located. We propose that this recombination pattern speeds up divergence of the X- and Y-linked pseudoautosomal regions between the parental species and results in their incompatibility in the male hybrids.

Mammalian Comparative Genomics Reveals Genetic and Epigenetic Features Associated with Genome Reshuffling in Rodentia

Understanding how mammalian genomes have been reshuffled through structural changes is fundamental to the dynamics of its composition, evolutionary relationships between species and, in the long run, speciation. In this work, we reveal the evolutionary genomic landscape in Rodentia, the most diverse and speciose mammalian order, by whole-genome comparisons of six rodent species and six representative outgroup mammalian species. The reconstruction of the evolutionary breakpoint regions across rodent phylogeny shows an increased rate of genome reshuffling that is approximately two orders of magnitude greater than in other mammalian species here considered. We identified novel lineage and clade-specific breakpoint regions within Rodentia and analyzed their gene content, recombination rates and their relationship with constitutive lamina genomic associated domains, DNase I hypersensitivity sites and chromatin modifications. We detected an accumulation of protein-coding genes in evolutionary breakpoint regions, especially genes implicated in reproduction and pheromone detection and mating. Moreover, we found an association of the evolutionary breakpoint regions with active chromatin state landscapes, most probably related to gene enrichment. Our results have two important implications for understanding the mechanisms that govern and constrain mammalian genome evolution. The first is that the presence of genes related to species-specific phenotypes in evolutionary breakpoint regions reinforces the adaptive value of genome reshuffling. Second, that chromatin conformation, an aspect that has been often overlooked in comparative genomic studies, might play a role in modeling the genomic distribution of evolutionary breakpoints.

Mapping Causal Variants with Single-Nucleotide Resolution Reveals Biochemical Drivers of Phenotypic Change

Understanding the sequence determinants that give rise to diversity among individuals and species is the central challenge of genetics. However, despite ever greater numbers of sequenced genomes, most genome-wide association studies cannot distinguish causal variants from linked passenger mutations spanning many genes. We report that this inherent challenge can be overcome in model organisms. By pushing the advantages of inbred crossing to its practical limit in Saccharomyces cerevisiae, we improved the statistical resolution of linkage analysis to single nucleotides. This "super-resolution" approach allowed us to map 370 causal variants across 26 quantitative traits. Missense, synonymous, and cis-regulatory mutations collectively gave rise to phenotypic diversity, providing mechanistic insight into the basis of evolutionary divergence. Our data also systematically unmasked complex genetic architectures, revealing that multiple closely linked driver mutations frequently act on the same quantitative trait. Single-nucleotide mapping thus complements traditional deletion and overexpression screening paradigms and opens new frontiers in quantitative genetics.

Broad-scale recombination pattern in the primitive bird Rhea americana (Ratites, Palaeognathae)

Birds have genomic and chromosomal features that make them an attractive group to analyze the evolution of recombination rate and the distribution of crossing over. Yet, analyses are biased towards certain species, especially domestic poultry and passerines. Here we analyze for the first time the recombination rate and crossover distribution in the primitive ratite bird, Rhea americana (Rheiformes, Palaeognathae). Using a cytogenetic approach for in situ mapping of crossovers we found that the total genetic map is 3050 cM with a global recombination rate of 2.1 cM/Mb for female rheas. In the five largest macrobivalents there were 3 or more crossovers in most bivalents. Recombination rates for macrobivalents ranges between 1.8-2.1 cM/Mb and the physical length of their synaptonemal complexes is highly predictive of their genetic lengths. The crossover rate at the pseudoautosomal region is 2.1 cM/Mb, similar to those of autosomal pairs 5 and 6 and only slightly higher compared to other macroautosomes. It is suggested that the presence of multiple crossovers on the largest macrobivalents is a feature common to many avian groups, irrespective of their position throughout phylogeny. These data provide new insights to analyze the heterogeneous recombination landscape of birds.

Chromosome Synapsis and Recombination in Male Hybrids between Two Chromosome Races of the Common Shrew (Sorex araneus L., Soricidae, Eulipotyphla)

Hybrid zones between chromosome races of the common shrew (Sorex araneus) provide exceptional models to study the potential role of chromosome rearrangements in the initial steps of speciation. The Novosibirsk and Tomsk races differ by a series of Robertsonian fusions with monobrachial homology. They form a narrow hybrid zone and generate hybrids with both simple (chain of three chromosomes) and complex (chain of eight or nine) synaptic configurations. Using immunolocalisation of the meiotic proteins, we examined chromosome pairing and recombination in males from the hybrid zone. Homozygotes and simple heterozygotes for Robertsonian fusions showed a low frequency of synaptic aberrations (<10%). The carriers of complex synaptic configurations showed multiple pairing abnormalities, which might lead to reduced fertility. The recombination frequency in the proximal regions of most chromosomes of all karyotypes was much lower than in the other regions. The strong suppression of recombination in the pericentromeric regions and co-segregation of race specific chromosomes involved in the long chains would be expected to lead to linkage disequilibrium between genes located there. Genic differentiation, together with the high frequency of pairing aberrations in male carriers of the long chains, might contribute to maintenance of the narrow hybrid zone.

Landscaping Crossover Interference Across a Genome

The evolutionary success of eukaryotic organisms crucially depends on the capacity to produce genetic diversity through reciprocal exchanges of each chromosome pair, or crossovers (COs), during meiosis. It has been recognized that COs arise more evenly across a given chromosome than at random. This phenomenon, termed CO interference, occurs pervasively in eukaryotes and may confer a selective advantage. We describe here a multipoint linkage analysis procedure for segregating families to quantify the strength of CO interference over the genome, and extend this procedure to illustrate the landscape of CO interference in natural populations. We further discuss the crucial role of CO interference in amplifying and maintaining genetic diversity through sex-, stress-, and age-induced differentiation.

Recombination correlates with synaptonemal complex length and chromatin loop size in bovids-insights into mammalian meiotic chromosomal organization

Homologous chromosomes exchange genetic information through recombination during meiosis, a process that increases genetic diversity, and is fundamental to sexual reproduction. In an attempt to shed light on the dynamics of mammalian recombination and its implications for genome organization, we have studied the recombination characteristics of 112 individuals belonging to 28 different species in the family Bovidae. In particular, we analyzed the distribution of RAD51 and MLH1 foci during the meiotic prophase I that serve, respectively, as proxies for double-strand breaks (DSBs) which form in early stages of meiosis and for crossovers. In addition, synaptonemal complex length and meiotic DNA loop size were estimated to explore how genome organization determines DSBs and crossover patterns. We show that although the number of meiotic DSBs per cell and recombination rates observed vary between individuals of the same species, these are correlated with diploid number as well as with synaptonemal complex and DNA loop sizes. Our results illustrate that genome packaging, DSB frequencies, and crossover rates tend to be correlated, while meiotic chromosomal axis length and DNA loop size are inversely correlated in mammals. Moreover, axis length, DSB frequency, and crossover frequencies all covary, suggesting that these correlations are established in the early stages of meiosis.

Recombination rate variation in mice from an isolated island

Recombination rate is a heritable trait that varies among individuals. Despite the major impact of recombination rate on patterns of genetic diversity and the efficacy of selection, natural variation in this phenotype remains poorly characterized. We present a comparison of genetic maps, sampling 1212 meioses, from a unique population of wild house mice (Mus musculus domesticus) that recently colonized remote Gough Island. Crosses to a mainland reference strain (WSB/EiJ) reveal pervasive variation in recombination rate among Gough Island mice, including subchromosomal intervals spanning up to 28% of the genome. In spite of this high level of polymorphism, the genomewide recombination rate does not significantly vary. In general, we find that recombination rate varies more when measured in smaller genomic intervals. Using the current standard genetic map of the laboratory mouse to polarize intervals with divergent recombination rates, we infer that the majority of evolutionary change occurred in one of the two tested lines of Gough Island mice. Our results confirm that natural populations harbour a high level of recombination rate polymorphism and highlight the disparities in recombination rate evolution across genomic scales.

Variation and Evolution of the Meiotic Requirement for Crossing Over in Mammals

The segregation of homologous chromosomes at the first meiotic division is dependent on the presence of at least one well-positioned crossover per chromosome. In some mammalian species, however, the genomic distribution of crossovers is consistent with a more stringent baseline requirement of one crossover per chromosome arm. Given that the meiotic requirement for crossing over defines the minimum frequency of recombination necessary for the production of viable gametes, determining the chromosomal scale of this constraint is essential for defining crossover profiles predisposed to aneuploidy and understanding the parameters that shape patterns of recombination rate evolution across species. Here, I use cytogenetic methods for in situ imaging of crossovers in karyotypically diverse house mice (Mus musculus domesticus) and voles (genus Microtus) to test how chromosome number and configuration constrain the distribution of crossovers in a genome. I show that the global distribution of crossovers in house mice is thresholded by a minimum of one crossover per chromosome arm, whereas the crossover landscape in voles is defined by a more relaxed requirement of one crossover per chromosome. I extend these findings in an evolutionary metaanalysis of published recombination and karyotype data for 112 mammalian species and demonstrate that the physical scale of the genomic crossover distribution has undergone multiple independent shifts from one crossover per chromosome arm to one per chromosome during mammalian evolution. Together, these results indicate that the chromosomal scale constraint on crossover rates is itself a trait that evolves among species, a finding that casts light on an important source of crossover rate variation in mammals.

Genome analysis and avirulence gene cloning using a high-density RADseq linkage map of the flax rust fungus, Melampsora lini

BACKGROUND

Rust fungi are an important group of plant pathogens that cause devastating losses in agricultural, silvicultural and natural ecosystems. Plants can be protected from rust disease by resistance genes encoding receptors that trigger a highly effective defence response upon recognition of specific pathogen avirulence proteins. Identifying avirulence genes is crucial for understanding how virulence evolves in the field.

RESULTS

To facilitate avirulence gene cloning in the flax rust fungus, Melampsora lini, we constructed a high-density genetic linkage map using single nucleotide polymorphisms detected in restriction site-associated DNA sequencing (RADseq) data. The map comprises 13,412 RADseq markers in 27 linkage groups that together span 5860 cM and contain 2756 recombination bins. The marker sequences were used to anchor 68.9 % of the M. lini genome assembly onto the genetic map. The map and anchored assembly were then used to: 1) show that M. lini has a high overall meiotic recombination rate, but recombination distribution is uneven and large coldspots exist; 2) show that substantial genome rearrangements have occurred in spontaneous loss-of-avirulence mutants; and 3) identify the AvrL2 and AvrM14 avirulence genes by map-based cloning. AvrM14 is a dual-specificity avirulence gene that encodes a predicted nudix hydrolase. AvrL2 is located in the region of the M. lini genome with the lowest recombination rate and encodes a small, highly-charged proline-rich protein.

CONCLUSIONS

The M. lini high-density linkage map has greatly advanced our understanding of virulence mechanisms in this pathogen by providing novel insights into genome variability and enabling identification of two new avirulence genes.

Unique sex chromosome systems in Ellobius: How do male XX chromosomes recombine and undergo pachytene chromatin inactivation?

Most mammalian species have heteromorphic sex chromosomes in males, except for a few enigmatic groups such as the mole voles Ellobius, which do not have the Y chromosome and Sry gene. The Ellobius (XX ♀♂) system of sex chromosomes has no analogues among other animals. The structure and meiotic behaviour of the two X chromosomes were investigated for males of the sibling species Ellobius talpinus and Ellobius tancrei. Their sex chromosomes, despite their identical G-structure, demonstrate short synaptic fragments and crossover-associated MLH1 foci in both telomeric regions only. The chromatin undergoes modifications in the meiotic sex chromosomes. SUMO-1 marks a small nucleolus-like body of the meiotic XX. ATR and ubiH2A are localized in the asynaptic area and the histone γH2AFX covers the entire XX bivalent. The distribution of some markers of chromatin inactivation differentiates sex chromosomes of mole voles from those of other mammals. Sex chromosomes of both studied species have identical recombination and meiotic inactivation patterns. In Ellobius, similar chromosome morphology masks the functional heteromorphism of the male sex chromosomes, which can be seen at meiosis.

A High-Resolution SNP Array-Based Linkage Map Anchors a New Domestic Cat Draft Genome Assembly and Provides Detailed Patterns of Recombination

High-resolution genetic and physical maps are invaluable tools for building accurate genome assemblies, and interpreting results of genome-wide association studies (GWAS). Previous genetic and physical maps anchored good quality draft assemblies of the domestic cat genome, enabling the discovery of numerous genes underlying hereditary disease and phenotypes of interest to the biomedical science and breeding communities. However, these maps lacked sufficient marker density to order thousands of shorter scaffolds in earlier assemblies, which instead relied heavily on comparative mapping with related species. A high-resolution map would aid in validating and ordering chromosome scaffolds from existing and new genome assemblies. Here, we describe a high-resolution genetic linkage map of the domestic cat genome based on genotyping 453 domestic cats from several multi-generational pedigrees on the Illumina 63K SNP array. The final maps include 58,055 SNP markers placed relative to 6637 markers with unique positions, distributed across all autosomes and the X chromosome. Our final sex-averaged maps span a total autosomal length of 4464 cM, the longest described linkage map for any mammal, confirming length estimates from a previous microsatellite-based map. The linkage map was used to order and orient the scaffolds from a substantially more contiguous domestic cat genome assembly (Felis catus v8.0), which incorporated ∼20 × coverage of Illumina fragment reads. The new genome assembly shows substantial improvements in contiguity, with a nearly fourfold increase in N50 scaffold size to 18 Mb. We use this map to report probable structural errors in previous maps and assemblies, and to describe features of the recombination landscape, including a massive (∼50 Mb) recombination desert (of virtually zero recombination) on the X chromosome that parallels a similar desert on the porcine X chromosome in both size and physical location.

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

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

Meiotic behaviour of evolutionary sex-autosome translocations in Bovidae

The recurrent occurrence of sex-autosome translocations during mammalian evolution suggests common mechanisms enabling a precise control of meiotic synapsis, recombination and inactivation of sex chromosomes. We used immunofluorescence and FISH to study the meiotic behaviour of sex chromosomes in six species of Bovidae with evolutionary sex-autosome translocations (Tragelaphus strepsiceros, Taurotragus oryx, Tragelaphus imberbis, Tragelaphus spekii, Gazella leptoceros and Nanger dama ruficollis). The autosomal regions of fused sex chromosomes showed normal synapsis with their homologous counterparts. Synapsis in the pseudoautosomal region (PAR) leads to the formation of characteristic bivalent (in T. imberbis and T. spekii with X;BTA13/Y;BTA13), trivalent (in T. strepsiceros and T. oryx with X/Y;BTA13 and G. leptoceros with X;BTA5/Y) and quadrivalent (in N. dama ruficollis with X;BTA5/Y;BTA16) structures at pachynema. However, when compared with other mammals, the number of pachynema lacking MLH1 foci in the PAR was relatively high, especially in T. imberbis and T. spekii, species with both sex chromosomes involved in sex autosome translocations. Meiotic transcriptional inactivation of the sex-autosome translocations assessed by γH2AX staining was restricted to their gonosomal regions. Despite intraspecies differences, the evolutionary fixation of sex-autosome translocations among bovids appears to involve general mechanisms ensuring sex chromosome pairing, synapsis, recombination and inactivation.

Recombination patterns reveal information about centromere location on linkage maps

Linkage mapping is often used to identify genes associated with phenotypic traits and for aiding genome assemblies. Still, many emerging maps do not locate centromeres - an essential component of the genomic landscape. Here, we demonstrate that for genomes with strong chiasma interference, approximate centromere placement is possible by phasing the same data used to generate linkage maps. Assuming one obligate crossover per chromosome arm, information about centromere location can be revealed by tracking the accumulated recombination frequency along linkage groups, similar to half-tetrad analyses. We validate the method on a linkage map for sockeye salmon (Oncorhynchus nerka) with known centromeric regions. Further tests suggest that the method will work well in other salmonids and other eukaryotes. However, the method performed weakly when applied to a male linkage map (rainbow trout; O. mykiss) characterized by low and unevenly distributed recombination - a general feature of male meiosis in many species. Further, a high frequency of double crossovers along chromosome arms in barley reduced resolution for locating centromeric regions on most linkage groups. Despite these limitations, our method should work well for high-density maps in species with strong recombination interference and will enrich many existing and future mapping resources.

Variation of Meiotic Recombination Rates and MLH1 Foci Distribution in Spermatocytes of Cattle, Sheep and Goats

Despite similar genome sizes, a great variability in recombination rates is observed in mammals. We used antibodies against SYCP3, MLH1 and centromeres to compare crossover frequency, position along chromosome arms and the effect of crossover interference in spermatocytes of 4 species from the family Bovidae (Bos taurus, 2n = 60, tribe Bovini; Ovis aries, 2n = 54, Capra hircus, 2n = 60 and Ammotragus lervia, 2n = 58, tribe Caprini). Despite significant individual variability, our results also show significant differences in both recombination rates and the total length of autosomal synaptonemal complexes (SC) between cattle (47.53 MLH1 foci/cell, 244.59 µm) and members of the tribe Caprini (61.83 MLH1 foci, 296.19 µm) which can be explained by the length of time that has passed since their evolutionary divergence. Sheep displayed the highest number of MLH1 foci per cell and recombination density, although they have a lower diploid chromosome number caused by centric fusions corresponding to cattle chromosomes 1;3, 2;8 and 5;11. However, the proportion of MLH1 foci observed on the fused chromosomes in sheep (26.14%) was significantly lower than on the orthologous acrocentrics in cattle (27.6%) and goats (28.2%), and their distribution along the SC arms differed significantly. The reduced recombination rate in metacentrics is probably caused by interference acting across the centromere.

Linkage Mapping Reveals Strong Chiasma Interference in Sockeye Salmon: Implications for Interpreting Genomic Data

Meiotic recombination is fundamental for generating new genetic variation and for securing proper disjunction. Further, recombination plays an essential role during the rediploidization process of polyploid-origin genomes because crossovers between pairs of homeologous chromosomes retain duplicated regions. A better understanding of how recombination affects genome evolution is crucial for interpreting genomic data; unfortunately, current knowledge mainly originates from a few model species. Salmonid fishes provide a valuable system for studying the effects of recombination in nonmodel species. Salmonid females generally produce thousands of embryos, providing large families for conducting inheritance studies. Further, salmonid genomes are currently rediploidizing after a whole genome duplication and can serve as models for studying the role of homeologous crossovers on genome evolution. Here, we present a detailed interrogation of recombination patterns in sockeye salmon (Oncorhynchus nerka). First, we use RAD sequencing of haploid and diploid gynogenetic families to construct a dense linkage map that includes paralogous loci and location of centromeres. We find a nonrandom distribution of paralogs that mainly cluster in extended regions distally located on 11 different chromosomes, consistent with ongoing homeologous recombination in these regions. We also estimate the strength of interference across each chromosome; results reveal strong interference and crossovers are mostly limited to one per arm. Interference was further shown to continue across centromeres, but metacentric chromosomes generally had at least one crossover on each arm. We discuss the relevance of these findings for both mapping and population genomic studies.

Effect of species-specific differences in chromosome morphology on chromatin compaction and the frequency and distribution of RAD51 and MLH1 foci in two bovid species: cattle (Bos taurus) and the common eland (Taurotragus oryx)

Meiotic recombination between homologous chromosomes is crucial for their correct segregation into gametes and for generating diversity. We compared the frequency and distribution of MLH1 foci and RAD51 foci, synaptonemal complex (SC) length and DNA loop size in two related Bovidae species that share chromosome arm homology but show an extreme difference in their diploid chromosome number: cattle (Bos taurus, 2n = 60) and the common eland (Taurotragus oryx, 2nmale = 31). Compared to cattle, significantly fewer MLH1 foci per cell were observed in the common eland, which can be attributed to the lower number of initial double-strand breaks (DSBs) detected as RAD51 foci in leptonema. Despite the significantly shorter total autosomal SC length and longer DNA loop size of the common eland bi-armed chromosomes compared to those of bovine acrocentrics, the overall crossover density in the common eland was still lower than in cattle, probably due to the reduction in the number of MLH1 foci in the proximal regions of the bi-armed chromosomes. The formation of centric fusions during karyotype evolution of the common eland accompanied by meiotic chromatin compaction has greater implications in the reduction in the number of DSBs in leptonema than in the decrease of MLH1 foci number in pachynema.

Genomic correlates of recombination rate and its variability across eight recombination maps in the western honey bee (Apis mellifera L.)

Background

Meiotic recombination has traditionally been explained based on the structural requirement to stabilize homologous chromosome pairs to ensure their proper meiotic segregation. Competing hypotheses seek to explain the emerging findings of significant heterogeneity in recombination rates within and between genomes, but intraspecific comparisons of genome-wide recombination patterns are rare. The honey bee (Apis mellifera) exhibits the highest rate of genomic recombination among multicellular animals with about five cross-over events per chromatid.

Results

Here, we present a comparative analysis of recombination rates across eight genetic linkage maps of the honey bee genome to investigate which genomic sequence features are correlated with recombination rate and with its variation across the eight data sets, ranging in average marker spacing ranging from 1 Mbp to 120 kbp. Overall, we found that GC content explained best the variation in local recombination rate along chromosomes at the analyzed 100 kbp scale. In contrast, variation among the different maps was correlated to the abundance of microsatellites and several specific tri- and tetra-nucleotides.

Conclusions

The combined evidence from eight medium-scale recombination maps of the honey bee genome suggests that recombination rate variation in this highly recombining genome might be due to the DNA configuration instead of distinct sequence motifs. However, more fine-scale analyses are needed. The empirical basis of eight differing genetic maps allowed for robust conclusions about the correlates of the local recombination rates and enabled the study of the relation between DNA features and variability in local recombination rates, which is particularly relevant in the honey bee genome with its exceptionally high recombination rate.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1281-2) contains supplementary material, which is available to authorized users.

Sex-specific recombination maps for individual macrochromosomes in the Japanese quail (Coturnix japonica)

Meiotic recombination in the Japanese quail was directly studied by immunolocalization of mutL homolog 1 (MLH1), a mismatch repair protein of mature recombination nodules. In total, 15,862 crossovers were scored along the autosomal synaptonemal complexes in 308 meiotic nuclei from males and females. Crossover frequencies calculated from MLH1 foci show wide similitude between males and females with slightly higher number of foci in females. From this analysis, we predict that the sex-averaged map length of the Japanese quail is 2580 cM, with a genome-wide recombination rate of 1.9 cM/Mb. MLH1 focus mapping along the six largest bivalents showed few intersex differences in the distribution of crossovers along with variant patterns in metacentric and acrocentric macrobivalents. These results provide valuable information to complement linkage map analysis in the species while providing insight into our understanding of the mechanisms of crossover distribution along chromosome arms.

Genetic recombination variation in wild Robertsonian mice: on the role of chromosomal fusions and Prdm9 allelic background

Despite the existence of formal models to explain how chromosomal rearrangements can be fixed in a population in the presence of gene flow, few empirical data are available regarding the mechanisms by which genome shuffling contributes to speciation, especially in mammals. In order to shed light on this intriguing evolutionary process, here we present a detailed empirical study that shows how Robertsonian (Rb) fusions alter the chromosomal distribution of recombination events during the formation of the germline in a Rb system of the western house mouse (Mus musculus domesticus). Our results indicate that both the total number of meiotic crossovers and the chromosomal distribution of recombination events are reduced in mice with Rb fusions and that this can be related to alterations in epigenetic signatures for heterochromatinization. Furthermore, we detected novel house mouse Prdm9 allelic variants in the Rb system. Remarkably, mean recombination rates were positively correlated with a decrease in the number of ZnF domains in the Prdm9 gene. The suggestion that recombination can be modulated by both chromosomal reorganizations and genetic determinants that control the formation of double-stranded breaks during meiosis opens new avenues for understanding the role of recombination in chromosomal speciation.

Rapid and inexpensive whole-genome genotyping-by-sequencing for crossover localization and fine-scale genetic mapping

The reshuffling of existing genetic variation during meiosis is important both during evolution and in breeding. The reassortment of genetic variants relies on the formation of crossovers (COs) between homologous chromosomes. The pattern of genome-wide CO distributions can be rapidly and precisely established by the short-read sequencing of individuals from F₂ populations, which in turn are useful for quantitative trait locus (QTL) mapping. Although sequencing costs have decreased precipitously in recent years, the costs of library preparation for hundreds of individuals have remained high. To enable rapid and inexpensive CO detection and QTL mapping using low-coverage whole-genome sequencing of large mapping populations, we have developed a new method for library preparation along with Trained Individual GenomE Reconstruction, a probabilistic method for genotype and CO predictions for recombinant individuals. In an example case with hundreds of F₂ individuals from two Arabidopsis thaliana accessions, we resolved most CO breakpoints to within 2 kb and reduced a major flowering time QTL to a 9-kb interval. In addition, an extended region of unusually low recombination revealed a 1.8-Mb inversion polymorphism on the long arm of chromosome 4. We observed no significant differences in the frequency and distribution of COs between F₂ individuals with and without a functional copy of the DNA helicase gene RECQ4A. In summary, we present a new, cost-efficient method for large-scale, high-precision genotyping-by-sequencing.