Heterochromatic siRNAs Research Articles

SRNAminer: A multifunctional toolkit for next-generation sequencing small RNA data mining in plants

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or “golden” set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA.

Coding and noncoding transcriptomes of NODULIN HOMEOBOX (NDX)-deficient Arabidopsis inflorescence

Arabidopsis NODULIN HOMEOBOX (NDX) is a plant-specific transcriptional regulator whose role in small RNA biogenesis and heterochromatin homeostasis has recently been described. Here we extend our previous transcriptomic analysis to the flowering stage of development. We performed mRNA-seq and small RNA-seq measurements on inflorescence samples of wild-type and ndx1-4 mutant (WiscDsLox344A04) Arabidopsis plants. We identified specific groups of differentially expressed genes and noncoding heterochromatic siRNA (hetsiRNA) loci/regions whose transcriptional activity was significantly changed in the absence of NDX. In addition, data obtained from inflorescence were compared with seedling transcriptomics data, which revealed development-specific changes in gene expression profiles. Overall, we provide a comprehensive data source on the coding and noncoding transcriptomes of NDX-deficient Arabidopsis flowers to serve as a basis for further research on NDX function.

Dynamically expressed small RNAs, substantially driven by genomic structural variants, contribute to transcriptomic changes during tomato domestication.

Tomato has undergone extensive selections during domestication. Recent progress has shown that genomic structural variants (SVs) have contributed to gene expression dynamics during tomato domestication, resulting in changes of important traits. Here, we performed comprehensive analyses of small RNAs (sRNAs) from nine representative tomato accessions. We demonstrate that SVs substantially contribute to the dynamic expression of the three major classes of plant sRNAs: microRNAs (miRNAs), phased secondary short interfering RNAs (phasiRNAs), and 24-nucleotide heterochromatic siRNAs (hc-siRNAs). Changes in the abundance of phasiRNAs and 24-nucleotide hc-siRNAs likely contribute to the alteration of mRNA gene expression in cis during tomato domestication, particularly for genes associated with biotic and abiotic stress tolerance. We also observe that miRNA expression dynamics are associated with imprecise processing, alternative miRNA-miRNA* selections, and SVs. SVs mainly affect the expression of less-conserved miRNAs that do not have established regulatory functions or low abundant members in highly expressed miRNA families. Our data highlight different selection pressures on miRNAs compared to phasiRNAs and 24-nucleotide hc-siRNAs. Our findings provide insights into plant sRNA evolution as well as SV-based gene regulation during crop domestication. Furthermore, our dataset provides a rich resource for mining the sRNA regulatory network in tomato.

Functional specialization of monocot DCL3 and DCL5 proteins through the evolution of the PAZ domain.

Monocot DICER-LIKE3 (DCL3) and DCL5 produce distinct 24-nt small interfering RNAs (siRNAs), heterochromatic siRNAs (hc-siRNAs) and phased secondary siRNAs (phasiRNAs), respectively. The former small RNAs are linked to silencing of transposable elements and heterochromatic repeats, and the latter to reproductive processes. It is assumed that these DCLs evolved from an ancient ‘eudicot-type’ DCL3 ancestor, which may have produced both types of siRNAs. However, how functional differentiation was achieved after gene duplication remains elusive. Here, we find that monocot DCL3 and DCL5 exhibit biochemically distinct preferences for 5′ phosphates and 3′ overhangs, consistent with the structural properties of their in vivo double-stranded RNA substrates. Importantly, these distinct substrate specificities are determined by the PAZ domains of DCL3 and DCL5, which have accumulated mutations during the course of evolution. These data explain the mechanism by which these DCLs cleave their cognate substrates from a fixed end, ensuring the production of functional siRNAs. Our study also indicates how plants have diversified and optimized RNA silencing mechanisms during evolution.

The evolutionary history of small RNAs in Solanaceae.

The Solanaceae or “nightshade” family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family’s small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities.

Elevated transcription of transposable elements is accompanied by het-siRNA-driven de novo DNA methylation in grapevine embryogenic callus

BackgroundSomatic variation is a valuable source of trait diversity in clonally propagated crops. In grapevine, which has been clonally propagated worldwide for centuries, important phenotypes such as white berry colour are the result of genetic changes caused by transposable elements. Additionally, epiallele formation may play a role in determining geo-specific (‘terroir’) differences in grapes and thus ultimately in wine. This genomic plasticity might be co-opted for crop improvement via somatic embryogenesis, but that depends on a species-specific understanding of the epigenetic regulation of transposable element (TE) expression and silencing in these cultures. For this reason, we used whole-genome bisulphite sequencing, mRNA sequencing and small RNA sequencing to study the epigenetic status and expression of TEs in embryogenic callus, in comparison with leaf tissue.ResultsWe found that compared with leaf tissue, grapevine embryogenic callus cultures accumulate relatively high genome-wide CHH methylation, particularly across heterochromatic regions. This de novo methylation is associated with an abundance of transcripts from highly replicated TE families, as well as corresponding 24 nt heterochromatic siRNAs. Methylation in the TE-specific CHG context was relatively low over TEs located within genes, and the expression of TE loci within genes was highly correlated with the expression of those genes.ConclusionsThis multi-‘omics analysis of grapevine embryogenic callus in comparison with leaf tissues reveals a high level of genome-wide transcription of TEs accompanied by RNA-dependent DNA methylation of these sequences in trans. This provides insight into the genomic conditions underlying somaclonal variation and epiallele formation in plants regenerated from embryogenic cultures, which is an important consideration when using these tissues for plant propagation and genetic improvement.

Dri1 mediates heterochromatin assembly via RNAi and histone deacetylation.

Heterochromatin, a transcriptionally silenced chromatin domain, is important for genome stability and gene expression. Histone 3 lysine 9 methylation (H3K9me) and histone hypoacetylation are conserved epigenetic hallmarks of heterochromatin. In fission yeast, RNA interference (RNAi) plays a key role in H3K9 methylation and heterochromatin silencing. However, how RNAi machinery and histone deacetylases (HDACs) are coordinated to ensure proper heterochromatin assembly is still unclear. Previously, we showed that Dpb4, a conserved DNA polymerase epsilon subunit, plays a key role in the recruitment of HDACs to heterochromatin during S phase. Here, we identified a novel RNA-binding protein Dri1 that interacts with Dpb4. GFP-tagged Dri1 forms distinct foci mostly in the nucleus, showing a high degree of colocalization with Swi6/Heterochromatin Protein 1. Deletion of dri1+ leads to defects in silencing, H3K9me, and heterochromatic siRNA generation. We also showed that Dri1 physically associates with heterochromatic transcripts, and is required for the recruitment of the RNA-induced transcriptional silencing (RITS) complex via interacting with the complex. Furthermore, loss of Dri1 decreases the association of the Sir2 HDAC with heterochromatin. We further demonstrated that the C-terminus of Dri1 that includes an intrinsically disordered (IDR) region and three zinc fingers is crucial for its role in silencing. Together, our evidences suggest that Dri1 facilitates heterochromatin assembly via the RNAi pathway and HDAC.

SRNAanno\u2014a database repository of uniformly annotated small RNAs in plants

Small RNAs (sRNAs) are essential regulatory molecules, and there are three major sRNA classes in plants: microRNAs (miRNAs), phased small interfering RNAs (phased siRNAs or phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). Excluding miRNAs, the other two classes are not well annotated or available in public databases for most sequenced plant genomes. We performed a comprehensive sRNA annotation of 143 plant species that have fully sequenced genomes and next-generation sequencing sRNA data publicly available. The results are available via an online repository called sRNAanno (www.plantsRNAs.org). Compared with other public plant sRNA databases, we obtained was much more miRNA annotations, which are more complete and reliable because of the consistent and highly stringent criteria used in our miRNA annotations. sRNAanno also provides free access to genomic information for >22,721 PHAS loci and >22 million hc-siRNA loci annotated from these 143 plant species. Both miRNA and PHAS loci can be easily browsed to view their main features, and a collection of archetypal trans-acting siRNA 3 (TAS3) genes were annotated separately for quick access. To facilitate the ease of sRNA annotation, sRNAanno provides free service for sRNA annotations to the community. In summary, the sRNAanno database is a great resource to facilitate genomic and genetic research on plant small RNAs.

The explant developmental stage profoundly impacts small RNA-mediated regulation at the dedifferentiation step of maize somatic embryogenesis

Maize somatic embryogenesis (SE) requires the induction of embryogenic callus and establishment of proliferation before plant regeneration. The molecular mechanisms underlying callus embryogenic potential are not well understood. Here we explored the role of small RNAs (sRNAs) and the accumulation of their target transcripts in maize SE at the dedifferentiation step using VS-535 zygotic embryos collected at distinct developmental stages and displaying contrasting in vitro embryogenic potential and morphology. MicroRNAs (miRNAs), trans-acting siRNAs (tasiRNAs), heterochromatic siRNAs (hc-siRNAs) populations and their RNA targets were analyzed by high-throughput sequencing. Abundances of specific miRNAs, tasiRNAs and targets were validated by qRT-PCR. Unique accumulation patterns were found for immature embryo at 15 Days After Pollination (DAP) and for the callus induction from this explant, as compared to 23 DAP and mature embryos. miR156, miR164, miR166, tasiARFs and the 24 nt hc-siRNAs displayed the most strikingly different patterns between explants and during dedifferentiation. According to their role in auxin responses and developmental cues, we conclude that sRNA-target regulation operating within the 15 DAP immature embryo explant provides key molecular hints as to why this stage is relevant for callus induction with successful proliferation and plant regeneration.

Several phased siRNA annotation methods can frequently misidentify 24 nucleotide siRNA-dominated PHAS loci.

Small RNAs regulate key physiological functions in land plants. Small RNAs can be divided into two categories: microRNAs (miRNAs) and short interfering RNAs (siRNAs); siRNAs are further subdivided into transposon/repetitive region‐localized heterochromatic siRNAs and phased siRNAs (phasiRNAs). PhasiRNAs are produced from the miRNA‐mediated cleavage of a Pol II RNA transcript; the miRNA cleavage site provides a defined starting point from which phasiRNAs are produced in a distinctly phased pattern. 21–22 nucleotide (nt)‐dominated phasiRNA‐producing loci (PHAS) are well represented in all land plants to date. In contrast, 24 nt‐dominated PHAS loci are known to be encoded only in monocots and are generally restricted to male reproductive tissues. Currently, only one miRNA (miR2275) is known to trigger the production of these 24 nt‐dominated PHAS loci. In this study, we use stringent methodologies in order to examine whether or not 24 nt‐dominated PHAS loci also exist in Arabidopsis thaliana. We find that highly expressed heterochromatic siRNAs were consistently misidentified as 24 nt‐dominated PHAS loci using multiple PHAS‐detecting algorithms. We also find that MIR2275 is not found in A. thaliana, and it seems to have been lost in the last common ancestor of Brassicales. Altogether, our research highlights the potential issues with widely used PHAS‐detecting algorithms which may lead to false positives when trying to annotate new PHAS, especially 24 nt‐dominated loci.

Characterization of siRNAs clusters in Arabidopsis thaliana galls induced by the root-knot nematode Meloidogyne incognita

BackgroundRoot-knot nematodes (RKN), genus Meloidogyne, are plant parasitic worms that have the ability to transform root vascular cylinder cells into hypertrophied, multinucleate and metabolically over-active feeding cells. Redifferentiation into feeding cells is the result of a massive transcriptional reprogramming of root cells targeted by RKN. Since RKN are able to induce similar feeding cells in roots of thousands of plant species, these worms are thought to manipulate essential and conserved plant molecular pathways.ResultsSmall non-coding RNAs of uninfected roots and infected root galls induced by M. incognita from Arabidopsis thaliana were sequenced by high throughput sequencing. SiRNA populations were analysed by using the Shortstack algorithm. We identified siRNA clusters that are differentially expressed in infected roots and evidenced an over-representation of the 23–24 nt siRNAs in infected tissue. This size corresponds to heterochromatic siRNAs (hc-siRNAs) which are known to regulate expression of transposons and genes at the transcriptional level, mainly by inducing DNA methylation.ConclusionsCorrelation of siRNA clusters expression profile with transcriptomic data identified several protein coding genes that are candidates to be regulated by siRNAs at the transcriptional level by RNA directed DNA methylation (RdDM) pathway either directly or indirectly via silencing of neighbouring transposable elements.

Evolution and Diversification of Small RNA Pathways in Flowering Plants.

Small regulatory RNAs guide gene silencing at the DNA or RNA level through repression of complementary sequences. The two main forms of small RNAs are microRNA (miRNA) and small interfering RNA (siRNAs), which are generated from the processing of different forms of double-stranded RNA (dsRNA) precursors. These two forms of small regulatory RNAs function in distinct but overlapping gene silencing pathways in plants. Gene silencing pathways in eukaryotes evolved from an ancient prokaryotic mechanism involved in genome defense against invasive genetic elements, but has since diversified to also play a crucial role in regulation of endogenous gene expression. Here, we review the biogenesis of the different forms of small RNAs in plants, including miRNAs, phased, secondary siRNAs (phasiRNAs) and heterochromatic siRNAs (hetsiRNAs), with a focus on their functions in genome defense, transcriptional and post-transcriptional gene silencing, RNA-directed DNA methylation, trans-chromosomal methylation and paramutation. We also discuss the important role that gene duplication has played in the functional diversification of gene silencing pathways in plants, and we highlight recently discovered components of gene silencing pathways in plants.

Ambient temperature regulates the expression of a small set of sRNAs influencing plant development through NF-YA2 and YUC2.

Plants substantially alter their developmental programme upon changes in the ambient temperature. The 21-24nt small RNAs (sRNAs) are important gene expression regulators, which play a major role in development and adaptation. However, little is known about how the different sRNA classes respond to changes in the ambient temperature. We profiled the sRNA populations in four different tissues of Arabidopsis thaliana plants grown at 15°C, 21°C, and 27°C. We found that only a small fraction (0.6%) of the sRNA loci are ambient temperature-controlled. We identified thermoresponsive microRNAs and identified their target genes using degradome libraries. We verified that the target of the thermoregulated miR169, NF-YA2, is also ambient temperature-regulated. NF-YA2, as the component of the conserved transcriptional regulator NF-Y complex, binds the promoter of the flowering time regulator FT and the auxin biosynthesis gene YUC2. Other differentially expressed loci include thermoresponsive phased siRNA loci that target various auxin pathway genes and tRNA fragments. Furthermore, a temperature-dependent 24-nt heterochromatic siRNA locus in the promoter of YUC2 may contribute to the epigenetic regulation of auxin homeostasis. This holistic approach facilitated a better understanding of the role of different sRNA classes in ambient temperature adaptation of plants.

MicroRNAs in Marchantia polymorpha.

Contents Summary 409 I. Introduction 409 II. RNA silencing machinery in Marchantia polymorpha 410 III. miRNA prediction by integrating omics approach 410 IV. miRNAs and their targets in Marchantia polymorpha 410 V. Mpo-miR390-mediated MpTAS3 tasiRNA biogenesis and potential tasiARF target MpARF2 414 VI. Artificial miRNA and CRISPR-CAS9 edited MIR gene in Marchantia polymorpha 414 VII. Conclusions 415 Acknowledgements 415 References 415 SUMMARY: The liverwort Marchantia polymorpha occupies an important phylogenetic position for comparative studies of land plant gene regulation. Multiple gene regulatory pathways mediated by small RNAs, including microRNAs (miRNAs), trans-acting short-interfering RNAs, and heterochromatic siRNAs often associated with RNA-dependent DNA methylation, have been characterized in flowering plants. Genes for essential components for all of these small RNA-mediated gene silencing pathways are found in M.polymorpha as well as the moss Phsycomitrella patens, indicating that these pathways existed in the ancestral land plant. However, only seven miRNAs are conserved across land plants, with both ancestral and novel targets identified in M.polymorpha. There is little or no evidence that any of these conserved miRNAs are present in algae. As with other plants investigated, most miRNAs in M.polypmorpha exhibit lineage-specific evolution. Application of artificial miRNA and CRISPR-Cas9 technologies in genetic studies of M.polymorpha provide avenues to further investigate miRNA biology.

AGO4 is specifically required for heterochromatic siRNA accumulation at Pol V-dependent loci in Arabidopsis thaliana.

In plants, 24 nucleotide long heterochromatic siRNAs (het-siRNAs) transcriptionally regulate gene expression by RNA-directed DNA methylation (RdDM). The biogenesis of most het-siRNAs depends on the plant-specific RNA polymerase IV (Pol IV), and ARGONAUTE4 (AGO4) is a major het-siRNA effector protein. Through genome-wide analysis of sRNA-seq data sets, we found that AGO4 is required for the accumulation of a small subset of het-siRNAs. The accumulation of AGO4-dependent het-siRNAs also requires several factors known to participate in the effector portion of the RdDM pathway, including RNA POLYMERASE V (POL V), DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2) and SAWADEE HOMEODOMAIN HOMOLOGUE 1 (SHH1). Like many AGO proteins, AGO4 is an endonuclease that can 'slice' RNAs. We found that a slicing-defective AGO4 was unable to fully recover AGO4-dependent het-siRNA accumulation from ago4 mutant plants. Collectively, our data suggest that AGO4-dependent siRNAs are secondary siRNAs dependent on the prior activity of the RdDM pathway at certain loci.

Profiling New Small RNA Sequences.

Small RNAs are key molecules in RNA silencing pathways that exert the sequence-specific regulation of gene expression and chromatin modifications in many eukaryotes. In plants, endogenous small RNAs, including microRNAs (miRNAs), trans-acting short interfering RNAs (tasiRNAs), and heterochromatic siRNAs (hc-siRNAs), play an important role in switching or orchestrating biological processes during the development and at the onset of stress responses. These endogenous and exogenous small RNAs are mainly 20-24 nucleotides in length. In addition, viral genome-derived siRNAs of similar lengths are produced during viral infection, and they exhibit anti-viral defense activity in RNA silencing pathway.Here, we introduce a method to isolate and characterize small RNA molecules possibly applicable to a wide range of plant resources and tissues. After purification from total RNAs, small RNAs were subjected to Illumina sequencing analysis using compatible reagents kits. Following the sample preparation protocol, small RNAs are ligated first at the 3'- and then at the 5'-end to the respective RNA adapters followed by reverse transcription with a set of primers to produce cDNAs with Index sequences at ends. After PCR amplification, cDNAs are subjected (after gel purification) to RNA-seq analysis. This method could be applied to isolate small RNAs from different sources and characterize small RNA profiles to compare different sets of samples, e.g., wild-type and mutant plants, plants under different stress environments, and virus-infected plants because the starting RNA material is free of contaminated starch or similar material which would block further analysis.

Improved Placement of Multi-mapping Small RNAs.

High-throughput sequencing of small RNAs (sRNA-seq) is a popular method used to discover and annotate microRNAs (miRNAs), endogenous short interfering RNAs (siRNAs), and Piwi-associated RNAs (piRNAs). One of the key steps in sRNA-seq data analysis is alignment to a reference genome. sRNA-seq libraries often have a high proportion of reads that align to multiple genomic locations, which makes determining their true origins difficult. Commonly used sRNA-seq alignment methods result in either very low precision (choosing an alignment at random), or sensitivity (ignoring multi-mapping reads). Here, we describe and test an sRNA-seq alignment strategy that uses local genomic context to guide decisions on proper placements of multi-mapped sRNA-seq reads. Tests using simulated sRNA-seq data demonstrated that this local-weighting method outperforms other alignment strategies using three different plant genomes. Experimental analyses with real sRNA-seq data also indicate superior performance of local-weighting methods for both plant miRNAs and heterochromatic siRNAs. The local-weighting methods we have developed are implemented as part of the sRNA-seq analysis program ShortStack, which is freely available under a general public license. Improved genome alignments of sRNA-seq data should increase the quality of downstream analyses and genome annotation efforts.

Genome-wide analysis of single non-templated nucleotides in plant endogenous siRNAs and miRNAs.

Plant small RNAs are subject to various modifications. Previous reports revealed widespread 3′ modifications (truncations and non-templated tailing) of plant miRNAs when the 2′-O-methyltransferase HEN1 is absent. However, non-templated nucleotides in plant heterochromatic siRNAs have not been deeply studied, especially in wild-type plants. We systematically studied non-templated nucleotide patterns in plant small RNAs by analyzing small RNA sequencing libraries from Arabidopsis, tomato, Medicago, rice, maize and Physcomitrella. Elevated rates of non-templated nucleotides were observed at the 3′ ends of both miRNAs and endogenous siRNAs from wild-type specimens of all species. ‘Off-sized’ small RNAs, such as 25 and 23 nt siRNAs arising from loci dominated by 24 nt siRNAs, often had very high rates of 3′-non-templated nucleotides. The same pattern was observed in all species that we studied. Further analysis of 24 nt siRNA clusters in Arabidopsis revealed distinct patterns of 3′-non-templated nucleotides of 23 nt siRNAs arising from heterochromatic siRNA loci. This pattern of non-templated 3′ nucleotides on 23 nt siRNAs is not affected by loss of known small RNA 3′-end modifying enzymes, and may result from modifications added to longer heterochromatic siRNA precursors.

The expanding world of small RNAs in plants.

Plant genomes encode various small RNAs that function in distinct, yet overlapping, genetic and epigenetic silencing pathways. However, the abundance and diversity of small-RNA classes varies among plant species, suggesting coevolution between environmental adaptations and gene-silencing mechanisms. Biogenesis of small RNAs in plants is well understood, but we are just beginning to uncover their intricate regulation and activity. Here, we discuss the biogenesis of plant small RNAs, such as microRNAs, secondary siRNAs and heterochromatic siRNAs, and their diverse cellular and developmental functions, including in reproductive transitions, genomic imprinting and paramutation. We also discuss the diversification of small-RNA-directed silencing pathways through the expansion of RNA-dependent RNA polymerases, DICER proteins and ARGONAUTE proteins.

Comprehensive Annotation of Physcomitrella patens Small RNA Loci Reveals That the Heterochromatic Short Interfering RNA Pathway Is Largely Conserved in Land Plants.

Many plant small RNAs are sequence-specific negative regulators of target mRNAs and/or chromatin. In angiosperms, the two most abundant endogenous small RNA populations are usually 21-nucleotide microRNAs (miRNAs) and 24-nucleotide heterochromatic short interfering RNAs (siRNAs). Heterochromatic siRNAs are derived from repetitive regions and reinforce DNA methylation at targeted loci. The existence and extent of heterochromatic siRNAs in other land plant lineages has been unclear. Using small RNA-sequencing (RNA-seq) of the moss Physcomitrella patens, we identified 1090 loci that produce mostly 23- to 24-nucleotide siRNAs. These loci are mostly in intergenic regions with dense DNA methylation. Accumulation of siRNAs from these loci depends upon P. patens homologs of DICER-LIKE3 (DCL3), RNA-DEPENDENT RNA POLYMERASE2, and the largest subunit of DNA-DEPENDENT RNA POLYMERASE IV, with the largest subunit of a Pol V homolog contributing to expression at a smaller subset of the loci. A MINIMAL DICER-LIKE (mDCL) gene, which lacks the N-terminal helicase domain typical of DCL proteins, is specifically required for 23-nucleotide siRNA accumulation. We conclude that heterochromatic siRNAs, and their biogenesis pathways, are largely identical between angiosperms and P. patens, with the notable exception of the P. patens-specific use of mDCL to produce 23-nucleotide siRNAs.