Computational Biology and Machine Learning Approaches to Identify the Rubber Tree (Hevea brasiliensis L.) Genome Encoded Potential MicroRNAs Targeting Rubber Tree Virus 1

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.