Diversity and determinants of recombination landscapes in flowering plants

Genome sequence of the wild species, Spinacia tetrandra, including a phased sequence of the extensive sex-linked region, revealing partial degeneration in evolutionary strata with unusual properties

Genetic degeneration is a striking feature of Y chromosomes, often involving losses of many genes carried on the X chromosome. However, the time course of gene losses remains unclear. Sex chromosomes of plants evolved more recently than animals' highly degenerated ones, making them ideal for studying degeneration timing. To investigate Spinacia sex chromosome evolution and the time course of degeneration, we compared genome sequences of cultivated Spinacia oleracea, with a small Y-linked region on Chr4, with its two wild relatives. In spinach and its closest relative Spinacia turkestanica, the Y duplication region (YDR) introduced a male-determining factor into Chr4's low-recombining pericentromeric region. In other words, a turnover event occurred in these species' recent common ancestor. The homologous Chr4 of the more distantly related S. tetrandra has a c. 133 Mb completely sex-linked and partially degenerated region, possibly reflecting the ancestral state. Sequence divergence analysis suggests that two 'evolutionary strata' evolved shortly before the two Spinacia lineages split. Consistent with the turnover hypothesis, the YDR of the other two Spinacia species is not within the S. tetrandra older stratum. We discuss the unexpected findings in S. tetrandra that genetic degeneration, genomic rearrangements, and repetitive sequence density are all greatest in the younger stratum.

A well-annotated genome of Apium graveolens var. dulce cv. Challenger, a celery with resistance to Fusarium oxysporum f. sp. apii race 2

Celery (Apium graveolens var. dulce) production can be limited by the fungal pathogen Fusarium oxysporum f. sp. apii (Foa), particularly at temperatures above 22°C. Because celery has a narrow genetic base, an intraspecific admixture of Apium graveolens was developed into cv. Challenger, which is resistant to Foa race 2, the causal agent of Fusarium yellows, but susceptible to Foa race 4, a relatively unrelated causal agent of Fusarium wilt. We assembled a high-quality, chromosome-level physical map of Challenger with 40 464 RNA-based, protein-coding gene models in 3.3 Gbp and anchored it with a genetic map. Although there is high gene density and higher recombination at the ends of the chromosomes, an average of 56% of the genes/chromosome are in lower recombination zones (<0.025 cM/Mb). We identified Challenger's nucleotide-binding and leucine-rich repeat receptors (NLRs) and pattern recognition receptors (PRRs), the two gene families that encode most resistance (R) genes. In three treatment groups (mock-infested or infested with either Foa race 2 or race 4), 243 NLRs and 445 PRRs were quantified in the celery crowns via Quant-Seq 3' mRNA-Seq (Tag-Seq). We compared the genomes of Challenger with that of the previously published cv. Ventura, which is moderately susceptible to Foa race 2. We present a toolbox for genome-assisted breeding for celery that includes annotated gene models, a protocol for genotype-by-sequencing, documentation of the expression of NLRs and PRRs, and a straightforward strategy for introgressing selected NLR superclusters, 83% of which are in higher recombination regions.

Out of Africa: The genomic footprints of Vietnamese Robusta coffee

Vietnam is the main producer of Robusta (Coffea canephora) coffee, but faces several future agronomic challenges. These may be addressed through breeding for improved cultivars and more sustainable cropping systems. For such efforts to be successful and efficient, locally available genetic resources must be understood. Indeed, while C. canephora exhibits high genetic diversity in its native tropical African forests, only a part of it contributed to the worldwide diffusion of Robusta. Here we traced the African origins of Robusta accessions cultivated in the Central Highlands of Vietnam. A total of 126 Robusta accessions from the Vietnam coffee germplasm collection were characterized, including historical, elite and local cultivated clones. Their genetic diversity and origins were inferred through comparisons with wild reference samples using a new set of 261 genome-wide SNPs. A core set of 45 accessions that maximize the genetic distance and allelic richness were identified for conservation and breeding priorities. Full genome sequencing of these individuals helped to closely trace the origins of chromosomal segments back to different, geographically-structured wild African genetic groups. All Vietnamese Robusta accessions displayed Congo Basin (ER group) origins, albeit to various extents. However, we also uncovered contribution from several other genetic groups, variously from the Guinean region (D), the central African Atlantic coast (AG), and Eastern CAR/Uganda (OB), in 31 hybrid individuals. These source groups have been widely used in crossbreeding to develop elite clones. In addition, using whole-genome sequencing data, we also identified various admixture patterns at the chromosome level among the hybrids, which might provide valuable information for selecting breeding materials.

Genetic Isolation among Four Lineages of Silene nutans

Speciation is the process leading to the emergence of new species. While being usually progressive, it can sometimes be fast with rapid emergence of reproductive barriers leading to high level of reproductive isolation. Some reproductive barriers might leave signatures in the genome, through elevated level of genetic differentiation at specific loci. Similar signatures might also be the results of linked selection acting in low recombination regions. Nottingham catchfly (Silene nutans) is a Caryophyllaceae species composed of four genetically differentiated lineages for which strong and asymmetric levels of reproductive isolation have been identified. Using population transcriptomic data from several individuals of the four lineages, we inferred the best evo-demographic scenario leading to the current reproductive isolation of these four lineages. We also tested whether loci exhibiting high level of genetic differentiation represented barrier loci or were located in low recombination regions, evolving under strong influence of linked selection. Overall, the four lineages of S. nutans have diverged in strict isolation, likely during the different glacial period, through migration in distinct glacial refugia. Speciation between these four lineages appeared to be particularly fast, likely due to fast evolving plastid genome accelerating plastid-nuclear co-evolution and the probability of plastid-nuclear incompatibilities in inter-lineage hybrids.

Mixed Outcomes in Recombination Rates After Domestication: Revisiting Theory and Data

The process of domestication has altered many phenotypes. Selection on these phenotypes has long been hypothesised to indirectly select for increases in the genome-wide recombination rate. This hypothesis is potentially consistent with theory on the evolution of the recombination rate, but empirical support has been unclear. We review relevant theory, lab-based experiments, and data comparing recombination rates in wild progenitors and their domesticated counterparts. We utilise population sequencing data and a deep learning method to infer genome-wide recombination rates for new comparisons of chicken/red junglefowl, sheep/mouflon, and goat/bezoar. We find evidence of increased recombination in domesticated goats compared to bezoars but more mixed results in chicken and generally decreased recombination in domesticated sheep compared to mouflon. Our results add to a growing body of literature in plants and animals that finds no consistent evidence of an increase in genome-wide recombination with domestication.

Genomic Divergence Shaped the Genetic Regulation of Meiotic Homologous Recombination in Brassica Allopolyploids

The tight regulation of meiotic recombination between homologs is disrupted in Brassica AAC allotriploids, a genomic configuration that may have facilitated the formation of rapeseed (Brassica napus L.) ∼7,500 years ago. Indeed, the presence of the haploid C genome induces supernumerary crossovers between homologous A chromosomes with dramatically reshaped distribution. However, the genetic mechanisms driving this phenomenon and their divergence between nascent and established lineages remain unclear. To address these concerns, we generated hybrids carrying additional C chromosomes derived either from an established lineage of the allotetraploid B. napus or from its diploid progenitor B. oleracea. We then assessed recombination variation across twelve populations by mapping male meiotic crossovers using single nucleotide polymorphism markers evenly distributed across the sequenced A genome. Our findings reveal that the C09 chromosome of B. oleracea is responsible for the formation of additional crossovers near pericentromeric regions. Interestingly, its counterpart from an established lineage of B. napus shows no significant effect on its own, despite having a similar content of meiotic genes. However, we showed that the B. napus C09 chromosome influences crossover formation through inter-chromosomal epistatic interactions with other specific C chromosomes. These results provide new insights into the genetic regulation of homologous recombination in Brassica and emphasize the role of genomic divergence since the formation of the allopolyploid B. napus.

Independent genetic basis of meiotic crossover positioning and interference in domestic pigs

Meiotic crossover patterning shows huge variation within and between chromosomes, individuals, and species, yet the molecular and evolutionary causes and consequences of this variation remain poorly understood. A key step is to understand the genetic architecture of the crossover rate, positioning, and interference to determine if these factors are governed by common or distinct genetic processes. Here, we investigate individual variation in autosomal crossover count, crossover position (measured as both intra-chromosomal shuffling and distance to telomere), and crossover interference in a large breeding population of domestic pigs (N = 82,474 gametes). We show that all traits are heritable in females at the gamete (h2 = 0.07-0.11) and individual mean levels (h2 = 0.08-0.41). In females, crossover count, and interference are strongly associated with RNF212, but crossover positioning is associated with SYCP2, MEI4, and PRDM9. Our results show that crossover positioning and rate/interference are driven by distinct genetic processes in female pigs and have the capacity to evolve independently.

Hi-reComb: constructing recombination maps from bulk gamete Hi-C sequencing

Recombination is central to genetics and to evolution of sexually reproducing organisms. However, obtaining accurate estimates of recombination rates, and of how they vary along chromosomes, continues to be challenging. To advance our ability to estimate recombination rates, we present Hi-reComb, a new method and software for estimation of recombination maps from bulk gamete chromosome conformation capture sequencing (Hi-C). Simulations show that Hi-reComb produces robust, accurate recombination landscapes. With empirical data from sperm of five fish species we show the advantages of this approach, including joint assessment of recombination maps and large structural variants, map comparisons using bootstrap, and workflows with trio phasing vs. Hi-C phasing. With off-the-shelf library construction and a straightforward rapid workflow, our approach will facilitate routine recombination landscape estimation for a broad range of studies and model organisms in genetics and evolutionary biology. Hi-reComb is open-source and freely available at https://github.com/millanek/Hi-reComb.

Recombination and the role of pseudo-overdominance in polyploid evolution

Natural selection is an imperfect force that can under some conditions fail to prevent the buildup of deleterious mutations. Small population sizes and the lack of recombination are two such scenarios that reduce the efficiency of selection. Under these conditions, the disconnect between deleterious genetic load and individual fitness due to the masking of recessive deleterious mutations in heterozygous individuals may result in the emergence of pseudo-overdominance, wherein the buildup of haplotypes with complementary sets of deleterious mutations results in apparent heterozygote advantage and an increase in linked neutral diversity. In polyploids, the presence of additional allelic copies magnifies this masking effect and may therefore increase the probability of pseudo-overdominance. Here, we simulate the evolution of small diploid and autotetraploid populations to identify the conditions that support the evolution of pseudo-overdominance. We discover that pseudo-overdominance evolves under a much wider range of parameters in autotetraploids than in diploids with identical population sizes, and that in many parts of parameter space there is an inverse relationship between fitness and recombination rate. These results imply that pseudo-overdominance may be more common than previously thought. We conclude by discussing the current evidence for pseudo-overdominance in species with polyploid histories, as well as its implications in agriculture due to the prevalence of polyploidy in crops.

Rapid evolution of recombination landscapes during the divergence of cichlid ecotypes in Lake Masoko

Variation of recombination rate along the genome is of crucial importance to rapid adaptation and organismal diversification. Many unknowns remain regarding how and why recombination landscapes evolve in nature. Here, we reconstruct recombination maps based on linkage disequilibrium and use subsampling and simulations to derive a new measure of recombination landscape evolution: the Population Recombination Divergence Index (PRDI). Using PRDI, we show that fine-scale recombination landscapes differ substantially between two cichlid fish ecotypes of Astatotilapia calliptera that diverged only ~2,500 generations ago. Perhaps surprisingly, recombination landscape differences are not driven by divergence in terms of allele frequency (FST) and nucleotide diversity (Δ(π)): although there is some association, we observe positive PRDI in regions where FST and Δ(π) are zero. We found a stronger association between the evolution of recombination and 47 large haplotype blocks that are polymorphic in Lake Masoko, cover 21% of the genome, and appear to include multiple inversions. Among haplotype blocks, there is a strong and clear association between the degree of recombination divergence and differences between ecotypes in heterozygosity, consistent with recombination suppression in heterozygotes. Overall, our work provides a holistic view of changes in population recombination landscapes during the early stages of speciation with gene flow.

Centromere diversity and its evolutionary impacts on plant karyotypes and plant reproduction

Karyotype changes are a formidable evolutionary force by directly impacting cross-incompatibility, gene dosage, genetic linkage, chromosome segregation, and meiotic recombination landscape. These changes often arise spontaneously and are commonly detected within plant lineages, even between closely related accessions. One element that can influence drastic karyotype changes after only one (or few) plant generations is the alteration of the centromere position, number, distribution, or even its strength. Here, we briefly explore how these different centromere configurations can directly result in karyotype rearrangements, impacting plant reproduction and meiotic recombination.

Comparison of Recombination Rate, Reference Bias, and Unique Pangenomic Haplotypes in Cannabis sativa Using Seven De Novo Genome Assemblies

Genomic characterization of Cannabis sativa has accelerated rapidly in the last decade as sequencing costs have decreased and public and private interest in the species has increased. Here, we present seven new chromosome-level haplotype-phased genomes of C. sativa. All of these genotypes were alive at the time of publication, and several have numerous years of associated phenotype data. We performed a k-mer-based pangenome analysis to contextualize these assemblies within over 200 existing assemblies. This allowed us to identify unique haplotypes and genomic diversity among Cannabis sativa genotypes. We leveraged linkage maps constructed from F2 progeny of two of the assembled genotypes to characterize the recombination rate across the genome showing strong periphery-biased recombination. Lastly, we re-aligned a bulk segregant analysis dataset for the major-effect flowering locus Early1 to several of the new assemblies to evaluate the impact of reference bias on the mapping results and narrow the locus to a smaller region of the chromosome. These new assemblies, combined with the continued propagation of the genotypes, will contribute to the growing body of genomic resources for C. sativa to accelerate future research efforts.

The recombination landscape of the barn owl, from families to populations

Homologous recombination is a meiotic process that generates diversity along the genome and interacts with all evolutionary forces. Despite its importance, studies of recombination landscapes are lacking due to methodological limitations and limited data. Frequently used approaches include linkage mapping based on familial data that provides sex-specific broad-scale estimates of realized recombination and inferences based on population linkage disequilibrium that reveal a more fine-scale resolution of the recombination landscape, albeit dependent on the effective population size and the selective forces acting on the population. In this study, we use a combination of these 2 methods to elucidate the recombination landscape for the Afro-European barn owl (Tyto alba). We find subtle differences in crossover placement between sexes that lead to differential effective shuffling of alleles. Linkage disequilibrium-based estimates of recombination are concordant with family-based estimates and identify large variation in recombination rates within and among linkage groups. Larger chromosomes show variation in recombination rates, while smaller chromosomes have a universally high rate that shapes the diversity landscape. We find that recombination rates are correlated with gene content, genetic diversity, and GC content. We find no conclusive differences in the recombination landscapes between populations. Overall, this comprehensive analysis enhances our understanding of recombination dynamics, genomic architecture, and sex-specific variation in the barn owl, contributing valuable insights to the broader field of avian genomics.

Selection Can Favor a Recombination Landscape That Limits Polygenic Adaptation

Modifiers of recombination rates have been described but the selective pressures acting on them and their effect on adaptation to novel environments remain unclear. We performed experimental evolution in the nematode Caenorhabditis elegans using alternative rec-1 alleles modifying the position of meiotic crossovers along chromosomes without detectable direct fitness effects. We show that adaptation to a novel environment is impaired by the allele that decreases recombination rates in the genomic regions containing fitness variation. However, the allele that impairs adaptation is indirectly favored by selection, because it increases recombination rates and reduces the associations among beneficial and deleterious variation located in its chromosomal vicinity. These results validate theoretical expectations about the evolution of recombination but suggest that genome-wide polygenic adaptation is of little consequence to indirect selection on recombination rate modifiers.

The plant early recombinosome: a high security complex to break DNA during meiosis

KEY MESSAGE

The formacion of numerous unpredictable DNA Double Strand Breaks (DSBs) on chromosomes iniciates meiotic recombination. In this perspective, we propose a 'multi-key lock' model to secure the risky but necesary breaks as well as a 'one per pair of cromatids' model for the topoisomerase-like early recombinosome. During meiosis, homologous chromosomes recombine at few sites of crossing-overs (COs) to ensure correct segregation. The initiation of meiotic recombination involves the formation of DNA double strand breaks (DSBs) during prophase I. Too many DSBs are dangerous for genome integrity: if these DSBs are not properly repaired, it could potentially lead to chromosomal fragmentation. Too few DSBs are also problematic: if the obligate CO cannot form between bivalents, catastrophic unequal segregation of univalents lead to the formation of sterile aneuploid spores. Research on the regulation of the formation of these necessary but risky DSBs has recently advanced in yeast, mammals and plants. DNA DSBs are created by the enzymatic activity of the early recombinosome, a topoisomerase-like complex containing SPO11. This opinion paper reviews recent insights on the regulation of the SPO11 cofactors necessary for the introduction of temporally and spatially controlled DSBs. We propose that a 'multi-key-lock' model for each subunit of the early recombinosome complex is required to secure the formation of DSBs. We also discuss the hypothetical implications that the established topoisomerase-like nature of the SPO11 core-complex can have in creating DSB in only one of the two replicated chromatids of early prophase I meiotic chromosomes. This hypothetical 'one per pair of chromatids' DSB formation model could optimize the faithful repair of the self-inflicted DSBs. Each DSB could use three potential intact homologous DNA sequences as repair template: one from the sister chromatid and the two others from the homologous chromosomes.

Analysis of the winter oilseed rape recombination landscape suggests maternal-paternal bias

Recombination, the reciprocal exchange of DNA between homologous chromosomes, is a mandatory step necessary for meiosis progression. Crossovers between homologous chromosomes generate new combinations of alleles and maintain genetic diversity. Due to genetic, epigenetic, and environmental factors, the recombination landscape is highly heterogeneous along the chromosomes and it also differs between populations and between sexes. Here, we investigated recombination characteristics across the 19 chromosomes of the model allopolyploid crop species oilseed rape (Brassica napus L.), using two unique multiparental populations derived from two genetically divergent founder pools, each of which comprised 50 genetically diverse founder accessions. A fully balanced, pairwise chain-crossing scheme was utilized to create each of the two populations. A total of 3213 individuals, spanning five successive generations, were genotyped using a 15K SNP array. We observed uneven distribution of recombination along chromosomes, with some genomic regions undergoing substantially more frequent recombination in both populations. In both populations, maternal recombination events were more frequent than paternal recombination. This study provides unique insight into the recombination landscape at chromosomal level and reveals a maternal-paternal bias for recombination number with implications for breeding.

Mating systems and recombination landscape strongly shape genetic diversity and selection in wheat relatives

How and why genetic diversity varies among species is a long-standing question in evolutionary biology. Life history traits have been shown to explain a large part of observed diversity. Among them, mating systems have one of the strongest impacts on genetic diversity, with selfing species usually exhibiting much lower diversity than outcrossing relatives. Theory predicts that a high rate of selfing amplifies selection at linked sites, reducing genetic diversity genome-wide, but frequent bottlenecks and rapid population turn-over could also explain low genetic diversity in selfers. However, how linked selection varies with mating systems and whether it is sufficient to explain the observed difference between selfers and outcrossers has never been tested. Here, we used the Aegilops/Triticum grass species, a group characterized by contrasted mating systems (from obligate outcrossing to high selfing) and marked recombination rate variation across the genome, to quantify the effects of mating system and linked selection on patterns of neutral and selected polymorphism. By analyzing phenotypic and transcriptomic data of 13 species, we show that selfing strongly affects genetic diversity and the efficacy of selection by amplifying the intensity of linked selection genome-wide. In particular, signatures of adaptation were only found in the highly recombining regions in outcrossing species. These results bear implications for the evolution of mating systems and, more generally, for our understanding of the fundamental drivers of genetic diversity.

Genetic background affects the strength of crossover interference in house mice

Meiotic recombination is required for faithful chromosome segregation in most sexually reproducing organisms and shapes the distribution of genetic variation in populations. Both the overall rate and the spatial distribution of crossovers vary within and between species. Adjacent crossovers on the same chromosome tend to be spaced more evenly than expected at random, a phenomenon known as crossover interference. Although interference has been observed in many taxa, the factors that influence the strength of interference are not well understood. We used house mice (Mus musculus), a well-established model system for understanding recombination, to study the effects of genetics and age on recombination rate and interference in the male germline. We analyzed crossover positions in 503 progeny from reciprocal F1 hybrids between inbred strains representing the three major subspecies of house mice. Consistent with previous studies, autosomal alleles from M. m. musculus tend to increase recombination rate, while inheriting a M. m. musculus X chromosome decreases recombination rate. Old males transmit an average of 0.6 more crossovers per meiosis (5.0%) than young males, though the effect varies across genetic backgrounds. We show that the strength of crossover interference depends on genotype, providing a rare demonstration that interference evolves over short timescales. Differences between reciprocal F1s suggest that X-linked factors modulate the strength of interference. Our findings motivate additional comparisons of interference among recently diverged species and further examination of the role of paternal age in determining the number and positioning of crossovers.

Mutation potentiates migration swamping in polygenic local adaptation

Locally adapted traits can exhibit a wide range of genetic architectures, from pronounced divergence at a few loci to small frequency divergence at many loci. The type of architecture that evolves depends strongly on the migration rate, as weakly selected loci experience swamping and do not make lasting contributions to divergence. Simulations from previous studies showed that even when mutations are strongly selected and should resist migration swamping, the architecture of adaptation can collapse and become transient at high mutation rates. Here, we use an analytical two-population model to study how this transition in genetic architecture depends upon population size, strength of selection, and parameters describing the mutation process. To do this, we develop a mathematical theory based on the diffusion approximation to predict the threshold mutation rate above which the transition occurs. We find that this performs well across a wide range of parameter space, based on comparisons with individual-based simulations. The threshold mutation rate depends most strongly on the average effect size of mutations, weakly on the strength of selection, and marginally on the population size. Across a wide range of the parameter space, we observe that the transition to a transient architecture occurs when the trait-wide mutation rate is 10-3-10-2, suggesting that this phenomenon is potentially relevant to complex traits with a large mutational target. On the other hand, based on the apparent stability of genetic architecture in many classic examples of local adaptation, our theory suggests that per-trait mutation rates are often relatively low.

Alternating between even and odd ploidy levels switches on and off the recombination control, even near the centromeres

Meiotic recombination is a key biological process in plant evolution and breeding, as it generates genetic diversity in each generation through the formation of crossovers (COs). However, due to their importance in genome stability, COs are highly regulated in frequency and distribution. We previously demonstrated that this strict regulation of COs can be modified, both in terms of CO frequency and distribution, in allotriploid Brassica hybrids (2n = 3x = 29; AAC) resulting from a cross between Brassica napus (2n = 4x = 38; AACC) and Brassica rapa (2n = 2x = 20; AA). Using the recently updated B. napus genome now including pericentromeres, we demonstrated that COs occur in these cold regions in allotriploids, as close as 375 kb from the centromere. Reverse transcription quantitative PCR (RT-qPCR) of various meiotic genes indicated that Class I COs are likely involved in the increased recombination frequency observed in allotriploids. We also demonstrated that this modified recombination landscape can be maintained via successive generations of allotriploidy (odd ploidy level). This deregulated meiotic behavior reverts to strict regulation in allotetraploid (even ploidy level) progeny in the second generation. Overall, we provide an easy way to manipulate tight recombination control in a polyploid crop.

Relationship Between Genetic and Phenotypic Variations in Natural Populations of Perennial and Biennial Sagebrush

Plant responses to environmental heterogeneity depend on life-history traits, which could relate to phenotypical and genetic characteristics. To elucidate this relationship, we examined the variation in population genetics and functional traits of short- and long-lived Artemisia species that are co-occurring in the steppes of Mongolia. Mongolian steppes represent stressful and water-limited habitats, demanding phenotypic modifications in the short term and/or genetic adaptation in the long term. However, detailed knowledge is missing about both plant phenotypic and genetic differentiation, and their interrelationships in temperate grasslands. Here, we investigated 21 populations of the widely distributed subshrub Artemisia frigida and the herbaceous biennial Artemisia scoparia. Genetic variation was assessed with newly developed simple sequence repeats (SSRs) markers. Functional trait data were collected from each individual, and data on environmental variables was collected for each population. We detected significantly higher genetic diversity in the biennial species (H E = 0.86) compared with the perennial (H E = 0.79). For both species, the largest share of genetic variation was partitioned within populations (96%). Population genetic structure in the biennial A. scoparia was weak, while the perennial A. frigida showed some spatial genetic structure, which was impacted by geographical factors, soil nutrients, and precipitation amount. Morphology-related functional traits (i.e., plant height) were predominantly associated with environmental variables rather than with genetic variation, whereas physiology-related trait (i.e., specific leaf area [SLA]) was partly genetically determined.

Diversity in Recombination Hotspot Characteristics and Gene Structure Shape Fine-Scale Recombination Patterns in Plant Genomes

During the meiosis of many eukaryote species, crossovers tend to occur within narrow regions called recombination hotspots. In plants, it is generally thought that gene regulatory sequences, especially promoters and 5' to 3' untranslated regions, are enriched in hotspots, but this has been characterized in a handful of species only. We also lack a clear description of fine-scale variation in recombination rates within genic regions and little is known about hotspot position and intensity in plants. To address this question, we constructed fine-scale recombination maps from genetic polymorphism data and inferred recombination hotspots in 11 plant species. We detected gradients of recombination in genic regions in most species, yet gradients varied in intensity and shape depending on specific hotspot locations and gene structure. To further characterize recombination gradients, we decomposed them according to gene structure by rank and number of exons. We generalized the previously observed pattern that recombination hotspots are organized around the boundaries of coding sequences, especially 5' promoters. However, our results also provided new insight into the relative importance of the 3' end of genes in some species and the possible location of hotspots away from genic regions in some species. Variation among species seemed driven more by hotspot location among and within genes than by differences in size or intensity among species. Our results shed light on the variation in recombination rates at a very fine scale, revealing the diversity and complexity of genic recombination gradients emerging from the interaction between hotspot location and gene structure.

Genetic background affects the strength of crossover interference in house mice

Meiotic recombination is required for faithful chromosome segregation in most sexually reproducing organisms and shapes the distribution of genetic variation in populations. Both the overall rate and the spatial distribution of crossovers vary within and between species. Adjacent crossovers on the same chromosome tend to be spaced more evenly than expected at random, a phenomenon known as crossover interference. Although interference has been observed in many taxa, the factors that influence the strength of interference are not well understood. We used house mice (Mus musculus), a well-established model system for understanding recombination, to study the effects of genetics and age on recombination rate and interference in the male germline. We analyzed crossover positions in 503 progeny from reciprocal F1 hybrids between inbred strains representing the three major subspecies of house mice. Consistent with previous studies, autosomal alleles from M. m. musculus tend to increase recombination rate, while inheriting a M. m. musculus X chromosome decreases recombination rate. Old males transmit an average of 0.6 more crossovers per meiosis (5.0%) than young males, though the effect varies across genetic backgrounds. We show that the strength of crossover interference depends on genotype, providing a rare demonstration that interference evolves over short timescales. Differences between reciprocal F1s suggest that X-linked factors modulate the strength of interference. Our findings motivate additional comparisons of interference among recently diverged species and further examination of the role of paternal age in determining the number and positioning of crossovers.

Understanding the Genetic Basis of Variation in Meiotic Recombination: Past, Present, and Future

Meiotic recombination is a fundamental feature of sexually reproducing species. It is often required for proper chromosome segregation and plays important role in adaptation and the maintenance of genetic diversity. The molecular mechanisms of recombination are remarkably conserved across eukaryotes, yet meiotic genes and proteins show substantial variation in their sequence and function, even between closely related species. Furthermore, the rate and distribution of recombination shows a huge diversity within and between chromosomes, individuals, sexes, populations, and species. This variation has implications for many molecular and evolutionary processes, yet how and why this diversity has evolved is not well understood. A key step in understanding trait evolution is to determine its genetic basis-that is, the number, effect sizes, and distribution of loci underpinning variation. In this perspective, I discuss past and current knowledge on the genetic basis of variation in recombination rate and distribution, explore its evolutionary implications, and present open questions for future research.

Insights into spinach domestication from genome sequences of two wild spinach progenitors, Spinacia turkestanica and Spinacia tetrandra

Cultivated spinach (Spinacia oleracea) is a dioecious species. We report high-quality genome sequences for its two closest wild relatives, Spinacia turkestanica and Spinacia tetrandra, which are also dioecious, and are used to study the genetics of spinach domestication. Using a combination of genomic approaches, we assembled genomes of both these species and analyzed them in comparison with the previously assembled S. oleracea genome. These species diverged c. 6.3 million years ago (Ma), while cultivated spinach split from S. turkestanica 0.8 Ma. In all three species, all six chromosomes include very large gene-poor, repeat-rich regions, which, in S. oleracea, are pericentromeric regions with very low recombination rates in both male and female genetic maps. We describe population genomic evidence that the similar regions in the wild species also recombine rarely. We characterized 282 structural variants (SVs) that have been selected during domestication. These regions include genes associated with leaf margin type and flowering time. We also describe evidence that the downy mildew resistance loci of cultivated spinach are derived from introgression from both wild spinach species. Collectively, this study reveals the genome architecture of spinach assemblies and highlights the importance of SVs during the domestication of cultivated spinach.

Identification of a major QTL conferring resistance to wheat yellow mosaic virus derived from the winter wheat 'Hokkai 240' on chromosome 2AS

Wheat yellow mosaic disease is a soilborne disease caused by wheat yellow mosaic virus (WYMV). Symptoms include yellow mosaic coloring of leaves, stunting, and growth inhibition. Here we conducted a detailed analysis of resistance to this virus in winter wheat 'Hokkai 240' by carrying out inoculation tests of WYMV and conducting field tests. The resistance level observed in 'Hokkai 240' was compared with those in varieties harboring known resistance genes. In the inoculation tests, 'Hokkai 240' showed resistance to WYMV Pathotypes I and II and partial resistance to Pathotype III. This result was contrary to the sensitive responses to the three pathotypes exhibited by the variety harboring resistance gene on chromosome 2DL. In fields infected with WYMV Pathotypes II and III, 'Hokkai 240' plants exhibited few disease symptoms and little proliferation of the virus. By analyzing the quantitative trait loci (QTLs) in recombinant inbred lines from a cross between 'Hokkai 240' and 'Nanbukomugi', a single major QTL, Q.Ymhk, from 'Hokkai 240', which had significant effects on Pathotypes II and III of the virus, was detected in the proximity of snp4212 and snp4215 mapped on chromosome 2AS. These results indicate that Q.Ymhk may be useful for developing broad resistance to WYMV in wheat breeding programs.

Chromosome segregation during spermatogenesis occurs through a unique center-kinetic mechanism in holocentric moth species

Precise regulation of chromosome dynamics in the germline is essential for reproductive success across species. Yet, the mechanisms underlying meiotic chromosomal events such as homolog pairing and chromosome segregation are not fully understood in many species. Here, we employ Oligopaint DNA FISH to investigate mechanisms of meiotic homolog pairing and chromosome segregation in the holocentric pantry moth, Plodia interpunctella, and compare our findings to new and previous studies in the silkworm moth, Bombyx mori, which diverged from P. interpunctella over 100 million years ago. We find that pairing in both Bombyx and Plodia spermatogenesis is initiated at gene-rich chromosome ends. Additionally, both species form rod shaped cruciform-like bivalents at metaphase I. However, unlike the telomere-oriented chromosome segregation mechanism observed in Bombyx, Plodia can orient bivalents in multiple different ways at metaphase I. Surprisingly, in both species we find that kinetochores consistently assemble at non-telomeric loci toward the center of chromosomes regardless of where chromosome centers are located in the bivalent. Additionally, sister kinetochores do not seem to be paired in these species. Instead, four distinct kinetochores are easily observed at metaphase I. Despite this, we find clear end-on microtubule attachments and not lateral microtubule attachments co-orienting these separated kinetochores. These findings challenge the classical view of segregation where paired, poleward-facing kinetochores are required for accurate homolog separation in meiosis I. Our studies here highlight the importance of exploring fundamental processes in non-model systems, as employing novel organisms can lead to the discovery of novel biology.

Why should we study plant sex chromosomes?

Understanding plant sex chromosomes involves studying interactions between developmental and physiological genetics, genome evolution, and evolutionary ecology. We focus on areas of overlap between these. Ideas about how species with separate sexes (dioecious species, in plant terminology) can evolve are even more relevant to plants than to most animal taxa because dioecy has evolved many times from ancestral functionally hermaphroditic populations, often recently. One aim of studying plant sex chromosomes is to discover how separate males and females evolved from ancestors with no such genetic sex-determining polymorphism, and the diversity in the genetic control of maleness vs femaleness. Different systems share some interesting features, and their differences help to understand why completely sex-linked regions may evolve. In some dioecious plants, the sex-determining genome regions are physically small. In others, regions without crossing over have evolved sometimes extensive regions with properties very similar to those of the familiar animal sex chromosomes. The differences also affect the evolutionary changes possible when the environment (or pollination environment, for angiosperms) changes, as dioecy is an ecologically risky strategy for sessile organisms. Dioecious plants have repeatedly reverted to cosexuality, and hermaphroditic strains of fruit crops such as papaya and grapes are desired by plant breeders. Sex-linked regions are predicted to become enriched in genes with sex differences in expression, especially when higher expression benefits one sex function but harms the other. Such trade-offs may be important for understanding other plant developmental and physiological processes and have direct applications in plant breeding.

Divergence time shapes gene reuse during repeated adaptation

When diverse lineages repeatedly adapt to similar environmental challenges, the extent to which the same genes are involved (gene reuse) varies across systems. We propose that divergence time among lineages is a key factor driving this variability: as lineages diverge, the extent of gene reuse should decrease due to reductions in allele sharing, functional differentiation among genes, and restructuring of genome architecture. Indeed, we show that many genomic studies of repeated adaptation find that more recently diverged lineages exhibit higher gene reuse during repeated adaptation, but the relationship becomes less clear at older divergence time scales. Thus, future research should explore the factors shaping gene reuse and their interplay across broad divergence time scales for a deeper understanding of evolutionary repeatability.

rec-1 loss of function increases recombination in the central gene clusters at the expense of autosomal pairing centers

Meiotic control of crossover (CO) number and position is critical for homologous chromosome segregation and organismal fertility, recombination of parental genotypes, and the generation of novel genetic combinations. We here characterize the recombination rate landscape of a rec-1 loss of function modifier of CO position in Caenorhabditis elegans, one of the first ever modifiers discovered. By averaging CO position across hermaphrodite and male meioses and by genotyping 203 single-nucleotide variants covering about 95% of the genome, we find that the characteristic chromosomal arm-center recombination rate domain structure is lost in the loss of function rec-1 mutant. The rec-1 loss of function mutant smooths the recombination rate landscape but is insufficient to eliminate the nonuniform position of CO. Lower recombination rates in the rec-1 mutant are particularly found in the autosomal arm domains containing the pairing centers. We further find that the rec-1 mutant is of little consequence for organismal fertility and egg viability and thus for rates of autosomal nondisjunction. It nonetheless increases X chromosome nondisjunction rates and thus male appearance. Our findings question the maintenance of recombination rate heritability and genetic diversity among C. elegans natural populations, and they further suggest that manipulating genetic modifiers of CO position will help find quantitative trait loci located in low-recombining genomic regions normally refractory to discovery.

The genome-wide meiotic recombination landscape in ciliates and its implications for crossover regulation and genome evolution

Meiotic recombination is essential for sexual reproduction and its regulation has been extensively studied in many taxa. However, genome-wide recombination landscape has not been reported in ciliates and it remains unknown how it is affected by the unique features of ciliates: the synaptonemal complex (SC)-independent meiosis and the nuclear dimorphism. Here, we show the recombination landscape in the model ciliate Tetrahymena thermophila by analyzing single-nucleotide polymorphism datasets from 38 hybrid progeny. We detect 1021 crossover (CO) events (35.8 per meiosis), corresponding to an overall CO rate of 9.9 cM/Mb. However, gene conversion by non-crossover is rare (1.03 per meiosis) and not biased towards G or C alleles. Consistent with the reported roles of SC in CO interference, we find no obvious sign of CO interference. CO tends to occur within germ-soma common genomic regions and many of the 44 identified CO hotspots localize at the centromeric or subtelomeric regions. Gene ontology analyses show that CO hotspots are strongly associated with genes responding to environmental changes. We discuss these results with respect to how nuclear dimorphism has potentially driven the formation of the observed recombination landscape to facilitate environmental adaptation and the sharing of machinery among meiotic and somatic recombination.

Meiotic recombination dynamics in plants with repeat-based holocentromeres shed light on the primary drivers of crossover patterning

Centromeres strongly affect (epi)genomic architecture and meiotic recombination dynamics, influencing the overall distribution and frequency of crossovers. Here we show how recombination is regulated and distributed in the holocentric plant Rhynchospora breviuscula, a species with diffused centromeres. Combining immunocytochemistry, chromatin analysis and high-throughput single-pollen sequencing, we discovered that crossover frequency is distally biased, in sharp contrast to the diffused distribution of hundreds of centromeric units and (epi)genomic features. Remarkably, we found that crossovers were abolished inside centromeric units but not in their proximity, indicating the absence of a canonical centromere effect. We further propose that telomere-led synapsis of homologues is the feature that best explains the observed recombination landscape. Our results hint at the primary influence of mechanistic features of meiotic pairing and synapsis rather than (epi)genomic features and centromere organization in determining the distally biased crossover distribution in R. breviuscula, whereas centromeres and (epi)genetic properties only affect crossover positioning locally.

Evolution of a plant sex chromosome driven by expanding pericentromeric recombination suppression

Recombination suppression around sex-determining gene(s) is a key step in evolution of sex chromosomes, but it is not well understood how it evolves. Recently evolved sex-linked regions offer an opportunity to understand the mechanisms of recombination cessation. This paper analyses such a region on Silene latifolia (Caryophyllaceae) sex chromosomes, where recombination was suppressed in the last 120 thousand years ("stratum 3"). Locating the boundaries of the stratum 3 in S. latifolia genome sequence revealed that this region is far larger than assumed previously-it is about 14 Mb long and includes 202 annotated genes. A gradient of X:Y divergence detected in the stratum 3, with divergence increasing proximally, indicates gradual recombination cessation, possibly caused by expansion of pericentromeric recombination suppression (PRS) into the pseudoautosomal region. Expansion of PRS was also the likely cause for the formation of the older stratum 2 on S. latifolia sex chromosomes. The role of PRS in sex chromosome evolution has been underappreciated, but it may be a significant factor, especially in the species with large chromosomes where PRS is often extensive.

Pericentromeric recombination suppression and the 'large X effect' in plants

X chromosome was reported to be a major contributor to isolation between closely related species-the 'large X' effect (LXE). The causes of LXE are not clear, but the leading theory is that it is caused by recessive species incompatibilities exposed in the phenotype due to the hemizygosity of X-linked genes in the heterogametic sex. However, the LXE was also reported in species with relatively recently evolved sex chromosomes where Y chromosome is not completely degenerate and X-linked genes are not hemizygous, such as the plant Silene latifolia. Recent genome sequencing and detailed genetic mapping in this species revealed a massive (> 330 Mb) non- or rarely-recombining pericentromeric region on the X chromosome (Xpr) that comprises ~ 90% of the chromosome and over 13% of the entire genome. If any of the Xpr genes are involved in species incompatibilities, this would oppose interspecific gene flow for other genes tightly linked in the Xpr. Here we test the hypothesis that the previously reported LXE in S. latifolia is caused by the lack of recombination on most of the X chromosome. Based on genome-wide analysis of DNA polymorphism and gene expression in S. latifolia and its close cross-compatible relative S. dioica, we report that the rarely-recombining regions represent a significant barrier for interspecific gene flow. We found little evidence for any additional factors contributing to the LXE, suggesting that extensive pericentromeric recombination suppression on the X-chromosome is the major if not the only cause of the LXE in S. latifolia and S. dioica.

Epigenetic regulation during meiosis and crossover

Meiosis is a distinctive type of cell division that reorganizes genetic material between generations. The initial stages of meiosis consist of several crucial steps which include double strand break, homologous chromosome pairing, break repair and crossover. Crossover frequency varies depending on the position on the chromosome, higher at euchromatin region and rare at heterochromatin, centromeres, telomeres and ribosomal DNA. Crossover positioning is dependent on various factors, especially epigenetic modifications. DNA methylation, histone post-translational modifications, histone variants and non-coding RNAs are most probably playing an important role in positioning of crossovers on a chromosomal level as well as hotspot level. DNA methylation negatively regulates crossover frequency and its effect is visible in centromeres, pericentromeres and heterochromatin regions. Pericentromeric chromatin and heterochromatin mark studies have been a centre of attraction in meiosis. Crossover hotspots are associated with euchromatin regions having specific chromatin modifications such as H3K4me3, H2A.Z. and H3 acetylation. This review will provide the current understanding of the epigenetic role in plants during meiotic recombination, chromosome synapsis, double strand break and hotspots with special attention to euchromatin and heterochromatin marks. Further, the role of epigenetic modifications in regulating meiosis and crossover in other organisms is also discussed.

Meiotic recombination is confirmed to be unusually high in the fission yeast Schizosaccharomyces pombe

In most eukaryotes, meiotic crossovers (COs) are limited to 1-3 per chromosome, and are prevented from occurring close to one another by CO interference. The fission yeast Schizosaccharomyces pombe, an exception to these general rules, was reported to have the highest CO number per chromosome and no or weak interference. However, global CO frequency was indirectly estimated, calling for confirmation. Here, we used an innovative strategy to determine COs genome-wide in S. pombe. We confirmed weak CO interference, acting at physical distances compatible with the patterning of recombination precursors. We revealed a slight co-variation in CO number between chromosomes, suggesting that a limiting pro-CO factor varies between meiocytes. CO number per chromosome varies proportionally with chromosome size, with the three chromosomes having, on average, 15.9, 12.5, and 7.0 COs, respectively. This reinforces S. pombe's status as the eukaryote with the highest CO number per chromosome described to date.

Telomere-to-telomere haplotype-resolved reference genome reveals subgenome divergence and disease resistance in triploid Cavendish banana

Banana is one of the most important crops of the world. Cavendish-type bananas, which have a monospecific Musa acuminata origin (AAA), account for around half of the global banana production, thereby are of great significance for human societies. However, until now, the high-quality haplotype-resolved reference genome was still undecoded for banana cultivars. Here, we reported the telomere-to-telomere (T2T) and haplotype-resolved reference genome of 'Baxijiao' (Cavendish) consisting of three haploid assemblies. The sizes of the three haploid assemblies were estimated to be 477.16 Mb, 477.18 Mb, and 469.57 Mb, respectively. Although with monospecific origins, the three haploid assemblies showed great differences with low levels of sequence collinearity. Several large reciprocal translocations were identified among chromosomes 1, 4, and 7. An expansion of gene families that might affect fruit quality and aroma was detected, such as those belonging to sucrose/disaccharide/oligosaccharide catabolic processes, sucrose metabolic process, starch metabolic process, and aromatic compound biosynthetic process. Besides, an expansion of gene families related to anther and pollen development was observed, which could be associated with parthenocarpy and sterility of the Cavendish cultivar. Finally, much fewer resistance genes were identified in 'Baxijiao' than in M. acuminata, particularly in the gene clusters in chromosomes 3 and 10, providing potential targets to explore for molecular analysis of disease resistance in banana. This T2T haplotype-resolved reference genome will thus be a valuable genetic resource for biological studies, molecular breeding, and genetic improvement of banana.

Nascent evolution of recombination rate differences as a consequence of chromosomal rearrangements

Reshuffling of genetic variation occurs both by independent assortment of chromosomes and by homologous recombination. Such reshuffling can generate novel allele combinations and break linkage between advantageous and deleterious variants which increases both the potential and the efficacy of natural selection. Here we used high-density linkage maps to characterize global and regional recombination rate variation in two populations of the wood white butterfly (Leptidea sinapis) that differ considerably in their karyotype as a consequence of at least 27 chromosome fissions and fusions. The recombination data were compared to estimates of genetic diversity and measures of selection to assess the relationship between chromosomal rearrangements, crossing over, maintenance of genetic diversity and adaptation. Our data show that the recombination rate is influenced by both chromosome size and number, but that the difference in the number of crossovers between karyotypes is reduced as a consequence of a higher frequency of double crossovers in larger chromosomes. As expected from effects of selection on linked sites, we observed an overall positive association between recombination rate and genetic diversity in both populations. Our results also revealed a significant effect of chromosomal rearrangements on the rate of intergenic diversity change between populations, but limited effects on polymorphisms in coding sequence. We conclude that chromosomal rearrangements can have considerable effects on the recombination landscape and consequently influence both maintenance of genetic diversity and efficiency of selection in natural populations.

The origin and evolution of sex chromosomes, revealed by sequencing of the Silene latifolia female genome

White campion (Silene latifolia, Caryophyllaceae) was the first vascular plant where sex chromosomes were discovered. This species is a classic model for studies on plant sex chromosomes due to presence of large, clearly distinguishable X and Y chromosomes that originated de novo about 11 million years ago (mya), but lack of genomic resources for this relatively large genome (∼2.8 Gb) remains a significant hurdle. Here we report S. latifolia female genome assembly integrated with sex-specific genetic maps of this species, focusing on sex chromosomes and their evolution. The analysis reveals a highly heterogeneous recombination landscape with strong reduction in recombination rate in the central parts of all chromosomes. Recombination on the X chromosome in female meiosis primarily occurs at the very ends, and over 85% of the X chromosome length is located in a massive (∼330 Mb) gene-poor, rarely recombining pericentromeric region (Xpr). The results indicate that the non-recombining region on the Y chromosome (NRY) initially evolved in a relatively small (∼15 Mb), actively recombining region at the end of the q-arm, possibly as a result of inversion on the nascent X chromosome. The NRY expanded about 6 mya via linkage between the Xpr and the sex-determining region, which may have been caused by expanding pericentromeric recombination suppression on the X chromosome. These findings shed light on the origin of sex chromosomes in S. latifolia and yield genomic resources to assist ongoing and future investigations into sex chromosome evolution.

Patterns, mechanisms, and consequences of homoeologous exchange in allopolyploid angiosperms: a genomic and epigenomic perspective

Allopolyploids result from hybridization between different evolutionary lineages coupled with genome doubling. Homoeologous chromosomes (chromosomes with common shared ancestry) may undergo recombination immediately after allopolyploid formation and continue over successive generations. The outcome of this meiotic pairing behavior is dynamic and complex. Homoeologous exchanges (HEs) may lead to the formation of unbalanced gametes, reduced fertility, and selective disadvantage. By contrast, HEs could act as sources of novel evolutionary substrates, shifting the relative dosage of parental gene copies, generating novel phenotypic diversity, and helping the establishment of neo-allopolyploids. However, HE patterns vary among lineages, across generations, and even within individual genomes and chromosomes. The causes and consequences of this variation are not fully understood, though interest in this evolutionary phenomenon has increased in the last decade. Recent technological advances show promise in uncovering the mechanistic basis of HEs. Here, we describe recent observations of the common patterns among allopolyploid angiosperm lineages, underlying genomic and epigenomic features, and consequences of HEs. We identify critical research gaps and discuss future directions with far-reaching implications in understanding allopolyploid evolution and applying them to the development of important phenotypic traits of polyploid crops.

The genomics of linkage drag in inbred lines of sunflower

Crop wild relatives represent valuable sources of alleles for crop improvement, including adaptation to climate change and emerging diseases. However, introgressions from wild relatives might have deleterious effects on desirable traits, including yield, due to linkage drag. Here, we analyzed the genomic and phenotypic impacts of wild introgressions in inbred lines of cultivated sunflower to estimate the impacts of linkage drag. First, we generated reference sequences for seven cultivated and one wild sunflower genotype, as well as improved assemblies for two additional cultivars. Next, relying on previously generated sequences from wild donor species, we identified introgressions in the cultivated reference sequences, as well as the sequence and structural variants they contain. We then used a ridge-regression best linear unbiased prediction (BLUP) model to test the effects of the introgressions on phenotypic traits in the cultivated sunflower association mapping population. We found that introgression has introduced substantial sequence and structural variation into the cultivated sunflower gene pool, including >3,000 new genes. While introgressions reduced genetic load at protein-coding sequences, they mostly had negative impacts on yield and quality traits. Introgressions found at high frequency in the cultivated gene pool had larger effects than low-frequency introgressions, suggesting that the former likely were targeted by artificial selection. Also, introgressions from more distantly related species were more likely to be maladaptive than those from the wild progenitor of cultivated sunflower. Thus, breeding efforts should focus, as far as possible, on closely related and fully compatible wild relatives.