A novel DCL2-dependent miRNA pathway in tomato affects susceptibility to RNA viruses

Wang et al. show that Dicer-like 2 (DCL2) is the major Dicer in tomato defense against tobacco mosaic virus (TMV) and potato virus X (PVX) and that it is involved in the biogenesis of endogenous 22-nt sRNA.

Tomato Dicer-like2 (slDCL2) is a key component of resistance pathways against potato virus X (PVX) and tobacco mosaic virus (TMV). It is also required for production of endogenous small RNAs, including miR6026 and other noncanonical microRNAs (miRNAs). The slDCL2 mRNAs are targets of these slDCL2-dependent RNAs in a feedback loop that was disrupted by target mimic RNAs of miR6026. In lines expressing these RNAs, there was correspondingly enhanced resistance against PVX and TMV. These findings illustrate a novel miRNA pathway in plants and a crop protection strategy in which miRNA target mimicry elevates expression of defense-related mRNAs.

Nucleosomes and DNA methylation shape meiotic DSB frequency in Arabidopsis thaliana transposons and gene regulatory regions

Meiotic recombination initiates from DNA double-strand breaks (DSBs) generated by SPO11 topoisomerase-like complexes. Meiotic DSB frequency varies extensively along eukaryotic chromosomes, with hotspots controlled by chromatin and DNA sequence. To map meiotic DSBs throughout a plant genome, we purified and sequenced Arabidopsis thaliana SPO11-1-oligonucleotides. SPO11-1-oligos are elevated in gene promoters, terminators, and introns, which is driven by AT-sequence richness that excludes nucleosomes and allows SPO11-1 access. A positive relationship was observed between SPO11-1-oligos and crossovers genome-wide, although fine-scale correlations were weaker. This may reflect the influence of interhomolog polymorphism on crossover formation, downstream from DSB formation. Although H3K4me3 is enriched in proximity to SPO11-1-oligo hotspots at gene 5' ends, H3K4me3 levels do not correlate with DSBs. Repetitive transposons are thought to be recombination silenced during meiosis, to prevent nonallelic interactions and genome instability. Unexpectedly, we found high SPO11-1-oligo levels in nucleosome-depleted Helitron/Pogo/Tc1/Mariner DNA transposons, whereas retrotransposons were coldspots. High SPO11-1-oligo transposons are enriched within gene regulatory regions and in proximity to immunity genes, suggesting a role as recombination enhancers. As transposon mobility in plant genomes is restricted by DNA methylation, we used the met1 DNA methyltransferase mutant to investigate the role of heterochromatin in SPO11-1-oligo distributions. Epigenetic activation of meiotic DSBs in proximity to centromeres and transposons occurred in met1 mutants, coincident with reduced nucleosome occupancy, gain of transcription, and H3K4me3. Together, our work reveals a complex relationship between chromatin and meiotic DSBs within A. thaliana genes and transposons, with significance for the diversity and evolution of plant genomes.

Nucleosomes and DNA methylation shape meiotic DSB frequency in Arabidopsis thaliana transposons and gene regulatory regions

Meiotic recombination initiates from DNA double-strand breaks (DSBs) generated by SPO11 topoisomerase-like complexes. Meiotic DSB frequency varies extensively along eukaryotic chromosomes, with hotspots controlled by chromatin and DNA sequence. To map meiotic DSBs throughout a plant genome, we purified and sequenced Arabidopsis thaliana SPO11-1-oligonucleotides. SPO11-1-oligos are elevated in gene promoters, terminators, and introns, which is driven by AT-sequence richness that excludes nucleosomes and allows SPO11-1 access. A positive relationship was observed between SPO11-1-oligos and crossovers genome-wide, although fine-scale correlations were weaker. This may reflect the influence of interhomolog polymorphism on crossover formation, downstream from DSB formation. Although H3K4me3 is enriched in proximity to SPO11-1-oligo hotspots at gene 5' ends, H3K4me3 levels do not correlate with DSBs. Repetitive transposons are thought to be recombination silenced during meiosis, to prevent nonallelic interactions and genome instability. Unexpectedly, we found high SPO11-1-oligo levels in nucleosome-depleted Helitron/Pogo/Tc1/Mariner DNA transposons, whereas retrotransposons were coldspots. High SPO11-1-oligo transposons are enriched within gene regulatory regions and in proximity to immunity genes, suggesting a role as recombination enhancers. As transposon mobility in plant genomes is restricted by DNA methylation, we used the met1 DNA methyltransferase mutant to investigate the role of heterochromatin in SPO11-1-oligo distributions. Epigenetic activation of meiotic DSBs in proximity to centromeres and transposons occurred in met1 mutants, coincident with reduced nucleosome occupancy, gain of transcription, and H3K4me3. Together, our work reveals a complex relationship between chromatin and meiotic DSBs within A. thaliana genes and transposons, with significance for the diversity and evolution of plant genomes.

Regulation of genome-wide DNA methylation by mobile small RNAs

Contents Summary 540 I. Introduction 540 II. There are different types of sRNA mobility 541 III. Mechanisms of sRNA movement 541 IV. Long-distance, shoot-root, mobile siRNAs influence DNA methylation in recipient tissues 541 V. Classes of interactions between shoot-root mobile siRNAs and DNA methylation 542 VI. Loci targeted directly and indirectly by shoot-root mobile siRNAs are associated with different histone modifications 543 VII. Is mobile siRNA-regulated DNA methylation important in specific tissues or under specific conditions? 543 VIII. Mobile sRNAs can be used to modify plant traits 544 IX. Conclusions 544 Acknowledgements 544 References 544 SUMMARY: RNA-directed DNA methylation (RdDM) at cytosine residues regulates gene expression, silences transposable elements and influences genome stability. The mechanisms responsible for RdDM are guided to target loci by small RNAs (sRNAs) that can move within plants cell to cell and long distance. Here we discuss recent advances in the understanding of interactions between mobile sRNAs and DNA methylation. We describe the mechanisms of sRNA movement, the differences between known classes of mobile sRNA-DNA methylation interactions and the limits of current knowledge. Finally, we discuss potential applications of mobile sRNAs in modifying plant traits.

Methods for discovering genomic loci exhibiting complex patterns of differential methylation

BACKGROUND

Cytosine methylation is widespread in most eukaryotic genomes and is known to play a substantial role in various regulatory pathways. Unmethylated cytosines may be converted to uracil through the addition of sodium bisulphite, allowing genome-wide quantification of cytosine methylation via high-throughput sequencing. The data thus acquired allows the discovery of methylation 'loci'; contiguous regions of methylation consistently methylated across biological replicates. The mapping of these loci allows for associations with other genomic factors to be identified, and for analyses of differential methylation to take place.

RESULTS

The segmentSeq R package is extended to identify methylation loci from high-throughput sequencing data from multiple experimental conditions. A statistical model is then developed that accounts for biological replication and variable rates of non-conversion of cytosines in each sample to compute posterior likelihoods of methylation at each locus within an empirical Bayesian framework. The same model is used as a basis for analysis of differential methylation between multiple experimental conditions with the baySeq R package. We demonstrate the capability of this method to analyse complex data sets in an analysis of data derived from multiple Dicer-like mutants in Arabidopsis. This reveals several novel behaviours at distinct sets of loci in response to loss of one or more of the Dicer-like proteins that indicate an antagonistic relationship between the Dicer-like proteins at at least some methylation loci. Finally, we show in simulation studies that this approach can be significantly more powerful in the detection of differential methylation than many existing methods in data derived from both mammalian and plant systems.

CONCLUSIONS

The methods developed here make it possible to analyse high-throughput sequencing of the methylome of any given organism under a diverse set of experimental conditions. The methods are able to identify methylation loci and evaluate the likelihood that a region is truly methylated under any given experimental condition, allowing for downstream analyses that characterise differences between methylated and non-methylated regions of the genome. Futhermore, diverse patterns of differential methylation may also be characterised from these data.

Mobile small RNAs and their role in regulating cytosine methylation of DNA

Small (s)RNAs of 21 to 24 nucleotides are associated with RNA silencing and methylation of DNA cytosine residues. All sizes can move from cell-to-cell and long distance in plants, directing RNA silencing in destination cells. Twenty-four nucleotide sRNAs are the predominant long-distance mobile species. Thousands move from shoot to root, where they target methylation of transposable elements both directly and indirectly. We derive several classes of interaction between small RNAs and methylation and use these to explore the mechanisms of methylation and gene expression that associate with mobile sRNA signaling.

Recombination Rate Heterogeneity within Arabidopsis Disease Resistance Genes

Meiotic crossover frequency varies extensively along chromosomes and is typically concentrated in hotspots. As recombination increases genetic diversity, hotspots are predicted to occur at immunity genes, where variation may be beneficial. A major component of plant immunity is recognition of pathogen Avirulence (Avr) effectors by resistance (R) genes that encode NBS-LRR domain proteins. Therefore, we sought to test whether NBS-LRR genes would overlap with meiotic crossover hotspots using experimental genetics in Arabidopsis thaliana. NBS-LRR genes tend to physically cluster in plant genomes; for example, in Arabidopsis most are located in large clusters on the south arms of chromosomes 1 and 5. We experimentally mapped 1,439 crossovers within these clusters and observed NBS-LRR gene associated hotspots, which were also detected as historical hotspots via analysis of linkage disequilibrium. However, we also observed NBS-LRR gene coldspots, which in some cases correlate with structural heterozygosity. To study recombination at the fine-scale we used high-throughput sequencing to analyze ~1,000 crossovers within the RESISTANCE TO ALBUGO CANDIDA1 (RAC1) R gene hotspot. This revealed elevated intragenic crossovers, overlapping nucleosome-occupied exons that encode the TIR, NBS and LRR domains. The highest RAC1 recombination frequency was promoter-proximal and overlapped CTT-repeat DNA sequence motifs, which have previously been associated with plant crossover hotspots. Additionally, we show a significant influence of natural genetic variation on NBS-LRR cluster recombination rates, using crosses between Arabidopsis ecotypes. In conclusion, we show that a subset of NBS-LRR genes are strong hotspots, whereas others are coldspots. This reveals a complex recombination landscape in Arabidopsis NBS-LRR genes, which we propose results from varying coevolutionary pressures exerted by host-pathogen relationships, and is influenced by structural heterozygosity.

DNA methylation epigenetically silences crossover hot spots and controls chromosomal domains of meiotic recombination in Arabidopsis

Yelina et al. show that RNA-directed DNA methylation is sufficient to locally silence Arabidopsis euchromatic crossover hot spots and is associated with increased nucleosome density and H3K9me2. This work demonstrates that DNA methylation plays a key role in establishing domains of meiotic recombination along chromosomes.

During meiosis, homologous chromosomes undergo crossover recombination, which is typically concentrated in narrow hot spots that are controlled by genetic and epigenetic information. Arabidopsis chromosomes are highly DNA methylated in the repetitive centromeres, which are also crossover-suppressed. Here we demonstrate that RNA-directed DNA methylation is sufficient to locally silence Arabidopsis euchromatic crossover hot spots and is associated with increased nucleosome density and H3K9me2. However, loss of CG DNA methylation maintenance in met1 triggers epigenetic crossover remodeling at the chromosome scale, with pericentromeric decreases and euchromatic increases in recombination. We used recombination mutants that alter interfering and noninterfering crossover repair pathways (fancm and zip4) to demonstrate that remodeling primarily involves redistribution of interfering crossovers. Using whole-genome bisulfite sequencing, we show that crossover remodeling is driven by loss of CG methylation within the centromeric regions. Using cytogenetics, we profiled meiotic DNA double-strand break (DSB) foci in met1 and found them unchanged relative to wild type. We propose that met1 chromosome structure is altered, causing centromere-proximal DSBs to be inhibited from maturation into interfering crossovers. These data demonstrate that DNA methylation is sufficient to silence crossover hot spots and plays a key role in establishing domains of meiotic recombination along chromosomes.

The use of duplex-specific nuclease in ribosome profiling and a user-friendly software package for Ribo-seq data analysis

Ribosome profiling is a technique that permits genome-wide, quantitative analysis of translation and has found broad application in recent years. Here we describe a modified profiling protocol and software package designed to benefit more broadly the translation community in terms of simplicity and utility. The protocol, applicable to diverse organisms, including organelles, is based largely on previously published profiling methodologies, but uses duplex-specific nuclease (DSN) as a convenient, species-independent way to reduce rRNA contamination. We show that DSN-based depletion compares favorably with other commonly used rRNA depletion strategies and introduces little bias. The profiling protocol typically produces high levels of triplet periodicity, facilitating the detection of coding sequences, including upstream, downstream, and overlapping open reading frames (ORFs) and an alternative ribosome conformation evident during termination of protein synthesis. In addition, we provide a software package that presents a set of methods for parsing ribosomal profiling data from multiple samples, aligning reads to coding sequences, inferring alternative ORFs, and plotting average and transcript-specific aspects of the data. Methods are also provided for extracting the data in a form suitable for differential analysis of translation and translational efficiency.

Sex Chromosome Dosage Compensation in Heliconius Butterflies: Global yet Still Incomplete?

The evolution of heterogametic sex chromosomes is often—but not always—accompanied by the evolution of dosage compensating mechanisms that mitigate the impact of sex-specific gene dosage on levels of gene expression. One emerging view of this process is that such mechanisms may only evolve in male-heterogametic (XY) species but not in female-heterogametic (ZW) species, which will consequently exhibit “incomplete” sex chromosome dosage compensation. However, recent results suggest that at least some Lepidoptera (moths and butterflies) may prove to be an exception to this prediction. Studies in bombycoid moths indicate the presence of a chromosome-wide epigenetic mechanism that effectively balances Z chromosome gene expression between the sexes by reducing Z-linked expression in males. In contrast, strong sex chromosome dosage effects without any reduction in male Z-linked expression were previously reported in a pyralid moth, suggesting a lack of any such dosage compensating mechanism. Here we report an analysis of sex chromosome dosage compensation in Heliconius butterflies, sampling multiple individuals for several different adult tissues (head, abdomen, leg, mouth, and antennae). Methodologically, we introduce a novel application of linear mixed-effects models to assess dosage compensation, offering a unified statistical framework that can estimate effects specific to chromosome, to sex, and their interactions (i.e., a dosage effect). Our results show substantially reduced Z-linked expression relative to autosomes in both sexes, as previously observed in bombycoid moths. This observation is consistent with an increasing body of evidence that some lepidopteran species possess an epigenetic dosage compensating mechanism that reduces Z chromosome expression in males to levels comparable with females. However, this mechanism appears to be imperfect in Heliconius, resulting in a modest dosage effect that produces an average 5–20% increase in male expression relative to females on the Z chromosome, depending on the tissue. Thus our results in Heliconius reflect a mixture of previous patterns reported for Lepidoptera. In Heliconius, a moderate pattern of incomplete dosage compensation persists apparently despite the presence of an epigenetic dosage compensating mechanism. The chromosomal distributions of sex-biased genes show an excess of male-biased and a dearth of female-biased genes on the Z chromosome relative to autosomes, consistent with predictions of sexually antagonistic evolution.

Arabidopsis meiotic crossover hot spots overlap with H2A.Z nucleosomes at gene promoters

PRDM9 directs human meiotic crossover hot spots to intergenic sequence motifs, whereas budding yeast hot spots overlap regions of low nucleosome density (LND) in gene promoters. To investigate hot spots in plants, which lack PRDM9, we used coalescent analysis of genetic variation in Arabidopsis thaliana. Crossovers increased toward gene promoters and terminators, and hot spots were associated with active chromatin modifications, including H2A.Z, histone H3 Lys4 trimethylation (H3K4me3), LND and low DNA methylation. Hot spot-enriched A-rich and CTT-repeat DNA motifs occurred upstream and downstream, respectively, of transcriptional start sites. Crossovers were asymmetric around promoters and were most frequent over CTT-repeat motifs and H2A.Z nucleosomes. Pollen typing, segregation and cytogenetic analysis showed decreased numbers of crossovers in the arp6 H2A.Z deposition mutant at multiple scales. During meiosis, H2A.Z forms overlapping chromosomal foci with the DMC1 and RAD51 recombinases. As arp6 reduced the number of DMC1 or RAD51 foci, H2A.Z may promote the formation or processing of meiotic DNA double-strand breaks. We propose that gene chromatin ancestrally designates hot spots within eukaryotes and PRDM9 is a derived state within vertebrates.

High-throughput sequencing of cytosine methylation in plant DNA

: Cytosine methylation is a significant and widespread regulatory factor in plant systems. Methods for the high-throughput sequencing of methylation have allowed a greatly improved characterisation of the methylome. Here we discuss currently available methods for generation and analysis of high-throughput sequencing of methylation data. We also discuss the results previously acquired through sequencing plant methylomes, and highlight remaining challenges in this field.

Empirical Bayesian analysis of paired high-throughput sequencing data with a beta-binomial distribution

BACKGROUND

Pairing of samples arises naturally in many genomic experiments; for example, gene expression in tumour and normal tissue from the same patients. Methods for analysing high-throughput sequencing data from such experiments are required to identify differential expression, both within paired samples and between pairs under different experimental conditions.

RESULTS

We develop an empirical Bayesian method based on the beta-binomial distribution to model paired data from high-throughput sequencing experiments. We examine the performance of this method on simulated and real data in a variety of scenarios. Our methods are implemented as part of the RbaySeq package (versions 1.11.6 and greater) available from Bioconductor (http://www.bioconductor.org).

CONCLUSIONS

We compare our approach to alternatives based on generalised linear modelling approaches and show that our method offers significant gains in performance on simulated data. In testing on real data from oral squamous cell carcinoma patients, we discover greater enrichment of previously identified head and neck squamous cell carcinoma associated gene sets than has previously been achieved through a generalised linear modelling approach, suggesting that similar gains in performance may be found in real data. Our methods thus show real and substantial improvements in analyses of high-throughput sequencing data from paired samples.

Metastable differentially methylated regions within Arabidopsis inbred populations are associated with modified expression of non-coding transcripts

Individual plants within a population may vary at both genetic and epigenetic levels. The rate of genetic divergence and its underlying mechanisms is well understood. Less is known about the factors contributing to epigenetic divergence among isogenic populations except that, despite the presence of mechanisms that faithfully maintain epigenetic marks, epigenetic differences are more frequent than genetic variation. Epigenetically divergent stretches of isogenic DNA sequence are called epialleles. Currently, it is not clear why certain regions exhibit variable epigenetic status. We identified and characterised two long RNA transcripts with altered expression and DNA methylation in an ago5 mutant. However, further investigation revealed that these changes were not dependent upon AGO5. Rather, the variable transcription of these loci in Arabidopsis mutant and wild-type populations corresponds to spontaneous differential methylated regions (DMRs) or epialleles. These two DMRs are delineated by RNAs which are highly expressed when the DMR is hypomethylated. Furthermore, they control the expression of 5′ transcriptional start site mRNA variants of nearby protein coding genes. Our data support the recent observations that meiotically stable DMRs exist within inbred populations. We further demonstrate that DMR boundaries can be defined by putative non-coding promoter-associated transcripts.

Identifying small interfering RNA loci from high-throughput sequencing data

MOTIVATION

Small interfering RNAs (siRNAs) are produced from much longer sequences of double-stranded RNA precursors through cleavage by Dicer or a Dicer-like protein. These small RNAs play a key role in genetic and epigenetic regulation; however, a full understanding of the mechanisms by which they operate depends on the characterization of the precursors from which they are derived.

RESULTS

High-throughput sequencing of small RNA populations allows the locations of the double-stranded RNA precursors to be inferred. We have developed methods to analyse small RNA sequencing data from multiple biological sources, taking into account replicate information, to identify robust sets of siRNA precursors. Our methods show good performance on both a set of small RNA sequencing data in Arabidopsis thaliana and simulated datasets.

AVAILABILITY

Our methods are available as the Bioconductor (www.bioconductor.org) package segmentSeq (version 1.5.6 and above).

Getting a full dose? Reconsidering sex chromosome dosage compensation in the silkworm, Bombyx mori

Dosage compensation—equalizing gene expression levels in response to differences in gene dose or copy number—is classically considered to play a critical role in the evolution of heteromorphic sex chromosomes. As the X and Y diverge through degradation and gene loss on the Y (or the W in female-heterogametic ZW taxa), it is expected that dosage compensation will evolve to correct for sex-specific differences in gene dose. Although this is observed in some organisms, recent genome-wide expression studies in other taxa have revealed striking exceptions. In particular, reports that both birds and the silkworm moth (Bombyx mori) lack dosage compensation have spurred speculation that this is the rule for all female-heterogametic taxa. Here, we revisit the issue of dosage compensation in silkworm by replicating and extending the previous analysis. Contrary to previous reports, our efforts reveal a pattern typically associated with dosage compensated taxa: the global male:female expression ratio does not differ between the Z and autosomes. We believe the previous report of unequal male:female ratios on the Z reflects artifacts of microarray normalization in conjunction with not testing a major assumption that the male:female global expression ratio was unbiased for autosomal loci. However, we also find that the global Z chromosome expression is significantly reduced relative to autosomes, a pattern not expected in dosage compensated taxa. This combination of male:female parity with an overall reduction in expression for sex-linked loci is not consistent with the prevailing evolutionary theory of sex chromosome evolution and dosage compensation.

baySeq: empirical Bayesian methods for identifying differential expression in sequence count data

Background

High throughput sequencing has become an important technology for studying expression levels in many types of genomic, and particularly transcriptomic, data. One key way of analysing such data is to look for elements of the data which display particular patterns of differential expression in order to take these forward for further analysis and validation.

Results

We propose a framework for defining patterns of differential expression and develop a novel algorithm, baySeq, which uses an empirical Bayes approach to detect these patterns of differential expression within a set of sequencing samples. The method assumes a negative binomial distribution for the data and derives an empirically determined prior distribution from the entire dataset. We examine the performance of the method on real and simulated data.

Conclusions

Our method performs at least as well, and often better, than existing methods for analyses of pairwise differential expression in both real and simulated data. When we compare methods for the analysis of data from experimental designs involving multiple sample groups, our method again shows substantial gains in performance. We believe that this approach thus represents an important step forward for the analysis of count data from sequencing experiments.

Small silencing RNAs in plants are mobile and direct epigenetic modification in recipient cells

A silencing signal in plants with an RNA specificity determinant moves through plasmodesmata and the phloem. To identify the mobile RNA, we grafted Arabidopsis thaliana shoots to roots that would be a recipient for the silencing signal. Using mutants that block small RNA (sRNA) biogenesis in either source or recipient tissue, we found that transgene-derived sRNA as well as a substantial proportion of the endogenous sRNA had moved across the graft union, and we provide evidence that 24-nucleotide mobile sRNAs direct epigenetic modifications in the genome of the recipient cells. Mobile sRNA thus represents a mechanism for transmitting the specification of epigenetic modification and could affect genome defense and responses to external stimuli that have persistent effects in plants.

The Arabidopsis RNA-directed DNA methylation argonautes functionally diverge based on their expression and interaction with target loci

Argonaute (AGO) effectors of RNA silencing bind small RNA (sRNA) molecules and mediate mRNA cleavage, translational repression, or epigenetic DNA modification. In many organisms, these targeting mechanisms are devolved to different products of AGO multigene families. To investigate the basis of AGO functional diversification, we characterized three closely related Arabidopsis thaliana AGOs (AGO4, AGO6, and AGO9) implicated in RNA-directed DNA methylation. All three AGOs bound 5' adenosine 24-nucleotide sRNAs, but each exhibited different preferences for sRNAs from different heterochromatin-associated loci. This difference was reduced when AGO6 and AGO9 were expressed from the AGO4 promoter, indicating that the functional diversification was partially due to differential expression of the corresponding genes. However, the AGO4-directed pattern of sRNA accumulation and DNA methylation was not fully recapitulated with AGO6 or AGO9 expressed from the AGO4 promoter. Here, we show that sRNA length and 5' nucleotide do not account for the observed functional diversification of these AGOs. Instead, the selectivity of sRNA binding is determined by the coincident expression of the AGO and sRNA-generating loci, and epigenetic modification is influenced by interactions between the AGO protein and the different target loci. These findings highlight the importance of tissue specificity and AGO-associated proteins in influencing epigenetic modifications.

Analysis of small RNA in fission yeast; centromeric siRNAs are potentially generated through a structured RNA

The formation of heterochromatin at the centromeres in fission yeast depends on transcription of the outer repeats. These transcripts are processed into siRNAs that target homologous loci for heterochromatin formation. Here, high throughput sequencing of small RNA provides a comprehensive analysis of centromere-derived small RNAs. We found that the centromeric small RNAs are Dcr1 dependent, carry 5'-monophosphates and are associated with Ago1. The majority of centromeric small RNAs originate from two remarkably well-conserved sequences that are present in all centromeres. The high degree of similarity suggests that this non-coding sequence in itself may be of importance. Consistent with this, secondary structure-probing experiments indicate that this centromeric RNA is partially double-stranded and is processed by Dicer in vitro. We further demonstrate the existence of small centromeric RNA in rdp1Delta cells. Our data suggest a pathway for siRNA generation that is distinct from the well-documented model involving RITS/RDRC. We propose that primary transcripts fold into hairpin-like structures that may be processed by Dcr1 into siRNAs, and that these siRNAs may initiate heterochromatin formation independent of RDRC activity.