Plant microRNAs and their role in defense against viruses: a bioinformatics approach
© Pérez-Quintero et al; licensee BioMed Central Ltd. 2010
Received: 24 March 2010
Accepted: 1 July 2010
Published: 1 July 2010
microRNAs (miRNAs) are non-coding short RNAs that regulate gene expression in eukaryotes by translational inhibition or cleavage of complementary mRNAs. In plants, miRNAs are known to target mostly transcription factors and are implicated in diverse aspects of plant growth and development. A role has been suggested for the miRNA pathway in antiviral defense in plants. In this work, a bioinformatics approach was taken to test whether plant miRNAs from six species could have antiviral activity by targeting the genomes of plant infecting viruses.
All plants showed a repertoire of miRNAs with potential for targeting viral genomes. The viruses were targeted by abundant and conserved miRNA families in regions coding for cylindrical inclusion proteins, capsid proteins, and nuclear inclusion body proteins. The parameters for our predicted miRNA:target pairings in the viral genomes were similar to those for validated targets in the plant genomes, indicating that our predicted pairings might behave in-vivo as natural miRNa-target pairings. Our screening was compared with negative controls comprising randomly generated miRNAs, animal miRNAs, and genomes of animal-infecting viruses. We found that plant miRNAs target plant viruses more efficiently than any other sequences, but also, miRNAs can either preferentially target plant-infecting viruses or target any virus without preference.
Our results show a strong potential for antiviral activity of plant miRNAs and suggest that the miRNA pathway may be a support mechanism to the siRNA pathway in antiviral defense.
RNA silencing is a conserved defense mechanism that plants and other eukaryotes use to protect their genomes against aberrant nucleic acids. This process uses short RNAs (20-30 nt) to recognize and manipulate complementary nucleic acids [1, 2]. At least five classes of these small regulatory RNAs have been characterized, including microRNAs (miRNAs), small interference RNAs (siRNAs), transacting siRNAs (ta-siRNAs), natural antisense siRNAs (nat-siRNAs) and, in metazoans, the Piwi-interacting RNAs [3, 4]. miRNAs and siRNAs are chemically indistinguishable and participate in partially overlapping pathways; both are derived from double-stranded RNA (dsRNA) and are then processed into 21-22 nt single stranded molecules by Dicer or a Dicer-like enzyme; later, they are incorporated into the RNA-induced silencing complex (RISC) to guide the cleavage or translational repression of the complementary strand [1, 5]. The main differences between miRNAs and siRNAs lie in their biogenesis and in their target molecules. siRNAs are generally derived from endogenous aberrant dsRNAs or from exogenous agents such as viruses, and silence the same molecule from which they originated. miRNAs, instead, originate from nuclear genes and act in trans, silencing mRNAs from other genes [6, 7].
In plants miRNAs were described first in Arabidopsis [8, 9], and later in other species. To date, there over 2200 plant miRNAs from over 30 species available at the miRBase . Most of these miRNAs target transcription factors and thus are implicated in diverse aspects of plant growth and development [11, 12].
In addition to regulate the endogenous expression of some genes, miRNAs could have a direct role in viral defense. This has been shown for various cases in animal-infecting viruses. For example miR-32 restricts the replication of the primate foamy virus type 1, miR-122 targets the hepatitis C virus and at least four miRNAs expressed in T-cells impair HIV replication [2, 13–15]. Also an important role for miRNAs in antiviral defense in humans has been suggested through bioinformatics . Likewise, animal-infecting viruses can encode miRNAs to regulate both the viral life cycle and the interaction between viruses and their hosts [17, 18].
Whereas siRNAs are known to play an important and direct role in antiviral defense in plants [19, 20], so far, there has not been proof of naturally occurring plant microRNAs with antiviral activity. It has been shown, using genetically modified viruses and plants, that complementarity between a plant miRNA and a virus genome is enough for antiviral activity. Transgenic tobacco and Arabidopsis plants displayed resistance against Cucumber mosaic virus (CMV), Turnip yellow mosaic virus (TYMV) and Turnip mosaic virus (TuMV) when expressing artificial miRNAS directed against regions in the viruses' genomes [21–23]. Also, inserting the target sequence of host plant's miRNAs in the virus genome can impair virus infectivity; however, the virus can escape rapidly of the miRNA action by mutations .
It has been suggested that virtually any endogenous small RNA could hold an intrinsic, albeit fortuitous, antiviral potential (by random complementarity) that is independent of its cellular function [15, 24–26]. Also, several sequences of 20-25 nt located within Arabidopsis intergenic regions share perfect or near perfect complementarity with a variety of plant virus genomes, but have not been validated as miRNAs yet . There are also a large number of non-conserved RNAs with unknown targets ("orphan" miRNAs) that could have an antiviral role and constitute a reservoir of defensive molecules due to their complementarity to invading viral genomes .
In this work, we present a bioinformatics approach to explore the possibility of endogenous plant miRNAs having a role in antiviral defense by targeting the genomes of plant-infecting viruses and the results are considered in the context of the evolution of plant-virus interactions.
Statistics for plant microRNAs vs plant viruses' genomes screenings
Number of miRNAs
Number of viruses (total size)
20 (149. 9 Kb)
20 (146. 1 Kb)
20 (223.7 Kb)
20 (166. 2 Kb)
16 (166. 4)
23 (188. 8 Kb)
Matching percentage (Hits/Sample size) × 100
Positive miRNAs b (%)
Positive viruses c (%)
To validate our hypotheses that plant-infecting viruses are more likely to be targeted by plant miRNAs than by other sequences and that plant miRNAs preferentially target plant-infecting viruses over other sequences, we conducted the following analyses. We created a group of negative controls to screen for miRNA targets in the following cases i) animal miRNAs vs plant virus genomes, ii) random generated miRNAs vs plant virus genomes, iii) randomized plant miRNAs vs plant virus genomes and iv) plant miRNAs vs animal virus genomes. As a positive control we screened the set of plant miRNAs against their validated target sequences.
miRANDA scoring values for each screening
Plant miRNAs vs Plant viruses
Animal miRNAs vs Plant viruses
Rand1 vs Plant viruses
Rand2 vs Plant viruses
Plant miRNAs vs Animal viruses
Plant miRNAs vs Control Targets
Free-folding Energy (kcal/mol)
The miRanda score of the positive control was significantly higher than the score of the plant miRNA vs plant viruses screening, while the miRanda score for three of the four negative controls was significantly lower. However, all miRANDA scores are above the threshold of what is considered necessary for biological activity. We should also take into account that this parameter gives a high weight to pairing in the 5'region of the miRNA which is not as crucial for plant miRNAs activity as for animal miRNAs (Table 2) .
Matching percentages of plant microRNAs vs plant viruses screenings compared to negative controls
Plant microRNAs vs Plant viruses screenings
Animal microRNAs vs Plant viruses
Rand1 vs Plant viruses
Rand2 vs Plant viruses
Plant microRNAs vs Animal viruses
A. thaliana- m vs
A. thaliana - viruses
↓ < 0.0001
↓ < 0.0001
↓ < 0.0001
↑ < 0.0001
G. max microRNAs vs
G.max - viruses
↓ < 0.0001
↓ < 0.0001
↓ < 0.0001
O. sativa microRNAs vs
↓ < 0.0001
↓ < 0.0001
↓ < 0.0001
S. bicolor microRNAs vs
S.bicolor - viruses
↓ < 0.0001
↓ < 0.0001
↓ < 0.0001
↑ < 0.0001
↓ < 0.0001
↓ < 0.0001
↓ < 0.0001
↓ < 0.0001
Z. mays microRNAs vs
Z. mays viruses
↓ < 0.0001
↓ < 0.0001
↓ < 0.0001
↑ < 0.0001
Using a bioinformatics approach we found that plant miRNAs potentially target genomic regions in plant-infecting viruses. To validate our results we carried out several positive and negative controls and these showed that the genomes of plant viruses are preferentially targeted by their host's miRNAs but were not conclusive regarding the specificity of plant miRNAs for the genomes of plant viruses. A similar trend has been found using a bioinformatics approach with animal miRNAs vs animal viruses , where the miRNA pathway has been proved to have antiviral role in Metazoans [2, 13–15]. This suggests that our predicted pairings could also have a biological function, although an experimental biological validation is necessary. It is possible that some of the viral targets found in this study are the result of purely fortuitous matches as has been suggested by various authors [15, 26, 27, 30]. Even if these pairings are the result of chance instead of selection, it is possible that given the right physiological circumstances (e.g. high expression of the miRNAs, lack of silencing suppressor in the virus) these miRNAs would efficiently silence the predicted targets. This hypothesis is supported by studies showing that artificial miRNAs can mediate antiviral defense in plants and that complementarity with the target is enough to produce resistance [21–24]. Also, plants defective in miRNA-silencing have shown to be more susceptible to some viruses .
It was reported that human miRNAs were more likely to target the genomes of human-infecting viruses over non-host's viruses . Such specificity could not be demonstrated for plant miRNAs in the present study. However, a large amount of the targets we found for plant miRNAs in the genomes of animal viruses are in non-coding regions and are therefore unlikely to impair viral activity (Figure 4). Additionally, some predicted targets of plant miRNAs were found both in plant and animal viruses (e.g. capsid genes) which may indicate a preference to target conserved regions in viruses. Finally, it is possible that the genomes of plant-infecting viruses are undergoing rapid evolution to avoid targeting by plant miRNAs, therefore giving lower matching percentages than expected. This is plausible since it has been shown that viruses can rapidly evolve to escape miRNA targeting in plants .
To identify possible plant miRNAs in the viruses genomes we used strict parameters based on experimentally valid miRNA:target pairings to ensure potential biological activity. Even considering the inherent difficulties of the computational prediction of miRNA targets, which often results large number of false positive targets , it is possible that our conservative approach has underestimated the number of candidate targets. Increasing evidence has shown that miRNA-mediated silencing in plants can occur in relaxed miRNA:target pairings, mainly leading to translational arrest instead of mRNA cleavage, although the mechanisms are not fully understood [17, 33–36]. Once the criteria for miRNA-mediated translational arrest in plants are fully understood, new approaches searching for plant miRNA targets in viral genomes may be necessary.
We found that miRNAs from deeply conserved and highly expressed families (e. g. families 156, 395, 159, 166, 160) have more potential targets in the viruses' genomes. This could suggest a way in which abundant plant miRNAs are selected to have multiplicity of functions including pathogen defense. This is supported by the fact that these families have multiple targets within the plant genomes , and some of them have been shown to be differentially expressed in response to stresses. For example, miRNAs 395 and 399 are responsive to abiotic stress (phosphorus and sulfates starvation) [37, 38], and miRNAs 156, 159 and 160 are responsive to viral infections [39–42].
By contrast, the more phylogenetically restricted families (e.g. families 495, 414, 818, 854, 1861), may be participating in more specific plant-virus interactions. Indeed, in some plants there is a large diversity of non-conserved and "young" miRNAs with still unknown targets that could be potentially employed against viral sequences [43, 44]. The lack of potential antiviral activity for some microRNA families could also be the result of them being expressed at very low levels or in a tissue or cell-specific manner, thus being less likely to play a significant role in antiviral defense.
It is also important to consider some arguments that do not support a putative function of plant miRNAs as an effective option for antiviral defense. First, most viruses encode for silencing repressors, which could directly interfere with the miRNA machinery [27, 45, 46]. Second, viral genomes evolve much faster than host miRNAs [11, 24]. Third, the miRNA signal is neither systemic nor quickly amplified . Nevertheless, using miRNAs to protect against virus might be an advantageous preemptive measure (a plant would be resistant to viruses that has never encountered before) benefitting of their ability to pair with multiple targets .
The apparent inadequacy of miRNAs as an antiviral defense mechanism may indicate that their role is not as direct as siRNAs. On one side miRNAs may simply act as a support mechanism for siRNAs. On the other side, the targets found here may be a reflection of a virus adaptation phenomenon in which they take advantage of the host miRNAs to suppress their own replication to evade immune elimination and establish in this way a persistent infection as has been suggested by Mahajan et al., . In this case the role of miRNAs would be to reach an equilibrated host-virus interaction .
Also, these results can be discussed in the context of the hypothesis proposed by Lu et al., , which states that early in plant evolution miRNAs played an important role in anti-viral defense and then novel functions evolved after the requirements of survival were satisfied [26, 33]. At this initial time, plant miRNAs may have been crucial for shaping the host ranges of several virus groups. Then, some of these "antiviral miRNAs" might have been selected to regulate endogenous genes after fortuitous matching. Both the rapid evolution of viruses and the necessity of precise gene regulation could have worked as selective pressures towards the modern miRNA pathway since the requirement for a high degree of complementarity between plant miRNAs and their targets can act as a stabilizer, preventing sequence drift even over long periods of evolutionary time . Many miRNAs might have been originated from invading viral sequences, a pathway for miRNA evolution that has been suggested previously for plants . Additionally, bioinformatics evidence suggests a transition from viral sequence to siRNA to miRNA gene in plants . Our candidate targets may be an indication of these virus-derived miRNAs, especially those found for phylogenetically restricted miRNA families with unknown genomic targets.
Our work presents initial evidence for the suspected potential of antiviral activity mediated by plant miRNAs, which is likely to have played a role in early plant evolution and in shaping host ranges for plant infecting viruses.
miRNA sequences from six plants (Arabidopsis thaliana, Glycine max, Oryza sativa, Sorghum bicolor, Vitis vinifera, and Zea mays) were downloaded from the miRBASE . These species were selected for having at least 60 available sequenced miRNAs as late as March 2009, and for being hosts of at least 10 plant-infecting viruses with fully sequenced genomes. For comparisons, miRNAs from eight Metazoan species (Coenorhabditis elegans, Drosophila melanogaster, Dario rerio, Gallus gallus, Homo sapiens, Mus musculus, Ornithorhynchus anatinus, and Pan troglodytes) were selected and 50 miRNAs sequences for each animal were also downloaded from the miRBASE .
Complete genome sequences for plant and animal-infecting viruses were obtained from Genbank . Host ranges and related information for plant viruses were consulted using the Description of Plant Virus Database, DPVWeb  and the Plant virus Database, VIDE . For animal-infecting viruses we used the International Committee on Taxonomy of Viruses (ICTV) database .
Two sets of random miRNAS were made, one using a Perl script generating random 21-nucleotide sequences , and another one by randomizing the plant miRNA sequences with the Bioedit software , doing 1000 random swap operations.
Information for verified targets of plant miRNAs were obtained from Tarbase , the Arabidopsis thaliana small RNA project , the Arabidospsis MPSS Database , the Plant microRNA Database , or primary literature [11, 59, 60]. The corresponding sequences were downloaded from Genbank .
Targets for each set of miRNAs were searched in viral genomes using a modified version of miRanda (v September 2008) . This software uses a scoring system based on the complementarities of nucleotides, similar to the Smith-Waterman algorithm. The scoring matrix used for this analysis also allows G = U 'wobble' pairs, which are important for the accurate detection of RNA:RNA duplexes. The algorithm uses folding routines from the Vienna 1.3 RNA secondary structure programming library . Although miRanda was originally designed to search for miRNA targets in animals, it is versatile enough to be modified and has been used to search for targets in viruses and plants, and has proven to be an efficient method [30, 62]. The miRanda screenings were repeated several times using randomly generated subsets of either the miRNA or the viral genome sets.
MiRanda screenings were made using different combinations of miRNAs and viral genomes. The main one was plant miRNAs against plant viruses' genomes. This was compared with four other control screenings: (i) animal miRNAs vs plant viruses, (ii) random 21 nt sequences (Rand1) vs plant viruses (iii) randomized plant miRNAs (Rand2) vs plant viruses and (iv) plant miRNAs vs animal-infecting viruses. As a positive control, the plant miRNAs were screened against 190 sequences corresponding to verified miRNA targets in the plant genomes.
The criteria to consider a sequence as a putative miRNA target were: four or fewer mismatches overall, only one or none mismatches in the 5' region of the miRNA (positions 1 to 12), no more than two consecutives mismatches in positions 13 to 21, no mismatches in positions 10 and 11. Additionally, the miRNA:target pair should have low free-energy of bonding (maximum -20 kcal/mol). These criteria are based on experimental work and have been extensively used for miRNA target prediction in various plants [23, 63, 64].
Four miRanda parameters obtained in the different screenings were used to compare and validate the predicted targets. These parameters were: a) the free-folding energy of the miRNA:target pair, which is commonly used as a measure for miRNA target prediction and indicates the stability of the miRNA:target duplex and the likeliness of correct matching and cleavage; b) the percentage identity, which indicates how many bases are complementary between the miRNAs and the target; c) the Z-score, which is based on a distribution of the shuffled alignment score; a high Z-score means that the alignment is least likely to be the result of chance; and d) the miRanda score, which weights all the others parameters and also each base pair in the alignment based on complementarity and position; it represents a measure of the number of mismatches and their distribution (mismatches in the 5' end of the target are given a higher penalization) .
The main variable used to compare the screenings was the matching percentage = [Number of candidates/(Size of the virus' genome (kb) × Number of miRNAs)] × 100, which is the percentage of the screened sample that resulted in target candidates. For statistical analysis, the Shapiro Normality test and Wilcoxon tests were performed with the software R .
To compare the targeted regions in the viral genomes, the number of hits in each region was divided by the average size in kilobases of this region in the various viruses' genomes.
We are grateful to Liliana López Klein for their support in the statistical analyses. Thanks to Ivan Acosta for fruitful discussions and for critically reading the manuscript and for editing. This research was supported partially by the DIB (Dirección de Investigaciones, sede Bogotá), Universidad Nacional de Colombia.
- Baulcombe D: RNAi silencing in plants. Nature. 2004, 431: 356-363. 10.1038/nature02874.PubMedView Article
- Obbard DJ, Gordon KHJ, Buck AH, Jiggins FM: The evolution of RNAi as a defence against viruses and transposable elements. Phil Trans R Soc B. 2009, 364: 99-115. 10.1098/rstb.2008.0168.PubMedPubMed CentralView Article
- Meins F, Si-Ammour A, Blevins T: RNA silencing systems and their relevance to plant development. Annu Rev Cell Dev Biol. 2005, 21: 297-318. 10.1146/annurev.cellbio.21.122303.114706. 4PubMedView Article
- Vaucheret H: Post-transcriptional small RNA pathways in plants: Mechanisms and regulations. Genes Dev. 2006, 20: 759-771. 10.1101/gad.1410506.PubMedView Article
- Lindbo J, Dougherty W: Plant Pathology and RNAi: A Brief History. Annu Rev Phytopathol. 2005, 43: 191-204. 10.1146/annurev.phyto.43.040204.140228.PubMedView Article
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.PubMedView Article
- Tang GL, Reinhart BJ, Bartel DP, Zamore PD: A biochemical framework for RNA silencing in plants. Genes Dev. 2003, 17: 49-63. 10.1101/gad.1048103.PubMedPubMed CentralView Article
- Llave C, Kasschau K, Rector M, Carrington J: Endogenous and Silencing-Associated Small RNAs in Plants. The Plant Cell. 2002, 14: 1605-1619. 10.1105/tpc.003210.PubMedPubMed CentralView Article
- Reinhart B, Weinstein E, Rhoades M, Bartel B, Bartel D: MicroRNAs in plants. Gene Dev. 2002, 16: 1616-1626. 10.1101/gad.1004402.PubMedPubMed CentralView Article
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, D154-D158. 38 Database
- Dugas D, Bartel B: MicroRNA regulation of gene expression in plants. Curr Opin Plant Biol. 2004, 7: 512-52. 10.1016/j.pbi.2004.07.011.PubMedView Article
- Mallory A, Vaucheret H: Functions of microRNAs and related small RNAs in plants. Nat Genet. 2006, 38 (Suppl): S31-S36. 10.1038/ng1791.PubMedView Article
- Berkhout B, Haasnoot J: The interplay between virus infection and the cellular RNA interference Machinery. FEBS Letters. 2006, 580: 2896-2902. 10.1016/j.febslet.2006.02.070.PubMedView Article
- Bushati N, Cohen S: microRNA functions. Annu Rev Cell Dev Biol. 2007, 23: 175-205. 10.1146/annurev.cellbio.23.090506.123406.PubMedView Article
- Lecellier CH, Dunoyer P, Arar K, Lehmann-Che J, Eyquem S, Himber C, Saïb A, Voinnet O: A cellular microRNA mediates antiviral defense in human cells. Science. 2005, 308 (5721): 480-1. 10.1126/science.1108784.View Article
- Watanabe Y, Kishi A, Yachie N, Kanai A, Tomita M: Computational analysis of microRNA-mediated antiviral defense in humans. FEBS Letters. 2007, 581: 4603-4610. 10.1016/j.febslet.2007.08.049.PubMedView Article
- Nair V, Zavolan M: Virus-encoded microRNAs: novel regulators of gene expression. Trends Microbiol. 2006, 14 (4): 169-175. 10.1016/j.tim.2006.02.007.PubMedView Article
- Scaria V, Hariharan M, Pillai B, Maiti S, Brahmachari S: Host-virus genome interactions: macro roles for microRNAs. Cell Microbiol. 2007, 9 (12): 2784-2794V. 10.1111/j.1462-5822.2007.01050.x.PubMedView Article
- Ding SW, Li H, Lu R, Li F, Li WX: RNA silencing: a conserved antiviral immunity of plants and animals. Virus Res. 2004, 102: 109-115. 10.1016/j.virusres.2004.01.021.PubMedView Article
- Mahmood-ur-Rahman AI, Husnain T, Riazuddin S: RNA interference: the story of gene silencing in plants and humans. Biotechnol Adv. 2008, 26 (3): 202-9. 10.1016/j.biotechadv.2007.12.002.PubMedView Article
- Niu Q, Lin S, Reyes J, Chen K, Wu H, Ye S, Chua N: Expression of artificial microRNAs in transgenic Arabidopsis thaliana confers virus resistance. Nature Biotechnol. 2006, 24 (11): 1420-1427. 10.1038/nbt1255.View Article
- Qu J, Ye J, Fang R: Artificial MicroRNA-Mediated Virus Resistance in Plants. J Virol. 2007, 81: 6690-6699. 10.1128/JVI.02457-06.PubMedPubMed CentralView Article
- Lin S, Wu H, Elena S, Chen K, Niu Q, Ye S, Chen C, Chua N: Molecular Evolution of a Viral Non-Coding Sequence under the Selective Pressure of amiRNA-Mediated Silencing. PLoS Pathog. 2009, 5 (2): e1000312-10.1371/journal.ppat.1000312.PubMedPubMed CentralView Article
- Simón-Mateo C, García J: MicroRNA-Guided Processing Impairs Plum Pox Virus Replication, but the Virus Readily Evolves To Escape This Silencing Mechanism. J Virol. 2006, 80 (5): 2429-2436. 10.1128/JVI.80.5.2429-2436.2006.PubMedPubMed CentralView Article
- Dunoyer P, Voinnet O: The complex interplay between plant viruses and host RNA-silencing pathways. Curr Opin Plant Biol. 2005, 8: 415-423. 10.1016/j.pbi.2005.05.012.PubMedView Article
- Lu Y, Gan Q, Chi X, Qin S: Roles of microRNA in plant defense and virus offense interaction. Plant Cell Rep. 2008, 27: 1571-1579. 10.1007/s00299-008-0584-z.PubMedView Article
- Llave C: MicroRNAs: more than a role in plant development?. Mol plant pathol. 2004, 5 (4): 361-366. 10.1111/j.1364-3703.2004.00227.x.PubMedView Article
- Büchen-Osmond C: Taxonomy and Classification of Viruses. Manual of Clinical Microbiology. 2003, ASM Press, Washington DC, 2: 1217-1226. 8
- Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS: MicroRNA targets in Drosophila. Genome Biol. 2003, 5 (1): R1-10.1186/gb-2003-5-1-r1.PubMedPubMed CentralView Article
- Hsu P, Lin L, Hsu S, Hsu J, Huang H: ViTa: prediction of host microRNAs targets on viruses. Nucl Acids Res. 2007, D381-D385. 10.1093/nar/gkl1009. 35 Database
- MacDiarmid R: RNA silencing in productive virus infections. Annu Rev Phytopathol. 2005, 43: 523-44. 10.1146/annurev.phyto.43.040204.140204.PubMedView Article
- Barbato C, Arisi I, Frizzo M, Brandi R, Da Sacco L, Masotti A: Computational Challenges in miRNA Target Predictions: To Be or Not to Be a True Target?. J Biomed Biotechnol. 2009, 9-Article ID 803069.
- Bartel B, Bartel D: MicroRNAs: At the Root of Plant Development?. Plant Physiol. 2003, 132: 709-717. 10.1104/pp.103.023630.PubMedPubMed CentralView Article
- Brennecke J, Stark A, Russell R, Cohen S: Principles of MicroRNA-Target Recognition. PLoS Biol. 2005, 3 (3): e85-10.1371/journal.pbio.0030085.PubMedPubMed CentralView Article
- Brodersen P, Voinnet O: Revisiting the principles of microRNA target recognition and mode of action. Nat Rev Mol Cell Bio. 10: 141-148. 10.1038/nrm2619.
- Mallory A, Bouche N: MicroRNA-directed regulation: to cleave or not to cleave. Trends Plant Sci. 2008, 13 (7): 359-367. 10.1016/j.tplants.2008.03.007.PubMedView Article
- Fujii H, Chiou TJ, Lin SI, Aung K, Zhu JK: A miRNA involved in phosphate starvation response in Arabidopsis. Curr Biol. 2005, 15: 2038-2043. 10.1016/j.cub.2005.10.016.PubMedView Article
- Jones-Rhoades MW, Bartel DP: Computational identification of plant micro-459 RNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14 (460): 787-799. 10.1016/j.molcel.2004.05.027.PubMedView Article
- Bazzini AA, Hopp HE, Beachy RN, Asurmendi S: Infection and coaccumulation of tobacco mosaic virus proteins alter microRNA levels, correlating with symptom and plant development. P Natl Acad Sci USA. 2007, 104 (29): 12157-12162. 10.1073/pnas.0705114104.View Article
- Feng J, Wang K, Liu X, Chen S, Chen J: The quantification of tomato microRNAs response to viral infection by stem-loop real-time RT-PCR. Gene. 2009, 437 (1-2): 14-21. 10.1016/j.gene.2009.01.017.PubMedView Article
- Naqvi A, Choudhury N, Mohd Q, Haq R, Mukherjee S: MicroRNAs as biomarkers in Tomato Leaf Curl Virus (ToLCV) disease. Nucl acids S. 2008, 52: 507-508. 10.1093/nass/nrn257.View Article
- Tagami Y, Inaba N, Kutsuna N, Kurihara Y, Watanabe Y: Specific enrichment of miRNAs in Arabidopsis thaliana infected with tobacco mosaic virus. DNA Res. 2007, 14: 227-233. 10.1093/dnares/dsm022.PubMedPubMed CentralView Article
- Axtell MJ: Evolution of microRNAs and their targets: are all microRNAs biologically relevant. Biochim Biophys Acta. 2008, 1779 (11): 725-34.PubMedView Article
- Sunkar R, Jagadeeswaran G: In silico identification of conserved microRNAs in large number of diverse plant species. BMC Plant Biol. 2008, 8: 37-10.1186/1471-2229-8-37.PubMedPubMed CentralView Article
- Duan C, Wang C, Huishan G: Regulation of microRNA on plant development and viral infection. Chinese Sci Bull. 2006, 51 (3): 269--278. 10.1007/s11434-006-0269-3.View Article
- Zhang B, Pan X, Cobb GP, Anderson T: Plant microRNA: A small regulatory molecule with big impact. Dev Biol. 2006, 289: 3-16. 10.1016/j.ydbio.2005.10.036.PubMedView Article
- Mahajan V, Drake A, Chen J: Virus-specific host miRNAs: antiviral defenses or promoters of persistent infection?. Trends Immunol. 30 (1):
- Shabalina A, Koonin E: Origins and evolution of eukaryotic RNA interference. Trends Ecol and Evol. 2008, 23 (10): 578-587. 10.1016/j.tree.2008.06.005.View Article
- Piriyapongs J, Jordan IK: Dual coding of siRNAs and miRNAs by plant transposable elements. RNA. 2008, 14: 814-821. 10.1261/rna.916708.View Article
- Benson D, Karsch-Mizrachi I, Lipman D, Ostell J, Wheele D: GenBank. Nucleic Acids Res. 2005, D34-D38. 33 Database
- Adams MJ, Antoniw JF: DPVweb: a comprehensive database of plant and fungal virus genes and genomes. Nucleic Acids Res. 2006, D382-D385. 10.1093/nar/gkj023. 34 Database
- Brunt AA, Crabtree K, Dallwitz MJ, Gibbs AJ, Watson L, Zurcher EJ, (eds): Plant Viruses Online: Descriptions and Lists from the VIDE Database. (1996 onwards), [http://www.agls.uidaho.edu/ebi/vdie//].
- Stajich JE, Block D, Boulez K, Brenner SE, et al: The Bioperl toolkit: Perl modules for the life sciences. Genome Res. 2002, 10: 1611-8. 10.1101/gr.361602.View Article
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Export Nucl acids S. 1999, 41: 95-98.
- Papadopoulos GL, Reczko M, Simossis VA, Sethupathy P, Hatzigeorgiou AG: The database of experimentally supported targets: a functional update of TarBase. Nucl Acids Res. 2009, D155-8. 10.1093/nar/gkn809. 37 Database
- Backman TW, Sullivan CM, Cumbie JS, Miller ZA, Chapman EJ, Fahlgren N, Givan SA, Carrington JC, Kasschau KD: Update of ASRP: the Arabidopsis Small RNA Project database. Nucleic Acids Res. 2008, 36: D982-5. 10.1093/nar/gkm997.PubMedPubMed CentralView Article
- Meyers BC, Lee DK, Vu TH, Tej SS, Edberg SB, Matvienko M, Tindell LD: Arabidopsis MPSS: an online resource for quantitative expression analysis. Plant Physiol. 2004, 135: 801-813. 10.1104/pp.104.039495.PubMedPubMed CentralView Article
- Zhang Z, Yu J, Li D, Zhang Z, Lio F, Zhou X, Wang T, Ling Y, Su Z: PMRD: plant microRNA database. Nucleic Acids Res. 2010, D806-D813. 10.1093/nar/gkp818. 38 Database
- Fahlgreen N, Howell M, Kasschau K, Chapman E, Sullivan CM, Crumbie JS, Givan SA, Law TF, Grant SR, Dangl JF, Carrington JC: High-Throughput Sequencing of Arabidopsis microRNAs: Evidence for Frequent Birth and Death of MIRNA Genes. PLoS ONE. 2007, 2: e219. 10.1371/journal.pone.0000219.View Article
- Garcia D: A miRacle in plant development: role of microRNAs in cell differentiation and patterning. Semin Cell Dev Biol. 2008, 19 (6): 586-95. 10.1016/j.semcdb.2008.07.013.PubMedView Article
- Hofacker IL, Fontana W, Stadler PF, Bonhoeffer S, Tacker P, Schuster P: Fast Folding and Comparison of RNA Secondary Structures. Monatshefte f Chemie. 1994, 125: 167-188. 10.1007/BF00818163.View Article
- Maziere P, Enright A: Prediction of microRNA targets. Drug Discov Today. 2007, 12 (11/12): 452-458. 10.1016/j.drudis.2007.04.002.PubMedView Article
- Zhang B, Pan X, Wang Q, Cobba G, Anderson T: Computational identification of microRNAs and their targets. Comput Biol Chem. 2006, 30: 395-407. 10.1016/j.compbiolchem.2006.08.006.PubMedView Article
- Schwab R, Palatnik J, Riester M, Schimmer C, Schmid M, Weigell D: Specific Effects of MicroRNA on the Plant Transcriptome. Dev Cell. 2005, 8: 517-527. 10.1016/j.devcel.2005.01.018.PubMedView Article
- R Development Core Team: R: A language and environment for statistical computing R Foundation for Statistical Computing, Vienna, Austria. 2005, [http://www.R-project.org]
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.