Tapping panel dryness (TPD), a complex physiological syndrome associated with the rubber tree (Hevea brasiliensis Muell. Arg.), causes cessation of latex drainage upon tapping and thus threatens rubber production. Rubber tree virus 1 (RTV1) is a novel positive-sense single-stranded RNA virus from the Betaflexiviridae (genus Capillovirus), which has been established to cause TPD. MicroRNAs (miRNAs) play an important role in the interplay between viruses and host cells. In this study, we identified the rubber tree genome-encoded miRNAs and their therapeutic targets against RTV1. We applied computational algorithms to predict target binding sites of rubber tree miRNAs potentially targeting RTV1 RNA genome. Mature rubber-tree miRNAs are retrieved from the miRBase database and are used for hybridization of the RTV1 genome. A total of eleven common rubber-tree miRNAs were identified based on consensus genomic positions. The consensus of four algorithms predicted the hybridization sites of the hbr-miR396a and hbr-miR398 at common genomic loci (6676 and 1840), respectively. A miRNA-regulatory network of rubber tree was constructed with the RTV1— ORFs using Circos, is illustrated to analyze therapeutic targets. Overall, this study provides the first computational evidence of the reliable miRNA–mRNA interaction between specific rubber tree miRNAs and RTV1 genomic RNA transcript. Therefore, the predicted data offer valuable evidence for the development of RTV1-resistant rubber tree in the future. Our work suggests that similar computational host miRNA prediction strategies are warranted for identification of the miRNA targets in the other viral genomes.

Tapping panel dryness (TPD) syndrome is a complex disease of Rubber tree (Hevea brasiliensis L.) which causes cessation of latex drainage upon tapping of rubber tree. Rubber tree virus (RTV1) was identified as a novel pathogen associated with rubber tree and a potential causal agent of TPD. RTV1 is a monopartite RNA virus that is linear, non-enveloped and has a single-stranded (ss) positive RNA genome of approximately 6081 nucleotides and is composed of two major open reading frames (ORFs), ORF1 (polyprotein), and ORF2 (movement protein. This study aimed to investigate the possibility of rubber genome encoded tree microRNAs (miRNAs) as novel therapeutic targets against RTV1 using in silico algorithms. Mature rubber tree miRNAs are retrieved from the miRBase database and are used for hybridization of RTV1 using five different five different computational algorithms including miRanda, RNA22, RNAhybrid and psRNATarget. A total of eleven common rubber tree miRNAs were identified based on consensus genomic positions. The consensus of four algorithms predicted the hybridization sites of hbr-miR396a and hbr-miR398 at common locus positions 6676, 1840 respectively. To validate the prediction, secondary structures of the consensual rubber tree miRNAs and free energy of duplex binding were calculated using the RNAfold and RNAcofold algorithms respectively. We created a plot between rubber tree miRNAs and RTV1 ORFs by using Circos algorithm. In this study, we predicted eleven consensual rubber tree miRNAs. Among these miRNAs, hbr-miR398 was identified as the most effectual miRNA that may target the ORF1 gene of the RTV1 genome. The predicted data will be important in the development of rubber trees resistant to RTV1.