1. Introduction
Terpenoids, a diverse class of compounds produced by plants, encompass important substances such as gibberellins, carotenoids, and chlorophyll. These compounds, also known as isoprene derivatives, are crucial for various aspects of plant growth and development [
1,
2]. The fundamental structural unit of terpenoids is isoprene (C5H8), and they are classified according to the number of these isoprene units they possess. This classification includes semiterpenes (C5), monoterpenes (C10), sesquiterpenes (C15), diterpenes (C20), triterpenes (C30), and so on [
3]. Isoprene in plants is vital for maintaining membrane fluidity, facilitating respiration, and enhancing photosynthesis. In turn, it supports overall plant growth and regulation. These specialized metabolites also play a significant role in mediating the plant’s defense and adaptive responses to environmental challenges. They are integral to allelopathic interactions and plant–pathogen relationships, serving to protect plants from herbivores and disease-causing agents, and the volatiles they produce influence the behavior of insects [
4]. Additionally, terpenoids are known to attract pollinators and seed-dispersing animals, thereby enhancing the survival and propagation of plant species [
5]. The economic significance of isoprenoids extends to various industries, as they serve as crucial components in the production of pharmaceuticals, nutritional supplements, flavorings, pigments, agrochemicals, and disinfectants, and numerous other applications. This underscores the importance of understanding and harnessing the potential of these bioactive compounds for both ecological and commercial purposes [
6].
The biosynthesis of plant terpenoids is a complex process that can be segmented into three principal stages. The initial stage focuses on the synthesis of the fundamental precursors, namely isoprene pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP). The second stage involves the generation of direct precursor compounds, which are then used in the third stage for the synthesis and modification of a diverse array of terpenoid structures [
5]. IPP and DMAPP are indispensable precursors in the manufacture of terpenoids, synthesized primarily through two key metabolic pathways. The first metabolic pathway is the mevalonate (MVA) pathway, which operates in the cytoplasm and commences with acetyl-CoA. The second is the 2-C-methyl-D-erythritol-4-phosphate (MEP) pathway, which takes place in the plastids [
7]. The MVA pathway is prevalent across diverse organisms, including archaea, certain Gram-positive bacteria, yeast, and animals. Conversely, the MEP pathway serves as the primary route for most Gram-negative bacteria, cyanobacteria, and green algae [
8,
9,
10,
11]. Many organisms rely solely on one of these pathways, whereas higher plants and specific algae, such as
Cyanidium caldarium and
Ochromonas danica, utilize both the MVA and MEP pathways [
12,
13,
14]. A key enzyme in the MEP pathway is 4-(cytidine 5′-diphospho) -2-C-methyl-D-erythritol kinase (CMK), which catalyzes the transformation of 4-(cytidine 5′-diphospho)-2-C-methyl-D-erythritol (CDE-MP) to 4-(cytidine 5′-diphospho) -2-C-methyl-D-erythritol 2-phosphate (CDE-ME2P) [
15,
16]. CMK shares significant sequence homology with other kinases, particularly homoserine and mevalvalate kinases. Through studies employing virus-induced gene silencing (VIGS) targeted at the tobacco
CMK gene, its pivotal role in chloroplast development has been uncovered, resulting in observable reductions in both the number and size of chloroplasts, along with decreased levels of metabolites originating from the MEP pathway, including chlorophyll and carotenoids [
17]. In
Ginkgo biloba, two homologous
CMK genes,
GbCMK1 and
GbCMK2, have been successfully cloned and characterized. Their existence suggests a sophisticated regulatory mechanism within the MEP pathway, allowing
G. biloba to efficiently distribute the precursors IPP and DMAPP between its primary and secondary metabolic processes. This highlights the plant’s employment of a multi-faceted regulatory strategy involving CMK isoenzymes to fine-tune its secondary metabolism [
18].
P. massoniana belongs to the subgenus of Pinaceae, a family known as the largest within the gymnosperms, characterized by its possession of divascular bundles. This expansive family comprises approximately 10 genera and around 230 species globally, with China alone contributing to approximately 10 genera and roughly 113 species [
19]. Members of the Pinaceae family are of paramount importance due to their utility as sources of solid wood and high-quality pulp. Additionally, several species within this family are valuable for the extraction and refinement of various chemical raw materials, such as turpentine. Some species are edible or possess medicinal properties, while others are utilized in ornamental gardening [
19,
20].
P. massoniana particularly stands out for its production of rosin, which is abundant in monoterpenes and diterpenes. This resin is harvested not just for extracting turpentine and rosin but also functions as a biofuel, thereby carrying substantial ecological and economic importance [
21]. China accounts for over one-third of the global turpentine production, positioning it as a leading producer. As a result, future research will prioritize understanding the relationship between the components of resin production and pine wood nematodes in the development of high-yielding turpentine varieties [
22]. Methyl-D-erythritol-4- phosphate (MEP) pathway serves as a crucial and indispensable metabolic route for terpenoid biosynthesis in
P. massoniana, with CMK being the sole kinase involved in this pathway [
23]. However, despite its significance, there have been no dedicated studies to date exploring the transcriptional expression and regulatory mechanisms of the
CMK gene in
P. massoniana. In the present study, we elucidate the function of
PmCMK through comprehensive analyses, including tissue-specific expression profiling, heterologous expression, and transient transformation, leveraging the previously identified open reading frame (ORF) and promoter sequences of
PmCMK. These endeavors are not only instrumental in advancing our understanding of the regulatory pathways governing terpenoid biosynthesis but also contribute significantly to the breeding of high-yielding resin-producing strains of
P. massoniana.
2. Materials and Methods
2.1. Plant Material
The experimental materials comprised two-year-old potted seedlings of P. massoniana, originating from the seed orchard located at the Baisha State-owned Forest Farm in Shanghang, Fujian Province, China. These materials were used to clone the ORF and promoter of PmCMK. The samples of P. massoniana with high resin production were sourced from Gaolan Village in Xindong Town, Gaozhou City, Guangdong Province, while the samples with low resin production were collected from Zhong Village in Pintai Town, Yunan County, Guangdong Province. They were used to measure the expression of PmCMK. For the study, tobacco (Nicotiana × sanderae) was cultivated in a controlled laboratory setting, with a photoperiod of 16 h of light alternating with 8 h of darkness, maintaining a constant temperature of 24 °C and a light intensity of 100 μmol m−2 s−1.
2.2. Cloning of the PmCMK and Bioinformatics Analysis of the Protein Sequence
RNA was extracted from two-year-old potted seedlings of
P. massoniana using the FastPure Universal Plant TotalRNA Isolation Kit (Nanjing Vazyme Biotech Co., Ltd., Nanjing, China). First-strand cDNA was synthesized using the One-step gDNA Removal and cDNA Synthesis SuperMix kit (Yeasen Biotech, Shanghai, China). The ORF of the
PmCMK gene (1218 bp) was cloned from the material of
P. massoniana utilizing the sequence from NCBI (GenBank accession number: MW892445.1). Primers CMK-ORF-F and CMK-ORF-R (
Supplementary Table S1) were designed according to the
PmCMK sequence. A 50 μL PCR system was used. The process involved an initial pre-denaturation step at 98 °C for 1 min, followed by 35 cycles of denaturation at 98 °C for 10 s, annealing at 58 °C for 10 sec and extension at 72 °C for 1 min. Finally, a final extension step was performed at 72 °C for 5 min. After cloning, the pEASY-Blunt vector was connected, transferred into
Escherichia coli, and the bacterial solution was sent for sequencing (Tsingke Biotechnology Co., Ltd., Beijing, China).
The sequences of high homology to
PmCMK were obtained by searching the NCBI database with Blastp (
Supplementary Table S2). The analysis and prediction of the PmCMK protein was conducted with Jalview 2.11.3.3 software. The phylogenetic tree was constructed utilizing MEGA 5.1 software.
2.3. Relative Expression Levels of PmCMK in High- and Low-Resin-Producing Materials in P. massoniana
The tender shoots of annual seedlings from high-resin-producing and low-resin-producing
P. massoniana were collected separately. Additionally, seeds from both high- and low-resin-producing
P. massoniana were cultivated under a photoperiod of 16 h of light alternating with 8 h of darkness, maintaining a constant temperature of 24 °C and a light intensity of 100 μmol m
−2 s
−1 for two months, and the entire plants were then harvested. Expression levels of
PmCMK were examined in four types of materials: the tender shoots of annual seedlings from high-resin-producing
P. massoniana, the tender shoots of annual seedlings from low-resin-producing
P. massoniana, the entire two-month-old plants from high-resin-producing
P. massoniana, and the entire two-month-old plants from low-resin-producing
P. massoniana. The four materials are shown in
Supplementary Figure S1. We employed the FastPure Universal Plant TotalRNA Isolation Kit (Nanjing Vazyme Biotech Co., Ltd.) to extract RNA from the samples in high- and low-resin-producing materials.
The differential expression of the PmCMK gene across various tissues was examined through reverse transcription quantitative polymerase chain reaction (RT-qPCR). RT-qPCR was performed using the SYBR Green method. The 10 µL reaction mixture consisted of 1 µL of 20-fold diluted cDNA, 5 µL of SYBR Green Real-time PCR Master Mix, 0.4 µL each of 10 µM primers, and 3.2 µL of ddH2O. The composition of the reaction system and the amplification protocol adhered to the instructions in the Hieff UNICON Universal Blue qPCR SYBR Green Master Mix manual. The relative transcript abundance was determined using the 2−ΔΔCt method. Each qRT-PCR result was conducted with three biological replicates and three technical replicates. Biological replicates were three P. massoniana samples from different single plants with the same growth conditions and the same parts. The three technical replicates were experimental replicates. Data were organized and analyzed using Excel 2021, and graphs were plotted in GraphPad Prism 8.0.2.
2.4. Promoter Cloning and Transient Expression Analysis
DNA was extracted from
P. massoniana using the FastPure Plant DNA Isolation Mini Kit-BOX2 (Nanjing Vazyme Biotech Co., Ltd.). Based on the sequence of the
PmCMK promoter (1654 bp), a set of primers was designed for amplification, utilizing the genomic DNA (gDNA) of
P. massoniana as the template. To predict the cis-acting elements within the cloned promoter sequences, the PlantCARE database was utilized. Primers, incorporating restriction sites for HindIII and BamHI, were designed based on the pBI121 vector, which contains GUS tags and the promoter sequence. Transient expression vectors, pBI121-PmCMK pro::GUS, were constructed and subsequently transformed into
Agrobacterium GV3101 cells, along with the empty plasmid vector pBI121-CaMV::GUS as a control. After that, the leaves of a 20-day-old tissue culture of tobacco seedlings under aseptic conditions were transformed using the
Agrobacterium-mediated method [
5]. After the transformation, the plants underwent a 2-day dark treatment at 28 °C. Subsequently, the leaves of different groups were sprayed with 100 μM indole-3-acetic acid (IAA), 100 μM salicylic acid (SA), 100 μM methyl jasmonate (MeJA), and 20% polyethylene glycol 6000 (PEG6000). Each experiment involved three biological replicates and three technical replicates. The leaves were then stained with GUS dye, incubated overnight at 37 °C, and decolorized using 75% ethanol to remove chlorophyll. Finally, the GUS staining was visualized under a stereo-fluorescence microscope (M205FA, Leica, Wetzlar, Germany) for photographic documentation. This process allowed for a qualitative assessment of GUS expression in response to various treatments, providing valuable insights into the regulatory mechanisms of the
PmCMK promoter.
2.5. Genetically Modified Tobacco
In this experiment, an
Agrobacterium-mediated transformation method was employed for the infection of wild safflower tobacco. Initially,
Agrobacterium was cultured in 1 mL of Luria–Bertani (LB) medium and incubated with shaking at 28 °C for one day. The bacterial culture was then transferred to 5 mL of fresh LB medium supplemented with 50 mg/mL kanamycin (Kan) and incubated overnight with shaking at 28 °C. The following day, 1 mL of the overnight bacterial culture was aseptically transferred to 50 mL of LB medium containing 50 mg/mL Kan and 25 mg/mL rifampicin (Rif), and the culture was shaken at 28 °C until the optical density at 600 nm (OD
600) reached a range of 0.4 to 0.6. The infiltration solution was prepared by ddH
2O containing 10 mM MgCl
2, 10 mM MES, and 200 μM AS, and the ratio of the infiltration solution used to resuspend the bacterial suspension to the original volume of the bacterial suspension was 1:1. Leaves from two-month-old tobacco plants were trimmed into 1 × 1 cm pieces and immersed in the infection solution, shaking gently for approximately 10 min. The inoculation method referred to the method used by Heidari [
24]. Afterward, the infection solution was decanted, and the leaves were blotted dry before being placed on Murashige and Skoog (MS) medium containing 2.0 mg/L 6-benzylaminopurine (6-BA) and 0.5 mg/L IAA to facilitate the infection process. The plants were then subjected to two days of dark incubation followed by a photoperiod of 16 h of light and 8 h of darkness with a constant temperature of 24 °C and a light intensity of 125 μmol m
−2 s
−1. Five days post-infection, the tobacco leaves were transferred to an MS medium containing 2.0 mg/L 6-BA, 0.5 mg/L IAA, 500 mg/L cefotaxime (Cef), and 50 mg/L Kan to select for transformed cells. Under the same conditions as the previous stage, approximately 30 days later, small green regenerated shoots began to emerge around the inoculated leaves. These shoots were excised and cultured in a root induction medium (1/2 MS supplemented with 25 mg/L Kan and 100 mg/L Cef) to encourage root development. Once the tissue culture seedlings had developed a robust root system, they were transplanted into soil to continue their growth and development. The tobacco seedlings were genetically analyzed at the DNA level to confirm the presence of the transgene and to quantify the number of successfully transformed seedlings. This process ensured the successful generation of transgenic tobacco plants for further study and application.
2.6. The Contents of Chlorophyll a, Chlorophyll b, Carotenoid, and CMK
Fresh transgenic tobacco leaves were selected as the sample. The contents of chlorophyll a, chlorophyll b, and carotenoid were determined by the reagents in the plant chloroplast pigments’ (chlorophyll a, b, and carotenoids) detection kit (Shanghai Hengyuan Biotechnology Co., Ltd., Shanghai, China). CMK activity in transgenic tobacco was determined by the plant CMK ELISA kit (Shanghai Hengyuan Biotechnology Co., Ltd.).
2.7. Relative Expression Levels of Genes Related to Terpenoid Synthesis Pathways in Transgenic Tobacco
The relative expression levels of genes associated with the terpenoid synthesis pathway were determined using the RT-qPCR method. The selected tobacco genes were as follows:
DXS and
DXR from the MEP pathway,
HMGR and
MK from the MVA pathway,
GGPPS from the GGPP pathway,
PSY and
ZDS from the lycopene pathway,
β-LCY prior to the downstream β-carotene generation, and
BCH2 from the zeaxanthin pathway. The primer sequences are provided in
Supplementary Table S1. The gene name and GenBank ID of genes related to the terpene synthesis pathway are provided in
Supplementary Table S3. The methodology and reaction system employed for RT-qPCR were identical to those described in
Section 2.3.
2.8. Statistical Analysis
Each experimental result was derived from three biological replicates and three technical replicates. Data were organized and calculated using Excel, and GraphPad Prism 8.0.2 was utilized for plotting. * represents the p-value of the statistical test, where * indicates p < 0.05, ** indicates p < 0.01, *** indicates p < 0.001, and **** indicates p < 0.0001.
4. Discussion
Terpenoids constitute a class of natural compounds that are abundantly found in higher plants, primarily derived from isoprene (C5H8). Terpenoids generated by plants, such as gibberellins, carotenoids, and chlorophyll, play important roles in plant growth and development [
2]. The metabolites derived from terpenoid compounds can specifically function in mediating interactions between tissues or defense and adaptive responses pertaining to the environment [
25]. The MEP pathway, as a crucial metabolic route for the synthesis of terpenoids in plants, represents a significant component in terpenoid research. This article focused on the
PmCMK gene, which has yet to be investigated in
P. massoniana, aiming to clone and elucidate its function, thereby contributing novel insights to the study of the MEP pathway in
P. massoniana. In this paper,
PmCMK gene was cloned, and a bioinformatics analysis was performed. The results of the conserved domain analysis showed that
PmCMK had the GHMP kinase domain of the GHMP family. This family comprises unique ATP-dependent enzymes, wherein GHMP kinases feature a conserved Gly/Ser-rich region that plays a crucial role in ATP binding, isoprenoid and amino acid biosynthesis, as well as carbohydrate metabolism [
26,
27,
28]. In this paper, materials with a high and low resin production were collected, and the relative expression levels were measured. The results indicated that in both the shoots and the entire plants, the expression of
PmCMK was notably more pronounced in high-resin seedlings as compared to their low-resin counterparts (
Figure 3), which also verified the important role of
PmCMK in the resin production pathway of
P. massoniana. Previous identification results have established a positive correlation between the expression level of the
PmGGPPS gene and resin production, with the expression of
PmGGPPS largely aligning with the varying trends in resin production, both high and low [
29]. This provides a certain reference for further research on the expression of terpene-related genes in high- and low-resin-producing
P. massoniana seedlings in the future.
A correlation between MeJA and terpene metabolism has been established in numerous studies [
30,
31]. MeJA treatment increases the production of monoterpenes and diterpenes of wood in
Picea abies, along with a corresponding increase in the activities of their respective terpenoid synthases, indicating that MeJA can enhance the transcription of genes responsible for resin acid biosynthesis [
32,
33]. Previous studies have demonstrated hormonal effects on terpene production and
TPS gene expression, with MeJA serving as a crucial component of terpenes that play a vital role in defense mechanisms in kiwifruit (
Actinidia chinensis Planch.). An exogenous application of SA to kiwifruit produced significant damage on the surface of kiwifruit berries, whereas no notable phenotypic changes were observed after MeJA treatment. Interestingly, when kiwifruit berries were treated with a combination of MeJA and SA, the SA-induced damage on the fruit surface was significantly reduced, indicating that MeJA is effective in mitigating SA-induced damage [
34]. Moreover, previous research has revealed disparities in the expression of
PmCMK in
P. massoniana under SA and MeJA treatments, which is consistent with the identification of cis-acting elements within the
PmCMK promoter presented in this paper [
28].
In this study, the content of CMK, chlorophyll, and carotenoid in the overexpression lines was observed to be higher than in the wild type, suggesting that the overexpression of
PmCMK affects the synthesis of these terpenoids. Furthermore, we discovered an upsurge in the expression levels of
DXS and
DXR in the MEP pathway of the transgenic tobacco. Based on previous research on
DXS and
DXR in the terpenoid synthesis pathway, DXS stands out as a highly regulated enzyme, acting as a crucial rate-limiting step in the MEP pathway across various plant species, and is an effective target for regulating the synthesis of terpenoids in plants [
35,
36].
PmCMK increased the expression levels of
DXS and
DXR, which greatly affected the rate of the MEP pathway, and thus affected the synthesis and regulation of terpenoids. In transgenic plants, the expression of
HMGR in MVA pathway also increases, while the expression of
MK decreases, indicating that genes in the MEP pathway have a certain impact on genes in the MVA pathway. Furthermore, the expression levels of
PSY,
ZDS,
β-LCY of the β-carotene pathway and
BCH2 of the zeaxanthin pathway in transgenic lines were higher than those of the wild type. This suggests that
PmCMK genes have a greater influence on the downstream terpenoid processing pathways, and that the increase in the expression of these genes also corresponds to the increase in the expression of chlorophyll and carotene in the previous transgenic tobacco. This indicates that
PmCMK affects the production of terpenoids to a certain degree, and this effect is complex, as various genes exert different effects during distinct stages of the terpenoid synthesis pathway. This is consistent with the results shown in the study of the tomato fruit’s MEP pathway, where
SlWRKY35 activation positively regulated carotenoid biosynthesis. In that study, they found that
SlWRKY35 overexpression could increase the expression of most genes in the carotenoid synthesis pathway. However, compared to the downstream pathway, the expression of genes in the upstream carotenoid synthesis pathway was notably higher. The overexpression of
SlWRKY35 had a more profound effect on the expression of MEP pathway genes, particularly
DXS1 in tomatoes. This underscores the complexity of the terpenoid synthesis pathway [
37].