Supplementary data are available at Bioinformatics online.

Meiotic recombination initiates from DNA double-strand breaks (DSBs) generated by SPO11 topoisomerase-like complexes. Meiotic DSB frequency varies extensively along eukaryotic chromosomes, with hotspots controlled by chromatin and DNA sequence. To map meiotic DSBs throughout a plant genome, we purified and sequenced SPO11-1-oligonucleotides. SPO11-1-oligos are elevated in gene promoters, terminators, and introns, which is driven by AT-sequence richness that excludes nucleosomes and allows SPO11-1 access. A positive relationship was observed between SPO11-1-oligos and crossovers genome-wide, although fine-scale correlations were weaker. This may reflect the influence of interhomolog polymorphism on crossover formation, downstream from DSB formation. Although H3K4me3 is enriched in proximity to SPO11-1-oligo hotspots at gene 5' ends, H3K4me3 levels do not correlate with DSBs. Repetitive transposons are thought to be recombination silenced during meiosis, to prevent nonallelic interactions and genome instability. Unexpectedly, we found high SPO11-1-oligo levels in nucleosome-depleted Helitron/Pogo/Tc1/Mariner DNA transposons, whereas retrotransposons were coldspots. High SPO11-1-oligo transposons are enriched within gene regulatory regions and in proximity to immunity genes, suggesting a role as recombination enhancers. As transposon mobility in plant genomes is restricted by DNA methylation, we used the DNA methyltransferase mutant to investigate the role of heterochromatin in SPO11-1-oligo distributions. Epigenetic activation of meiotic DSBs in proximity to centromeres and transposons occurred in mutants, coincident with reduced nucleosome occupancy, gain of transcription, and H3K4me3. Together, our work reveals a complex relationship between chromatin and meiotic DSBs within genes and transposons, with significance for the diversity and evolution of plant genomes.

When defining bacterial populations through whole genome sequencing (WGS) the samples often have unknown evolutionary histories. With the increased use of next generation WGS in routine diagnostics, surveillance and epidemiology a vast amount of short read data is available, with phylogenetic trees (dendograms) used to visualise the relationships and similarities between samples. Standard reference and assembly based methods can take substantial amounts of time to generate these phylogenetic relationships, with the computation time often exceeding the time to sequence the samples in the first place. Faster methods (Ondov et al. 2016; Wood and Salzberg 2014) can loosely classify samples into known taxonomic categories, however the loss of granularity means the relationships between samples is reduced. This can be the difference between ruling a sample in or out of an outbreak, which is a clinically important finding for genomic epidemiologists. Other methods (Boratyn et al. 2014) are closed source which prevents independent scrutiny. Saf-fronTree utilises the k-mer profiles between samples to rapidly construct a tree, directly from raw reads in FASTQ format or contigs in FASTA format. It supports NGS data (such as Illumina), 3rd generation long read data (Pacbio/Nanopore) and assembled sequences (FASTA). Firstly, a k-mer count database is constructed for each sample using KMC (Kokot, Długosz, and Deorowicz 2017). Next, the intersection of the k-mer databases is found for each pair of samples, with the number of k-mers in common recorded in a distance matrix. Finally, the distance matrix is used to construct a UPGMA tree (Sokal and Michener 1958) in Newick format. This tree method was chosen as it is fast, however the final result is lower quality than slower methods which perform ancestral sequence reconstructions (Stamatakis 2014). The computational complexity of the algorithm is O(Nˆ2), so is best suited to datasets of less than 50 samples. This can give rapid insights into small datasets in minutes, rather than hours. SaffonTree provides better granular-ity than MLST as it uses more of the underlying genome, can operate at low depth of coverage, is reference free, species agnostic, and has a low memory requirement.

During meiosis, homologous chromosomes undergo crossover recombination, which creates genetic diversity and balances homolog segregation. Despite these critical functions, crossover frequency varies extensively within and between species. Although natural crossover recombination modifier loci have been detected in plants, causal genes have remained elusive. Using natural Arabidopsis thaliana accessions, we identified two major recombination quantitative trait loci (rQTLs) that explain 56.9% of crossover variation in Col×Ler F 2 populations. We mapped rQTL1 to semidominant polymorphisms in HEI10, which encodes a conserved ubiquitin E3 ligase that regulates crossovers. Null hei10 mutants are haploinsufficient, and, using genome-wide mapping and immunocytology, we show that transformation of additional HEI10 copies is sufficient to more than double euchromatic crossovers. However, heterochromatic centromeres remained recombination-suppressed. The strongest HEI10-mediated crossover increases occur in subtelomeric euchromatin, which is reminiscent of sex differences in Arabidopsis recombination. Our work reveals that HEI10 naturally limits Arabidopsis crossovers and has the potential to influence the response to selection.

Histone variants have been proposed to act as determinants for post-translational modifications (PTM) with widespread regulatory functions. In this report, we identify a histone-modifying enzyme that selectively methylates the replication-dependent histone H3 variant H3.1. The crystal structure of the SET domain of the histone H3 lysine 27 (H3K27) methyltransferase ARABIDOPSIS TRITHORAX-RELATED PROTEIN 5 (ATXR5) in complex with a H3.1 peptide shows that ATXR5 contains a bipartite catalytic domain that specifically “reads” alanine 31 of H3.1. Variation at position 31 between H3.1 and replication-independent H3.3 is conserved in plants and animals, and threonine 31 in H3.3 is responsible for inhibiting the activity of ATXR5 and its paralog ATXR6. Our results suggest a simple model for the mitotic inheritance of the heterochromatic mark H3K27me1 and the protection of H3.3-enriched genes against heterochromatization during DNA replication.

Eukaryotic centromeres contain the kinetochore, which connects chromosomes to the spindle allowing segregation. During meiosis, centromeres are suppressed for inter-homolog crossover, as recombination in these regions can cause chromosome missegregation and aneuploidy. Plant centromeres are surrounded by transposon-dense pericentromeric heterochromatin that is epigenetically silenced by histone 3 lysine 9 dimethylation (H3K9me2), and DNA methylation in CG and non-CG sequence contexts. However, the role of these chromatin modifications in control of meiotic recombination in the pericentromeres is not fully understood. Here, we show that disruption of Arabidopsis thaliana H3K9me2 and non-CG DNA methylation pathways, for example, via mutation of the H3K9 methyltransferase genes KYP/SUVH4 SUVH5 SUVH6, or the CHG DNA methyltransferase gene CMT3, increases meiotic recombination in proximity to the centromeres. Using immunocytological detection of MLH1 foci and genotyping by sequencing of recombinant plants, we observe that H3K9me2 and non-CG DNA methylation pathway mutants show increased pericentromeric crossovers. Increased pericentromeric recombination in H3K9me2/non-CG mutants occurs in hybrid and inbred backgrounds and likely involves contributions from both the interfering and noninterfering crossover repair pathways. We also show that meiotic DNA double-strand breaks (DSBs) increase in H3K9me2/non-CG mutants within the pericentromeres, via purification and sequencing of SPO11-1-oligonucleotides. Therefore, H3K9me2 and non-CG DNA methylation exert a repressive effect on both meiotic DSB and crossover formation in plant pericentromeric heterochromatin. Our results may account for selection of enhancer trap Dissociation (Ds) transposons into the CMT3 gene by recombination with proximal transposon launch-pads.

During meiosis, homologous chromosomes undergo reciprocal crossovers, which generate genetic diversity and underpin classical crop improvement. Meiotic recombination initiates from DNA double-strand breaks (DSBs), which are processed into single-stranded DNA that can invade a homologous chromosome. The resulting joint molecules can ultimately be resolved as crossovers. In , competing pathways balance the repair of ∼100-200 meiotic DSBs into ∼10 crossovers per meiosis, with the excess DSBs repaired as noncrossovers. To bias DSB repair toward crossovers, we simultaneously increased dosage of the procrossover E3 ligase gene and introduced mutations in the anticrossovers helicase genes and As and increase interfering and noninterfering crossover pathways, respectively, they combine additively to yield a massive meiotic recombination increase. Interestingly, we also show that increased dosage increases crossover coincidence, which indicates an effect on interference. We also show that patterns of interhomolog polymorphism and heterochromatin drive recombination increases distally towards the subtelomeres in both and backgrounds, while the centromeres remain crossover suppressed. These results provide a genetic framework for engineering meiotic recombination landscapes in plant genomes.

Meiotic crossover frequency varies extensively along chromosomes and is typically concentrated in hotspots. As recombination increases genetic diversity, hotspots are predicted to occur at immunity genes, where variation may be beneficial. A major component of plant immunity is recognition of pathogen Avirulence (Avr) effectors by resistance (R) genes that encode NBS-LRR domain proteins. Therefore, we sought to test whether NBS-LRR genes would overlap with meiotic crossover hotspots using experimental genetics in Arabidopsis thaliana. NBS-LRR genes tend to physically cluster in plant genomes; for example, in Arabidopsis most are located in large clusters on the south arms of chromosomes 1 and 5. We experimentally mapped 1,439 crossovers within these clusters and observed NBS-LRR gene associated hotspots, which were also detected as historical hotspots via analysis of linkage disequilibrium. However, we also observed NBS-LRR gene coldspots, which in some cases correlate with structural heterozygosity. To study recombination at the fine-scale we used high-throughput sequencing to analyze ~1,000 crossovers within the RESISTANCE TO ALBUGO CANDIDA1 (RAC1) R gene hotspot. This revealed elevated intragenic crossovers, overlapping nucleosome-occupied exons that encode the TIR, NBS and LRR domains. The highest RAC1 recombination frequency was promoter-proximal and overlapped CTT-repeat DNA sequence motifs, which have previously been associated with plant crossover hotspots. Additionally, we show a significant influence of natural genetic variation on NBS-LRR cluster recombination rates, using crosses between Arabidopsis ecotypes. In conclusion, we show that a subset of NBS-LRR genes are strong hotspots, whereas others are coldspots. This reveals a complex recombination landscape in Arabidopsis NBS-LRR genes, which we propose results from varying coevolutionary pressures exerted by host-pathogen relationships, and is influenced by structural heterozygosity.