Solanesol accumulation in plants is affected by genetic and environmental factors [2; 6; 13]. Solanesol content varies widely among tobacco varieties [2; 14]. To identify tobacco stocks with high solanesol content, we had determined the solanesol content of 168 Chinese flue-cured tobacco germplasm resources planted in four geographical regions of China in 2014 and 2015; we found that their contents were within the range of 0.70-4.13% [14]. The'major gene plus polygene' mixed inheritance model showed that solanesol content in tobacco leaves is controlled by two pairs of 'co-dominant major gene plus additive dominant polygene', with major gene transmission rates of 33.61% and 53.15% and polygene transmission rates of 30.01% and 13.64% [15]. After infection with tobacco mosaic virus (TMV) and _Pseudomonas syringae_ pv. _tabaci_, resistant varieties showed significantly higher solanesol content, while susceptible varieties exhibited no significant increases in solanesol content, indicating that solanesol plays an important role in the response of tobacco to biotic stress [16]. Previous studies have shown that drought, shade, long-wavelength light, and rare earth element treatments all lead to increased solanesol content [2; 13]. However, there have been few studies on the influence of environmental temperature on solanesol accumulation in tobacco leaves. Transcriptome analysis has been successfully used to study the molecular mechanisms of plant responses to environmental changes [17; 18; 19]. In a broad sense, the transcriptome refers to the collection of all RNAs--including messenger RNA (mRNA), ribosomal RNA, transfer RNA, and non-coding RNA--transcribed by a cell or tissue of an organism in a specific state; in a narrow sense, it refers to the collection of all mRNAs [20]. Oligo(dT) magnetic beads can be used to enrich all mRNAs transcribed by a specific tissue or cell at any specific time or location that can be sequenced using the Illumina HiSeq platform [21]. In this way, mRNA sequences and their abundances can be rapidly and comprehensively analyzed, and new genes and new transcripts can be identified [22]. Transcriptome sequencing can be used to identify differentially expressed genes (DEGs) related to different phenotypes under different treatment conditions [17; 18; 19]. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway predictions can be used for enrichment analysis and for categorizing DEGs [23]. The most active genes among DEGs can be identified to guide further research and assist in elucidating biological mechanisms [24]. Solanesol biosynthesis occurs in chloroplasts of higher plants via the 2-C-methyl-D-erythritol 4-phosphate (MEP) pathway [2; 5; 6]. Moreover, chloroplasts are organelles for plant photosynthesis, that use chlorophyll to convert light energy into chemical energy, transforming CO\({}_{2}\) and water into sugar [25]. In our previous study, 30 candidate genes involved in solanesol biosynthesis were identified, including 1-deoxy-d-xylulose 5-phosphate synthase (_DXS_), 1-deoxy-d-xylulose 5-phosphate reductoisomerase (_DXR_), 2-C-methyl-d-erythritol 4-phosphate cytidyltransferase (_IspD_), 4-diphosphocytidyl-2-C-methyl-d-erythritol kinase (_IspE_), 2-C-methyl-d-erythritol 2,4-cyclo-diphosphate synthase (_IspF_), 1-hydroxy-2-methyl-2-(_E_)-butenyl 4-diphosphate synthase(_IspG_), 1-hydroxy-2-methyl-2-(_E_)-butenyl 4-diphosphate reductase (_IspH_), isopentenyl diphosphate isomerase (_IPI_), and solanesyl diphosphate synthase (_SPS_) genes [5]. In _Arabidopsis thaliana_, _AtSPS1_ overexpression lightened the damage induced by photooxidative stress, and _AtSPS1_ and _AtSPS2_ gene silencing induced photosystem II photoinhibition [26; 27].