Transcriptome Analysis of Nicotiana Benthamiana Infected by Tobacco Curly Shoot Virus
Ke Li,Gentu Wu,Mingjun Li,Mingge Ma,Jiang Du,Miao Sun,Xianchao Sun,Ling Qing
DOI: https://doi.org/10.1186/s12985-018-1044-1
IF: 5.916
2018-01-01
Virology Journal
Abstract:Background: Tobacco curly shoot virus (TbCSV) is a monopartite begomovirus associated with betasatellite (Tobacco curly shoot betasatellite, TbCSB), which causes serious leaf curl disease on tomato and tobacco in China. It is interesting that TbCSV induced severe upward leaf curling in Nicotiana benthamiana, but in the presence of TbCSB, symptoms changed to be downward leaf curling. However, the mechanism of interactions between viral pathogenicity, host defense, viral-betasatellite interactions and virus-host interactions remains unclear. Methods: In this study, RNA-seq was used to analyze differentially expressed genes (DEGs) in N. benthamiana plants infected by TbCSV (Y35A) and TbCSV together with TbCSB (Y35AB) respectively. Results: Through mapping to N. benthamiana reference genome, 59,814 unigenes were identified. Transcriptome analysis revealed that a total of 4081 and 3196 DEGs were identified in Y35AB vs CK (control check) and Y35A vs CK, respectively. Both GO and KEGG analyses were conducted to classify the DEGs. Ten of the top 15 GO terms were enriched in both DEGs of Y35AB vs CK and Y35A vs CK, and these enriched GO terms mainly classified into three categories including biological process, cellular component and molecular function. KEGG pathway analysis indicated that 118 and 111 pathways were identified in Y35AB vs CK and Y35A vs CK, respectively, of which nine and six pathways were significantly enriched. Three major pathways in Y35AB vs CK involved in metabolic pathways, carbon metabolism and photosynthesis, while those in Y35A vs CK were related to Ribosome, Glyoxylate and dicarboxylate metabolism and DNA replication. We observed that 8 PR genes were significantly up-regulated and 44 LRR-RLK genes were significantly differentially expressed in Y35A treatment or in Y35AB treatment. In addition, 7 and 13 genes were identified to be significantly changed in biosynthesis and signal transduction pathway of brassinosteroid (BR) and jasmonic acid (JA) respectively. Conclusions: These results presented here would be particularly useful to further elucidate the response of the host plant against virus infection.