Cloning and Functional Analysis of NtHDR in ‘Jinzhanyintai’ of Narcissus tazetta var. chinensis M.Roem

1. Introduction

Floral fragrance consists of a diverse range of volatile organic compounds (VOCs) [1]. In defense against pathogens, parasites, and herbivores, as well as for attracting pollinators, plants release a vast amounts of phytogenic VOCs [2]. Terpenes represent a pivotal class of plant volatiles, and their metabolic pathways have been extensively elucidated. The precursors of terpenoids, including Isopentenyl pyrophosphate (IPP) and dimethylallyl diphosphate (DMAPP), are synthesized through two distinct routes: the cytoplasmic mevalonate pathway and the plastid-localized methylerythritol 4-phosphate pathway (MEP). Subsequently, terpenoids are generated via the catalytic actions of various terpene synthases [3,4,5,6]. Notably, the MEP pathway predominantly governs the biosynthesis of monoterpenes [7], which undergo a seven-step enzymatic reaction, where 1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate reductase (HDR) serves as the final step of the reaction catalase. This last catalytic step precisely regulates the production of IPP and DMAPP mixtures [8,9,10]. Variations in HDR expression have been demonstrated to correlate significantly with alterations in terpene metabolite levels. Studies involving Solanum lycopersicum and Arabidopsis thaliana seeds demonstrated that increased HDR expression led to a corresponding increase in tetra-terpenoid carotenoid content [11]. Furthermore, anatomical investigations of HDR-silenced Nicotiana tabacum leaves have revealed the structural disarray in fenestrated tissue accompanied by a reduction in plastid numbers [12]. In the study of metabolic engineering for the valuable terpenoid tanshinone [13], HDR overexpression was proved to indirectly or directly facilitate the accumulation of these terpenoids. In summary, HDR genes can exert indispensable regulatory roles in plant growth and development, with their overexpression resulting in substantial elevation of indirect precursor levels. This can contribute to the upregulation of biosynthesis pathways for numerous downstream products within the MEP pathway [14].

Chinese Narcissus (Narcissus tazetta var. chinensis M.Roem.), a renowned traditional flower in China, has been cultivated for millennia and is esteemed for its lush green leaves, graceful blossoms, enchanting fragrance, and exceptional ornamental value [15]. Among its significant cultivars, 'Jingzhanyintai' holds a special place. Extensive previous investigations [16,17,18,19,20,21,22] have revealed that the primary components of the floral fragrance volatiles in Chinese Narcissus 'Jingzhanyintai' encompass terpenoids and phenylpropanoids, including monoterpene ocimene, linalool, methyl benzoate, and ethyl benzoate. In recent years, although research on key genes involved in the metabolic pathway of Chinese Narcissus floral fragrance has become an increasing concentration [17,23,24,25,26,27,28], studies that specifically aim to address the HDR remain unreported. To address this knowledge gap, we conducted a comprehensive study where the HDR gene was successfully cloned from the petals of Chinese Narcissus 'Jingzhanyintai' using RT-PCR, leveraging transcriptome data. Subsequent subcellular localization experiments provided insights into its specific function. The expression pattern of the NtHDR was thoroughly analyzed through quantitative real-time PCR analysis. Furthermore, to achieve functional validation, the NtHDR was applied into the generation of Nicotiana benthamiana via transformation, coupled with the employment of the solid-phase microextraction coupled with gas chromatography–mass spectrometry (SPME-GC-MS) technique. This integrated approach established a solid foundation for delving deeper into the floral metabolism and molecular regulatory mechanisms underlying the captivating flower fragrance of Chinese Narcissus.

2. Materials and Methods

2.1. Plant Materials

The plant material utilized in this study comprised Chinese Narcissus cultivar 'Jinzhanyintai' which was purchased from Zhangzhou, Fujian, China. The plants were cultivated within a controlled artificial climate chamber, maintaining a temperature of 20°C for 14 h of light exposure, as well as a temperature of 16°C for 8 h of darkness, alongside a humidity level of 70%. Figure 1 visually presents the varying stages of ‘Jinzhanyintai’ flowering. At different time points, the petals, corona, leaves, bulbs, and roots were promptly subjected to snap-freezing in liquid nitrogen and subsequently stored in a -80°C refrigerator for subsequent analysis.

N. benthamiana was grown in a light incubator at 22°C for 22 h of light and 2 h of dark at 50–70% humidity for transient transformation assays.

2.2. Total RNA Extraction and cDNA Synthesis

The total RNA extraction from 'Jinzhanyintai' petals, corona, bulbs, roots, and leaves was performed through the application of the Quick RNA Isolation Kit (0416-50GX, Huayueyang, Beijing, China), according to the manufacturer’s instructions. The purity and concentration of the obtained RNA samples were subsequently assessed using NanoDrop2000. For cDNA synthesis, a final reaction volume of 20 µL from total RNA (1 μg) was employed, with the utilization of Oligo(dT)₁₅ primers, and AMV Reverse Transcriptase (A3500, Promega, Madison, WI, USA) according to manufacturer’s specifications.

2.3. Cloning and Expression Vector Construction of NtHDR

Using the transcriptome sequencing database as a reference, a homology blast analysis was conducted employing the HDR protein sequence from the model plant Arabidopsis thaliana. Subsequently, specific primers were designed following careful screening of the HDR homologous gene sequence (Table 1). Gene coding sequence (CDS) amplification was conducted with applying PrimeSTAR^® HS (Premix) (Takara, Kusatsu, Japan). The PCR amplification system and protocol were as follows: PrimeSTAR^® HS (Premix, Kusatsu, Japan) 25 µL, 4 µL cDNA template, 2 µL of a 10 μmol·L^-1 mixture of upstream and downstream primers and supplemented with ddH₂O to a total volume of 50 µL. The amplification procedure was as follows: 95°C for 3 min; 95°C for 15 s, 56℃ for 15 s, 72°C for 2 min (35 cycles); 72°C for 5 min; and 4°C +∞. Following PCR, the agarose gel containing the target fragment was excised and purified according to the instructions of the DNA Recovery Kit (DP209, Tiangen, Beijing, China). The purified fragment was then ligated into pCE2 TA/Blunt-Zero Vector (Vazyme, Nanjing, China) and transformed into Escherichia coli DH5α receptor state (CB101, Tiangen, Beijing, China).Positive monoclonal colonies were selected and sent for sequencing to Genewiz. Subsequently, NtHDR was ligated to the pBI121 expression vector through the In-Fusion Cloning, with XmaI and XbaI selected as the enzyme cleavage sites. The enzyme cleavage system and timing were operated according to the instructions provided by BioLab. The homologous amplification primers applied are detailed in Table 2. The PCR amplification system and protocol remained consistent with the aforementioned steps, with the amplification template being the plasmid obtained after correct sequencing in the previous step. The PCR products were recovered by agarose gel purification. The homologous recombination system was prepared following the guidelines outlined in the ClonExpress^® II One Step Cloning Kit (Vazyme, Nanjing, China): 200 ng of the digested pBI121 vector, 60 ng of PCR product, 4 µL of 5×CEII Buffer, 2 µL of ExnaseII, and 20 µL of ddH₂O were mixed and incubated at 37°C for 30 min. The resulting mixture was then transferred into E. coli DH5α receptor state. Positive monoclonal colonies were selected and submitted to Genewiz for sequencing.

2.4. NtHDR Bioinformatics Analysis

The structural domain identification of the proteins was analyzed through adopting the NCBI Conserved Domain Search online tool (https://www.ncbi.nlm.nih.gov/). The amino acid composition, protein molecular weight, theoretical isoelectric point (pI), and stability were predicted through utilizing the online software ProtParam tool (http://web.expasy.org/protparam/). The hydrophobic property and charge distribution of the protein were analyzed via the ProtScale tool (https://web.expasy.org/protscale/). To analyze the transmembrane domain of the NtTPS protein, the TMHMM server 2.0 software (http://www.cbs.dtu.dk/services/TMHMM) was employed. The presence of a signal peptide in NtTPS was predicted with the SignalP4.1 Server (http://www.cbs.dtu.dk/services/SignalP/). The tertiary protein structure predictions were conducted through the application of SWISS-MODEL (https://swissmodel.expasy.org/) platforms, respectively. The sequence similarity of comparison was conducted to analyze the NtHDR sequence results using BLAST from NCBI [32]. The Neighbor-Joining (NJ) method was employed to construct the phylogenetic tree using MEGA 7.0 software with 1000 bootstrap replicates [33,34].

2.5. Fluorescence Quantification of NtHDR in ‘Jinzhanyintai’

Specific amplification primers were designed through applying Primer Premier 6 (Table 1), with primers for reference gene (Actin) selected based on previous studies [35]. Moreover, following the manufacturer’s recommendation, a LightCycler^® 480 fluorescent qPCR instrument and TB Green Premix ExTaqII were utilized for conducting the experiments and the qRT-PCR assays, respectively. Furthermore, each experiment was performed with three technical replicates. The relative expression of the NtHDR gene was calculated using the 2^-△△Ct method[36]. Correlation analysis between the relative expression of the NtHDR gene in various periods of Narcissus flowering and in diverse plant parts was performed with SPSS 25.

2.6. Subcellular Localization Assay of NtHDR Proteins in Tobacco

Subcellular localization studies of NtHDR were conducted through the agrobacterium-mediated transient expression technique. The plant overexpression vector plasmid pBI121-NtHDR-GFP was transferred into GV3101 agrobacterium tumefaciens competent cells (Shanghai Weidi Biotechnology, Shanghai, China) via heat stimulation. The resulting single colonies were cultured for 60 h and subjected to colony PCR for identification. Positive single colonies were selected and transferred to 20 mL of LB liquid medium containing 50 mg/L kanamycin and 25 mg/L rifampicin, followed by being incubated at 28℃ with shaking at 200 rpm until the OD₆₀₀ value reached approximately 0.8 [37]. The bacteria were collected, suspended in a formulation containing 10 mM MES-KOH, 10 mM MgCl₂, and 150 µM acetosyringone, and left for 2–3 h. The suspension was injected into N. benthamiana leaves and incubated under low light conditions for 48 h. Subsequently, chloroplast autofluorescence was observed and imaged with a LSM 880 (ZEISS) laser confocal microscope, with an excitation wavelength of 633 nm and an emission wavelength of 684 nm. Additionally, GFP fluorescence was observed and imaged with an excitation wavelength of 488 nm and an emission wavelength of 533 nm [38,39] (Zhang et al., 2022; Liu, 2017). The plant expression vector plasmid pBI121-GFP was utilized as a positive control and treated as pBI121-NtHDR-GFP.

2.7. Agrobacterium-Mediated Transformation and Identification of N. benthamiana

The plant overexpression vector plasmid pBI121-NtHDR was transferred into GV3101 agrobacterium tumefaciens competent cells through heat stimulation, followed by colony PCR to confirm successful transformation. The subsequent culture procedures were the same as described in Section 4.6. The bacteria were collected, suspended in the MS liquid medium, and allowed to incubate for 2–3 h.

Under aseptic conditions in the ultra-clean table, healthy leaves from well-cultured sterile N. benthamiana seedlings were carefully excised into 0.8 cm × 0.8 cm sized leaves (avoiding the veins as much as possible) using sterilized knives. These leaves were then placed onto a culture medium containing MS+6-BA 2.0 mg/L+NAA 0.2 mg/L and kept in the dark for a period of 2 d. Pre-cultured N. benthamiana leaves were subsequently immersed in a prepared agrobacterium solution for approximately 10 min, during which the solution was continuously agitated to ensure optimal contact with the leaves. After completion, excess bacterial solution was removed using sterile filter paper, and the leaves were rinsed 2–3 times with sterile water. Subsequently, the leaves were placed on sterile filter paper to dry. A new pre-culture medium was prepared by placing three layers of sterile filter paper in a flat position, and the pre-treated leaves were positioned backside up on the filter paper layers. This setup was maintained in the dark for an additional 2 d. Once adventitious buds started to grow, the leaves were transferred to a normal light condition and the culture medium was changed approximately once every 10 d. The adventitious buds presenting excellent growth were excised and transferred to a differentiation culture medium containing MS+NAA 0.2 mg/L + Kan 100 mg/L + Tim 300 mg/L. For rooting, the adventitious buds were cultured in a 1/2MS+Kan 100 mg/L + Tim 300 mg/L medium until they developed 2–3 leaves and reached sufficient height. Once the root system became adequately grown, the cap of the culture bottle was loosened, and sterile water was added for a period of 3–4 d under low light conditions. Prior to transplanting the seedlings into small plastic pots, the medium on the roots was carefully cleaned, and the seedlings were covered with transparent glass or plastic containers for a period of 10 d after transplantation. This step aimed to maintain a high ambient humidity, allowing the seedlings to gradually adapt to the new environment.

Transgenic N. benthamiana DNA was extracted following the protocol outlined in the DNA Extraction Kit (MF735-plus, Juhemei, Beijing, China) to ensure high-quality DNA samples Subsequently, to validate the expression of NtHDR at the protein level in the transgenic plants, petals were carefully selected from plants that had been identified as positive for DNA. Moreover, RNA extraction was performed with the RNA Extraction Kit, followed by reverse transcription to synthesize cDNA. The PCR identification process utilized specific primers, namely 35SF and NtHDR-R (Table 2), to amplify the target sequence.

2.8. Analysis of Volatile Components in Transgenic N. benthamiana Flowers

Entirely open flowers of transgenic and wild-type N. benthamiana were selected with three replicates performed for each assay). The samples were prepared by mixing approximately 0.2 g of randomly selected N. benthamiana flowers in a clean 10 mL round-bottomed headspace vial at room temperature (25±2°C). The vial was sealed using a blue Teflon silicone spacer, leaving about 2/3 of the space at the top. Analysis was performed using a 50/30 μm DVB/CAR/PDMS solid-phase microextraction head and HP-5MS (30 m × 250 μm × 0.25 μm) flexible quartz capillary column with a GC-MS model Agilent 7890A-5975C. The headspace solid-phase microextraction parameters were adjusted based on the work of Sun Miao [40] with appropriate modifications. The aging temperature and time of the extraction head were set at 250°C for 1 h. The extraction process was conducted at 50°C for 30 min, followed by a desorption time of 3 min. The injection port was of the auto injection type. Helium gas was used as the carrier gas, with a column flow rate of 1 mL/min in the pulsed non-split injection mode. The heating temperature of the injection port was maintained at 250°C. For MS analysis, the EI ionization mode was employed, with an electron energy of 70 eV. The ion source and the four-stage rod temperature were set at 230°C and 150°C, respectively. The full scan acquisition mode covered a scan range of 50–600 amu, and the EMV mode was set to the relative value. The equilibration time for the gas phase column chamber was 0.25 min. The ramp-up procedure was as follows: an initial temperature of 55°C held for 3 min, followed by a ramp-up rate of 3°C/min to 148°C for 1 min, then a ramp-up rate of 10°C/min to 210°C for 3 min, and finally a post-run temperature of 55°C for 1 min.

The identification of compounds involved a two-step approach: retention index method and qualitative analysis of volatiles in N. benthamiana flowers with the application of the NIST14 spectral library. In addition, manual spectral analysis was performed, taking into account the mass spectra and retention time indexes reported in relevant published literature. This comprehensive analysis facilitated the identification of specific volatile components.

3. Results

3.1. Characterization of NtHDR

The ORF of NtHDR spanned 1380 bp, encoding a protein consisting of 460 amino acid residues. The NtHDR can be classified under the Isoprenoid_Biosyn_C1 superfamily. The analysis with the ProtParam online tool revealed that the molecular weight of NtHDR was approximately 51.38kDa. The theoretical pI of the protein was calculated to be 5.58. Furthermore, the protein instability coefficient was determined to be 32.35, indicating its relatively unstable nature. The mean hydrophilicity (GRAVY) value of -0.372 suggested that NtHDR possessed a hydrophilic property (Figure 2C). Prediction analyses conducted through the application of Signa lP4.1 Server and TMHMM indicated the absence of signal peptide sites and transmembrane regions in the NtHDR protein, thereby categorizing it as a non-secreted and non-transmembrane protein (Figure 2A,B). To visualize the protein’s 3D structure, SWISS-MODEL was utilized, resulting in a model exhibiting approximately 81.21% similarity to the reference model (Figure 2D).

3.2. Phylogenetic Analysis and Sequence Alignment of NtHDR

The functional domain analysis revealed that NtHDR was attributed to the lytB-ispH superfamily of genes, with the encoded protein exhibiting a highly conserved region. Further analysis of NtHDR sequence homology indicated the highest amino acid sequence similarity with A. officinalis at 88.35%. Comparatively, Z. officinale, Elaeis guineensis, Musa acuminata subsp. malaccensis, and Ananas comosus presented 83.97%, 83.3%, 81.03%, and 80.99% gene similarity, respectively (Figure 3). These findings suggested a significant homology of HDR across various species, with strong conservation observed at the C-terminus of the encoded protein. Phylogenetic tree analysis, performed using MEGA7.0 (Figure 4), demonstrated the relationship between HDR genes of 'Jinzhanyintai' and other species. NtHDR was found to be closely associated with A. officinalis (XP_020241009.1), Z. officinale (XP_042395642.1), Dendrobium officinaleKimura (AGH62555.1), Oncidium flexuosum (ACJ83116.1), E. guineensis (XP_010909277.1), Phoenix dactylifera (XP_008799160.1), M. acuminata subsp. Malaccensis (XP_009412299.1), and A. comosus (XP_020104485.1), as well as other monocotyledons, forming a closely clustered branch. Notably, the closest relatives were observed to be A. officinalis and Z. officinale.

3.3. NtHDR Expression Specificity Analysis

Quantitative real-time fluorescence analysis revealed the expression pattern of NtHDR in various tissues, including petals, corona, leaves, bulbs, and roots, indicating certain variations among them (Figure 5A). Notably, NtHDR exhibited higher expression levels in petals and corona while it presented lower expressions in other tissues, with a ranking order of petals > corona > leaves > roots > bulbs. Moreover, examining the expression level of NtHDR in varying flowering stages revealed interesting trends. In petals, the expression level demonstrated an increasing and then decreasing trend during blooming (Figure 5B), while in the corona, the expression level presented a continuous increasing trend throughout the flowering process (Figure 5C). Overall, the expression pattern of NtHDR demonstrated a regular association with the flowering stage, and the pronounced increase in expression level observed in the corona suggested a potential involvement of NtHDR in the metabolism of floral terpenoids.

3.4. NtHDR Subcellular Localization Analysis

The pBI121-NtHDR-GFP plasmid-containing Agrobacterium GV3101 strain was infiltrated into N. benthamiana leaves, while the pBI121-GFP plasmid-containing Agrobacterium GV3101 strain served as a control. Subsequent observation by utilizing laser confocal microscopy revealed that the fusion protein of the target gene HDR and GFP generated green fluorescence exclusively on chloroplasts (Figure 6). In the control group, a distinct and bright green fluorescence signal of GFP was identified within the nucleus and cytoplasm, which indicated the localization of NtHDR on chloroplasts. This observation aligned with the knowledge that the MEP pathway, a metabolic pathway involved in floral metabolism, occurred within cellular chloroplasts, as predicted by the website and previous knowledge.

3.5. Stable Expression of NtHDR in N. benthamiana

A PCR-based screening method was employed to screen a total of 10 transgenic N. benthamiana plants overexpressing NtHDR, with wild-type plants serving as controls for subsequent investigations (Figure 7). Comprehensive observations, including plant size, leaf morphology, flower color, petal size, and florescence, were conducted on both NtHDR overexpressing N. benthamiana plants and wild-type plants. However, no significant differences were detected between the two groups. The breeding process of transgenic N. benthamiana plants is illustrated in Figure 8 below.

3.6. Determination and Analysis of Floral Volatiles of Transgenic N. benthamiana

Volatile compounds from the flowers of transgenic N. benthamiana plants and wild-type plants were collected using headspace solid-phase microextraction and subsequently analyzed using GC-MS. The analysis revealed discernible differences in the volatile components between the two plant groups(Table 2). Specifically, the transgenic N. benthamiana plants expressing NtHDR exhibited the release of linalool (RT=13.297), benzyl alcohol (RT=10.690), and phenethanol (RT=14.151) as evident in the chromatogram (Figure 9-B). Conversely, these three volatile components were absent in the wild-type plants (Figure 9-A).

Table 2. Main volatile components in fresh flower of different strain N.benthamiana.

Main volatile components	Retention time/min	Relative content/%
		A Wild-type N. benthamiana			B Transgenic N. benthamiana
		CK1	CK2	CK3	H1	H2	H3
Benzyl alcohol	10.690	—	—	—	—	1.288	1.687
Phenylacetaldehyde	10.832	1.682	1.001	—	2.652	4.247	1.105
Linalool	13.297	—	—	—	0.551	1.086	0.777
Phenethyl alcohol	14.151	—	—	—	—	1.812	—
Decamethylcyclopentasiloxane	15.963	10.766	10.488	7.715	16.885	12.974	11.461
2-Isobutyl-3-methoxypyrazine	16.893	0.870	—	—	1.003	1.940	—
Cinnamic aldehyde	20.923	3.203	—	0.233	4.047	—	4.133

Table 2. Main volatile components in fresh flower of different strain N.benthamiana.

Main volatile components	Retention time/min	Relative content/%
		A Wild-type N. benthamiana			B Transgenic N. benthamiana
		CK1	CK2	CK3	H1	H2	H3
Benzyl alcohol	10.690	—	—	—	—	1.288	1.687
Phenylacetaldehyde	10.832	1.682	1.001	—	2.652	4.247	1.105
Linalool	13.297	—	—	—	0.551	1.086	0.777
Phenethyl alcohol	14.151	—	—	—	—	1.812	—
Decamethylcyclopentasiloxane	15.963	10.766	10.488	7.715	16.885	12.974	11.461
2-Isobutyl-3-methoxypyrazine	16.893	0.870	—	—	1.003	1.940	—
Cinnamic aldehyde	20.923	3.203	—	0.233	4.047	—	4.133

4. Discussion

In this study, we cloned an HDR enzyme gene from Chinese Narcissus 'Jinzhanyintai' and successfully introduced it into N. benthamiana plants to achieve stable transformation. The floral fragrance components of the transgenic plants were investigated using SPME-GC-MS. Our analysis revealed the presence of the monoterpene volatile linalool, as well as the phenylpropanoids benzyl alcohol and phenylethanol, in the NtHDR overexpression plants, whereas these components were undetectable in the wild-type plants. In addition, our findings align with previous research conducted by Pang Hongdong [29] who transferred the SAMT enzyme gene serving as a key component in the metabolic pathway of phenylpropanoids/phenylcyclic compounds in Rhamnus phyllanthus into N. tabacum. Pang Hongdong observed higher levels of hexenal, benzyl alcohol, and linalool, while methyl benzoate was not detected. The author hypothesized that the introduction of CpSAMT may have disrupted the synthesis pathway of methyl benzoate in N. benthamiana, thereby leading to more accumulation of benzyl alcohol. Hence, considering the outcomes of our study, we speculated that the production of linalool, benzyl alcohol, and phenylethanol may suggest an unknown association between the phenylpropanoid/phenylcyclic compound metabolic pathway and the terpene metabolic pathway. However, due to the incomplete understanding of the phenylpropanoid/phenylcyclic compound metabolic pathways, further investigations are necessary to elucidate the specific mechanisms underlying their interplay with terpene metabolic pathways. Furthermore, Zhang Tengxun [30] transferred the DXS enzyme gene, an essential component of the MEP pathway in Lilium brownii, into N. tabacum and observed significantly higher linalool content in the transgenic flowers compared to wild-type plants. Similarly, in our experiment, the introduction of the key enzyme gene HDR from the MEP pathway into N. benthamiana resulted in increased linalool content in the transgenic plants’ flowers. Chen Jindou [31] transferred the putative phenylpropane transcription factor CpODORANT into N. tabacum and suggested that the monoterpene ocimene was only present in the control and wild-type N. tabacum flowers, whereas it could not be detected in the flowers of transgenic CpODORANT N. tabacum. Conversely, a small amount of sesquiterpenes geraniene and germacrene D were produced. This observation may help explain the complexity of the floral aroma regulation mechanism in N. tabacum. Based on the aforementioned findings, we inferred that the overexpression of NtHDR in N. benthamiana plants led to increased accumulation of IPP and DMAPP, the floral precursor substances in the MEP pathway. Consequently, this upregulation affected the downstream formation of the monoterpene compound linalool, as well as the phenylpropanoid metabolic pathway in N. benthamiana, resulting in the production of benzyl alcohol and phenylethanol. Overall, our findings provide preliminary evidence supporting the significant regulatory role of the HDR gene in the floral metabolic pathway at the transgenic level.

5. Conclusions

In this study, HDR, a key enzyme gene located at the end of the MEP pathway, was cloned from the petals of Nar. tazetta var．chinensis 'Jinzhanyintai' and named NtHDR. NtHDR played a role in the floral biosynthesis process and was capable of promoting the synthesis of terpene volatiles.