Early transcriptional response of terpenoid metabolism to Colletotrichum gloeosporioides in a resistant wild strawberry Fragaria nilgerrensis Author Mehmood, Nasir * & College of Horticulture and the Fujian provincial Key Laboratory of Plant Functional Biology, Fujian Agriculture and Forestry University, Fuzhou, 350002, China Author Yuan, Yuan * & College of Horticulture and the Fujian provincial Key Laboratory of Plant Functional Biology, Fujian Agriculture and Forestry University, Fuzhou, 350002, China Author Ali, Mohammed Egyptian Deserts Gene Bank, Department of Genetic Resources, Desert Research Center, Egypt Author Ali, Muhammad College of Agriculture and Biotechnology, Zhejiang University, Hangzhou, 310058, China Author Iftikhar, Junaid * & College of Horticulture and the Fujian provincial Key Laboratory of Plant Functional Biology, Fujian Agriculture and Forestry University, Fuzhou, 350002, China Author Cheng, Chunzhen * & College of Horticulture and the Fujian provincial Key Laboratory of Plant Functional Biology, Fujian Agriculture and Forestry University, Fuzhou, 350002, China Author Lyu, Meiling * & College of Horticulture and the Fujian provincial Key Laboratory of Plant Functional Biology, Fujian Agriculture and Forestry University, Fuzhou, 350002, China Author Wu, Binghua * & College of Horticulture and the Fujian provincial Key Laboratory of Plant Functional Biology, Fujian Agriculture and Forestry University, Fuzhou, 350002, China text Phytochemistry 2021 112590 2021-01-31 181 1 12 http://dx.doi.org/10.1016/j.phytochem.2020.112590 journal article 10.1016/j.phytochem.2020.112590 1873-3700 8290677 2.3. Isoprenoid biosynthesis genes from F. nilgerrensis leaf transcriptome Fruits of wild and cultivated varieties of strawberry contains numerous types of terpenoids as the major components of the aroma essential oils, represented mainly by myrtenol, α- pinene, β- pinene, myrtenyl acetate, sabinene, β- myrcene, α- phellandrene, β- phellandrene, dihydromyrcenol, α- terpinolene, α- terpineol, linalool, nerolidol, ()-limonene and (+)-limonene ( Aharoni et al., 2004 ). However, the composition of terpenoids in leave are rarely reported. Our GC/MS profile unveils that the leaf contained low amount (less than 2% of total extracts) of terpenoids in total. However, upon challenged with C. gloeosporioides , rapid induction in terpenoid biosynthesis is quite impressed. Thus, we played particular attention to the fungus-elicited transcriptional regulation of terpenoid biosynthesis network in the leaf based on our transcriptome data. Plant terpenoids biosynthesis are initiated from the plastidial MEP (methylerythritol phosphate) and the cytosolic/peroxisomal MVA (mevalonic acid) pathways that generate the building blocks isopentyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) for diverse downstream conversions ( Fig. 4 ). Therefore, we first identified unigenes coding for enzymes in various steps of terpenoid biosynthesis pathways in the transcriptome, such as isopentyl diphosphate synthases (IPPS), dimethylallyl diphosphate synthases (DMADPS), geranyl diphosphate synthases (GPPS), farnesyl diphosphate synthases (FPPS) and geranylgeranyl diphosphate synthases (GGPPS). In addition, most unigenes associated with biosynthesis of isoprenoid in both MEP and MVA pathways were found, including FnDXR (1-deoxy-D-xylulose-5- phosphate reductoisomerase), FnISPF (2-C-methyl-D-erythritol-2,4- cyclodiphosphate synthase), FnHDS1 ((E)-4-hydroxy-3-methylbut- 2- enyl-diphosphate synthase 1), FnHDR (4-hydroxy-3-methylbut-2-enyl diphosphate reductase), FnIDI1 (isopentenyl-diphosphate deltaisomerase 1), FnAACT1 (acetyl-CoA C-acetyltransferase 1), FnHMGS (hydroxylmethylglutaryl-CoA synthase), FnHMGR4 and 5 (hydroxymethylglutaryl-CoA reductase (NADPH) 4 and 5), FnMDC (diphosphomevalonate decarboxylase). In total 57 unigenes related to MEP, MVA and other terpenoid backbone biosynthesis pathway were identified and assigned gene names based on homology (Table S4 and Fig. 4 ). Terpene synthase (TPS) catalysis the final step for various terpenoid molecules. In total 59 homologous TPS unigenes from our dataset were identified and determined on the basis of TPS sequence resemblances that were discovered in the reference of canonical annotation database (Table S5 and Fig. 4 ). Fifteen of these unigenes were annotated as monoterpenes, limonene and pinene degradation biosynthesis, such as (3 S)-linalool synthase, (+)-neomenthol dehydrogenase, ()-alphaterpineol synthase, ()-camphene/tricyclene synthase and aldehyde dehydrogenase (NAD + ). Twenty-three unigenes were annotated as sesquiterpene biosynthesis, such as farnesyl-diphosphate farnesyltransferase, NAD + -dependent farnesol dehydrogenase, alpha-farnesene synthase, ()-germacrene D synthase and beta-amyrin synthase. Also, twelve unigenes were annotated as diterpene biosynthesis, such as ent-copalyl diphosphate synthase, ent-kaurene synthase, ent-kaurene oxidase, ent-kaurenoic acid monooxygenase, gibberellin-44 dioxygenase, gibberellin 3 beta-dioxygenase and gibberellin 2 beta-dioxygenase. Finally, nine unigenes related to squalene monooxygenase were annotated as triterpene biosynthesis (Table S5 and Fig. 4 ). Fig. 2. Global differentially-expressed-genes (DEG) in leaves of F. nilgerrensis after challenged by C. gloeosporioides . (A) DEG heatmap depicts clustering of the early upregulated genes (up to 18-h post challenge) and the late respondents ( > 24 h) in the upper and lower parts, respectively, with top GO terms indicated. (B) Major GO term enrichment in an early responsive gene cluster 7 showing that terpenoid biosynthesis genes are co-expressed with antipathogen immune-responsive pathway genes. Each time-point contains 2 biological replications and the expression data was zero-normalized before hierarchical clustering is performed (Pearson uncentered metric, average linkage ordering) using Cluster 3.0 version 1.58. Heatmap was examined and generated using Java Tree- View version 1.1.6r4. Gene Ontology (GO) enrichment analysis of the DEGs is implemented by the topGO R package-based Kolmogorov–Smirnov test. Based on the transcript sequences, we also detected simple sequence repeat (SSR) motifs in the identified unigenes associating with terpenoid biosynthesis, e.g. Hmgcr 4, FnDXS 3, FnDXR, FnMCT, FnIDI 1, FnFPPS 1, FnGGPPS 1, FnCMT, FnFNTB 2, FnSPS 1 , and unigenes related to (3 S)- linalool synthase, aldehyde dehydrogenase (NAD + ), 3, 4, 6, ent-copalyl diphosphate synthase, ent-kaurene synthase 1, ent-karene oxidase, ent-kaurenoic acid monooxygenase 2, gibberellin-3-beta dioxygenase, farnesyl-diphosphate, farnesyltransferase, squalene monooxygenase 2, and beta-amyrin synthase 3 (Table S6). The identified homologous genes provide an overview for further deciphering terpenoid biosynthesis in leaves of the N. nilgerrensis and may be a reference for other strawberry cultivars as well ( Fig. 4 ). Next, we examined the expression profiles of these unigenes during the timecourse after fungal challenge. In the general terpenoid biosynthesis pathways, more unigenes from the MVA, MEP and those related to diphosphate synthases (PPS) showed up-regulation as early at 3 hpi with various fold-change levels, while only a few unigenes were repressed, especially c27310. graph_c0 ( FnDHDDS4 ) and c15242. graph_c0 ( FnDHDDS1 ) ( Fig. 5A ). For TPS -like unigenes, ten unigenes ( c26859. graph_c0 , c23001. graph_c0 , c26382. graph_c2 , c24450. graph_c0 , c34218. graph_c0 , c34395. graph_c0 , c39471. graph_c0 , c31932. graph_c0 , c24276. graph_c0 and c27886. graph_c0 ) were up-regulated remarkedly at 3–6 hpi, and further fourteen unigenes ( c24788. graph_c1 , c27601. graph_c0 , c9433. graph_c0 , c26113. graph_c2 , c8446. graph_c0 , c26141. graph_c0 , c13323. graph_c0 , c11811. graph_c0 , c9937. graph_c0 , 18669. graph_c0 , c28365. graph_c0 , c1733. graph_c0 , c45499. graph_c0 , c19725. graph_c0 and c19236. graph_c0 ) were induced significantly at 12–18 hpi ( Fig. 5B ). These results indicate that in response to C. gloeosporioides challenge, terpenoid biosynthesis genes are among the most active transcripts of the specialised metabolism during the early stage.