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.

BACKGROUND: Vancomycin is frequently used as a last line of defence against infections due to multidrug-resistant Staphylococcus aureus (S. aureus). A recent finding described the acquisition of vancomycin-resistant S. aureus strains by the integration of an enterococcal plasmid containing the vanA operon into the S. aureus chromosome via homologous recombination involving a specific integration site called locus L2.
METHODS: To characterise all mechanisms of acquisition of vanA, this study analysed the 15 706 S. aureus genomes to look for vanA and described its genetic environment.
RESULTS: A complete vanA operon was found in 25 S. aureus strains isolated from 12 patients, including nine co-isolated with vancomycin-resistant Enterococcus strains. VanA was found within transposon Tn1546-like elements on 17 plasmids and eight chromosomes. VanA might be acquired through conjugation of enterococcal and staphylococcal plasmids, transposition of Tn1546 carrying vanA and plasmid integration into the chromosome. Further, L2 was detected in 2087 genomes (13.3%) of S. aureus strains across different continents. Six potential chromosomal hotspots for integration of the entire vanA-containing enterococcal plasmid were identified by homologous recombination via L2.
CONCLUSIONS: These findings suggest that the recently described scenario in a New York patient could be reproduced anywhere. Surveillance of this possibility is mandatory, especially in patients with vancomycin-resistant Enterococcus infection or colonisation.

BACKGROUND: GGC and GCC short tandem repeats (STRs) are of various evolutionary, biological, and pathological implications. However, the fundamental two-repeats (dyads) of these STRs are widely unexplored.
RESULTS: On a genome-wide scale, we mapped (GGC)2 and (GCC)2 dyads in human, and found monumental colonies (distance between each dyad < 500 bp) of extraordinary density, and in some instances periodicity. The largest (GCC)2 and (GGC)2 colonies were intergenic, homogeneous, and human-specific, consisting of 219 (GCC)2 on chromosome 2 (probability < 1.545E-219) and 70 (GGC)2 on chromosome 9 (probability = 1.809E-148). We also found that several colonies were shared in other great apes, and directionally increased in density and complexity in human, such as a colony of 99 (GCC)2 on chromosome 20, that specifically expanded in great apes, and reached maximum complexity in human (probability 1.545E-220). Numerous other colonies of evolutionary relevance in human were detected in other largely overlooked regions of the genome, such as chromosome Y and pseudogenes. Several of the genes containing or nearest to those colonies were divergently expressed in human.
CONCLUSIONS: In conclusion, (GCC)2 and (GGC)2 form unprecedented genomic colonies that coincide with the evolution of human and other great apes. The extent of the genomic rearrangements leading to those colonies support overlooked recombination hotspots, shared across great apes. The identified colonies deserve to be studied in mechanistic, evolutionary, and functional platforms.

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

The genetic variability of SARS-CoV-2 (genus Betacoronavirus, family Coronaviridae) has been scrutinized since its first detection in December 2019. Although the role of structural variants, particularly deletions, in virus evolution is little explored, these genome changes are extremely frequent. They are associated with relevant processes, including immune escape and attenuation. Deletions commonly occur in accessory ORFs and might even lead to the complete loss of one or more ORFs. This scenario poses an interesting question about the origin and spreading of extreme structural rearrangements that persist without compromising virus viability. Here, we analyze the genome of SARS-CoV-2 in late 2021 in Uruguay and identify a Delta lineage (AY.20) that experienced a large deletion (872 nucleotides according to the reference Wuhan strain) that removes the 7a, 7b, and 8 ORFs. Deleted viruses coexist with wild-type (without deletion) AY.20 and AY.43 strains. The Uruguayan deletion is like those identified in Delta strains from Poland and Japan but occurs in a different Delta clade. Besides providing proof of the circulation of this large deletion in America, we infer that the 872-deletion arises by the consecutive occurrence of a 6-nucleotide deletion, characteristic of delta strains, and an 866-nucleotide deletion that arose independently in the AY.20 Uruguayan lineage. The largest deletion occurs adjacent to transcription regulatory sequences needed to synthesize the nested set of subgenomic mRNAs that serve as templates for transcription. Our findings support the role of transcription sequences as a hotspot for copy-choice recombination and highlight the remarkable dynamic of SARS-CoV-2 genomes.

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

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

A salmon louse (Lepeophtheirus salmonis salmonis) genetic linkage map was constructed to serve as a genomic resource for future investigations into the biology of this important marine parasitic copepod species, and to provide insights into the inheritance patterns of genetic markers in this species. SNP genotyping of 8 families confirmed the presence of 15 linkage groups based upon the assignment of 93,773 markers. Progeny sample size weight adjusted map sizes in males (with the exception of SL12 and SL15) ranged in size from 96.50 cM (SL11) to 134.61 cM (SL06), and total combined map steps or bins ranged from 143 (SL09) to 203 (SL13). The SL12 male map was the smallest linkage group with a weight-averaged size of 3.05 cM with 6 recombination bins. Male:female specific recombination rate differences are 10.49:1 and represent one of the largest reported sex-specific differences for any animal species. Recombination ratio differences (M:F) ranged from 1.0 (SL12) to 29:1 (SL15). The number of markers exhibiting normal Mendelian segregation within the sex linkage group SL15 was extremely low (N = 80) in comparison to other linkage groups genotyped [range: 1459 (SL12)-10206 markers (SL05)]. Re-evaluation of Mendelian inheritance patterns of markers unassigned to any mapping parent according to hemizygous segregation patterns (models presented) identified matches for many of these markers to hemizygous patterns. The greatest proportion of these markers assigned to SL15 (N increased to 574). Inclusion of the hemizygous markers revised SL15 sex-specific recombination rate differences to 28:1. Recombination hot- and coldspots were identified across all linkage groups with all linkage groups possessing multiple peaks. Nine of 13 linkage groups evaluated possessed adjacent domains with hot-coldspot transitional zones. The most common pattern was for one end of the linkage to show elevated recombination in addition to internal regions. For SL01 and SL06, however, a terminal region with high recombination was not evident while a central domain possessing extremely high-recombination levels was present. High levels of recombination were weakly coupled to higher levels of SNP variation within domains, but this association was very strong for the central domains of SL01 and SL06. From the pooled paternal half-sib lots (several virgin females placed with 1 male), only 1 or two surviving family lots were obtained. Surviving families possessed parents where both the male and female possessed either inherently low or high recombination rates. This study provides insight into the organization of the sea louse genome, and describes large differences in recombination rate that exist among individuals of the same sex, and between the sexes. These differences in recombination rate may be coupled to the capabilities of this species to adapt to environmental and pharmaceutical treatments, given that family survivorship appears to be enhanced when parents have similar recombination levels.

The programmed formation of hundreds of DNA double-strand breaks (DSBs) is essential for proper meiosis and fertility. In mice and humans, the location of these breaks is determined by the meiosis-specific protein PRDM9, through the DNA-binding specificity of its zinc-finger domain. PRDM9 also has methyltransferase activity. Here, we show that this activity is required for H3K4me3 and H3K36me3 deposition and for DSB formation at PRDM9-binding sites. By analyzing mice that express two PRDM9 variants with distinct DNA-binding specificities, we show that each variant generates its own set of H3K4me3 marks independently from the other variant. Altogether, we reveal several basic principles of PRDM9-dependent DSB site determination, in which an excess of sites are designated through PRDM9 binding and subsequent histone methylation, from which a subset is selected for DSB formation.

Traditional plant breeding relies on meiotic recombination for mixing of parental alleles to create novel allele combinations. Detailed analysis of recombination patterns in model organisms shows that recombination is tightly regulated within the genome, but frequencies vary extensively along chromosomes. Despite being a model organism for fruit developmental studies, high-resolution recombination patterns are lacking in tomato. In this study, we developed a novel methodology to use low-coverage resequencing to identify genome-wide recombination patterns and applied this methodology on 60 tomato Recombinant Inbred Lines (RILs). Our methodology identifies polymorphic markers from the low-coverage resequencing population data and utilizes the same data to locate the recombination breakpoints in individuals by using a variable sliding window. We identified 1,445 recombination sites comprising 112 recombination prone regions enriched for AT-rich DNA motifs. Furthermore, the recombination prone regions in tomato preferably occurred in gene promoters over intergenic regions, an observation consistent with Arabidopsis thaliana, Zea mays and Mimulus guttatus. Overall, our cost effective method and findings enhance the understanding of meiotic recombination in tomato and suggest evolutionarily conserved recombination associated genomic features.