Next Article in Journal
Ecological Restoration in Eastern Canada Using Four Early-Successional Species on Severely Degraded Sites Using a Factorial of Site-Preparation Treatments: Growth and Biomass over Two Years
Previous Article in Journal
Simulation of Radial Growth of Mongolian Pine (Pinus sylvestris var. mongolica) and Dahurian Larch (Larix gmelinii) Using the Vaganov–Shashkin Model in the Greater Khingan Range, Northeast China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Functional Characterization of Terpene Synthases from Masson Pine (Pinus massoniana) under Feeding of Monochamus alternatus Adults

1
Co-Innovation Center for the Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China
2
College of Forestry, Nanjing Forestry University, Nanjing 210037, China
3
School of Landscape Architecture, Jiangsu Vocational College of Agriculture and Forestry, Zhenjiang 212499, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Forests 2024, 15(2), 244; https://doi.org/10.3390/f15020244
Submission received: 19 December 2023 / Revised: 19 January 2024 / Accepted: 25 January 2024 / Published: 27 January 2024
(This article belongs to the Section Forest Health)

Abstract

:
Conifers have evolved sophisticated terpenoid defenses for protection against herbivores and pathogens. Pinus massoniana Lamb. is the most widely distributed pioneer afforestation and resin tree species in China, but is seriously harmed by pine wilt disease. Monochamus alternatus is the main vector of pine wilt disease in China. Monoterpenes, sesquiterpenes and diterpenes, the main secondary defensive compounds of P. massoniana, are catalyzed by different terpene synthases (TPSs), which participate in the important defense pathways against external biotic and abiotic stresses. Here, we aimed to identify the terpene synthases (TPSs) in P. massoniana, responding to the feeding of M. alternatus, and to characterize the functions and products of the mono-TPSs. We identified six differentially expressed TPS genes in the P. massoniana fed upon by M. alternatus, including four mono-TPS and two sesqui-TPS genes. The functions of the four mono-TPSs were verified by analysis of the main product and by-products of these mono-TPSs. (+)-α-Pinene, (−)-α-pinene, and limonene were the major products of TPS (+)-α-pinene, TPS (−)-α-pinene, and TPS limonene, respectively, but TPS (−)-β-pinene only catalyzed a trace amount of (−)-β-pinene in the products. Our findings shed light on the potential relationships between the structure of terpene synthases and their corresponding products.

1. Introduction

Conifers rely on constitutive and induced defense systems, primarily based on terpenoids, against various biotic and abiotic stresses [1,2,3]. The coniferous terpenoid defense system consists of monoterpenoids, sesquiterpenoids, diterpenoids, and the corresponding resin ducts and cells [4,5]. Conifers can actively change the composition of terpenoids to defend against insect herbivores, affecting insect behavior and performance or attracting natural enemies [6,7,8]. For example, a methyl jasmonate-stimulated feeding treatment affected the release of terpene volatiles from five different conifers [9]. Methyl jasmonate treatment changes the chemical composition of Pinus sylvestris and influences the meal properties of Hylobius abietis [10]. Many herbivorous insects have evolved to use host-produced terpenoids, particularly volatile monoterpenoids and sesquiterpenoids, as critical chemical signals for host recognition and location. For instance, Dendroctonus frontalis and D. terebrans have olfactory and behavioral responses to volatiles (α-pinene, β-pinene, myrcene, limonene, and 4-allylanisole) released from the resins excreted by their hosts [11]. Furthermore, Xyleborus glabratus is attracted to sesquiterpenes, such as α-copaene, α-cubebene, α-humulene, and calamenene, emitted from the cambium of their lauraceous hosts [12].
The components and composition of constitutive terpenoids vary under different biotic and abiotic stresses [13], which is directly related to the distinct gene expression levels of terpene synthases (TPSs) under different conditions [14,15]. In Picea abies and Pinus massoniana (Pm), terpenoid biosynthesis pathways are activated under feeding stress or methyl jasmonate induction. This activation leads to upregulated gene expression of mono-TPS PaJF67 and Pm TPSs, which finally enhances the biosynthesis of some terpenoids, such as α-pinene, β-pinene, and (+)-3-Carene [16,17,18]. Identification and characterization of TPS genes for genetic improvement can improve the genetic improvement of plants and enhance resistance against invaders. For example, transgenic soybean with overexpression of GmAFS showed significantly improved resistance to soybean cyst nematode [19]; expression of the 7-epizingiberene synthase and Z-Z-farnesyl-diphosphate synthases in the glandular trichomes of cultivated tomato resulted in the resistance of tomatoes to several herbivore pests [20].
In plants, TPSs catalyze the formation of many monoterpenes, sesquiterpenes, and diterpenes with different structures, using precursors [21]. Monoterpenes, sesquiterpenes, and diterpenes are synthesized from geranyl diphosphate (GPP), farnesyl diphosphate (FPP), and geranylgeranyl diphosphate, respectively. They are, in turn, synthesized from different numbers of C5s by isoprenyl transferase. Terpenoids usually exist as enantiomers but often with one enantiomer as the dominant product [22]. Terpene synthases can produce one or more products [23]. The production of sesquiterpenes is usually lower than that of monoterpenes and diterpenes but has the highest structural diversity in coniferous trees. Due to difficulty in the structural identification of sesquiterpenes and the diversity of sesquiterpene synthase products, the physiology and biochemistry of sesquiterpene synthases have barely been studied in coniferous trees. Phylogenetic analysis based on protein sequences has divided gymnosperm and some angiosperm TPSs into seven subfamilies: TPS-a, TPS-b, TPS-c, TPS-e/f, TPS-g, and TPS-h [24,25,26,27,28]. The TPS-d subfamily is entirely composed of TPSs from gymnosperms that are involved in primary and secondary metabolism [21]. Some TPS genes in P. abies, P.sitchensis, P. glauca, Abies grandis, Pseudotsuga menziesii, P. taeda, P. contorta, and P. banksiana have been identified, cloned, and functionally verified [29,30,31,32,33,34].
Monochamus alternatus is a wood-boring insect, widely distributed from East to South Asia, including China, Japan, Korea, Laos, Taiwan, Vietnam, and some other parts of the world, including Canada, Denmark, and the United States [35]. The larvae of this beetle feed and pupate in the phloem and xylem of host plants, the newly emerged adults mate after pre-mating feeding, and the female adults spawn in the ovipositing incisions [36] (Figure 1). M. alternatus can spread pine wood nematodes (Bursaphelenchus xylophilus), which are the causal agent of pine wilt disease, during adult feeding [37,38]. P. massoniana is an important tree species widely distributed in southern China because of its high tolerance of acidic soils and high oleoresin yield, with a distribution area of about 20 million square kilometers in China. However, P. massoniana is the preferred host of M. alternatus and is seriously damaged by pine wilt disease. Since the disease was first discovered in China in 1982, it has killed billions of pine trees, causing huge economic losses, and is the main threat to the healthy pine forests widely distributed in China [39].
In our previous study, the terpene-based defense of P. massoniana was induced by M. alternatus adult feeding. Transcriptome sequencing and paraffin sections showed that most TPSs of P. massoniana showed increased expression, which was followed by the formation of traumatic resin ducts and release of terpenoids [18,40]. Among the released volatile terpenoids, monoterpenoids and sesquiterpenoids had the highest contents, including α-pinene, β-phellandrene, β-pinene, β-myrcene, camphene, caryophyllene, and longifolene. In addition, many diterpenes are released. We also proved that methyl jasmonate induction significantly promoted the release of volatile terpenoids from P. massoniana [41]. At the same time, the pre-induction of methyl jasmonate significantly reduced the feeding area of the M. alternatus adult [41]. These findings indicate that the terpenoids are important secondary metabolites of P. massoniana in response to M. alternatus adult feeding, and TPSs are involved in the biosynthesis of terpenoids. Terpenoids play a major role in M. alternatus host and mate localization, as well as in chemical communication among the pine wood nematode, M. alternatus, and P. massoniana [42,43,44]. Some TPSs in P. massoniana have been identified, such as Pm GGPPS2, which affects the synthesis of diterpenoids [45], as well as Pm TPS4 and Pm TPS21, which regulate the synthesis of α-pinene and longifolene and participate in the defense against pine wood nematodes [46]. Nevertheless, the links between the structures of TPSs and their products in P. massoniana remain unclear.
In this study, we first modelled the structures of six TPSs in P. massoniana that had been found to be significantly induced by M. alternatus adult feeding, including four mono-TPS genes, TPS (+)-α-pinene, TPS (−)-α-pinene, TPS (−)-β-pinene, and TPS limonene, and two sesqui-TPS genes, TPS selinene and TPS longifolene, after cloning and analyzing their sequences. Prokaryotic expression was then performed for monoterpene synthesis to determine their final products. The results from this study will help reveal the potential correlations between the structures of TPSs and their corresponding terpene products, which will provide a theoretical basis for the transgenic breeding of M. alternatus resistant Masson pines to control pine wilt disease.

2. Materials and Methods

2.1. Plant and Insect Materials

A two-year-old P. massoniana seedling from Pingxiang City, Jiangxi Province (26°76′ N, 114°28′ E) in China was planted individually in a greenhouse (26 ± 0.5 °C, relative humidity = 70 ± 5%, 16:8 h light/dark photoperiod) with nutrient medium (humus/turf soil/perlite, 3:1:1 by volume) and weekly deep watering [40]. The third-instar larvae of M. alternatus were collected from P. massoniana host trees in Quanjiao County, Anhui Province, China (41°31′ N, 117°74′ E). Subsequently, they were reared individually on an artificial diet at a stationary temperature of 26 ± 0.5 °C, a relative humidity of 60 ± 5%, and a 16:8 h light/dark photoperiod, as described by Chen Rui-Xu et al. [40]. Newly emerged adults were collected and individually kept in 250-mL glass flasks, provisioned with fresh P. massoniana twigs as food in climate-controlled rooms (26 ± 0.5 °C, relative humidity = 60 ± 5%, 16:8 h light/dark photoperiod). Adult beetles (mixture of male and female), visually estimated to be of similar physiological conditions, were chosen, and starved for 24 h before use. Finally, the plant materials were authenticated by Prof. Kongshu Ji from Nanjing Forestry University, China. The voucher specimen (NF1001839) was deposited at Dendrological Herbarium, Nanjing Forestry University (Institution code from Chinese Virtual Herbarium).

2.2. Sampling

The feeding treatment was administered according to Chen Rui-Xu et al. [18]. Healthy P. massoniana seedlings under similar growth conditions were selected and divided into treatment and control groups. The starved adults of M. alternatus were caged individually on a seedling stem using a wire mesh (25 cm × 25 cm, hole size = 4 mm2) surrounding the stem of P. massoniana seedlings. After feeding for 24 h, the wire mesh and insects were removed (N = 3). Subsequently, the seedlings of the treatment and control groups were sampled. The stems of P. massoniana, 5 cm above and below the feeding center of M. alternatus, were cut into small sections (1–2 cm long), transferred to ribonuclease-free centrifuge tubes, and quickly frozen in liquid nitrogen. The stems were stored at −80 °C.

2.3. RNA Extraction and Gene Sequence Analysis

Pm TPS gene sequences and corresponding gene-specific primers were derived from the previously determined M. alternatus feeding transcriptome [18] (BioProject accession number PRJNA328403). Four monoterpene synthase genes, unigene 26781, unigene 36618, unigene 26791, and unigene 18558, and two sesquiterpene synthase genes, unigene 380 and unigene 12384, were selected for amplification of the coding sequence (CDS). The total RNA was extracted from the sections of P. massoniana using a Miniprep RNA Purification Kit (TIANGEN, Beijing, China) and reverse-transcribed to synthesize complementary DNA (cDNA) using a 5 × All-In-One RT Master Mix (Accurate Biology, China). The SMARTer® Rapid Amplification of cDNA Ends 5′/3′ Kit (Clontech Laboratories, Inc., Mountain View, CA, USA) was used for 5′ or 3′ rapid amplification of cDNA via polymerase chain reaction to obtain the complete sequences. The amplified products were recovered using the MiniBEST Agarose Gel DNA Extraction Kit Ver.4.0 (Takara Bio Inc., Otsu, Japan). The amplified sequences were ligated to a pCE2 TA/Blunt-Zero vector (Vazyme, Nanjing, China) to identify full-length sequences. The full-length sequences of Pm TPS (+)-α-pinene, Pm TPS (−)-α-pinene, Pm TPS (−)-β-pinene, Pm TPS limonene, Pm TPS selinene, and TPS longifolene were submitted to GenBank (accession numbers MN652001.1, MN652002.1, MN652003.1, MN652005.1, MN652006.1, MN652007.1). The basic local alignment search tool at the National Center for Biotechnology Information was used to identify sequences homologous to Pm TPSs (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 13 January 2022).

2.4. Bioinformatics Analysis

The protein sequence alignment was conducted using DNAMAN 6.0 software (Lynnon Biosoft, San Ramon, CA, USA. Version number 6.0.3.99). The online platform ExPASY https://web.expasy.org/compute_pi/, accessed on 1 December 2021) was used to calculate isoelectric point, molecular weight, and hydrophobic/hydrophilic prediction. In addition, the online tool TargetP 1.1 (http://www.cbs.dtu.dk/services/TargetP/, accessed on 1 December 2021) was used to predict the amino-terminal (N-terminal) signal peptide sequence. The coding sequences of 85 monoterpene and sesquiterpene synthases from pine trees were downloaded from the NCBI database (Table S1) and aligned using Clustal W software (https://www.genome.jp/tools-bin/clustalw, accessed on 13 December 2021). MEGA 7.0 software (Version 7.0.26) was used to construct the phylogenetic tree as per the maximum likelihood method with 500 bootstrap replicates [47].
SOPMA online software (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html, accessed on 26 December 2021) was used to analyze and predict the secondary structures of Pm TPS proteins. In addition, the online software MEME (Multiple EM for Motif Elicitation) (https://meme-suite.org/meme/tools/meme, accessed on 11 January 2023) was used to predict and analyze the motifs of Pm TPS genes. The motif width was set to 2–300, and the number of motifs was set to 20. The amino acid sequences of 91 terpene synthases in the phylogenetic tree were used for motif analysis.

2.5. Homology Modeling and Molecular Docking

Homology models of six Pm TPSs were built using the Swiss-Model online website (https://swissmodel.expasy.org, accessed on 11 January 2022) based on the tertiary structure of α-bisabolene synthase from A. grandis (protein data bank (PDB)no. 3sae.1. A) and taxadiene synthase from Taxus brevifolia (PDB no. 3p5r.1. A). Modeled structures were confirmed to be of high quality, with Ramachandran plot statistics exceeding 90% using PROCHECK [48]. In addition, the substrates GPP and FPP were docked with the model structures of Pm mono-TPSs and sesqui-TPSs, respectively, using SYBYL-X software version 2.0 (Certara, Princeton, NJ, USA).

2.6. Expression and Purification of Recombinant Pm Mono-TPSs

Six Pm TPS-coding sequences were amplified as described above, and a homologous recombination insertion into the pET-32a (+) vector was employed. Homologous recombinant plasmids were verified via sequencing and transferred into BL21 (DE3) competent cells. Luria–Bertani medium (1 L) containing 100 μg/mL Amp was inoculated with five individual colonies and cultured until an optical density at 600 nm (OD600) of at least 0.8 was reached at 37 °C and 220 rpm. Isopropyl-β-D-thiogalactopyranoside was added to a final concentration of 0.5 mM and induced to express at 20 °C and 100 rpm for 12 h. Cultures were collected via centrifugation at 13,523× g for 15 min. The bacterial cells were collected and resuspended in 100 mL resuspension buffer (50 mM PB, 300 mM NaCl, 20 mM imidazole, pH 7.4). Phenylmethylsulfonyl fluoride was added to a final concentration of 1 mM and ultrasonically broken in an ice bath (power 200 W, work 3 s, pause 5 s, for 10 min per 30 mL of resuspended bacterial cells). After centrifugation at 13,523× g for 30 min, the supernatant and precipitate were collected (dissolved in 8 M urea). Subsequently, 20 μL of each was added to 5 μL loading buffer, kept in a 100 °C water bath for 10 min, and then centrifuged for a short time. Sodium dodecyl-sulfate polyacrylamide gel electrophoresis was performed on 7.5% gel (100 V, 120 min). The protein gel was stained with Coomassie brilliant blue and photographed. Recombinant Pm mono-TPS proteins were purified to homogeneity by Ni2+-affinity chromatography and then desalted at 4 °C for 16 h in a dialysis bag (14 kDa) with a buffer composed of 25 mM 4-(2-hydroxyethyl)-1-piperazineethanesulphonic acid, 100 mM NaCl, and 1 mM dithiothreitol (DTT). The dialysis target protein was frozen at −80 °C for 4 h and treated in freeze-drying apparatus for more than 24 h until the protein was completely dried. Dried terpene synthase protein was diluted to approximately 5 mg/mL using ddH2O.

2.7. TPS Enzyme Assays and Gas Chromatography–Flame Ionization Detection/Gas Chromatography–Mass Spectrometry (GC–FID/GC–MS) Analysis

The TPS enzyme activities were tested according to a previously published procedure (Martin et al. 2002), with modifications. Before the assay, the frozen proteins were placed at 37 °C until thawing. The diluted proteins were pre-equilibrated with assay buffers (25 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (pH 7.4), 100 mM KC1, 10 mM MgCl2, 5% (v/v) glycerol, and 5 mM DTT). Enzyme activity was assessed using 1 mL of the pre-equilibrated proteins with the addition of 25 μL GPP (50 mM, with 50% v/v methanol, perform nitrogen purge before use). All enzyme assays were overlaid with 800 µL of hexane (Aladdin Biological Technology Co., Ltd., Shanghai, China) and 400 ng/mL isobutyl benzene (Aladdin Biological Technology Co., Ltd.) to collect the released volatiles, then incubated at 37 °C, 80 rpm for 1 h, 3 h, and 5 h. The upper layer hexane (200 µL) was taken into a sample vial for GC–FID/GC–MS analysis. Chromatographic grade n-hexane, (+)-α-pinene, (−)-α-pinene, (−)-β-pinene, R-(+)-limonene, linalool, and myrcene standard were purchased from Aladdin Biological Technology Co., Ltd.
The Pm mono-TPS products were analyzed on a GC–FID fitted with an HP-5 column (30 m × 0.25 mm × 0.25 μm, Agilent Technologies Inc., Santa Clara, CA, USA). The injection volume was 1 μL of the hexane extract. The oven temperature was programmed from an initial temperature of 50 °C for 3 min to 250 °C at increments of 10 °C min−1. The GC–MS analysis was performed using the same program as that for GC–FID. The MS conditions were obtained as follows: an ion source, 230 °C; electron energy, 70 eV; and a scan range, 50–500 mass units; helium was used as the carrier gas, 1.0 mL/min. Molecule identification was performed using the National Institute of Standards and Technology 14 library. The results were confirmed by comparing the retention times of target monoterpenes with authentic standards, including (+)-α-pinene, (−)-α-pinene, (−)-β-pinene, R-(+)-limonene, linalool, and myrcene (Aladdin Biological Technology Co., Ltd.). If available, identification was confirmed by the Kovats index or authentic standards [49]. The C7-C40 saturated alkanes were purchased from Sigma-Aldrich (St. Louis, MO, USA).

3. Results

3.1. Sequence, Phylogenetic, and Motif Analysis of Six Pm TPSs

To characterize the terpene synthase genes involved in the M. alternatus feeding induced defense in P. massoniana, six genes, Pm TPS (+)-α-pinene, Pm TPS (−)-α-pinene, Pm TPS (−)-β-pinene, Pm TPS limonene, Pm TPS selinene, and Pm TPS longifolene, which were differentially expressed after M. alternatus inoculation, were identified in the P. massoniana transcriptome [38]. The BLAST alignment showed hits with known TPS genes based on the presence of conserved sequence characteristics shared by all TPSs. The six predicted proteins encoded by the above Pm TPSs ranged between 580–635 amino acids in size, with predicted pIs of 5.71, 5.53, 5.74, 5.44, 5.6, and 5.24, respectively. All Pm mono-TPSs had signal peptide sequences at the N-terminus, ranging from 36–48 aa in length, and localized to chloroplasts (Table S2). In contrast, Pm sesqui-TPSs did not contain signal peptide sequences (Table S2). Hydrophobic/hydrophilic analyses showed that all six Pm TPSs were hydrophilic (Figure S1).
As shown in Figure 2a, the mono-TPSs Pm TPS (+)-α-pinene, Pm TPS (−)-α-pinene, Pm TPS (−)-β-pinene, Pm TPS limonene and the sesqui-TPS Pm TPS selinene contain the arginine-tryptophan motif RRX8W, which is conserved among most mono-TPSs and some sesqui-TPSs, near the N-terminus. The RRX8W motif was missing in Pm TPS longifolene. The highly conserved aspartate-rich motif DDXXD and RXR motif were present in these six Pm TPSs. However, another conserved motif, DTE, was not found in P. massoniana. This is crucial for chelating divalent cations, typically Mg2+, in the carboxyl-terminal (C-terminal) domain [50]. Phylogenetic analysis revealed that four monoterpene and two sesquiterpene synthases of P. massoniana were clustered into two subfamilies (TPS-d1 and TPS-d2) with monoterpene and sesquiterpene synthases from other pine plants, except for P. Abies TPS E (E-α-farnesene), P. engelmannii × P. glauca TPS E (E-α-farnesene synthase/E-β-ocimene), and P. teada TPS (α-farnesense), which formed into one clade as monoterpene synthases, in line with previous studies [25,29,51,52] (Figure 2b). Among them, Pm TPS (+)-α-pinene grouped with α-pinene of P. kesiya. Furthermore, these were grouped with TPS α-pinene of P. contorta and P. banksiana. The Pm TPS limonene was grouped with P. sitchensis TPS (−)-limonene and P. abies TPS (−)-limonene. The Pm TPS (−)-β-pinene clustered with TPS (−)-β-pinene and TPS α-terpineol from several Pinus spp. Two sesqui-TPSs, Pm TPS longifolene and P. sylvestris TPS longifolene, were grouped, whereas Pm TPS selinene was grouped with P. sylvestris TPS caryophyllene/humulene (cary/humu), and slightly distant from the A. grandies δ-selinene synthase, annotated as monofunctional but heterologous.
The motif types and permutations of P. massoniana (Figure 2c) and its homologous TPSs (Table S3) were analyzed using the online software MEME (Version 5.5.5). The number and species of the 20 predicted motifs ranged from 11 to 50 aa in length. All Pm mono-TPSs contained 20 motifs, except for Pm TPS (−)-β-pinene. Thirteen motifs (motifs 1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 12, 14, and 17) were found in all Pm TPSs. The conserved domains DDXXD, RRX8W, and RXR were found in motifs 1, 11, and 3, respectively. Motifs 15, 18, 19, and 20 were found only in the Pm mono-TPSs but not in the Pm sesqui-TPSs. Most of the Pm mono-TPSs contained all 20 motif structures in the order 18-11-12-20-7-13-8-19-9-5-17-3-15-1-6-4-2-14-10-16. Motifs 1, 5, 7, 8, 11, and 12 appeared in all sequences. None of the Pm sesqui-TPSs contained motifs 18 or 20. The motifs of TPS D1 and TPS D2 had their characteristics. In particular, the clade of four farnesene and three bisabolene synthases had a long, non-conservative region between motifs 12 and 7, which was not found at the corresponding sites in other clades (Figure S2).

3.2. Prediction of the Secondary Structure of Pm TPSs

The secondary structure of Pm TPS amino acid sequences was predicted using SOPMA software (Figure 3a). The main secondary structures of the six Pm TPSs were α-helices, accounting for 63.27%, 63.91%, 63.04%, 63.62%, 68.95%, and 68.28%, respectively (Table S4). The α-helix content of sesqui-TPSs was slightly higher than that of mono-TPSs. The random coil contents were 26.66%, 26.22%, 29.25%, 30.56%, 29.45%, 24.01%, and 23.79%, respectively. The content of β-sheet/corner (2.92%–4.14%) and extension chain (3.18%–5.52%) was the lowest. The β-sheet/turn and extended strands are interspersed in the α-helix structure, and random coils connect all the structures. Owing to the presence of a signal peptide, the α-helix content in the first 100 amino acids of Pm mono-TPSs was less than that of Pm sesqui-TPSs. All Pm mono-TPSs had a higher content of extended strands than Pm sesqui-TPSs, except for Pm TPS (−)-β-pinene.
SMART software was used to analyze the functional domains of Pm TPSs (Figure 3b). Two Pfam domains, terpene-synth and Terpene-synth-C, have been characterized. The Terpene-synth domain of Pm mono-TPSs starts from the N-terminal 73–75 aa region and ends at the 264–270 aa region, while the terpene-synth-C domain starts from the N-terminal 302–312 aa region and ends at the 492–577 aa region. In Pm sesqui-TPSs, the terpene-synth domain starts from 18/27 and ends at 211, while the terpene-synth-C domain starts from 247/257 and ends at 512/523. In addition, there is a low-complexity region in the N-terminus of Pm TPS (+)-α-pinene, as shown in purple. The sequences were SIVPSMSMSSTTSVS and SMSMSSTTSVS, both serine-rich.

3.3. Homology Modeling and Molecular Docking of Pm TPSs

The predicted Pm TPS structures, as shown in Figure 4, revealed that Pm TPSs comprised two structural domains, the N- and C-terminal, with each domain consisting of several α-helices, β-sheets/turns, and random coils/loops. The C-terminal domain forms 9–17 α-helices from helix A to helix K. The structures of helix C, helix D2, and helix H1 in Pm TPS (+)-α-pinene, Pm TPS (−)-β-pinene, and Pm TPS limonene had 2–3 fewer amino acid residues than those in Pm TPS (−)-α-pinene, Pm TPS selinene, and Pm TPS longifolene. Helix J is presented as a complete long helix (Pm TPS (+)-α-pinene, Pm TPS (−)-β-pinene, and Pm TPS limonene) or equally divided into two short helices (Pm TPS (−)-α-pinene, Pm TPS selinene, and Pm TPS longifolene), according to the template. Other α-helices were conserved in all Pm TPSs.
The results of Ramachandran plot evaluation (Figure S3) revealed that the structural score of Pm TPSs obtained via this modeling were of a high-quality model: 90.7%–94.8% of residues were found in the most favored region, 4.1%–8.2% of residues were found in additionally allowed regions, and 0.4%–1.5% of residues were found in generously allowed and disallowed regions. This justified the stability without any steric hindrance and reliance on the predicted model with acceptable qualities. All TPSs formed active pockets at the C-terminus, where GPP and FPP molecules form hydrogen bonds with different amino acid residues. The docking experiment showed that GPP and FPP could bind to monoterpene and sesquiterpene synthases, respectively. The confidence scores of the six terpene synthases were 4, 4, 4, 4, 3, and 2 for Pm TPS (+)-α-pinene, Pm TPS (−)-α-pinene, Pm TPS (−)-β-pinene, Pm TPS limonene, Pm TPS selinene, and Pm TPS longifolene, respectively. The total scores were 7.23, 7.24, 6.31, 7.68, 7.32, and 5.98 for Pm TPS (+)-α-pinene, Pm TPS (−)-α-pinene, Pm TPS (−)-β-pinene, Pm TPS limonene, Pm TPS selinene, and Pm TPS longifolene, respectively.
Docking simulation results are shown in Figure 5. The side chains of aspartic acid (Asp) 484 and tyrosine (Tyr) 555 formed three hydrogen bonds and one hydrogen bond with the diphosphate group of GPP, respectively, in the substrate-binding pocket in Pm TPS (+)-α-pinene. The amino acids Asp 484 and Tyr 555 are distributed in helices H2 and J, respectively. For Pm TPS (−)-α-pinene, four residues in the binding pocket form hydrogen bonds with GPP: the amino acid glutamine (Glu) 410 in the middle of helix F, the amino acid Glu 489 in the H2-α1 loop, the amino acid serine (Ser) 432 in the G1-G2 loop, and the amino acid arginine (Arg) 473 in the middle of helix H2. The amino acid Arg 473, Glu 410, Glu 489, and Ser 432 formed three, two, two, and one hydrogen bonds with the diphosphate group of GPP, respectively. For Pm TPS (−)-β-pinene, only hydrogen bonds were formed between the side chain of the amino acid threonine (Thr) 559, which is distributed in helix K, and the diphosphate group of GPP. For Pm TPS limonene, five residues in the binding pocket formed hydrogen bonds with GPP: cysteine (Cys) 343 and Asn 346 in helix D, the amino acid Tyr 319 in helix C, the amino acid Tyr 425 in helix F, and the amino acid Ser 451 near the end of helix G1. Each of the above amino acids formed one hydrogen bond with the diphosphate group of GPP.
The side chain of two closely related amino acids, Asp 321 and Asp 325 in helix D, formed a hydrogen bond with the diphosphate group of GPP in the substrate-binding pocket of Pm TPS selinene. For Pm TPS longifolene, four residues in the binding pocket formed hydrogen bonds with FPP: Asp 331 in helix D, the amino acid Ser 433 near the end of helix G1, and Arg 473 and lysine (Lys) 484 in helix H2. The amino acid Arg 473, Asp 331, Ser 433, and Lys 484 formed three, two, one, and one hydrogen bonds with the diphosphate group of FPP, respectively.

3.4. Heterologous Expression and Functional Verification of Pm Mono-TPSs

All cDNAs of Pm TPSs cloned into the pET32 (+) vector with thioredoxin and histidine (HIS) tags were successfully expressed in E. coli DH5α competent cells, and recombinant proteins were then purified with His 60 Ni Superflow Resin and Gravity Columns (Figure S4). Following an in vitro enzyme activity assay, Pm TPS (+)-α-pinene was confirmed to encode a mono-TPS enzyme that catalyzes the formation of mostly (+)-α-pinene, β-pinene, linalool, and a few minor R-(+)-limonene, as detected by GC–FID and GC–MS analysis (Figure 6, Figures S5 and S6, Table S5). Pm TPS (−)-α-pinene catalyzed the formation of mostly (−)-α-pinene, as well as myrcene, R-(+)-limonene, and linalool. Pm TPS β-pinene catalyzed the formation of mostly β-pinene and linalool, as well as a few (+)-α-pinene and R-(+)-limonene. Pm TPS limonene catalyzed the formation of mostly R-(+)-limonene and a few (−)-α-pinene and linalool. The GC–MS analysis of the reaction products catalyzed by Pm TPS (+)-α-pinene with GPP as a substrate showed that the reaction product concentration of (−)-α-pinene was 120 μg/mL after 1 h of reaction, whereas that of Pm TPS limonene, Pm TPS (+)-α-pinene, and Pm TPS β-pinene were 6, 1, and 1.36 μg/mL, respectively (Table S5). With the extension of the incubation time, the main product of Pm mono-TPSs varied and increased during 1–3 h of the in vitro reaction, except for Pm TPS (−)-β-pinene. Within 3–5 h of incubation, the concentrations of the main products of Pm TPS (−)-α-pinene decreased from 125.6 to 114 μg/mL, whereas the concentrations of Pm TPS (+)-α-pinene and Pm TPS limonene increased from 2.43 and 6.67 μg/mL to 2.6 and 8.57 μg/mL. The by-products of Pm TPS (+)-α-pinene, linalool, and (−)-β-pinene were detected after 1 h and 3 h of incubation, respectively, and gradually increased to 2.333 and 0.393 μg/mL after 5 h of incubation. Limonene was detected after 5 h of incubation. Linalool and (−)-α-pinene, the by-products of Pm TPS (−)-β-pinene, were detected after 1 h of incubation, but the content of linalool exceeded that of the main product β-pinene after 3 h of incubation. The content of (−)-α-pinene gradually increased to 0.67 μg/mL after 5 h of incubation but was lower than that of the main product (−)-β-pinene at 1.03 μg/mL. Limonene was detected after incubation for 5 h. The levels of limonene and (+) -α-pinene, by-products of Pm TPS (−)-α-pinene, accounted for approximately 5% of the main product and were detected after 1 h of incubation. The levels of (−)-α-pinene and linalool, two by-products of Pm TPS limonene, were lower than those of Pm TPS (−)-α-pinene by-products, and were detected after 1 and 3 h of incubation, respectively.

4. Discussion

In this study, we identified six TPS genes in P. massoniana, including four mono-TPS and two sesqui-TPS genes. We believe that these genes are involved in the response of P. massoniana to M. alternatus invasion. In vitro enzyme activity assay also verified that Pm mono-TPSs can catalyze substrate GPP to produce volatile terpenes. In our previous study [18,38,39], we proved that exogenous MeJA treating and M. alternatus adult feeding could induce a resin-based defense of P. massoniana, including secondary resin duct formation, induced resin terpenoid accumulation, and gene transcriptional changes related to terpenoid synthesis in metabolic pathways (Figure 7). TPSs catalyze complex carbocation cascade reactions on the prenyl diphosphate substrate, resulting in cyclic or linear terpene backbones [53], which determines the formation and composition of oleoresin terpenes. Therefore, studying the relationship between TPS structure, action site and production is helpful to elucidate the mechanism behind the synthesis of resin with complex composition, as well as the understanding, development and application of terpene synthase in vitro terpene biosynthesis.
The total length of Pm TPSs was similar to that reported in Pinaceae [54]. Each Pm mono-TPSs contained a peptide sequence of approximately 50 amino acid residues, similar to SaTPS1 identified from Santalum album [55]. Studies have shown that a lack of the N-terminal peptide FaNES1 leads to changes in subcellular localization, resulting in different products [56]. The Pm sesqui-TPSs had no peptide sequence, suggesting that it may function in the cytoplasm, whereas the sesqui-TPSs from Solanum habrochaites (santalene/bergamotene synthases) are located in chloroplasts, indicating that sesquiterpenes can also be synthesized in plastids [57]. Among the three characteristically conserved domains of TPS, RxR and DDxxD were found in all the sequences. The RRx8W domain was found not only in monoterpenes but also in Pm TPS selinene. Monoterpene synthases lacking the RRx8W motif are responsible for the production of acyclic mono-TPSs, a characteristic of the TPS-g subfamily [58]. Previous studies on the functional verification of menthol limonene synthase have shown that the tandem arginine sequence RR in RRx8W might catalyze the transformation of GPP to an intermediate product [59], suggesting that Pm TPS selinene can bind the substrate GPP. Additionally, in A. grandis, RRx8W was also found in the N-terminus of δ-selinene synthase, and its heterologous expression could produce 52 sesquiterpenes [32]. The conserved domain NSE/DTE was not found in Pm TPSs nor the reported TPS sequences of P. taeda, P. contorta, P. banksiana, and P. nigra subsp. laricio [22,34,60], suggesting that the conserved domain NSE/DTE may not be conserved in Pinus.
In the phylogenetic tree, we found that the amino acid sequences of TPSs in congeners were more similar to those of TPSs with the same main products in species from different genera (Figure 2b). For example, the distance between Pm TPS seli (Pm TPS selinene) and A. grandis TPS d-seli, both of which produce selinene, was greater than the distance between Pm TPS seli and Pinus sylvestris TPS cary/humu, which are congeners but have different main products. This feature is also found in other species. The similarity between caryophyllene synthase from Oryza sativa, Arabidopsis thaliana, Artemisia annua, and Cucumis sativus is less than 40%. In contrast, caryophyllene synthase from Oryza sativa has 40%–42% similarity with the synthases that produce other sesquiterpenes from Zea mays [61]. No special motifs were found in the sesquiterpene sequences used to construct the phylogenetic tree. None of the Pm sesqui-TPSs contained motifs 18 or 20. Picea abies TPS E, E-alpha-far, Pinus taeda TPS alpha-far, and Picea engelmannii × Picea glauca TPS (E, E)-alpha-far/(E)-beta-oci in TPS D1 contained an additional motif 15 compared to the other sesqui-TPSs. The absence of these motifs suggests a potential difference in the evolution of mono- and sesqui-TPSs. Furthermore, in the phylogenetic tree, motif 19 was absent in all sesqui-TPS sequences of Pinus, including P. massoniana. We hypothesize that this may provide evidence for the functional differentiation of TPSs occurring before species differentiation.
Previous study has shown that the catalytic function of TPSs is affected by kinetics rather than thermodynamics, meaning that the types of intermediates initially formed in the reaction can affect the formation of subsequent products [62]. The formation of multiple centromeres in the active center can generate a variety of substrates, leading to the formation of by-products in fewer reaction steps [63]. Molecular docking simulations showed that the substrate-binding pockets of sesqui-TPSs were larger than those of mono-TPSs. The substrate of sesqui-TPS is FPP (C15), which contains three double bonds and is more complex than GPP (C10), which contains two double bonds [64]. Because of the increase in the chain length of the substrate, the number of sesquiterpene intermediates is much greater than that of monoterpene intermediates, which enriches the diversity of product structures [65]. Larger pockets may facilitate the entry of substrates and the formation of products. Secondary metabolic enzymes tend to produce mechanistic elasticity during evolution, giving rise to new enzymes with a wider range of substrates and a variety of new products [66].
Some residues are responsible for catalytic elasticity within and around the active sites of TPSs. Changing these residues can regulate the products or catalytic rates [67,68,69]. For example, Xu Jin-Kun et al. [70] reported that Tyr residues (Y573) in S-limonene synthase may be involved in the earlier steps of the reaction, probably by controlling the conformation of the helical linalyl diphosphate intermediate. The replacement of Ile in A. thaliana ent-kaurene synthase with Thr can convert it into a pimaradiene synthase, but the catalytic efficiency is reduced by four times [71]. The replacement of nine amino acids in Nicotiana tabacum can accomplish conversion between the catalytic specificities of 5-epi-aristolochene synthase and premnaspirodiene synthase [68]. Molecular docking between Pm TPSs and the corresponding substrates showed that there were seven kinds of residues in Pm mono-TPSs participating in the catalytic process, Asp, Arg, Tyr, Glu, Trp, Thr, and Cys, which were distributed between helix C and helix K. Asp, Try, Arg, and Lys were found in Pm sesqui-TPSs, distributed between helix D and helix H2, similar to reports from A. grandis [31]. Mn2+ is a common catalyst in biological catalytic reactions, often combining with Asp [72], which may be a reason for the elevated catalytic efficiency of some TPSs in the presence of Mn2+ [73]. Thus, adding Mn2+ may accelerate the production of Pm (+)-α-pinene, Pm selinene, and Pm longifolene. At a pH ≥ 7, Asp and Glu are deprotonated, facilitating the formation of a carbonium ion intermediate. Therefore, when using Pm mono-TPSs in the catalytic reaction, it is critical to maintain a pH above 7 for stable production. The surfaces of the aromatic groups in Tyr and Trp are positively charged, which can stabilize the carbonium ion intermediate [62]. Because Pm sesqui-TPSs have fewer aromatic groups than Pm mono-TPSs, there may be more binding possibilities for the structures of the intermediates and products in Pm sesqui-TPSs than that in Pm mono-TPSs.
During expression vector construction, a purified tag was inserted at the N-terminus to avoid the influence of redundant tag sequences at the C-terminal active site during protein activity tests. The non-removal of peptide sequences in this study did not affect protein expression and in vitro activity tests, as in other studies [22,61,74,75,76,77]. Previous studies have shown that 1–1.5 h is the optimal reaction duration for collecting enzymatic reaction products [22,34]. In this study, most of the main products were produced within 1 h, whereas longer reaction durations only significantly increased the number of by-products. Therefore, under a consistent metal ion concentration and incubation temperature, it is recommended to complete the reaction of TPSs within 1–1.5 h. The products of mono-TPSs mainly possess chirality, with several trace products that have a related stereo-structure [78]. Heterologous expression of monoterpene synthases reported in Pinus is usually named after their main products, such as (+)-α-pinene synthase and (−)-α-pinene synthase from P. taeda, (+)-α-pinene synthases, (−)-α-pinene synthases, and (−)-β-pinene synthase from P. contorta and P. banksiana, as well as (−)-α-pinene synthase and β-pinene synthase from Pinus elliottii [22,34,79]. In this study, the main product profiles of the Pm TPSs were similar to those reported above but at different concentrations. Meanwhile, as observed in the TPSs of P. taeda [22], most Pm TPSs only generated certain enantiomers.
Four Pm mono-TPSs were heterologously expressed. Among them, Pm (−)-α-pinene produced the largest number of products and formed the largest number of hydrogen bonds with GPP. The varied amounts of by-product linalool produced by Pm mono-TPSs suggest that Pm (−)-α-pinene may have a higher catalytic specificity than Pm (+)-α-pinene and Pm (−)-β-pinene. Therefore, the number of active sites and hydrogen bonds between the substrate and enzyme may be positively related to product specificity. Chen Rui-Xu et al. [40] analyzed the volatiles produced by P. massoniana fed on by M. alternatus adults. The ratio of β- pinenes to α-pinenes in the volatiles was approximately 1:100, which is very close to the ratio of these two products in in vitro activity tests. However, the chemical composition of the products in the in vitro activity tests differed from that of the volatiles. This may be due to the catalytic functions of downstream enzymes, such as P450 enzymes, which in plants catalyze the transformation of some TPS products to other compounds [80]. For example, limonene in Mentha canadensis can be converted into (−)-carvone or (−)-isomentadiene by two P450 enzymes [81]. Notably, although the limonene synthase gene was upregulated in P. massoniana fed upon by M. alternatus adults, and limonene was the main heterologous expression product, this compound was not detected in volatiles. The limonene produced may have been converted to pinocarvone by downstream P450 enzymes [40].

5. Conclusions

Terpene synthases catalyze complex multistep reactions that generate hundreds of structurally diverse terpenoids of biological and commercial importance. In a previous study, we showed how defense mechanisms are triggered in P. massoniana through a complex series of changes, including traumatic resin duct formation, which induce terpenoid accumulation, key enzyme activation, and TPS gene expression [18,38]. In this study, we cloned six TPS genes from P. massoniana involved in these protective responses. We then subjected these genes to bioinformatic analysis and in vitro activity tests. Phylogenetic analysis indicated that the evolutionary distance between TPSs of the same species was greater than that between the same products within P. massoniana. Furthermore, we uncovered the structural features that may be responsible for the functional attributes of TPSs through molecular docking, verifying our findings via in vitro activity tests of Pm TPSs. Studying the catalytic mechanism of TPS could aid in understanding the structure–function relationship of TPSs in P. massoniana and provide a theoretical basis for resistant tree breeding through genetic modification.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/f15020244/s1, Figure S1. Hydrophobicity or hydrophilicity prediction of six amino acid sequences of Pm TPSs. Figure S2. Motif analysis of 91 terpene synthase sequences. Figure S3. Ramachandran plot of P. massoniana terpene synthases. Figure S4. SDS–PAGE of Pm mono-TPSs. Figure S5. The standard of four monoterpenes identified by GC–MS. Figure S6. Separation between the peaks of (+)-α-pinene and (−)-α-pinene in both authentic standards and the catalyzed productions. Table S1. Coding sequences of 85 monoterpene and sesquiterpene synthases from pine trees. Table S2. Features of TPS encoding sequences in Pinus massoniana. Table S3. Consensus sequence identified by MEME in TPSs of Pinus massoniana. Table S4. The secondary features of TPS encoding sequences in Pinus massoniana. Table S5. Terpene production content of Pm mono-TPSs.

Author Contributions

Methodology, R.C.; Formal analysis, Q.W.; Data curation, Q.W. and R.C.; Writing—original draft, Q.W.; Writing—review and editing, R.C., T.X. and D.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 32001322), the Natural Science Foundation for Colleges and Universities in Jiangsu Province (No. 20KJB220004), and the Science Foundation of Jiangsu Vocational college of Agriculture and Forestry (Nos. 2020kj002 and 2021kj90).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All relevant data is contained within the article. The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding authors.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

terpene synthase (TPS)
Pinus massoniana(Pm)
geranyl diphosphate (GPP)
farnesyl diphosphate (FPP)
gas chromatography-flame ionization detection (GC/FID)
gas chromatography-mass spectrometry (GC/MS)
protein data bank (PDB)
aspartic acid (Asp)
arginine (Arg)
tyrosine (Tyr)
glutamic acid (Glu)
tryptophan (Trp)
threonine (Thr)
cysteine (Cys)
lysine (Lys)
serin (Ser)
histidine (HIS)

References

  1. Tomlin, E.S.; Antonejevic, E.; Alfaro, R.I.; Borden, J.H. Changes in volatile terpene and diterpene resin acid composition of resistant and susceptible white spruce leaders exposed to simulated white pine weevil damage. Tree Physiol. 2000, 20, 1087–1095. [Google Scholar] [CrossRef]
  2. Trapp, S.; Croteau, R. Defensive resin biosynthesis in conifers. Annu. Rev. Plant Physiol. Plant Mol. Biol. 2001, 52, 689–724. [Google Scholar] [CrossRef]
  3. Bohlmann, J. Pine terpenoid defences in the mountain pine beetle epidemic and in other conifer pest interactions: Specialized enemies are eating holes into a diverse, dynamic and durable defence system. Tree Physiol. 2012, 32, 943–945. [Google Scholar] [CrossRef]
  4. Mumm, R.; Hilker, M. Direct and indirect chemical defence of pine against folivorous insects. Trends Plant Sci. 2006, 11, 351–358. [Google Scholar] [CrossRef]
  5. Zulak, K.G.; Bohlmann, J. Terpenoid Biosynthesis and Specialized Vascular Cells of Conifer Defense. J. Integr. Plant Biol. 2010, 52, 86–97. [Google Scholar] [CrossRef] [PubMed]
  6. Celedon, J.M.; Bohlmann, J. Oleoresin defenses in conifers: Chemical diversity, terpene synthases and limitations of oleoresin defense under climate change. New Phytol. 2019, 224, 1444–1463. [Google Scholar] [CrossRef]
  7. Miller, B.; Madilao, L.L.; Ralph, S.; Bohlmann, J. Insect-induced conifer defense. White pine weevil and methyl jasmonate induce traumatic resinosis, de novo formed volatile emissions, and accumulation of terpenoid synthase and putative octadecanoid pathway transcripts in Sitka spruce. Plant Physiol. 2005, 137, 369–382. [Google Scholar] [CrossRef]
  8. Schnee, C.; Köllner, T.G.; Held, M.; Turlings, T.C.; Gershenzon, J.; Degenhardt, J. The products of a single maize sesquiterpene synthase form a volatile defense signal that attracts natural enemies of maize herbivores. Proc. Natl. Acad. Sci. USA 2006, 103, 1129–1134. [Google Scholar] [CrossRef]
  9. Faiola, C.L.; Jobson, B.T.; VanReken, T.M. Impacts of simulated herbivory on volatile organic compound emission profiles from coniferous plants. Biogeosciences 2015, 12, 527–547. [Google Scholar] [CrossRef]
  10. Lundborg, L.; Fedderwitz, F.; Björklund, N.; Nordlander, G.; Borg-Karlson, A.K. Induced defenses change the chemical composition of pine seedlings and influence meal properties of the pine weevil Hylobius abietis. Phytochemistry 2016, 130, 99–105. [Google Scholar] [CrossRef] [PubMed]
  11. Munro, H.L.; Gandhi, K.J.; Barnes, B.F.; Montes, C.R.; Nowak, J.T.; Shepherd, W.P.; Villari, C.; Sullivan, B.T. Electrophysiological and behavioral responses Dendroctonus frontalis and D. terebrans (Coleoptera: Curculionidae) to resin odors of host pines (Pinus spp.). Chemoecology 2020, 30, 215–231. [Google Scholar] [CrossRef]
  12. Martini, X.; Hughes, M.A.; Smith, J.A.; Stelinski, L.L. Attraction of redbay ambrosia beetle, Xyleborus Glabratus, to leaf volatiles of its host plants in North America. J. Chem. Ecol. 2015, 41, 613–621. [Google Scholar] [CrossRef]
  13. Foti, V.; Araniti, F.; Manti, F.; Alicandri, E.; Giuffrè, A.M.; Bonsignore, C.P.; Castiglione, E.; Sorgonà, A.; Covino, S.; Paolacci, A.R.; et al. Profiling volatile terpenoids from calabrian pine stands infested by the pine processionary moth. Plants 2020, 9, 1362. [Google Scholar] [CrossRef]
  14. Huang, K.; Chen, R.; Chen, H.; Zhu, H.; Hao, D. Methyl jasmonate induced traumatic resin ducts differentiation in the sapling of Pinus masssoniana Lamb. J. Northeast. For. Univ. 2019, 47, 119–124+128. (In Chinese) [Google Scholar] [CrossRef]
  15. Fedderwitz, F.; Nordlander, G.; Ninkovic, V.; Björklund, N. Effects of jasmonate-induced resistance in conifer plants on the feeding behaviour of a bark-chewing insect, Hylobius abietis. J. Pest. Sci. 2016, 89, 97–105. [Google Scholar] [CrossRef]
  16. Fäldt, J.; Martin, D.; Miller, B.; Rawat, S.; Bohlmann, J. Traumatic resin defense in Norway spruce (Picea abies): Methyl jasmonate-induced terpene synthase gene expression, and cDNA cloning and functional characterization of (+)-3-carene synthase. Plant Mol. Biol. 2003, 51, 119–133. [Google Scholar] [CrossRef]
  17. Huber, D.P.; Philippe, R.N.; Godard, K.A.; Sturrock, R.N.; Bohlmann, J. Characterization of four terpene synthase cDNAs from methyl jasmonate-induced Douglas-fir, Pseudotsuga menziesii. Phytochemistry 2005, 66, 1427–1439. [Google Scholar] [CrossRef]
  18. Chen, R.; Li, Y.; Huang, K.; Li, H.; Chen, H.; Xu, T.; Hao, D. Time-course transcriptomic study of phenolic metabolism and P450 enzymes in Pinus massoniana Lamb. after feeding by Monochamus alternatus Hope. Scand. J. For. Res. 2019, 34, 569–576. [Google Scholar] [CrossRef]
  19. Lin, J.; Wang, D.; Chen, X.; Köllner, T.; Mazarei, M.; Guo, H.; Pantalone, V.R.; Arelli, P.; Stewart, C.N., Jr.; Wang, N.; et al. An (E, E)-α-farnesene synthase gene of soybean has a role in defence against nematodes and is involved in synthesizing insect-induced volatiles. Plant Biotechnol. J. 2017, 15, 510–519. [Google Scholar] [CrossRef] [PubMed]
  20. Bleeker, P.M.; Mirabella, R.; Diergaarde, P.J.; VanDoorn, A.; Tissier, A.; Kant, M.R.; Prins, M.; de Vos, M.; Haring, M.A.; Schuurink, R.C. Improved herbivore resistance in cultivated tomato with the sesquiterpene biosynthetic pathway from a wild relative. Proc. Natl. Acad. Sci. USA 2012, 109, 20124–20129. [Google Scholar] [CrossRef]
  21. Keeling, C.I.; Bohlmann, J. Genes, enzymes and chemicals of terpenoid diversity in the constitutive and induced defence of conifers against insects and pathogens. New Phytol. 2006, 170, 657–675. [Google Scholar] [CrossRef] [PubMed]
  22. Phillips, M.A.; Wildung, M.R.; Williams, D.C.; Hyatt, D.C.; Croteau, R. cDNA isolation, functional expression, and characterization of (+)-alpha-pinene synthase and (−)-alpha-pinene synthase from loblolly pine (Pinus taeda): Stereocontrol in pinene biosynthesis. Arch. Biochem. Biophys. 2003, 411, 267–276. [Google Scholar] [CrossRef] [PubMed]
  23. Martin, D.M.; Fäldt, J.; Bohlmann, J. Functional Characterization of Nine Norway Spruce TPS Genes and Evolution of Gymnosperm Terpene Synthases of the TPS-d Subfamily. Plant Physiol. 2004, 135, 1908–1927. [Google Scholar] [CrossRef] [PubMed]
  24. Bohlmann, J.; Meyer-Gauen, G.; Croteau, R. Plant terpenoid synthases: Molecular biology and phylogenetic analysis. Proc. Natl. Acad. Sci. USA 1998, 95, 4126–4133. [Google Scholar] [CrossRef] [PubMed]
  25. Warren, R.L.; Keeling, C.I.; Yuen, M.M.; Raymond, A.; Taylor, G.A.; Vandervalk, B.P.; Mohamadi, H.; Paulino, D.; Chiu, R.; Jackman, S.D.; et al. Improved white spruce (Picea glauca) genome assemblies and annotation of large gene families of conifer terpenoid and phenolic defense metabolism. Plant J. 2015, 83, 189–212. [Google Scholar] [CrossRef] [PubMed]
  26. Luchi, N.; Ma, R.; Capretti, P.; Bonello, P. Systemic induction of traumatic resin ducts and resin flow in Austrian pine by wounding and inoculation with Sphaeropsis sapinea and Diplodia scrobiculata. Planta 2005, 221, 75–84. [Google Scholar] [CrossRef] [PubMed]
  27. Bonello, P.; Gordon, T.R.; Storer, A.J. Systemic induced resistance in Monterey pine. For. Path. 2001, 31, 99–106. [Google Scholar] [CrossRef]
  28. Chen, F.; Tholl, D.; Bohlmann, J.; Pichersky, E. The family of terpene synthases in plants: A mid-size family of genes for specialized metabolism that is highly diversified throughout the kingdom. Plant J. 2011, 66, 212–229. [Google Scholar] [CrossRef]
  29. Keeling, C.I.; Weisshaar, S.; Ralph, S.G.; Jancsik, S.; Hamberger, B.; Dullat, H.K.; Bohlmann, J. Transcriptome mining, functional characterization, and phylogeny of a large terpene synthase gene family in spruce (Picea spp.). BMC Plant Biol. 2011, 11, 43. [Google Scholar] [CrossRef]
  30. Vogel, B.S.; Wildung, M.R.; Vogel, G.; Croteau, R. Abietadiene synthase from grand fir (Abies grandis). cDNA isolation, characterization, and bacterial expression of a bifunctional diterpene cyclase involved in resin acid biosynthesis. J. Biol. Chem. 1996, 271, 23262–23268. [Google Scholar] [CrossRef]
  31. Bohlmann, J.; Steele, C.L.; Croteau, R. Monoterpene synthases from grand fir (Abies grandis): cDNA isolation, characterization, and functional expression of myrcene synthase, (−)-(4S)-limonene synthase, and (−)-(1S,5S)-pinene synthase. J. Biol. Chem. 1997, 272, 21784–21792. [Google Scholar] [CrossRef]
  32. Steele, C.L.; Crock, J.; Bohlmann, J.; Croteau, R. Sesquiterpene Synthases from Grand Fir (Abies grandis). Comparison of constitutive and wound-induced activities, and cDNA isolation, characterization, and bacterial expression of delta-selinene synthase and gamma-humulene synthase. J. Biol. Chem. 1998, 273, 2078–2089. [Google Scholar] [CrossRef] [PubMed]
  33. Martin, D.M.; Bohlmann, J. Identification of Vitis vinifera (−)-alpha-terpineol synthase by in silico screening of full-length cDNA ESTs and functional characterization of recombinant terpene synthase. Phytochemistry 2004, 65, 1223–1229. [Google Scholar] [CrossRef] [PubMed]
  34. Hall, D.E.; Yuen, M.M.; Jancsik, S.; Quesada, A.L.; Dullat, H.K.; Li, M.; Henderson, H.; Arango-Velez, A.; Liao, N.Y.; Docking, R.T.; et al. Transcriptome resources and functional characterization of monoterpene synthases for two host species of the mountain pine beetle, lodgepole pine (Pinus contorta) and jack pine (Pinus banksiana). BMC Plant Biol. 2013, 13, 80. [Google Scholar] [CrossRef] [PubMed]
  35. Lee, S.; Oh, D.; Lee, S.; Chung, S.; Dong-Soon, K. Subspecific Synonym of Monochamus alternatus (Coleoptera: Cerambycidae): Population Genetics and Morphological Reassessment. J. Econ. Entomol. 2022, 115, 1987–1994. [Google Scholar] [CrossRef] [PubMed]
  36. Sugimoto, H.; Fukuta, M.; Togashi, K. Yearly changes in dispersal and life-history traits of Monochamus alternatus Hope with reference to its outbreak. J. Appl. Entomol. 2020, 144, 459–467. [Google Scholar] [CrossRef]
  37. Linit, M.J. Nemtaode-vector relationships in the pine wilt disease system. J. Nematol. 1988, 20, 227–235. [Google Scholar] [PubMed]
  38. Akbulut, S.; Stamps, W.T. Insect vectors of the pinewood nematode: A review of the biology and ecology of Monochamus species. For. Pathol. 2012, 42, 89–99. [Google Scholar] [CrossRef]
  39. Ye, J. Epidemic status of pine wilt disease in China and its prevention and control techniques and counter measures. Sci. Silvae Sin. 2019, 55, 1–10. (In Chinese) [Google Scholar]
  40. Chen, R.; He, X.; Chen, J.; Gu, T.; Liu, P.; Xu, T.; Teale, S.A.; Hao, D. Traumatic resin duct development, terpenoid formation, and related synthase gene expression in Pinus massoniana under feeding pressure of Monochamus alternatus. J. Plant Growth Regul. 2019, 38, 897–908. [Google Scholar] [CrossRef]
  41. Chen, R.; Huang, K.; Pan, S.; Xu, T.; Tan, J.; Hao, D. Jasmonate induced terpene-based defense in Pinus massoniana depresses Monochamus alternatus adult feeding. Pest. Manag. Sci. 2021, 77, 731–740. [Google Scholar] [CrossRef]
  42. Hao, D.; Ma, F.; Wang, Y.; Dai, H.; Zhang, Y. Electroantennogram and behavioural responses of Monochamus alternatus to volatiles from Pinus massoniana. J. Appl. Entomol. 2007, 44, 541–544. [Google Scholar] [CrossRef]
  43. Zhao, L.; Wei, W.; Kang, L.; Sun, J. Chemotaxis of the pinewood nematode, Bursaphelenchus xylophilus, to volatiles associated with host pine, Pinus massoniana, and its vector Monochamus alternatus. J. Chem. Ecol. 2007, 33, 1207–1216. [Google Scholar] [CrossRef]
  44. Yang, R.; Li, D.; Yi, S.; Wei, Y.; Wang, M. Odorant-binding protein 19 in Monochamus alternatus involved in the recognition of a volatile strongly emitted from ovipositing host pines. Insect Sci. 2023, 1–13. [Google Scholar] [CrossRef]
  45. Chen, X. Studies on Extraction and Expression of Lipid-Producing Genes in Masson’s Pine. Ph.D. Thesis, Chinese Academy of Forestry, Beijing, China, 2018. [Google Scholar]
  46. Liu, B.; Liu, Q.; Zhou, Z.; Yin, H.; Xie, Y.; Wei, Y. Two terpene synthases in resistant Pinus massoniana contribute to defence against Bursaphelenchus xylophilus. Plant Cell Environ. 2021, 44, 257–274. [Google Scholar] [CrossRef] [PubMed]
  47. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef]
  48. Laskowski, R.A.; Macarthur, M.W.; Moss, D.S.; Thornton, J.M. PROCHECK: A program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 1993, 26, 283–291. [Google Scholar] [CrossRef]
  49. Kovats, E. Gas chromatographic characterization of organic substances in the retention index system. Adv. Chromatogr. 1965, 1, 229–247. [Google Scholar]
  50. Marrero, P.F.; Poulter, C.D.; Edwards, P.A. Effects of site-directed mutagenesis of the highly conserved aspartate residues in domain II of farnesyl diphosphate synthase activity. J. Biol. Chem. 1992, 267, 21873–21878. [Google Scholar] [CrossRef]
  51. McCaskill, D.; Croteau, R. Prospects for the bioengineering of isoprenoid biosynthesis. Biochem. Eng. Biotechnol. 1997, 55, 107–146. [Google Scholar] [CrossRef]
  52. Kyte, J.; Doolittle, R.F. A simple method for displaying the hydropathic character of a protein. J. Mol. Biol. 1982, 157, 105–132. [Google Scholar] [CrossRef]
  53. Lanier, E.R.; Andersen, T.B.; Hamberger, B. Plant terpene specialized metabolism: Complex networks or simple linear pathways? Plant J. 2023, 114, 1178–1201. [Google Scholar] [CrossRef]
  54. Alicandri, E.; Paolacci, A.R.; Osadolor, S.; Sorgonà, A.; Badiani, M.; Ciaffi, M. On the evolution and functional diversity of terpene synthases in the Pinus species: A Review. J. Mol. Evol. 2020, 88, 253–283. [Google Scholar] [CrossRef]
  55. Zhang, X.; Niu, M.; Teixeira da Silva, J.A.; Zhang, Y.; Yuan, Y.; Jia, Y.; Xiao, Y.; Li, Y.; Fang, L.; Zeng, S.; et al. Identification and functional characterization of three new terpene synthase genes involved in chemical defense and abiotic stresses in Santalum album. BMC Plant Biol. 2019, 19, 115. [Google Scholar] [CrossRef]
  56. Aharoni, A.; Giri, A.P.; Verstappen, F.W.; Bertea, C.M.; Sevenier, R.; Sun, Z.; Jongsma, M.A.; Schwab, W.; Bouwmeester, H.J. Gain and loss of fruit flavor compounds produced by wild and cultivated strawberry species. Plant Cell 2004, 16, 3110–3131. [Google Scholar] [CrossRef]
  57. Williams, D.C.; McGarvey, D.J.; Katahira, E.J.; Croteau, R. Truncation of limonene synthase preprotein provides a fully active ‘pseudomature’ form of this monoterpene cyclase and reveals the function of the amino-terminal arginine pair. Biochemistry 1998, 37, 12213–12220. [Google Scholar] [CrossRef]
  58. Abbas, F.; Ke, Y.; Yu, R.; Fan, Y. Functional characterization and expression analysis of two terpene synthases involved in floral scent formation in Lili’m ‘Sibeia’. Planta 2019, 249, 71–93. [Google Scholar] [CrossRef]
  59. Spezio, M.; Wilson, D.B.; Karplus, P.A. Crystal structure of the catalytic domain of a thermophilic endocellulase. Biochemistry 1993, 32, 9906–9916. [Google Scholar] [CrossRef] [PubMed]
  60. Alicandri, E.; Covino, S.; Sebastiani, B.; Paolacci, A.R.; Badiani, M.; Sorgonà, A.; Ciaffi, M. Monoterpene Synthase Genes and Monoterpene Profiles in Pinus nigra subsp. laricio. Plants 2022, 11, 449. [Google Scholar] [CrossRef] [PubMed]
  61. Wang, L.; Fang, X.; Yang, C.; Li, J.; Chen, X. Biosynthesis and regulation of secondary terpenoid metabolism in plants. Sci. Sin. Vitae. 2013, 43, 1030–1046. (In Chinese) [Google Scholar] [CrossRef]
  62. Christianson, D.W. Structural biology and chemistry of the terpenoid cyclases. Chem. Rev. 2006, 106, 3412–3442. [Google Scholar] [CrossRef]
  63. Vedula, L.S.; Rynkiewicz, M.J.; Pyun, H.J.; Coates, R.M.; Cane, D.E.; Christianson, D.W. Molecular recognition of the substrate diphosphate group governs product diversity in trichodiene synthase mutants. Biochemistry 2005, 44, 6153–6163. [Google Scholar] [CrossRef] [PubMed]
  64. Durairaj, J.; Di Girolamo, A.; Bouwmeester, H.J.; de Ridder, D.; Beekwilder, J.; van Dijk, A.D. An analysis of characterized plant sesquiterpene synthases. Phytochemistry 2019, 158, 157–165. [Google Scholar] [CrossRef] [PubMed]
  65. Tian, B.; Poulter, C.D.; Jacobson, M.P. Defining the product chemical space of monoterpenoid synthases. PLoS Comput. Biol. 2016, 12, e1005053. [Google Scholar] [CrossRef] [PubMed]
  66. Weng, J.K.; Philippe, R.N.; Noel, J.P. The rise of chemodiversity in plants. Science 2012, 336, 1667–1670. [Google Scholar] [CrossRef] [PubMed]
  67. Starks, C.M.; Back, K.; Chappell, J.; Noel, J.P. Structural basis for cyclic terpene biosynthesis by tobacco 5-epi-aristolochene synthase. Science 1997, 277, 1815–1820. [Google Scholar] [CrossRef] [PubMed]
  68. Greenhagen, B.T.; O’Maille, P.E.; Noel, J.P.; Chappell, J. Identifying and manipulating structural determinates linking catalytic specificities in terpene synthases. Proc. Natl. Acad. Sci. USA 2006, 103, 9826–9831. [Google Scholar] [CrossRef] [PubMed]
  69. Faylo, J.L.; van Eeuwen, T.; Kim, H.J.; Gorbea Colón, J.J.; Garcia, B.A.; Murakami, K.; Christianson, D.W. Structural insight on assembly-line catalysis in terpene biosynthesis. Nat. Commun. 2021, 12, 3487. [Google Scholar] [CrossRef] [PubMed]
  70. Xu, J.; Xu, J.; Ai, Y.; Farid, R.A.; Tong, L.; Yang, D. Mutational analysis and dynamic simulation of S-limonene synthase reveal the importance of Y573: Insight into the cyclization mechanism in monoterpene synthases. Arch. Biochem. Biophys. 2018, 638, 27–34. [Google Scholar] [CrossRef]
  71. Xu, M.; Wilderman, P.R.; Peters, R.J. Following evolution’s lead to a single residue switch for diterpene synthase product outcome. Proc. Natl. Acad. Sci. USA 2007, 104, 7397–7401. [Google Scholar] [CrossRef]
  72. Qin, T. Amino acid residues in protein structure and function. Life Sci. Instrum. 2011, 9, 41–43. (In Chinese) [Google Scholar]
  73. Trindade, H.; Sena, I.; Figueiredo, A.C. Characterization of α-pinene synthase gene in Pinus pinaster and P. pinea in vitro cultures and differential gene expression following Bursaphelenchus xylophilus inoculation. Acta Physiol. Plant. 2016, 38, 143. [Google Scholar] [CrossRef]
  74. Bohlmann, J.; Phillips, M.; Ramachandiran, V.; Katoh, S.; Croteau, R. cDNA cloning, characterization, and functional expression of four new monoterpene synthase members of the Tpsd gene family from grand fir (Abies grandis). Arch. Biochem. Biophys. 1999, 368, 232–243. [Google Scholar] [CrossRef] [PubMed]
  75. Williams, D.C.; Wildung, M.R.; Jin, A.Q.; Dalal, D.; Oliver, J.S.; Coates, R.M.; Croteau, R. Heterologous expression and characterization of a “Pseudomature” form of taxadiene synthase involved in paclitaxel (Taxol) biosynthesis and evaluation of a potential intermediate and inhibitors of the multistep diterpene cyclization reaction. Arch. Biochem. Biophys. 2000, 379, 137–146. [Google Scholar] [CrossRef] [PubMed]
  76. Peters, R.J.; Flory, J.E.; Jetter, R.; Ravn, M.M.; Lee, H.J.; Coates, R.M.; Croteau, R.B. Abietadiene synthase from grand fir (Abies grandis): Characterization and mechanism of action of the “pseudomature” recombinant enzyme. Biochemistry 2000, 39, 15592–15602. [Google Scholar] [CrossRef]
  77. Borghi, M.; Xie, D.Y. Cloning and characterization of a monoterpene synthase gene from flowers of Camelina sativa. Planta 2018, 247, 443–457. [Google Scholar] [CrossRef]
  78. Christianson, D.W. Structural and chemical biology of terpenoid cyclases. Chem. Rev. 2017, 117, 11570–11648. [Google Scholar] [CrossRef]
  79. Diao, S.; Zhang, Y.; Luan, Q.; Ding, X.; Sun, J.; Jiang, J. Identification of TPS-d subfamily genes and functional characterization of three monoterpene synthases in slash pine. Ind. Crops Prod. 2022, 188, 115609. [Google Scholar] [CrossRef]
  80. Li, J.; Luo, X.; Zhao, P.; Zeng, Y. Post-modification enzymes involved in the biosynthesis of plant terpenoids. Plant Divers. 2009, 31, 461–468. [Google Scholar] [CrossRef]
  81. Ringer, K.L.; Davis, E.M.; Croteau, R. Monoterpene Metabolism. Cloning, Expression, and Characterization of (−)-Isopiperitenol/(−)-Carveol Dehydrogenase of Peppermint and Spearmint. Plant Physiol. 2005, 137, 863–872. [Google Scholar] [CrossRef]
Figure 1. Egg, larva, pupa, and male adult of M. alternatus.
Figure 1. Egg, larva, pupa, and male adult of M. alternatus.
Forests 15 00244 g001
Figure 2. (a) Amino acid sequence alignment of six terpene synthases (TPSs), Pm TPS (+)-α-pinene, Pm TPS (−)-α-pinene, Pm TPS (−)-β-pinene, Pm TPS limonene, Pm TPS selinene, and Pm TPS longifolene. The nucleotide sequences were aligned using DNAMAN 6.0. The conserved terpene synthase motifs, RRX8W, RXR, DDXXD, are underlined. Amino acid residues are shaded in light gray, gray, and black, showing 50, 75 and 100% identity, respectively, while dashes represent the gaps used for optimal alignment. (b) Phylogenetic tree of six Pm TPSs together with TPSs of other Pinus spp. The alignment of nucleotide sequences was performed using Clustal X 2.1, and the tree was constructed using the maximum likelihood method, with 500 replicates for bootstrapping, using MEGA 7 software (Version 7.0.26). All Pm TPSs are classified into TPS-d1and TPS-d2 subfamilies based on the taxonomic distribution of TPS families. The accession numbers of all genes are given in Supplementary Table S1. Missing data and positions with gaps were eliminated. (c) Motif analysis of mono- and sesqui-TPSs in P. massoniana.
Figure 2. (a) Amino acid sequence alignment of six terpene synthases (TPSs), Pm TPS (+)-α-pinene, Pm TPS (−)-α-pinene, Pm TPS (−)-β-pinene, Pm TPS limonene, Pm TPS selinene, and Pm TPS longifolene. The nucleotide sequences were aligned using DNAMAN 6.0. The conserved terpene synthase motifs, RRX8W, RXR, DDXXD, are underlined. Amino acid residues are shaded in light gray, gray, and black, showing 50, 75 and 100% identity, respectively, while dashes represent the gaps used for optimal alignment. (b) Phylogenetic tree of six Pm TPSs together with TPSs of other Pinus spp. The alignment of nucleotide sequences was performed using Clustal X 2.1, and the tree was constructed using the maximum likelihood method, with 500 replicates for bootstrapping, using MEGA 7 software (Version 7.0.26). All Pm TPSs are classified into TPS-d1and TPS-d2 subfamilies based on the taxonomic distribution of TPS families. The accession numbers of all genes are given in Supplementary Table S1. Missing data and positions with gaps were eliminated. (c) Motif analysis of mono- and sesqui-TPSs in P. massoniana.
Forests 15 00244 g002
Figure 3. (a) Predicted secondary structure of Pm TPSs by SOMPA. Numbers represent amino acid residues. Blue indicates residues part of an α-helix, green indicates those part of a β-sheet/turn, red indicates extended strands, and purple indicates residues in random coils; (b) Domain analysis of Pm TPSs using the SMART database, which predicted the presence of two domains that actively participate in the terpene biosynthesis reaction. The numbers on domain boxes represent amino acid residue positions, and the E-value is indicated in parentheses.
Figure 3. (a) Predicted secondary structure of Pm TPSs by SOMPA. Numbers represent amino acid residues. Blue indicates residues part of an α-helix, green indicates those part of a β-sheet/turn, red indicates extended strands, and purple indicates residues in random coils; (b) Domain analysis of Pm TPSs using the SMART database, which predicted the presence of two domains that actively participate in the terpene biosynthesis reaction. The numbers on domain boxes represent amino acid residue positions, and the E-value is indicated in parentheses.
Forests 15 00244 g003
Figure 4. (a) Predicted three-dimensional structure and (b) amino acid sequence alignment of the six terpenoid synthases (TPS) from P. massoniana. The protein resolution method is X-ray. (a) “N” represents the N-terminal domain on the left side of the dash line, while “C” represents the C-terminal domain. The α-helix segments are blue and ordered in white in alphabetical order, β-sheet/turns are in green, and random coil/loops are in white; (b) shows amino acid sequence alignment, α-helical segments are marked and underlined red (* indicate fully aligned amino acids). The homology modeling templates of the six terpene synthases (Pm TPS (+)-α-pinene, Pm TPS (−)-α-pinene, Pm TPS (−)-β-pinene, Pm TPS limonene, Pm TPS selinene, Pm TPS longifolene) are respectively α-bisabolene synthase (PDB no. 3sae.1.A), taxadiene synthase (PDB no. 3p5r.1.A), α-bisabolene synthase, α-bisabolene synthase, taxadiene synthase, and taxadiene synthase. The consistency between the modeling templates and Pm TPSs was greater than 30%.
Figure 4. (a) Predicted three-dimensional structure and (b) amino acid sequence alignment of the six terpenoid synthases (TPS) from P. massoniana. The protein resolution method is X-ray. (a) “N” represents the N-terminal domain on the left side of the dash line, while “C” represents the C-terminal domain. The α-helix segments are blue and ordered in white in alphabetical order, β-sheet/turns are in green, and random coil/loops are in white; (b) shows amino acid sequence alignment, α-helical segments are marked and underlined red (* indicate fully aligned amino acids). The homology modeling templates of the six terpene synthases (Pm TPS (+)-α-pinene, Pm TPS (−)-α-pinene, Pm TPS (−)-β-pinene, Pm TPS limonene, Pm TPS selinene, Pm TPS longifolene) are respectively α-bisabolene synthase (PDB no. 3sae.1.A), taxadiene synthase (PDB no. 3p5r.1.A), α-bisabolene synthase, α-bisabolene synthase, taxadiene synthase, and taxadiene synthase. The consistency between the modeling templates and Pm TPSs was greater than 30%.
Forests 15 00244 g004
Figure 5. Molecular docking of Pm TPSs and substrates. (Left) Contact surface between terpene synthase and substrate molecules. (Right) Specific amino acids forming hydrogen bonds with substrate. Dashed lines indicate possible hydrogen bonds.
Figure 5. Molecular docking of Pm TPSs and substrates. (Left) Contact surface between terpene synthase and substrate molecules. (Right) Specific amino acids forming hydrogen bonds with substrate. Dashed lines indicate possible hydrogen bonds.
Forests 15 00244 g005
Figure 6. Representative GC–FID traces showing products of recombinant Pm mono-TPSs. The labeled numbers represent corresponding products. 1: (+)-α-pinene, 2: (−)-α-pinene, 3: (−)-β-pinene, 4: limonene, 5: myrcene, 6: linalool, IS: internal standard isobutyl-benzene.
Figure 6. Representative GC–FID traces showing products of recombinant Pm mono-TPSs. The labeled numbers represent corresponding products. 1: (+)-α-pinene, 2: (−)-α-pinene, 3: (−)-β-pinene, 4: limonene, 5: myrcene, 6: linalool, IS: internal standard isobutyl-benzene.
Forests 15 00244 g006
Figure 7. P. massoniana resin-based defense induced by M. alternatus feeding/MeJA treatment.
Figure 7. P. massoniana resin-based defense induced by M. alternatus feeding/MeJA treatment.
Forests 15 00244 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wen, Q.; Chen, R.; Xu, T.; Hao, D. Functional Characterization of Terpene Synthases from Masson Pine (Pinus massoniana) under Feeding of Monochamus alternatus Adults. Forests 2024, 15, 244. https://doi.org/10.3390/f15020244

AMA Style

Wen Q, Chen R, Xu T, Hao D. Functional Characterization of Terpene Synthases from Masson Pine (Pinus massoniana) under Feeding of Monochamus alternatus Adults. Forests. 2024; 15(2):244. https://doi.org/10.3390/f15020244

Chicago/Turabian Style

Wen, Quanmin, Ruixu Chen, Tian Xu, and Dejun Hao. 2024. "Functional Characterization of Terpene Synthases from Masson Pine (Pinus massoniana) under Feeding of Monochamus alternatus Adults" Forests 15, no. 2: 244. https://doi.org/10.3390/f15020244

APA Style

Wen, Q., Chen, R., Xu, T., & Hao, D. (2024). Functional Characterization of Terpene Synthases from Masson Pine (Pinus massoniana) under Feeding of Monochamus alternatus Adults. Forests, 15(2), 244. https://doi.org/10.3390/f15020244

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop