In Silico Identification of Cassava Genome-Encoded MicroRNAs with Predicted Potential for Targeting the ICMV-Kerala Begomoviral Pathogen of Cassava

Muhammad Aleem Ashraf; Babar Ali; Judith K Brown; Imran Shahid; Naitong Yu

doi:10.3390/v15020486

In Silico Identification of Cassava Genome-Encoded MicroRNAs with Predicted Potential for Targeting the ICMV-Kerala Begomoviral Pathogen of Cassava

Viruses. 2023 Feb 9;15(2):486. doi: 10.3390/v15020486.

Authors

Muhammad Aleem Ashraf^{1

2}, Babar Ali², Judith K Brown³, Imran Shahid⁴, Naitong Yu¹

Affiliations

¹ Institute of Tropical Biosciences and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou 571101, China.
² Institute of Biological Sciences, Faculty of Natural and Applied Sciences, Khwaja Fareed University of Engineering and Information Technology, Rahim Yar Khan 64200, Pakistan.
³ School of Plant Sciences, The University of Arizona, Tucson, AZ 85721, USA.
⁴ Department of Pharmacology and Toxicology, Faculty of Medicine, Umm Al-Qura University, Makkah 21955, Saudi Arabia.

Abstract

Cassava mosaic disease (CMD) is caused by several divergent species belonging to the genus Begomovirus (Geminiviridae) transmitted by the whitefly Bemisia tabaci cryptic species group. In India and other parts of Asia, the Indian cassava mosaic virus-Kerala (ICMV-Ker) is an emergent begomovirus of cassava causing damage that results in reduced yield loss and tuber quality. Double-stranded RNA-mediated interference (RNAi) is an evolutionary conserved mechanism in eukaryotes and highly effective, innate defense system to inhibit plant viral replication and/or translation. The objective of this study was to identify and characterize cassava genome-encoded microRNAs (mes-miRNA) that are predicted to target ICMV-Ker ssDNA-encoded mRNAs, based on four in silico algorithms: miRanda, RNA22, Tapirhybrid, and psRNA. The goal is to deploy the predicted miRNAs to trigger RNAi and develop cassava plants with resistance to ICMV-Ker. Experimentally validated mature cassava miRNA sequences (n = 175) were downloaded from the miRBase biological database and aligned with the ICMV-Ker genome. The miRNAs were evaluated for base-pairing with the cassava miRNA seed regions and to complementary binding sites within target viral mRNAs. Among the 175 locus-derived mes-miRNAs evaluated, one cassava miRNA homolog, mes-miR1446a, was identified to have a predicted miRNA target binding site, at position 2053 of the ICMV-Ker genome. To predict whether the cassava miRNA might bind predicted ICMV-Ker mRNA target(s) that could disrupt viral infection of cassava plants, a cassava locus-derived miRNA-mRNA regulatory network was constructed using Circos software. The in silico-predicted cassava locus-derived mes-miRNA-mRNA network corroborated interactions between cassava mature miRNAs and the ICMV-Ker genome that warrant in vivo analysis, which could lead to the development of ICMV-Ker resistant cassava plants.

Keywords: Indian cassava mosaic virus-Kerala; R-language; RNA interference; begomovirus; computational miRNA algorithms; non-coding RNA; virus resistance.

Publication types

Research Support, Non-U.S. Gov't

MeSH terms

Animals
Begomovirus* / genetics
Lepidoptera*
Manihot* / genetics
MicroRNAs* / genetics
RNA Interference

Substances

MicroRNAs

Supplementary concepts

Indian cassava mosaic virus

Grants and funding

This research was funded by Hainan Provincial Natural Science Foundation (321RC640), National Key R&D Program of China (2019YFD1000500), Khwaja Fareed University of Engineering and Information Technology and University of Arizona under a MoU. The APC was funded by NTY.