Abstract
Weak flavor in flue-cured tobacco compromises quality, making aroma enhancement crucial. Sesame cake fertilizer (SF) has the potential for improving tobacco aroma, but its effects on aroma components and mechanism remain unclear. Here, a four-year field experiment was conducted in Southern tobacco region of China to compare SF with conventional fertilization (CK). We investigated how SF influences soil quality, root and leaf development, aroma precursor accumulation, and volatile aroma composition to enhance tobacco aroma quality. Our results indicated that SF improved soil structure, pH, organic matter, and cation exchange capacity (CEC). It promoted root growth, dry matter accumulation, and root activity and CEC. Leaf photosynthesis and plastid pigment content increased due to enhanced chloroplast ultrastructure. SF also boosted glandular trichome density and secretion, leading to higher aroma precursor accumulation, particularly cembratriene-diol. After curing, 17 of 18 differential volatile aroma substances were upregulated, including carotenoids, cembratriendid alkyl degradation products, esters, terpenes, and heterocyclic compounds. Further, SF significantly increased sensory quality of flue-cured tobacco by promoting aroma quality and volume, and electronic nose analysis also confirmed this. Therefore, SF improves tobacco aroma by enhancing soil health, root growth, and leaf precursor formation. The “soil-root-leaf-differential aroma substances” framework highlights its role in increasing carotenoid and cembratriene-diol content, contributing to higher volatile aroma concentrations. This study highlights the potential of SF as a sustainable agricultural product for improving soil health and tobacco quality.
Highlights
Sesame cake fertilizer (SF) enhances soil quality and promote tobacco root growth.
SF increased tobacco plastids by improving chloroplast ultrastructure in leaves.
SF boosted aroma precursors by enhancing glandular trichome density and secretion.
SF enriched tobacco aroma components and enhanced sensory quality.
1 Introduction
Tobacco (Nicotiana tabacum L.) is an important industrial crop worldwide, and flue-cured tobacco has the largest planting area (World Health Organization, 2017; ). Based on ecological, sensory, chemical, and metabolic characteristics, flue-cured tobacco is classified into eight types in China (). Fujian Province, the main flue-cured tobacco producing area in southern tobacco region of China, yields high-quality fresh-sweet and honey-sweet tobacco leaves (). However, over the past few decades, excessive chemical fertilizer use to boost yield has reduced soil quality and negatively affected crop yield and quality, including tobacco (; Zhang et al., 2018a). Soils in the Fujian tobacco-growing area are mainly acidic red and yellow soils, and prolonged excessive chemical fertilizer use has reduced organic matter, increased acidification and compaction, and imbalanced soil nutrients, resulting in uncoordinated chemical composition and weakened tobacco leaf aroma (; Zhang et al., 2018a; ). Therefore, exploring effective fertilization application strategies to improve soil and tobacco quality is imperative.
Organic fertilizers have proven effective as soil amendments and have been widely studied in crops (Wei et al., 2016; Yang et al., 2025) such as wheat (), maize (), rice (Wang et al., 2024a), and tobacco (Yan and Liu, 2020). They are typically derived from animal manure, plant residues, and microbial sources rich in organic matter and beneficial microorganisms (Thuriès et al., 2002; ). Organic fertilizers improve soil structure and quality, regulate soil microbial communities, and enhance crop yield and quality (; Yang et al., 2025). Animal manure-based organic fertilizers can improve soil structure and fertility and effectively control pests and diseases (Rowen et al., 2019). Replacing 15% of chemical nitrogen fertilizer with bioorganic fertilizer considerably improved tobacco yield and quality (Wang et al., 2024b). Combining biochar and high-carbon fertilizer with chemical fertilizers improves soil carbon and nitrogen pools, increases plant growth-promoting bacteria, enhances the aroma of tobacco leaves, and improves tobacco quality (Yan et al., 2020). Cow manure combined with various botanical oil meal organic fertilizers, including soybean, rapeseed, peanut bran, and sesame cake, improved tobacco yield and quality; however, only the combination with sesame cake substantially increased the pH of neutral tobacco soil (). Collectively, the effects of different organic fertilizers vary in terms of soil properties and flue-cured tobacco styles.
Sesame cake fertilizer is a high-quality organic amendment that enhances tobacco quality, particularly leaf aroma (Wu et al., 2006). This effect is attributed to its rich organic matter, essential nutrients, and unique volatile components (e.g., sesamol and sesamin) that release aromatic substances during fermentation (; ). Aroma is central to tobacco style characteristics and quality, driven by specific aroma components (). Tobacco aroma formation is influenced by the growth environment, including soil type, fertility, and pH (Yan and Liu, 2020; ; Yan et al., 2025). Soil physicochemical conditions directly affect tobacco root growth and nutrient absorption, thereby influencing leaf development and aroma substance synthesis (Yan et al., 2020; ). Glandular trichome secretions on tobacco leaves are a primary source of aromatic substances that determine tobacco aroma after curing (). In our previous research, we found that applying 600 kg ha-1 of sesame cake fertilizer effectively increased flue-cured tobacco sensory quality by improving aroma in Fujian (). However, the mechanism underlying the increase in tobacco aroma and sensory quality after sesame cake fertilizer application remains unclear.
Therefore, a four-year fertilization trial in tobacco was conducted to (a) elucidate the effects of sesame cake fertilizer on soil physicochemical properties, tobacco growth, and flue-cured quality; and (b) systematically explain how it improves tobacco leaf quality using a “soil-root-leaf-differential aroma substances” framework.
2 Materials and methods
2.1 Site description and experimental design
The field experiment began in 2021 at Dongliu (25°20’N, 115°90’E), Fujian Province, China, with sampling primarily conducted in 2023 and 2024. This area belongs to the subtropical monsoon temperate climatic zone, with average precipitation from 151 to 236 mm and average temperature ranging from 18.3 to 20.6 °C during the tobacco growing season (Supplementary Figure S1). The study site soil is classified as red loam. The initial chemical characteristics of the soil were as follows: soil total nitrogen (TN) 1.13 g kg−1, available phosphorus (AP) 207.85 mg kg−1, available potassium (AK) 202.19 mg kg−1, organic matter (OM) 18.36 g kg−1, and pH 4.52. The tobacco variety used in this study was Yunyan 87, with plants transplanted in February and harvested in July.
The field experiment was conducted during the tobacco growing seasons from 2021 to 2024 using a randomized block design with two treatments: (1) control (CK), conventional fertilizer; (2) sesame cake fertilizer (SF), applied 600 kg ha-1 sesame cake fertilizer on top of the conventional fertilizer (Supplementary Table S1). The sesame cake fertilizer supplied by Zhumadian Kangbo Huixin Agricultural Science and Technology Co. China. Both CK and SF treatments received equivalent total inputs of nitrogen (150 kg N ha-1), phosphorus (120 kg P ha-1), and potassium (510 kg K ha-1). Organic fertilizer was applied once as the base fertilizer. Inorganic fertilizers were applied partly as a base dose (60%) and partly as topdressing (40%) across three times. Each treatment had three replicates for a total of six plots, each plot was 12 m long and 10 m wide. Ridge cultivation was conducted with a row spacing of 110 cm and a hole spacing of 45 cm. All plots received the identical irrigation, herbicide, and pesticide in accordance with local technical specifications for high-quality tobacco production.
To accurately observe changes in glandular trichomes, we conducted a pot experiment at the Qilin North Farm of South China Agricultural University, Guangzhou Province, China (23°09’N, 113°15’E) from March to June 2024. Two treatments (CK and SF) were set up with a total of 32 pots (16 pots per treatment). The tobacco variety Yunyan 87 was grown in plastic pots (39 cm diameter, 34 cm height) filled with 30 kg of clay soil from the top layer (0–40 cm) of the field, with one plant per pot. Before transplanting, the fertilizer was thoroughly mixed with the soil in each pot and incubated outdoors for seven days. The fertilizer application rate was calculated as the ratio of 1100 plants ha-1 in the field (Supplementary Table S1).
2.2 Sampling and measurements
2.2.1 Soil properties
Soil samples were collected from a 0–40 cm depth in each plot using the five-point sampling method at the tobacco maturity stage in 2023 and 2024. Each sample was dried naturally and sieved through a 50-mesh sieve to measure the relevant indexes. Total carbon (TC) and TN were measured using a C/N elemental analyzer (Vario MAX, Elementar, Germany; Yan et al., 2020). Total phosphorus (TP) was quantified using an ascorbic acid colorimetric assay after digestion with H2SO4-HClO4. The total potassium (TK) was measured using a flame photometer. Available nitrogen (AN), AP, and AK were measured using alkali diffusion, molybdate-ascorbic acid, and flame atomic absorption spectroscopy, respectively. Soil pH was determined using the glass electrode method with water extraction (water-to-soil mass ratio of 2.5:1). Soil electrical conductivity (EC) was determined using a glass electrode (DDS-11A, Shanghai, China) in a 1:5 soil: water suspension (Zhang et al., 2024a). The soil organic matter (OM) was calculated using the potassium dichromate volumetric method with external heating. The soil cation exchange capacity (CEC) was determined using cobalt hexachloride hexaammonium cobalt leaching spectrophotometry (HJ 889-2017). Soil bulk density, porosity, and water content were quantified using the ring-knife, soil aggregate structure, and oven drying methods, respectively (Zou et al., 2018; Wang et al., 2024a).
2.2.2 Tobacco root traits and dry matter accumulation
Intact tobacco roots were collected from three randomly sampled plants per plot at the vigorous growth (55 d after transplanting) and maturity (110 d after transplanting) stages. Each root sample was gently cleaned, and the lateral roots were separated from the main root. Roots were separately arranged and scanned with a flatbed scanner (Microtek ScanMaker i850), and root images were analyzed using WinRHIZO software to quantify the total root length (RL), surface area (RSA), and mean diameter (RMD; ). Root activity (RA) was determined using the triamcinolone tetrazolium chloride method (Zhao et al., 2023). Root EC and CEC were measured using a conductivity meter and the drenching method, respectively (Zhang et al., 2024a). Three plants from each plot were sampled at the same stage and dissected into roots, stems, and leaves. All separated organs were cleaned and oven-dried at 80 °C to a constant weight, and their dry weights were recorded.
2.2.3 Tobacco leaf areas, function, and ultrastructure
Six plants per plot were randomly selected to measure plant height, leaf length, and leaf width at both the vigorous growth and maturity stages. Leaf area was calculated as average leaf length × average leaf width × 0.6345 (). The net photosynthetic rate (Pn), transpiration rate (Tr), stomatal conductance (Gs), and intercellular CO2 molar fraction (Ci) of the middle leaves (positions 9–12 from the top) of three randomly selected plants per plot were measured between 09:00 and 11:00 am on sunny days using a Targas-1 portable photosynthesis tester (PP Systems, USA; Tor-Ngern et al., 2021).
For pigment analysis and chloroplast ultrastructure, 0.2 g of fresh leaf tissue was collected from the same vein as the central leaf, extracted with 20 mL of 95% ethanol, and incubated in the dark for 24 h. Chlorophyll a (Chla), chlorophyll b (Chlb), and carotenoid (Car) contents were determined by measuring absorbance at 665, 649, and 470 nm, respectively (Wu et al., 2014). The same position of the central leaf blade, between the 6th and 8th branch veins, was selected; tissues were cut into 1 mm × 1 mm pieces and immediately fixed in 5% glutaraldehyde solution. To observe the chloroplast ultrastructure, the fixed tissues were further treated with 1% osmium tetroxide, dehydrated in graded acetone (30–100%), and embedded in epoxy resin. Ultrathin sections were then prepared using an ultrathin microtome (Leica EM UC6, Wetzlar, Germany), double-stained with uranyl acetate and lead citrate, and observed and photographed using a transmission electron microscope (JEM-1400 Plus, Japan), with ten fields of view for each treatment (Xin et al., 2024).
2.2.4 Tobacco leaf glandular trichomes and their secretion content
At the vigorous growth and maturity stages, fresh middle leaves (9th position from the top) from three randomly selected plants per treatment, avoiding the main veins, were cut and divided into two samples. One sample was immersed in a 0.2% aqueous rhodamine B solution for 20 min, rinsed three times with distilled water, and dried naturally under ventilated conditions. A 5 mm × 6 mm block of tissue was excised and observed with a super-depth-of-field microscope (VHX-500, Japan) to assess glandular trichome morphology from ten fields of view. The other sample was immediately stained with Sudan III solution, heated, dripped with dilute glycerol, and placed under a stereomicroscope (×30) to record glandular trichome density (number per millimeter) based on the average of 10 different fields of view (1.064 mm2 individual field area; Supplementary Figure S2, Zhang et al., 2018b).
Twenty leaf discs (11 cm diameter) per treatment were selected from the middle leaves to determine the leaf glandular trichome secretion content. The leaf discs were immersed in 1000 mL dichloromethane for 2 s, repeated eight times, and then 1 mL of internal standard (2.020 mg mL-1 sucrose octaacetate and 2.542 mg mL-1 heptadecan-1-ol alcohol mixture) was added. The extract was concentrated using a rotary evaporator and dried with a nitrogen blower. Gas chromatography–mass spectrometry (GC–MS) was conducted as previously described (Wang et al., 2021).
2.2.5 Volatile aroma substances content of flue-cured tobacco
GC–MS analysis of volatile organic compounds (VOCs) was performed using an Agilent 7697A-8890-5977B system (Agilent, Palo Alto, CA, USA). A VF-WAXms capillary column (25 m × 0.25 mm × 0.2 µm, Agilent CP9204) was used to separate the VOCs. VOCs were extracted through headspace solid-phase microextraction (HS–SPME) as follows: Six replicates per treatment of de-stemmed and crushed flue-cured middle tobacco leaves (C3F) were collected in 20 mL headspace sample bottles. Subsequently, 2.0 µL of internal standard (50 µg/mL n-Pentadecane-d32) was added to the headspace vials. The mixture was sealed and preheated in a water bath at 80 °C for 20 min to achieve headspace equilibrium. The solid-phase microextraction (SPME) fibers were then inserted into the headspace for 20 min to adsorb VOCs. A 1 µL aliquot was injected into the GC–MS system in split mode (10:1) for analysis. The samples were separated with a VF-WAXms (25 m × 0.25 mm × 0.2 µm) capillary column with the following GC column temperature program: held at 40 °C for 2 min, increased to 100 °C at 5 °C/min, then to 230 °C at 15 °C/min and held for 5 min, followed by an additional 2 min at 230 °C. Helium (99.999%) was used as the carrier gas at a flow rate of 2 mL min-1. Mass spectrometry was conducted in electron impact ionization mode at 70 eV, scanning from 50 to 500 m/z at 3.2 scan/s, with the ion source maintained at 280 °C. VOCs were identified by comparing the mass spectra of all detected metabolites with those in the NIST (version 2017), MS-DIAL (version 2021), and other public databases. Data were analyzed using the Majorbio cloud platform (https://cloud.majorbio.com; Zhang et al., 2024b).
2.2.6 Sensory quality and electronic nose of flue-cured tobacco
De-stemmed and shredded flue-cured middle tobacco leaves (C3F), with three replicates per treatment, were selected to evaluate the sensory quality and analyzed using an electronic nose. Sensory evaluation was performed by five experts from China Tobacco Guangdong Industrial Co., Ltd. Following the Guangdong tobacco sensory standard, evaluation indices included aroma quality (20%), aroma volume (35%), offensive taste (20%), irritancy (10%), and taste (15%), with a maximum total score of nine points. The overall odor profile was obtained using an electronic nose (PEN3, Airsense GmbH, Schwerin, Germany). Approximately 2.5 g of flue-cured tobacco sample was added to 20 mL of brown solution and equilibrated at room temperature for 1 h. Each sample was analyzed for 90 s at an injection rate of 300 mL min-1, and the average sensor response between 65 and 67 s was used for further analysis. The electronic nose sensor array comprised 10 metal-oxide semiconductor chemosensors to detect different odor characteristics (Supplementary Table S2; Zhang et al., 2024b).
2.3 Statistical analysis
Statistical analyses were performed using SPSS (version 26.0; SPSS Inc., Chicago, IL, USA). Significant differences were determined using student’s t-test at a 0.05 probability level. Orthogonal partial least squares discriminant analysis was performed using the SIMCA 14.1 (Umetrics, Ummea, Sweden). Different aromatic substances were screened based on variable importance in projection (VIP) > 1 and P < 0.05. The aroma compound descriptions were obtained from the flavor.net websites at http://www.flavornet.org and http://www.thegoodsc (Zhang et al., 2021; Wen et al., 2023). The Mantel test was performed using a Chiplot (https://www.chiplot.online), and the structural equation modeling (SEM) was constructed using the R studio “lavaan” package.
3 Results
3.1 Soil physical and chemical properties
The soil properties were considerably different after the application of SF (Table 1; Supplementary Table S3). For soil physical properties, bulk density and fractal dimension (D) under SF treatment significantly decreased by 8.20% and 6.02% on average across the year compared with CK. Water content, total porosity, proportion of large aggregates (> 0.25 mm) (R > 0.25), geometric mean diameter (GMD), and mean weight diameter (MWD) of the SF treatment on average across the year significantly increased by 22.39%, 4.63%, 8.14%, 15.17%, and 20.57%, respectively, compared with CK. For soil chemical properties, all chemical indices improved after the application of sesame cake fertilizer. Compared with CK-treated soil, the TC, TN, TK, AP, AK, OM, pH, EC, and CEC in SF-treated soil on average across the year were significantly increased by 17.23%, 10.58%, 12.20%, 15.86%, 5.43%, 19.05%, 6.39%, 22.28%, and 30.88%, respectively.
Table 1
| Properties | 2023 | 2024 | ||
|---|---|---|---|---|
| CK | SF | CK | SF | |
| TC (g·kg-1) | 10.47 ± 0.07b | 12.39 ± 0.14a | 10.62 ± 0.18b | 12.33 ± 0.05a |
| TN (g·kg-1) | 1.68 ± 0.01b | 1.81 ± 0.01a | 1.65 ± 0.01b | 1.87 ± 0.01a |
| TP (g·kg-1) | 0.93 ± 0.00a | 0.94 ± 0.00a | 0.94 ± 0.02a | 1.01 ± 0.02a |
| TK (g·kg-1) | 46.54 ± 0.13b | 52.15 ± 0.24a | 46.84 ± 0.42b | 52.61 ± 0.25a |
| AN (mg·kg-1) | 139.21 ± 0.89a | 143.78 ± 1.09a | 141.88 ± 1.64a | 146.90 ± 1.07a |
| AP (mg·kg-1) | 287.47 ± 1.26b | 319.79 ± 1.53a | 292.65 ± 1.15b | 352.34 ± 0.71a |
| AK (mg·kg-1) | 221.34 ± 1.06b | 232.65 ± 1.19a | 230.94 ± 1.24b | 244.18 ± 1.19a |
| рH | 4.60 ± 0.01b | 4.88 ± 0.01a | 4.59 ± 0.09b | 4.89 ± 0.01a |
| EC (ms·cm-1) | 0.17 ± 0.01b | 0.21 ± 0.01a | 0.16 ± 0.00b | 0.19 ± 0.01a |
| OM (g·kg-1) | 23.68 ± 0.12b | 27.55 ± 0.24a | 25.12 ± 1.04b | 30.55 ± 0.67a |
| CEC (cmol(+)·kg-1) | 12.59 ± 0.25b | 15.54 ± 0.16a | 7.34 ± 0.81b | 10.54 ± 0.36a |
| Bulk density (g·cm-3) | 1.01 ± 0.00a | 0.94 ± 0.00b | 1.03 ± 0.03a | 0.93 ± 0.02b |
| Water content (%) | 0.23 ± 0.01b | 0.28 ± 0.01a | 0.25 ± 0.01b | 0.30 ± 0.01a |
| Total porosity (%) | 0.54 ± 0.00b | 0.57 ± 0.00a | 0.56 ± 0.01b | 0.59 ± 0.01a |
| Capillary porosity (%) | 0.50 ± 0.00a | 0.51 ± 0.00a | 0.52 ± 0.01a | 0.53 ± 0.01a |
| Non-Capillary porosity (%) | 0.05 ± 0.00a | 0.06 ± 0.00a | 0.04 ± 0.01a | 0.06 ± 0.01a |
| R>0.25 mm (%) | 64.13 ± 0.14b | 69.63 ± 0.24a | 66.14 ± 1.42b | 71.26 ± 0.99a |
| GMD (mm) | 1.16 ± 0.01b | 1.26 ± 0.09a | 1.21 ± 0.04b | 1.43 ± 0.01a |
| MWD (mm) | 0.53 ± 0.01b | 0.67 ± 0.01a | 0.58 ± 0.03b | 0.71 ± 0.01a |
| D | 2.46 ± 0.01a | 2.24 ± 0.01b | 2.54 ± 0.01a | 2.46 ± 0.01b |
Soil physical and chemical properties at tobacco maturity stage in 2023 and 2024 after applying sesame cake fertilizer.
CK, conventional fertilization treatment; SF, sesame cake fertilizer treatment; TC, total carbon content; TN, total nitrogen content; TP, total phosphorus content; TK, total potassium content; AN, available nitrogen content; AP, available phosphorus content; AK, available potassium content; EC, electrical conductivity; OM, organic matter; CEC, cation exchange capacity; R > 0.25mm (%), the proportion of large aggregates > 0.25mm; GMD, geometric mean diameter; MWD, mean weight diameter; D, fractal dimension. Values are presented as mean ± standard error (n=3). Different letters denote significant differences according to Student’s t-test at P < 0.05.
3.2 The root system traits and dry matter accumulation
At the vigorous growth stage, the year significantly affected RSA, RMD, and EC; fertilizer application significantly affected RMD and RA; while the year and fertilizer interaction only significantly affected RSA (Table 2). Compared with CK, the RMD significantly increased by 13.92% in 2024 after applying SF. On average across the year, RA was significantly increased by 54.32%. At the maturity stage, the year did not significantly affect any root traits; fertilizer application significantly affected RL, RMD, EC, and CEC, and the year and fertilizer interaction only significantly affected RAS. On average, across the years, RL, RSA, RMD, EC, and CEC significantly increased by 42.34%, 43.42%, 12.41%, 25.99%, and 28.01%, respectively, in the SF-treated soil compared to CK-treated soil. RA did not significantly differ between the CK and SF treatments.
Table 2
| Year | Growth stage | Treatment | Root physical traits | Root chemical traits | ||||
|---|---|---|---|---|---|---|---|---|
| RL (cm) | RSA (cm²) | RMD (mm) | RA (ug.g.h) | EC (ms cm-1) | CEC (cmol (+) kg-1) | |||
| 2023 | Vigorous growth | CK | 5478.19 ± 421.60a | 1389.89 ± 126.29a | 0.76 ± 0.01a | 552.45 ± 15.03b | 0.47 ± 0.03a | 33.33 ± 4.41a |
| SF | 6330.45 ± 181.72a | 1734.51 ± 78.11a | 0.80 ± 0.01a | 913.56 ± 83.26a | 0.37 ± 0.03a | 31.00 ± 3.06a | ||
| Maturity | CK | 21081.99 ± 1674.37b | 5006.32 ± 481.10b | 1.04 ± 0.01b | 224.62 ± 18.55a | 0.51 ± 0.02b | 17.33 ± 0.88b | |
| SF | 28004.47 ± 1850.68a | 6768.82 ± 93.92a | 1.13 ± 0.03a | 302.22 ± 27.11a | 0.68 ± 0.05a | 22.33 ± 0.88a | ||
| 2024 | Vigorous growth | CK | 5416.78 ± 146.94a | 1874.25 ± 37.5a | 0.79 ± 0.02b | 516.24 ± 15.00b | 0.91 ± 0.04a | 24.33 ± 0.91a |
| SF | 6213.25 ± 734.99a | 2186.28 ± 244.2a | 0.90 ± 0.02a | 739.64 ± 69.22a | 0.89 ± 0.05a | 27.17 ± 1.02a | ||
| Maturity | CK | 18157.92 ± 2079.32b | 5587.98 ± 699.78b | 0.99 ± 0.04b | 275.90 ± 34.22a | 0.59 ± 0.01b | 21.12 ± 1.28b | |
| SF | 27572.41 ± 2466.13a | 8473.53 ± 613.24a | 1.15 ± 0.02a | 319.03 ± 20.52a | 0.70 ± 0.02a | 26.86 ± 1.14a | ||
| ANOVA factor | Vigorous growth | Y | ns | ** | ** | ns | *** | ns |
| F | ns | ns | ** | *** | ns | ns | ||
| Y×F | ns | * | ns | ns | ns | ns | ||
| Maturity | Y | ns | ns | ns | ns | ns | ns | |
| F | ** | ns | ** | ns | *** | *** | ||
| Y×F | ns | *** | ns | ns | ns | ns | ||
Tobacco root traits at the vigorous growth and maturity stage in 2023 and 2024 after applying sesame cake fertilizer.
CK, conventional fertilization treatment; SF, sesame cake fertilizer treatment; RL, root length; RSA, root surface area; RMD, root mean diameter; RA, root activity; EC, electrical conductivity; CEC, cation exchange capacity. Y denotes the year, F represents the fertilizer, Y×F denote their interaction effect. Values are presented as mean ± standard error (n=3). Different letters denote significant differences according to Student’s t-test at P < 0.05. ns, *, **, and *** indicate non-significance, significance at P < 0.05, < 0.01 and < 0.001 probability levels, respectively.
Compared with the CK, the root dry weight in the SF treatment significantly increased at the vigorous growth and maturity stages in both years compared to CK (Figures 1A, C, E, F). The stem and leaf dry weights did not significantly differ after SF application, except that the SF treatment was significantly higher than CK during the vigorous growth in 2024 (Figures 1B, D–F).
Figure 1
3.3 Leaf areas, photosynthesis, and ultrastructure
For agronomic traits, the year significantly affected plant height, leaf width, and leaf area at the vigorous growth stage and only affected leaf length at the maturity stage. Fertilizer application significantly affected leaf length, width, and area at the vigorous growth stage and only significantly affected leaf width at the maturity stage. The interaction between year and fertilizer did not significantly affect plant height or leaf area at either growth stage (Table 3). The plant height and leaf length were not significantly different after the application of SF, while the leaf width and leaf area were significantly increased in the SF treatment at the vigorous growth stage compared to CK but did not differ at the maturity stage.
Table 3
| Year | Growth stage | Treatment | Agronomic traits | Photosynthesis traits | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| Height (cm) | Leaf length (cm) | Leaf width (cm) | Leaf area (cm²) | Pn (μmol CO2 m-2 s-1) | Gs (mol m-2 s-1) | Ci (ppm) | Tr (mmol H2O m-2 s-1) | |||
| 2023 | Vigorous growth | CK | 66.33 ± 1.86a | 67.33 ± 0.88a | 31.83 ± 0.17b | 1493.17 ± 24.06b | 8.12 ± 0.34b | 118.83 ± 6.32b | 253.33 ± 3.12b | 1.00 ± 0.04b |
| SF | 70.67 ± 2.33a | 70.00 ± 0.58a | 36.50 ± 0.29a | 1667.71 ± 13.98a | 11.02 ± 0.33a | 228.83 ± 12.99a | 294.67 ± 2.73a | 1.22 ± 0.02a | ||
| Maturity | CK | 112.67 ± 0.88a | 83.17 ± 2.52a | 31.50 ± 1.80a | 1783.87 ± 88.25a | 0.86 ± 0.05b | 21.33 ± 2.64a | 297.0 ± 9.89a | 0.31 ± 0.04a | |
| SF | 115.00 ± 1.73a | 87.67 ± 1.86a | 32.50 ± 2.18a | 1875.84 ± 106.07a | 1.30 ± 0.08a | 20.33 ± 2.06a | 264.0 ± 12.16a | 0.39 ± 0.05a | ||
| 2024 | Vigorous growth | CK | 86.67 ± 1.86a | 69.33 ± 1.67a | 28.33 ± 1.33b | 1247.85 ± 78.35b | 9.80 ± 0.39b | 123.20 ± 3.61b | 206.67 ± 8.19b | 1.39 ± 0.03b |
| SF | 87.33 ± 1.45a | 72.33 ± 1.20a | 30.67 ± 1.20a | 1497.63 ± 34.34a | 12.80 ± 0.24a | 180.00 ± 14.25a | 268.33 ± 2.30a | 3.45 ± 0.04a | ||
| Maturity | CK | 112.67 ± 2.08a | 89.33 ± 0.33a | 31.67 ± 0.33a | 1795.00 ± 22.97a | 0.78 ± 0.05b | 17.33 ± 3.31a | 282.17 ± 19.13a | 0.25 ± 0.06a | |
| SF | 116.33 ± 0.88a | 91.33 ± 0.67a | 30.67 ± 0.88a | 1777.87 ± 63.20a | 1.29 ± 0.06a | 27.33 ± 3.78a | 307.33 ± 8.46a | 0.37 ± 0.02a | ||
| ANOVA factor | Vigorous growth | Y | *** | ns | ** | ** | *** | ns | *** | *** |
| F | ns | * | *** | ** | *** | *** | *** | *** | ||
| Y×F | ns | ns | ns | ns | *** | ns | ns | * | ||
| Maturity | Y | ns | * | ns | ns | ns | ns | ns | ns | |
| F | ns | ns | *** | ns | * | ns | *** | ns | ||
| Y×F | ns | ns | ns | ns | ns | ns | ns | * | ||
Agronomic and photosynthesis traits of tobacco leave during the vigorous growth and maturity stage in 2023 and 2024 after applying sesame cake fertilizer.
CK, conventional fertilization treatment; SF, sesame cake fertilizer treatment; Pn, net photosynthetic rate; Gs, stomatal conductance; Ci, intercellular CO2 concentration; Tr, transpiration rate. Y denotes the year, F represents the fertilizer, Y×F denote their interaction effect. Values are presented as mean ± standard error (n=6). Different letters denote significant differences according to Student’s t-test at P < 0.05. ns, *, **, and *** indicate non-significance, significance at P < 0.05, < 0.01 and < 0.001 probability levels, respectively.
For photosynthesis traits, the year significantly affected the Pn, Ci, and Tr at the vigorous growth stage but did not affect any of the parameters at the maturity stage. Fertilizer significantly affected all parameters at the vigorous growth stage and affected both Pn and Ci at the maturity stage. The year and fertilizer interaction significantly affected both Pn and Tr at the vigorous growth stage, and only Tr at the maturity stage. Pn, Gs, Ci, and Tr were all significantly increased in the SF treatment compared to the CK treatment during vigorous growth in both years. Moreover, only Pn was significantly increased at the maturity stage in the two years; on average, across the year, the Pn in SF was 58.27% higher than in CK (Table 3).
The ultrastructure of the leaf chloroplasts was further observed; it showed an increase in chloroplast size with leaf aging, and eventually, the outer membrane ruptured to release mature starch granules (Figure 2). Compared to CK, at the vigorous growth stage, a higher number of chloroplasts and starch granules, more intact leaf cell walls, and thicker stacked lamellae were observed in the SF treatment. Furthermore, at the maturity stage, larger chloroplast volume, starch granule number, and laxer-stacked lamellae were observed in the SF treatment. These results indicate that tobacco leaves developed a well-organized chloroplast lamellar structure with an increased number of chloroplasts and starch granules after SF application.
Figure 2
3.4 Glandular trichomes and plastid pigment
Three types of tobacco leaf glandular trichomes, long glandular trichomes (LGTs), short glandular trichomes (SGTs), and one type of non-secretory glandular trichome, exist (). Aqueous rhodamine B staining of the secretory glandular trichomes (Song et al., 2022) revealed that the LGTs stain more intensely than SGTs in both treatments (Supplementary Figure S2). The recorded trichome numbers showed that the number of LGTs per millimeter in the SF treatment was significantly 27.27% higher than that in the CK treatment across the vigorous growth stage, and the numbers of LGTs, SGTs, and total glandular trichomes (TGTs) in the SF treatment were all significantly increased by 41.91%, 82.04%, and 56.17%, respectively, at the maturity stage compared to the CK treatment (Figures 3A, B; Supplementary Figure S3).
Figure 3
The total glandular trichome secretion content was significantly increased by 37.37% and 8.86% in the SF treatment at the vigorous growth and maturity stages, respectively, compared to that in the CK treatment (Figure 3C). At the vigorous growth stage, the secreted contents of cembratriene-diol, sucrose esters, and total alkanes in the SF treatment were significantly higher than those in the CK treatment. At the maturity stage, all secretory components were significantly higher in the SF treatment than in the CK treatment, except for total alkanes (Supplementary Table S4).
During the vigorous growth stage, both year and fertilizer application significantly affected the contents of chlorophyll a, b and carotenoid in tobacco leaves at the vigorous growth stage, and their interaction did not significantly affect all the plastid pigment contents (Table 4). The chlorophyll and carotenoid contents in the SF treatment significantly increased by 30.46% and 32.00% in 2023 and by 18.75% and 15.33% in 2024, respectively, compared to the respective CK treatment. At the maturity stage, only the year significantly affected carotenoid content.
Table 4
| Year | Growth stage | Treatment | Chlorophyll a (mg L-1) | Chlorophyll b (mg L-1) | Chlorophyll (mg g-1) | Carotenoid (mg L-1) |
|---|---|---|---|---|---|---|
| 2023 | Vigorous Growth Stage | CK | 10.73 ± 1.15a | 4.37 ± 0.37b | 1.51 ± 0.15b | 1.75 ± 0.20b |
| SF | 13.83 ± 0.40a | 5.85 ± 0.06a | 1.97 ± 0.04a | 2.31 ± 0.05a | ||
| Maturity Stage | CK | 3.29 ± 0.63a | 1.96 ± 0.29a | 0.52 ± 0.09a | 0.59 ± 0.02a | |
| SF | 3.49 ± 0.24a | 2.12 ± 0.12a | 0.56 ± 0.03a | 0.63 ± 0.05a | ||
| 2024 | Vigorous Growth Stage | CK | 18.32 ± 0.59b | 5.63 ± 0.13b | 2.40 ± 0.07b | 4.11 ± 0.19b |
| SF | 21.48 ± 0.31a | 7.04 ± 0.33a | 2.85 ± 0.06a | 4.74 ± 0.09a | ||
| Maturity Stage | CK | 3.27 ± 0.22a | 1.58 ± 0.07a | 0.49 ± 0.03a | 0.96 ± 0.07a | |
| SF | 2.17 ± 0.29b | 1.80 ± 0.11a | 0.40 ± 0.04a | 0.84 ± 0.09a | ||
| ANOVA factor | Vigorous growth | Y | *** | *** | *** | *** |
| F | ** | *** | *** | * | ||
| Y×F | ns | ns | ns | ns | ||
| Maturity | Y | ns | ns | ns | *** | |
| F | ns | ns | ns | ns | ||
| Y×F | ns | ns | ns | ns |
Plastid pigment content in the middle-leaf of tobacco in 2023 and 2024 after applying sesame cake fertilizer.
CK, conventional fertilization treatment; SF, sesame cake fertilizer treatment; Y denotes the year, F represents the fertilizer, Y×F denote their interaction effect. Values are presented as mean ± standard error (n=3). Different letters denote significant differences according to Student’s t-test at P < 0.05. ns, *, **, and *** indicate non-significance, significance at P < 0.05, < 0.01 and < 0.001 probability levels, respectively.
3.5 Volatile aroma substances in flue-cured tobacco
A total of 108 volatile substances were detected in flue-cured tobacco using HS-SPME-GC-MS (Figure 4A; Supplementary Figure S4, Supplementary Table S5), and 51 volatile substances with aromatic characteristics were identified (Figure 4B). The aroma substances were mainly divided into neutral aroma substances, including 17 carotenoid degradation products, 3 aromatic amino acid cleavage substances, 5 Maillard reaction components, 2 cembratriendid alkyl degradation products, and 1 chlorophyll degradation products, as well as other aroma substances, including 3 terpenes, 5 alkanes, 4 heterocyclic compounds, 3 organic acids, and 8 esters (Supplementary Table S6). Compared with CK treatment, all the aroma substances in the SF treatment were increased, except for the Maillard reaction components (Figure 4B).
Figure 4
OPLS-DA was used to analyze the effects of SF on the aromatic substances in flue-cured tobacco. The cumulative interpretation rates (R2X and R2Y) of the model were 0.629 and 0.995, respectively, and the overall predictive ability (Q2) was 0.915. The intercept between the Q2 regression curve and the Y-axis was less than zero, indicating no overfitting in this model (Figure 4C). The coordinates of the two treatments were clearly separated, indicating that the application of SF had a significant effect on the aroma composition of flue-cured tobacco.
Based on VIP>1 and P<0.05 of the first principal component in the OPLS-DA model, 19 differential aroma substances from the 51 substances were screened (Supplementary Table S7). The differential aromatic substances were upregulated in the SF treatment, except for diethyl phthalate (Figure 4D). Further, based on aroma classification and flavor description, seven differential volatile substances with the carotenoid degradation products as aroma precursors mainly provided sweet, floral, woody, and fruity flavors; one differential volatile substance with the aromatic amino acid cleavage products, Maillard reaction components, and cembratriendid alkyl degradation products as aroma precursors mainly provided sweet and mild flavors; three differential volatile substances with esters as aroma precursors mainly provided sweet, bitter, floral, and mild flavors; two differential volatile substances with terpenes and alkanes as aroma precursors mainly provided floral, woody, fruity, waxy, and mild flavors; and one differential volatile substance with the organic acid and heterocyclic compound as aroma precursors mainly provided sour and bitter flavors (Figure 4E; Supplementary Table S7; Shi and Liu, 1998; ).
3.6 Sensory quality and electronic nose in the flue-cured tobacco
The total sensory quality score of flue-cured tobacco in the SF treatment was significantly higher by 5.76% and 16.03% than in the CK treatment in 2023 and 2024, respectively. This was mainly caused by the significantly increased aroma quality and volume, which were increased by 13.7% and 17.91%, respectively, in the SF treatment compared to CK on average (Figure 5A).
Figure 5
The PCA analysis of the electronic nose showed that PC1 contributed 80.3% and PC2 contributed 12.7% of the total variance (93.0%) (Figure 5B). The SF treatment was located in the upper right quadrant and positively correlated with W3S, W6S, W5S, W1S, W2W, and W1W, which are related to the long-chain alkanes, hydrogen, broad-range oxynitride, methane, compounds with aromatic ring structures, sulfur-containing organic compounds, sulfur, organic compounds, and terpenes. The SF treatment was also located on the negative y-axis and negatively correlated with W2S. CK was located in the upper left quadrant and was negatively correlated with W1C, W5C, and W3C, which are associated with aromatic components, hydrocarbons, and ammonia (Supplementary Figure S5, Supplementary Table S2). Therefore, the odor distributions of the SF and CK treatments were significantly different as identified by the electronic nose.
Furthermore, the relationship between sensory quality and differential aroma substances was analyzed (Figure 5C). Aroma quality, aroma volume, and total sensory quality score were all significantly positively correlated with solanone, nicotine, and phytyl acetate, and significantly negatively correlated with diethyl phthalate. Aftertaste was significantly positively correlated with the majority of differential aroma substances, and offensive taste was significantly positively correlated only with delta-damascone, ethyl oleate, and safranal. Irritancy was not significantly associated with differential aroma substances.
3.7 Relationship between tobacco growth and sensory quality with the application of sesame cake fertilizer
Positive correlations were found between the different soil indicators, except for BD and D (Figure 6A). Differential aroma substances were significantly positively correlated to TC, TN, OM (p<0.01 and r≥0.5), TK, WC, R>0.25 (p<0.05 and r≥0.5), and TP (p<0.05 and 0.25<r<0.5). Root growth was significantly and positively correlated with OM (p<0.05, r≥0.5), TN and R>0.25 (p<0.05 and 0.25<r<0.5), and TC (p<0.05, r<0.25).
Figure 6
A SEM was constructed based on the indicators with significant correlations to further analyze the relationships (Figure 6B; Supplementary Figure S6). The model fit was assessed using the Chi-square test (χ²), with the following criteria indicating good structural equation modeling (SEM) fit: χ²/df < 3, CFI > 0.9, NFI > 0.9 and RMSEA < 0.06. Root CEC (0.13), leaf carotenoid (0.94) and cembratriene-diol (0.24) had a significant positive effect on differential aroma substances, and a further significant positive effect on sensory quality (0.34).
4 Discussion
4.1 Sesame cake fertilizer improved the soil quality and tobacco root growth
As sesame cake fertilizer (SF) was applied, the soil quality tended to increase (Table 1). This could be attributed to the rich nutrients in the SF, such as organic matter, proteins, and nutrient elements (). These components synergistically improve the soil structure and enrich the nutrient pool of the soil (Vann et al., 2013; ). The significantly increased soil OM after the application of SF promoted the formation of large soil aggregates, thereby enhancing soil aggregate structure and improving total porosity (Zou et al., 2018; Sun et al., 2024; Uddin et al., 2025). Ruan et al. (2025) found that long-term organic fertilization increased the porosity by 500–1000 μm compared to inorganic fertilization. The enhanced soil aggregation and biological activities, including earthworm activity and root proliferation, further increased soil organic carbon content, which aligns with the trends observed in our study (Table 1). These changes reduce soil bulk density, enhance soil aeration, and increase soil water storage capacity and fertility retention (; Uddin et al., 2025). Improved soil aeration also promoted tobacco root growth, resulting in increased RL, RSA, and RMD, especially during the maturity stage (Table 2; ). A similar response has been observed in rice () and maize (Song et al., 2020), indicating that organic fertilizers can promote root growth and expansion.
The increased content of soil nitrogen, phosphorus, and potassium after the application of SF provided additional nutrients for root growth, which was also verified by the higher CEC of the soil and roots in the SF treatment (Tables 1, 2; ). The higher ability to absorb and exchange nutrient ions in the soil further maintains soil fertility, and the increased root CEC and RA indicate enhanced nutrient absorption and utilization from the soil, promoting root growth and dry matter accumulation (Figure 1, ; Singh et al., 2015; ). Mature root systems maintain efficient material turnover even as growth decelerates, thereby sustaining functional effectiveness through optimized resource cycling (; ). The improvement in soil properties and effective plant growth is also derived from the increased biodiversity and interactions within microbial taxa, which boost soil metabolism (; ). Furthermore, the application of SF affects soil pH. The soil pH in Fujian Province is acidic and negatively affects tobacco growth; however, it significantly improved after the application of SF (Table 1; ; Wang et al., 2019). This phenomenon could be attributed to the release of base cations during the decomposition of SF, followed by displacement reactions with acidic components in the soil, effectively neutralizing the soil acidity (). Taken together, the improved soil quality after SF application established a healthy soil microenvironment that further promoted tobacco root growth and nutrient absorption and utilization.
4.2 Sesame cake fertilizer promotes the formation of aroma precursors in leaves
In tobacco leaves, plastid pigments are not only key to photosynthesis but also serve as essential aroma precursors, as they decompose and transform into aroma substances during the tobacco curing process (; ). The application of SF increased the chlorophyll and carotenoid content, which is consistent with previous studies on rice (Table 4, ). The increased plastid pigment content could be attributed to a higher number of chloroplasts and starch granules and a well-developed chloroplast lamellar structure in tobacco leaves in the SF treatment (Figure 2, Zhang et al., 2024c). The increased carotenoid can degrade more volatile compounds such as β-damascone and megastigmatrienone, which imparts sweet and floral flavor to tobacco leaves, further enhancing the aroma quality and volume of cured tobacco (Shi and Liu, 1998). Our results also confirmed the promotion of the degradation pathway and increased sensory quality (Figures 4E, 5). Additionally, higher plastid pigment content increases the leaf net photosynthetic rate and promotes photosynthate supply, thus increasing leaf area primarily by widening the leaf width (Table 3, ; ).
Tobacco aroma substances are positively correlated with glandular trichome density on the leaf, as trichome secretion serves as a critical source of aroma precursors (). In the present study, a significant increase in both long and short glandular trichome density was observed after SF application, particularly at the maturity stage (Figure 3B). This may be attributed to the well-developed chloroplast structure in tobacco glandular trichomes and thylakoid granule lamella, which promote glandular trichome development during early growth stages, while sufficient photosynthate supply enhances trichome secretion accumulation (; ). Accordingly, glandular trichome secretion was also considerably increased by SF application (Figure 3C). The higher cembratriene-diol content not only enhanced the antimicrobial, insecticidal, and antioxidant activities of the tobacco leaves but also degraded into valuable aroma substances, mainly solanone, during the curing process, corresponding with the increased solanone content observed in the flue-cured tobacco (Figures 3C, 4D; Wahlberg and Enzell, 1987).
4.3 Sesame cake fertilizer promoted sensory quality through increased aroma substance content in flue-cured tobacco
Aroma composition and volume are critical for the style characteristics and quality of flue-cured tobacco and are determined by the content and proportion of various aromatic substances (). The application of SF considerably increased the sensory quality, mainly owing to a higher score in aroma quality and volume (Figure 5). This increase resulted from the substantially enhanced aroma substances, which is consistent with previous studies on organic fertilizer application (Figure 4; Yan et al., 2020; Yan and Liu, 2020). In this study, SF mainly increased compounds derived from carotenoid as precursors. Carotenoids can undergo oxidative cleavage to form terpenoid compounds such as β-Ionone, β-Ionone epoxide, and megastigmatrienone, which are primarily responsible for imparting sweet, floral, woody, and fruity flavors (Shi and Liu, 1998). These substances exhibit low odor thresholds and minimal irritation, resulting in a higher efficiency in aroma contribution compared to ordinary volatile components (Shi et al., 2023). Additionally, the increased cembratriene-diol secreted by glandular trichomes facilitated oxidative degradation to generate solanone (Figure 3C; Supplementary Table S4, Wahlberg and Enzell, 1987; Xu et al., 2024). Solanone synergizes with β-Ionone, megastigmatrienone, and other terpenoid substances, effectively balancing the harshness of smoke and enhancing the persistence and subtlety of aroma (; Yin et al., 2019). The electronic nose analysis also reflected changes in volatile aroma substances; the sensors contributing the most were W1W (sulfurs, organic components, terpenes), W2W (aromatic ring components, organosulfur), and W5S (broad range, oxynitride; Figure 5B). These findings are consistent with the classification results of aroma substances detected by HS-SPME-GC-MS (Figure 4), wherein compounds derived from carotenoids and cembratriene-diol as precursors belonged to terpenes, while those derived from phenylalanine as precursors belonged to aromatic ring components (; Wang et al., 2018; ). Therefore, the e-nose results verified the aroma enhancement by SF, demonstrating that it is an effective method for rapidly estimating the impact of organic fertilizer on tobacco aroma (Yang et al., 2015).
4.4 Comprehensive analysis of the effect of sesame cake fertilizer on tobacco quality
Based on a comprehensive analysis of soil, root, and leaf growth, and aroma substances, there was a significant relationship between soil quality, root growth, and tobacco aroma substances (Figure 6A). This indicates that SF in the soil enhances soil fertility and improves soil physical and chemical properties, which directly or indirectly interact with the roots, further affecting secondary metabolite production in the leaves (Singh et al., 2015; ; Uddin et al., 2025). Subsequently, the well-developed chloroplast structure with an increased carotenoid content, as well as the higher cembratriene-diol content, accumulated more aroma substance precursors, which contributed to the increased aroma substance content in the flue-cured tobacco and ultimately improved sensory quality (Figure 6B; Yan et al., 2020; ; Song et al., 2022).
5 Conclusion
Sesame cake fertilizer application improved soil quality, increased soil and root cation exchange capacity, and promoted root growth and activity. It also enhanced leaf photosynthetic efficiency and plastid pigment content, enriched glandular trichome secretion, and increased the accumulation of key aroma precursors. Furthermore, the concentration of volatile aroma substances—including carotenoids and cembratriendid alkyl degradation products, esters, terpenes, and heterocyclic compounds—as well as the final enriched aroma quality and volume of flue-cured tobacco, were considerably enhanced. This study revealed that SF promoted plant growth and improved the aromatic quality of cured tobacco by enhancing the “soil-root-leaf-differential aroma substance” interaction, providing an effective organic fertilizer application strategy for flue-cured tobacco production.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
CQ: Data curation, Visualization, Investigation, Software, Writing – original draft. JL: Project administration, Supervision, Conceptualization, Writing – original draft, Investigation, Resources. XZ: Software, Writing – original draft, Formal analysis. WY: Visualization, Writing – original draft, Project administration. XyZ: Supervision, Project administration, Writing – original draft. ZL: Supervision, Project administration, Writing – review & editing. DT: Supervision, Writing – original draft, Investigation. SD: Supervision, Writing – review & editing. JC: Writing – review & editing, Supervision. YW: Funding acquisition, Formal analysis, Writing – review & editing, Project administration, Conceptualization, Supervision.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This study was supported by the project of Guangdong China Tobacco Industry Co., LTD., scientific research project (2023440000340013), Longyan Company of Fujian Provincial Tobacco Corporation scientific research project (2022350824240020, 2023350800240041), Guangdong Basic and Applied Basic Research Foundation (2022A1515110353), Double first-class discipline promotion project (2021B10564001).
Acknowledgments
The authors want to thank the staff of South China Agricultural University, China Tobacco Guangdong Industrial Company Limited, and Longyan Company of Fujian Provincial Tobacco Corporation for their excellent work.
Conflict of interest
Authors JL, XZ, ZL, XyZ were employed by China Tobacco Guangdong Industrial Company Limited. Author DT was employed by Longyan Company of Fujian Provincial Tobacco Corporation.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors declare that this study received funding from Guangdong China Tobacco Industry Co., LTD. and from Longyan Company of Fujian Provincial Tobacco Corporation. The funders had the following involvement in the study: study design, data collection, decision to publish.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1654657/full#supplementary-material
Glossary
- AN
Available Nitrogen
- AP
Available Phosphorus
- AV
Aroma Volume
- BD
Bulk Density
- Car
Carotenoid
- Cem
Cembratriene-diol
- Ch
Chloroplast
- Chla
Chlorophyll a
- Chlb
Chlorophyll b
- CEC
Cation Exchange Capacity
- Ci
Intercellular CO2 Concentration
- CK
Conventional Fertilizer (Control)
- D
Fractal Dimension
- EC
Electrical Conductivity
- GMD
Geometric Mean Diameter
- Gs
Stomatal Conductance
- HS-SPME
Headspace Solid-Phase Microextraction
- LA
Leaf Area
- LGTs
Long Glandular Trichomes
- MS
Mass Spectrometry
- MWD
Mean Weight Diameter
- NAS
Neutral Aroma Substances
- OM
Organic Matter
- P
Porosity
- Pn
Net photosynthetic rate
- RA
Root Activity
- RMD
Root Mean Diameter
- RL
Root Length
- RSA
Root Surface Area
- SEM
Structural Equation Modeling
- SGTs
Short Glandular Trichomes
- SF
Sesame Cake Fertilizer
- SPME
Solid-Phase Microextraction
- TC
Total Carbon
- TN
Total Nitrogen
- TP
Total Phosphorus
- TK
Total Potassium
- Tr
Transpiration Rate
- VIP
Variable Importance in Projection
References
1
AppauA.DropeJ.WitoelarF.ChavezJ. J.LencuchaR. (2019). Why do farmers grow tobacco? A qualitative exploration of farmers perspectives in Indonesia and Philippines. Int. J. Enviro. Res. Public Health16, 2330. doi: 10.3390/ijerph16132330
2
BanožićM.JokićS.AčkarĐ.BlažićM.ŠubarićD. (2020). Carbohydrates—key players in tobacco aroma formation and quality determination. Molecules25, 1734. doi: 10.3390/molecules25071734
3
BrarB. S.SinghK.DheriG. S.KumarB. (2013). Carbon sequestration and soil carbon pools in a rice–wheat cropping system: effect of long-term use of inorganic fertilizers and organic manure. Soil Till. Res.128, 30–36. doi: 10.1016/j.still.2012.10.001
4
BuschM.HainA. S. A. R. (2002). Functional analysis of the early steps of carotenoid biosynthesis in tobacco. Plant Physiol.128, 439–453. doi: 10.1104/pp.128.2.439
5
CaiK.XiangZ.PanW.ZhaoH.RenZ.LeiB.et al. (2013). Identification and quantitation of glycosidically bound aroma compounds in three tobacco types by gas chromatography–mass spectrometry. J. Chromatogr. A.1311, 149–156. doi: 10.1016/j.chroma.2013.08.051
6
ChenJ.LiY.HeX.JiaoF.ZouC. (2021b). Influences of different curing methods on chemical compositions in different types of tobaccos. Ind. Crop Prod.167, 113534. doi: 10.1016/j.indcrop.2021.113534
7
ChenY.JiangZ.OuJ.LiuF.CaiG.TanK.et al. (2024). Nitrogen substitution practice improves soil quality of red soil (Ultisols) in South China by affecting soil properties and microbial community composition. Soil Till. Res.240, 106089. doi: 10.1016/j.still.2024.106089
8
ChenY.LvX.QinY.ZhangD.ZhangC.SongZ.et al. (2023). Effects of different botanical oil meal mixed with cow manure organic fertilizers on soil microbial community and function and tobacco yield and quality. Front. Microbiol.14. doi: 10.3389/fmicb.2023.1191059
9
ChenD.WangM.WangG.ZhouY.YangX.LiJ.et al. (2022). Functional organic fertilizers can alleviate tobacco (Nicotiana tabacum L.) continuous cropping obstacle via ameliorating soil physicochemical properties and bacterial community structure. Front. Bioeng. Biotechnol.10. doi: 10.3389/fbioe.2022.1023693
10
ChenM.ZhangS.LiuL.WuL.DingX. (2021a). Combined organic amendments and mineral fertilizer application increase rice yield by improving soil structure, P availability and root growth in saline-alkaline soil. Soil Till. Res.212, 105060. doi: 10.1016/j.still.2021.105060
11
DasP. P.SinghK. R.NagpureG.MansooriA.SinghR. P.GhaziI. A.et al. (2022). Plant-soil-microbes: a tripartite interaction for nutrient acquisition and better plant growth for sustainable agricultural practices. Environ. Res.214, 113821. doi: 10.1016/j.envres.2022.113821
12
EifediyiE. K.AhamefuleH. E.RemisonS. U.AliyuT. H.AkanbiN. (2017). Effects of neem seed cake and NPK fertilizer on the growth and yield of sesame (Sesamum indicum L.). Cercetari Agronomice Moldova.50, 57–72. doi: 10.1515/cerce-2017-0015
13
FengW.Sánchez-RodríguezA. R.BilyeraN.WangJ.WangX.HanY.et al. (2024). Mechanisms of biochar-based organic fertilizers enhancing maize yield on a Chinese Chernozem: Root traits, soil quality and soil microorganisms. Environ. Technol. Innov.36, 103756. doi: 10.1016/j.eti.2024.103756
14
GaoZ.LiJ. Z.LiuS. T.ChenY. (2023). Within-leaf chloroplasts and nitrogen allocation to thylakoids in relation to photosynthesis during grain filling in maize. Plant Physiol.Bioch.196, 830–840. doi: 10.1016/j.plaphy.2023.02.034
15
GengZ.HeP.GaoH.LiuJ.QiuJ.CaiB. (2023). Aroma precursors of cigars from different tobacco parts and origins, and their correlations with sensory characteristics. Front. Plant Sci.14. doi: 10.3389/fpls.2023.1264739
16
GörgüçA.BircanC.MehmetF. M. (2019). Sesame bran as an unexploited by-product: Effect of enzyme and ultrasound-assisted extraction on the recovery of protein and antioxidant compounds. Food Chem.283, 637–645. doi: 10.1016/j.foodchem.2019.01.077
17
GouZ.ZhaiZ.ZhangQ.LiY.FengX.ZhangY.et al. (2024). Characteristics of organic carbon in tobacco-growing soils with different pH. Bangladesh J. Bot.53, 745–755. doi: 10.3329/bjb.v53i30.76614
18
HeintzeS. G. (1960). Studies on cation-exchange capacities of roots. Plant Soil.13, 365–383. doi: 10.1007/bf01394648
19
HuW.ZhangY.RongX.ZhouX.FeiJ.PengJ.et al. (2024). Biochar and organic fertilizer applications enhance soil functional microbial abundance and agroecosystem multifunctionality. Biochar6, 3. doi: 10.1007/s42773-023-00296-w
20
HuW.ZhangY.XiangminR.FeiJ.PengJ.LuoG. (2023). Coupling amendment of biochar and organic fertilizers increases maize yield and phosphorus uptake by regulating soil phosphatase activity and phosphorus-acquiring microbiota. Agr. Ecosyst. Enviro.355, 11. doi: 10.1016/j.agee.2023.108582
21
JabborovaD.MaH.Bellingrath-KimuraS. D.WirthS. (2021). Impacts of biochar on basil (Ocimum basilicum) growth, root morphological traits, plant biochemical and physiological properties and soil enzymatic activities. Sci. Hortic.290, 110518. doi: 10.1016/j.scienta.2021.110518
22
JacobsD. F.TimmerV. R. (2005). Fertilizer-induced changes in rhizosphere electrical conductivity: relation to forest tree seedling root system growth and function. New For.30, 147–166. doi: 10.1007/s11056-005-6572-z
23
KeeneC. K.WagnerG. J. (1985). Direct demonstration of duvatrienediol biosynthesis in glandular heads of tobacco trichomes. Plant Physiol.79, 1026–1032. doi: 10.1104/pp.79.4.1026
24
KudoyarovaG.VeselovaS.HartungW.FarhutdinovR.VeselovD.SharipovaG. (2011). Involvement of root ABA and hydraulic conductivity in the control of water relations in wheat plants exposed to increased evaporative demand. Planta233, 87–94. doi: 10.1007/s00425-010-1286-7
25
LiJ.MaZ.DaiH.LiH.QiuJ.PangX. (2024b). Application of PLSR in correlating sensory and chemical properties of middle flue-cured tobacco leaves with honey-sweet and burnt flavour. Heliyon10, e29547. doi: 10.1016/j.heliyon.2024.e29547
26
LiY.YangS.DongY. (2024a). Response of agronomic and physiological traits of flue-cured tobacco (Nicotiana tabacum L.) to plant stem cell soil amendments. Physiol. Mol. Plant Pathol.131, 102290. doi: 10.1016/j.pmpp.2024.102290
27
LiangZ.JiH.ZhangH.WengM.CuiH. (2009). Morphology and structure of chloroplast in flue-cured tobacco trichomes after applying fertilizer. Acta Bot.Boreal.Occident. Sin.29, 0291–0295. doi: 10.3321/j.issn:1000-4025.2009.02.013
28
LiuL.LiH.ZhuS.GaoY.ZhengX.XuY. (2021b). The response of agronomic characters and rice yield to organic fertilization in subtropical China: A three-level meta-analysis. Field Crop Res.263, 108049. doi: 10.1016/j.fcr.2020.108049
29
LiuX.LiuH.ZhangY.LiuC.LiuY.LiZ.et al. (2023). Organic amendments alter microbiota assembly to stimulate soil metabolism for improving soil quality in wheat-maize rotation system. J. Environ. Manage.339, 117927. doi: 10.1016/j.jenvman.2023.117927
30
LiuJ.ShuA.SongW.ShiW.GaoZ. (2021a). Long-term organic fertilizer substitution increases rice yield by improving soil properties and regulating soil bacteria. Geoderma404, 115287. doi: 10.1016/j.geoderma.2021.115287
31
LiuA.YuanK.LiQ.LiuS.LiY.TaoM.et al. (2022). Metabolomics and proteomics revealed the synthesis difference of aroma precursors in tobacco leaves at various growth stages. Plant Physiol. Biochem.192, 308–319. doi: 10.1016/j.plaphy.2022.10.016
32
LuoG.LiL.FrimanV. P.GuoJ.GuoS.ShenQ.et al. (2018). Organic amendments increase crop yields by improving microbe-mediated soil functioning of agroecosystems: A meta-analysis. Soil Biol. Biochem.124, 105–115. doi: 10.1016/j.soilbio.2018.06.002
33
LuoD. S.WangB.QiaoX. Y. (2019). Explanation of national regionalization of leaves style of flue-cured tobacco. Acta Tabacaria Sin.25, 1–9. doi: 10.16472/j.Chinatobacco.2019.218
34
MeybergM.KrohnS.BrümmerB.KristenU. (1991). Ultrastructure and secretion of glandular trichomes of tobacco leaves. Flora185, 357–363. doi: 10.1016/s0367-2530(17)30495-4
35
MiaoY.StewartB. A.ZhangF. (2011). Long-term experiments for sustainable nutrient management in China. A review. Agron. Sustain. Dev.31, 397–414. doi: 10.1051/agro/2010034
36
MulvaneyR. L.KhanS. A.EllsworthT. R. (2009). Synthetic nitrogen fertilizers deplete soil nitrogen: a global dilemma for sustainable cereal production. J. Environ. Qual.38, 2295–2314. doi: 10.2134/jeq2008.0527
37
NouskaC.DeligeorgakiM.KyrkouC.MichaelidouA. M.MoschakisT.BiliaderisC. G.et al. (2024). Structural and physicochemical properties of sesame cake protein isolates obtained by different extraction methods. Food Hydrocoll.151, 109757. doi: 10.1016/j.foodhyd.2024.109757
38
OguchiR. (2010). Does the photosynthetic light-acclimation need change in leaf anatomy? Plant Cell Environ.26, 505–512. doi: 10.1046/j.1365-3040.2003.00981.x
39
QiC.ZhangX.LiD.TongD.YeW.LiH.et al. (2025). Effects of equal replacement of chemical fertilisers by sesame and peanut cake fertilisers on carbon and nitrogen metabolism and nitrogen utilisation of flue-cured tobacco. Jiangsu Agric. Sci. 53, 85–93. doi: 10.15889/j.issn.1002-1302.2025.12.012
40
QianX. J.XiaoJ.HongY. F.ShiG. Q.TanX.ZengW. L.et al. (2019). Spatial and temporal variability of soil pH in Longyan tobacco planting areas. Chin. J. Trop. Crops.40, 2061–2067. doi: 10.3969/j.issn.1000-2561.2019.10.021
41
QinZ.ChangY. L.ChenZ. M.WangY. G.FanW.GuL. B.et al. (2024). A novel strategy for preparing lignan-rich sesame oil from cold-pressed sesame seed cake by combining enzyme-assisted treatment and subcritical fluid extraction. Ind. Crop Prod.218, 11. doi: 10.1016/j.indcrop.2024.119041
42
RamosF. T.DoresE. F. D. C.WeberO. L. D. S.BeberD. C.CampeloJ. H.JrMaiaJ. C. D. S. (2018). Soil organic matter doubles the cation exchange capacity of tropical soil under no‐till farming in Brazil. J. Sci. Food. Agric.98, 3595–3602. doi: 10.1002/jsfa.8881
43
ReichertJ. M.PellegriniA.RodriguesM. F. (2019). Tobacco growth, yield and quality affected by soil constraints on steeplands. Ind. Crop Prod.128, 512–526. doi: 10.1016/j.indcrop.2018.11.037
44
RowenE.TookerJ. F.BlubaughC. K. (2019). Managing fertility with animal waste to promote arthropod pest suppression. Biol. Control.134, 130–140. doi: 10.1016/j.biocontrol.2019.04.012
45
RuanR.ZhangZ.LanT.WangY.LiW.ChenH.et al. (2025). The role of soil pore structure on nitrate release from soil organic matter and applied fertilizer under three fertilization regimes. Soil Till. Res.248, 106396. doi: 10.1016/j.still.2024.106396
46
ShiH. Z.LiuG. S. (1998). Tobacco aromatics (China: China Agriculture Press).
47
ShiJ.WuW.ZhangY.BaldermannS.PengQ.WangJ.et al. (2023). Comprehensive analysis of carotenoids constituents in purple-coloured leaves and carotenoid-derived aroma differences after processing into green, black, and white tea. Lwt173, 114286. doi: 10.1016/j.lwt.2022.114286
48
SinghP.RaiR. K.SumanA.SrivastavaT. K.SinghK. P.AryaN.et al. (2015). Soil-root interface changes in sugarcane plant and ratoon crops under subtropical conditions: implications for dry-matter accumulation. Commun. Soil Sci. Plant Anal.46, 454–475. doi: 10.1080/00103624.2014.997385
49
SongX.RazaviB. S.LudwigB.ZamanianK.ZangH.KuzyakovY.et al. (2020). Combined biochar and nitrogen application stimulates enzyme activity and root plasticity. Sci. Total Environ.735, 139393. doi: 10.1016/j.scitotenv.2020.139393
50
SongZ.WangP.ChenX.PengY.CaiB.SongJ.et al. (2022). Melatonin alleviates cadmium toxicity and abiotic stress by promoting glandular trichome development and antioxidant capacity in Nicotiana tabacum. Ecotoxicol. Environ. Saf.236, 113437. doi: 10.1016/j.ecoenv.2022.113437
51
SunX.NiuL.ZhangM.ZhangH.LiuH.ZhaoM.et al. (2024). Application of carbon-based nutrient fertilizer improved soil fertility and seed yield of Paeonia ostii ‘Feng Dan’. Ind. Crop Prod.212, 118348. doi: 10.1016/j.indcrop.2024.118348
52
ThurièsL.PansuM.Larré-LarrouyM. C.FellerC. (2002). Biochemical composition and mineralization kinetics of organic inputs in a sandy soil. Soil Biol. Biochem.34, 239–250. doi: 10.1016/s0038-0717(01)00178-x
53
Tor-NgernP.Chart-AsaC.ChanthornW.RodtassanaC.HasselquistN. J. (2021). Variation of leaf-level gas exchange rates and leaf functional traits of dominant trees across three successional stages in a southeast Asian tropical forest. For. Ecol. Manage.489, 119101. doi: 10.1016/j.foreco.2021.119101
54
UddinM. K.SahaB. K.WongV. N.PattiA. F. (2025). Organo-mineral fertilizer to sustain soil health and crop yield for reducing environmental impact: a comprehensive review. Eur. J. Agron.162, 127433. doi: 10.1016/j.eja.2024.127433
55
VannM. C.FisherL. R.JordanD. L.DavidS. W.HardyD. H.StewartA. M. (2013). Potassium rate and application effect on flue-cured tobacco. Agron. J.105, 304. doi: 10.2134/agronj2012.0259
56
WahlbergI.EnzellC. R. (1987). Tobacco isoprenoids. Nat. Prod. Rep.4, 237. doi: 10.1039/np9870400237
57
WangB.DengX.WangR.ZongguoX.TongW.MaE.et al. (2024b). Bio-organic substitution in tobacco (Nicotiana tabacum L) cultivation: Optimum strategy to lower carbon footprint and boost net ecosystem economic benefit. J. Environ. Manage.370, 122654. doi: 10.1016/j.jenvman.2024.122654
58
WangS.SunN.ZhangS.LongdozB.WellensJ.MeersmansJ.et al. (2024a). Soil organic carbon storage impacts on crop yields in rice-based cropping systems under different long-term fertilisation. Eur. J. Agron.161, 127357. doi: 10.1016/j.eja.2024.127357
59
WangJ.WangH.FuY.HuangT.WangX. (2021). Genetic variance and transcriptional regulation modulate terpenoid biosynthesis in trichomes of nicotiana tabacum under drought. Ind. Crop Prod.167, 113501. doi: 10.1016/j.indcrop.2021.113501
60
WangH.XuJ.LiuX.ZhangD.LiL.LiW.et al. (2019). Effects of long-term application of organic fertilizer on improving organic matter content and retarding acidity in red soil from China. Soil Tillage Res.195, 104382. doi: 10.1016/j.still.2019.104382
61
WangJ. A.YangG. H.LiC. X. (2018). Zonal distribution of neutral aroma components in flue-cured tobacco leaves. Phytochem. Lett.24, 125–130. doi: 10.1016/j.phytol.2018.01.014
62
WeiW.YanY.CaoJ.ChristieP.ZhangF.FanM. (2016). Effects of combined application of organic amendments and fertilizers on crop yield and soil organic matter: an integrated analysis of long-term experiments. Agric. Ecosyst. Environ.225, 86–92. doi: 10.1016/j.agee.2016.04.004
63
WenS.SunL.ZhangS.ChenZ.ChenR.LiZ.et al. (2023). The formation mechanism of aroma quality of green and yellow teas based on GC-MS/MS metabolomics. Food Res. Int.172, 113137. doi: 10.1016/j.foodres.2023.113137
64
World Health Organization (WHO) (2017). Tobacco and its environmental impact: an overview. Geneva: World Health Organization. Available online at: http://apps.who.int/iris.
65
WuQ.WuX.ZhangX.JiangC.XiaoB.ZhangY.et al. (2014). Mapping of two white stem genes in tetraploid common tobacco (Nicotiana tabacum L.). Mol. Breed.34, 1065–1074. doi: 10.1007/s11032-014-0097-0
66
WuX. P.ZhongX. M.QinY. Q.LiuS. (2006). Effects of proportional application of sesame seed cake fertilizers and chemical fertilizers on the aroma quality of flue-cured tobacco leaves. Acta Agronomica Sinica.32, 1554–1559. doi: 10.3321/j.issn:0496-3490.2006.10.021
67
XinL.TangM.ZhangL.HuangW.WangX.GaoY. (2024). Effects of saline-fresh water rotation irrigation on photosynthetic characteristics and leaf ultrastructure of tomato plants in a greenhouse. Agric. Water Manage.292, 108671. doi: 10.1016/j.agwat.2024.108671
68
XuM.DuY.HouX.ZhangZ.YanN. (2024). Chemical structures, biosynthesis, bioactivities, and utilisation values for the diterpenes produced in tobacco trichomes. Phytochemistry223, 114117. doi: 10.1016/j.phytochem.2024.114117
69
YanS.LiuG. (2020). Effect of increasing soil carbon content on tobacco aroma and soil microorganisms. Phytochem. Lett.36, 42–48. doi: 10.1016/j.phytol.2020.01.011
70
YanS.WangP.CaiX.WangC.ZwietenL. V.WangH.et al. (2025). Biochar-based fertilizer enhanced tobacco yield and quality by improving soil quality and soil microbial community. Environ. Technol. Innov.37, 103964. doi: 10.1016/j.eti.2024.103964
71
YanS.ZhaoJ.RenT.LiuG. (2020). Correlation between soil microbial communities and tobacco aroma in the presence of different fertilizers. Ind. Crop Prod.151, 112454. doi: 10.1016/j.indcrop.2020.112454
72
YangJ.AnH.LiJ.ZhengJ.LongJ. (2015). Study on analysis of volatile components in cigarette made through different drying processing by electronic nose. Chem. Res. Application.27, 246–250. doi: 10.3969/j.issn.1004-1656.2015.03.002
73
YangJ.RenY.JiaM.HuangS.GuoT.LiuB.et al. (2025). Improving soil quality and crop yield of fluvo-aquic soils through long-term organic-inorganic fertilizer combination: promoting microbial community optimization and nutrient utilization. Environ. Technol. Innov.37, 104050. doi: 10.1016/j.eti.2025.104050
74
YinF.KarangwaE.SongS.DuhoranimanaE.LinS.CuiH.et al. (2019). Contribution of tobacco composition compounds to characteristic aroma of Chinese faint-scent cigarettes through chromatography analysis and partial least squares regression. J. Chromatogr. B.1105, 217–227. doi: 10.1016/j.jchromb.2018.12.001
75
ZhangW.CaoJ.LiZ.LiQ.LaiX.SunL.et al. (2021). HS-SPME and GC/MS volatile component analysis of yinghong no. 9 dark tea during the pile fermentation process. Food Chem.357, 129654. doi: 10.1016/j.foodchem.2021.129654
76
ZhangM.GuoD.WuG.HanP.ShiY.ZhengT.et al. (2024b). Analysis of volatile compound metabolic profiles during the fermentation of filler tobacco leaves through integrated E-nose, GC–MS, GC-IMS, and sensory evaluation. J. Chromatogr. A.1737, 465472. doi: 10.1016/j.chroma.2024.465472
77
ZhangK.HanX.FuY.KhanZ.ZhangB.BiJ.et al. (2024c). Biochar coating promoted rice growth under drought stress through modulating photosynthetic apparatus, chloroplast ultrastructure, stomatal traits and ROS homeostasis. Plant Physiol. Bioch.216, 109145. doi: 10.1016/j.plaphy.2024.109145
78
ZhangS. T.SongX. N.LiN.ZhangK.LiuG. S.LiX. D.et al. (2018a). Influence of high-carbon basal fertiliser on the structure and composition of a soil microbial community under tobacco cultivation. Res. Microbiol.169, 115–126. doi: 10.1016/j.resmic.2017.10.004
79
ZhangB.ZhangZ.BaiX.LiL.WuJ.LiuY.et al. (2024a). Long-term rice cultivation enhances root development and yields by improving the structural properties of soil aggregates in saline–alkaline environments. Environ. Technol. Innov.36, 103848. doi: 10.1016/j.eti.2024.103848
80
ZhangH. Y.ZhangS. T.YangY. X.JiaH. F.CuiH. (2018b). Metabolic flux engineering of cembratrienol production in both the glandular trichome and leaf mesophyll in nicotiana tabacum. Plant Cell Physiol.3, 566–574. doi: 10.1093/pcp/pcy004
81
ZhaoS.WangW.ChenX.GaoY.WuX.DingM.et al. (2023). Graphene oxide affected root growth, anatomy, and nutrient uptake in alfalfa. Ecotoxicol. Environ. Saf.250, 114483. doi: 10.1016/j.ecoenv.2022.114483
82
ZouC. M.LiY.HuangW.ZhaoG. K.PuG.r.SuJ. E.et al. (2018). Rotation and manure amendment increase soil macro-aggregates and associated carbon and nitrogen stocks in flue-cured tobacco production. Geoderma325, 49–58. doi: 10.1016/j.geoderma.2018.03.017
Summary
Keywords
flue-cured tobacco, sesame cake fertilizer, root growth, leaf glandular trichomes, aromatic substances, sensory quality
Citation
Qi C, Lei J, Ye W, Zhang X, Li Z, Zhong X, Tong D, Deng S, Chen J and Wang Y (2025) Sesame cake fertilizer improves tobacco aroma quality by boosting root growth and leaf aroma precursor formation. Front. Plant Sci. 16:1654657. doi: 10.3389/fpls.2025.1654657
Received
26 June 2025
Accepted
30 July 2025
Published
20 August 2025
Volume
16 - 2025
Edited by
Anoop Kumar Srivastava, Central Citrus Research Institute (ICAR), India
Reviewed by
Heiplanmi Rymbai, ICAR Research Complex for NEH Region, India
Dahe Qiao, Guizhou Academy of Agricultural Sciences (CAAS), China
Updates
Copyright
© 2025 Qi, Lei, Ye, Zhang, Li, Zhong, Tong, Deng, Chen and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jianjun Chen, chenjianjun@scau.edu.cn; Yuanyuan Wang, wangyuanyuan@scau.edu.cn
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.