MicroRNA profiling of the whitefly Bemisia tabaci Middle East-Aisa Minor I following the acquisition of Tomato yellow leaf curl China virus

Bi Wang, Lanlan Wang, Fangyuan Chen, Xiuling Yang, Ming Ding, Zhongkai Zhang, Shu-Sheng Liu, Xiao-Wei Wang, Xueping Zhou

The begomoviruses are the largest and most economically important group of plant viruses exclusively vectored by whitefly (Bemisia tabaci) in a circulative, persistent manner. During this process, begomoviruses and whitefly vectors have developed close relationships and complex interactions. However, the molecular mechanisms underlying these interactions remain largely unknown, and the microRNA profiles for viruliferous and nonviruliferous whiteflies have not been studied. Sequences of Argonaute 1(Ago1) and Dicer 1 (Dcr1) genes were cloned from B. tabaci MEAM1 cDNAs. Subsequently, deep sequencing of small RNA libraries from uninfected and Tomato yellow leaf curl China virus (TYLCCNV)-infected whiteflies was performed. The conserved and novel miRNAs were identified using the release of miRBase Version 19.0 and the prediction software miRDeep2, respectively. The sequencing results of selected deregulated and novel miRNAs were further confirmed using quantitative reverse transcription-PCR. Moreover, the previously published B. tabaci MEAM1 transcriptome database and the miRNA target prediction algorithm miRanda 3.1 were utilized to predict potential targets for miRNAs. Gene Ontology (GO) analysis was also used to classify the potential enriched functional groups of their putative targets. Ago1 and Dcr1orthologs with conserved domains were identified from B. tabaci MEAM1. BLASTn searches and sequence analysis identified 112 and 136 conserved miRNAs from nonviruliferous and viruliferous whitefly libraries respectively, and a comparison of the conserved miRNAs of viruliferous and nonviruliferous whiteflies revealed 15 up- and 9 down-regulated conserved miRNAs. 7 novel miRNA candidates with secondary pre-miRNA hairpin structures were also identified. Potential targets of conserved and novel miRNAs were predicted using GO analysis, for the targets of up- and down-regulated miRNAs, eight and nine GO terms were significantly enriched. We identified Ago1 and Dcr1 orthologs from whiteflies, which indicated that miRNA-mediated silencing is present in whiteflies. Our comparative analysis of miRNAs from TYLCCNV viruliferous and nonviruliferous whiteflies revealed the relevance of deregulated miRNAs for the post-transcriptional gene regulation in these whiteflies. The potential targets of all expressed miRNAs were also predicted. These results will help to acquire a better understanding of the molecular mechanism underlying the complex interactions between begomoviruses and whiteflies.