Cattle Sex-Specific Recombination and Genetic Control from a Large Pedigree Analysis

Meiotic recombination is an essential biological process that generates genetic diversity and ensures proper segregation of chromosomes during meiosis. From a large USDA dairy cattle pedigree with over half a million genotyped animals, we extracted 186,927 three-generation families, identified over 8.5 million maternal and paternal recombination events, and constructed sex-specific recombination maps for 59,309 autosomal SNPs. The recombination map spans for 25.5 Morgans in males and 23.2 Morgans in females, for a total studied region of 2,516 Mb (986 kb/cM in males and 1,085 kb/cM in females). The male map is 10% longer than the female map and the sex difference is most pronounced in the subtelomeric regions. We identified 1,792 male and 1,885 female putative recombination hotspots, with 720 hotspots shared between sexes. These hotspots encompass 3% of the genome but account for 25% of the genome-wide recombination events in both sexes. During the past forty years, males showed a decreasing trend in recombination rate that coincided with the artificial selection for milk production. Sex-specific GWAS analyses identified PRDM9 and CPLX1 to have significant effects on genome-wide recombination rate in both sexes. Two novel loci, NEK9 and REC114, were associated with recombination rate in both sexes, whereas three loci, MSH4, SMC3 and CEP55, affected recombination rate in females only. Among the multiple PRDM9 paralogues on the bovine genome, our GWAS of recombination hotspot usage together with linkage analysis identified the PRDM9 paralogue on chromosome 1 to be associated in the U.S. Holstein data. Given the largest sample size ever reported for such studies, our results reveal new insights into the understanding of cattle and mammalian recombination.     

A pronounced sex difference when two linked loci of the Japanese brown frogRana japonica are recombined

To investigate the sex differences of crossing-over inRana japonica, the recombination frequencies between two linked loci,LDH-B andMPI, were examined in offspring obtained from 29 matings involving either females or males heterozygous at both loci. No recombinants were found among 303 offspring bred from 6 heterozygous males, whereas 160 were parental and 138 were recombinants among 298 offspring bred from 9 heterozygous females. The recombination frequency was 46.3%. These results show that the sex difference when two loci were recombined was marked and significant inRana japonica.     

Genetic linkage mapping of allozyme loci in even- and odd-year pink salmon (Oncorhynchus gorbuscha)

We constructed genetic linkage maps of allozyme loci in even- and odd-year pink salmon (Oncorhynchus gorbuscha), using the total of 320 families (each female was crossed with two different males, and 80 females and 160 males were used for each of even year and odd year). The maps include eight linkage groups involving 22 loci. We observed substantial variation in recombination frequencies among different families within broodline and between sexes within broodlines. In the linkage analysis between sAAT-3* and sMDH-B1,2*, two even-year families and one odd-year family exhibited evidence of association, but two even-year and one odd-year families did not. Recombination rate tends to be reduced in males in pink salmon. The ratio of recombination rate (female/male), which ranged from 1.7 to infinity, averaged 2.8 in the even-year crosses and 3.2 in the odd-year crosses. The linkage groups (LG) I and II involving sAAT and mAH loci, which probably duplicated in the recent tetraploidization event, and the orders of loci in the LGs I (sAAT-3* --> mAH-4*) and II (mAH-3* --> sAAT-4*) were reversed, suggesting the possible paracentric inversion during salmonid evolution after the duplication.     

Variation in recombination rate across the genome: evidence and implications

Recent data from humans and other species provide convincing evidence of variation in recombination rate in different genomic regions. Comparison of physical and genetic maps reveals variation on a scale of megabases, with substantial differences between sexes. Recombination is often suppressed near centromeres and elevated near telomeres, but neither of these observations is true for all chromosomes. In humans, patterns of linkage disequilibrium and experimental measures of recombination from sperm-typing reveal dramatic hotspots of recombination on a scale of kilobases. Genome-wide variation in the amount of crossing-over may be due to variation in the density of hotspots, the intensity of hotspots, or both. Theoretical models of selection and linkage predict that genetic variation will be reduced in regions of low recombination, and this prediction is supported by data from several species. Heterogeneity in rates of crossing-over provides both an opportunity and a challenge for identifying disease genes: as associations occur in blocks, genomic regions containing disease loci may be identified with relatively few markers, yet identifying the causal mutations is unlikely to be achieved through associations alone.     

A coalescent model of recombination hotspots

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

Causes and consequences of crossing-over evidenced via a high-resolution recombinational landscape of the honey bee

BackgroundSocial hymenoptera, the honey bee (Apis mellifera) in particular, have ultra-high crossover rates and a large degree of intra-genomic variation in crossover rates. Aligned with haploid genomics of males, this makes them a potential model for examining the causes and consequences of crossing over. To address why social insects have such high crossing-over rates and the consequences of this, we constructed a high-resolution recombination atlas by sequencing 55 individuals from three colonies with an average marker density of 314 bp/marker.ResultsWe find crossing over to be especially high in proximity to genes upregulated in worker brains, but see no evidence for a coupling with immune-related functioning. We detect only a low rate of non-crossover gene conversion, contrary to current evidence. This is in striking contrast to the ultrahigh crossing-over rate, almost double that previously estimated from lower resolution data. We robustly recover the predicted intragenomic correlations between crossing over and both population level diversity and GC content, which could be best explained as indirect and direct consequences of crossing over, respectively.ConclusionsOur data are consistent with the view that diversification of worker behavior, but not immune function, is a driver of the high crossing-over rate in bees. While we see both high diversity and high GC content associated with high crossing-over rates, our estimate of the low non-crossover rate demonstrates that high non-crossover rates are not a necessary consequence of high recombination rates.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0566-0) contains supplementary material, which is available to authorized users.     

Prospects for association mapping in classical inbred mouse strains

The collection of classical inbred mouse strains displays heritable variation in a large number of complex traits. Many generations of historical recombination have contributed to the panel of classical strain genomes, raising the possibility that quantitative trait loci could be located with high resolution by correlating strain genotypes and phenotypes. Although this association mapping framework has been successful in several empirical applications, its expected performance remains unclear. We used computer simulations based on a publicly available, dense single-nucleotide polymorphism (SNP) map to measure the power and false-positive rate of association mapping on a genomic scale across 30 commonly used classical inbred strains. Expected power is (i) often low for phenotypic effect sizes that are realistic for complex traits, (ii) highly variable across the genome, and (iii) correlated with linkage disequilibrium, aspects of the allele frequency distribution, and haplotype characteristics, as predicted by theory. Simulations also demonstrate clear potential for spurious associations to be generated by unequal relatedness among the strains. These findings suggest that association mapping in the classical strains is best applied in combination with other procedures, such as QTL mapping.     

Lack of correlation between crossing-over and chromosome distance between inversions in Drosophila ananassae

Two linked inversions, AL and ZE, located in the opposite limbs of the second chromosome of Drosophila ananassae are separated from each other by nearly 32% of the total length of the second chromosome. Crossing-over between these inversions when heterozygous was studied in females and males by the salivary-gland smear technique using karyotypically homozygous stocks. The results of recombination experiments show that there is a strong suppression of recombination between inversions when heterozygous, in spite of a large euchromatic distance available for crossing-over between them. Thus there is no correlation between chromosome distance and crossing-over between heterozygous inversions in the second chromosome of D. ananassae when studied cytologically.     

Genetic Variants in REC8, RNF212, and PRDM9 Influence Male Recombination in Cattle

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

Whole‐genome resequencing of extreme phenotypes in collared flycatchers highlights the difficulty of detecting quantitative trait loci in natural populations

Dissecting the genetic basis of phenotypic variation in natural populations is a long‐standing goal in evolutionary biology. One open question is whether quantitative traits are determined only by large numbers of genes with small effects, or whether variation also exists in large‐effect loci. We conducted genomewide association analyses of forehead patch size (a sexually selected trait) on 81 whole‐genome‐resequenced male collared flycatchers with extreme phenotypes, and on 415 males sampled independent of patch size and genotyped with a 50K SNP chip. No SNPs were genomewide statistically significantly associated with patch size. Simulation‐based power analyses suggest that the power to detect large‐effect loci responsible for 10% of phenotypic variance was <0.5 in the genome resequencing analysis, and <0.1 in the SNP chip analysis. Reducing the recombination by two‐thirds relative to collared flycatchers modestly increased power. Tripling sample size increased power to >0.8 for resequencing of extreme phenotypes (N = 243), but power remained <0.2 for the 50K SNP chip analysis (N = 1245). At least 1 million SNPs were necessary to achieve power >0.8 when analysing 415 randomly sampled phenotypes. However, power of the 50K SNP chip to detect large‐effect loci was nearly 0.8 in simulations with a small effective population size of 1500. These results suggest that reliably detecting large‐effect trait loci in large natural populations will often require thousands of individuals and near complete sampling of the genome. Encouragingly, far fewer individuals and loci will often be sufficient to reliably detect large‐effect loci in small populations with widespread strong linkage disequilibrium.     

A primary genetic linkage map for human chromosome 12

A primary genetic map for human chromosome 12 has been constructed from data on 23 restriction fragment length polymorphic systems collected in 38 normal families with large sibships. Linkage analysis of the genotypic data has ordered 16 loci into a continuous genetic map of 111 cM in males and 258 cM in females. Although most of the genetic map reflects a higher rate of recombination in females relative to males, significantly more frequent recombination was observed in males than in females in intervals between loci on the distal portion of the short arm of the chromosome. The mapping data shown here will serve as a first step toward a high-resolution genetic map for human chromosome 12.     

Spontaneous recombination in males ofDrosophila bipectinata

Spontaneous recombination in males ofDrosophila bipectinata was tested in five wild type laboratory stocks of different geographic origins by using sepia eye and black body colour double recessive mutant stock. The results indicate thatDrosophila bipectinata exhibits spontaneous male recombination. Further, recombination occurs at low rate and there is interstrain variation with respect to the rate of male crossing-over. This is the first report of spontaneous recombination in males ofDrosophila bipectinata.     

A site specific increase in recombination in Drosophila ananassae

The e65 pi; bri ru stock of Drosophila ananassae produced an extremely high rate of recombination in males when made heterozygous with any one of the wild type stocks. We analyzed and characterized the genetic factors which caused this phenomenon. We show that the second chromosome of the e65 pi; bri ru stock carries an enhancer of male recombination. The enhancer, En(2)-ep, is located between Om(2C) and Arc. The enhancement of meiotic recombination both in males and females was also observed at the specific region between Om(2C) and Arc on 2L. The magnitude of increased recombination was 30-40 fold in males and 13-30 fold in females. The relation between the hotspot of recombination in both sexes and the enhancer of male recombination is discussed.     

Evolution of the genomic recombination rate in murid rodents

Although very closely related species can differ in their fine-scale patterns of recombination hotspots, variation in the average genomic recombination rate among recently diverged taxa has rarely been surveyed. We measured recombination rates in eight species that collectively represent several temporal scales of divergence within a single rodent family, Muridae. We used a cytological approach that enables in situ visualization of crossovers at meiosis to quantify recombination rates in multiple males from each rodent group. We uncovered large differences in genomic recombination rate between rodent species, which were independent of karyotypic variation. The divergence in genomic recombination rate that we document is not proportional to DNA sequence divergence, suggesting that recombination has evolved at variable rates along the murid phylogeny. Additionally, we document significant variation in genomic recombination rate both within and between subspecies of house mice. Recombination rates estimated in F(1) hybrids reveal evidence for sex-linked loci contributing to the evolution of recombination in house mice. Our results provide one of the first detailed portraits of genomic-scale recombination rate variation within a single mammalian family and demonstrate that the low recombination rates in laboratory mice and rats reflect a more general reduction in recombination rate across murid rodents.     

The influence of recombination on SNP diversity in chickens

The association between recombination rate and diversity, but not divergence is considered to be driven mainly by natural selection: fixation of positively selected variants and associated hitchhiking effects and/or background selection eliminating deleterious alleles. In the present study, we investigated the relationship between recombination rate, SNP diversity and interspecies divergence for 29 loci in chickens. We found that recombination rate is positively correlated with nucleotide diversity but is not correlated with interspecies divergence. It appears that variation in recombination rate explains over 30% of the variation in levels of diversity among 29 loci. Our data suggested that natural selection is a main factor in shaping SNP diversity in chickens. Since SNP diversity is significantly lower at Z-linked than at autosomal loci, we argued that genetic hitchhiking might be more important than background selection in producing the observed correlation.     

The beta -globin recombinational hotspot reduces the effects of strong selection around HbC, a recently arisen mutation providing resistance to malaria

Recombination is expected to reduce the effect of selection on the extent of linkage disequilibrium (LD), but the impact that recombinational hotspots have on sites linked to selected mutations has not been investigated. We empirically determine chromosomal linkage phase for 5.2 kb spanning the beta -globin gene and hotspot. We estimate that the HbC mutation, which is positively selected because of malaria, originated <5,000 years ago and that selection coefficients are 0.04-0.09. Despite strong selection and the recent origin of the HbC allele, recombination (crossing-over or gene conversion) is observed within 1 kb 5' of the selected site on more than one-third of the HbC chromosomes sampled. The rapid decay in LD upstream of the HbC allele demonstrates the large effect the ss-globin hotspot has in mitigating the effects of positive selection on linked variation.     

Genetic analysis of genome-scale recombination rate evolution in house mice

The rate of meiotic recombination varies markedly between species and among individuals. Classical genetic experiments demonstrated a heritable component to population variation in recombination rate, and specific sequence variants that contribute to recombination rate differences between individuals have recently been identified. Despite these advances, the genetic basis of species divergence in recombination rate remains unexplored. Using a cytological assay that allows direct in situ imaging of recombination events in spermatocytes, we report a large (∼30%) difference in global recombination rate between males of two closely related house mouse subspecies (Mus musculus musculus and M. m. castaneus). To characterize the genetic basis of this recombination rate divergence, we generated an F2 panel of inter-subspecific hybrid males (n = 276) from an intercross between wild-derived inbred strains CAST/EiJ (M. m. castaneus) and PWD/PhJ (M. m. musculus). We uncover considerable heritable variation for recombination rate among males from this mapping population. Much of the F2 variance for recombination rate and a substantial portion of the difference in recombination rate between the parental strains is explained by eight moderate- to large-effect quantitative trait loci, including two transgressive loci on the X chromosome. In contrast to the rapid evolution observed in males, female CAST/EiJ and PWD/PhJ animals show minimal divergence in recombination rate (∼5%). The existence of loci on the X chromosome suggests a genetic mechanism to explain this male-biased evolution. Our results provide an initial map of the genetic changes underlying subspecies differences in genome-scale recombination rate and underscore the power of the house mouse system for understanding the evolution of this trait.     

Evaluating human autosomal loci for sexually antagonistic viability selection in two large biobanks

Sex and sexual differentiation are pervasive across the tree of life. Because females and males often have substantially different functional requirements, we expect selection to differ between the sexes. Recent studies in diverse species, including humans, suggest that sexually antagonistic viability selection creates allele frequency differences between the sexes at many different loci. However, theory and population-level simulations indicate that sex-specific differences in viability would need to be very large to produce and maintain reported levels of between-sex allelic differentiation. We address this contradiction between theoretical predictions and empirical observations by evaluating evidence for sexually antagonistic viability selection on autosomal loci in humans using the largest cohort to date (UK Biobank, n = 487,999) along with a second large, independent cohort (BioVU, n = 93,864). We performed association tests between genetically ascertained sex and autosomal loci. Although we found dozens of genome-wide significant associations, none replicated across cohorts. Moreover, closer inspection revealed that all associations are likely due to cross-hybridization with sex chromosome regions during genotyping. We report loci with potential for mis-hybridization found on commonly used genotyping platforms that should be carefully considered in future genetic studies of sex-specific differences. Despite being well powered to detect allele frequency differences of up to 0.8% between the sexes, we do not detect clear evidence for this signature of sexually antagonistic viability selection on autosomal variation. These findings suggest a lack of strong ongoing sexually antagonistic viability selection acting on single locus autosomal variation in humans.