High‐resolution mapping of the recombination landscape of the phytopathogen Fusarium graminearum suggests two‐speed genome evolution

Recombination is a major evolutionary force, increasing genetic diversity and permitting efficient coevolution of fungal pathogen(s) with their host(s). The ascomycete Fusarium graminearum is a devastating pathogen of cereal crops, and can contaminate food and feed with harmful mycotoxins. Previous studies have suggested a high adaptive potential of this pathogen, illustrated by an increase in pathogenicity and resistance to fungicides. In this study, we provide the first detailed picture of the crossover events occurring during meiosis and discuss the role of recombination in pathogen evolution. An experimental recombinant population (n = 88) was created and genotyped using 1306 polymorphic markers obtained from restriction site‐associated DNA sequencing (RAD‐seq) and aligned to the reference genome. The construction of a high‐density linkage map, anchoring 99% of the total length of the reference genome, allowed the identification of 1451 putative crossovers, positioned at a median resolution of 24 kb. The majority of crossovers (87.2%) occurred in a relatively small portion of the genome (30%). All chromosomes demonstrated recombination‐active sections, which had a near 15‐fold higher crossover rate than non‐active recombinant sections. The recombination rate showed a strong positive correlation with nucleotide diversity, and recombination‐active regions were enriched for genes with a putative role in host–pathogen interaction, as well as putative diversifying genes. Our results confirm the preliminary analysis observed in other F. graminearum strains and suggest a conserved ‘two‐speed’ recombination landscape. The consequences with regard to the evolutionary potential of this major fungal pathogen are also discussed.     

Whole‐genome patterns of linkage disequilibrium across flycatcher populations clarify the causes and consequences of fine‐scale recombination rate variation in birds

Recombination rate is heterogeneous across the genome of various species and so are genetic diversity and differentiation as a consequence of linked selection. However, we still lack a clear picture of the underlying mechanisms for regulating recombination. Here we estimated fine‐scale population recombination rate based on the patterns of linkage disequilibrium across the genomes of multiple populations of two closely related flycatcher species (Ficedula albicollis and F. hypoleuca). This revealed an overall conservation of the recombination landscape between these species at the scale of 200 kb, but we also identified differences in the local rate of recombination despite their recent divergence (<1 million years). Genetic diversity and differentiation were associated with recombination rate in a lineage‐specific manner, indicating differences in the extent of linked selection between species. We detected 400–3,085 recombination hotspots per population. Location of hotspots was conserved between species, but the intensity of hotspot activity varied between species. Recombination hotspots were primarily associated with CpG islands (CGIs), regardless of whether CGIs were at promoter regions or away from genes. Recombination hotspots were also associated with specific transposable elements (TEs), but this association appears indirect due to shared preferences of the transposition machinery and the recombination machinery for accessible open chromatin regions. Our results suggest that CGIs are a major determinant of the localization of recombination hotspots, and we propose that both the distribution of TEs and fine‐scale variation in recombination rate may be associated with the evolution of the epigenetic landscape.     

Crossover homeostasis in yeast meiosis

Crossovers produced by homologous recombination promote accurate chromosome segregation in meiosis and are controlled such that at least one forms per chromosome pair and multiple crossovers are widely spaced. Recombination initiates with an excess number of double-strand breaks made by Spo11 protein. Thus, crossover control involves a decision by which some breaks give crossovers while others follow a predominantly noncrossover pathway(s). To understand this decision, we examined recombination when breaks are reduced in yeast spo11 hypomorphs. We find that crossover levels tend to be maintained at the expense of noncrossovers and that genomic loci differ in expression of this "crossover homeostasis." These findings define a previously unsuspected manifestation of crossover control, i.e., that the crossover/noncrossover ratio can change to maintain crossovers. Our results distinguish between existing models of crossover control and support the hypothesis that an obligate crossover is a genetically programmed event tied to crossover interference.     

The Impact of Recombination Hotspots on Genome Evolution of a Fungal Plant Pathogen

Recombination has an impact on genome evolution by maintaining chromosomal integrity, affecting the efficacy of selection, and increasing genetic variability in populations. Recombination rates are a key determinant of the coevolutionary dynamics between hosts and their pathogens. Historic recombination events created devastating new pathogens, but the impact of ongoing recombination in sexual pathogens is poorly understood. Many fungal pathogens of plants undergo regular sexual cycles, and sex is considered to be a major factor contributing to virulence. We generated a recombination map at kilobase-scale resolution for the haploid plant pathogenic fungus Zymoseptoria tritici. To account for intraspecific variation in recombination rates, we constructed genetic maps from two independent crosses. We localized a total of 10,287 crossover events in 441 progeny and found that recombination rates were highly heterogeneous within and among chromosomes. Recombination rates on large chromosomes were inversely correlated with chromosome length. Short accessory chromosomes often lacked evidence for crossovers between parental chromosomes. Recombination was concentrated in narrow hotspots that were preferentially located close to telomeres. Hotspots were only partially conserved between the two crosses, suggesting that hotspots are short-lived and may vary according to genomic background. Genes located in hotspot regions were enriched in genes encoding secreted proteins. Population resequencing showed that chromosomal regions with high recombination rates were strongly correlated with regions of low linkage disequilibrium. Hence, genes in pathogen recombination hotspots are likely to evolve faster in natural populations and may represent a greater threat to the host.     

The last bastion? X chromosome genotyping of Anopheles gambiae species pair males from a hybrid zone reveals complex recombination within the major candidate ‘genomic island of speciation’

Speciation with gene flow may be aided by reduced recombination helping to build linkage between genes involved in the early stages of reproductive isolation. Reduced recombination on chromosome X has been implicated in speciation within the Anopheles gambiae complex, species of which represent the major Afrotropical malaria vectors. The most recently diverged, morphologically indistinguishable, species pair, A. gambiae and Anopheles coluzzii, ubiquitously displays a ‘genomic island of divergence’ spanning over 4 Mb from chromosome X centromere, which represents a particularly promising candidate region for reproductive isolation genes, in addition to containing the diagnostic markers used to distinguish the species. Very low recombination makes the island intractable for experimental recombination studies, but an extreme hybrid zone in Guinea Bissau offers the opportunity for natural investigation of X‐island recombination. SNP analysis of chromosome X hemizygous males revealed: (i) strong divergence in the X‐island despite a lack of autosomal divergence; (ii) individuals with multiple‐recombinant genotypes, including likely double crossovers and localized gene conversion; (iii) recombination‐driven discontinuity both within and between the molecular species markers, suggesting that the utility of the diagnostics is undermined under high hybridization. The largely, but incompletely protected nature of the X centromeric genomic island is consistent with a primary candidate area for accumulation of adaptive variants driving speciation with gene flow, while permitting some selective shuffling and removal of genetic variation.     

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 Red Queen theory of recombination hotspots

Recombination hotspots are small chromosomal regions, where meiotic crossover events happen with high frequency. Recombination is initiated by a double‐strand break (DSB) that requires the intervention of the molecular repair mechanism. The DSB repair mechanism may result in the exchange of homologous chromosomes (crossover) and the conversion of the allelic sequence that breaks into the one that does not break (biased gene conversion). Biased gene conversion results in a transmission advantage for the allele that does not break, thus preventing recombination and rendering recombination hotspots transient. How is it possible that recombination hotspots persist over evolutionary time (maintaining the average chromosomal crossover rate) when they are self‐destructive? This fundamental question is known as the recombination hotspot paradox and has attracted much attention in recent years. Yet, that attention has not translated into a fully satisfactory answer. No existing model adequately explains all aspects of the recombination hotspot paradox. Here, we formulate an intragenomic conflict model resulting in Red Queen dynamics that fully accounts for all empirical observations regarding the molecular mechanisms of recombination hotspots, the nonrandom targeting of the recombination machinery to hotspots and the evolutionary dynamics of hotspot turnover.     

Engineering meiotic recombination pathways in rice

In the last 15 years, outstanding progress has been made in understanding the function of meiotic genes in the model dicot and monocot plants Arabidopsis and rice (Oryza sativa L.), respectively. This knowledge allowed to modulate meiotic recombination in Arabidopsis and, more recently, in rice. For instance, the overall frequency of crossovers (COs) has been stimulated 2.3‐ and 3.2‐fold through the inactivation of the rice FANCM and RECQ4 DNA helicases, respectively, two genes involved in the repair of DNA double‐strand breaks (DSBs) as noncrossovers (NCOs) of the Class II crossover pathway. Differently, the programmed induction of DSBs and COs at desired sites is currently explored by guiding the SPO11‐1 topoisomerase‐like transesterase, initiating meiotic recombination in all eukaryotes, to specific target regions of the rice genome. Furthermore, the inactivation of 3 meiosis‐specific genes, namely PAIR1, OsREC8 and OsOSD1, in the Mitosis instead of Meiosis (MiMe) mutant turned rice meiosis into mitosis, thereby abolishing recombination and achieving the first component of apomixis, apomeiosis. The successful translation of Arabidopsis results into a crop further allowed the implementation of two breakthrough strategies that triggered parthenogenesis from the MiMe unreduced clonal egg cell and completed the second component of diplosporous apomixis. Here, we review the most recent advances in and future prospects of the manipulation of meiotic recombination in rice and potentially other major crops, all essential for global food security.     

Studies on crossover-specific mutants and the distribution of crossing over in Drosophila females

In Drosophila females, the majority of recombination events do not become crossovers and those that do occur are nonrandomly distributed. Furthermore, a group of Drosophila mutants specifically reduce crossing over, suggesting that crossovers depend on different gene products than noncrossovers. In mei-218 mutants, crossing over is reduced by approximately 90% while noncrossovers and the initiation of recombination remain unchanged. Importantly, the residual crossovers have a more random distribution than wild-type. It has been proposed that mei-218 has a role in establishing the crossover distribution by determining which recombination sites become crossovers. Surprisingly, a diverse group of genes, including those required for double strand break (DSB) formation or repair, have an effect on crossover distribution. Not all of these mutants, however, have a crossover-specific defect like mei-218 and it is not understood why some crossover-defective mutants alter the distribution of crossovers. Intragenic recombination experiments suggest that mei-218 is required for a molecular transition of the recombination intermediate late in the DSB repair pathway. We propose that the changes in crossover distribution in some crossover-defective mutants are a secondary consequence of the crossover reductions. This may be the activation of a regulatory system that ensures at least one crossover per chromosome, and which compensates for an absence of crossovers by attempting to generate them at random locations.     

Srs2 and Sgs1-Top3 suppress crossovers during double-strand break repair in yeast

Very few gene conversions in mitotic cells are associated with crossovers, suggesting that these events are regulated. This may be important for the maintenance of genetic stability. We have analyzed the relationship between homologous recombination and crossing-over in haploid budding yeast and identified factors involved in the regulation of crossover outcomes. Gene conversions unaccompanied by a crossover appear 30 min before conversions accompanied by exchange, indicating that there are two different repair mechanisms in mitotic cells. Crossovers are rare (5%), but deleting the BLM/WRN homolog, SGS1, or the SRS2 helicase increases crossovers 2- to 3-fold. Overexpressing SRS2 nearly eliminates crossovers, whereas overexpression of RAD51 in srs2Delta cells almost completely eliminates the noncrossover recombination pathway. We suggest Sgs1 and its associated topoisomerase Top3 remove double Holliday junction intermediates from a crossover-producing repair pathway, thereby reducing crossovers. Srs2 promotes the noncrossover synthesis-dependent strand-annealing (SDSA) pathway, apparently by regulating Rad51 binding during strand exchange.     

Multiple Pathways Suppress Non-Allelic Homologous Recombination during Meiosis in Saccharomyces cerevisiae

Recombination during meiosis in the form of crossover events promotes the segregation of homologous chromosomes by providing the only physical linkage between these chromosomes. Recombination occurs not only between allelic sites but also between non-allelic (ectopic) sites. Ectopic recombination is often suppressed to prevent non-productive linkages. In this study, we examined the effects of various mutations in genes involved in meiotic recombination on ectopic recombination during meiosis. RAD24, a DNA damage checkpoint clamp-loader gene, suppressed ectopic recombination in wild type in the same pathway as RAD51. In the absence of RAD24, a meiosis-specific recA homolog, DMC1, suppressed the recombination. Homology search and strand exchange in ectopic recombination occurred when either the RAD51 or the DMC1 recA homolog was absent, but was promoted by RAD52. Unexpectedly, the zip1 mutant, which is defective in chromosome synapsis, showed a decrease, rather than an increase, in ectopic recombination. Our results provide evidence for two types of ectopic recombination: one that occurs in wild-type cells and a second that occurs predominantly when the checkpoint pathway is inactivated.     

Rapid evolutionary dynamics in a 2.8‐Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi‐gene families, the DNA sequence of a 2.8‐Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance‐like gene families, was analyzed in the diploid grass Aegilops tauschii, the D‐genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non‐syntenic genes and only 13 ancestral genes are shared among these grass species. These non‐syntenic genes, including prolamin and resistance‐like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non‐syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin‐related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance‐like genes and subsequent differential expansion of each R gene family. The high frequency of non‐syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8‐Mb region with nine‐fold higher than the genome‐wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.     

Replication protein A (AtRPA1a) is required for class I crossover formation but is dispensable for meiotic DNA break repair

Replication protein A (RPA) is involved in many aspects of DNA metabolism including meiotic recombination. Many species possess a single RPA1 gene but Arabidopsis possesses five RPA1 paralogues. This feature has enabled us to gain further insight into the meiotic role of RPA1. Proteomic analysis implicated one of the AtRPA1 family (AtRPA1a) in meiosis. Immunofluorescence studies confirmed that AtRPA1a is associated with meiotic chromosomes from leptotene through to early pachytene. Analysis of an Atrpa1a mutant revealed that AtRPA1a is not essential at early stages in the recombination pathway. DNA double‐strand breaks are repaired in Atrpa1a, but the mutant is defective in the formation of crossovers, exhibiting a 60% reduction in chiasma frequency. Consistent with this, localization of recombination proteins AtRAD51 and AtMSH4 appears normal, whereas the numbers of AtMLH1 and AtMLH3 foci at pachytene are significantly reduced. This suggests that the defect in Atrpa1a is manifested at the stage of second‐end capture. Analysis of Atrpa1a/Atmsh4 and Atrpa1a/Atmlh3 double mutants indicates that loss of AtRPA1a predominantly affects the formation of class I, interference‐dependent crossovers.     

The dynamic loss and gain of introns during the evolution of the Brassicaceae

Sequence comparison allows the detailed analysis of evolution at the nucleotide and amino acid levels, but much less information is known about the structural evolution of genes, i.e. how the number, length and distribution of introns change over time. We constructed a parsimonious model for the evolutionary rate of intron loss (IL) and intron gain (IG) within the Brassicaceae and found that IL/IG has been highly dynamic, with substantial differences between and even within lineages. The divergence of the Brassicaceae lineages I and II marked a dramatic change in the IL rate, with the common ancestor of lineage I losing introns three times more rapidly than the common ancestor of lineage II. Our data also indicate a subsequent declining trend in the rate of IL, although in Arabidopsis thaliana introns continue to be lost at approximately the ancestral rate. Variations in the rate of IL/IG within lineage II have been even more remarkable. Brassica rapa appears to have lost introns approximately 15 times more rapidly than the common ancestor of B. rapa and Schenkiella parvula, and approximately 25 times more rapidly than its sister species Eutrema salsugineum. Microhomology was detected at the splice sites of several dynamic introns suggesting that the non‐homologous end‐joining and double‐strand break repair is a common pathway underlying IL/IG in these species. We also detected molecular signatures typical of mRNA‐mediated IL, but only in B. rapa.     

Clonal reproduction shapes evolution in the lizard malaria parasite Plasmodium floridense

The preponderant clonal evolution hypothesis (PCE) predicts that frequent clonal reproduction (sex between two clones) in many pathogens capable of sexual recombination results in strong linkage disequilibrium and the presence of discrete genetic subdivisions characterized by occasional gene flow. We expand on the PCE and predict that higher rates of clonal reproduction will result in: (1) morphologically cryptic species that exhibit (2) low within‐species variation and (3) recent between‐species divergence. We tested these predictions in the Caribbean lizard malaria parasite Plasmodium floridense using 63 single‐infection samples in lizards collected from across the parasite's range, and sequenced them at two mitochondrial, one apicoplast, and five nuclear genes. We identified 11 provisionally cryptic species within P. floridense, each of which exhibits low intraspecific variation and recent divergence times between species (some diverged approximately 110,000 years ago). Our results are consistent with the hypothesis that clonal reproduction can profoundly affect diversification of species capable of sexual recombination, and suggest that clonal reproduction may have led to a large number of unrecognized pathogen species. The factors that may influence the rates of clonal reproduction among pathogens are unclear, and we discuss how prevalence and virulence may relate to clonal reproduction.     

SHORT INDELS ARE SUBJECT TO INSERTION‐BIASED GENE CONVERSION

Recombination between homologous loci is accompanied by formation of heteroduplexes. Repairing mismatches in heteroduplexes often leads to single nucleotide substitutions in a process known as gene conversion. Gene conversion was shown to be GC‐biased in different organisms; that is, a W(A or T)→S(G or C) substitution is more likely in this process than a S→W substitution. Here, we show that the insertion/deletion ratio for short noncoding indels that reach fixation between species is positively correlated with the recombination rate in Drosophila melanogaster, Homo sapiens, and Saccharomyces cerevisiae. This correlation is both due to an increase of the fixation rate of insertions and decrease of the fixation rate of deletions in regions of high recombination. Whole‐genome data on indel polymorphism and divergence in D. melanogaster rule out mutation biases and selection as the cause of this trend, pointing to insertion‐biased gene conversion as the most likely explanation. The bias toward insertions is the strongest for single‐nucleotide indels, and decreases with indel length. In regions of high recombination rate this bias leads to an up to ～5‐fold excess of fixed short insertions over deletions, and substantially affects the evolution of DNA segments.     

RECOMBINATION AND MIGRATION RATES IN NATURAL POPULATIONS OF BACILLUS SUBTILIS AND BACILLUS MOJAVENSIS

We have investigated the rates of recombination and migration in native populations of two closely related, naturally competent Bacillus species. Native soil isolates of Bacillus subtilis and Bacillus mojavensis were obtained from three continents and, within North America, from populations at a range of geographical distances from one another. The rate of recombination within populations of each species was estimated from restriction-site data for three genes. Recombination was shown to occur within each species at about the same rate as neutral mutation, whatever the geographical scale or phylogenetic scale over which strains were sampled. The rate of migration between populations was estimated by a cladistic analysis and was shown to be high (i.e., Nm > 1), even among populations on different continents. The level of migration within each species is sufficient to prevent neutral geographical divergence within species.     

Phosphorylation of the Synaptonemal Complex Protein Zip1 Regulates the Crossover/Noncrossover Decision during Yeast Meiosis

Interhomolog crossovers promote proper chromosome segregation during meiosis and are formed by the regulated repair of programmed double-strand breaks. This regulation requires components of the synaptonemal complex (SC), a proteinaceous structure formed between homologous chromosomes. In yeast, SC formation requires the “ZMM” genes, which encode a functionally diverse set of proteins, including the transverse filament protein, Zip1. In wild-type meiosis, Zmm proteins promote the biased resolution of recombination intermediates into crossovers that are distributed throughout the genome by interference. In contrast, noncrossovers are formed primarily through synthesis-dependent strand annealing mediated by the Sgs1 helicase. This work identifies a conserved region on the C terminus of Zip1 (called Zip1 4S), whose phosphorylation is required for the ZMM pathway of crossover formation. Zip1 4S phosphorylation is promoted both by double-strand breaks (DSBs) and the meiosis-specific kinase, MEK1/MRE4, demonstrating a role for MEK1 in the regulation of interhomolog crossover formation, as well as interhomolog bias. Failure to phosphorylate Zip1 4S results in meiotic prophase arrest, specifically in the absence of SGS1. This gain of function meiotic arrest phenotype is suppressed by spo11Δ, suggesting that it is due to unrepaired breaks triggering the meiotic recombination checkpoint. Epistasis experiments combining deletions of individual ZMM genes with sgs1-md zip1-4A indicate that Zip1 4S phosphorylation functions prior to the other ZMMs. These results suggest that phosphorylation of Zip1 at DSBs commits those breaks to repair via the ZMM pathway and provides a mechanism by which the crossover/noncrossover decision can be dynamically regulated during yeast meiosis.