Recombination Hotspots Research Articles

Data Preprocessing Methods for Selective Sweep Detection using Convolutional Neural Networks.

The identification of positive selection has been framed as a classification task, with Convolutional Neural Networks (CNNs) already outperforming summary statistics and likelihood-based approaches in accuracy. Despite the prevalence of CNN-based methods that manipulate the pixels of images representing raw genomic data as a preprocessing step to improve classification accuracy, the efficacy of these pixel-rearrangement techniques remains inadequately examined, particularly in the presence of confounding factors like population bottlenecks, migration and recombination hotspots. We introduce a set of pixel rearrangement algorithms aimed at enhancing CNN classification accuracy in detecting selective sweeps. These algorithms are employed to assess the performance of four CNN models for selective sweep detection. Our findings illustrate that the judicious application of rearrangement algorithms notably enhances the overall classification accuracy of a CNN across various datasets simulating confounding factors. We observed that sorting the columns of the genomic matrices has higher on CNN performance than rearranging the sequences. To some extent, these rearrangement algorithms are more robust to misspecified demographic models compared with the utilization of the default preprocessing algorithm as suggested by the respective authors of each CNN architecture. We provide the data rearrangement algorithms as a distinct package available for download at: https://github.com/Zhaohq96/Genetic-data-rearrangement.

Characterization and distribution of de novo mutations in the zebra finch

Germline de novo mutations (DNMs) provide the raw material for evolution. The DNM rate varies considerably between species, sexes and chromosomes. Here, we identify DNMs in the zebra finch (Taeniopygia guttata) across 16 parent-offspring trios using two genome assemblies of different quality. Using an independent genotyping assay, we validate 82% of the 150 candidate DNMs. DNM rates are consistent between both assemblies, with estimates of 6.14 × 10–9 and 6.36 × 10–9 per site per generation. We observe a strong paternal bias in DNM rates (male-to-female ratio ɑ ≈ 4), but this bias is in transition mutations only, leading to a transition-to-transversion ratio of 3.18 and 3.57. Finally, we find that DNMs tend to be randomly distributed across chromosomes, not associated with recombination hotspots or genic regions. However, the sex chromosome chrZ shows a roughly fourfold increased DNM rate compared to autosomes, which is more than the expected increase due to chrZ spending two-thirds of its time in males. Overall, our results further enhance our understanding of DNMs in passerine songbirds.

Squash vein yellowing virus from California emerged in the Middle East via intragenic and intergeneric recombination events in the hypervariable potyvirus P1 and ipomovirus P1a genes.

We present the complete sequence of the genomic RNA of an isolate of squash vein yellowing virus (Ipomovirus cucurbitavenaflavi) from California (SqVYV-CA) and show it is a recombinant virus with a highly divergent 5' UTR and proximal P1a gene. The evolution of SqVYV-CA involved an intrageneric event between unknown potyviruses, related to isolates of papaya ringspot virus (Potyvirus papayanuli) from the Old World, and an intergeneric event between this recombinant potyvirus (minor parent) and an isolate of SqVYV from Israel (SqVYV-IL) (major parent). These events occurred in mixed infections and in the potyvirus P1 and ipomovirus P1a recombination hotspots and resulted in SqVYV-CA having a potyvirus 5' UTR and chimeric P1-P1a gene/protein and the remainder of the genome from SqVYV-IL. The SqVYV-CA chimeric P1-P1a gene is under positive selection, and the protein is intrinsically disordered and may localize to the nucleus via nuclear localization signals in the P1 part. The C-terminal SqVYV-IL P1a part also diverged but retained the conserved serine protease motif. Furthermore, substantial divergence in SqVYV isolates from the Middle East was associated with genetic drift and a long evolutionary history in this region. The finding that the host range and symptomatology in cucurbits of SqVYV-CA is similar to those of SqVYV from Florida and SqVYV-IL, indicated that the recombinant part of the genome had no obvious effect on the virus-host interaction. A divergent part of the P1 sequence of the SqVYV-CA P1-P1a gene was used to develop a primer pair and RT-PCR test for specific detection of SqVYV-CA. This test was used to detect spread of SqVYV-CA to a new production area of California in 2021 and 2022. Together, these results demonstrate (i) a high level of genetic diversity exists among isolates of SqVYV and involved intra- and intergeneric recombination and genetic drift (mutation), (ii) evidence that SqVYV originated in the Middle East and that there were independent introductions into the New World and (iii) the remarkable genetic flexibility of the 5' proximal genes of these viruses.

Feline bocaviruses found in Thailand have undergone genetic recombination for their evolutions

Feline bocaviruses (FBoVs) have been discovered for a decade and are often detected in feces, possibly associated with diarrhea in cats. Studies on FBoV evolution remain limited and have mainly focused on prevalence and genetic characterization. Although genetic recombination serves as a potential mechanism in bocavirus evolution, research on this process for FBoVs has been scarce. In this study, we characterized 19 complete coding sequences of FBoVs obtained from Thai cats, revealing that FBoV-1, -2, and -3 were endemic in Thailand. Genetic characterizations showed that most Thai FBoVs were closely related to previously detected strains in Thailand and China. Recombination analyses indicated intragenic, intraspecies recombination in all FBoV species, with recombination breakpoints commonly found in the NP1 and VP1/2 genes, highlighting these genes may be hotspots for FBoV recombination. However, no interspecies recombination was detected. Selective pressure analysis of various FBoV genes revealed that these viruses underwent purifying selection. Although the VP1/2 gene of all FBoV species was under the strongest negative selection pressure, positive selection sites were only found in FBoV-1 and FBoV-3. This study is the first to identify natural recombination in FBoV-2 and FBoV-3 and provides evidence that genetic recombination is a potential driver of FBoV evolutions. Additionally, this study offers up-to-date information on the genetic characteristics, evolutionary dynamics, and selective pressure status of FBoVs, which should be continuously monitored.

Deep longitudinal lower respiratory tract microbiome profiling reveals genome-resolved functional and evolutionary dynamics in critical illness

The lower respiratory tract (LRT) microbiome impacts human health, especially among critically ill patients. However, comprehensive characterizations of the LRT microbiome remain challenging due to low microbial mass and host contamination. We develop a chelex100-based low-biomass microbial-enrichment method (CMEM) that enables deep metagenomic profiling of LRT samples to recover near-complete microbial genomes. We apply the method to 453 longitudinal LRT samples from 157 intensive care unit (ICU) patients in three geographically distant hospitals. We recover 120 high-quality metagenome-assembled genomes (MAGs) and associated plasmids without culturing. We detect divergent longitudinal microbiome dynamics and hospital-specific dominant opportunistic pathogens and resistomes in pneumonia patients. Diagnosed pneumonia and the ICU stay duration were associated with the abundance of specific antibiotic-resistance genes (ARGs). Moreover, CMEM can serve as a robust tool for genome-resolved analyses. MAG-based analyses reveal strain-specific resistome and virulome among opportunistic pathogen strains. Evolutionary analyses discover increased mobilome in prevailing opportunistic pathogens, highly conserved plasmids, and new recombination hotspots associated with conjugative elements and prophages. Integrative analysis with epidemiological data reveals frequent putative inter-patient strain transmissions in ICUs. In summary, we present a genome-resolved functional, transmission, and evolutionary landscape of the LRT microbiota in critically ill patients.

Multiple independent losses of crossover interference during yeast evolutionary history.

Meiotic recombination is essential for the accurate chromosome segregation and the generation of genetic diversity through crossover and gene conversion events. Although this process has been studied extensively in a few selected model species, understanding how its properties vary across species remains limited. For instance, the ancestral ZMM pathway that generates interference-dependent crossovers has undergone multiple losses throughout evolution, suggesting variations in the regulation of crossover formation. In this context, we first characterized the meiotic recombination landscape and properties of the Kluyveromyces lactis budding yeast. We then conducted a comprehensive analysis of 29,151 recombination events (19, 212 COs and 9, 939 NCOs) spanning 577 meioses in the five budding yeast species Saccharomyces cerevisiae, Saccharomyces paradoxus, Lachancea kluyveri, Lachancea waltii and K. lactis. Eventually, we found that the Saccharomyces yeasts displayed higher recombination rates compared to the non-Saccharomyces yeasts. In addition, bona fide crossover interference and associated crossover homeostasis were detected in the Saccharomyces species only, adding L. kluyveri and K. lactis to the list of budding yeast species that lost crossover interference. Finally, recombination hotspots, although highly conserved within the Saccharomyces yeasts are not conserved beyond the Saccharomyces genus. Overall, these results highlight great variability in the recombination landscape and properties through budding yeasts evolution.

Identification and functional analysis of recent IS transposition events in rhizobia

Rhizobia are alpha- and beta- Proteobacteria that, through the establishment of symbiotic interactions with leguminous plants, are able to fix atmospheric nitrogen as ammonium. The successful establishment of a symbiotic interaction is highly dependent on the availability of nitrogen sources in the soil, and on the specific rhizobia strain. Insertion sequences (ISs) are simple transposable genetic elements that can move to different locations within the host genome and are known to play an important evolutionary role, contributing to genome plasticity by acting as recombination hot-spots, and disrupting coding and regulatory sequences. Disruption of coding sequences may have occurred either in a common ancestor of the species or more recently. By means of ISComapare, we identified Differentially Located ISs (DLISs) in nearly related rhizobial strains of the genera Bradyrhizobium, Mesorhizobium, Rhizobium and Sinorhizobium. Our results revealed that recent IS transposition could have a role in adaptation by enabling the activation and inactivation of genes that could dynamically affect the competition and survival of rhizobia in the rhizosphere.Graphical

Haplotype-Aware Detection of SERPINA1 Variants by Nanopore Sequencing

Alpha-1 antitrypsin (AAT) is an acute-phase reactant with immunomodulatory properties that mainly inhibits neutrophil elastase. Low serum levels cause AAT deficiency (AATD), an underdiagnosed condition that predisposes to pulmonary and hepatic diseases. The SERPINA1 gene, which encodes AAT, contains more than 500 variants. PI*Z and PI*S alleles are the most diagnosed causes of AATD, but the role of the SERPINA1 haplotypes in AAT function remains unknown. SERPINA1 gene was PCR amplified from 94 asthma patients, using primers with tails for indexing. Sequencing libraries were loaded into a MinION-Mk1C, and MinKNOW was used for basecalling and demultiplexing. Nanofilt and Minimap2 were used for filtering and mapping/alignment. Variant calling/phasing were performed with PEPPER-Margin-DeepVariant. SERPINA1 gene was 100% covered for all samples, with a minimum sequencing depth of 500X. 75 single nucleotide variants and 4 indels were detected, with 45 and 2 of them highly polymorphic (MAF>0.1), respectively. Nine of the SNVs showed differences in allele frequencies when compared with the overall Spanish population. More than 90% of heterozygous SNVs were phased, yielding 91 and 58 different haplotypes for each SERPINA1 amplified region. Haplotype-based Linkage Disequilibrium (LD) analysis suggests that a recombination hotspot could generate variation in the SERPINA1 gene. The proposed workflow enables haplotype-aware genotyping of the SERPINA1 gene by nanopore sequencing, which will allow the development of novel AATD diagnostic strategies.

Novel crossover and recombination hotspots massively spread across primate genomes

BackgroundThe recombination landscape and subsequent natural selection have vast consequences forevolution and speciation. However, most of the crossover and recombination hotspots are yet to be discovered. We previously reported the relevance of C and G trinucleotide two-repeat units (CG-TTUs) in crossovers and recombination.MethodsOn a genome-wide scale, here we mapped all combinations of A and T trinucleotide two-repeat units (AT-TTUs) in human, consisting of AATAAT, ATAATA, ATTATT, TTATTA, TATTAT, and TAATAA. We also compared a number of the colonies formed by the AT-TTUs (distance between consecutive AT-TTUs < 500 bp) in several other primates and mouse.ResultsWe found that the majority of the AT-TTUs (> 96%) resided in approximately 1.4 million colonies, spread throughout the human genome. In comparison to the CG-TTU colonies, the AT-TTU colonies were significantly more abundant and larger in size. Pure units and overlapping units of the pure units were readily detectable in the same colonies, signifying that the units were the sites of unequal crossover. We discovered dynamic sharedness of several of the colonies across the primate species studied, which mainly reached maximum complexity and size in human.ConclusionsWe report novel crossover and recombination hotspots of the finest molecular resolution, massively spread and shared across the genomes of human and several other primates. With respect to crossovers and recombination, these genomes are far more dynamic than previously envisioned.

Genomic Underpinnings of Cytoplasmic Incompatibility: CIF Gene-Neighborhood Diversification Through Extensive Lateral Transfers and Recombination in Wolbachia.

Cytoplasmic incompatibility (CI), a non-Mendelian genetic phenomenon, involves the manipulation of host reproduction by Wolbachia, a maternally transmitted alphaproteobacterium. The underlying mechanism is centered around the CI Factor (CIF) system governed by two genes, cifA and cifB, where cifB induces embryonic lethality, and cifA counteracts it. Recent investigations have unveiled intriguing facets of this system, including diverse cifB variants, prophage association in specific strains, copy number variation, and rapid component divergence, hinting at a complex evolutionary history. We utilized comparative genomics to systematically classify CIF systems, analyze their locus structure and domain architectures, and reconstruct their diversification and evolutionary trajectories. Our new classification identifies ten distinct CIF types, featuring not just versions present in Wolbachia, but also other intracellular bacteria, and eukaryotic hosts. Significantly, our analysis of CIF loci reveals remarkable variability in gene composition and organization, encompassing an array of diverse endonucleases, variable toxin domains, deubiquitinating peptidases (DUBs), prophages, and transposons. We present compelling evidence that the components within the loci have been diversifying their sequences and domain architectures through extensive, independent lateral transfers and interlocus recombination involving gene conversion. The association with diverse transposons and prophages, coupled with selective pressures from host immunity, likely underpins the emergence of CIF loci as recombination hotspots. Our investigation also posits the origin of CifB-REase domains from mobile elements akin to CR (Crinkler-RHS-type) effectors and Tribolium Medea1 factor, which is linked to another non-Mendelian genetic phenomenon. This comprehensive genomic analysis offers novel insights into the molecular evolution and genomic foundations of Wolbachia-mediated host reproductive control.

Deciphering the evolutionary landscape of severe fever with thrombocytopenia syndrome virus across East Asia.

Severe fever with thrombocytopenia syndrome virus (SFTSV) poses a significant public health challenge in East Asia, necessitating a deeper understanding of its evolutionary dynamics to effectively manage its spread and pathogenicity. This study provides a comprehensive analysis of the genetic diversity, recombination patterns, and selection pressures across the SFTSV genome, utilizing an extensive dataset of 2041 sequences from various hosts and regions up to November 2023. Employing maximum likelihood and Bayesian evolutionary analysis by sampling trees (BEAST), we elucidated the phylogenetic relationships among nine distinct SFTSV genotypes (A, B1, B2, B3, B4, C, D, E, and F), revealing intricate patterns of viral evolution and genotype distribution across China, South Korea, and Japan. Furthermore, our analysis identified 34 potential reassortments, underscoring a dynamic genetic interplay among SFTSV strains. Genetic recombination was observed most frequently in the large segment and least in the small segment, with notable recombination hotspots characterized by stem-loop hairpin structures, indicative of a structural propensity for genetic recombination. Additionally, selection pressure analysis on critical viral genes indicated a predominant trend of negative selection, with specific sites within the RNA-dependent RNA polymerase and glycoprotein genes showing positive selection. These sites suggest evolutionary adaptations to host immune responses and environmental pressures. This study sheds light on the intricate evolutionary mechanisms shaping SFTSV, offering insights into its adaptive strategies and potential implications for vaccine development and therapeutic interventions.

Recombination of Porcine Reproductive and Respiratory Syndrome Virus: Features, Possible Mechanisms, and Future Directions.

Recombination is a pervasive phenomenon in RNA viruses and an important strategy for accelerating the evolution of RNA virus populations. Recombination in the porcine reproductive and respiratory syndrome virus (PRRSV) was first reported in 1999, and many case reports have been published in recent years. In this review, all the existing reports on PRRSV recombination events were collected, and the genotypes, parental strains, and locations of the recombination breakpoints have been summarized and analyzed. The results showed that the recombination pattern constantly changes; whether inter- or intra-lineage recombination, the recombination hotspots vary in different recombination patterns. The virulence of recombinant PRRSVs was higher than that of the parental strains, and the emergence of virulence reversion was caused by recombination after using MLV vaccines. This could be attributed to the enhanced adaptability of recombinant PRRSV for entry and replication, facilitating their rapid propagation. The aim of this paper was to identify common features of recombinant PRRSV strains, reduce the recombination risk, and provide a foundation for future research into the mechanism of PRRSV recombination.

Homologous Recombination Shapes the Architecture and Evolution of Bacterial Genomes.

Homologous recombination is a key evolutionary force that varies considerably across bacterial species. However, how the landscape of homologous recombination varies across genes and within individual genomes has only been studied in a few species. Here, we used Approximate Bayesian Computation to estimate the recombination rate along the genomes of 145 bacterial species. Our results show that homologous recombination varies greatly along bacterial genomes and shapes many aspects of genome architecture and evolution. The genomic landscape of recombination presents several key signatures: rates are highest near the origin of replication in most species, patterns of recombination generally appear symmetrical in both replichores (i.e. replicational halves of circular chromosomes) and most species have genomic hotpots of recombination. Furthermore, many closely related species share conserved landscapes of recombination across orthologs indicating that recombination landscapes are conserved over significant evolutionary distances. We show evidence that recombination drives the evolution of GC-content through increasing the effectiveness of selection and not through biased gene conversion, thereby contributing to an ongoing debate. Finally, we demonstrate that the rate of recombination varies across gene function and that many hotspots of recombination are associated with adaptive and mobile regions often encoding genes involved in pathogenicity.

The Prediction of Recombination Hotspot Based on Automated Machine Learning

Meiotic recombination plays a pivotal role in genetic evolution. Genetic variation induced by recombination is a crucial factor in generating biodiversity and a driving force for evolution. At present, the development of recombination hotspot prediction methods has encountered challenges related to insufficient feature extraction and limited generalization capabilities. This paper focused on the research of recombination hotspot prediction methods. We explored deep learning-based recombination hotspot prediction and scrutinized the shortcomings of prevalent models in addressing the challenge of recombination hotspot prediction. To addressing these deficiencies, an automated machine learning approach was utilized to construct recombination hotspot prediction model. The model combined sequence information with physicochemical properties by employing TF-IDF-Kmer and DNA composition components to acquire more effective feature data. Experimental results validate the effectiveness of the feature extraction method and automated machine learning technology used in this study. The final model was validated on three distinct datasets and yielded accuracy rates of 97.14%, 79.71%, and 98.73%, surpassing the current leading models by 2%, 2.56%, and 4%, respectively. In addition, we incorporated tools such as SHAP and AutoGluon to analyze the interpretability of black-box models, delved into the impact of individual features on the results, and investigated the reasons behind misclassification of samples. Finally, an application of recombination hotspot prediction was established to facilitate easy access to necessary information and tools for researchers. The research outcomes of this paper underscore the enormous potential of automated machine learning methods in gene sequence prediction.

Corrigendum: The evolutionary turnover of recombination hot spots contributes to speciation in mice.

Corrigendum: The evolutionary turnover of recombination hot spots contributes to speciation in mice.

Increased Positive Selection in Highly Recombining Genes Does not Necessarily Reflect an Evolutionary Advantage of Recombination.

It is commonly thought that the long-term advantage of meiotic recombination is to dissipate genetic linkage, allowing natural selection to act independently on different loci. It is thus theoretically expected that genes with higher recombination rates evolve under more effective selection. On the other hand, recombination is often associated with GC-biased gene conversion (gBGC), which theoretically interferes with selection by promoting the fixation of deleterious GC alleles. To test these predictions, several studies assessed whether selection was more effective in highly recombining genes (due to dissipation of genetic linkage) or less effective (due to gBGC), assuming a fixed distribution of fitness effects (DFE) for all genes. In this study, I directly derive the DFE from a gene's evolutionary history (shaped by mutation, selection, drift, and gBGC) under empirical fitness landscapes. I show that genes that have experienced high levels of gBGC are less fit and thus have more opportunities for beneficial mutations. Only a small decrease in the genome-wide intensity of gBGC leads to the fixation of these beneficial mutations, particularly in highly recombining genes. This results in increased positive selection in highly recombining genes that is not caused by more effective selection. Additionally, I show that the death of a recombination hotspot can lead to a higher dN/dS than its birth, but with substitution patterns biased towards AT, and only at selected positions. This shows that controlling for a substitution bias towards GC is therefore not sufficient to rule out the contribution of gBGC to signatures of accelerated evolution. Finally, although gBGC does not affect the fixation probability of GC-conservative mutations, I show that by altering the DFE, gBGC can also significantly affect nonsynonymous GC-conservative substitution patterns.

High prevalence of PRDM9-independent recombination hotspots in placental mammals

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

Effects of demographic history on recombination hotspots in soybean.

Recombination is the primary mechanism underlying genetic improvement in populations and allows plant breeders to create new allelic combinations for agronomic improvement. Soybean [Glycine max (L.) Merr.] has gone through multiple genetic bottlenecks that have significantly affected its genetic diversity, linkage disequilibrium, and altered allele frequencies. To investigate the impact of genetic bottlenecks on recombination hotspots in soybeans, historical recombination was studied in three soybean populations. The populations were wild soybean [Glycine soja (Sieb. and Zucc.)], landraces, and North American elite soybean cultivars that have been genotyped with the SoySNP50K BeadChip. While each population after a genetic bottleneck had an increased average haplotype block size, they did not have a significant difference in the number of hotspots between each population. Instead, the increase in observed haplotype block size is likely due to an elimination of individuals that contained historical recombination at hotspots which decreased the observed rate of recombination for the hotspot after each genetic bottleneck. Conversely, heterochromatic DNA which has an increased haplotype block size compared to euchromatic DNA had a significantly different number of hotspots but not a significant difference in the average hotspot recombination rate. Previously identified genomic motifs associated with hotspots were also associated with hotspots found in the historical populations suggesting a common mechanism. This characterization of historical recombination hotspots in soybeans provides further insights into the effect genetic bottlenecks and selection have on recombination hotspots.

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

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

A fine-scale genetic map of the Japanese population.

Genetic maps are fundamental resources for linkage and association studies. A fine-scale genetic map can be constructed by inferring historical recombination events from the genome-wide structure of linkage disequilibrium-a non-random association of alleles among loci-by using population-scale sequencing data. We constructed a fine-scale genetic map and identified recombination hotspots from 10 092 551 bi-allelic high-quality autosomal markers segregating among 150 unrelated Japanese individuals whose genotypes were determined by high-coverage (30×) whole-genome sequencing, and the genotype quality was carefully controlled by using their parents' and offspring's genotypes. The pedigree information was also utilized for haplotype phasing. The resulting genome-wide recombination rate profiles were concordant with those of the worldwide population on a broad scale, and the resolution was much improved. We identified 9487 recombination hotspots and confirmed the enrichment of previously known motifs in the hotspots. Moreover, we demonstrated that the Japanese genetic map improved the haplotype phasing and genotype imputation accuracy for the Japanese population. The construction of a population-specific genetic map will help make genetics research more accurate.