Ancient Origin and Recent Innovations of RNA Polymerase IV and V

An RdDM-independent function of Pol V transcripts in gene regulation and plant defence

RNA polymerase V (Pol V) and Pol IV are known to be specialized for RNA-directed DNA methylation (RdDM). Here we report that Pol V, but not Pol IV, regulates hundreds of genes including jasmonic acid-responsive genes and confers plant defence to Botrytis cinerea and Spodoptera exigua. About half of the Pol V-regulated genes are associated with Pol V transcripts (PVTs). We thus hypothesized that some PVTs could regulate gene expression in an RdDM-independent manner. To test this hypothesis, we studied three PVTs, PVT-ERF5a/b and PVT-ERF6, as models. PVT-ERF5a/b and PVT-ERF6 are transcribed from the upstream regions of ERF5 and ERF6 and positively regulate their transcription, thereby regulating plant defence. Such regulation involves PVT-dependent H3K4me3 deposition and requires the DRD1-DMS3-RDM1 complex that mediates Pol V recruitment to the target loci. These findings highlight an unprecedented role for PVTs in regulating gene transcription, apart from serving as scaffold RNAs to direct DNA methylation.

A critical role of a plant-specific TFIIB-related protein, BRP1, in salicylic acid-mediated immune response

An integral part of plant immunity is transcription reprogramming by concerted action of specific transcription factors that activate or repress genes through recruitment or release of RNA polymerase II (Pol II). Pol II is assembled into Pol II holoenzyme at the promoters through association with a group of general transcription factors including transcription factor IIB (TFIIB) to activate transcription. Unlike other eukaryotic organisms, plants have a large family of TFIIB-related proteins with 15 members in Arabidopsis including several plant-specific TFIIB-related proteins (BRPs). Molecular genetic analysis has revealed important roles of some BRPs in plant reproductive processes. In this study, we report that Arabidopsis knockout mutants for BRP1, the founding member of the BRP protein family, were normal in growth and development, but were hypersusceptible to the bacterial pathogen Psuedomonas syringae. The enhanced susceptibility of the brp1 mutants was associated with reduced expression of salicylic acid (SA) biosynthetic gene ISOCHORISMATE SYNTHASE 1 (ICS1) and SA-responsive PATHOGENESIS-RELATED (PR) genes. Pathogen-induced SA accumulation was reduced in the brp1 mutants and exogenous SA rescued the brp1 mutants for resistance to the bacterial pathogen. In uninfected plants, BRP1 was primarily associated with the plastids but pathogen infection induced its accumulation in the nucleus. BRP1 acted as a transcription activator in plant cells and binded to the promoter of ICS1. These results collectively indicate that BRP1 is a functionally specialized transcription factor that increasingly accumulates in the nucleus in response to pathogen infection to promote defense gene expression.

Charophytic Green Algae Encode Ancestral Polymerase IV/Polymerase V Subunits and a CLSY/DRD1 Homolog

In flowering plants, euchromatic transposons are transcriptionally silenced by RNA-directed DNA Methylation, a small RNA-guided de novo methylation pathway. RNA-directed DNA Methylation requires the activity of the RNA Polymerases IV and V, which produce small RNA precursors and noncoding targets of small RNAs, respectively. These polymerases are distinguished from Polymerase II by multiple plant-specific paralogous subunits. Most RNA-directed DNA Methylation components are present in all land plants, and some have been found in the charophytic green algae, a paraphyletic group that is sister to land plants. However, the evolutionary origin of key RNA-directed DNA Methylation components, including the two largest subunits of Polymerase IV and Polymerase V, remains unclear. Here, we show that multiple lineages of charophytic green algae encode a single-copy precursor of the largest subunits of Polymerase IV and Polymerase V, resolving the two presumed duplications in this gene family. We further demonstrate the presence of a Polymerase V-like C-terminal domain, suggesting that the earliest form of RNA-directed DNA Methylation utilized a single Polymerase V-like polymerase. Finally, we reveal that charophytic green algae encode a single CLSY/DRD1-type chromatin remodeling protein, further supporting the presence of a single specialized polymerase in charophytic green algae.

SPT6L, a newly discovered ancestral component of the plant RNA-directed DNA methylation pathway

RNA-directed DNA methylation (RdDM) is driven by small RNAs (sRNAs) complementary to the nascent transcript of RNA polymerase V (Pol V). sRNAs associated with ARGONAUTE (AGO) proteins are tethered to Pol V mainly by the AGO-hook domain of its subunit NRPE1. We found, by in silico analyses, that Pol V strongly colocalizes on chromatin with another AGO-hook protein, SPT6-like (SPT6L), which is a known essential transcription elongation factor of Pol II. Our phylogenetic analysis revealed that SPT6L acquired its AGO-binding capacity already in the most basal streptophyte algae, even before the emergence of Pol V, suggesting that SPT6L might be a driving force behind the RdDM evolution. Since its emergence, SPT6L with the AGO-hook represents the only conserved SPT6 homolog in Viridiplantae, implying that the same protein is involved in both Pol II and Pol V complexes. To better understand the role of SPT6L in the Pol V complex, we characterized genomic loci where these two colocalize and uncovered that DNA methylation there is more dynamic, driven by higher levels of sRNAs often from non-canonical RdDM pathways and more dependent on chromatin modifying and remodeling proteins like MORC. Pol V loci with SPT6L are highly depleted in helitrons but enriched in gene promoters for which locally and temporally precise methylation is necessary. In view of these results, we discuss potential roles of multiple AGO-hook domains present in the Pol V complex and speculate that SPT6L mediates de novo methylation of naïve loci by interconnecting Pol II and Pol V activities.

Molecular mechanisms of the RNA polymerases in plant RNA-directed DNA methylation

In plants, two atypical DNA-dependent RNA polymerases, RNA polymerase IV (Pol IV) and Pol V, and an RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) together produce noncoding RNAs (ncRNAs) to guide the plant-specific RNA-directed DNA methylation (RdDM). Although both Pol IV and Pol V have evolved from the canonical Pol II, they have adapted to different roles in RdDM. The mechanisms of their adaptation are key to understanding plant DNA methylation and the divergent evolution of polymerases. In this review, we summarize insights that have emerged from recent structural studies of Pol IV, Pol V, and RDR2 and discuss their structural features critical for efficient ncRNA production in RdDM.

CLSY docking to Pol IV requires a conserved domain critical for small RNA biogenesis and transposon silencing

Eukaryotes must balance the need for gene transcription by RNA polymerase II (Pol II) against the danger of mutations caused by transposable element (TE) proliferation. In plants, these gene expression and TE silencing activities are divided between different RNA polymerases. Specifically, RNA polymerase IV (Pol IV), which evolved from Pol II, transcribes TEs to generate small interfering RNAs (siRNAs) that guide DNA methylation and block TE transcription by Pol II. While the Pol IV complex is recruited to TEs via SNF2-like CLASSY (CLSY) proteins, how Pol IV partners with the CLSYs remains unknown. Here we identified a conserved CYC-YPMF motif that is specific to Pol IV and is positioned on the complex exterior. Furthermore, we found that this motif is essential for the co-purification of all four CLSYs with Pol IV, but that only one CLSY is present in any given Pol IV complex. These findings support a "one CLSY per Pol IV" model where the CYC-YPMF motif acts as a CLSY-docking site. Indeed, mutations in and around this motif phenocopy pol iv null mutants. Together, these findings provide structural and functional insights into a critical protein feature that distinguishes Pol IV from other RNA polymerases, allowing it to promote genome stability by targeting TEs for silencing.

RNA Polymerases IV and V Are Involved in Olive Fruit Development

Transcription is carried out in most eukaryotes by three multimeric complexes (RNA polymerases I, II and III). However, plants contain two additional RNA polymerases (IV and V), which have evolved from RNA polymerase II. RNA polymerases II, IV and V contain both common and specific subunits that may specialise some of their functions. In this study, we conducted a search for the genes that putatively code for the specific subunits of RNA polymerases IV and V, as well as those corresponding to RNA polymerase II in olive trees. Based on the homology with the genes of Arabidopsis thaliana, we identified 13 genes that putatively code for the specific subunits of polymerases IV and V, and 16 genes that code for the corresponding specific subunits of polymerase II in olives. The transcriptomic analysis by RNA-Seq revealed that the expression of the RNA polymerases IV and V genes was induced during the initial stages of fruit development. Given that RNA polymerases IV and V are involved in the transcription of long non-coding RNAs, we investigated their expression and observed relevant changes in the expression of this type of RNAs. Particularly, the expression of the intergenic and intronic long non-coding RNAs tended to increase in the early steps of fruit development, suggesting their potential role in this process. The positive correlation between the expression of RNA polymerases IV and V subunits and the expression of non-coding RNAs supports the hypothesis that RNA polymerases IV and V may play a role in fruit development through the synthesis of this type of RNAs.

Conserved and non-conserved RNA-target modules in plants: lessons for a better understanding of Marchantia development

A wide variety of functional regulatory non-coding RNAs (ncRNAs) have been identified as essential regulators of plant growth and development. Depending on their category, ncRNAs are not only involved in modulating target gene expression at the transcriptional and post-transcriptional levels but also are involved in processes like RNA splicing and RNA-directed DNA methylation. To fulfill their molecular roles properly, ncRNAs must be precisely processed by multiprotein complexes. In the case of small RNAs, DICER-LIKE (DCL) proteins play critical roles in the production of mature molecules. Land plant genomes contain at least four distinct classes of DCL family proteins (DCL1-DCL4), of which DCL1, DCL3 and DCL4 are also present in the genomes of bryophytes, indicating the early divergence of these genes. The liverwort Marchantia polymorpha has become an attractive model species for investigating the evolutionary history of regulatory ncRNAs and proteins that are responsible for ncRNA biogenesis. Recent studies on Marchantia have started to uncover the similarities and differences in ncRNA production and function between the basal lineage of bryophytes and other land plants. In this review, we summarize findings on the essential role of regulatory ncRNAs in Marchantia development. We provide a comprehensive overview of conserved ncRNA-target modules among M. polymorpha, the moss Physcomitrium patens and the dicot Arabidopsis thaliana, as well as Marchantia-specific modules. Based on functional studies and data from the literature, we propose new connections between regulatory pathways involved in Marchantia's vegetative and reproductive development and emphasize the need for further functional studies to understand the molecular mechanisms that control ncRNA-directed developmental processes.

DNA-dependent RNA polymerases in plants

DNA-dependent RNA polymerases (Pols) transfer the genetic information stored in genomic DNA to RNA in all organisms. In eukaryotes, the typical products of nuclear Pol I, Pol II, and Pol III are ribosomal RNAs, mRNAs, and transfer RNAs, respectively. Intriguingly, plants possess two additional Pols, Pol IV and Pol V, which produce small RNAs and long noncoding RNAs, respectively, mainly for silencing transposable elements. The five plant Pols share some subunits, but their distinct functions stem from unique subunits that interact with specific regulatory factors in their transcription cycles. Here, we summarize recent advances in our understanding of plant nucleus-localized Pols, including their evolution, function, structures, and transcription cycles.

RETINOBLASTOMA-RELATED interactions with key factors of the RNA-directed DNA methylation (RdDM) pathway and its influence on root development

MAIN CONCLUSION

Our study presents evidence for a novel mechanism for RBR function in transcriptional gene silencing by interacting with key players of the RdDM pathway in Arabidopsis and several plant clades. Transposable elements and other repetitive elements are silenced by the RNA-directed DNA methylation pathway (RdDM). In RdDM, POLIV-derived transcripts are converted into double-stranded RNA (dsRNA) by the activity of RDR2 and subsequently processed into 24 nucleotide short interfering RNAs (24-nt siRNAs) by DCL3. 24-nt siRNAs serve as guides to direct AGO4-siRNA complexes to chromatin-bound POLV-derived transcripts generated from the template/target DNA. The interaction between POLV, AGO4, DMS3, DRD1, RDM1 and DRM2 promotes DRM2-mediated de novo DNA methylation. The Arabidopsis Retinoblastoma protein homolog (RBR) is a master regulator of the cell cycle, stem cell maintenance, and development. We in silico predicted and explored experimentally the protein-protein interactions (PPIs) between RBR and members of the RdDM pathway. We found that the largest subunits of POLIV and POLV (NRPD1 and NRPE1), the shared second largest subunit of POLIV and POLV (NRPD/E2), RDR1, RDR2, DCL3, DRM2, and SUVR2 contain canonical and non-canonical RBR binding motifs and several of them are conserved since algae and bryophytes. We validated experimentally PPIs between Arabidopsis RBR and several of the RdDM pathway proteins. Moreover, seedlings from loss-of-function mutants in RdDM and RBR show similar phenotypes in the root apical meristem. We show that RdDM and SUVR2 targets are up-regulated in the 35S:AmiGO-RBR background.

Methylation in DNA, histone, and RNA during flowering under stress condition: A review

Flowering is the most critical transition period in the whole lifecycle of plants, and it is a highly sensitive period to stress. New combinations of temperature, drought stress, carbon dioxide and other abiotic/biotic conditions resulting from contemporary climate change affect the flowering process. Plants have evolved several strategies to deal with environmental stresses, including epigenetic modifications. Numerous studies show that environmental stresses trigger methylation/demethylation during flowering to preserve/accelerate plant lifecycle. What's more, histone and DNA methylation can be induced to respond to stresses, resulting in changes of flowering gene expression and enhancing stress tolerance in plants. Furthermore, RNA methylation may influence stress-regulated flowering by regulating mRNA stability and antioxidant mechanism. Our review presents the involvement of methylation in stress-repressed and stress-induced flowering. The crosstalk between methylation and small RNAs, phytohormones and exogenous substances (such as salicylic acid, nitric oxide) during flowering under different stresses were discussed. The latest regulatory evidence of RNA methylation in stress-regulated flowering was collected for the first time. Meanwhile, the limited evidences of methylation in biotic stress-induced flowering were summarized. Thus, the review provides insights into understanding of methylation mechanism in stress-regulated flowering and makes use for the development of regulating plant flowering at epigenetic level in the future.

Expansion and contraction of small RNA and methylation machinery throughout plant evolution

The revolution in sequencing has created a wealth of plant genomes that can be mined to understand the evolution of biological complexity. Complexity is often driven by gene duplication, which allows paralogs to specialize in an activity of the ancestral gene or acquire novel functions. Angiosperms encode a variety of gene silencing pathways that share related machinery for small RNA biosynthesis and function. Recent phylogenetic analysis of these gene families plots the expansion, specialization, and occasional contraction of this core machinery. This analysis reveals the ancient origin of RNA-directed DNA Methylation in early land plants, or possibly their algal ancestors, as well as ongoing duplications that evolve novel small RNA pathways.

The renaissance and enlightenment of Marchantia as a model system

The liverwort Marchantia polymorpha has been utilized as a model for biological studies since the 18th century. In the past few decades, there has been a Renaissance in its utilization in genomic and genetic approaches to investigating physiological, developmental, and evolutionary aspects of land plant biology. The reasons for its adoption are similar to those of other genetic models, e.g. simple cultivation, ready access via its worldwide distribution, ease of crossing, facile genetics, and more recently, efficient transformation, genome editing, and genomic resources. The haploid gametophyte dominant life cycle of M. polymorpha is conducive to forward genetic approaches. The lack of ancient whole-genome duplications within liverworts facilitates reverse genetic approaches, and possibly related to this genomic stability, liverworts possess sex chromosomes that evolved in the ancestral liverwort. As a representative of one of the three bryophyte lineages, its phylogenetic position allows comparative approaches to provide insights into ancestral land plants. Given the karyotype and genome stability within liverworts, the resources developed for M. polymorpha have facilitated the development of related species as models for biological processes lacking in M. polymorpha.

A null allele of the pol IV second subunit impacts stature and reproductive development in Oryza sativa

All eukaryotes possess three DNA-dependent RNA polymerases, Pols I-III, while land plants possess two additional polymerases, Pol IV and Pol V. Derived through duplication of Pol II subunits, Pol IV produces 24-nt short interfering RNAs that interact with Pol V transcripts to target de novo DNA methylation and silence transcription of transposons. Members of the grass family encode additional duplicated subunits of Pol IV and V, raising questions regarding the function of each paralog. In this study, we identify a null allele of the putative Pol IV second subunit, NRPD2, and demonstrate that NRPD2 is the sole subunit functioning with NRPD1 in small RNA production and CHH methylation in leaves. Homozygous nrpd2 mutants have neither gametophytic defects nor embryo lethality, although adult plants are dwarf and sterile.

The Mediator Complex: A Central Coordinator of Plant Adaptive Responses to Environmental Stresses

As sessile organisms, plants are constantly exposed to a variety of environmental stresses and have evolved adaptive mechanisms, including transcriptional reprogramming, in order to survive or acclimate under adverse conditions. Over the past several decades, a large number of gene-specific transcription factors have been identified in the transcriptional regulation of plant adaptive responses. The Mediator complex plays a key role in transducing signals from gene-specific transcription factors to the transcription machinery to activate or repress target gene expression. Since its first purification about 15 years ago, plant Mediator complex has been extensively analyzed for its composition and biological functions. Mutants of many plant Mediator subunits are not lethal but are compromised in growth, development and response to biotic and abiotic stress, underscoring a particularly important role in plant adaptive responses. Plant Mediator subunits also interact with partners other than transcription factors and components of the transcription machinery, indicating the complexity of the regulation of gene expression by plant Mediator complex. Here, we present a comprehensive discussion of recent analyses of the structure and function of plant Mediator complex, with a particular focus on its roles in plant adaptive responses to a wide spectrum of environmental stresses and associated biological processes.

Gene duplication and rate variation in the evolution of plastid ACCase and Clp genes in angiosperms

While the chloroplast (plastid) is known for its role in photosynthesis, it is also involved in many other metabolic pathways essential for plant survival. As such, plastids contain an extensive suite of enzymes required for non-photosynthetic processes. The evolution of the associated genes has been especially dynamic in flowering plants (angiosperms), including examples of gene duplication and extensive rate variation. We examined the role of ongoing gene duplication in two key plastid enzymes, the acetyl-CoA carboxylase (ACCase) and the caseinolytic protease (Clp), responsible for fatty acid biosynthesis and protein turnover, respectively. In plants, there are two ACCase complexes-a homomeric version present in the cytosol and a heteromeric version present in the plastid. Duplications of the nuclear-encoded homomeric ACCase gene and retargeting of one resultant protein to the plastid have been previously reported in multiple species. We find that these retargeted homomeric ACCase proteins exhibit elevated rates of sequence evolution, consistent with neofunctionalization and/or relaxation of selection. The plastid Clp complex catalytic core is composed of nine paralogous proteins that arose via ancient gene duplication in the cyanobacterial/plastid lineage. We show that further gene duplication occurred more recently in the nuclear-encoded core subunits of this complex, yielding additional paralogs in many species of angiosperms. Moreover, in six of eight cases, subunits that have undergone recent duplication display increased rates of sequence evolution relative to those that have remained single copy. We also compared substitution patterns between pairs of Clp core paralogs to gain insight into post-duplication evolutionary routes. These results show that gene duplication and rate variation continue to shape the plastid proteome.

Several Isoforms for Each Subunit Shared by RNA Polymerases are Differentially Expressed in the Cultivated Olive Tree (Olea europaea L.)

Plants contain five nuclear RNA polymerases, with RNA pols IV and V in addition to conserved eukaryotic RNA pols I, II, and III. These transcriptional complexes share five common subunits, which have been extensively analyzed only in yeasts. By taking advantage of the recently published olive tree cultivar (Olea europaea L. cv. Picual) genome, we performed a genome-wide analysis of the genomic composition corresponding to subunits common to RNA pols. The cultivated olive tree genome is quite complex and contains many genes with several copies. We also investigated, for the first time, gene expression patterns for subunits common to RNA pols using RNA-Seq under different economically and biologically relevant conditions for the cultivar "Picual": tissues/organs, biotic and abiotic stresses, and early development from seeds. Our results demonstrated the existence of a multigene family of subunits common to RNA pols, and a variable number of paralogs for each subunit in the olive cultivar "Picual." Furthermore, these isoforms display specific and differentiated expression profiles depending on the isoform and growth conditions, which may be relevant for their role in olive tree biology.

Evolution of CG Methylation Maintenance Machinery in Plants

Cytosine methylation is an epigenetic mark present in most eukaryotic genomes that contributes to the regulation of gene expression and the maintenance of genome stability. DNA methylation mostly occurs at CG sequences, where it is initially deposited by de novo DNA methyltransferases and propagated by maintenance DNA methyltransferases (DNMT) during DNA replication. In this review, we first summarize the mechanisms maintaining CG methylation in mammals that involve the DNA Methyltransferase 1 (DNMT1) enzyme and its cofactor, UHRF1 (Ubiquitin-like with PHD and RING Finger domain 1). We then discuss the evolutionary conservation and diversification of these two core factors in the plant kingdom and speculate on potential functions of novel homologues typically observed in land plants but not in mammals.

The visual story of data storage: From storage properties to user interfaces

About fifty times more data has been created than there are stars in the observable universe. Current trends in data creation and consumption mean that the devices and storage media we use will require more physical space. Novel data storage media such as DNA are considered a viable alternative. Yet, the introduction of new storage technologies should be accompanied by an evaluation of user requirements. To assess such needs, we designed and conducted a survey to rank different storage properties adapted for visualization. That is, accessibility, capacity, usage, mutability, lifespan, addressability, and typology. Withal, we reported different storage devices over time while ranking them by their properties. Our results indicated a timeline of three distinct periods: magnetic, optical and electronic, and alternative media. Moreover, by investigating user interfaces across different operating systems, we observed a predominant presence of bar charts and tree maps for the usage of a medium and its file directory hierarchy, respectively. Taken together with the results of our survey, this allowed us to create a customized user interface that includes data visualizations that can be toggled for both user groups: Experts and Public.

Atypical DNA methylation, sRNA-size distribution, and female gametogenesis in Utricularia gibba

The most studied DNA methylation pathway in plants is the RNA Directed DNA Methylation (RdDM), a conserved mechanism that involves the role of noncoding RNAs to control the expansion of the noncoding genome. Genome-wide DNA methylation levels have been reported to correlate with genome size. However, little is known about the catalog of noncoding RNAs and the impact on DNA methylation in small plant genomes with reduced noncoding regions. Because of the small length of intergenic regions in the compact genome of the carnivorous plant Utricularia gibba, we investigated its repertoire of noncoding RNA and DNA methylation landscape. Here, we report that, compared to other angiosperms, U. gibba has an unusual distribution of small RNAs and reduced global DNA methylation levels. DNA methylation was determined using a novel strategy based on long-read DNA sequencing with the Pacific Bioscience platform and confirmed by whole-genome bisulfite sequencing. Moreover, some key genes involved in the RdDM pathway may not represented by compensatory paralogs or comprise truncated proteins, for example, U. gibba DICER-LIKE 3 (DCL3), encoding a DICER endonuclease that produces 24-nt small-interfering RNAs, has lost key domains required for complete function. Our results unveil that a truncated DCL3 correlates with a decreased proportion of 24-nt small-interfering RNAs, low DNA methylation levels, and developmental abnormalities during female gametogenesis in U. gibba. Alterations in female gametogenesis are reminiscent of RdDM mutant phenotypes in Arabidopsis thaliana. It would be interesting to further study the biological implications of the DCL3 truncation in U. gibba, as it could represent an initial step in the evolution of RdDM pathway in compact genomes.

Intertwined evolution of plant epigenomes and genomes

DNA methylation is found across eukaryotes; however, plants have evolved patterns and pathways of DNA methylation that are distinct from animals and fungi. DNA methylation shapes the evolution of genomes through its direct roles in transposon silencing, gene expression, genome stability, and its impact on mutation rates. In return the diversity of DNA methylation across species is shaped by genome sequence evolution. Extensive diversification of key DNA methylation pathways has continued in plants through gene duplication and loss. Meanwhile, frequent movement of transposons has altered local DNA methylation patterns and the genes affected. Only recently has the diversity and evolutionary history of plant DNA methylation become evident with the availability of increasing genomic and epigenomic data. However, much remains unresolved regarding the evolutionary forces that have shaped the dynamics of the complex and intertwined history of plant genome and epigenome evolution.

Eukaryotic RNA Polymerases: The Many Ways to Transcribe a Gene

In eukaryotic cells, three nuclear RNA polymerases (RNA pols) carry out the transcription from DNA to RNA, and they all seem to have evolved from a single enzyme present in the common ancestor with archaea. The multiplicity of eukaryotic RNA pols allows each one to remain specialized in the synthesis of a subset of transcripts, which are different in the function, length, cell abundance, diversity, and promoter organization of the corresponding genes. We hypothesize that this specialization of RNA pols has conditioned the evolution of the regulatory mechanisms used to transcribe each gene subset to cope with environmental changes. We herein present the example of the homeostatic regulation of transcript levels versus changes in cell volume. We propose that the diversity and instability of messenger RNAs, transcribed by RNA polymerase II, have conditioned the appearance of regulatory mechanisms based on different gene promoter strength and mRNA stability. However, for the regulation of ribosomal RNA levels, which are very stable and transcribed mainly by RNA polymerase I from only one promoter, different mechanisms act based on gene copy variation, and a much simpler regulation of the synthesis rate.

Genome-wide analyses across Viridiplantae reveal the origin and diversification of small RNA pathway-related genes

Small RNAs play a major role in the post-transcriptional regulation of gene expression in eukaryotes. Despite the evolutionary importance of streptophyte algae, knowledge on small RNAs in this group of green algae is almost non-existent. We used genome and transcriptome data of 34 algal and plant species, and performed genome-wide analyses of small RNA (miRNA & siRNA) biosynthetic and degradation pathways. The results suggest that Viridiplantae started to evolve plant-like miRNA biogenesis and degradation after the divergence of the Mesostigmatophyceae in the streptophyte algae. We identified two major evolutionary transitions in small RNA metabolism in streptophyte algae; during the first transition, the origin of DCL-New, DCL1, AGO1/5/10 and AGO4/6/9 in the last common ancestor of Klebsormidiophyceae and all other streptophytes could be linked to abiotic stress responses and evolution of multicellularity in streptophytes. During the second transition, the evolution of DCL 2,3,4, and AGO 2,3,7 as well as DRB1 in the last common ancestor of Zygnematophyceae and embryophytes, suggests their possible contribution to pathogen defense and antibacterial immunity. Overall, the origin and diversification of DICER and AGO along with several other small RNA pathway-related genes among streptophyte algae suggested progressive adaptations of streptophyte algae during evolution to a subaerial environment.

Wang, Liang et al. conduct a genome-wide investigation into the origin of small RNA pathway-related genes in Viridiplantae. Their findings suggest that streptophyte algae progressively adapted to a subaerial environment through generation of these pathways.

Non-coding RNA polymerases that silence transposable elements and reprogram gene expression in plants

Multisubunit RNA polymerase (Pol) complexes are the core machinery for gene expression in eukaryotes. The enzymes Pol I, Pol II and Pol III transcribe distinct subsets of nuclear genes. This family of nuclear RNA polymerases expanded in terrestrial plants by the duplication of Pol II subunit genes. Two Pol II-related enzymes, Pol IV and Pol V, are highly specialized in the production of regulatory, non-coding RNAs. Pol IV and Pol V are the central players of RNA-directed DNA methylation (RdDM), an RNA interference pathway that represses transposable elements (TEs) and selected genes. Genetic and biochemical analyses of Pol IV/V subunits are now revealing how these enzymes evolved from ancestral Pol II to sustain non-coding RNA biogenesis in silent chromatin. Intriguingly, Pol IV-RdDM regulates genes that influence flowering time, reproductive development, stress responses and plant–pathogen interactions. Pol IV target genes vary among closely related taxa, indicating that these regulatory circuits are often species-specific. Data from crops like maize, rice, tomato and Brassicarapa suggest that dynamic repositioning of TEs, accompanied by Pol IV targeting to TE-proximal genes, leads to the reprogramming of plant gene expression over short evolutionary timescales.

RNA-directed DNA Methylation

RNA-directed DNA methylation (RdDM) is a biological process in which non-coding RNA molecules direct the addition of DNA methylation to specific DNA sequences. The RdDM pathway is unique to plants, although other mechanisms of RNA-directed chromatin modification have also been described in fungi and animals. To date, the RdDM pathway is best characterized within angiosperms (flowering plants), and particularly within the model plant Arabidopsis thaliana. However, conserved RdDM pathway components and associated small RNAs (sRNAs) have also been found in other groups of plants, such as gymnosperms and ferns. The RdDM pathway closely resembles other sRNA pathways, particularly the highly conserved RNAi pathway found in fungi, plants, and animals. Both the RdDM and RNAi pathways produce sRNAs and involve conserved Argonaute, Dicer and RNA-dependent RNA polymerase proteins.

RdDM has been implicated in a number of regulatory processes in plants. The DNA methylation added by RdDM is generally associated with transcriptional repression of the genetic sequences targeted by the pathway. Since DNA methylation patterns in plants are heritable, these changes can often be stably transmitted to progeny. As a result, one prominent role of RdDM is the stable, transgenerational suppression of transposable element (TE) activity. RdDM has also been linked to pathogen defense, abiotic stress responses, and the regulation of several key developmental transitions. Although the RdDM pathway has a number of important functions, RdDM-defective mutants in Arabidopsis thaliana are viable and can reproduce, which has enabled detailed genetic studies of the pathway. However, RdDM mutants can have a range of defects in different plant species, including lethality, altered reproductive phenotypes, TE upregulation and genome instability, and increased pathogen sensitivity. Overall, RdDM is an important pathway in plants that regulates a number of processes by establishing and reinforcing specific DNA methylation patterns, which can lead to transgenerational epigenetic effects on gene expression and phenotype.

Formation and Recognition of UV-Induced DNA Damage within Genome Complexity

Ultraviolet (UV) light is a natural genotoxic agent leading to the formation of photolesions endangering the genomic integrity and thereby the survival of living organisms. To prevent the mutagenetic effect of UV, several specific DNA repair mechanisms are mobilized to accurately maintain genome integrity at photodamaged sites within the complexity of genome structures. However, a fundamental gap remains to be filled in the identification and characterization of factors at the nexus of UV-induced DNA damage, DNA repair, and epigenetics. This review brings together the impact of the epigenomic context on the susceptibility of genomic regions to form photodamage and focuses on the mechanisms of photolesions recognition through the different DNA repair pathways.

Small RNA Function in Plants: From Chromatin to the Next Generation

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."

RNA-directed DNA Methylation and sexual reproduction: expanding beyond the seed

Two trends are changing our understanding of RNA-directed DNA methylation. In model systems like Arabidopsis, tissue-specific analysis of DNA methylation is uncovering dynamic changes in methylation during sexual reproduction and unraveling the contribution of maternal and paternal epigenomes to the developing embryo. These studies indicate that RNA-directed DNA Methylation might be important for mediating balance between maternal and paternal contributions to the endosperm. At the same time, researchers are moving beyond Arabidopsis to illuminate the ancestral role of RdDM in non-flowering plants that lack an endosperm, suggesting that RdDM might play a broader role in sexual reproduction.

Arabidopsis RNA Polymerase IV generates 21-22 nucleotide small RNAs that can participate in RNA-directed DNA methylation and may regulate genes

The plant-specific RNA Polymerase IV (Pol IV) transcribes heterochromatic regions, including many transposable elements (TEs), with the well-described role of generating 24 nucleotide (nt) small interfering RNAs (siRNAs). These siRNAs target DNA methylation back to TEs to reinforce the boundary between heterochromatin and euchromatin. In the male gametophytic phase of the plant life cycle, pollen, Pol IV switches to generating primarily 21-22 nt siRNAs, but the biogenesis and function of these siRNAs have been enigmatic. In contrast to being pollen-specific, we identified that Pol IV generates these 21-22 nt siRNAs in sporophytic tissues, likely from the same transcripts that are processed into the more abundant 24 nt siRNAs. The 21-22 nt forms are specifically generated by the combined activities of DICER proteins DCL2/DCL4 and can participate in RNA-directed DNA methylation. These 21-22 nt siRNAs are also loaded into ARGONAUTE1 (AGO1), which is known to function in post-transcriptional gene regulation. Like other plant siRNAs and microRNAs incorporated into AGO1, we find a signature of genic mRNA cleavage at the predicted target site of these siRNAs, suggesting that Pol IV-generated 21-22 nt siRNAs may function to regulate gene transcript abundance. Our data provide support for the existing model that in pollen Pol IV functions in gene regulation. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Construction and characterization of a knock-down RNA interference line of OsNRPD1 in rice (Oryza sativa ssp japonica cv Nipponbare)

In plants, RNA-directed DNA methylation (RdDM) is a silencing mechanism relying on the production of 24-nt small interfering RNAs (siRNAs) by RNA POLYMERASE IV (Pol IV) to trigger methylation and inactivation of transposable elements (TEs). We present the construction and characterization of osnrpd1, a knock-down RNA interference line of OsNRPD1 gene that encodes the largest subunit of Pol IV in rice (Oryza sativa ssp japonica cv Nipponbare). We show that osnrpd1 displays a lower accumulation of OsNRPD1 transcripts, associated with an overall reduction of 24-nt siRNAs and DNA methylation level in all three contexts, CG, CHG and CHH. We uncovered new insertions of known active TEs, the LTR retrotransposons Tos17 and Lullaby and the long interspersed nuclear element-type retrotransposon Karma. However, we did not observe any clear developmental phenotype, contrary to what was expected for a mutant severely affected in RdDM. In addition, despite the presence of many putatively functional TEs in the rice genome, we found no evidence of in planta global reactivation of transposition. This knock-down of OsNRPD1 likely led to a weakly affected line, with no effect on development and a limited effect on transposition. We discuss the possibility that a knock-out mutation of OsNRPD1 would cause sterility in rice. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

The NRPD1 N-terminus contains a Pol IV-specific motif that is critical for genome surveillance in Arabidopsis

RNA-guided surveillance systems constrain the activity of transposable elements (TEs) in host genomes. In plants, RNA polymerase IV (Pol IV) transcribes TEs into primary transcripts from which RDR2 synthesizes double-stranded RNA precursors for small interfering RNAs (siRNAs) that guide TE methylation and silencing. How the core subunits of Pol IV, homologs of RNA polymerase II subunits, diverged to support siRNA biogenesis in a TE-rich, repressive chromatin context is not well understood. Here we studied the N-terminus of Pol IV’s largest subunit, NRPD1. Arabidopsis lines harboring missense mutations in this N-terminus produce wild-type (WT) levels of NRPD1, which co-purifies with other Pol IV subunits and RDR2. Our in vitro transcription and genomic analyses reveal that the NRPD1 N-terminus is critical for robust Pol IV-dependent transcription, siRNA production and DNA methylation. However, residual RNA-directed DNA methylation observed in one mutant genotype indicates that Pol IV can operate uncoupled from the high siRNA levels typically observed in WT plants. This mutation disrupts a motif uniquely conserved in Pol IV, crippling the enzyme's ability to inhibit retrotransposon mobilization. We propose that the NRPD1 N-terminus motif evolved to regulate Pol IV function in genome surveillance.

Establishment, maintenance, and biological roles of non-CG methylation in plants

Cytosine DNA methylation is prevalent throughout eukaryotes and prokaryotes. While most commonly thought of as being localized to dinucleotide CpG sites, non-CG sites can also be modified. Such non-CG methylation is widespread in plants, occurring at trinucleotide CHG and CHH (H = A, T, or C) sequence contexts. The prevalence of non-CG methylation in plants is due to the plant-specific CHROMOMETHYLASE (CMT) and RNA-directed DNA Methylation (RdDM) pathways. These pathways have evolved through multiple rounds of gene duplication and gene loss, generating epigenomic variation both within and between species. They regulate both transposable elements and genes, ensure genome integrity, and ultimately influence development and environmental responses. In these capacities, non-CG methylation influence and shape plant genomes.

Tissue-specific transposon-associated small RNAs in the gymnosperm tree, Norway spruce

BACKGROUND

Small RNAs (sRNAs) are regulatory molecules impacting on gene expression and transposon activity. MicroRNAs (miRNAs) are responsible for tissue-specific and environmentally-induced gene repression. Short interfering RNAs (siRNA) are constitutively involved in transposon silencing across different type of tissues. The male gametophyte in angiosperms has a unique set of sRNAs compared to vegetative tissues, including phased siRNAs from intergenic or genic regions, or epigenetically activated siRNAs. This is contrasted by a lack of knowledge about the sRNA profile of the male gametophyte of gymnosperms.

RESULTS

Here, we isolated mature pollen from male cones of Norway spruce and investigated its sRNA profiles. While 21-nt sRNAs is the major size class of sRNAs in needles, in pollen 21-nt and 24-nt sRNAs are the most abundant size classes. Although the 24-nt sRNAs were exclusively derived from TEs in pollen, both 21-nt and 24-nt sRNAs were associated with TEs. We also investigated sRNAs from somatic embryonic callus, which has been reported to contain 24-nt sRNAs. Our data show that the 24-nt sRNA profiles are tissue-specific and differ between pollen and cell culture.

CONCLUSION

Our data reveal that gymnosperm pollen, like angiosperm pollen, has a unique sRNA profile, differing from vegetative leaf tissue. Thus, our results reveal that angiosperm and gymnosperm pollen produce new size classes not present in vegetative tissues; while in angiosperm pollen 21-nt sRNAs are generated, in the gymnosperm Norway spruce 24-nt sRNAs are generated. The tissue-specific production of distinct TE-derived sRNAs in angiosperms and gymnosperms provides insights into the diversification process of sRNAs in TE silencing pathways between the two groups of seed plants.

Tracing the origin and evolution history of methylation-related genes in plants

BACKGROUND

DNA methylation is a crucial epigenetic modification, which is involved in many biological processes, including gene expression regulation, embryonic development, cell differentiation and genomic imprinting etc. And it also involves many key regulatory genes in eukaryotes. By tracing the evolutionary history of methylation-related genes, we can understand the origin and expansion time of these genes, which helps to understand the evolutionary history of plants, and we can also understand the changes of DNA methylation patterns in different species. However, most studies on the evolution of methylation-related genes failed to be carried out for the whole DNA methylation pathway.

RESULTS

In this study, we conducted a comprehensive identification of 33 methylation-related genes in 77 species, and investigated gene origin and evolution throughout the plant kingdom. We found that the origin of genes responsible for methylation maintenance and demethylation evolved early, while most de novo methylation-related genes appeared late. The methylation-related genes were expanded by whole genome duplication and tandem replication, but were also accompanied by a large number of gene absence events in different species. The gene length and intron length varied a lot in different species, but exon structure and functional domains were relatively conserved. The phylogenetic relationships of methylation-related genes were traced to reveal the evolution history of DNA methylation in different species. The expression patterns of methylation-related genes have changed during the evolution of species, and the expression patterns of these genes in different species can be clustered into four categories.

CONCLUSIONS

The study describes a global characterization of DNA methylation-related genes in the plant kingdom. The similarities and differences in origin time, gene structure and phylogenetic relationship of these genes lead us to understand the evolutionary conservation and dynamics of DNA methylation in plants.

DNA methylation in Marchantia polymorpha

Methylation of DNA is an epigenetic mechanism for the control of gene expression. Alterations in the regulatory pathways involved in the establishment, perpetuation and removal of DNA methylation can lead to severe developmental alterations. Our understanding of the mechanistic aspects and relevance of DNA methylation comes from remarkable studies in well-established angiosperm plant models including maize and Arabidopsis. The study of plant models positioned at basal lineages opens exciting opportunities to expand our knowledge on the function and evolution of the components of DNA methylation. In this Tansley Insight, we summarize current progress in our understanding of the molecular basis and relevance of DNA methylation in the liverwort Marchantia polymorpha.

An Insect Counteradaptation against Host Plant Defenses Evolved through Concerted Neofunctionalization

Antagonistic chemical interactions between herbivorous insects and their host plants are often thought to coevolve in a stepwise process, with an evolutionary innovation on one side being countered by a corresponding advance on the other. Glucosinolate sulfatase (GSS) enzyme activity is essential for the Diamondback moth, Plutella xylostella, to overcome a highly diversified secondary metabolite-based host defense system in the Brassicales. GSS genes are located in an ancient cluster of arylsulfataselike genes, but the exact roles of gene copies and their evolutionary trajectories are unknown. Here, we combine a functional investigation of duplicated insect arylsulfatases with an analysis of associated nucleotide substitution patterns. We show that the Diamondback moth genome encodes three GSSs with distinct substrate spectra and distinct expression patterns in response to glucosinolates. Contrary to our expectations, early functional diversification of gene copies was not indicative of a coevolutionary arms race between host and herbivore. Instead, both copies of a duplicated arylsulfatase gene evolved concertedly in the context of an insect host shift to acquire novel detoxifying functions under positive selection, a pattern of duplicate gene retention that we call “concerted neofunctionalization.”

24-nt reproductive phasiRNAs are broadly present in angiosperms

Small RNAs are key regulators in plant growth and development. One subclass, phased siRNAs (phasiRNAs) require a trigger microRNA for their biogenesis. In grasses, two pathways yield abundant phasiRNAs during anther development; miR2275 triggers one class, 24-nt phasiRNAs, coincident with meiosis, while a second class of 21-nt phasiRNAs are present in premeiotic anthers. Here we report that the 24-nt phasiRNA pathway is widely present in flowering plants, indicating that 24-nt reproductive phasiRNAs likely originated with the evolutionary emergence of anthers. Deep comparative genomic analyses demonstrated that this miR2275/24-nt phasiRNA pathway is widely present in eudicots plants, however, it is absent in legumes and in the model plant Arabidopsis, demonstrating a dynamic evolutionary history of this pathway. In Solanaceae species, 24-nt phasiRNAs were observed, but the miR2275 trigger is missing and some loci displaying 12-nt phasing. Both the miR2275-triggered and Solanaceae 24-nt phasiRNAs are enriched in meiotic stages, implicating these phasiRNAs in anther and/or pollen development, a spatiotemporal pattern consistent in all angiosperm lineages that deploy them.

Four putative SWI2/SNF2 chromatin remodelers have dual roles in regulating DNA methylation in Arabidopsis

DNA methylation is a conserved epigenetic mark that is critical for many biological processes in plants and mammals. In Arabidopsis, the antagonistic activities of RNA-directed DNA methylation (RdDM) and ROS1-dependent active DNA demethylation are key for the dynamic regulation of locus-specific DNA methylation. However, the molecular factors that coordinate RdDM and active demethylation are largely unknown. Here we report that CLSY4 and its three paralogous SWI2/SNF2-type chromatin-remodeling proteins function in both RdDM and DNA demethylation in Arabidopsis. We initially identified CLSY4 in a genetic screen for DNA demethylation factors and subsequently demonstrated that it also is important in RdDM. Comprehensive genetic analyses using single and high order mutants of CLSY family proteins revealed their roles as double agents in the balance between methylation and demethylation reactions. The four CLSY proteins collectively are necessary for the canonical RdDM pathway; at the same time, each CLSY likely mediates DNA demethylation at specific loci where DNA methylation depends on RdDM. These results indicate that the four chromatin-remodeling proteins have dual functions in regulating genomic DNA methylation, and thus provide new insights into the dynamic regulation of DNA methylation in a model multicellular eukaryotic organism.

Modeling Interactions between Transposable Elements and the Plant Epigenetic Response: A Surprising Reliance on Element Retention

Transposable elements (TEs) compose the majority of angiosperm DNA. Plants counteract TE activity by silencing them epigenetically. One form of epigenetic silencing requires 21-22 nt small interfering RNAs that act to degrade TE mRNA and may also trigger DNA methylation. DNA methylation is reinforced by a second mechanism, the RNA-dependent DNA methylation (RdDM) pathway. RdDM relies on 24 nt small interfering RNAs and ultimately establishes TEs in a quiescent state. These host factors interact at a systems level, but there have been no system level analyses of their interactions. Here, we define a deterministic model that represents the propagation of active TEs, aspects of the host response and the accumulation of silenced TEs. We describe general properties of the model and also fit it to biological data in order to explore two questions. The first is why two overlapping pathways are maintained, given that both are likely energetically expensive. Under our model, RdDM silenced TEs effectively even when the initiation of silencing was weak. This relationship implies that only a small amount of RNAi is needed to initiate TE silencing, but reinforcement by RdDM is necessary to efficiently counter TE propagation. Second, we investigated the reliance of the host response on rates of TE deletion. The model predicted that low levels of deletion lead to few active TEs, suggesting that silencing is most efficient when methylated TEs are retained in the genome, thereby providing one explanation for the large size of plant genomes.

Evidence for a Unique DNA-Dependent RNA Polymerase in Cereal Crops

Gene duplication is an important driver for the evolution of new genes and protein functions. Duplication of DNA-dependent RNA polymerase (Pol) II subunits within plants led to the emergence of RNA Pol IV and V complexes, each of which possess unique functions necessary for RNA-directed DNA Methylation. Comprehensive identification of Pol V subunit orthologs across the monocot radiation revealed a duplication of the largest two subunits within the grasses (Poaceae), including critical cereal crops. These paralogous Pol subunits display sequence conservation within catalytic domains, but their carboxy terminal domains differ in length and character of the Ago-binding platform, suggesting unique functional interactions. Phylogenetic analysis of the catalytic region indicates positive selection on one paralog following duplication, consistent with retention via neofunctionalization. Positive selection on residue pairs that are predicted to interact between subunits suggests that paralogous subunits have evolved specific assembly partners. Additional Pol subunits as well as Pol-interacting proteins also possess grass-specific paralogs, supporting the hypothesis that a novel Pol complex with distinct function has evolved in the grass family, Poaceae.

Metagenomic sequencing suggests a diversity of RNA interference-like responses to viruses across multicellular eukaryotes

RNA interference (RNAi)-related pathways target viruses and transposable element (TE) transcripts in plants, fungi, and ecdysozoans (nematodes and arthropods), giving protection against infection and transmission. In each case, this produces abundant TE and virus-derived 20-30nt small RNAs, which provide a characteristic signature of RNAi-mediated defence. The broad phylogenetic distribution of the Argonaute and Dicer-family genes that mediate these pathways suggests that defensive RNAi is ancient, and probably shared by most animal (metazoan) phyla. Indeed, while vertebrates had been thought an exception, it has recently been argued that mammals also possess an antiviral RNAi pathway, although its immunological relevance is currently uncertain and the viral small RNAs (viRNAs) are not easily detectable. Here we use a metagenomic approach to test for the presence of viRNAs in five species from divergent animal phyla (Porifera, Cnidaria, Echinodermata, Mollusca, and Annelida), and in a brown alga-which represents an independent origin of multicellularity from plants, fungi, and animals. We use metagenomic RNA sequencing to identify around 80 virus-like contigs in these lineages, and small RNA sequencing to identify viRNAs derived from those viruses. We identified 21U small RNAs derived from an RNA virus in the brown alga, reminiscent of plant and fungal viRNAs, despite the deep divergence between these lineages. However, contrary to our expectations, we were unable to identify canonical (i.e. Drosophila- or nematode-like) viRNAs in any of the animals, despite the widespread presence of abundant micro-RNAs, and somatic transposon-derived piwi-interacting RNAs. We did identify a distinctive group of small RNAs derived from RNA viruses in the mollusc. However, unlike ecdysozoan viRNAs, these had a piRNA-like length distribution but lacked key signatures of piRNA biogenesis. We also identified primary piRNAs derived from putatively endogenous copies of DNA viruses in the cnidarian and the echinoderm, and an endogenous RNA virus in the mollusc. The absence of canonical virus-derived small RNAs from our samples may suggest that the majority of animal phyla lack an antiviral RNAi response. Alternatively, these phyla could possess an antiviral RNAi response resembling that reported for vertebrates, with cryptic viRNAs not detectable through simple metagenomic sequencing of wild-type individuals. In either case, our findings show that the antiviral RNAi responses of arthropods and nematodes, which are highly divergent from each other and from that of plants and fungi, are also highly diverged from the most likely ancestral metazoan state.

Functional Dissection of the Pol V Largest Subunit CTD in RNA-Directed DNA Methylation

Plant multisubunit RNA polymerase V (Pol V) transcription recruits Argonaute-small interfering RNA (siRNA) complexes that specify sites of RNA-directed DNA methylation (RdDM) for gene silencing. Pol V's largest subunit, NRPE1, evolved from the largest subunit of Pol II but has a distinctive C-terminal domain (CTD). We show that the Pol V CTD is dispensable for catalytic activity in vitro yet essential in vivo. One CTD subdomain (DeCL) is required for Pol V function at virtually all loci. Other CTD subdomains have locus-specific effects. In a yeast two-hybrid screen, the 3'→ 5' exoribonuclease RRP6L1 was identified as an interactor with the DeCL and glutamine-serine (QS)-rich subdomains located downstream of an Argonaute-binding subdomain. Experimental evidence indicates that RRP6L1 trims the 3' ends of Pol V transcripts sliced by Argonaute 4 (AGO4), suggesting a model whereby the CTD enables the spatial and temporal coordination of AGO4 and RRP6L1 RNA processing activities.

Identification and Evolutionary Characterization of ARGONAUTE-Binding Platforms

ARGONAUTE (AGO) proteins are eukaryotic RNA silencing effectors that interact with their binding partners via short peptide motifs known as AGO hooks. AGO hooks tend to cluster in one region of the protein to create an AGO-binding platform. In addition to the presence of AGO hooks, AGO-binding platforms are intrinsically disordered, contain tandem repeat arrays, and have weak sequence conservation even between close relatives. These characteristics make it difficult to identify and perform evolutionary analysis of these regions. Because of their weak sequence conservation, only a few AGO-binding platforms are characterized, and the evolution of these regions is only poorly understood. In this chapter we describe modules developed for computational identification and evolutionary analysis of AGO-binding platforms, with particular emphasis on understanding evolution of the tandem repeat arrays.

Revisiting Criteria for Plant MicroRNA Annotation in the Era of Big Data

MicroRNAs (miRNAs) are ∼21-nucleotide-long regulatory RNAs that arise from endonucleolytic processing of hairpin precursors. Many function as essential posttranscriptional regulators of target mRNAs and long noncoding RNAs. Alongside miRNAs, plants also produce large numbers of short interfering RNAs (siRNAs), which are distinguished from miRNAs primarily by their biogenesis (typically processed from long double-stranded RNA instead of single-stranded hairpins) and functions (typically via roles in transcriptional regulation instead of posttranscriptional regulation). Next-generation DNA sequencing methods have yielded extensive data sets of plant small RNAs, resulting in many miRNA annotations. However, it has become clear that many miRNA annotations are questionable. The sheer number of endogenous siRNAs compared with miRNAs has been a major factor in the erroneous annotation of siRNAs as miRNAs. Here, we provide updated criteria for the confident annotation of plant miRNAs, suitable for the era of "big data" from DNA sequencing. The updated criteria emphasize replication and the minimization of false positives, and they require next-generation sequencing of small RNAs. We argue that improved annotation systems are needed for miRNAs and all other classes of plant small RNAs. Finally, to illustrate the complexities of miRNA and siRNA annotation, we review the evolution and functions of miRNAs and siRNAs in plants.

Catalytic properties of RNA polymerases IV and V: accuracy, nucleotide incorporation and rNTP/dNTP discrimination

All eukaryotes have three essential nuclear multisubunit RNA polymerases, abbreviated as Pol I, Pol II and Pol III. Plants are remarkable in having two additional multisubunit RNA polymerases, Pol IV and Pol V, which synthesize noncoding RNAs that coordinate RNA-directed DNA methylation for silencing of transposons and a subset of genes. Based on their subunit compositions, Pols IV and V clearly evolved as specialized forms of Pol II, but their catalytic properties remain undefined. Here, we show that Pols IV and V differ from one another, and Pol II, in nucleotide incorporation rate, transcriptional accuracy and the ability to discriminate between ribonucleotides and deoxyribonucleotides. Pol IV transcription is considerably more error-prone than Pols II or V, which may be tolerable in its synthesis of short RNAs that serve as precursors for siRNAs targeting non-identical members of transposon families. By contrast, Pol V exhibits high fidelity transcription, similar to Pol II, suggesting a need for Pol V transcripts to faithfully reflect the DNA sequence of target loci to which siRNA–Argonaute silencing complexes are recruited.

Conservation and divergence of small RNA pathways and microRNAs in land plants

Background

As key regulators of gene expression in eukaryotes, small RNAs have been characterized in many seed plants, and pathways for their biogenesis, degradation, and action have been defined in model angiosperms. However, both small RNAs themselves and small RNA pathways are not well characterized in other land plants such as lycophytes and ferns, preventing a comprehensive evolutionary perspective on small RNAs in land plants.

Results

Using 25 representatives from major lineages of lycophytes and ferns, most of which lack sequenced genomes, we characterized small RNAs and small RNA pathways in these plants. We identified homologs of DICER-LIKE (DCL), ARGONAUTE (AGO), and other genes involved in small RNA pathways, predicted over 2600 conserved microRNA (miRNA) candidates, and performed phylogenetic analyses on small RNA pathways as well as miRNAs. Pathways underlying miRNA biogenesis, degradation, and activity were established in the common ancestor of land plants, but the 24-nucleotide siRNA pathway that guides DNA methylation is incomplete in sister species of seed plants, especially lycophytes. We show that the functional diversification of key gene families such as DCL and AGO as observed in angiosperms occurred early in land plants followed by parallel expansion of the AGO family in ferns and angiosperms. We uncovered a conserved AGO subfamily absent in angiosperms.

Conclusions

Our phylogenetic analyses of miRNAs in bryophytes, lycophytes, ferns, and angiosperms refine the time-of-origin for conserved miRNA families as well as small RNA machinery in land plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1291-2) contains supplementary material, which is available to authorized users.

Assessing the Gene Content of the Megagenome: Sugar Pine (Pinus lambertiana)

Sugar pine (Pinus lambertiana Douglas) is within the subgenus Strobus with an estimated genome size of 31 Gbp. Transcriptomic resources are of particular interest in conifers due to the challenges presented in their megagenomes for gene identification. In this study, we present the first comprehensive survey of the P. lambertiana transcriptome through deep sequencing of a variety of tissue types to generate more than 2.5 billion short reads. Third generation, long reads generated through PacBio Iso-Seq have been included for the first time in conifers to combat the challenges associated with de novo transcriptome assembly. A technology comparison is provided here to contribute to the otherwise scarce comparisons of second and third generation transcriptome sequencing approaches in plant species. In addition, the transcriptome reference was essential for gene model identification and quality assessment in the parallel project responsible for sequencing and assembly of the entire genome. In this study, the transcriptomic data were also used to address questions surrounding lineage-specific Dicer-like proteins in conifers. These proteins play a role in the control of transposable element proliferation and the related genome expansion in conifers.

The Argonaute-binding platform of NRPE1 evolves through modulation of intrinsically disordered repeats

Argonaute (Ago) proteins are important effectors in RNA silencing pathways, but they must interact with other machinery to trigger silencing. Ago hooks have emerged as a conserved motif responsible for interaction with Ago proteins, but little is known about the sequence surrounding Ago hooks that must restrict or enable interaction with specific Argonautes. Here we investigated the evolutionary dynamics of an Ago-binding platform in NRPE1, the largest subunit of RNA polymerase V. We compared NRPE1 sequences from > 50 species, including dense sampling of two plant lineages. This study demonstrates that the Ago-binding platform of NRPE1 retains Ago hooks, intrinsic disorder, and repetitive character while being highly labile at the sequence level. We reveal that loss of sequence conservation is the result of relaxed selection and frequent expansions and contractions of tandem repeat arrays. These factors allow a complete restructuring of the Ago-binding platform over 50-60 million yr. This evolutionary pattern is also detected in a second Ago-binding platform, suggesting it is a general mechanism. The presence of labile repeat arrays in all analyzed NRPE1 Ago-binding platforms indicates that selection maintains repetitive character, potentially to retain the ability to rapidly restructure the Ago-binding platform.

DNA methylome of the 20-gigabase Norway spruce genome

DNA methylation plays important roles in many biological processes, such as silencing of transposable elements, imprinting, and regulating gene expression. Many studies of DNA methylation have shown its essential roles in angiosperms (flowering plants). However, few studies have examined the roles and patterns of DNA methylation in gymnosperms. Here, we present genome-wide high coverage single-base resolution methylation maps of Norway spruce (Picea abies) from both needles and somatic embryogenesis culture cells via whole genome bisulfite sequencing. On average, DNA methylation levels of CG and CHG of Norway spruce were higher than most other plants studied. CHH methylation was found at a relatively low level; however, at least one copy of most of the RNA-directed DNA methylation pathway genes was found in Norway spruce, and CHH methylation was correlated with levels of siRNAs. In comparison with needles, somatic embryogenesis culture cells that are used for clonally propagating spruce trees showed lower levels of CG and CHG methylation but higher level of CHH methylation, suggesting that like in other species, these culture cells show abnormal methylation patterns.