MYB-related transcription factors control chloroplast biogenesis.

Chloroplast biogenesis is dependent on master regulators from the GOLDEN2-LIKE (GLK) family of transcription factors. However, glk mutants contain residual chlorophyll, indicating that other proteins must be involved. Here, we identify MYB-related transcription factors as regulators of chloroplast biogenesis in the liverwort Marchantia polymorpha and angiosperm Arabidopsis thaliana. In both species, double-mutant alleles in MYB-related genes show very limited chloroplast development, and photosynthesis gene expression is perturbed to a greater extent than in GLK mutants. Genes encoding enzymes of chlorophyll biosynthesis are controlled by MYB-related and GLK proteins, whereas those allowing CO2 fixation, photorespiration, and photosystem assembly and repair require MYB-related proteins. Regulation between the MYB-related and GLK transcription factors appears more extensive in A. thaliana than in M. polymorpha. Thus, MYB-related and GLK genes have overlapping as well as distinct targets. We conclude that MYB-related and GLK transcription factors orchestrate chloroplast development in land plants.

DNA methylation epigenetically silences crossover hot spots and controls chromosomal domains of meiotic recombination in Arabidopsis.

During meiosis, homologous chromosomes undergo crossover recombination, which is typically concentrated in narrow hot spots that are controlled by genetic and epigenetic information. Arabidopsis chromosomes are highly DNA methylated in the repetitive centromeres, which are also crossover-suppressed. Here we demonstrate that RNA-directed DNA methylation is sufficient to locally silence Arabidopsis euchromatic crossover hot spots and is associated with increased nucleosome density and H3K9me2. However, loss of CG DNA methylation maintenance in met1 triggers epigenetic crossover remodeling at the chromosome scale, with pericentromeric decreases and euchromatic increases in recombination. We used recombination mutants that alter interfering and noninterfering crossover repair pathways (fancm and zip4) to demonstrate that remodeling primarily involves redistribution of interfering crossovers. Using whole-genome bisulfite sequencing, we show that crossover remodeling is driven by loss of CG methylation within the centromeric regions. Using cytogenetics, we profiled meiotic DNA double-strand break (DSB) foci in met1 and found them unchanged relative to wild type. We propose that met1 chromosome structure is altered, causing centromere-proximal DSBs to be inhibited from maturation into interfering crossovers. These data demonstrate that DNA methylation is sufficient to silence crossover hot spots and plays a key role in establishing domains of meiotic recombination along chromosomes.

Nucleosomes and DNA methylation shape meiotic DSB frequency in Arabidopsis thaliana transposons and gene regulatory regions.

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 Arabidopsis thaliana 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 met1 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 met1 mutants, coincident with reduced nucleosome occupancy, gain of transcription, and H3K4me3. Together, our work reveals a complex relationship between chromatin and meiotic DSBs within A. thaliana genes and transposons, with significance for the diversity and evolution of plant genomes.

Recombination Rate Heterogeneity within Arabidopsis Disease Resistance Genes.

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.

Epigenetic remodeling of meiotic crossover frequency in Arabidopsis thaliana DNA methyltransferase mutants.

Meiosis is a specialized eukaryotic cell division that generates haploid gametes required for sexual reproduction. During meiosis, homologous chromosomes pair and undergo reciprocal genetic exchange, termed crossover (CO). Meiotic CO frequency varies along the physical length of chromosomes and is determined by hierarchical mechanisms, including epigenetic organization, for example methylation of the DNA and histones. Here we investigate the role of DNA methylation in determining patterns of CO frequency along Arabidopsis thaliana chromosomes. In A. thaliana the pericentromeric regions are repetitive, densely DNA methylated, and suppressed for both RNA polymerase-II transcription and CO frequency. DNA hypomethylated methyltransferase1 (met1) mutants show transcriptional reactivation of repetitive sequences in the pericentromeres, which we demonstrate is coupled to extensive remodeling of CO frequency. We observe elevated centromere-proximal COs in met1, coincident with pericentromeric decreases and distal increases. Importantly, total numbers of CO events are similar between wild type and met1, suggesting a role for interference and homeostasis in CO remodeling. To understand recombination distributions at a finer scale we generated CO frequency maps close to the telomere of chromosome 3 in wild type and demonstrate an elevated recombination topology in met1. Using a pollen-typing strategy we have identified an intergenic nucleosome-free CO hotspot 3a, and we demonstrate that it undergoes increased recombination activity in met1. We hypothesize that modulation of 3a activity is caused by CO remodeling driven by elevated centromeric COs. These data demonstrate how regional epigenetic organization can pattern recombination frequency along eukaryotic chromosomes.

Putative Arabidopsis THO/TREX mRNA export complex is involved in transgene and endogenous siRNA biosynthesis.

RNA silencing in plants and some animals has a non-cell-autonomous effect due to an RNA signal that moves between cells or organs. To identify unique factors involved in this process, we analyzed a group of Arabidopsis mutants with defective spread of RNA silencing from a transgene expressed specifically in the phloem. These mutants accumulated reduced amounts of small interfering (si)RNA from the transgene locus and from endogenous loci TAS1, TAS2, and an inverted repeat locus IR71. The defect in TAS1 and TAS2 siRNA biogenesis is in the processing of a long siRNA precursor. We mapped the mutations to a gene encoding the Arabidopsis homolog of a protein, TEX1, which is involved in intracellular transport of RNA in animals. TEX1 is a component of the THO/TREX complex, and we show that the Arabidopsis TEX1 interacts with other predicted components of a THO/TREX complex. Correspondingly, we found at least two other components of the Arabidopsis THO core complex that are involved in RNA silencing. To reconcile the effect of these mutations on transgene and endogenous gene siRNA, we propose a mechanism in which THO/TREX processes or transports a long RNA molecule so that it can be a template for secondary siRNA production.

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

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

Quantification and Sequencing of Crossover Recombinant Molecules from Arabidopsis Pollen DNA.

During meiosis, homologous chromosomes undergo recombination, which can result in formation of reciprocal crossover molecules. Crossover frequency is highly variable across the genome, typically occurring in narrow hotspots, which has a significant effect on patterns of genetic diversity. Here we describe methods to measure crossover frequency in plants at the hotspot scale (bp-kb), using allele-specific PCR amplification from genomic DNA extracted from the pollen of F1 heterozygous plants. We describe (1) titration methods that allow amplification, quantification and sequencing of single crossover molecules, (2) quantitative PCR methods to more rapidly measure crossover frequency, and (3) application of high-throughput sequencing for study of crossover distributions within hotspots. We provide detailed descriptions of key steps including pollen DNA extraction, prior identification of hotspot locations, allele-specific oligonucleotide design, and sequence analysis approaches. Together, these methods allow the rate and recombination topology of plant hotspots to be robustly measured and compared between varied genetic backgrounds and environmental conditions.

Juxtaposition of heterozygous and homozygous regions causes reciprocal crossover remodelling via interference during Arabidopsis meiosis.

During meiosis homologous chromosomes undergo crossover recombination. Sequence differences between homologs can locally inhibit crossovers. Despite this, nucleotide diversity and population-scaled recombination are positively correlated in eukaryote genomes. To investigate interactions between heterozygosity and recombination we crossed Arabidopsis lines carrying fluorescent crossover reporters to 32 diverse accessions and observed hybrids with significantly higher and lower crossovers than homozygotes. Using recombinant populations derived from these crosses we observed that heterozygous regions increase crossovers when juxtaposed with homozygous regions, which reciprocally decrease. Total crossovers measured by chiasmata were unchanged when heterozygosity was varied, consistent with homeostatic control. We tested the effects of heterozygosity in mutants where the balance of interfering and non-interfering crossover repair is altered. Crossover remodeling at homozygosity-heterozygosity junctions requires interference, and non-interfering repair is inefficient in heterozygous regions. As a consequence, heterozygous regions show stronger crossover interference. Our findings reveal how varying homolog polymorphism patterns can shape meiotic recombination.

High-throughput analysis of meiotic crossover frequency and interference via flow cytometry of fluorescent pollen in Arabidopsis thaliana.

During meiosis, reciprocal exchange between homologous chromosomes occurs as a result of crossovers (COs). CO frequency varies within genomes and is subject to genetic, epigenetic and environmental control. As robust measurement of COs is limited by their low numbers, typically 1-2 per chromosome, we adapted flow cytometry for use with Arabidopsis transgenic fluorescent protein-tagged lines (FTLs) that express eCFP, dsRed or eYFP fluorescent proteins in pollen. Segregation of genetically linked transgenes encoding fluorescent proteins of distinct colors can be used to detect COs. The fluorescence of up to 80,000 pollen grains per individual plant can be measured in 10-15 min using this protocol. A key element of CO control is interference, which inhibits closely spaced COs. We describe a three-color assay for the measurement of CO frequency in adjacent intervals and calculation of CO interference. We show that this protocol can be used to detect changes in CO frequency and interference in the fancm zip4 double mutant. By enabling high-throughput measurement of CO frequency and interference, these methods will facilitate genetic dissection of meiotic recombination control.

Arabidopsis meiotic crossover hot spots overlap with H2A.Z nucleosomes at gene promoters.

PRDM9 directs human meiotic crossover hot spots to intergenic sequence motifs, whereas budding yeast hot spots overlap regions of low nucleosome density (LND) in gene promoters. To investigate hot spots in plants, which lack PRDM9, we used coalescent analysis of genetic variation in Arabidopsis thaliana. Crossovers increased toward gene promoters and terminators, and hot spots were associated with active chromatin modifications, including H2A.Z, histone H3 Lys4 trimethylation (H3K4me3), LND and low DNA methylation. Hot spot-enriched A-rich and CTT-repeat DNA motifs occurred upstream and downstream, respectively, of transcriptional start sites. Crossovers were asymmetric around promoters and were most frequent over CTT-repeat motifs and H2A.Z nucleosomes. Pollen typing, segregation and cytogenetic analysis showed decreased numbers of crossovers in the arp6 H2A.Z deposition mutant at multiple scales. During meiosis, H2A.Z forms overlapping chromosomal foci with the DMC1 and RAD51 recombinases. As arp6 reduced the number of DMC1 or RAD51 foci, H2A.Z may promote the formation or processing of meiotic DNA double-strand breaks. We propose that gene chromatin ancestrally designates hot spots within eukaryotes and PRDM9 is a derived state within vertebrates.

SDE5, the putative homologue of a human mRNA export factor, is required for transgene silencing and accumulation of trans-acting endogenous siRNA.

Post-transcriptional gene silencing (PTGS) is a sequence-specific RNA degradation process conserved in fungi, plants and animals. The trigger of the mechanism is double-stranded RNA derived from transgenic or endogenous loci and formed by intra- or inter-molecular interactions of single-stranded RNAs or the action of RNA-dependent RNA polymerases (RDRs). Double-stranded RNA from various sources is processed by one of the four Dicer-like (DCL) proteins in Arabidopsis, and the resulting short RNAs enter into at least four different pathways, one of which involves the production of trans-acting short interfering RNAs (tasiRNAs). We report here a novel gene (SDE5) that is required for transgene silencing and the production of tasiRNAs. Mutation in SDE5 also results in hyper-susceptibility to cucumber mosaic virus but not turnip mosaic virus. However, like RDR6, SDE5 is not involved in inverted repeat-induced transgene silencing or the biogenesis of microRNAs and 24 nt siRNAs produced by DCL3. Based on these results, we propose that SDE5 acts together with RDR6 in generating double-stranded RNA from specific single-stranded RNAs. As the sequence of SDE5 has sequence features shared by TAP, a human mRNA export factor, we propose that its role could be in the transport of RNA molecules that are converted into a double-stranded form by RDR6.