Mycoreovirus genome rearrangements associated with RNA silencing deficiency

Mycoreovirus 1 (MyRV1) has 11 double-stranded RNA genome segments (S1 to S11) and confers hypovirulence to the chestnut blight fungus, Cryphonectria parasitica. MyRV1 genome rearrangements are frequently generated by a multifunctional protein, p29, encoded by a positive-strand RNA virus, Cryphonectria hypovirus 1. One of its functional roles is RNA silencing suppression. Here, we explored a possible link between MyRV1 genome rearrangements and the host RNA silencing pathway using wild-type (WT) and mutant strains of both MyRV1 and the host fungus. Host strains included deletion mutants of RNA silencing components such as dicer-like (dcl) and argonaute-like (agl) genes, while virus strains included an S4 internal deletion mutant MyRV1/S4ss. Consequently, intragenic rearrangements with nearly complete duplication of the three largest segments, i.e. S1, S2 and S3, were observed even more frequently in the RNA silencing-deficient strains Δdcl2 and Δagl2 infected with MyRV1/S4ss, but not with any other viral/host strain combinations. An interesting difference was noted between genome rearrangement events in the two host strains, i.e. generation of the rearrangement required prolonged culture for Δagl2 in comparison with Δdcl2. These results suggest a role for RNA silencing that suppresses genome rearrangements of a dsRNA virus.     

An antagonistic function for Arabidopsis DCL2 in development and a new function for DCL4 in generating viral siRNAs

Plants contain more DICER-LIKE (DCL) enzymes and double-stranded RNA binding (DRB) proteins than other eukaryotes, resulting in increased small RNA network complexities. Analyses of single, double, triple and quadruple dcl mutants exposed DCL1 as a sophisticated enzyme capable of producing both microRNAs (miRNAs) and siRNAs, unlike the three other DCLs, which only produce siRNAs. Depletion of siRNA-specific DCLs results in unbalanced small RNA levels, indicating a redeployment of DCL/DRB complexes. In particular, DCL2 antagonizes the production of miRNAs and siRNAs by DCL1 in certain circumstances and affects development deleteriously in dcl1 dcl4 and dcl1 dcl3 dcl4 mutant plants, whereas dcl1 dcl2 dcl3 dcl4 quadruple mutant plants are viable. We also show that viral siRNAs are produced by DCL4, and that DCL2 can substitute for DCL4 when this latter activity is reduced or inhibited by viruses, pointing to the competitiveness of DCL2. Given the complexity of the small RNA repertoire in plants, the implication of each DCL, in particular DCL2, in the production of small RNAs that have no known function will constitute one of the next challenges.     

Respective contributions of Arabidopsis DCL2 and DCL4 to RNA silencing

Dicer proteins are central to the different mechanisms involving RNA interference. Plants have evolved multiple DICER‐LIKE (DCL) copies, thus enabling functional diversification. In Arabidopsis, DCL2 and DCL4 process double‐stranded RNA into 22 and 21 nucleotide small interfering (si)RNAs, respectively, and have overlapping functions with regards to virus and transgene silencing. Nonetheless, some studies have reported that dcl2 or dcl4 single mutations are sometimes sufficient to hinder silencing. To better dissect the role of DCL2 and DCL4, we analyzed silencing kinetics and efficiencies using different transgenic systems in single and double mutant backgrounds. The results indicate that DCL2 stimulates transitivity and secondary siRNA production through DCL4 while being sufficient for silencing on its own. Notably, silencing of 35S‐driven transgenes functions more efficiently in dcl4 mutants, indicating that DCL4 mostly obscures DCL2 in wild‐type plants. Nonetheless, in a dcl4 mutant compromised in phloem‐originating silencing, ectopically expressed DCL2 allows restoration of silencing, suggesting that DCL2 is not, or poorly, expressed in phloem. Remarkably, this ectopic DCL2 contribution to phloem‐originating silencing is dependent on the activity of RNA‐DEPENDENT RNA POLYMERASE6. These results indicate that, despite differences in the silencing activity of their small RNA products, DCL2 and DCL4 mostly act redundantly yet hierarchically when present simultaneously.     

The interaction between DCL1 and HYL1 is important for efficient and precise processing of pri-miRNA in plant microRNA biogenesis

It has been reported that some double-stranded RNA (dsRNA) binding proteins interact with small RNA biogenesis-related RNase III enzymes. However, their biological significance is poorly understood. Here we examine the relationship between the Arabidopsis microRNA- (miRNA) producing enzyme DCL1 and the dsRNA binding protein HYL1. In the hyl1-2 mutant, the processing steps of miR163 biogenesis were partially impaired; increased accumulation of pri-miR163 and reduced accumulation of short pre-miR163 and mature miR163 as well as misplaced cleavages in the stem structure of pri-miR163 were detected. These misplaced cleavages were similar to those previously observed in the dcl1-9 mutant, in which the second double-stranded RNA binding domain of the protein was disrupted. An immunoprecipitation assay using Agrobacterium-mediated transient expression in Nicotiana benthamiana showed that HYL1 was able to form a complex with wild-type DCL1 protein, but not with the dcl1-9 mutant protein. We also examined miR164b and miR166a biogenesis in hyl1-2 and dcl1-9. Increased accumulation of pri-miRNAs and reduced accumulation of pre-miRNAs and mature miRNAs were detected. Misplaced cleavage on pri-miR164b was observed only in dcl1-9 but not in hyl1-2, whereas not on pri-miR166a in either mutant. These results indicate that HYL1 has a function in assisting efficient and precise cleavage of pri-miRNA through interaction with DCL1.     

A Turnip crinkle virus Isolate Lacking the CP Counter‐Defence Protein Gene Providing Protection Against the Wild‐Type Strain is Associated with Highly Localized RNA Silencing

Cross‐protection has been used successfully and commercially to control a range of virus diseases for which the selection of suitable mild strains of plant viruses is necessary. Turnip crinkle virus (TCV) is highly pathogenic on Arabidopsis plants and its silencing suppressor‐defective mutant, TCVΔCP, can induce highly localized RNA silencing which is differs from that of other protective strains. We found that TCVΔCP provides some protection against wild‐type TCV but lacks complete protection, and the relative locations of the protective virus and challenge virus affect the degree of cross‐protection. However, similar cross‐protection afforded by TCVΔCP is not observed in Nicotiana benthamiana plants. As expected, TCVΔCP pre‐infected Arabidopsis plants fail to protect against infection with the unrelated Cucumber mosaic virus, strain Fhy. It appears that cross‐protection afforded by TCVΔCP requires that the challenge virus be very similar in sequence, which is a characteristic of RNA silencing. In order to investigate whether the protection is associated with the highly localized RNA silencing, mutant plants involved in key silencing pathway genes of RNA silencing machinery, including dcl2, dcl4 and triple dcl2/dcl3/dcl4 mutants were used. The results demonstrate that cross‐protection afforded by TCVΔCP is dependent on host RNA silencing, and both DCL2 and DCL4 play important roles in this process.     

Mutations in the antiviral RNAi defense pathway modify Brome mosaic virus RNA recombinant profiles

RNA interference (RNAi) mechanism targets viral RNA for degradation. To test whether RNAi gene products contributed to viral RNA recombination, a series of Arabidopsis thaliana RNAi-defective mutants were infected with Brome mosaic virus (BMV) RNAs that have been engineered to support crossovers within the RNA3 segment. Single-cross RNA3-RNA1, RNA3-RNA2, and RNA3-RNA3 recombinants accumulated in both the wild-type (wt) and all knock-out lines at comparable frequencies. However, a reduced accumulation of novel 3' mosaic RNA3 recombinants was observed in ago1, dcl2, dcl4, and rdr6 lines but not in wt Col-0 or the dcl3 line. A BMV replicase mutant accumulated a low level of RNA3-RNA1 single-cross recombinants in Col-0 plants while, in a dcl2 dcl4 double mutant, the formation of both RNA3-RNA1 and mosaic recombinants was at a low level. A control infection in the cpr5-2 mutant, a more susceptible BMV Arabidopsis host, generated similar-to-Col-0 profiles of both single-cross and mosaic recombinants, indicating that recombinant profiles were, to some extent, independent of a viral replication rate. Also, the relative growth experiments revealed similar selection pressure for recombinants among the host lines. Thus, the altered recombinant RNA profiles have originated at the level of recombinant formation rather than because of altered selection. In conclusion, the viral replicase and the host RNAi gene products contribute in distinct ways to BMV RNA recombination. Our studies reveal that the antiviral RNAi mechanisms are utilized by plant RNA viruses to increase their variability, reminiscent of phenomena previously demonstrated in fungi.     

Rapid identification of cell-specific, internalizing RNA aptamers with bioinformatics analyses of a cell-based aptamer selection

BackgroundThe broad applicability of RNA aptamers as cell-specific delivery tools for therapeutic reagents depends on the ability to identify aptamer sequences that selectively access the cytoplasm of distinct cell types. Towards this end, we have developed a novel approach that combines a cell-based selection method (cell-internalization SELEX) with high-throughput sequencing (HTS) and bioinformatics analyses to rapidly identify cell-specific, internalization-competent RNA aptamers.

Methodology/Principal Findings

We demonstrate the utility of this approach by enriching for RNA aptamers capable of selective internalization into vascular smooth muscle cells (VSMCs). Several rounds of positive (VSMCs) and negative (endothelial cells; ECs) selection were performed to enrich for aptamer sequences that preferentially internalize into VSMCs. To identify candidate RNA aptamer sequences, HTS data from each round of selection were analyzed using bioinformatics methods: (1) metrics of selection enrichment; and (2) pairwise comparisons of sequence and structural similarity, termed edit and tree distance, respectively. Correlation analyses of experimentally validated aptamers or rounds revealed that the best cell-specific, internalizing aptamers are enriched as a result of the negative selection step performed against ECs.

Conclusions and Significance

We describe a novel approach that combines cell-internalization SELEX with HTS and bioinformatics analysis to identify cell-specific, cell-internalizing RNA aptamers. Our data highlight the importance of performing a pre-clear step against a non-target cell in order to select for cell-specific aptamers. We expect the extended use of this approach to enable the identification of aptamers to a multitude of different cell types, thereby facilitating the broad development of targeted cell therapies.     

Controlled RISC loading efficiency of miR168 defined by miRNA duplex structure adjusts ARGONAUTE1 homeostasis

Micro RNAs (miRNAs) are processed from precursor RNA molecules with precisely defined secondary stem-loop structures. ARGONAUTE1 (AGO1) is the main executor component of miRNA pathway and its expression is controlled via the auto-regulatory feedback loop activity of miR168 in plants. Previously we have shown that AGO1 loading of miR168 is strongly restricted leading to abundant cytoplasmic accumulation of AGO-unbound miR168. Here, we report, that intrinsic RNA secondary structure of MIR168a precursor not only defines the processing of miR168, but also precisely adjusts AGO1 loading efficiency determining the biologically active subset of miR168 pool. Our results show, that modification of miRNA duplex structure of MIR168a precursor fragment or expression from artificial precursors can alter the finely adjusted loading efficiency of miR168. In dcl1-9 mutant where, except for miR168, production of most miRNAs is severely reduced this mechanism ensures the elimination of unloaded AGO1 proteins via enhanced AGO1 loading of miR168. Based on this data, we propose a new competitive loading mechanism model for miR168 action: the miR168 surplus functions as a molecular buffer for controlled AGO1 loading continuously adjusting the amount of AGO1 protein in accordance with the changing size of the cellular miRNA pool.     

Recovery of dicer-like 1-late flowering phenotype by miR172 expressed by the noncanonical DCL4-dependent biogenesis pathway

MicroRNAs (miRNAs) act as down-regulators of gene expression, and play a dominant role in eukaryote development. In Arabidopsis thaliana, DICER-LIKE 1 (DCL1) is the main processor in miRNA biogenesis, and dcl1 mutants show various developmental defects at the early stage of embryogenesis or at gamete formation. However, miRNAs responsible for the respective developmental stages of the dcl1 defects have not been identified. Here, we developed a DCL1-independent miRNA expression system using the unique DCL4-dependent miRNA, miR839. By replacing the mature sequence in the miR839 precursor sequence with that of miR172, one of the most widely conserved miRNAs in angiosperms, we succeeded in expressing miR172 from a chimeric miR839 precursor in dcl1-7 plants and observed the repression of miR172 target gene expression. In parallel, the DCL4-dependent miR172 expression rescued the late flowering phenotype of dcl1-7 by acceleration of flowering. We established the DCL1-independent miRNA expression system, and revealed that the reduction of miR172 expression is responsible for the dcl1-7 late flowering phenotype.     

A novel non-radioactive primase–pyrophosphatase activity assay and its application to the discovery of inhibitors of Mycobacterium tuberculosis primase DnaG

Bacterial DNA primase DnaG synthesizes RNA primers required for chromosomal DNA replication. Biochemical assays measuring primase activity have been limited to monitoring formation of radioactively labelled primers because of the intrinsically low catalytic efficiency of DnaG. Furthermore, DnaG is prone to aggregation and proteolytic degradation. These factors have impeded discovery of DnaG inhibitors by high-throughput screening (HTS). In this study, we expressed and purified the previously uncharacterized primase DnaG from Mycobacterium tuberculosis (Mtb DnaG). By coupling the activity of Mtb DnaG to that of another essential enzyme, inorganic pyrophosphatase from M. tuberculosis (Mtb PPiase), we developed the first non-radioactive primase–pyrophosphatase assay. An extensive optimization of the assay enabled its efficient use in HTS (Z′ = 0.7 in the 384-well format). HTS of 2560 small molecules to search for inhibitory compounds yielded several hits, including suramin, doxorubicin and ellagic acid. We demonstrate that these three compounds inhibit Mtb DnaG. Both suramin and doxorubicin are potent (low-µM) DNA- and nucleotide triphosphate-competitive priming inhibitors that interact with more than one site on Mtb DnaG. This novel assay should be applicable to other primases and inefficient DNA/RNA polymerases, facilitating their characterization and inhibitor discovery.     

The Privileged Chemical Space Predictor (pcsp): A computer program that identifies privileged chemical space from screens of modularly assembled chemical libraries

Modularly assembled combinatorial libraries are often used to identify ligands that bind to and modulate the function of a protein or a nucleic acid. Much of the data from screening these compounds, however, is not efficiently utilized to define structure–activity relationships (SAR). If SAR data are accurately constructed, it can enable the design of more potent binders. Herein, we describe a computer program called Privileged Chemical Space Predictor (pcsp) that statistically determines SAR from high-throughput screening (HTS) data and then identifies features in small molecules that predispose them for binding a target. Features are scored for statistical significance and can be utilized to design improved second generation compounds or more target-focused libraries. The program’s utility is demonstrated through analysis of a modularly assembled peptoid library that previously was screened for binding to and inhibiting a group I intron RNA from the fungal pathogen Candida albicans.     

Differential sensitivity of Arabidopsis siRNA biogenesis mutants to genotoxic stress

Plant response to stress has been linked to different RNA-silencing processes and epigenetic mechanisms. Our recent results showed that Arabidopsis thaliana Dicer-like (DCL) mutants were impaired in transgenerational changes, recombination frequency and stress tolerance. We also found that transgenerational changes were dependent on changes in DNA methylation. Here, we hypothesized that plants deficient in the production of small RNAs would show an impaired abiotic stress response. To test this, we exposed A. thaliana dcl2, dcl3, dcl4, dcl2 dcl3 (d2d3), dcl2 dcl4 (d2d4), dcl2 dcl3 dcl4 (d2d3d4), nrpd1a, rdr2 and rdr6 mutants to methyl methane sulfonate (MMS). We found dcl4 and rdr6 to be more sensitive and dcl2, dcl3, d2d3 and rdr2 plants more resistant to MMS, as shown by fresh weight, root length and survival rate. The in vitro repair assay showed the lower ability of dcl2 and dcl3 to repair UV-damaged DNA. To summarize, we found that whereas mutants impaired in transactivating siRNA biogenesis were more sensitive to MMS, mutants impaired in natural antisense siRNA and heterochromatic siRNA biogeneses were more tolerant. Our data suggest that plant response to MMS is in part regulated through biogenesis of various siRNAs.