Qingfa Wu

Virus discovery by deep sequencing and assembly of virus-derived small silencing RNAs.

In response to infection, invertebrates process replicating viral RNA genomes into siRNAs of discrete sizes to guide virus clearance by RNA interference. Here, we show that viral siRNAs sequenced from fruit fly, mosquito, and nematode cells were all overlapping in sequence, suggesting a possibility of using siRNAs for viral genome assembly and virus discovery. To test this idea, we examined contigs assembled from published small RNA libraries and discovered five previously undescribed viruses from cultured Drosophila cells and adult mosquitoes, including three with a positive-strand RNA genome and two with a dsRNA genome. Notably, four of the identified viruses exhibited only low sequence similarities to known viruses, such that none could be assigned into an existing virus genus. We also report detection of virus-derived PIWI-interacting RNAs (piRNAs) in Drosophila melanogaster that have not been previously described in any other host species and demonstrate viral genome assembly from viral piRNAs in the absence of viral siRNAs. Thus, this study provides a powerful culture-independent approach for virus discovery in invertebrates by assembling viral genomes directly from host immune response products without prior virus enrichment or amplification. We propose that invertebrate viruses discovered by this approach may include previously undescribed human and vertebrate viral pathogens that are transmitted by arthropod vectors.

RNAi-mediated viral immunity requires amplification of virus-derived siRNAs in Arabidopsis thaliana

In diverse eukaryotic organisms, Dicer-processed, virus-derived small interfering RNAs direct antiviral immunity by RNA silencing or RNA interference. Here we show that in addition to core dicing and slicing components of RNAi, the RNAi-mediated viral immunity in Arabidopsis thaliana requires host RNA-directed RNA polymerase (RDR) 1 or RDR6 to produce viral secondary siRNAs following viral RNA replication-triggered biogenesis of primary siRNAs. We found that the two antiviral RDRs exhibited specificity in targeting the tripartite positive-strand RNA genome of cucumber mosaic virus (CMV). RDR1 preferentially amplified the 5′-terminal siRNAs of each of the three viral genomic RNAs, whereas an increased production of siRNAs targeting the 3′ half of RNA3 detected in rdr1 mutant plants appeared to be RDR6-dependent. However, siRNAs derived from a single-stranded 336-nucleotide satellite RNA of CMV were not amplified by either antiviral RDR, suggesting avoidance of the potent RDR-dependent silencing as a strategy for the molecular parasite of CMV to achieve preferential replication. Our work thus identifies a distinct mechanism for the amplification of immunity effectors, which together with the requirement for the biogenesis of endogenous siRNAs, may play a role in the emergence and expansion of eukaryotic RDRs.

The 21-Nucleotide, but Not 22-Nucleotide, Viral Secondary Small Interfering RNAs Direct Potent Antiviral Defense by Two Cooperative Argonautes in Arabidopsis thaliana

This work identifies cooperative action of ARGONAUTE1 and ARGONAUTE2 in virus resistance conferred by 21-nucleotide virus-derived small interfering RNAs (siRNAs). It also reveals that 22-nucleotide viral siRNAs do not guide efficient antiviral defense, demonstrating a qualitative difference between 21- and 22-nucleotide classes of siRNAs in RNA silencing. Arabidopsis thaliana defense against distinct positive-strand RNA viruses requires production of virus-derived secondary small interfering RNAs (siRNAs) by multiple RNA-dependent RNA polymerases. However, little is known about the biogenesis pathway and effector mechanism of viral secondary siRNAs. Here, we describe a mutant of Cucumber mosaic virus (CMV-Δ2b) that is silenced predominantly by the RNA-DEPENDENT RNA POLYMERASE6 (RDR6)-dependent viral secondary siRNA pathway. We show that production of the viral secondary siRNAs targeting CMV-Δ2b requires SUPPRESSOR OF GENE SILENCING3 and DICER-LIKE4 (DCL4) in addition to RDR6. Examination of 25 single, double, and triple mutants impaired in nine ARGONAUTE (AGO) genes combined with coimmunoprecipitation and deep sequencing identifies an essential function for AGO1 and AGO2 in defense against CMV-Δ2b, which act downstream the biogenesis of viral secondary siRNAs in a nonredundant and cooperative manner. Our findings also illustrate that dicing of the viral RNA precursors of primary and secondary siRNA is insufficient to confer virus resistance. Notably, although DCL2 is able to produce abundant viral secondary siRNAs in the absence of DCL4, the resultant 22-nucleotide viral siRNAs alone do not guide efficient silencing of CMV-Δ2b. Possible mechanisms for the observed qualitative difference in RNA silencing between 21- and 22-nucleotide secondary siRNAs are discussed.

The birth and death of microRNA genes in Drosophila

MicroRNAs (miRNAs) are small, endogenously expressed RNAs that regulate mRNAs post-transcriptionally. The class of miRNA genes, like other gene classes, should experience birth, death and persistence of its members. We carried out deep sequencing of miRNAs from three species of Drosophila, and obtained 107,000 sequences that map to no fewer than 300 loci that were not previously known. We observe a large class of miRNA genes that are evolutionarily young, with a rate of birth of 12 new genes per million years (Myr). Most of these new miRNAs originated from non-miRNA sequences. Among the new genes, we estimate that 96% disappeared quickly in the course of evolution; only 4% of new miRNA genes were retained by natural selection. Furthermore, only 60% of these retained genes became integrated into the transcriptome in the long run (60 Myr). This small fraction (2.5%) of surviving miRNAs may later on become moderately or highly expressed. Our results suggest that there is a high birth rate of new miRNA genes, accompanied by a comparably high death rate. The estimated net gain of long-lived miRNA genes, which is not strongly affected by either the depth or the breadth (number of tissues) of sequencing, is 0.3 genes per Myr in Drosophila.

Viral Suppressors of RNA-Based Viral Immunity: Host Targets

RNA-based viral immunity (RVI) operates in fungi, plants, and invertebrates to specifically destroy viral RNAs using the cellular RNA interference machinery. Discovery of diverse viral proteins as suppressors of RNA silencing provides strong validation for the effectiveness of RVI. Here we review recent studies that have revealed new mechanistic insights into plant and insect viral suppressors of RVI or suggested a role for RNA silencing suppression during mammalian viral infection.

Homology-independent discovery of replicating pathogenic circular RNAs by deep sequencing and a new computational algorithm

A common challenge in pathogen discovery by deep sequencing approaches is to recognize viral or subviral pathogens in samples of diseased tissue that share no significant homology with a known pathogen. Here we report a homology-independent approach for discovering viroids, a distinct class of free circular RNA subviral pathogens that encode no protein and are known to infect plants only. Our approach involves analyzing the sequences of the total small RNAs of the infected plants obtained by deep sequencing with a unique computational algorithm, progressive filtering of overlapping small RNAs (PFOR). Viroid infection triggers production of viroid-derived overlapping siRNAs that cover the entire genome with high densities. PFOR retains viroid-specific siRNAs for genome assembly by progressively eliminating nonoverlapping small RNAs and those that overlap but cannot be assembled into a direct repeat RNA, which is synthesized from circular or multimeric repeated-sequence templates during viroid replication. We show that viroids from the two known families are readily identified and their full-length sequences assembled by PFOR from small RNAs sequenced from infected plants. PFOR analysis of a grapevine library further identified a viroid-like circular RNA 375 nt long that shared no significant sequence homology with known molecules and encoded active hammerhead ribozymes in RNAs of both plus and minus polarities, which presumably self-cleave to release monomer from multimeric replicative intermediates. A potential application of the homology-independent approach for viroid discovery in plant and animal species where RNA replication triggers the biogenesis of siRNAs is discussed.

Poly A- Transcripts Expressed in HeLa Cells

Background Transcripts expressed in eukaryotes are classified as poly A+ transcripts or poly A- transcripts based on the presence or absence of the 3′ poly A tail. Most transcripts identified so far are poly A+ transcripts, whereas the poly A- transcripts remain largely unknown. Methodology/Principal Findings We developed the TRD (Total RNA Detection) system for transcript identification. The system detects the transcripts through the following steps: 1) depleting the abundant ribosomal and small-size transcripts; 2) synthesizing cDNA without regard to the status of the 3′ poly A tail; 3) applying the 454 sequencing technology for massive 3′ EST collection from the cDNA; and 4) determining the genome origins of the detected transcripts by mapping the sequences to the human genome reference sequences. Using this system, we characterized the cytoplasmic transcripts from HeLa cells. Of the 13,467 distinct 3′ ESTs analyzed, 24% are poly A-, 36% are poly A+, and 40% are bimorphic with poly A+ features but without the 3′ poly A tail. Most of the poly A- 3′ ESTs do not match known transcript sequences; they have a similar distribution pattern in the genome as the poly A+ and bimorphic 3′ ESTs, and their mapped intergenic regions are evolutionarily conserved. Experiments confirmed the authenticity of the detected poly A- transcripts. Conclusion/Significance Our study provides the first large-scale sequence evidence for the presence of poly A- transcripts in eukaryotes. The abundance of the poly A- transcripts highlights the need for comprehensive identification of these transcripts for decoding the transcriptome, annotating the genome and studying biological relevance of the poly A- transcripts.

Identification of Viruses and Viroids by Next-Generation Sequencing and Homology-Dependent and Homology-Independent Algorithms

A fast, accurate, and full indexing of viruses and viroids in a sample for the inspection and quarantine services and disease management is desirable but was unrealistic until recently. This article reviews the rapid and exciting recent progress in the use of next-generation sequencing (NGS) technologies for the identification of viruses and viroids in plants. A total of four viroids/viroid-like RNAs and 49 new plant RNA and DNA viruses from 18 known or unassigned virus families have been identified from plants since 2009. A comparison of enrichment strategies reveals that full indexing of RNA and DNA viruses as well as viroids in a plant sample at single-nucleotide resolution is made possible by one NGS run of total small RNAs, followed by data mining with homology-dependent and homology-independent computational algorithms. Major challenges in the application of NGS technologies to pathogen discovery are discussed.

A large quantity of novel human antisense transcripts detected by LongSAGE

MOTIVATION Taking advantage of the high sensitivity and specificity of LongSAGE tag for transcript detection and genome mapping, we analyzed the 632 813 unique human LongSAGE tags deposited in public databases to identify novel human antisense transcripts. RESULTS Our study identified 45 321 tags that match the antisense strand of 9804 known mRNA sequences, 6606 of which contain antisense ESTs and 3198 are mapped only by SAGE tags. Quantitative analysis showed that the detected antisense transcripts are present at levels lower than their counterpart sense transcripts. Experimental results confirmed the presence of antisense transcripts detected by the antisense tags. We also constructed an antisense tag database that can be used to identify the antisense SAGE tags originated from the antisense strand of known mRNA sequences included in the RefSeq database. CONCLUSIONS Our study highlights the benefits of exploring SAGE data for comprehensive identification of human antisense transcripts and demonstrates the prevalence of antisense transcripts in the human genome.

Discovery of Replicating Circular RNAs by RNA-Seq and Computational Algorithms

Replicating circular RNAs are independent plant pathogens known as viroids, or act to modulate the pathogenesis of plant and animal viruses as their satellite RNAs. The rate of discovery of these subviral pathogens was low over the past 40 years because the classical approaches are technical demanding and time-consuming. We previously described an approach for homology-independent discovery of replicating circular RNAs by analysing the total small RNA populations from samples of diseased tissues with a computational program known as progressive filtering of overlapping small RNAs (PFOR). However, PFOR written in PERL language is extremely slow and is unable to discover those subviral pathogens that do not trigger in vivo accumulation of extensively overlapping small RNAs. Moreover, PFOR is yet to identify a new viroid capable of initiating independent infection. Here we report the development of PFOR2 that adopted parallel programming in the C++ language and was 3 to 8 times faster than PFOR. A new computational program was further developed and incorporated into PFOR2 to allow the identification of circular RNAs by deep sequencing of long RNAs instead of small RNAs. PFOR2 analysis of the small RNA libraries from grapevine and apple plants led to the discovery of Grapevine latent viroid (GLVd) and Apple hammerhead viroid-like RNA (AHVd-like RNA), respectively. GLVd was proposed as a new species in the genus Apscaviroid, because it contained the typical structural elements found in this group of viroids and initiated independent infection in grapevine seedlings. AHVd-like RNA encoded a biologically active hammerhead ribozyme in both polarities, and was not specifically associated with any of the viruses found in apple plants. We propose that these computational algorithms have the potential to discover novel circular RNAs in plants, invertebrates and vertebrates regardless of whether they replicate and/or induce the in vivo accumulation of small RNAs.