The evolution of RNAi as a defence against viruses and transposable elements - PubMed (original) (raw)

Review

The evolution of RNAi as a defence against viruses and transposable elements

Darren J Obbard et al. Philos Trans R Soc Lond B Biol Sci. 2009.

Abstract

RNA interference (RNAi) is an important defence against viruses and transposable elements (TEs). RNAi not only protects against viruses by degrading viral RNA, but hosts and viruses can also use RNAi to manipulate each other's gene expression, and hosts can encode microRNAs that target viral sequences. In response, viruses have evolved a myriad of adaptations to suppress and evade RNAi. RNAi can also protect cells against TEs, both by degrading TE transcripts and by preventing TE expression through heterochromatin formation. The aim of our review is to summarize and evaluate the current data on the evolution of these RNAi defence mechanisms. To this end, we also extend a previous analysis of the evolution of genes of the RNAi pathways. Strikingly, we find that antiviral RNAi genes, anti-TE RNAi genes and viral suppressors of RNAi all evolve rapidly, suggestive of an evolutionary arms race between hosts and parasites. Over longer time scales, key RNAi genes are repeatedly duplicated or lost across the metazoan phylogeny, with important implications for RNAi as an immune defence.

PubMed Disclaimer

Figures

Figure 1

Figure 1

RNAi pathways in Drosophila. (a) The antiviral siRNA pathway. Dicer-2 cuts dsRNA into siRNAs, which are loaded into an Argonaute-containing RISC that targets RNA for degradation. (b) The miRNA pathway. Primary miRNAs are transcribed from the genome, processed by Drosha and Dicer-1 into mature miRNAs. These are loaded into the Argonaute-containing effector complex (RISC), which binds mRNAs and recruits additional factors to inhibit translation. (c) The ‘ping-pong’ model of TE silencing in the Drosophila germ line. Aubergine and Argonaute-3 (Piwi family Argonautes) alternately cleave sense (red) and antisense (blue) transcripts from TEs, guided by piRNAs generated in the other half of the cycle. Cleavage both inactivates the transcript and generates the 5′ end of a new piRNA. The new piRNA-precursor is bound by the partner Piwi family member and the 3′ end degraded and then modified by the addition of a methyl group (Me). The nuclear localization of Piwi suggests that it might mediate heterochromatin assembly (Klattenhoff & Theurkauf 2008). It is unknown how this occurs, but one possibility is that the active Piwi complex binds nascent TE transcripts to recruit heterochromatin factors (Grewal & Elgin 2007). Recently, a fourth pathway has been identified, which targets TE transcripts in both the soma and the germ line, using Dcr-2 and Ago-2 from the antiviral pathway, but Loqs from the miRNA pathway (Chung et al. 2008; Czech et al. 2008).

Figure 2

Figure 2

The evolutionary rate of VSRs in positive sense ssRNA plant viruses. The rate of protein evolution (_K_A, non-synonymous sites) relative to the rate of neutral evolution (_K_S, synonymous sites) for genes taken from 17 closely related pairs of ssRNA plant viruses. For viral suppressors of RNAi (VSRs) mean _K_A/_K_S=0.29 (_n_=20), and for other genes mean _K_A/_K_S=0.14 (_n_=72). In 14 of the 17 viruses, the rate of VSR evolution was higher than the average rate of the other genes (p<0.01, sign test). Where possible, genome pairs comprised two isolates of the same viral taxon (see the electronic supplementary material for accession details).

Figure 3

Figure 3

Evolutionary rate of RNAi-related genes in Drosophila. The rate of protein evolution (_K_A, non-synonymous sites) relative to the rate of neutral evolution (_K_S, synonymous sites) between D. melanogaster and D. simulans, plotted against aligned gene length. All 10 581 orthologous genes for which sequences are available are shown in grey, and 22 RNAi-pathway genes are shown as filled circles, squares and triangles. The smoothed median (red line) and the 75th, 85th and 95th percentiles (shaded background) are plotted for a sliding window (average 380 genes wide). RNAi-related genes evolve significantly more rapidly than the genome average (Wilcoxon signed-rank test, p<0.01). Six are in the fastest 5 per cent of genes (R2D2, Maelstrom, Krimper, Ago-2, Aubergine and Dcr-2), and a further 12 are in the upper 50 per cent. The piRNA pathway (purple circles) and viRNA pathway (red triangles) each evolve more rapidly than other genes (p<0.001 and p<0.01, respectively) but the miRNA pathway (green squares) does not (_p_=0.34). _K_A/KS was calculated by the method of Li et al. (1985). This analysis partially duplicates that of Obbard et al. (2006), but uses newer genome releases and an alternative estimator of _K_A/_K_S.

Figure 4

Figure 4

The rate of adaptive evolution in Drosophila RNAi genes. The estimated number of adaptive substitutions per codon (as opposed to neutral substitutions fixed by random genetic drift) for RNAi-pathway genes was estimated using the method of Smith & Eyre-Walker (2002). The analysis uses the D. simulans polymorphism data from Obbard et al. (; Ago-1, Ago-2, Dcr-1, Dcr-2, Loqs, R2D2) and Begun et al. (; all other genes), together with the fixed differences between these datasets and the D. melanogaster genome. Negative estimates (Ago-3 and Krimper) arise due to the relatively large number of amino acid polymorphisms in these genes. Genes with individually significant McDonald–Kreitman tests (McDonald & Kreitman 1991) are indicated by asterisks (*p<0.05, **p<0.01, ***p<0.001, no correction for multiple tests). Note that the differences in sample size mean that the power of the test varies greatly between genes (details available in the electronic supplementary material).

Figure 5

Figure 5

Evolutionary relationships of metazoans and their complements of selected RNAi genes. The presence of an RdRp and the number (where known) of SID and Dicer family members (Weinstock et al. 2006; Gordon & Waterhouse 2007) are indicated in table 1. The metazoan phylogeny is compiled from Bourlat et al. (2006), Delsuc et al. (2006) and Simionato et al. (2007). See figure 1 in the electronic supplementary material for sources.

Similar articles

Cited by

References

    1. Ahlquist P. RNA-dependent RNA polymerases, viruses, and RNA silencing. Science. 2002;296:1270–1273. doi:10.1126/science.1069132 - DOI - PubMed
    1. Alcami A., Koszinowski U.H. Viral mechanisms of immune evasion. Trends Microbiol. 2000;8:410–418. doi:10.1016/S0966-842X(00)01830-8 - DOI - PMC - PubMed
    1. Allen T.M., O'Connor D.H., Jing P.C., Dzuris J.L., Mothe B.R., Vogel T.U. Tat-specific cytotoxic T lymphocytes select for SIV escape variants during resolution of primary viraemia. Nature. 2000;407:386–390. doi:10.1038/35036559 - DOI - PubMed
    1. Anandalakshmi R., Pruss G.J., Ge X., Marathe R., Mallory A.C., Smith T.H. A viral suppressor of gene silencing in plants. Proc. Natl Acad. Sci. USA. 1998;95:13 079–13 084. doi:10.1073/pnas.95.22.13079 - DOI - PMC - PubMed
    1. Aravin A.A., Klenov M.S., Vagin V.V., Bantignies F., Cavalli G., Gvozdev V.A. Dissection of a natural RNA silencing process in the Drosophila melanogaster germ line. Mol. Cell. Biol. 2004;24:6742–6750. doi:10.1128/MCB.24.15.6742-6750.2004 - DOI - PMC - PubMed

Publication types

MeSH terms

Substances

LinkOut - more resources