RNAi is a conserved mechanism in which small interfering RNAs (siRNAs) guide the degradation of cognate RNAs, but also promote heterochromatin assembly at repetitive DNA elements such as centromeric repeats1,2. However, the full extent of RNAi functions and its endogenous targets have not been explored. Here we show that in the fission yeast Schizosaccharomyces pombe, RNAi and heterochromatin factors cooperate to silence diverse loci, including sexual differentiation genes, genes encoding transmembrane proteins, and retrotransposons that are also targeted by the exosome RNA degradation machinery. In the absence of the exosome, transcripts are processed preferentially by the RNAi machinery, revealing siRNA clusters and corresponding increase in heterochromatin modifications across large domains containing genes and retrotransposons. We show that the generation of siRNAs and heterochromatin assembly by RNAi is triggered by a mechanism involving the canonical poly(A) polymerase Pla1 and an associated RNA surveillance factor Red1, which also activate the exosome. Remarkably, siRNA production and heterochromatin modifications at these target loci are regulated by environmental growth conditions, and by developmental signals that induce gene expression during sexual differentiation. Our analyses uncover interplay between RNAi and the exosome that is conserved in higher eukaryotes, and show that differentiation signals modulate RNAi silencing to regulate developmental genes.

T4 RNA ligases are commonly used to attach adapters to RNAs, but large differences in ligation efficiency make detection and quantitation problematic. We developed a ligation selection strategy using random RNAs in combination with high-throughput sequencing to gain insight into the differences in efficiency of ligating pre-adenylated DNA adapters to RNA 3′-ends. After analyzing biases in RNA sequence, secondary structure and RNA-adapter cofold structure, we conclude that T4 RNA ligases do not show significant primary sequence preference in RNA substrates, but are biased against structural features within RNAs and adapters. Specifically, RNAs with less than three unstructured nucleotides at the 3′-end and RNAs that are predicted to cofold with an adapter in unfavorable structures are likely to be poorly ligated. The effect of RNA-adapter cofold structures on ligation is supported by experiments where the ligation efficiency of specific miRNAs was changed by designing adapters to alter cofold structure. In addition, we show that using adapters with randomized regions results in higher ligation efficiency and reduced ligation bias. We propose that using randomized adapters may improve RNA representation in experiments that include a 3′-adapter ligation step.

High-throughput sequencing (HTS) has become a powerful tool for the detection of and sequence characterization of microRNAs (miRNA) and other small RNAs (sRNA). Unfortunately, the use of HTS data to determine the relative quantity of different miRNAs in a sample has been shown to be inconsistent with quantitative PCR and Northern Blot results. Several recent studies have concluded that the major contributor to this inconsistency is bias introduced during the construction of sRNA libraries for HTS and that the bias is primarily derived from the adaptor ligation steps, specifically where single stranded adaptors are sequentially ligated to the 3’ and 5’-end of sRNAs using T4 RNA ligases. In this study we investigated the effects of ligation bias by using a pool of randomized ligation substrates, defined mixtures of miRNA sequences and several combinations of adaptors in HTS library construction. We show that like the 3’ adaptor ligation step, the 5’ adaptor ligation is also biased, not because of primary sequence, but instead due to secondary structures of the two ligation substrates. We find that multiple secondary structural factors influence final representation in HTS results. Our results provide insight about the nature of ligation bias and allowed us to design adaptors that reduce ligation bias and produce HTS results that more accurately reflect the actual concentrations of miRNAs in the defined starting material.

BackgroundMobile group II introns insert site-specifically into DNA target sites by a mechanism termed retrohoming in which the excised intron RNA reverse splices into a DNA strand and is reverse transcribed by the intron-encoded protein. Retrohoming is mediated by a ribonucleoprotein particle that contains the intron-encoded protein and excised intron RNA, with target specificity determined largely by base pairing of the intron RNA to the DNA target sequence. This feature enabled the development of mobile group II introns into bacterial gene targeting vectors (“targetrons”) with programmable target specificity. Thus far, however, efficient group II intron-based gene targeting reactions have not been demonstrated in eukaryotes.Methodology/Principal FindingsBy using a plasmid-based Xenopus laevis oocyte microinjection assay, we show that group II intron RNPs can integrate efficiently into target DNAs in a eukaryotic nucleus, but the reaction is limited by low Mg2+ concentrations. By supplying additional Mg2+, site-specific integration occurs in up to 38% of plasmid target sites. The integration products isolated from X. laevis nuclei are sensitive to restriction enzymes specific for double-stranded DNA, indicating second-strand synthesis via host enzymes. We also show that group II intron RNPs containing either lariat or linear intron RNA can introduce a double-strand break into a plasmid target site, thereby stimulating homologous recombination with a co-transformed DNA fragment at frequencies up to 4.8% of target sites. Chromatinization of the target DNA inhibits both types of targeting reactions, presumably by impeding RNP access. However, by using similar RNP microinjection methods, we show efficient Mg2+-dependent group II intron integration into plasmid target sites in zebrafish (Danio rerio) embryos and into plasmid and chromosomal target sites in Drosophila melanogster embryos, indicating that DNA replication can mitigate effects of chromatinization.Conclusions/SignificanceOur results provide an experimental foundation for the development of group II intron-based gene targeting methods for higher organisms.

We find that group II intron EcI5, a subclass CL/IIB1 intron from an Escherichia coli virulence plasmid, is highly active in retrohoming in E. coli. Both full-length EcI5 and an EcI5-DORF intron with the intron-encoded protein expressed separately from the same donor plasmid retrohome into a recipient plasmid target site at substantially higher frequencies than do similarly configured Lactococcus lactis Ll.LtrB introns. A comprehensive view of DNA target site recognition by EcI5 was obtained from selection experiments with donor and recipient plasmid libraries in which different recognition elements were randomized. These experiments suggest that EcI5, like other mobile group II introns, recognizes DNA target sequences by using both the intron-encoded protein and base-pairing of the intron RNA, with the latter involving EBS1, EBS2, and EBS3 sequences characteristic of class IIB introns. The intron-encoded protein appears to recognize a small number of bases flanking those recognized by the intron RNA, but their identity is different than in previously characterized group II introns. A computer algorithm based on the empirically determined DNA recognition rules enabled retargeting of EcI5 to integrate specifically at 10 different sites in the chromosomal lacZ gene at frequencies up to 98% without selection. Our findings provide insight into modes of DNA target site recognition used by mobile group II introns. More generally, they show how the diversity of mobile group II introns can be exploited to provide a large variety of different target specificities and potentially other useful properties for gene targeting.