Recombination Hotspots Research Articles

Turning coldspots into hotspots: targeted recruitment of axis protein Hop1 stimulates meiotic recombination in Saccharomyces cerevisiae.

The DNA double-strand breaks that initiate meiotic recombination are formed in the context of the meiotic chromosome axis, which in Saccharomyces cerevisiae contains a meiosis-specific cohesin isoform and the meiosis-specific proteins Hop1 and Red1. Hop1 and Red1 are important for double-strand break formation; double-strand break levels are reduced in their absence and their levels, which vary along the lengths of chromosomes, are positively correlated with double-strand break levels. How axis protein levels influence double-strand break formation and recombination remains unclear. To address this question, we developed a novel approach that uses a bacterial ParB-parS partition system to recruit axis proteins at high levels to inserts at recombination coldspots where Hop1 and Red1 levels are normally low. Recruiting Hop1 markedly increased double-strand breaks and homologous recombination at target loci, to levels equivalent to those observed at endogenous recombination hotspots. This local increase in double-strand breaks did not require Red1 or the meiosis-specific cohesin component Rec8, indicating that, of the axis proteins, Hop1 is sufficient to promote double-strand break formation. However, while most crossovers at endogenous recombination hotspots are formed by the meiosis-specific MutLγ resolvase, crossovers that formed at an insert locus were only modestly reduced in the absence of MutLγ, regardless of whether or not Hop1 was recruited to that locus. Thus, while local Hop1 levels determine local double-strand break levels, the recombination pathways that repair these breaks can be determined by other factors, raising the intriguing possibility that different recombination pathways operate in different parts of the genome.

Expressed genes and their new alleles identification during fibre elongation reveal the genetic factors underlying improvements of fibre length in cotton.

SummaryInterspecific breeding in cotton takes advantage of genetic recombination among desirable genes from different parental lines. However, the expression new alleles (ENAs) from crossovers within genic regions and their significance in fibre length (FL) improvement are currently not understood. Here, we generated resequencing genomes of 191 interspecific backcross inbred lines derived from CRI36 (Gossypium hirsutum) × Hai7124 (Gossypium barbadense) and 277 dynamic fibre transcriptomes to identify the ENAs and extremely expressed genes (eGenes) potentially influencing FL, and uncovered the dynamic regulatory network of fibre elongation. Of 35 420 eGenes in developing fibres, 10 366 ENAs were identified and preferentially distributed in chromosomes subtelomeric regions. In total, 1056–1255 ENAs showed transgressive expression in fibres at 5–15 dpa (days post‐anthesis) of some BILs, 520 of which were located in FL‐quantitative trait locus (QTLs) and GhFLA9 (recombination allele) was identified with a larger effect for FL than GhFLA9 of CRI36 allele. Using ENAs as a type of markers, we identified three novel FL‐QTLs. Additionally, 456 extremely eGenes were identified that were preferentially distributed in recombination hotspots. Importantly, 34 of them were significantly associated with FL. Gene expression quantitative trait locus analysis identified 1286, 1089 and 1059 eGenes that were colocalized with the FL trait at 5, 10 and 15 dpa, respectively. Finally, we verified the Ghir_D10G011050 gene linked to fibre elongation by the CRISPR‐cas9 system. This study provides the first glimpse into the occurrence, distribution and expression of the developing fibres genes (especially ENAs) in an introgression population, and their possible biological significance in FL.

Divergent Viruses Discovered in Swine Alter the Understanding of Evolutionary History and Genetic Diversity of the Respirovirus Genus and Related Porcine Parainfluenza Viruses.

ABSTRACTParamyxoviridae is a rapidly growing family of viruses, whose potential for cross-species transmission makes it difficult to predict the harm of newly emerging viruses to humans and animals. To better understand their diversity, evolutionary history, and co-evolution with their hosts, we analyzed a collection of porcine parainfluenza virus (PPIV) genomes to reconstruct the species classification basis and evolutionary history of the Respirovirus genus. We sequenced 17 complete genomes of porcine respirovirus 1 (also known as porcine parainfluenza virus 1; PPIV-1), thereby nearly tripling the number of currently available PPIV-1 genomes. We found that PPIV-1 was widely prevalent in China with two divergent lineages, PPIV-1a and PPIV-1b. We further provided evidence that a new species, porcine parainfluenza virus 2 (PPIV-2), had recently emerged in China. Our results pointed to a need for revising the current species demarcation criteria of the Respirovirus genus. In addition, we used PPIV-1 as an example to explore recombination and diversity of the Respirovirus genus. Interestingly, we only detected heterosubtypic recombination events between PPIV-1a and PPIV-1b with no intrasubtypic recombination events. The recombination hotspots highlighted a diverse geography-dependent genome structure of paramyxovirus infecting swine in China. Furthermore, we found no evidence of co-evolution between respirovirus and its host, indicating frequent cross-species transmission. In summary, our analyses showed that swine can be infected with a broad range of respiroviruses and recombination may serve as an important evolutionary mechanism for the Respirovirus genus’ greater diversity in genome structure than previously anticipated.IMPORTANCE Livestock have emerged as critically underrecognized sources of paramyxovirus diversity, including pigs serving as the source of Nipah virus (NiV) and swine parainfluenza virus type 3, and goats and bovines harboring highly divergent viral lineages. Here, we identified a new species of Respirovirus genus named PPIV-2 in swine and proposed to revise the species demarcation criteria of the Respirovirus genus. We found heterosubtypic recombination events and high genetic diversity in PPIV-1. Further, we showed that genetic recombination may have occurred in the Respirovirus genus which may be associated with host range expansion. The continued expansion of Respirovirus genus diversity in livestock with relatively high human contact rates requires enhanced surveillance and ongoing evaluation of emerging cross-species transmission threats.

Cryptic Genes for Interbacterial Antagonism Distinguish Rickettsia Species Infecting Blacklegged Ticks From Other Rickettsia Pathogens.

BackgroundThe genus Rickettsia (Alphaproteobacteria: Rickettsiales) encompasses numerous obligate intracellular species with predominantly ciliate and arthropod hosts. Notable species are pathogens transmitted to mammals by blood-feeding arthropods. Mammalian pathogenicity evolved from basal, non-pathogenic host-associations; however, some non-pathogens are closely related to pathogens. One such species, Rickettsia buchneri, is prevalent in the blacklegged tick, Ixodes scapularis. While I. scapularis transmits several pathogens to humans, it does not transmit Rickettsia pathogens. We hypothesize that R. buchneri established a mutualism with I. scapularis, blocking tick superinfection with Rickettsia pathogens.MethodsTo improve estimates for assessing R. buchneri infection frequency in blacklegged tick populations, we used comparative genomics to identify an R. buchneri gene (REIS_1424) not present in other Rickettsia species present throughout the I. scapularis geographic range. Bioinformatic and phylogenomics approaches were employed to propose a function for the hypothetical protein (263 aa) encoded by REIS_1424.ResultsREIS_1424 has few analogs in other Rickettsiales genomes and greatest similarity to non-Proteobacteria proteins. This cohort of proteins varies greatly in size and domain composition, possessing characteristics of Recombination hotspot (Rhs) and contact dependent growth inhibition (CDI) toxins, with similarity limited to proximal C-termini (~145 aa). This domain was named CDI-like/Rhs-like C-terminal toxin (CRCT). As such proteins are often found as toxin-antidote (TA) modules, we interrogated REIS_1423 (151 aa) as a putative antidote. Indeed, REIS_1423 is similar to proteins encoded upstream of CRCT domain-containing proteins. Accordingly, we named these proteins CDI-like/Rhs-like C-terminal toxin antidotes (CRCA). R. buchneri expressed both REIS_1423 and REIS_1424 in tick cell culture, and PCR assays showed specificity for R. buchneri over other rickettsiae and utility for positive detection in three tick populations. Finally, phylogenomics analyses uncovered divergent CRCT/CRCA modules in varying states of conservation; however, only R. buchneri and related Tamurae/Ixodes Group rickettsiae carry complete TA modules.ConclusionWe hypothesize that Rickettsia CRCT/CRCA modules circulate in the Rickettsia mobile gene pool, arming rickettsiae for battle over arthropod colonization. While its functional significance remains to be tested, R. buchneri CRCT/CRCA serves as a marker to positively identify infection and begin deciphering the role this endosymbiont plays in the biology of the blacklegged tick.

Coexpression of MEIOTIC-TOPOISOMERASE VIB-dCas9 with guide RNAs specific to a recombination hotspot is insufficient to increase crossover frequency in Arabidopsis.

During meiosis, homologous chromosomes pair and recombine, which can result in reciprocal crossovers that increase genetic diversity. Crossovers are unevenly distributed along eukaryote chromosomes and show repression in heterochromatin and the centromeres. Within the chromosome arms, crossovers are often concentrated in hotspots, which are typically in the kilobase range. The uneven distribution of crossovers along chromosomes, together with their low number per meiosis, creates a limitation during crop breeding, where recombination can be beneficial. Therefore, targeting crossovers to specific genome locations has the potential to accelerate crop improvement. In plants, meiotic crossovers are initiated by DNA double-strand breaks that are catalyzed by SPO11 complexes, which consist of 2 catalytic (SPO11-1 and SPO11-2) and 2 noncatalytic subunits (MTOPVIB). We used the model plant Arabidopsis thaliana to coexpress an MTOPVIB-dCas9 fusion protein with guide RNAs specific to the 3a crossover hotspot. We observed that this was insufficient to significantly change meiotic crossover frequency or pattern within 3a. We discuss the implications of our findings for targeting meiotic recombination within plant genomes.

Deep learning identifies and quantifies recombination hotspot determinants.

MotivationRecombination is one of the essential genetic processes for sexually reproducing organisms, which can happen more frequently in some regions, called recombination hotspots. Although several factors, such as PRDM9 binding motifs, are known to be related to the hotspots, their contributions to the recombination hotspots have not been quantified, and other determinants are yet to be elucidated. Here, we propose a computational method, RHSNet, based on deep learning and signal processing, to identify and quantify the hotspot determinants in a purely data-driven manner, utilizing datasets from various studies, populations, sexes and species.ResultsRHSNet can significantly outperform other sequence-based methods on multiple datasets across different species, sexes and studies. In addition to being able to identify hotspot regions and the well-known determinants accurately, more importantly, RHSNet can quantify the determinants that contribute significantly to the recombination hotspot formation in the relation between PRDM9 binding motif, histone modification and GC content. Further cross-sex, cross-population and cross-species studies suggest that the proposed method has the generalization power and potential to identify and quantify the evolutionary determinant motifs.Availability and implementation https://github.com/frankchen121212/RHSNet.Supplementary information Supplementary data are available at Bioinformatics online.

Phenotypic characterization and genome analysis of a novel Salmonella Typhimurium phage having unique tail fiber genes

Salmonella enterica serovar Typhimurium is a foodborne pathogen causing occasional outbreaks of enteric infections in humans. Salmonella has one of the largest pools of temperate phages in its genome that possess evolutionary significance for pathogen. In this study, we characterized a novel temperate phage Salmonella phage BIS20 (BIS20) with unique tail fiber genes. It belongs to the subfamily Peduovirinae genus Eganvirus and infects Salmonella Typhimurium strain (SE-BS17; Acc. NO MZ503545) of poultry origin. Phage BIS20 was viable only at biological pH and temperature ranges (pH7 and 37 °C). Despite being temperate BIS20 significantly slowed down the growth of host strain for 24 h as compared to control (P < 0.009). Phage BIS20 features 29,477-base pair (bp) linear DNA genome with 53% GC content and encodes for 37 putative ORFs. These ORFs have mosaic arrangement as indicated by its ORF similarity to various phages and prophages in NCBI. Genome analysis indicates its similarity to Salmonella enterica serovar Senftenberg prophage (SEStP) sequence (Nucleotide similarity 87.7%) and Escherichia virus 186 (~ 82.4% nucleotide similarity). Capsid genes were conserved however those associated with tail fiber formation and assembly were unique to all members of genus Eganvirus. We found strong evidence of recombination hotspot in tail fiber gene. Our study identifies BIS20 as a new species of genus Eganvirus temperate phages as its maximum nucleotide similarity is 82.4% with any phage in NCBI. Our findings may contribute to understanding of origin of new temperate phages.

Recombination and mutation shape variations in the major histocompatibility complex

The major histocompatibility complex (MHC) is closely associated with numerous diseases, but its high degree of polymorphism complicates the discovery of disease-associated variants. In principle, recombination and de novo mutations are two critical factors responsible for MHC polymorphisms. However, direct evidence for this hypothesis is lacking. Here, we report the generation of fine-scale MHC recombination and de novo mutation maps of ∼5 Mb by deep sequencing (> 100×) of the MHC genome for 17 MHC recombination and 30 non-recombination Han Chinese families (a total of 190 individuals). Recombination hotspots and Han-specific breakpoints are located in close proximity at haplotype block boundaries. The average MHC de novo mutation rate is higher than the genome-wide de novo mutation rate, particularly in MHC recombinant individuals. Notably, mutation and recombination generated polymorphisms are located within and outside linkage disequilibrium regions of the MHC, respectively, and evolution of the MHC locus was mainly controlled by positive selection. These findings provide insights on the evolutionary causes of the MHC diversity and may facilitate the identification of disease-associated genetic variants.

Banana bunchy top virus genetic diversity in Pakistan and association of diversity with recombination in its genomes

Banana Bunchy top virus (BBTV) is a multipartite circular single strand DNA virus that belongs to genus Babuvirus and family Nanoviridae. It causes significant crop losses worldwide and also in Pakistan. BBTV is present in Pakistan since 1988 however, till now only few (about twenty only) sequence of genomic components have been reported from the country. To have insights into current genetic diversity in Pakistan fifty-seven genomic components including five complete genomes (comprises of DNA-R, -U3, -S, -M, -C and -N components) were sequenced in this study. The genetic diversity analysis of populations from Pakistan showed that DNA-R is highly conserved followed by DNA-N, whereas DNA-U3 is highly diverse with the most diverse Common Region Stem-loop (CR-SL) in BBTV genome, a functional region, which previously been reported to have undergone recombination in Pakistani population. A Maximum Likelihood (ML) phylogenetic analysis of entire genomes of isolates by using sequence of all the components concatenated together with the reported genomes around the world revealed deeper insights about the origin of the disease in Pakistan. A comparison of the genetic diversity of Pakistani and entire BBTV populations around the world indicates that there exists a correlation between genetic diversity and recombination. Population genetics analysis indicated that the degree of selection pressure differs depending on the area and genomic component. A detailed analysis of recombination across various components and functional regions suggested that recombination is closely associated with the functional parts of BBTV genome showing high genetic diversity. Both genetic diversity and recombination analyses suggest that the CR-SL is a recombination hotspot in all BBTV genomes and among the six components DNA-U3 is the only recombined component that has extensively undergone inter and intragenomic recombination. Diversity analysis of recombinant regions results on average one and half fold increase and, in some cases up to four-fold increase due to recombination. These results suggest that recombination is significantly contributing to the genetic diversity of BBTV populations around the world.

PRDM9-directed recombination hotspots depleted near meiotically transcribed genes.

PRDM9 drives recombination hotspots in some mammals, including mice and apes. Non-functional orthologs of PRDM9 are present in a wide variety of vertebrates, but why it is functionally maintained in some lineages is not clear. One possible explanation is that PRDM9 plays a role in ensuring that meiosis is successful. During meiosis, available DNA may be a limiting resource given the tight packaging of chromosomes and could lead to competition between two key processes: meiotic transcription and recombination. Here we explore this potential competition and the role that PRDM9 might play in their interaction. Leveraging existing mouse genomic data, we use resampling schemes that simulate shuffled features along the genome and models that account for the rarity of features in the genome, to test if PRDM9 influences interactions between recombination hotspots and meiotic transcription in a whole genome framework. We also explored possible DNA sequence motifs associated to clusters of hotspots not tied to transcription or PRDM9. We find evidence of competition between meiotic transcription and recombination, with PRDM9 appearing to relocate recombination to avoid said conflict. We also find that retrotransposons may be playing a role in directing hotspots in the absence of other factors.

Canis familiaris (Great Dane domestic dog)

Previous assembly gaps that are now resolved in the Great Dane genome are enriched for sequences with an extreme GC content and are located near transcription start sites and recombination hot spots. The identified gene models include 2151 genes whose start positions map to a gap in the CanFam3.1 assembly. Due to their sequence composition, many of these sequences are recalcitrant to analysis using Illumina short-read sequencing. Several additional canine assemblies are now available that have improved continuity and gene annotations.

Exploring the Natural Origins of SARS-CoV-2 in the Light of Recombination.

The lack of an identifiable intermediate host species for the proximal animal ancestor of SARS-CoV-2, and the large geographical distance between Wuhan and where the closest evolutionary related coronaviruses circulating in horseshoe bats (members of the Sarbecovirus subgenus) have been identified, is fueling speculation on the natural origins of SARS-CoV-2. We performed a comprehensive phylogenetic study on SARS-CoV-2 and all the related bat and pangolin sarbecoviruses sampled so far. Determining the likely recombination events reveals a highly reticulate evolutionary history within this group of coronaviruses. Distribution of the inferred recombination events is nonrandom with evidence that Spike, the main target for humoral immunity, is beside a recombination hotspot likely driving antigenic shift events in the ancestry of bat sarbecoviruses. Coupled with the geographic ranges of their hosts and the sampling locations, across southern China, and into Southeast Asia, we confirm that horseshoe bats, Rhinolophus, are the likely reservoir species for the SARS-CoV-2 progenitor. By tracing the recombinant sequence patterns, we conclude that there has been relatively recent geographic movement and cocirculation of these viruses’ ancestors, extending across their bat host ranges in China and Southeast Asia over the last 100 years. We confirm that a direct proximal ancestor to SARS-CoV-2 has not yet been sampled, since the closest known relatives collected in Yunnan shared a common ancestor with SARS-CoV-2 approximately 40 years ago. Our analysis highlights the need for dramatically more wildlife sampling to: 1) pinpoint the exact origins of SARS-CoV-2’s animal progenitor, 2) the intermediate species that facilitated transmission from bats to humans (if there is one), and 3) survey the extent of the diversity in the related sarbecoviruses’ phylogeny that present high risk for future spillovers.

Genome Plasticity of African Swine Fever Virus: Implications for Diagnostics and Live-Attenuated Vaccines.

African swine fever (ASF) is a highly contagious transboundary viral hemorrhagic disease of domestic and wild pigs presenting a significant threat to the global swine industry. Following its introduction in Caucasus, Georgia, in 2007, the genome of the genotype II of African swine fever virus (ASFV) strain Georgia-07 and its derivatives accumulated significant mutations, resulting in the emergence of genetic variants within short epidemiological timescales as it spreads and infects different hosts in diverse ecosystems, causing outbreaks in Europe, South Asia, South East Asia and the Caribbean. This suggests that ASFV, with a comparatively large and complex DNA genome, is susceptible to genetic mutations including deletions and that although the virus is environmentally stable, it is genetically unstable. This has implications for the development of vaccines and diagnostic tests for disease detection and surveillance. Analysis of the ASFV genome revealed recombination hotspots, which in double-stranded DNA (dsDNA) viruses represent key drivers of genetic diversity. The ability of pox virus, a dsDNA virus with a genome complexity similar to ASFV, regaining virulence following the deletion of a virulence gene via gene amplification, coupled with the recent emergence and spread of live-attenuated ASFV vaccine strains causing disease and death in pigs in China, raise legitimate concerns around the use of live-attenuated ASFV vaccines in non-endemic regions to control the potential introduction. Further research into the risk of using live-attenuated ASFV in non-endemic regions is highly needed.

Construction of Genetic Linkage Maps From a Hybrid Family of Large Yellow Croaker (Larimichthys crocea).

Consensus and sex-specific genetic linkage maps for large yellow croaker (Larimichthys crocea) were constructed using samples from an F1 family produced by crossing a Daiqu female and a Mindong male. A total of 20,147 single nucleotide polymorphisms (SNPs) by restriction site associated DNA sequencing were assigned to 24 linkage groups (LGs). The total length of the consensus map was 1757.4 centimorgan (cM) with an average marker interval of 0.09 cM. The total length of female and male linkage map was 1533.1 cM and 1279.2 cM, respectively. The average female-to-male map length ratio was 1.2 ± 0.23. Collapsed markers in the genetic maps were re-ordered according to their relative positions in the ASM435267v1 genome assembly to produce integrated genetic linkage maps with 9885 SNPs distributed across the 24 LGs. The recombination pattern of most LGs showed sigmoidal patterns of recombination, with higher recombination in the middle and suppressed recombination at both ends, which corresponds with the presence of sub-telocentric and acrocentric chromosomes in the species. The average recombination rate in the integrated female and male maps was respectively 3.55 cM/Mb and 3.05 cM/Mb. In most LGs, higher recombination rates were found in the integrated female map, compared to the male map, except in LG12, LG16, LG21, LG22, and LG24. Recombination rate profiles within each LG differed between the male and the female, with distinct regions indicating potential recombination hotspots. Separate quantitative trait loci (QTL) and association analyses for growth related traits in 6 months fish were performed, however, no significant QTL was detected. The study indicates that there may be genetic differences between the two strains, which may have implications for the application of DNA-information in the further breeding schemes.

RecBCD enzyme and Chi recombination hotspots as determinants of self vs. non-self: Myths and mechanisms.

Bacteria face a challenge when DNA enters their cells by transformation, mating, or phage infection. Should they treat this DNA as an invasive foreigner and destroy it, or consider it one of their own and potentially benefit from incorporating new genes or alleles to gain useful functions? It is frequently stated that the short nucleotide sequence Chi (5' GCTGGTGG 3'), a hotspot of homologous genetic recombination recognized by Escherichia coli's RecBCD helicase-nuclease, allows E. coli to distinguish its DNA (self) from any other DNA (non-self) and to destroy non-self DNA, and that Chi is "over-represented" in the E. coli genome. We show here that these latter statements (dogmas) are not supported by available evidence. We note Chi's wide-spread occurrence and activity in distantly related bacterial species and phages. We illustrate multiple, highly non-random features of the genomes of E. coli and coliphage P1 that account for Chi's high frequency and genomic position, leading us to propose that P1 selects for Chi's enhancement of recombination, whereas E. coli selects for the preferred codons in Chi. We discuss other, previously described mechanisms for self vs. non-self determination involving RecBCD and for RecBCD's destruction of DNA that cannot recombine, whether foreign or domestic, with or without Chi.

Chromosome-Wide Characterization of Intragenic Crossover in Shiitake Mushroom, Lentinula edodes.

Meiotic crossover plays a critical role in generating genetic variations and is a central component of breeding. However, our understanding of crossover in mushroom-forming fungi is limited. Here, in Lentinula edodes, we characterized the chromosome-wide intragenic crossovers, by utilizing the single-nucleotide polymorphisms (SNPs) datasets of an F1 haploid progeny. A total of 884 intragenic crossovers were identified in 110 single-spore isolates, the majority of which were closer to transcript start sites. About 71.5% of the intragenic crossovers were clustered into 65 crossover hotspots. A 10 bp motif (GCTCTCGAAA) was significantly enriched in the hotspot regions. Crossover frequencies around mating-type A (MAT-A) loci were enhanced and formed a hotspot in L. edodes. Genome-wide quantitative trait loci (QTLs) mapping identified sixteen crossover-QTLs, contributing 8.5–29.1% of variations. Most of the detected crossover-QTLs were co-located with crossover hotspots. Both cis- and trans-QTLs contributed to the nonuniformity of crossover along chromosomes. On chr2, we identified a QTL hotspot that regulated local, global crossover variation and crossover hotspot in L. edodes. These findings and observations provide a comprehensive view of the crossover landscape in L. edodes, and advance our understandings of conservation and diversity of meiotic recombination in mushroom-forming fungi.

Evolutionary Dynamics of Type 2 Porcine Reproductive and Respiratory Syndrome Virus by Whole-Genome Analysis.

Porcine reproductive and respiratory syndrome virus (PRRSV), an important pathogen in the swine industry, is a genetically highly diverse RNA virus. However, the phylogenetic and genomic recombination properties of this virus are not yet fully understood. In this study, we performed an integrated analysis of all available whole-genome sequences of type 2 PRRSV (n = 901) to reveal its evolutionary dynamics. The results showed that there were three distinct phylogenetic lineages of PRRSV in their distribution patterns. We identified that sublineage 2.7 (L2.7), associated with a NADC30 cluster, had the highest substitution rate and higher viral genetic diversity, and inter-lineage recombination is observed more frequently in L2.7 PRRSV compared to other sublineages. Most inter-lineage recombination events detected are observed between L2.7 PRRSVs (as major parents) and L3.4 (a JXA1-R-related cluster)/L3.7 (a WUH3-related cluster) PRRSVs (as minor parents). Moreover, the recombination hotspots are located in the structural protein gene ORF2 and ORF4, or in the non-structural protein gene nsp7. In addition, a GM2-related cluster, L3.2, shows inconsistent recombination modes compared to those of L2.7, suggesting that it may have undergone extensive and unique recombination in their evolutionary history. We also identified several amino acids under positive selection in GP2, GP4 and GP5, the major glycoproteins of PRRSV, showing the driving force behind adaptive evolution. Taken together, our results provide new insights into the evolutionary dynamics of PPRSV that contribute to our understanding of the critical factors involved in its evolution and guide future efforts to develop effective preventive measures against PRRSV.

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

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

WGS- versus ORF5-Based Typing of PRRSV: A Belgian Case Study.

Porcine reproductive and respiratory syndrome virus (PRRSV) is the causative agent of one of the most widespread and economically devastating diseases in the swine industry. Typing circulating PRRSV strains by means of sequencing is crucial for developing adequate control strategies. Most genetic studies only target the highly variable open reading frame (ORF) 5, for which an extensive database is available. In this study, we performed whole-genome sequencing (WGS) on a collection of 124 PRRSV-1 positive serum samples that were collected over a 5-year period (2015–2019) in Belgium. Our results show that (nearly) complete PRRSV genomes can be obtained directly from serum samples with a high success rate. Analysis of the coding regions confirmed the exceptionally high genetic diversity, even among Belgian PRRSV-1 strains. To gain more insight into the added value of WGS, we performed phylogenetic cluster analyses on separate ORF datasets as well as on a single, concatenated dataset (CDS) containing all ORFs. A comparison between the CDS and ORF clustering schemes revealed numerous discrepancies. To explain these differences, we performed a large-scale recombination analysis, which allowed us to identify a large number of potential recombination events that were scattered across the genome. As PRRSV does not contain typical recombination hot-spots, typing PRRSV strains based on a single ORF is not recommended. Although the typing accuracy can be improved by including multiple regions, our results show that the full genetic diversity among PRRSV strains can only be captured by analysing (nearly) complete genomes. Finally, we also identified several vaccine-derived recombinant strains, which once more raises the question of the safety of these vaccines.

Genetic structure and evolutionary diversity of mating-type (MAT) loci in\xa0Hypsizygus marmoreus

The mating compatibility in fungi is generally governed by genes located within a single or two unlinked mating type (MAT) loci. Hypsizygus marmoreus is an edible mushroom in the order Agaricales with a tetrapolar system, which contains two unlinked MAT loci-homeodomain (HD) transcription factor genes and pheromone/pheromone receptor genes (P/R). In this study, we analyzed the genetic structure and diversity of MAT loci in tetrapolar system of H. marmoreus through sequencing of 54 heterokaryon and 8 homokaryon strains. Although within the HD loci, the gene order was conserved, the gene contents were variable, and the HD loci haplotypes were further classified into four types. By analyzing the structure, phylogeny, and the HD transmissibility based on the progeny of these four HD mating-type loci types, we found that they were heritable and tightly linked at the HD loci. The P/R loci genes were found to comprise three pheromone receptors, three pheromones, and two pheromone receptor-like genes. Intra- and inter-specific phylogenetic analyses of pheromone receptors revealed that the STE3 genes were divided into three groups, and we thus theorize that they diverged before speciation. Comparative analysis of the MAT regions among 73 Basidiomycete species indicated that the diversity of HD and P/R loci in Agaricales and Boletales may contribute to mating compatibility. The number of HD genes were not correlated with the tetrapolar or bipolar systems. In H. marmoreus, the expression levels of these genes at HD and P/R loci of compatible strains were found higher than in those of homonuclear/homokaryotic strains, indicating that these mating genes acted as switches for mating processes. Further collinear analysis of HD loci in interspecific species found that HD loci contains conserved recombination hotspots showing major rearrangements in Coprinopsis cinerea and Schizophyllum commune, suggesting different mechanisms for evolution of physically linked MAT loci in these groups. It seems likely that gene rearrangements are common in Agaricales fungi around HD loci. Together, our study provides insights into the genomic basis of mating compatibility in H. marmoreus.