DICER-LIKE1a autoregulation based on intronic microRNA processing is required for stress adaptation in Physcomitrium patens.

The Physcomitrium patens DICER-LIKE1a (PpDCL1a) mRNA encoding the essential Dicer protein for microRNA (miRNA) biogenesis harbors an intronic miRNA (miR1047). An autoregulatory mechanism to control PpDCL1a abundance that is based on competitive processing of the intronic miRNA and proper PpDCL1a mRNA splicing has previously been proposed. If intron splicing occurs first the mRNA can be translated into the functional PpDCL1a protein, whereas the processing of the intronic miRNA catalyzed by PpDCL1a itself, prior to pre-mRNA splicing, generates a truncated transcript unable to produce a functional protein. This proposed autoregulation of DCL1 has not been functionally analyzed in any plant species, and the existence of this autoregulatory control is expected to have a general impact on the overall miRNA biogenesis pathway and the transcriptome that is under miRNA control. We abolished PpDCL1a autoregulatory feedback control by the precise deletion of the MIR1047-containing intron. The generated line displayed hypersensitivity to salt stress and hyposensitivity to the plant hormone ABA, accompanied by the disturbed expression of miRNAs and mRNAs, revealed by transcriptome analyses. The feedback control together with the phenotypic abnormalities and molecular changes in the intron-less line can be rescued by the re-insertion of a modified intron harboring a sequence-unrelated artificial miRNA. Our findings indicate the physiological importance of miR1047-based feedback control of PpDCL1a transcript abundance, which controls the expression of miRNAs, and their cognate target RNAs during salt stress adaptation, and suggests a key role for this autoregulation in the molecular adaptation of land plants to terrestrial habitats.

Identification of Small RNAs During High Light Acclimation in Arabidopsis thaliana.

The biological significance of non-coding RNAs (ncRNAs) has been firmly established to be important for the regulation of genes involved in stress acclimation. Light plays an important role for the growth of plants providing the energy for photosynthesis; however, excessive light conditions can also cause substantial defects. Small RNAs (sRNAs) are a class of non-coding RNAs that regulate transcript levels of protein-coding genes and mediate epigenetic silencing. Next generation sequencing facilitates the identification of small non-coding RNA classes such as miRNAs (microRNAs) and small-interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), but changes in the ncRNA transcriptome in response to high light are poorly understood. We subjected Arabidopsis plants to high light conditions and performed a temporal in-depth study of the transcriptome data after 3 h, 6 h, and 2 days of high light treatment. We identified a large number of high light responsive miRNAs and sRNAs derived from NAT gene pairs, lncRNAs and TAS transcripts. We performed target predictions for differentially expressed miRNAs and correlated their expression levels through mRNA sequencing data. GO analysis of the targets revealed an overrepresentation of genes involved in transcriptional regulation. In A. thaliana, sRNA-mediated regulation of gene expression in response to high light treatment is mainly carried out by miRNAs and sRNAs derived from NAT gene pairs, and from lncRNAs. This study provides a deeper understanding of sRNA-dependent regulatory networks in high light acclimation.

PpGRAS12 acts as a positive regulator of meristem formation in Physcomitrium patens.

KEY MESSAGE: This study focused on the key regulatory function of Physcomitrium patens GRAS12 gene underlying an increasing plant complexity, an important step in plant terrestrialization and the evolutionary history of life. The miR171-GRAS module has been identified as a key player in meristem maintenance in angiosperms. PpGRAS12 is a member of the GRAS family and a validated target for miR171 in Physcomitrium (Physcomitrella) patens. Here we show a regulatory function of miR171 at the gametophytic vegetative growth stage and targeted deletion of the PpGRAS12 gene adversely affects sporophyte production since fewer sporophytes were produced in ΔPpGRAS12 knockout lines compared to wild type moss. Furthermore, highly specific and distinct growth arrests were observed in inducible PpGRAS12 overexpression lines at the protonema stage. Prominent phenotypic aberrations including the formation of multiple apical meristems at the gametophytic vegetative stage in response to elevated PpGRAS12 transcript levels were discovered via scanning electron microscopy. The production of multiple buds in the PpGRAS12 overexpression lines similar to ΔPpCLV1a/1b disruption mutants is accompanied by an upregulation of PpCLE and downregulation of PpCLV1, PpAPB, PpNOG1, PpDEK1, PpRPK2 suggesting that PpGRAS12 acts upstream of these genes and negatively regulates the proposed pathway to specify simplex meristem formation. As CLV signaling pathway components are not present in the chlorophytic or charophytic algae and arose with the earliest land plants, we identified a key regulatory function of PpGRAS12 underlying an increasing plant complexity, an important step in plant terrestrialization and the evolutionary history of life.

Identification of small non-coding RNAs responsive to GUN1 and GUN5 related retrograde signals in Arabidopsis thaliana.

Chloroplast perturbations activate retrograde signalling pathways, causing dynamic changes of gene expression. Besides transcriptional control of gene expression, different classes of small non-coding RNAs (sRNAs) act in gene expression control, but comprehensive analyses regarding their role in retrograde signalling are lacking. We performed sRNA profiling in response to norflurazon (NF), which provokes retrograde signals, in Arabidopsis thaliana wild type (WT) and the two retrograde signalling mutants gun1 and gun5. The RNA samples were also used for mRNA and long non-coding RNA profiling to link altered sRNA levels to changes in the expression of their cognate target RNAs. We identified 122 sRNAs from all known sRNA classes that were responsive to NF in the WT. Strikingly, 142 and 213 sRNAs were found to be differentially regulated in both mutants, indicating a retrograde control of these sRNAs. Concomitant with the changes in sRNA expression, we detected about 1500 differentially expressed mRNAs in the NF-treated WT and around 900 and 1400 mRNAs that were differentially regulated in the gun1 and gun5 mutants, with a high proportion (~30%) of genes encoding plastid proteins. Furthermore, around 20% of predicted miRNA targets code for plastid-localised proteins. Among the sRNA-target pairs, we identified pairs with an anticorrelated expression as well pairs showing other expressional relations, pointing to a role of sRNAs in balancing transcriptional changes upon retrograde signals. Based on the comprehensive changes in sRNA expression, we assume a considerable impact of sRNAs in retrograde-dependent transcriptional changes to adjust plastidic and nuclear gene expression.

Identification of small RNAs during cold acclimation in Arabidopsis thaliana.

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 cis- and 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 cis- and 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.

ABA-Induced Vegetative Diaspore Formation in Physcomitrella patens.

The phytohormone abscisic acid (ABA) is a pivotal regulator of gene expression in response to various environmental stresses such as desiccation, salt and cold causing major changes in plant development and physiology. Here we show that in the moss Physcomitrella patens exogenous application of ABA triggers the formation of vegetative diaspores (brachycytes or brood cells) that enable plant survival in unfavorable environmental conditions. Such diaspores are round-shaped cells characterized by the loss of the central vacuole, due to an increased starch and lipid storage preparing these cells for growth upon suitable environmental conditions. To gain insights into the gene regulation underlying these developmental and physiological changes, we analyzed early transcriptome changes after 30, 60, and 180 min of ABA application and identified 1,030 differentially expressed genes. Among these, several groups can be linked to specific morphological and physiological changes during diaspore formation, such as genes involved in cell wall modifications. Furthermore, almost all members of ABA-dependent signaling and regulation were transcriptionally induced. Network analysis of transcription-associated genes revealed a large overlap of our study with ABA-dependent regulation in response to dehydration, cold stress, and UV-B light, indicating a fundamental function of ABA in diverse stress responses in moss. We also studied the evolutionary conservation of ABA-dependent regulation between moss and the seed plant Arabidopsis thaliana pointing to an early evolution of ABA-mediated stress adaptation during the conquest of the terrestrial habitat by plants.

Principles of nucleation of H3K27 methylation during embryonic development.

During embryonic development, maintenance of cell identity and lineage commitment requires the Polycomb-group PRC2 complex, which catalyzes histone H3 lysine 27 trimethylation (H3K27me3). However, the developmental origins of this regulation are unknown. Here we show that H3K27me3 enrichment increases from blastula stages onward in embryos of the Western clawed frog (Xenopus tropicalis) within constrained domains strictly defined by sequence. Strikingly, although PRC2 also binds widely to active enhancers, H3K27me3 is only deposited at a small subset of these sites. Using a Support Vector Machine algorithm, these sequences can be predicted accurately on the basis of DNA sequence alone, with a sequence signature conserved between humans, frogs, and fish. These regions correspond to the subset of blastula-stage DNA methylation-free domains that are depleted for activating promoter motifs, and enriched for motifs of developmental factors. These results imply a genetic-default model in which a preexisting absence of DNA methylation is the major determinant of H3K27 methylation when not opposed by transcriptional activation. The sequence and motif signatures reveal the hierarchical and genetically inheritable features of epigenetic cross-talk that impose constraints on Polycomb regulation and guide H3K27 methylation during the exit of pluripotency.

DICER-LIKE3 activity in Physcomitrella patens DICER-LIKE4 mutants causes severe developmental dysfunction and sterility.

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.

Expression of artificial microRNAs in Physcomitrella patens.

MicroRNAs (miRNAs) are ∼21-nt-long small RNAs transcribed from endogenous MIR genes which form precursor RNAs with a characteristic hairpin structure. MiRNAs control the expression of cognate target genes by binding to reverse complementary sequences resulting in cleavage or translational inhibition of the target RNA. Artificial miRNAs (amiRNAs) can be generated by exchanging the miRNA/miRNA sequence of endogenous MIR precursor genes, while maintaining the general pattern of matches and mismatches in the foldback. Thus, for functional gene analysis, amiRNAs can be designed to target any gene of interest. During the last decade, the moss Physcomitrella patens emerged as a model plant for functional gene analysis based on its unique ability to integrate DNA into the nuclear genome by homologous recombination which allows for the generation targeted gene knockout mutants. In addition to this, we developed a protocol to express amiRNAs in P. patens that has particular advantages over the generation of knockout mutants and might be used to speed up reverse genetics approaches in this model species.

Gene function analysis by artificial microRNAs in Physcomitrella patens.

MicroRNAs (miRNAs) are ~21 nt long small RNAs transcribed from endogenous MIR genes which form precursor RNAs with a characteristic hairpin structure. miRNAs control the expression of cognate target genes by binding to reverse complementary sequences resulting in cleavage or translational inhibition of the target RNA. Artificial miRNAs (amiRNAs) can be generated by exchanging the miRNA/miRNA sequence of endogenous MIR precursor genes, while maintaining the general pattern of matches and mismatches in the foldback. Thus, for functional gene analysis amiRNAs can be designed to target any gene of interest. During the last decade the moss Physcomitrella patens emerged as a model plant for functional gene analysis based on its unique ability to integrate DNA into the nuclear genome by homologous recombination which allows for the generation of targeted gene knockout mutants. In addition to this, we developed a protocol to express amiRNAs in P. patens that has particular advantages over the generation of knockout mutants and might be used to speed up reverse genetics approaches in this model species.

Transcriptional control of gene expression by microRNAs.

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.

Physcomitrella patens DCL3 is required for 22-24 nt siRNA accumulation, suppression of retrotransposon-derived transcripts, and normal development.

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.