- Safety and Quality Institute of Agricultural Products, Heilongjiang Academy of Agricultural Sciences, Harbin, Heilongjiang, China
Introduction: Tricyclazole is widely employed as a pesticide for controlling rice blast. While effective, it simultaneously acts as a stressor on rice plants.
Methods: Rice seedlings were treated with tricyclazole for 7 days, and the root and shoot samples were collected for analysis. Four biomarkers of oxidative stress, including SOD, CAT, POD, and MDA were measured to evaluate the physiological responses. Integrated transcriptome and metabolome analysis was performed to reveal the mechanisms underlying the tricyclazole-induced stress response.
Results: Tricyclazole significantly increased the levels of SOD, CAT, POD, and MDA in the root and shoot of rice. A variety of differential metabolites and differentially expressed genes (DEGs) were identified, which exhibited diverse functionalities. Flavonoid biosynthesis, glutathione metabolism, and phenylpropanoid biosynthesis were three important metabolic pathways in response to tricyclazole. These pathways encompassed many important metabolites and regulatory genes involved in stress responses. In addition, D-mannose, hispidulin-8-C-glucoside, dicumarol, and 4-hydroxyderricin in the root, and S-Adenosyl-L-methionine, nicotinate D-ribonucleoside, luteolin-7-O-(3'-O-coumaroyl) sophorotrioside, and D-Arabinono-1,4-lactone in the shoot were key metabolites involved in self-detoxification and resistance development of rice against tricyclazole. At the molecular level, these metabolites could be regulated by various DEGs involved in the systemic acquired resistance and response to toxic substances.
Discussion: A coordinated molecular and metabolic reprogramming triggered by tricyclazole enables the rice plant to mitigate stress and develop resistance.
Introduction
Rice (Oryza sativa L.), belongs to the genus Oryza within the family Poaceae, is one of the most important cereal crops (Kandil et al., 2022). Globally, rice ranks as the third most produced agricultural commodity and serves as a staple food for approximately half of the human population (Gross and Zhao, 2014; Nazish et al., 2022). In agricultural production, adversity stress is the most important factor affecting crop yield and quality. Rice is particularly susceptible to a range of abiotic stresses, such as drought, cold, salinity, and heavy metal toxicity, as well as biotic stresses caused by pathogens, such as fungi, bacteria, and viruses (Ahmad, 2022). Given that stress generates effects at the phenotypic, physiological, biochemical, and molecular levels, studying the mechanisms of stress response has become a key focus of agricultural research.
Rice blast is a widespread and destructive disease affecting cultivated rice globally, posing a serious threat to production and food security (Fernandez et al., 2021). Under favourable environmental conditions, the yield loss attributed to rice blast can reach 30–50% in major rice-growing regions (Ashkani et al., 2015). This devastating fungal disease is primarily caused by Magnaporthe oryzae, a hemibiotrophic pathogen that infects rice by initially maintaining host cells in the early phase, and subsequently kills cells in the late stage (Syauqi et al., 2022). During infection and colonisation, Magnaporthe oryzae mainly targets immune responses, destroys host defenses, and enhances the susceptibility of rice plants to fungal invasion (Fernandez and Orth, 2018). So far, pesticide application is one of the effective strategies to control rice blast, among which tricyclazole (8-methyl-[1,2,4]triazolo[3,4-b][1,3]benzothiazole) is one of the most common ones. As a triazole fungicide, tricyclazole is effective to control rice blast by inhibiting spore germination and adherent spore formation (Li et al., 2023; Liu et al., 2024c). However, tricyclazole, especially at high doses, also causes certain physiological toxicity to rice, thereby affecting rice growth and quality. Elucidating the underlying mechanisms of rice’s response to tricyclazole provides a theoretical basis for developing resistant varieties and formulating targeted strategies to reduce phytotoxicity, thereby supporting sustainable rice yield and quality.
Pesticides are a form of stress in plants, and their absorption in leaves or roots can lead to adverse outcomes, such as tissue damage, growth inhibition, and quality decline. In terms of underlying mechanisms, pesticides often trigger molecular-level changes and in turn, disrupt crucial metabolic processes (Liu et al., 2024a). Transcriptome and metabolome sequencing offer great possibilities for elucidating the intrinsic mechanisms of plants’ responses to pesticide stress, among which the former identifies key regulatory genes, while the latter captures dynamic changes in metabolic profiles (Lv et al., 2024). So far, some studies have confirmed relevant applications in rice. As reported, butachlor affects carbohydrate metabolism, and the key differentially expressed genes (DEGs) are involved in starch and sucrose metabolism (Liu and Zhu, 2020). Chlorpyrifos affects amino acid metabolism, and the key DEGs are involved in aspartate and glutamate metabolism (Liu and Zhu, 2020). Ametryn affects the metabolisms of carbohydrates, amino acids, hormones and phenylpropanoids in rice, and relevant genes encoding enzymes for the catabolism of ametryn are up-regulated (Qiao et al., 2023). ZhiNengCong elevates the levels of carbohydrates, vitamins, nucleotides, and amino acids in rice with bacterial leaf streak, and up-regulates genes involving resistance, growth, and sugar/amino acids/lipid metabolisms (Lu et al., 2025). In the present study, transcriptome and metabolome sequencing were performed in rice plants treated with tricyclazole. The metabolomic and transcriptomic changes were analysed in the root and shoot samples. This study aims to elucidate the mechanisms of stress defense and resistance establishment in response to tricyclazole, and to identify pivotal metabolites and regulatory genes involved.
Materials and methods
Sample processing
Rice (Wuyou 4) was used for the hydroponic cultivation of seedlings. When the seedlings reached the stage of one-leaf and one-heart (about 10 days after germination), the medium containing 20 mg/L tricyclazole was refreshed to trigger a stress response (P group). The seedlings cultured in conventional medium without tricyclazole were used as an internal check (CK group). The CK and P groups both have three biological replicates, with 50 seedlings per replicate. After 7 days of continuous cultivation, the seedlings were collected and then divided into the root and shoot parts. Four biomarkers of oxidative stress, including SOD, CAT, POD, and MDA were measured in the root and shoot samples as previously described (Chen et al., 2025). For the above physiological assay, three independent experiments were conducted, with three biological replicates in each group. Another part of the root and shoot samples (a total of 4 groups; N = 3 per group) was used for transcriptome and metabolome sequencing as follows.
Transcriptome sequencing
RNAs were isolated from samples of different groups using CTAB (Sangon, Shanghai, China) combined with pBIOZOL (Bioer, Hangzhou, China), quantified by a fluorescence analyser (Qubit 4.0, Thermo Fisher Scientific, Waltham, MA, USA), and assessed by a high-throughput biological fragment analyser (Qsep400, Bioptic, Taiwan, China). Based on qualified RNAs, libraries were constructed using an RNA-seq library prep kit (Vazyme, Nanjing, China), quantified by a fluorescence analyser (Qubit 4.0), and assessed by a fragment analyser (Agilent 2100, Santa Clara, CA, USA). Finally, a DNA nanoball (about 300 copies) was prepared from a qualified library and sequenced on the MGI platform (BGI Genomics, Shenzhen, China).
Transcriptome analysis
Raw data were processed using fastp v0.23.2, and the reads with adapters, with N > 10% of the base number of the reads, and with Q ≤ 25 (low quality) for more than 50% of the base number of the reads were filtered out. The number of raw reads, clean reads, Q20, Q30, and GC content of each sample were calculated. Clean reads were mapped to the reference genome (IRGSP-1.0_genome) using HISAT2 v2.2.1. Fragments Per Kilobase of transcript per Million fragments mapped (FPKM) was calculated to reflect gene expression levels using featureCounts v2.0.3. Principal component analysis (PCA) was performed using prcomp (R v4.1.2) to reveal the inter-group and intra-group variabilities. For two-group comparison, DEGs were determined using DESeq2 v1.22.1, and the threshold was log2FC ≥ 1 and False Discovery Rate (FDR) < 0.05 (Benjamini-Hochberg test). Based on the isolated DEGs, clustered heatmaps and volcano plots were generated. Besides, DEGs were aligned with KOG databases using diamond v2.0.9 to obtain annotation results.
Metabolome sequencing
The samples were frozen in a lyophiliser (Scientz-100F, Scientz, Ningbo, China) and ground into powder using a grinder (MM 400, Retsch, Haan, Germany). Subsequently, 30 mg powder was added into 1500 μL pre-cooled 70% methanolic aqueous internal standard extract. The mixture was vortexed for 30 s with 30 min intervals for six times in total. After 3 min of centrifugation at 12,000 rpm, the supernatant was collected and filtered through a microporous membrane (0.22 μm). The filtrate was used for metabolome sequencing using an UPLC-ESI-MS/MS system (SCIEX, Shanghai, China).
Metabolome analysis
The MS data were qualitatively analysed based on the MWDB database (Myway, Shenzhen, China). The isotope signals, repeated signals containing K+, Na+, NH4+, and repeated signals themselves are fragment ions of other substances with higher molecular weight were removed during analysis. In addition, the metabolites were quantified using the multiple reaction monitoring mode. PCA was performed using prcomp (R v4.1.2) to reveal the inter-group and intra-group variabilities. For two-group comparison, differential metabolites were determined by VIP > 1, fold change ≥ 2/≤ 0.5, and P < 0.05. VIP values were generated from the OPLS-DA (Orthogonal Partial Least Squares-Discriminant Analysis) model using MetaboAnalystR. A permutation test (200 permutations) was performed to avoid overfitting (Supplementary Figure S1). Based on the isolated differential metabolites, clustered heatmaps, volcano plots, and scatter plot were generated. Besides, Z-score standardization was used to reveal the distribution of important differential metabolites between groups.
Integrated analysis of metabolome and transcriptome
Based on the DEGs and differential metabolites, the KEGG pathways annotated by both omics were identified. If both the metabolomic and transcriptomic data have enrichment results in MetMap pathway (Myway), the co-annotated MetMap was added to the results of KEGG. If no MetMap pathways were co-annotated, only the KEGG results were displayed. Besides, DEGs and differential metabolites with |Pearson correlation coefficient| > 0.8 and P < 0.05 were used to draw cornetwork plots using Cytoscape v3.8.
Targeted validation of key metabolites
D-mannose and S-adenosyl-L-methionine were selected for targeted validation based on the commercial availability of reliable reference standards and the existence of well-established quantitative detection methods.
The concentration of D-mannose in the root was measured by gas chromatography-mass spectrometry (GC-MS). Simply, freeze-dried samples were ground and 20 mg powder samples were immersed in 500 µL methanol:isopropanol:water (3:3:2, v/v/v). After vortexing, sonication, and centrifugation, 50 µL supernatants were mixed with 20 μL internal standards. For accharide analysis, the internal standards were [²H6]-D-ribose (IR-25184; isoreag, Shanghai, China) and adonitol (R017943, Shanghai, China), prepared as a 1000 µg/mL solution in methanol. After dried under a nitrogen stream, the residue was derivatized by reaction with 100 μL methoxyamine hydrochloride in pyridine at 37 °C for 2 h, followed by reaction with 100 μL BSTFA at 37 °C for 30 min. Finally, 50 μL derivatization solution was diluted to 1 mL with n-hexane, filtered through a 0.22 μm membrane, and used for GC-MS analysis on the Agilent 8890-5977B platform (Agilent, Santa Clara, CA, USA). Based on the quantitative ion (m/z, 160) and qualitative ion (m/z, 319), a selected ion monitoring detection mode was employed to ensure the resolution of the target compound. The retention time (~7.54 min) and ion ratio (~102) of D-mannitol in the test samples matched those in the mixed standards, confirming accurate identification. For quality control, three quality control samples (pooled samples) were inserted between the mixed standards and the test samples during loading. The method performance characteristics were as follows: inter-day precision (n = 3), 9.34%; intra-day precision (n = 3), 9.10%; recovery, 93.07%. The calibration curve was y = 0.535674x – 2.162182E-004.
The concentration of S-adenosyl-L-methionine in the shoot was measured by liquid chromatography-mass spectrometry (LC-MS). Simply, 50 mg thawed samples were extracted in 500 μL methanol containing 10 μL internal standard for 30 min at -20 °C. The internal standard was [2H3]-SAM (SAM-d3) (485-80-3; SCRSTANDARD, Shanghai, China), prepared as a 500 ng/mL solution in methanol. After centrifugation, 250 μL supernatant was collected for LC-MS analysis on QTRAP 6500+ System (Sciex, Framingham, MA, USA). The ion mode was [M+H]+. Quantitative analysis was performed using the multiple reaction monitoring mode. The retention time of S-adenosyl-L-methionine (~5.68 min) in the test samples matched that in the mixed standards. For quality control, three quality control samples (pooled samples) were inserted between the mixed standards and the test samples during loading. The method performance characteristics were as follows: inter-day precision (N = 3), 6.36%; intra-day precision (N = 3), 5.87%; recovery, 128.80%. The calibration curve was y = 0.00943x - 0.03064.
Statistical analyses
Physiological and verification experiments were performed with three independent biological replicates (N = 3 per group). Data were expressed as mean ± standard deviation. Statistical analysis was conducted using GraphPad Prism 7.0 (GraphPad, San Diego, CA, USA). Differences between two groups were assessed by independent sample t-test. A p-value < 0.05 was considered statistically significant.
Results
The physiological responses of rice plant to the treatment of tricyclazole
After the treatment of tricyclazole, the physiological changes in rice plant were evaluated. As shown in Figure 1A, the fresh weights of the root and shoot were both slightly reduced by the treatment of tricyclazole, but the differences were not significant. Subsequently, the oxidative stress was analysed based on four biomarkers. In the root, the levels of SOD, CAT, POD, and MDA were significantly higher in the P-R group than in the CK-R group (P < 0.05). Consistently, the levels of SOD, CAT, POD, and MDA in the shoot were also significantly higher in the P-S group than in the CK-S group (P < 0.05) (Figures 1B-E).
Figure 1. Physiological changes in the root and shoot of rice plant after the treatment of tricyclazole. (A) Fresh weight. (B) SOD. (C) CAT. (D) POD. (E) MDA. *P < 0.05, **P < 0.01 vs. CK-R; #P < 0.05, ##P < 0.01 vs. CK-S.
Changes of transcriptomes in the root and shoot of rice plant after the treatment of tricyclazole
High-throughput sequencing was performed to reveal the transcriptomic differences in rice plant after the treatment of tricyclazole. Both the RNA samples and the sequencing data are of high quality (Supplementary Table S1 and 2). PCA displayed a high intra-group repeatability, evidenced by the fact that all intra-group replicates were clustered within a confidence ellipse. In addition to a slight overlap between the CK-R and P-R (treatment) groups, the four groups exhibited separated confidence ellipses, which indicated potential differences at the transcription level (Figure 2A). Cluster heatmap showed the presence of a large number of DEGs between the CK and P groups in the root and shoot samples (Figure 2B). According to the criteria of log2FC ≥ 1 and FDR < 0.05, 138 up-regulated and 713 down-regulated DEGs were determined in the root, and 585 up-regulated and 457 down-regulated DEGs were determined in the shoot between the CK and P groups (Figures 2C, D). KOG annotation revealed a diverse range of functional roles in these DEGs, with distinct profiles observed between the root and shoot (Figures 2E, F).
Figure 2. Transcriptomic differences in the root and shoot of rice plant after the treatment of tricyclazole. (A) PCA; different colored dots represent samples from different groups. (B) Cluster heatmap; red represents up-regulation and blue represents down-regulation. (C) Volcanic diagram in the root. (D) Volcanic diagram in the shoot; green points represent down-regulated DEGs, red points represent up-regulated DEGs, and gray points represent genes without significant differences. (E) Column diagram of KOG annotation in the root. (F) Column diagram of KOG annotation in the shoot.
Changes of metabolomes in the root and shoot of rice plant after the treatment of tricyclazole
UPLC-MS/MS was further performed to reveal the metabolomic differences in rice plant after the treatment of tricyclazole. A total of 3058 metabolites were detected before filtering. PCA showed that the intra-group replicates were all within a confidence ellipse, indicating a good repeatability. A clear separation among the four groups indicates the potential differences at the metabolic level (Figure 3A). Between the CK and P groups in the root and shoot, a large number of differential metabolites were observed on the cluster heatmap (Figure 3B). According to the criteria of VIP > 1, fold change ≥ 2/≤ 0.5, and P < 0.05, 63 up-regulated and 363 down-regulated metabolites were determined in the root, and 141 up-regulated and 340 down-regulated metabolites were determined in the shoot between the CK and P groups (Figures 3C, D). These differential metabolites have a wide range of categories, and some metabolites within the same category may also have differences between the root and shoot (Figures 3E, F).
Figure 3. Metabolomic differences in the root and shoot of rice plant after the treatment of tricyclazole. (A) PCA; different colored dots represent samples from different groups. (B) Cluster heatmap; red represents up-regulation and green represents down-regulation. (C) Volcanic diagram in the root. (D) Volcanic diagram in the shoot; green points represent down-regulated metabolites, red points represent up-regulated metabolites, and gray points represent metabolites without significant differences. (E) Scatter plot in the root. (F) Scatter plot in the shoot; each point represents a metabolite, and different colours indicate different classifications.
Functional enrichment of differential metabolites and DEGs
Bubble plots were constructed using KEGG pathways co-enriched by two omics datasets. The top 25 terms were listed in Figure 4. In the root, the differential metabolites and DEGs were mainly co-enriched in the flavonoid biosynthesis, flavone and flavonol biosynthesis, phenylpropanoid biosynthesis, glutathione metabolism, biosynthesis of secondary metabolites, biosynthesis of amino acids, etc. In addition to the flavonoid biosynthesis, glutathione metabolism, phenylpropanoid biosynthesis, and biosynthesis of secondary metabolites, the differential metabolites and DEGs in the shoot were also co-enriched in zeatin biosynthesis, luteolin aglycones biosynthesis, diterpenoid biosynthesis, tryptophan metabolism, chrysoeriol aglycones biosynthesis, etc.
Figure 4. Bubble plots of co-enriched KEGG by differential metabolites and DEGs in the root and shoot. Red, yellow, and blue represent a change in the significance of enrichment from high, medium, to low; the shape of bubbles represents different omics, and the size of bubbles represents the number of differential metabolites and DEGs.
Key metabolic pathways in response to tricyclazole
Based on the key metabolic pathways (flavonoid biosynthesis, glutathione metabolism, phenylpropanoid biosynthesis) enriched in both the root and shoot samples, a metabolic diagram was constructed to display relevant metabolites and regulatory genes (Figure 5). In flavonoid biosynthesis, tricyclazole increased sakuranetin and apigenin levels in the root, which might be mediated by Os04t0175900–01 and Os04t0101400-01, respectively. In phenylpropanoid biosynthesis, tricyclazole induced the accumulation of 4-Coumaryl alcohol in the root, possibly mediated by Os04t0612700–01 and Os04t0229100-01. In addition, 4-Coumaroylshikimate and p-Coumaroyl quinic acid were differential metabolites involved in both the flavonoid biosynthesis and phenylpropanoid biosynthesis, and their accumulation in the root might be positively regulated by Os06t0185300-01. Besides, glutathione disulfide is a key metabolite in glutathione metabolism. Under tricyclazole treatment, the up-regulation of Os02t0664000-01, Os02t0664050-00, and Os07t0638400–01 might contribute to the accumulation of glutathione disulfide in the shoot.
Figure 5. A schematic diagram depicting key metabolic pathways, including flavonoid biosynthesis, glutathione metabolism, and phenylpropanoid biosynthesis. Genes and metabolites are colored red or green to indicate up- or down-regulation, based on their fold changes in the root and shoot.
Identification of key metabolites and relevant genes involved in stress response
Referring to Z-score, the top 50 differential metabolites in the root and shoot were displayed in Figure 6. These metabolites may be important participants in rice’s response to tricyclazole. For example, D-mannose, hispidulin-8-C-glucoside, dicumarol, and 4-hydroxyderricin were all increased by the treatment of tricyclazole in the root. In the shoot, tricyclazole treatment increased S-Adenosyl-L-methionine, nicotinate D-ribonucleoside, luteolin-7-O-(3’-O-coumaroyl)sophorotrioside, and D-Arabinono-1,4-lactone. Targeted GC-MS subsequently validated that the concentration of D-mannose in the root was significantly higher in the P-R group than in the CK-R group (65.89 ± 0.52 vs. 41.72 ± 1.49 μg/g, P < 0.001). Targeted LC-MS validated that the concentration of S-Adenosyl-L-methionine in the shoot was significantly higher in the P-S group than in the CK-S group (72.64 ± 1.97 vs. 54.20 ± 1.08 μg/g, P < 0.01) (Figure 7). Potential molecular networks for these metabolites were further constructed. As shown in Figures 8, nine and eight DEGs implicated in the systemic acquired resistance were revealed in the root and shoot, respectively. Additionally, 13 and 20 DEGs implicated in the response to toxic substances were revealed in the root and shoot samples, respectively. These DEGs may be involved in the response to tricyclazole by mediating stress response and resistance development.
Figure 6. The Z-value plots of differential metabolites (top 50) in the root and shoot. Green and yellow dots represent the values of metabolites in the P and CK groups, respectively.
Figure 7. Validation of D-mannose and S-adenosyl-L-methionine concentrations by targeted GC -MS and LC-MS, respectively. ***P < 0.001 vs. CK-R; ##P < 0.01 vs. CK-S.
Figure 8. Correlation network plots constructed by DEGs implicated in the systemic acquired resistance and the response to toxic substances, and the screened metabolites.
Discussion
Tricyclazole is a pesticide widely used for the treatment of rice blast (Li et al., 2023; Liu et al., 2024c). Although tricyclazole effectively inhibits Magnaporthe oryzae infection, it simultaneously triggers defense responses in rice, imposing a form of stress on the plant. In the present study, the treatment of tricyclazole significantly increased the levels of SOD, CAT, POD, and MDA in the root and shoot of rice plants. These results indicate that tricyclazole induces a state of oxidative stress in rice plants, which in turn activates a systemic antioxidant defense response. Subsequently, metabolome and transcriptome sequencing were performed to reveal the mechanisms underlying the responses of rice plant to tricyclazole. As a result, a variety of differential metabolites and DEGs were isolated from the root and shoot samples. These factors may be key responders to tricyclazole and participate in a coordinated defense response in rice plants.
In order to determine which metabolic pathways are affected by tricyclazole, co-enrichment analysis was performed on the differential metabolites and DEGs. Integrated analysis showed that the flavonoid biosynthesis, glutathione metabolism, and phenylpropanoid biosynthesis were key metabolic pathways enriched in both the root and shoot samples. As reported, phenylpropanoids are a group of secondary metabolites that are associated with tolerances to drought stress (Shim et al., 2025), soluble aluminum (Wang et al., 2023), and ultraviolet B radiation (Liu et al., 2024b) in rice. As a branch of phenylpropanoid, flavonoids are also important in response to abiotic and biotic stresses in plants (Ahmed et al., 2021; Chen et al., 2023). In rice, flavonoids help to combat drought and ultraviolet radiation stress (Jan et al., 2022), salt and heat stress (Jan et al., 2021), and thiocyanate pollution (Tian et al., 2024). Similarly, glutathione is also an essential metabolite in plants that mediates the responses to both abiotic and biotic stresses by virtue of its roles in detoxification, antioxidant defense, and redox homeostasis (Noctor et al., 2012). Evidence has determined that glutathione helps rice resist complete submergence, chromium toxicity, and salt stress. The above evidence supports that these metabolic pathways may play a central role in rice’s response to tricyclazole. Furthermore, to better elucidate the involved metabolites and genes, a metabolic diagram was constructed. This diagram highlights several important metabolites and regulatory genes potentially implicated in tricyclazole stress response. For example, naringenin 7-O-methyltransferase encoded by Os04t0175900–01 increases the accumulation of sakuranetin, a crucial defense metabolite against environmental stress (Miyamoto, 2025). Os04t0101400-01, which encodes a flavone synthase II, is involved in the synthesis of apigenin, contributing to antioxidant activity and stress tolerance in plants (Zhu et al., 2024). Os04t0612700–01 and Os04t0229100–01 both encode cinnamyl-alcohol dehydrogenase, a key enzyme involved in stress defense in plants, are associated with the accumulation of 4-Coumaryl alcohol (Park et al., 2018; Bagniewska-Zadworna et al., 2014). In addition, glutathione peroxidase encoded by Os02t0664000-01, Os02t0664050-00, and Os07t0638400–01 can eliminate peroxides by oxidizing reduced glutathione to glutathione disulfide (Averill-Bates, 2023). To sum up, tricyclazole may activate the stress response of rice plants by affecting the flavonoid, glutathione, and phenylpropanoid metabolic pathways at both the transcriptome and metabolome levels.
Further screening was conducted on potential metabolites involved in stress response. Some metabolites may be critical in responding to tricyclazole-induced stress. For example, Hispidulin 8-C-(2- Osyringoyl-b-glucopyranoside) is involved in antioxidant defence of Haberlea rhodopensis to cold acclimation and relevant desiccation (Georgieva et al., 2021). Ascorbic acid is an important participant in scavenging ROS and blocking photoinhibition in plants under stress (Celi et al., 2023; Yang et al., 2023a). Nicotinate D-ribonucleoside, the ribonucleoside form of nicotinic acid, serves as a key intermediate in the salvage pathway of nicotinamide adenine dinucleotide (NAD) biosynthesis in plants. NAD-mediated ROS plays a significant role in controlling plant growth and conferring stress resistance (Hong et al., 2020). Elevated levels of nicotinate D-ribonucleoside have been observed in Cu(II)-exposed yeast (Farrés et al., 2016), microbial inoculated rhizospheric soil (Lai et al., 2024), and downy mildew-infected sunflower (Zhu et al., 2023). Luteolin-7-O-(3’-O-coumaroyl) isophorotrioside may also be implicated in stress responses by acting as a derivative of luteolin. Luteolin is involved in plants’ resistance to Sclerotinia sclerotiorum (Robison et al., 2018), Mythimna separata (Qi et al., 2016), nitrogen deficiency (Kováčik and Klejdus, 2014), and Cd (Zhu et al., 2024). D-Arabinono-1,4-lactone serves as a direct precursor in the biosynthetic pathway of L-ascorbic acid (AsA). As one of the most pivotal antioxidants in plants, AsA contributes to plant resistance against a range of abiotic stresses (Wu et al., 2024b). Bao et al. reported that D-Arabinono-1,4-lactone enhances the tolerances to drought and chilling in tobacco and stylo (Bao et al., 2016). In the present study, consistent changes in the above metabolites were observed by metabolomics analysis in tricyclazole-treated rice, either in the root or in the shoot. The changes in these metabolites reflect a metabolic reprogramming, likely contributing to the activation of defence mechanisms. Besides, dicumarol and 4-hydroxyderricin are associated with antibacterial activity. 4-Hydroxyderricin is a chalcone that exhibits an antimicrobial activity against various pathogens, including viruses, bacteria, fungi, and protozoa (Caesar et al., 2018; Adhikari et al., 2025). Beyond its antimicrobial property, dicumarol also exhibits insecticidal activity and the ability to induce phytoalexin production in plants (Li et al., 2014; Rehman et al., 2017). These two metabolites may contribute to the antifungal property of tricyclazole against Magnaporthe oryzae in rice plants. It should be noted that the functional inferences for the above metabolites are inferred from studies in other plant species or under other stress conditions. Although these comparisons provide reasonable hypotheses for our findings, the specific physiological functions of these metabolites in tricyclazole-treated rice remain putative and require further validation through targeted studies.
D-mannose and S-adenosyl-L-methionine are also two important metabolites potentially involved in the response to tricyclazole. Consistent with the metabolomics results, verification experiments further confirmed the accumulation of D-mannose in the root and S-adenosyl-L-methionine in the shoot following tricyclazole treatment. D-mannose is a monosaccharide that elevated under salt stress in tomato (Wu et al., 2024a) and under cold stress in Dendrobium huoshanense (Wu et al., 2023). The role of D-Mannose in plant stress tolerance may be closely associated with ascorbic acid biosynthesis. As a key precursor, D-mannose is utilized in the D-mannose/L-galactose pathway to synthesize ascorbic acid (Mishra et al., 2024). The resulting ascorbic acid can scavenge ROS and enhance plant stress tolerance by regulating the antioxidant pool (Duan et al., 2016; Chen et al., 2020). The accumulation of D-mannose under tricyclazole treatment suggests the activation of an antioxidant defense mechanism centered on ascorbic acid biosynthesis. On the other hand, S-Adenosyl-L-methionine serves as a methyl donor in methylation reactions, which participates in the stress responses to Ralstonia solanacearum and alkali (Gong et al., 2014; Wang et al., 2019), drought (Peleg et al., 2011), freezing (Kou et al., 2018), and salt (Ma et al., 2017). Mechanistically, the involvement of S-adenosyl-L-methionine in stress responses is primarily linked to the ethylene, polyamine, and glutathione biosynthetic pathways, all of which are activated in plants under stress conditions (Yang et al., 2023b; Noctor et al., 2012). Specifically, 1-aminocyclopropane-1-carboxylic acid synthase catalyzes the cleavage of S-Adenosyl-L-methionine to produce 1-aminocyclopropane-1-carboxylic acid, which is then converted into ethylene by 1-aminocyclopropane-1-carboxylic acid oxidase. Ethylene can coordinate plant growth and stress management by interacting with auxin, abscisic acid, salicylic acid, jasmonic acid, and gibberellins (Shukla, 2025). After decarboxylation, S-adenosyl-L-methionine also acts as an aminopropyl donor in polyamine synthesis. Under stress, polyamines positively influence plant growth, ion balance, water retention, photosynthesis, and ROS production (Napieraj et al., 2023). As discussed earlier, glutathione metabolism is a key metabolic pathway in response to tricyclazole. In transsulfuration, the sulfur atom from S-adenosyl-L-methionine can be converted into cysteine through a series of enzymatic steps, making it a precursor for glutathione (Lu, 2000). Thus, the accumulation of S-adenosyl-L-methionine may enhance tricyclazole tolerance by activating the ethylene, polyamine, and glutathione biosynthesis.
Moreover, the molecular basis of the above differential metabolites was further explored by constructing correlation networks. The results showed that these metabolites were linked to a variety of DEGs implicated in the systemic acquired resistance and response to toxic substances. For example, NPR1 is an activator of systemic acquired resistance involved in plant immunity (Zavaliev et al., 2020). NPR1 plays a crucial role in enhancing plant resistance to biotic and abiotic stresses by modulating the interaction between salicylic acid and other defense and growth hormones (Zavaliev and Dong, 2023). Salicylic acid glucosyltransferase also acts an important regulator in plant defence response through its regulatory role in salicylic acid (Kobayashi et al., 2020). In the molecular networks, Os03t0667100-02, which encodes NPR1, was correlated with D-mannose and S-adenosyl-L-methionine. Os09t0517900–01 and Os04t0206500-01, both encoding pathogen-inducible salicylic acid glucosyltransferase, were linked to D-mannose and S-adenosyl-L-methionine, respectively. These connections may be direct, or indirectly bridged via salicylic acid, with GDP-D-mannose/ascorbic acid and ethylene as key linkers (Miao et al., 2024; Mortimer et al., 2013; Shukla, 2025). Additionally, some other DEGs in the networks may also be key regulators in stress responses, such as Os02t0314600-01 (aspartyl protease family protein) and Os02t0483500-01 (BAHD acyltransferase), which correlate with D-mannose, as well as Os07t0673200-01 (E3 ubiquitin-protein ligase BAH), Os01t0175850-00 (UDP-glucosyltransferase), and Os01t0369700-03 (glutathione S-transferase), which are linked to S-adenosyl-L-methionine. These defense-related genes may act through modulation of metabolic networks to ultimately coordinate the plant’s defense response.
The possible mechanisms by which tricyclazole induces defence responses in rice plant are as follows. The root, acting as the direct site of tricyclazole perception, rapidly initiates transcriptional reprogramming. This reprogramming triggers a series of metabolic events, characterised by the accumulation of D-mannose, hispidulin-8-C-glucoside, dicumarol, 4-hydroxyderricin, etc. Among these, the first two metabolites exhibit antioxidant potential, and the latter two exhibit antimicrobial activity. Their enrichment suggests that the root system reinforces its chemical defence barrier to mitigate the potential adverse effects of tricyclazole and also to counteract potential secondary infections. As the response is transduced to the shoot via systemic signaling, it continues to trigger coordinated transcriptional reprogramming and associated metabolic adjustments. The accumulation of key metabolites, including S-adenosyl-L-methionine, nicotinate D-ribonucleoside, flavonoid luteolin-7-O-(3’-O-coumaroyl)sophorotrioside, and D-arabino-1,4-lactone, collectively illustrates the metabolic investment of the shoot system in antioxidant defence mechanisms. These metabolic adaptations may contribute to the adaptive detoxification of tricyclazole in rice plant, thereby helping to counteract the increase in oxidative stress (Figure 9).
This study remain has some limitations. First, the use of a single concentration and exposure time may not fully capture the dynamic effects of tricyclazole. Second, validation experiments relying on a single time point may miss critical early peaks or dynamic patterns and therefore limit mechanistic interpretation. Third, the proposed mechanisms remain preliminary, and future studies are needed to confirm the specific roles of key genes and metabolites identified.
Conclusion
In conclusion, tricyclazole activated oxidative stress in the roots and shoots of rice plants. Flavonoid biosynthesis, glutathione metabolism, and phenylpropanoid biosynthesis were the main metabolic pathways in response to tricyclazole. D-mannose, hispidulin-8-C-glucoside, dicumarol, 4-hydroxyderricin (root), S-Adenosyl-L-methionine, nicotinate D-ribonucleoside, luteolin-7-O-(3’-O-coumaroyl) sophorotrioside, and D-Arabinono-1,4-lactone (shoot) might be key metabolites involved in self-detoxification and resistance development in tricyclazole-treated rice. In addition, transcriptomic reprogramming orchestrated the metabolic changes under tricyclazole treatment, and some associated DEGs were implicated in the systemic acquired resistance and response to toxic substances.
Data availability statement
The data presented in this study are deposited in the SRA repository, accession number: PRJNA1359757 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1359757).
Author contributions
WH: Conceptualization, Visualization, Funding acquisition, Writing – review & editing, Formal analysis, Writing – original draft, Investigation, Methodology, Supervision, Data curation. DS: Investigation, Data curation, Writing – review & editing. AC: Data curation, Methodology, Writing – review & editing. GC: Writing – review & editing, Investigation, Data curation. FL: Methodology, Writing – review & editing, Investigation. XZ: Data curation, Formal analysis, Writing – review & editing. JD: Writing – review & editing, Formal analysis, Data curation. JL: Investigation, Writing – review & editing, Visualization. HR: Writing – review & editing, Data curation, Visualization. WG: Writing – review & editing, Investigation, Formal analysis. BL: Data curation, Conceptualization, Supervision, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study was funded by Research Business Expense Project of Heilongjiang Provincial Research Institute (CZKYF2025-1-B017), and Key Laboratory of Quality and Safety of Cereals and Their Products, State Administration for Market Regulation (GWZL2025-002).
Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declared that generative AI was not 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.2026.1723722/full#supplementary-material
References
Adhikari, S., Nath, P., Deb, V. K., Das, N., Banerjee, A., Pathak, S., et al. (2025). Pharmacological potential of natural chalcones: a recent studies and future perspective. Front. Pharmacol. 16, 1570385. doi: 10.3389/fphar.2025.1570385
Ahmad, M. (2022). Genomics and transcriptomics to protect rice (Oryza sativa. L.) from abiotic stressors: -pathways to achieving zero hunger. Front. Plant Sci. 13, 1002596. doi: 10.3389/fpls.2022.1002596
Ahmed, U., Rao, M. J., Qi, C., Xie, Q., Noushahi, H. A., Yaseen, M., et al. (2021). Expression profiling of flavonoid biosynthesis genes and secondary metabolites accumulation in populus under drought stress. Molecules 26, 5546. doi: 10.3390/molecules26185546
Ashkani, S., Rafii, M. Y., Shabanimofrad, M., Miah, G., Sahebi, M., Azizi, P., et al. (2015). Molecular breeding strategy and challenges towards improvement of blast disease resistance in rice crop. Front. Plant Sci. 6, 886. doi: 10.3389/fpls.2015.00886
Averill-Bates, D. A. (2023). The antioxidant glutathione. Vitam. Horm. 121, 109–141. doi: 10.1016/bs.vh.2022.09.002
Bagniewska-Zadworna, A., Barakat, A., Lakomy, P., Smoliński, D. J., and Zadworny, M. (2014). Lignin and lignans in plant defence: insight from expression profiling of cinnamyl alcohol dehydrogenase genes during development and following fungal infection in Populus. Plant Sci. 229, 111–121. doi: 10.1016/j.plantsci.2014.08.015
Bao, G., Zhuo, C., Qian, C., Xiao, T., Guo, Z., and Lu, S. (2016). Co-expression of NCED and ALO improves vitamin C level and tolerance to drought and chilling in transgenic tobacco and stylo plants. Plant Biotechnol. J. 14, 206–214. doi: 10.1111/pbi.12374
Caesar, L. K., Kellogg, J. J., Kvalheim, O. M., Cech, R. A., and Cech, N. B. (2018). Integration of biochemometrics and molecular networking to identify antimicrobials in angelica keiskei. Planta Med. 84, 721–728. doi: 10.1055/a-0590-5223
Celi, G. E. A., Gratão, P. L., Lanza, M., and Reis, A. R. D. (2023). Physiological and biochemical roles of ascorbic acid on mitigation of abiotic stresses in plants. Plant Physiol. Biochem. 202, 16. doi: 10.1016/j.plaphy.2023.107970
Chen, W., Hu, T., Ye, J., Wang, B., Liu, G., Wang, Y., et al. (2020). A CCAAT-binding factor, SlNFYA10, negatively regulates ascorbate accumulation by modulating the D-mannose/L-galactose pathway in tomato. Hortic. Res. 7, 200. doi: 10.1038/s41438-020-00418-6
Chen, S., Wang, X., Cheng, Y., Gao, H., and Chen, X. (2023). A review of classification, biosynthesis, biological activities and potential applications of flavonoids. Molecules 28, 4982. doi: 10.3390/molecules28134982
Chen, G., Xu, J., Wang, J., Huang, W., Liu, F., Dong, J., et al. (2025). Integrated analysis of transcriptome and metabolome revealed clomazone biodegradation in maize seedlings. Ecotoxicol. Environ. Saf. 298, 14. doi: 10.1016/j.ecoenv.2025.118302
Duan, W., Ren, J., Li, Y., Liu, T., Song, X., Chen, Z., et al. (2016). Conservation and Expression Patterns Divergence of Ascorbic Acid d-mannose/l-galactose Pathway Genes in Brassica rapa. Front. Plant Sci. 7, 778. doi: 10.3389/fpls.2016.00778
Farrés, M., Piña, B., and Tauler, R. (2016). LC-MS based metabolomics and chemometrics study of the toxic effects of copper on Saccharomyces cerevisiae. Metallomics 8, 790–798. doi: 10.1039/C6MT00021E
Fernandez, J., Lopez, V., Kinch, L., Pfeifer, M. A., Gray, H., Garcia, N., et al. (2021). Role of two metacaspases in development and pathogenicity of the rice blast fungus magnaporthe oryzae. mBio 12, 03471–03420. doi: 10.1128/mBio.03471-20
Fernandez, J. and Orth, K. (2018). Rise of a cereal killer: the biology of magnaporthe oryzae biotrophic growth. Trends Microbiol. 26, 582–597. doi: 10.1016/j.tim.2017.12.007
Georgieva, K., Mihailova, G., Gigova, L., Dagnon, S., Simova-Stoilova, L., and Velitchkova, M. (2021). The role of antioxidant defense in freezing tolerance of resurrection plant Haberlea rhodopensis. Physiol. Mol. Biol. Plants 27, 1119–1133. doi: 10.1007/s12298-021-00998-0
Gong, B., Li, X., Vandenlangenberg, K. M., Wen, D., Sun, S., Wei, M., et al. (2014). Overexpression of S-adenosyl-L-methionine synthetase increased tomato tolerance to alkali stress through polyamine metabolism. Plant Biotechnol. J. 12, 694–708. doi: 10.1111/pbi.12173
Gross, B. L. and Zhao, Z. (2014). Archaeological and genetic insights into the origins of domesticated rice. Proc. Natl. Acad. Sci. United States America 111, 6190–6197. doi: 10.1073/pnas.1308942110
Hong, Y., Wang, Z., Shi, H., Yao, J., Liu, X., Wang, F., et al. (2020). Reciprocal regulation between nicotinamide adenine dinucleotide metabolism and abscisic acid and stress response pathways in Arabidopsis. PloS Genet. 16, e1008892. doi: 10.1371/journal.pgen.1008892
Jan, R., Khan, M. A., Asaf, S., Lubna, Waqas, M., Park, J. R., et al. (2022). Drought and UV radiation stress tolerance in rice is improved by overaccumulation of non-enzymatic antioxidant flavonoids. Antioxidants 11, 917. doi: 10.3390/antiox11050917
Jan, R., Kim, N., Lee, S. H., Khan, M. A., Asaf, S., Lubn, et al. (2021). Enhanced flavonoid accumulation reduces combined salt and heat stress through regulation of transcriptional and hormonal mechanisms. Front. Plant Sci. 12, 796956. doi: 10.3389/fpls.2021.796956
Kandil, E. E., El-Banna, A., Tabl, D. M. M., Mackled, M. I., Ghareeb, R. Y., Al-Huqail, A. A., et al. (2022). Zinc nutrition responses to agronomic and yield traits, kernel quality, and pollen viability in rice (Oryza sativa L.). Front. Plant Sci. 13, 791066. doi: 10.3389/fpls.2022.791066
Kobayashi, Y., Fukuzawa, N., Hyodo, A., Kim, H., Mashiyama, S., Ogihara, T., et al. (2020). Role of salicylic acid glucosyltransferase in balancing growth and defence for optimum plant fitness. Mol. Plant Pathol. 21, 429–442. doi: 10.1111/mpp.12906
Kou, S., Chen, L., Tu, W., Scossa, F., Wang, Y., Liu, J., et al. (2018). The arginine decarboxylase gene ADC1, associated to the putrescine pathway, plays an important role in potato cold-acclimated freezing tolerance as revealed by transcriptome and metabolome analyses. Plant J. 96, 1283–1298. doi: 10.1111/tpj.14126
Kováčik, J. and Klejdus, B. (2014). Induction of phenolic metabolites and physiological changes in chamomile plants in relation to nitrogen nutrition. Food Chem. 142, 334–341. doi: 10.1016/j.foodchem.2013.07.074
Lai, X., Duan, W., Zhang, W., Peng, Z., Wang, X., Wang, H., et al. (2024). Integrative analysis of microbiome and metabolome revealed the effect of microbial inoculant on microbial community diversity and function in rhizospheric soil under tobacco monoculture. Microbiol. Spectr. 12, 04046–04023. doi: 10.1128/spectrum.04046-23
Li, C., Chen, Y., Huang, L., Zhang, Y., Cao, N., Guo, X., et al. (2023). Potential toxicity and dietary risk of tricyclazole to Chinese mitten crab (Eriocheir sinensis) in the rice-crab co-culture model. Environ. pollut. 316, 26. doi: 10.1016/j.envpol.2022.120514
Li, J., Hou, Z., Chen, G. H., Li, F., Zhou, Y., Xue, X. Y., et al. (2014). Synthesis, antibacterial activities, and theoretical studies of dicoumarols. Org. Biomol. Chem. 12, 5528–5535. doi: 10.1039/C4OB00772G
Liu, Y., Guo, L., Liu, L., Xu, L., Kuang, H., Xu, X., et al. (2024c). A paper-based lateral flow immunochromatographic sensor for the detection of tricyclazole in rice. Food Chem. 459, 14. doi: 10.1016/j.foodchem.2024.140434
Liu, N., Huang, J., Liu, X., Wu, J., and Huang, M. (2024a). Pesticide-induced metabolic disruptions in crops: A global perspective at the molecular level. Sci. Total Environ. 957, 2. doi: 10.1016/j.scitotenv.2024.177665
Liu, X., Xie, Z., Xin, J., Yuan, S., Liu, S., Sun, Y., et al. (2024b). OsbZIP18 is a positive regulator of phenylpropanoid and flavonoid biosynthesis under UV-B radiation in rice. Plants 13, 498. doi: 10.3390/plants13040498
Liu, N. and Zhu, L. (2020). Metabolomic and transcriptomic investigation of metabolic perturbations in oryza sativa L. Triggered by three pesticides. Environ. Sci. Technol. 54, 6115–6124. doi: 10.1021/acs.est.0c00425
Lu, S. C. (2000). S-adenosylmethionine. Int. J. Biochem. Cell Biol. 32, 391–395. doi: 10.1016/S1357-2725(99)00139-9
Lu, C., Yue, Y., Liu, B., Jiang, Y., Li, Z., Latif, M. Z., et al. (2025). A multi-omics approach reveals the effects of bio- and chemical- pesticides on rice yield and quality under disease stress. Rice 18, 025–00822. doi: 10.1186/s12284-025-00822-3
Lv, Z., Jiang, S., Kong, S., Zhang, X., Yue, J., Zhao, W., et al. (2024). Advances in single-cell transcriptome sequencing and spatial transcriptome sequencing in plants. Plants 13, 1679. doi: 10.3390/plants13121679
Ma, C., Wang, Y., Gu, D., Nan, J., Chen, S., and Li, H. (2017). Overexpression of S-adenosyl-l-methionine synthetase 2 from sugar beet M14 increased arabidopsis tolerance to salt and oxidative stress. Int. J. Mol. Sci. 18, 847. doi: 10.3390/ijms18040847
Miao, Y., Li, H., Pan, J., Zhou, B., He, T., Wu, Y., et al. (2024). Salicylic acid modulates secondary metabolism and enhanced colchicine accumulation in long yellow daylily (Hemerocallis citrina). AoB Plants 16, 29. doi: 10.1093/aobpla/plae029
Mishra, S., Sharma, A., and Srivastava, A. K. (2024). Ascorbic acid: a metabolite switch for designing stress-smart crops. Crit. Rev. Biotechnol. 44, 1350–1366. doi: 10.1080/07388551.2023.2286428
Miyamoto, K. (2025). Regulatory mechanisms and evolutionary insights of phytoalexin biosynthesis in rice. Biosci. Biotechnol. Biochem. 7, 166. doi: 10.1093/bbb/zbaf166
Mortimer, J. C., Yu, X., Albrecht, S., Sicilia, F., Huichalaf, M., Ampuero, D., et al. (2013). Abnormal glycosphingolipid mannosylation triggers salicylic acid-mediated responses in Arabidopsis. Plant Cell. 25, 1881–1894. doi: 10.1105/tpc.113.111500
Napieraj, N., Janicka, M., and Reda, M. (2023). Interactions of polyamines and phytohormones in plant response to abiotic stress. Plants (Basel) 12, 1159. doi: 10.3390/plants12051159
Nazish, T., Arshad, M., Jan, S. U., Javaid, A., Khan, M. H., Naeem, M. A., et al. (2022). Transporters and transcription factors gene families involved in improving nitrogen use efficiency (NUE) and assimilation in rice (Oryza sativa L.). Transgenic Res. 31, 23–42. doi: 10.1007/s11248-021-00284-5
Noctor, G., Mhamdi, A., Chaouch, S., Han, Y., Neukermans, J., Marquez-Garcia, B., et al. (2012). Glutathione in plants: an integrated overview. Plant Cell Environ. 35, 454–484. doi: 10.1111/j.1365-3040.2011.02400.x
Park, H. L., Kim, T. L., Bhoo, S. H., Lee, T. H., Lee, S. W., and Cho, M. H. (2018). Biochemical characterization of the rice cinnamyl alcohol dehydrogenase gene family. Molecules 23, 2659. doi: 10.3390/molecules23102659
Peleg, Z., Reguera, M., Tumimbang, E., Walia, H., and Blumwald, E. (2011). Cytokinin-mediated source/sink modifications improve drought tolerance and increase grain yield in rice under water-stress. Plant Biotechnol. J. 9, 747–758. doi: 10.1111/j.1467-7652.2010.00584.x
Qi, J., Sun, G., Wang, L., Zhao, C., Hettenhausen, C., Schuman, M. C., et al. (2016). Oral secretions from Mythimna separata insects specifically induce defence responses in maize as revealed by high-dimensional biological data. Plant Cell Environ. 39, 1749–1766. doi: 10.1111/pce.12735
Qiao, Y., Zhang, N., Liu, J., and Yang, H. (2023). Interpretation of ametryn biodegradation in rice based on joint analyses of transcriptome, metabolome and chemo-characterization. J. Hazard Mater. 445, 1. doi: 10.1016/j.jhazmat.2022.130526
Rehman, S., Ikram, M., and Subhan, F. (2017). Synthesis of new dicoumarol based zinc compounds and their invitro antimicrobial studies. Rev. la Sociedad Química Mexico 59, 137–142. doi: 10.29356/jmcs.v59i2.27
Robison, F. M., Turner, M. F., Jahn, C. E., Schwartz, H. F., Prenni, J. E., Brick, M. A., et al. (2018). Common bean varieties demonstrate differential physiological and metabolic responses to the pathogenic fungus Sclerotinia sclerotiorum. Plant Cell Environ. 41, 2141–2154. doi: 10.1111/pce.13176
Shim, Y., Kim, B., Choi, Y., Cho, S. H., Kim, Y., Kim, S. H., et al. (2025). Rice OsDof12 enhances tolerance to drought stress by activating the phenylpropanoid pathway. Plant J. 121, 30. doi: 10.1111/tpj.17175
Shukla, D. (2025). Ethylene signaling in plant development and stress adaptation. Methods Mol. Biol. 2945, 205–216. doi: 10.1007/978-1-0716-4650-2
Syauqi, J., Chen, R. K., Cheng, A. H., Wu, Y. F., Chung, C. L., Lin, C. C., et al. (2022). Surveillance of rice blast resistance effectiveness and emerging virulent isolates in Taiwan. Plant Dis. 106, 3187–3197. doi: 10.1094/PDIS-12-21-2806-RE
Tian, P., Feng, Y. X., and Li, Y. H. (2024). SOS! Hydrogen sulfide enhances the flavonoid early warning system in rice plants to cope with thiocyanate pollution. Toxics 12, 591. doi: 10.3390/toxics12080591
Wang, G., Kong, J., Cui, D., Zhao, H., Niu, Y., Xu, M., et al. (2019). Resistance against Ralstonia solanacearum in tomato depends on the methionine cycle and the γ-aminobutyric acid metabolic pathway. Plant J. 97, 1032–1047. doi: 10.1111/tpj.14175
Wang, J., Su, C., Cui, Z., Huang, L., Gu, S., Jiang, S., et al. (2023). Transcriptomics and metabolomics reveal tolerance new mechanism of rice roots to Al stress. Front. Genet. 13, 1063984. doi: 10.3389/fgene.2022.1063984
Wu, F., Chen, Z., Xu, X., Xue, X., Zhang, Y., and Sui, N. (2024a). Halotolerant Bacillus sp. strain RA coordinates myo-inositol metabolism to confer salt tolerance to tomato. J. Integr. Plant Biol. 66, 1871–1885. doi: 10.1111/jipb.13733
Wu, P., Li, B., Liu, Y., Bian, Z., Xiong, J., Wang, Y., et al. (2024b). Multiple physiological and biochemical functions of ascorbic acid in plant growth, development, and abiotic stress response. Int. J. Mol. Sci. 25, 1832. doi: 10.3390/ijms25031832
Wu, L., Meng, F., Su, X., Chen, N., Peng, D., and Xing, S. (2023). Transcriptomic responses to cold stress in Dendrobium huoshanense C.Z. Tang et S.J. Cheng. Physiol. Mol. Biol. Plants 29, 1633–1646. doi: 10.1007/s12298-023-01385-7
Yang, L., Wang, X., Zhao, F., Zhang, X., Li, W., Huang, J., et al. (2023b). Roles of S-adenosylmethionine and its derivatives in salt tolerance of cotton. Int. J. Mol. Sci. 24, 9517. doi: 10.3390/ijms24119517
Yang, D. Y., Zhuang, K. Y., and Ma, N. N. (2023a). Overexpression of SlGGP-LIKE gene enhanced the resistance of tomato to salt stress. Protoplasma 260, 625–635. doi: 10.1007/s00709-022-01800-y
Zavaliev, R. and Dong, X. (2023). NPR1, a key immune regulator for plant survival under biotic and abiotic stresses. Mol. Cell. 84, 131–141. doi: 10.1016/j.molcel.2023.11.018
Zavaliev, R., Mohan, R., Chen, T., and Dong, X. (2020). Formation of NPR1 condensates promotes cell survival during the plant immune response. Cell 182, 1093–1108. doi: 10.1016/j.cell.2020.07.016
Zhu, K., Bao, H., Han, S., Zhao, Y., Wang, M., Zhao, R., et al. (2023). Analysis of key differential metabolites in sunflower after downy mildew infection. Romanian Agric. Res. 40, 65. doi: 10.59665/rar4007
Keywords: metabolome, rice, stress, transcriptome, tricyclazole
Citation: Huang W, Shi D, Cheng A, Chen G, Liu F, Zhang X, Dong J, Lan J, Ren H, Guo W and Liu B (2026) The responses of rice plant to tricyclazole at the transcriptome and metabolome levels. Front. Plant Sci. 17:1723722. doi: 10.3389/fpls.2026.1723722
Received: 13 October 2025; Accepted: 19 January 2026; Revised: 13 January 2026;
Published: 03 February 2026.
Edited by:
Can Wang, Southwest Jiaotong University, ChinaReviewed by:
Shumei Fang, Heilongjiang Bayi Agricultural University, ChinaAkshay Patil, Anand Agricultural University, India
Sifu Li, Hunan Academy of Agricultural Sciences (CAAS), China
Copyright © 2026 Huang, Shi, Cheng, Chen, Liu, Zhang, Dong, Lan, Ren, Guo and Liu. 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: Wengong Huang, aHVhbmd3ZW5nb25nMTczNkAxNjMuY29t; Baohai Liu, c2hzbGJoQDEyNi5jb20=
Dongmei Shi