Genetic differentiation and intrinsic genomic features explain variation in recombination hotspots among cocoa tree populations.

BACKGROUND: Recombination plays an important evolutionary role by breaking up haplotypes and shuffling genetic variation. This process impacts the ability of selection to eliminate deleterious mutations or increase the frequency of beneficial mutations in a population. To understand the role of recombination generating and maintaining haplotypic variation in a population, we can construct fine-scale recombination maps. Such maps have been used to study a variety of model organisms and proven to be informative of how selection and demographics shape species-wide variation. Here we present a fine-scale recombination map for ten populations of Theobroma cacao - a non-model, long-lived, woody crop. We use this map to elucidate the dynamics of recombination rates in distinct populations of the same species, one of which is domesticated. RESULTS: Mean recombination rates in range between 2.5 and 8.6 cM/Mb for most populations of T. cacao with the exception of the domesticated Criollo (525 cM/Mb) and Guianna, a more recently established population (46.5 cM/Mb). We found little overlap in the location of hotspots of recombination across populations. We also found that hotspot regions contained fewer known retroelement sequences than expected and were overrepresented near transcription start and termination sites. We find mutations in FIGL-1, a protein shown to downregulate cross-over frequency in Arabidopsis, statistically associated to higher recombination rates in domesticated Criollo. CONCLUSIONS: We generated fine-scale recombination maps for ten populations of Theobroma cacao and used them to understand what processes are associated with population-level variation in this species. Our results provide support to the hypothesis of increased recombination rates in domesticated plants (Criollo population). We propose a testable mechanistic hypothesis for the change in recombination rate in domesticated populations in the form of mutations to a previously identified recombination-suppressing protein. Finally, we establish a number of possible correlates of recombination hotspots that help explain general patterns of recombination in this species.

The recombination landscape of introgression in yeast.

Meiotic recombination is an evolutionary force that acts by breaking up genomic linkage, increasing the efficacy of selection. Recombination is initiated with a double-strand break which is resolved via a crossover, which involves the reciprocal exchange of genetic material between homologous chromosomes, or a non-crossover, which results in small tracts of non-reciprocal exchange of genetic material. Crossover and non-crossover rates vary between species, populations, individuals, and across the genome. In recent years, recombination rate has been associated with the distribution of ancestry derived from past interspecific hybridization (introgression) in a variety of species. We explore this interaction of recombination and introgression by sequencing spores and detecting crossovers and non-crossovers from two crosses of the yeast Saccharomyces uvarum. One cross is between strains which each contain introgression from their sister species, S. eubayanus, while the other cross has no introgression present. We find that the recombination landscape is significantly different between S. uvarum crosses, and that some of these differences can be explained by the presence of introgression in one cross. Crossovers are reduced and non-crossovers are increased in heterozygous introgression compared to syntenic regions in the cross without introgression. This translates to reduced allele shuffling within introgressed regions, and an overall reduction of shuffling on most chromosomes with introgression compared to the syntenic regions and chromosomes without introgression. Our results suggest that hybridization can significantly influence the recombination landscape, and that the reduction in allele shuffling contributes to the initial purging of introgression in the generations following a hybridization event.

Temperature affects recombination rate plasticity and meiotic success between thermotolerant and cold tolerant yeast species.

Meiosis is required for the formation of gametes in all sexually reproducing species and the process is well conserved across the tree of life. However, meiosis is sensitive to a variety of external factors, which can impact chromosome pairing, recombination, and fertility. For example, the optimal temperature for successful meiosis varies between species of plants and animals. This suggests that meiosis is temperature sensitive, and that natural selection may act on variation in meiotic success as organisms adapt to different environmental conditions. To understand how temperature alters the successful completion of meiosis, we utilized two species of the budding yeast Saccharomyces with different temperature preferences: thermotolerant Saccharomyces cerevisiae and cold tolerant Saccharomyces uvarum. We surveyed three metrics of meiosis: sporulation efficiency, spore viability, and recombination rate in multiple strains of each species. As per our predictions, the proportion of cells that complete meiosis and form spores is temperature sensitive, with thermotolerant S. cerevisiae having a higher temperature threshold for successful meiosis than cold tolerant S. uvarum. We confirmed previous observations that S. cerevisiae recombination rate varies between strains and across genomic regions, and add new results that S. uvarum has higher recombination rates than S. cerevisiae. We find that temperature significantly influences recombination rate plasticity in S. cerevisiae and S. uvarum, in agreement with studies in animals and plants. Overall, these results suggest that meiotic thermal sensitivity is associated with organismal thermal tolerance, and may even result in temporal reproductive isolation as populations diverge in thermal profiles.

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.

Local ancestry analysis reveals genomic convergence in extremophile fishes.

The molecular basis of convergent phenotypes is often unknown. However, convergence at a genomic level is predicted when there are large population sizes, gene flow among diverging lineages or strong genetic constraints. We used whole-genome resequencing to investigate genomic convergence in fishes ( Poecilia spp.) that have repeatedly colonized hydrogen sulfide (H2S)-rich environments in Mexico. We identified genomic similarities in both single nucleotide polymorphisms (SNPs) and structural variants (SVs) among independently derived sulfide spring populations, with approximately 1.2% of the genome being shared among sulfidic ecotypes. We compared these convergent genomic regions to candidate genes for H2S adaptation identified from transcriptomic analyses and found that a significant proportion of these candidate genes (8%) were also in regions where sulfidic individuals had similar SNPs, while only 1.7% were in regions where sulfidic individuals had similar SVs. Those candidate genes included genes involved in sulfide detoxification, the electron transport chain (the main toxicity target of H2S) and other processes putatively important for adaptation to sulfidic environments. Regional genomic similarity across independent populations exposed to the same source of selection is consistent with selection on standing variation or introgression of adaptive alleles across divergent lineages. However, combined with previous analyses, our data also support that adaptive changes in mitochondrially encoded subunits arose independently via selection on de novo mutations. Pressing questions remain on what conditions ultimately facilitate the independent rise of adaptive alleles at the same loci in separate populations, and thus, the degree to which evolution is repeatable or predictable. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.