[1] [1] GOODING G J, RUFTY R C. Distribution, incidence, and strainsof viruses in burley tobacco in North Carolina［J］. Plant Disease, 1987, 71: 38-40.

[2] [2] HEMENWAY C, FANG R X, KANIEWSKI W K, et al. Analysis of the mechanism of protection in transgenic plants expressing the potato virus X coat protein or its antisense RNA［J］. Embo Journal, 1988, 7: 1273-1280.

[3] [3] LUCAS G B. Diseases of tobacco, 3rd edition［M］. Raleigh: Biological Consulting Associates, 1975, 407-410.

[4] [4] BRANTLEY J D, HUNT A G. The N-terminal protein of the polyprotein encoded by the potyvirus tobacco vein mottling virus is an rna-binding protein［J］. Journal of General Virology, 1993, 74: 1157-1162.

[6] [6] PIRONE T P. Effect of tobacco vein mottling virus on yield of burley tobacco cultivars［J］. Tobacco International, 1974, 176(19): 65-66.

[7] [7] PIRONE T P, DAVIS D L. Modification of the chemical composi-tion of burley tobacco by infection with tobacco vein mottling virus［J］. Tobacco Science and Technology, 1977, 21: 83-84.

[8] [8] PIRONE T P, GOODING G J. Effect of tobacco vein mottling virus on field-growth burley tobacco varieties［J］. Plant Disease Reporter, 1973, 57: 845-847.

[9] [9] SUN P B. AUSTIN P, TOZSER J, et al. Structural determinants of tobacco vein mottling virus protease substrate specificity［J］. Protein Science, 2010, 19(11): 2240-2251.

[10] [10] PIRONE T P, GOODING G J, SMILEY J H. Tobacco vein mottling virus on burley tobacco in Kentucky［J］. Plant Disease Reporter 1973, 57: 841-844.

[11] [11] ATREYA C D, PIRONE T P. Mutational analysis of the helper component-proteinase gene of a potyvirus: effects of amino-acid substitutions, deletions, and gene replacement on virulence and aphid transmissibility［J］. Proceedings of the National Academy of Sciences of the United States of America, 1993, 90(24): 11919-11923.

[12] [12] ANSORGE W J. Next-generation DNA sequencing techniques［J］. New Biotechnology, 2009, 25(4): 195-203.

[13] [13] SCHUSTER S C. Next-generation sequencing transforms today’s biology［J］. Nature Methods, 2008, 5(1): 16-18.

[14] [14] CHEN T, LV Y, ZHAO T, et al. Comparative transcriptome profiling of a resistant vs. susceptible tomato (Solanum lycopersicum) cultivar in response to infection by tomato yellow leaf curl virus［J］. PLoS One, 2013, 8(11): e80816.

[15] [15] FABRE F, MONTARRY J, COVILLE J, et al. Modelling the evolutionary dynamics of viruses within their hosts: a case study using high-throughput sequencing［J］. PLoS Pathogens, 2012, 8(4): e1002654.

[16] [16] FAN H, SUN H, WANG Y, et al. Deep sequencing-based transcriptome profiling reveals comprehensive insights into the responses of Nicotiana benthamiana to beet necrotic yellow vein virus infections containing or lacking RNA4［J］. PLoS One, 2014, 9(1): e85284.

[17] [17] HU Q, HOLLUNDER J, NIEHL A, et al. Specific impact of tobamovirus infection on the Arabidopsis small RNA profile［J］. PLoS One, 2011, 6(5): e19549.

[18] [18] LU J, DU Z X, KONG J, et al. Transcriptome analysis of Nicotiana tabacum infected by cucumber mosaic virus during systemic symptom development［J］. PLoS One, 2012, 7: e43447.

[19] [19] ZHANG Y, PEI X, ZHANG C, et al. De novo foliar transcriptome of Chenopodium amaranticolor and analysis of its gene expression during virus-induced hypersensitive response［J］. PLoS One, 2012, 7(9): e45953.

[20] [20] ZHENG W, MA L, ZHAO J, et al. Comparative transcriptome analysis of two rice varieties in response to rice stripe virus and small brown planthoppers during early interaction［J］. PLoS One, 2013, 8(12): e82126.

[25] [25] MARTIN M. Cutadapt removes adapter sequences from high-throughput sequencing reads［J］. Embnet Journal, 2011, 17: 10-12.

[26] [26] DAEHWAN K, BEN L, STEVEN L, et al. A fast spliced aligner with low memory requirements［J］. Nature Methods, 2015, 12: 357-360.

[27] [27] MENG J. A protocol for RNA methylation differential analysis with MeRIP-Seq data and exomePeak R/Bioconductor package［J］. Methods, 2014, 69(3): 274-281.

[28] [28] BAILEY T L, MIKAEL B, BUSKE F A, et al. MEME SUITE: tools for motif discovery and searching［J］. Nucleic Acids Research, 2009, 37(Web Server issue): W202-208.

[29] [29] HEINZ S. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and b cell identities［J］. Molecular Cell, 2010, 38(4): 576.

[30] [30] YU G. An R/Bioconductor package for ChIP peakannotation, comparison and visualization［J］. Bioinformatics, 2015, 31(14): 2382-2383.

[31] [31] PERTEA M. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads［J］. Nature Biotechnology, 2015, 33: 290-295.

[32] [32] ROBINSON M D. A bioconductor package for differential expression analysis of digital gene expression data［J］. Bioinformatics, 2010, 26(1): 139.