Search

Duckweeds are amongst the fastest growing of higher plants, making them attractive high-biomass targets for biofuel feedstock production. Their fronds have high rates of fatty acid synthesis to meet the demand for new membranes, but triacylglycerols (TAG) only accumulate to very low levels. Here we report on the engineering of Lemna japonica for the synthesis and accumulation of TAG in its fronds. This was achieved by expression of an estradiol-inducible cyan fluorescent protein-Arabidopsis WRINKLED1 fusion protein (CFP-AtWRI1), strong constitutive expression of a mouse diacylglycerol:acyl-CoA acyltransferase2 (MmDGAT), and a sesame oleosin variant (SiOLE(*)). Individual expression of each gene increased TAG accumulation by 1- to 7-fold relative to controls, while expression of pairs of these genes increased TAG by 7- to 45-fold. In uninduced transgenics containing all three genes, TAG accumulation increased by 45-fold to 3.6% of dry weight (DW) without severely impacting growth, and by 108-fold to 8.7% of DW after incubation on medium containing 100 μm estradiol for 4 days. TAG accumulation was accompanied by an increase in total fatty acids of up to three-fold to approximately 15% of DW. Lipid droplets from fronds of all transgenic lines were visible by confocal microscopy of BODIPY-stained fronds. At a conservative 12 tonnes (dry matter) per acre and 10% (DW) TAG, duckweed could produce 350 gallons of oil/acre/year, approximately seven-fold the yield of soybean, and similar to that of oil palm. These findings provide the foundation for optimizing TAG accumulation in duckweed and present a new opportunity for producing biofuels and lipidic bioproducts.

Epigenetic modifications that arise during plant and animal development, such as DNA and histone modification, are mostly reset during gamete formation, but some are inherited from the germline including those marking imprinted genes. Small RNAs guide these epigenetic modifications, and some are also inherited by the next generation. In C. elegans, inherited small RNA precursors have poly (UG) tails, but how inherited small RNAs are distinguished in other animals and plants is unknown. Pseudouridine (Ψ) is the most abundant RNA modification but has not been explored in small RNAs. Here, we develop novel assays to detect Ψ in short RNA sequences, demonstrating its presence in mouse and Arabidopsis microRNAs and their precursors. We also detect substantial enrichment in germline small RNAs, namely epigenetically activated siRNAs (easiRNAs) in Arabidopsis pollen, and piwi-interacting piRNAs in mouse testis. In pollen, pseudouridylated easiRNAs are localized to sperm cells, and we found that PAUSED/HEN5 (PSD), the plant homolog of Exportin-t, interacts genetically with Ψ and is required for transport of easiRNAs into sperm cells from the vegetative nucleus. We further show that Exportin-t is required for the triploid block: chromosome dosage-dependent seed lethality that is epigenetically inherited from pollen. Thus, Ψ has a conserved role in marking inherited small RNAs in the germline.

Plant genomes are largely comprised of retrotransposons which can replicate through ‘copy and paste' mechanisms. Long terminal repeat (LTR) retrotransposons are the major class of retrotransposons in plant species, and importantly they broadly affect the expression of nearby genes. Although most LTR retrotransposons are non-functional, active retrotranspositions have been reported in plant species or mutants under normal growth condition and environmental stresses. With the well-defined reference genome and numerous mutant alleles, Arabidopsis studies have significantly expanded our understanding of retrotransposon regulation. Active LTR retrotransposon loci produce virus-like particles to perform reverse transcription, and their complementary DNA can be inserted into new genomic loci. Due to the detrimental consequences of retrotransposition, plants like animals, have developed transcriptional and post-transcriptional silencing mechanisms. Recently several different genome-wide techniques have been developed to understand LTR retrotransposition in Arabidopsis and different plant species. Transposome, methylome, transcriptome, translatome and small RNA sequencing data have revealed how host silencing mechanisms can affect multiple steps of retrotransposition. These recent advances shed light on future mechanistic studies of retrotransposition as well as retrotransposon diversity.

Epigenetic reprogramming occurs during gametogenesis as well as during embryogenesis to reset the genome for early development. In flowering plants, many heterochromatic marks are maintained in sperm, but asymmetric DNA methylation is mostly lost. Asymmetric DNA methylation is dependent on small RNA but the re-establishment of silencing in embryo is not well understood. Here we demonstrate that small RNAs direct the histone H3 lysine 9 dimethylation during Arabidopsis thaliana embryonic development, together with asymmetric DNA methylation. This de novo silencing mechanism depends on the catalytic domain of SUVH9, a Su(Var)3-9 homolog thought to be catalytically inactive.

Epigenetic marks are reprogrammed in the gametes to reset genomic potential in the next generation. In mammals, paternal chromatin is extensively reprogrammed through the global erasure of DNA methylation and the exchange of histones with protamines1,2. Precisely how the paternal epigenome is reprogrammed in flowering plants has remained unclear since DNA is not demethylated and histones are retained in sperm3,4. Here, we describe a multi-layered mechanism by which H3K27me3 is globally lost from histone-based sperm chromatin in Arabidopsis. This mechanism involves the silencing of H3K27me3 writers, activity of H3K27me3 erasers and deposition of a sperm-specific histone, H3.10 (ref. 5), which we show is immune to lysine 27 methylation. The loss of H3K27me3 facilitates the transcription of genes essential for spermatogenesis and pre-configures sperm with a chromatin state that forecasts gene expression in the next generation. Thus, plants have evolved a specific mechanism to simultaneously differentiate male gametes and reprogram the paternal epigenome.

The selection and firing of DNA replication origins play key roles in ensuring that eukaryotes accurately replicate their genomes. This process is not well documented in plants due in large measure to difficulties in working with plant systems. We developed a new functional assay to label and map very early replicating loci that must, by definition, include at least a subset of replication origins. Arabidopsis (Arabidopsis thaliana) cells were briefly labeled with 5-ethynyl-2′-deoxy-uridine, and nuclei were subjected to two-parameter flow sorting. We identified more than 5500 loci as initiation regions (IRs), the first regions to replicate in very early S phase. These were classified as strong or weak IRs based on the strength of their replication signals. Strong initiation regions were evenly spaced along chromosomal arms and depleted in centromeres, while weak initiation regions were enriched in centromeric regions. IRs are AT-rich sequences flanked by more GC-rich regions and located predominantly in intergenic regions. Nuclease sensitivity assays indicated that IRs are associated with accessible chromatin. Based on these observations, initiation of plant DNA replication shows some similarity to, but is also distinct from, initiation in other well-studied eukaryotic systems.

Although paramutation has been well-studied at a few hallmark loci involved in anthocyanin biosynthesis in maize, the cellular and molecular mechanisms underlying the phenomenon remain largely unknown. Previously described actors of paramutation encode components of the RNA-directed DNA-methylation (RdDM) pathway that participate in the biogenesis of 24-nucleotide small interfering RNAs (24-nt siRNAs) and long non-coding RNAs. In this study, we uncover an ARGONAUTE (AGO) protein as an effector of the RdDM pathway that is in charge of guiding 24-nt siRNAs to their DNA target to create de novo DNA methylation. We combined immunoprecipitation, small RNA sequencing and reverse genetics to, first, validate AGO104 as a member of the RdDM effector complex and, then, investigate its role in paramutation. We found that AGO104 binds 24-nt siRNAs involved in RdDM, including those required for paramutation at the b1 locus. We also show that the ago104-5 mutation causes a partial reversion of the paramutation phenotype at the b1 locus, revealed by intermediate pigmentation levels in stem tissues. Therefore, our results place AGO104 as a new member of the RdDM effector complex that plays a role in paramutation at the b1 locus in maize.

Centromeres attach chromosomes to spindle microtubules during cell division and, despite this conserved role, show paradoxically rapid evolution and are typified by complex repeats. We used ultra-long-read sequencing to generate the Col-CEN Arabidopsis thaliana genome assembly that resolves all five centromeres. The centromeres consist of megabase-scale tandemly repeated satellite arrays, which support high CENH3 occupancy and are densely DNA methylated, with satellite variants private to each chromosome. CENH3 preferentially occupies satellites with least divergence and greatest higher-order repetition. The centromeres are invaded by ATHILA retrotransposons, which disrupt genetic and epigenetic organization of the centromeres. Crossover recombination is suppressed within the centromeres, yet low levels of meiotic DSBs occur that are regulated by DNA methylation. We propose that Arabidopsis centromeres are evolving via cycles of satellite homogenization and retrotransposon-driven diversification.One-sentence summaryLong read sequencing and assembly of the Arabidopsis centromeres reveals their genetic and epigenetic topography.

Centromeres attach chromosomes to spindle microtubules during cell division and, despite this conserved role, show paradoxically rapid evolution and are typified by complex repeats. We used long-read sequencing to generate the Col-CEN

Paramutation is an exception among eukaryotes, in which epigenetic information is conserved through mitosis and meiosis. It has been studied for over 70 years in maize, but the mechanisms involved are largely unknown. Previously described actors of paramutation encode components of the RNA-dependent DNA-methylation (RdDM) pathway all involved in the biogenesis of 24-nt small RNAs. However, no actor of paramutation have been identified in the effector complex of RdDM. Here, through a combination of reverse genetics, immunolocalization and immunoprecipitation (siRNA-IP) we found that ARGONAUTE104 (AGO104), AGO105 and AGO119 are members of the RdDM effector complex in maize and bind siRNAs produced from the tandem repeats required for paramutation at theb1locus. We also showed that AGO104 is an effector of theb1paramutation in maize.Author summaryReprogramming of epigenetic information has been described in both plants and mammals. Here, we show that maizeARGONAUTE (AGO) AGO104andAGO105/AGO119, respectively the close homologs ofA. thaliana AGO9andAGO4, are required to enable paramutation at theb1locus in maize. Paramutation is an epigenetic phenomenon that is stable over many generations (both mitotically and meiotically). A classic example is thebooster1(b1) gene in maize, where the weakly expressedBooster’(B’) allele stably decreases the expression of theBooster-Intense(B-I) allele, and changes it into a newB’allele. This newB’allele will in turn changeB-Iinto newB’in subsequent crosses. Previous research demonstrated that paramutation requires several proteins involved in the biosynthesis of small interfering RNAs (siRNAs) all related to the RNA-dependent DNA-methylation (RdDM) pathway. Yet, few members of the RdDM were functionally identified in maize. Here, we identify two new members of the maize RdDM pathway, and provide evidence that they are also involved in paramutation at theb1locus.

The innate and adaptive immune response are regulated by biological clocks, and circulating lymphocytes are lowest at sunrise. Accordingly, severity of disease in mouse models is highly dependent on the time of day of viral infection. Here, we explore whether circadian immunity contributes significantly to seasonality of respiratory viruses, including influenza and SARS-CoV-2. Susceptibility-Infection-Recovery-Susceptibility (SIRS) models of influenza and SIRS-derived models of COVID-19 suggest that local sunrise time is a better predictor of the basic reproductive number (R0) than climate, even when day length is taken into account. Moreover, these models predict a window of susceptibility when local sunrise time corresponds to the morning commute and contact rate is expected to be high. Counterfactual modeling suggests that retaining daylight savings time in the fall would reduce the length of this window, and substantially reduce seasonal waves of respiratory infections.

Transcriptional silencing by RNAi paradoxically relies on transcription, but how the transition from transcription to silencing is achieved has remained unclear. The Cryptic Loci Regulator complex (CLRC) inSchizosaccharomyces pombeis a cullin-ring E3 ligase required for silencing that is recruited by RNAi. We found that the E2 ubiquitin conjugating enzyme Ubc4 interacts with CLRC and mono-ubiquitinates the histone H3K9 methyltransferase Clr4SUV39H1, promoting the transition from co-transcriptional gene silencing (H3K9me2) to transcriptional gene silencing (H3K9me3). Ubiquitination of Clr4 occurs in an intrinsically disordered region (IDR), which undergoes robust liquid-liquid phase separation (LLPS), along with Swi6HP1the effector of transcriptional gene silencing. Phase separation of Clr4 and Swi6 is exquisitely sensitive to non-coding RNA (ncRNA), which promotes dimerization, chromatin association, and di-, but not tri-methylation instead. Ubc4-CLRC also targets the transcriptional co-activator Bdf2BRD4, down-regulating centromeric transcription and small RNA production. The deubiquitinase Ubp3 counteracts both activities.

Small RNA molecules can target a particular virus, gene, or transposable element (TE) with a high degree of specificity. Their ability to move from cell to cell and recognize targets in trans also allows building networks capable of regulating a large number of related targets at once. In the case of epigenetic silencing, small RNA may use the widespread distribution of TEs in eukaryotic genomes to coordinate many loci across developmental and generational time. Here, we discuss the intriguing role of plant small RNA in targeting transposons and repeats in pollen and seeds. Epigenetic reprogramming in the germline and early seed development provides a mechanism to control genome dosage, imprinted gene expression, and incompatible hybridizations via the "triploid block."

Rootless plants in the genus Wolffia are some of the fastest growing known plants on Earth. Wolffia have a reduced body plan, primarily multiplying through a budding type of asexual reproduction. Here, we generated draft reference genomes for Wolffia australiana (Benth.) Hartog & Plas, which has the smallest genome size in the genus at 357 Mb and has a reduced set of predicted protein-coding genes at about 15,000. Comparison between multiple high-quality draft genome sequences from W. australiana clones confirmed loss of several hundred genes that are highly conserved among flowering plants, including genes involved in root developmental and light signaling pathways. Wolffia has also lost most of the conserved nucleotide-binding leucine-rich repeat (NLR) genes that are known to be involved in innate immunity, as well as those involved in terpene biosynthesis, while having a significant overrepresentation of genes in the sphingolipid pathways that may signify an alternative defense system. Diurnal expression analysis revealed that only 13% of Wolffia genes are expressed in a time-of-day (TOD) fashion, which is less than the typical ∼40% found in several model plants under the same condition. In contrast to the model plants Arabidopsis and rice, many of the pathways associated with multicellular and developmental processes are not under TOD control in W. australiana, where genes that cycle the conditions tested predominantly have carbon processing and chloroplast-related functions. The Wolffia genome and TOD expression data set thus provide insight into the interplay between a streamlined plant body plan and optimized growth.

Wolffia is the fastest growing plant genus on Earth with a recorded doubling time of less than a day. Wolffia has a dramatically reduced body plan, primarily growing through a continuous, budding-type asexual reproduction with no obvious phase transition. Most plants are bound by the 24-hour light-dark cycle with the majority of processes such as gene expression partitioned or phased to a specific time-of-day (TOD). However, the role that TOD information and the circadian clock plays in facilitating the growth of a fast-growing plant is unknown. Here we generated draft reference genomes for Wolffia australiana (Benth.) Hartog & Plas to monitor gene expression over a two-day time course under light-dark cycles. Wolffia australiana has the smallest genome size in the genus at 357 Mb and has a dramatically reduced gene set at 15,312 with a specific loss of root (WOX5), vascular (CASP), circadian (TOC1), and light-signaling (NPH3) genes. Remarkably, it has also lost all but one of the NLR genes that are known to be involved in innate immunity. In addition, only 13% of its genes cycle, which is far less than in other plants, with an overrepresentation of genes associated with carbon processing and chloroplast-related functions. Despite having a focused set of cycling genes, TOD cis-elements are conserved in W. australiana, consistent with the overall conservation of transcriptional networks. In contrast to the model plants Arabidopsis thaliana and Oryza sativa, the reduction in cycling genes correlates with fewer pathways under TOD control in Wolffia, which could reflect a release of functional gating. Since TOD networks and the circadian clock work to gate activities to specific times of day, this minimization of regulation may enable Wolffia to grow continuously with optimal economy. Wolffia is an ideal model to study the transcriptional control of growth and the findings presented here could serve as a template for plant improvement.

In Arabidopsis, LTR retrotransposons are activated by mutations in the chromatin gene DECREASE in DNA METHYLATION 1 (DDM1), giving rise to 21- to 22-nt epigenetically activated siRNA (easiRNA) that depend on RNA DEPENDENT RNA POLYMERASE 6 (RDR6). We purified virus-like particles (VLPs) from ddm1 and ddm1rdr6 mutants in which genomic RNA is reverse transcribed into complementary DNA. High-throughput short-read and long-read sequencing of VLP DNA (VLP DNA-seq) revealed a comprehensive catalog of active LTR retrotransposons without the need for mapping transposition, as well as independent of genomic copy number. Linear replication intermediates of the functionally intact COPIA element EVADE revealed multiple central polypurine tracts (cPPTs), a feature shared with HIV in which cPPTs promote nuclear localization. For one member of the ATCOPIA52 subfamily (SISYPHUS), cPPT intermediates were not observed, but abundant circular DNA indicated transposon “suicide” by auto-integration within the VLP. easiRNA targeted EVADE genomic RNA, polysome association of GYPSY (ATHILA) subgenomic RNA, and transcription via histone H3 lysine-9 dimethylation. VLP DNA-seq provides a comprehensive landscape of LTR retrotransposons and their control at transcriptional, post-transcriptional, and reverse transcriptional levels.

MaizeCODE is a project aimed at identifying and analyzing functional elements in the maize genome. In its initial phase, MaizeCODE assayed up to five tissues from four maize strains (B73, NC350, W22, TIL11) by RNA-Seq, Chip-Seq, RAMPAGE, and small RNA sequencing. To facilitate reproducible science and provide both human and machine access to the MaizeCODE data, we enhanced SciApps, a cloud-based portal, for analysis and distribution of both raw data and analysis results. Based on the SciApps workflow platform, we generated new components to support the complete cycle of MaizeCODE data management. These include publicly accessible scientific workflows for the reproducible and shareable analysis of various functional data, a RESTful API for batch processing and distribution of data and metadata, a searchable data page that lists each MaizeCODE experiment as a reproducible workflow, and integrated JBrowse genome browser tracks linked with workflows and metadata. The SciApps portal is a flexible platform that allows the integration of new analysis tools, workflows, and genomic data from multiple projects. Through metadata and a ready-to-compute cloud-based platform, the portal experience improves access to the MaizeCODE data and facilitates its analysis.

RNA interference (RNAi), a cellular process through which small RNAs target and regulate complementary RNA transcripts, has well-characterized roles in post-transcriptional gene regulation and transposon repression. Recent studies have revealed additional conserved roles for RNAi proteins, such as Argonaute and Dicer, in chromosome function. By guiding chromatin modification, RNAi components promote chromosome segregation during both mitosis and meiosis and regulate chromosomal and genomic dosage response. Small RNAs and the RNAi machinery also participate in the resolution of DNA damage. Interestingly, many of these lesser-studied functions seem to be more strongly conserved across eukaryotes than are well-characterized functions such as the processing of microRNAs. These findings have implications for the evolution of RNAi since the last eukaryotic common ancestor, and they provide a more complete view of the functions of RNAi.

In Arabidopsis, LTR-retrotransposons are activated by mutations in the chromatin remodeler DECREASE in DNA METHYLATION 1 (DDM1), giving rise to 21-22nt epigenetically activated siRNAs (easiRNAs) that depend on RNA DEPENDENT RNA POLYMERASE 6 (RDR6). We purified virus-like-particles (VLPs) from ddm1 and ddm1rdr6 mutants in which genomic RNA is reverse transcribed into complementary DNA. Next generation short-read and long-read sequencing of VLP DNA (VLP DNA-seq) revealed a comprehensive catalog of active LTR-retrotransposons without the need for mapping transposition, and independent of genomic copy number. Linear replication intermediates of a functionally intact copia element EVADE revealed multiple central polypurine tracts (cPPT), a feature shared with HIV where cPPT promote nuclear localization. For one member of the ATCOPIA52 subfamily (SISYPHUS), cPPT intermediates were not observed, but abundant circular DNA indicated transposon “suicide” by auto-integration within the VLP. easiRNA targeted EVADE genomic RNA, polysome association of GYPSY (ATHILA) subgenomic RNA, and transcription via histone H3 lysine-9 dimethylation. VLP DNA-seq provides a comprehensive landscape of LTR-retrotransposons, and their control at transcriptional, post-transcriptional and reverse transcriptional levels.

Plant cells undergo two types of cell cycles–the mitotic cycle in which DNA replication is coupled to mitosis, and the endocycle in which DNA replication occurs in the absence of cell division. To investigate DNA replication programs in these two types of cell cycles, we pulse labeled intact root tips of maize (Zea mays) with 5-ethynyl-2’-deoxyuridine (EdU) and used flow sorting of nuclei to examine DNA replication timing (RT) during the transition from a mitotic cycle to an endocycle. Comparison of the sequence-based RT profiles showed that most regions of the maize genome replicate at the same time during S phase in mitotic and endocycling cells, despite the need to replicate twice as much DNA in the endocycle and the fact that endocycling is typically associated with cell differentiation. However, regions collectively corresponding to 2% of the genome displayed significant changes in timing between the two types of cell cycles. The majority of these regions are small with a median size of 135 kb, shift to a later RT in the endocycle, and are enriched for genes expressed in the root tip. We found larger regions that shifted RT in centromeres of seven of the ten maize chromosomes. These regions covered the majority of the previously defined functional centromere, which ranged between 1 and 2 Mb in size in the reference genome. They replicate mainly during mid S phase in mitotic cells but primarily in late S phase of the endocycle. In contrast, the immediately adjacent pericentromere sequences are primarily late replicating in both cell cycles. Analysis of CENH3 enrichment levels in 8C vs 2C nuclei suggested that there is only a partial replacement of CENH3 nucleosomes after endocycle replication is complete. The shift to later replication of centromeres and possible reduction in CENH3 enrichment after endocycle replication is consistent with a hypothesis that centromeres are inactivated when their function is no longer needed., Author summary In traditional cell division, or mitosis, a cell’s genetic material is duplicated and then split between two daughter cells. In contrast, in some specialized cell types, the DNA is duplicated a second time without an intervening division step, resulting in cells that carry twice as much DNA. This phenomenon, which is called the endocycle, is common during plant development. At each step, DNA replication follows an ordered program in which highly compacted DNA is unraveled and replicated in sections at different times during the synthesis (S) phase. In plants, it is unclear whether traditional and endocycle programs are the same, especially since endocycling cells are typically in the process of differentiation. Using root tips of maize, we found that in comparison to replication in the mitotic cell cycle, there is a small portion of the genome whose replication in the endocycle is shifted in time, usually to later in S phase. Some of these regions are scattered around the genome and mostly coincide with active genes. However, the most prominent shifts occur in centromeres. The shift to later replication in centromeres is noteworthy because they orchestrate the process of separating duplicated chromosomes into daughter cells, a function that is not needed in the endocycle.

Unlike in Arabidopsis, loss of function of Polymerase IV in the Brassicaceae species Capsella rubella causes microspore arrest, revealing an important role for Polymerase IV in pollen development.In Arabidopsis (Arabidopsis thaliana), DNA-dependent RNA polymerase IV (Pol IV) is required for the formation of transposable element (TE)-derived small RNA transcripts. These transcripts are processed by DICER-LIKE3 into 24-nucleotide small interfering RNAs (siRNAs) that guide RNA-directed DNA methylation. In the pollen grain, Pol IV is also required for the accumulation of 21/22-nucleotide epigenetically activated siRNAs, which likely silence TEs via post-transcriptional mechanisms. Despite this proposed role of Pol IV, its loss of function in Arabidopsis does not cause a discernible pollen defect. Here, we show that the knockout of NRPD1, encoding the largest subunit of Pol IV, in the Brassicaceae species Capsella (Capsella rubella), caused postmeiotic arrest of pollen development at the microspore stage. As in Arabidopsis, all TE-derived siRNAs were depleted in Capsella nrpd1 microspores. In the wild-type background, the same TEs produced 21/22-nucleotide and 24-nucleotide siRNAs; these processes required Pol IV activity. Arrest of Capsella nrpd1 microspores was accompanied by the deregulation of genes targeted by Pol IV-dependent siRNAs. TEs were much closer to genes in Capsella compared with Arabidopsis, perhaps explaining the essential role of Pol IV in pollen development in Capsella. Our discovery that Pol IV is functionally required in Capsella microspores emphasizes the relevance of investigating different plant models.

tRNA fragments (tRFs) are a class of small, regulatory RNAs with diverse functions. 3'-Derived tRFs perfectly match long terminal repeat (LTR)-retroelements which use the 3'-end of tRNAs to prime reverse transcription. Recent work has shown that tRFs target LTR-retroviruses and -transposons for the RNA interference (RNAi) pathway and also inhibit mobility by blocking reverse transcription. The highly conserved tRNA primer binding site (PBS) in LTR-retroelements is a unique target for 3'-tRFs to recognize and block abundant but diverse LTR-retrotransposons that become transcriptionally active during epigenetic reprogramming in development and disease. 3'-tRFs are processed from full-length tRNAs under so far unknown conditions and potentially protect many cell types. tRFs appear to be an ancient link between RNAi, transposons, and genome stability.

Eukaryotes use a temporally regulated process, known as the replication timing program, to ensure that their genomes are fully and accurately duplicated during S phase. Replication timing programs are predictive of genomic features and activity and are considered to be functional readouts of chromatin organization. Although replication timing programs have been described for yeast and animal systems, much less is known about the temporal regulation of plant DNA replication or its relationship to genome sequence and chromatin structure. We used the thymidine analog, 5-ethynyl-2'-deoxyuridine, in combination with flow sorting and Repli-Seq to describe, at high-resolution, the genome-wide replication timing program for Arabidopsis (Arabidopsis thaliana) Col-0 suspension cells. We identified genomic regions that replicate predominantly during early, mid, and late S phase, and correlated these regions with genomic features and with data for chromatin state, accessibility, and long-distance interaction. Arabidopsis chromosome arms tend to replicate early while pericentromeric regions replicate late. Early and mid-replicating regions are gene-rich and predominantly euchromatic, while late regions are rich in transposable elements and primarily heterochromatic. However, the distribution of chromatin states across the different times is complex, with each replication time corresponding to a mixture of states. Early and mid-replicating sequences interact with each other and not with late sequences, but early regions are more accessible than mid regions. The replication timing program in Arabidopsis reflects a bipartite genomic organization with early/mid-replicating regions and late regions forming separate, noninteracting compartments. The temporal order of DNA replication within the early/mid compartment may be modulated largely by chromatin accessibility.

Sultana et al. (2019) and Flasch et al. (2019) determined integration patterns of human LINE-1 (long interspersed element-1) retrotransposons highlighting their interaction with DNA replication guided by their 5'-TTTT/AA-3' integration motif and nucleotide biases in the genome.

Transposable elements are mobile genetic elements that are prevalent in plant genomes and are silenced by epigenetic modification. Different epigenetic modification pathways play distinct roles in the control of transposable element transcription, replication and recombination. The Arabidopsis genome contains families of all of the major transposable element classes, which are differentially enriched in particular genomic regions. Whole genome sequencing and DNA methylation profiling of hundreds of natural Arabidopsis accessions has revealed that transposable elements exhibit significant intraspecific genetic and epigenetic variation, and that genetic variation often underlies epigenetic variation. Together, epigenetic modification and the forces of selection define the scope within which transposable elements can contribute to, and control, genome evolution.

In Arabidopsis (Arabidopsis thaliana), DNA-dependent RNA polymerase IV (Pol IV) is required for the formation of transposable element (TE)-derived small RNA transcripts. These transcripts are processed by DICER-LIKE3 into 24-nucleotide small interfering RNAs (siRNAs) that guide RNA-directed DNA methylation. In the pollen grain, Pol IV is also required for the accumulation of 21/22-nucleotide epigenetically activated siRNAs, which likely silence TEs via post-transcriptional mechanisms. Despite this proposed role of Pol IV, its loss of function in Arabidopsis does not cause a discernible pollen defect. Here, we show that the knockout of NRPD1, encoding the largest subunit of Pol IV, in the Brassicaceae species Capsella (Capsella rubella), caused postmeiotic arrest of pollen development at the microspore stage. As in Arabidopsis, all TE-derived siRNAs were depleted in Capsella nrpd1 microspores. In the wild-type background, the same TEs produced 21/22-nucleotide and 24-nucleotide siRNAs; these processes required Pol IV activity. Arrest of Capsella nrpd1 microspores was accompanied by the deregulation of genes targeted by Pol IV-dependent siRNAs. TEs were much closer to genes in Capsella compared with Arabidopsis, perhaps explaining the essential role of Pol IV in pollen development in Capsella. Our discovery that Pol IV is functionally required in Capsella microspores emphasizes the relevance of investigating different plant models.

In Arabidopsis thaliana, the DNA-dependent RNA polymerase IV (Pol IV) is required for the formation of transposable element (TE)-derived small RNA (sRNA) transcripts. These transcripts are processed by DICER-LIKE 3 into 24-nt small interfering RNAs (siRNAs) that guide RNA-dependent DNA methylation. In the pollen grain, Pol IV is also required for the accumulation of 21/22-nt epigenetically-activated siRNAs (easiRNAs) that likely silence TEs by post-transcriptional mechanisms. Despite this proposed functional role, loss of Pol IV function in Arabidopsis does not cause a discernable pollen defect. Here, we show that loss of NRPD1, encoding the largest subunit of Pol IV in the Brassicaceae Capsella rubella, causes post-meiotic arrest of pollen development at the microspore stage. As in Arabidopsis, all TE-derived siRNAs were depleted in Capsella nrpd1 microspores. In wild-type background, we found that the same TEs produced 21/22-nt and 24-nt siRNAs, leading us to propose that Pol IV is generating the direct precursors for 21-24-nt siRNAs, which are targeted by different DICERs. Arrest of Capsella nrpd1 microspores was accompanied by deregulation of genes targeted by Pol IV-dependent siRNAs. The distance of TEs to genes was much closer in Capsella rubella compared to Arabidopsis thaliana, providing a possible explanation for the essential role of Pol IV for pollen development in Capsella. Our study in Capsella uncovers a functional requirement of Pol IV in microspores, emphasizing the relevance of investigating different plant models.One-sentence summaryLoss of Polymerase IV function in Capsella rubella causes microspore arrest, revealing an important functional role of Polymerase IV during pollen development.The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantcell.org) is: Claudia Kohler (Claudia.Kohler@slu.se)

Heterochromatic regions of the genome are epigenetically regulated to maintain a heritable '"silent state"'. In fission yeast and other organisms, epigenetic silencing is guided by nascent transcripts, which are targeted by the RNA interference pathway. The key effector complex of the RNA interference pathway consists of small interfering RNA molecules (siRNAs) associated with Argonaute, assembled into the RNA-induced transcriptional silencing (RITS) complex. This review focuses on our current understanding of how RITS promotes heterochromatin formation, and in particular on the role of Argonaute-containing complexes in many other functions such as quelling, release of RNA polymerases, cellular quiescence and genome defense.

International audience; 5-methyl cytosine is widespread in plant genomes in both CG and non-CG contexts. During replication, hemi-methylation on parental DNA strands guides symmetric CG methylation on nascent strands, but non-CG methylation requires modified histones and small RNA guides. Here, we used immortalized Arabidopsis cell suspensions to sort replicating nuclei and determine genome-wide cytosine methylation dynamics during the plant cell cycle. We find that symmetric mCG and mCHG are selectively retained in actively dividing cells in culture, whereas mCHH is depleted. mCG becomes transiently asymmetric during S phase but is rapidly restored in G2, whereas mCHG remains asymmetric throughout the cell cycle. Hundreds of loci gain ectopic CHG methylation, as well as 24-nt small interfering RNAs (siRNAs) and histone H3 lysine dimethylation (H3K9me2), without gaining CHH methylation. This suggests that spontaneous epialleles that arise in plant cell cultures are stably maintained by siRNA and H3K9me2 independent of the canonical RNA-directed DNA methylation (RdDM) pathway. In contrast, loci that fail to produce siRNA may be targeted for demethylation when the cell cycle arrests. Comparative analysis with methylomes of various tissues and cell types suggests that loss of small-RNA-directed non-CG methylation during DNA replication promotes germline reprogramming and epigenetic variation in plants propagated as clones.

Summary Chromosomal inversions can provide windows onto the cytogenetic, molecular, evolutionary and demographic histories of a species. Here we investigate a paracentric 1.17‐Mb inversion on chromosome 4 of Arabidopsis thaliana with nucleotide precision of its borders. The inversion is created by Vandal transposon activity, splitting an F‐box and relocating a pericentric heterochromatin segment in juxtaposition with euchromatin without affecting the epigenetic landscape. Examination of the RegMap panel and the 1001 Arabidopsis genomes revealed more than 170 inversion accessions in Europe and North America. The SNP patterns revealed historical recombinations from which we infer diverse haplotype patterns, ancient introgression events and phylogenetic relationships. We find a robust association between the inversion and fecundity under drought. We also find linkage disequilibrium between the inverted region and the early flowering Col‐FRIGIDA allele. Finally, SNP analysis elucidates the origin of the inversion to South‐Eastern Europe approximately 5000 years ago and the FRI‐Col allele to North‐West Europe, and reveals the spreading of a single haplotype to North America during the 17th to 19th century. The ‘American haplotype’ was identified from several European localities, potentially due to return migration., Significance Statement Structural rearrangements redefine chromosomes, shape genome diversity and can have profound effects on selection, adaptation and spread. Here we elucidate the history of a paracentric inversion in Arabidopsis thaliana, including its origin a few thousand years ago, its maintenance under certain environmental conditions and its migration patterns, from Europe to North America and back.

DNA methylation is an epigenetic mark that ensures silencing of transposable elements (TEs) and affects gene expression in many organisms. The function of different DNA methylation regulatory pathways has been largely characterized in the model plantArabidopsis thaliana. However, far less is known about DNA methylation regulation and functions in basal land plants. Here we focus on the liverwortMarchantia polymorpha, an emerging model species that represents a basal lineage of land plants. We identified MpMET, theM. polymorphaorthologue of theMETHYLTRANSFERASE 1(MET1) gene required for maintenance of methylation at CG sites in angiosperms. We generated Mpmetmutants using the CRISPR/Cas9 system, which showed a significant loss of CG methylation and severe morphological changes and developmental defects. The mutants developed many adventitious shoot-like structures, suggesting that MpMETis required for maintaining differentiated cellular identities in the gametophyte. Numerous TEs were up-regulated, even though non-CG methylation was highly increased at TEs in the Mpmetmutants. Closer inspection of CHG methylation revealed features unique toM. polymorpha. Methylation of CCG sites inM. polymorphadoes not depend onMET1, unlike inA. thalianaandPhyscomitrella patens. Furthermore, unlikeA. thaliana,M. polymorphashows higher methylation level at CAG sites than at other CHG contexts and CAG/CTG sites are mostly methylated asymmetrically. Interestingly, CAG and CTG methylation reached comparable levels and symmetry upon loss of CG methylation. Our results highlight the diversity of non-CG methylation regulatory mechanisms in plants.

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.

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 DNA replication, chromatin is disrupted ahead of the replication fork, and epigenetic information must be restored behind the fork. How epigenetic marks are inherited through DNA replication remains poorly understood. Histone H3 lysine 9 (H3K9) methylation and histone hypoacetylation are conserved hallmarks of heterochromatin. We previously showed that the inheritance of H3K9 methylation during DNA replication depends on the catalytic subunit of DNA polymerase epsilon, Cdc20. Here we show that the histone-fold subunit of Pol epsilon, Dpb4, interacts an uncharacterized small histone-fold protein, SPCC16C4.22, to form a heterodimer in fission yeast. We demonstrate that SPCC16C4.22 is nonessential for viability and corresponds to the true ortholog of Dpb3. We further show that the Dpb3-Dpb4 dimer associates with histone deacetylases, chromatin remodelers, and histones and plays a crucial role in the inheritance of histone hypoacetylation in heterochromatin. We solve the 1.9-A crystal structure of Dpb3-Dpb4 and reveal that they form the H2A-H2B-like dimer. Disruption of Dpb3-Dpb4 dimerization results in loss of heterochromatin silencing. Our findings reveal a link between histone deacetylation and H3K9 methylation and suggest a mechanism for how two processes are coordinated during replication. We propose that the Dpb3-Dpb4 heterodimer together with Cdc20 serves as a platform for the recruitment of chromatin modifiers and remodelers that mediate heterochromatin assembly during DNA replication, and ensure the faithful inheritance of epigenetic marks in heterochromatin.

Chromosome dosage plays a significant role in reproductive isolation and speciation in both plants and animals, but underlying mechanisms are largely obscure1. Transposable elements can promote hybridity through maternal small RNA2, and have been postulated to regulate dosage response via neighboring imprinted genes3,4. Here, we show that a highly conserved microRNA in plants, miR845, targets the tRNAMetprimer-binding site (PBS) of LTR-retrotransposons inArabidopsispollen, and triggers the accumulation of 21 to 22-nucleotide small RNA in a dose dependent fashion via RNA polymerase IV. We show that these epigenetically activated small-interfering RNAs (easiRNAs) mediate hybridization barriers between diploid seed parents and tetraploid pollen parents (“the triploid block”), and that natural variation for miR845 may account for “endosperm balance” allowing formation of triploid seeds. Targeting the PBS with small RNA is a common mechanism for transposon control in mammals and plants, and provides a uniquely sensitive means to monitor chromosome dosage and imprinting in the developing seed.