Full-Length Transcriptome Characterization and Functional Analysis of Pathogenesis-Related Proteins in Lilium Oriental Hybrid ‘Sorbonne’ Infected with Botrytis elliptica

Gray mold (Botrytis elliptica) causes a deleterious fungal disease that decreases the ornamental value and yield of lilies. Lilium oriental hybrid ‘Sorbonne’ is a variety that is resistant to gray mold. Understanding the mechanism of resistance against B. elliptica infection in ‘Sorbonne’ can provide a basis for the genetic improvement in lily plants. In this study, a PacBio Sequel II system was used to sequence the full-length transcriptome of Lilium ‘Sorbonne’ after inoculation with B. elliptica. A total of 46.64 Gb subreads and 19,102 isoforms with an average length of 1598 bp were obtained. A prediction analysis revealed 263 lncRNAs, and 805 transcription factors, 4478 simple sequence repeats, and 17,752 coding sequences were identified. Pathogenesis-related proteins (PR), which may play important roles in resistance against B. elliptica infection, were identified based on the full-length transcriptome data and previously obtained second-generation transcriptome data. Nine non-redundant potential LhSorPR proteins were identified and assigned to two groups that were composed of two LhSorPR4 and seven LhSorPR10 proteins based on their genetic relatedness. The real-time quantitative reverse transcription PCR (qRT-PCR) results showed that the patterns of expression of nine differentially expressed PR genes under B. elliptica stress were basically consistent with the results of transcriptome sequencing. The pattern of expression of LhSorPR4s and LhSorPR10s genes in different tissues was analyzed, and the expression of each gene varied. Furthermore, we verified the function of LhSorPR4-2 gene in Lilium. The expression of LhSorPR4-2 was induced by phytohormones such as methyl jasmonate, salicylic acid, and ethephon. Moreover, the promoter region of LhSorPR4-2 was characterized by several functional domains associated with phytohormones and stress response. The overexpression of LhSorPR4-2 gene in ‘Sorbonne’ increased the resistance of the lily plant to B. elliptica and correlated with high chitinase activity. This study provides a full-length transcript database and functionally analyzed the resistance of PR gene to B. elliptica in Lilium, thereby introducing the candidate gene LhSorPR4-2 to breed resistance in Lilium.