Full-Length Transcriptome Sequencing-Based Analysis of var. in Response to Venom.

is a major international quarantine pest that recently emerged in northeast China to specifically invade conifers. During female oviposition, venom is injected into the host together with its symbiotic fungus to alter the normal physiology and weaken or even kill the tree. In China, the Mongolian pine ( var. ), an important wind-proof and sand-fixing species, is the unique host of . To explore the interplay between venom and Mongolian pine, we performed a transcriptome comparative analysis of a 10-year-old Mongolian pine after wounding and inoculation with venom. The analysis was performed at 12 h, 24 h and 72 h. PacBio ISO-seq was used and integrated with RNA-seq to construct an accurate full-length transcriptomic database. We obtained 52,963 high-precision unigenes, consisting of 48,654 (91.86%) unigenes that were BLASTed to known sequences in the public database and 4309 unigenes without any annotation information, which were presumed to be new genes. The number of differentially expressed genes (DEGs) increased with the treatment time, and the DEGs were most abundant at 72 h. A total of 706 inoculation-specific DEGs (475 upregulated and 231 downregulated) and 387 wounding-specific DEGs (183 upregulated and 204 downregulated) were identified compared with the control. Under venom stress, we identified 6 DEGs associated with reactive oxygen species (ROS) and 20 resistance genes in Mongolian pine. Overall, 52 transcription factors (TFs) were found under venom stress, 45 of which belonged to the AP2/ERF TF family and were upregulated. A total of 13 genes related to the photosystem, 3 genes related photo-regulation, and 9 TFs were identified under wounding stress. In conclusion, several novel putative genes were found in Mongolian pine by PacBio ISO seq. Meanwhile, we also identified various genes that were resistant to venom, such as GAPDH, GPX, CAT, FL2, CERK1, and HSP83A, etc.