MicroRNAs (miRNAs) control gene expression in animals and plants. Like another class of small RNAs, siRNAs, they affect gene expression posttranscriptionally. While siRNAs in addition act in transcriptional gene silencing, a role of miRNAs in transcriptional regulation has been less clear. We show here that in moss Physcomitrella patens mutants without a DICER-LIKE1b gene, maturation of miRNAs is normal but cleavage of target RNAs is abolished and levels of these transcripts are drastically reduced. These mutants accumulate miRNA:target-RNA duplexes and show hypermethylation of the genes encoding target RNAs, leading to gene silencing. This pathway occurs also in the wild-type upon hormone treatment. We propose that initiation of epigenetic silencing by DNA methylation depends on the ratio of the miRNA and its target RNA.

BackgroundDevelopment of eukaryotic organisms is controlled by transcription factors that trigger specific and global changes in gene expression programs. In plants, MADS-domain transcription factors act as master regulators of developmental switches and organ specification. However, the mechanisms by which these factors dynamically regulate the expression of their target genes at different developmental stages are still poorly understood.ResultsWe characterized the relationship of chromatin accessibility, gene expression, and DNA binding of two MADS-domain proteins at different stages of Arabidopsis flower development. Dynamic changes in APETALA1 and SEPALLATA3 DNA binding correlated with changes in gene expression, and many of the target genes could be associated with the developmental stage in which they are transcriptionally controlled. We also observe dynamic changes in chromatin accessibility during flower development. Remarkably, DNA binding of APETALA1 and SEPALLATA3 is largely independent of the accessibility status of their binding regions and it can precede increases in DNA accessibility. These results suggest that APETALA1 and SEPALLATA3 may modulate chromatin accessibility, thereby facilitating access of other transcriptional regulators to their target genes.ConclusionsOur findings indicate that different homeotic factors regulate partly overlapping, yet also distinctive sets of target genes in a partly stage-specific fashion. By combining the information from DNA-binding and gene expression data, we are able to propose models of stage-specific regulatory interactions, thereby addressing dynamics of regulatory networks throughout flower development. Furthermore, MADS-domain TFs may regulate gene expression by alternative strategies, one of which is modulation of chromatin accessibility.

Endogenous 24 nt short interfering RNAs (siRNAs), derived mostly from intergenic and repetitive genomic regions, constitute a major class of endogenous small RNAs in flowering plants. Accumulation of Arabidopsis thaliana 24 nt siRNAs requires the Dicer family member DCL3, and clear homologs of DCL3 exist in both flowering and non-flowering plants. However, the absence of a conspicuous 24 nt peak in the total RNA populations of several non-flowering plants has raised the question of whether this class of siRNAs might, in contrast to the ancient 21 nt microRNAs (miRNAs) and 21–22 nt trans-acting siRNAs (tasiRNAs), be an angiosperm-specific innovation. Analysis of non-miRNA, non-tasiRNA hotspots of small RNA production within the genome of the moss Physcomitrella patens revealed multiple loci that consistently produced a mixture of 21–24 nt siRNAs with a peak at 23 nt. These Pp23SR loci were significantly enriched in transposon content, depleted in overlap with annotated genes, and typified by dense concentrations of the 5-methyl cytosine (5 mC) DNA modification. Deep sequencing of small RNAs from two independent Ppdcl3 mutants showed that the P. patens DCL3 homolog is required for the accumulation of 22–24 nt siRNAs, but not 21 nt siRNAs, at Pp23SR loci. The 21 nt component of Pp23SR-derived siRNAs was also unaffected by a mutation in the RNA-dependent RNA polymerase mutant Pprdr6. Transcriptome-wide, Ppdcl3 mutants failed to accumulate 22–24 nt small RNAs from repetitive regions while transcripts from two abundant families of long terminal repeat (LTR) retrotransposon-associated reverse transcriptases were up-regulated. Ppdcl3 mutants also displayed an acceleration of leafy gametophore production, suggesting that repetitive siRNAs may play a role in the development of P. patens. We conclude that intergenic/repeat-derived siRNAs are indeed a broadly conserved, distinct class of small regulatory RNAs within land plants.

SUMMARYThe moss Physcomitrella patens is an important model organism for studying plant evolution, development, physiology and biotechnology. Here we have generated microarray gene expression data covering the principal developmental stages, culture forms and some environmental/stress conditions. Example analyses of developmental stages and growth conditions as well as abiotic stress treatments demonstrate that (i) growth stage is dominant over culture conditions, (ii) liquid culture is not stressful for the plant, (iii) low pH might aid protoplastation by reduced expression of cell wall structure genes, (iv) largely the same gene pool mediates response to dehydration and rehydration, and (v) AP2/EREBP transcription factors play important roles in stress response reactions. With regard to the AP2 gene family, phylogenetic analysis and comparison with Arabidopsis thaliana shows commonalities as well as uniquely expressed family members under drought, light perturbations and protoplastation. Gene expression profiles for P. patens are available for the scientific community via the easy-to-use tool at https://www.genevestigator.com. By providing large-scale expression profiles, the usability of this model organism is further enhanced, for example by enabling selection of control genes for quantitative real-time PCR. Now, gene expression levels across a broad range of conditions can be accessed online for P. patens.

Trans-acting small interfering RNAs (ta-siRNAs) are plant-specific siRNAs released from TAS precursor transcripts after microRNA-dependent cleavage, conversion into double-stranded RNA, and Dicer-dependent phased processing. Like microRNAs (miRNAs), ta-siRNAs direct site-specific cleavage of target RNAs at sites of extensive complementarity. Here, we show that the DICER-LIKE 4 protein of Physcomitrella patens (PpDCL4) is essential for the biogenesis of 21 nucleotide (nt) ta-siRNAs. In ΔPpDCL4 mutants, off-sized 23 and 24-nt ta-siRNAs accumulated as the result of PpDCL3 activity. ΔPpDCL4 mutants display severe abnormalities throughout Physcomitrella development, including sterility, that were fully reversed in ΔPpDCL3/ΔPpDCL4 double-mutant plants. Therefore, PpDCL3 activity, not loss of PpDCL4 function per se, is the cause of the ΔPpDCL4 phenotypes. Additionally, we describe several new Physcomitrella trans-acting siRNA loci, three of which belong to a new family, TAS6. TAS6 loci are typified by sliced miR156 and miR529 target sites and are in close proximity to PpTAS3 loci.

RNA interference (RNAi) is a mechanism that regulates genes by either transcriptional (TGS) or posttranscriptional gene silencing (PTGS), required for genome maintenance and proper development of an organism. Small non-coding RNAs are the key players in RNAi and have been intensively studied in eukaryotes. In plants, several classes of small RNAs with specific sizes and dedicated functions have evolved. The major classes of small RNAs include microRNAs (miRNAs) and small interfering RNAs (siRNAs), which differ in their biogenesis. miRNAs are synthesized from a short hairpin structure while siRNAs are derived from long double-stranded RNAs (dsRNA). Both miRNA and siRNAs control the expression of cognate target RNAs by binding to reverse complementary sequences mediating cleavage or translational inhibition of the target RNA. They also act on the DNA and cause epigenetic changes such as DNA methylation and histone modifications. In the last years, the analysis of plant RNAi pathways was extended to the bryophyte Physcomitrella patens, a non-flowering, non-vascular ancient land plant that diverged from the lineage of seed plants approximately 450 million years ago. Based on a number of characteristic features and its phylogenetic key position in land plant evolution P. patens emerged as a plant model species to address basic as well as applied topics in plant biology. Here we summarize the current knowledge on the role of RNAi in P. patens that shows functional overlap with RNAi pathways from seed plants, and also unique features specific to this species.

ORCID IDs: 0000-0003-2740-5991 (N.M.); 0000-0003-2973-9514 (J.M.)The seedling-lethal Arabidopsis thaliana high chlorophyll fluorescence145 (hcf145) mutation leads to reduced stability of the plastid tricistronic psaA-psaB-rps14 mRNA and photosystem I (PSI) deficiency. Here, we genetically mapped the HCF145 gene, which encodes a plant-specific, chloroplast-localized, modular protein containing two homologous domains related to the polyketide cyclase family comprising 37 annotated Arabidopsis proteins of unknown function. Two further highly conserved and previously uncharacterized tandem repeat motifs at the C terminus, herein designated the transcript binding motif repeat (TMR) domains, confer sequence-specific RNA binding capability to HCF145. Homologous TMR motifs are often found as multiple repeats in quite diverse proteins of green and red algae and in the cyanobacterium Microcoleus sp PCC 7113 with unknown function. HCF145 represents the only TMR protein found in vascular plants. Detailed analysis of hcf145 mutants in Arabidopsis and Physcomitrella patens as well as in vivo and in vitro RNA binding assays indicate that HCF145 has been recruited in embryophyta for the stabilization of the psaA-psaB-rps14 mRNA via specific binding to its 59 untranslated region. The polyketide cyclase-related motifs support association of the TMRs to the psaA RNA, presumably pointing to a regulatory role in adjusting PSI levels according to the requirements of the plant cell.

Background: Cold stress causes dynamic changes in gene expression that are partially caused by small non-coding RNAs since they regulate protein coding transcripts and act in epigenetic gene silencing pathways. Thus, a detailed analysis of transcriptional changes of small RNAs (sRNAs) belonging to all known sRNA classes such as microRNAs (miRNA) and small interfering RNA (siRNAs) in response to cold contributes to an understanding of cold-related transcriptome changes. Result: We subjected A. thaliana plants to cold acclimation conditions (4°C) and analyzed the sRNA transcriptomes after 3 h, 6 h and 2 d. We found 93 cold responsive differentially expressed miRNAs and only 14 of these were previously shown to be cold responsive. We performed miRNA target prediction for all differentially expressed miRNAs and a GO analysis revealed the overrepresentation of miRNA-targeted transcripts that code for proteins acting in transcriptional regulation. We also identified a large number of differentially expressed cisand trans-nat-siRNAs, as well as sRNAs that are derived from long non-coding RNAs. By combining the results of sRNA and mRNA profiling with miRNA target predictions and publicly available information on transcription factors, we reconstructed a cold-specific, miRNA and transcription factor dependent gene regulatory network. We verified the validity of links in the network by testing its ability to predict target gene expression under cold acclimation. Conclusion: In A. thaliana, miRNAs and sRNAs derived from cisand trans-NAT gene pairs and sRNAs derived from lncRNAs play an important role in regulating gene expression in cold acclimation conditions. This study provides a fundamental database to deepen our knowledge and understanding of regulatory networks in cold acclimation.