Full-Length Transcriptome Sequencing-Based Analysis of Pinus sylvestris var. mongolica in Response to Sirex noctilio Venom

Abstract

Simple SummarySirex noctilio, as a devastating international forestry quarantine pest whose venom can cause a series of physiological changes in the host plants, such as needle wilting, yellowing, decreased transpiration rate and increased respiration rate, etc. In this study, a full-length reference transcript of Pinus sylvestris var. mongolica was constructed by combining second- and third-generation transcriptome sequencing technologies. We also identified the specific expression genes and transcription factors of P. sylvestris var. mongolica under S. noctilio venom and wounding stress. S. noctilio venom mainly induced the expression of genes related to ROS, GAPDH and GPX, and mechanical damage mainly induced the photosynthesis−related genes. The results provide a better understanding of the molecular regulation of pine trees in response to S. noctilio venom.Sirex noctilio 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 Pinus physiology and weaken or even kill the tree. In China, the Mongolian pine (Pinus sylvestris var. mongolica), an important wind-proof and sand-fixing species, is the unique host of S. noctilio. To explore the interplay between S. noctilio venom and Mongolian pine, we performed a transcriptome comparative analysis of a 10-year-old Mongolian pine after wounding and inoculation with S. noctilio 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 S. noctilio venom, such as GAPDH, GPX, CAT, FL2, CERK1, and HSP83A, etc.