Combined metabolome and transcriptome analyses reveal that growing under Red shade affects secondary metabolite content in Huangjinya green tea

Shading treatments impact the tea (Camellia sinensis L.) quality. The sunlight sensitive varieties can be grown under shading nets for better growth and secondary metabolite content. Here, we studied the responses of a sunlight sensitive green tea variety "Huangjinya" by growing under colored shading nets (red, yellow, blue, and black (75% and 95%) shading rates) to find out the most suitable color of the shading net. Red shading was the most promising treatment as it positively affected the weight and length of 100 one-bud-three leaves and reduced the degree and rate of new shoots burn compared to control (natural sunlight). We then explored the comparative metabolomic changes in response to red shading by using UPLC-ESI-MS/MS system. The amino acids and derivatives, flavonoids, and alkaloids were downaccumulated whereas lipids, organic acids, and lignans were upaccumulated in Red shade grown tea samples. The red shading nets caused a decreased catechin, epicatechin, dopamine, and L-tyramine contents but increased caffeine content. We then employed transcriptome sequencing to find key changes in expressions of related genes and pathways. Notably, key genes associated with the phenylpropanoid and flavonoid biosynthesis pathways exhibited complex regulation. These expression changes suggested a potential trend of polymerization or condensation of simple molecules like catechin or pelargonidin into larger molecules like glucoside or proanthocyanidins. Here, Red shading net triggered higher expression of genes enriched in lipid biosynthesis and jasmonic acid biosynthesis, suggesting an interplay of fatty acids and JA in improving tea performance. These findings contribute to the metabolic responses of Huangjinya tea to red shading nets which might have implications for flavor and health benefits. Our data provide a foundation for further exploration and optimization of cultivation practices for this unique tea variety.

Morphological and molecular response mechanisms of the root system of different Hemarthria compressa species to submergence stress

Introduction

The severity of flood disasters is increasing due to climate change, resulting in a significant reduction in the yield and quality of forage crops worldwide. This poses a serious threat to the development of agriculture and livestock. Hemarthria compressa is an important high-quality forage grass in southern China. In recent years, frequent flooding has caused varying degrees of impacts on H. compressa and their ecological environment.

Methods

In this study, we evaluated differences in flooding tolerance between the root systems of the experimental materials GY (Guang Yi, flood-tolerant) and N1291 (N201801291, flood-sensitive). We measured their morphological indexes after 7 d, 14 d, and 21 d of submergence stress and sequenced their transcriptomes at 8 h and 24 h, with 0 h as the control.

Results

During submergence stress, the number of adventitious roots and root length of both GY and N1291 tended to increase, but the overall growth of GY was significantly higher than that of N1291. RNA-seq analysis revealed that 6046 and 7493 DEGs were identified in GY-8h and GY-24h, respectively, and 9198 and 4236 DEGs in N1291-8h and N1291-24h, respectively, compared with the control. The GO and KEGG enrichment analysis results indicated the GO terms mainly enriched among the DEGs were oxidation-reduction process, obsolete peroxidase reaction, and other antioxidant-related terms. The KEGG pathways that were most significantly enriched were phenylpropanoid biosynthesis, plant hormone signal transduction etc. The genes of transcription factor families, such as C2H2, bHLH and bZIP, were highly expressed in the H. compressa after submergence, which might be closely related to the submergence adaptive response mechanisms of H. compressa.

Discussion

This study provides basic data for analyzing the molecular and morphological mechanisms of H. compressa in response to submergence stress, and also provides theoretical support for the subsequent improvement of submergence tolerance traits of H. compressa.

Enhancing transcriptome analysis in medicinal plants: multiple unigene sets in Astragalus membranaceus

Astragalus membranaceus is a medicinal plant mainly used in East Asia and contains abundant secondary metabolites. Despite the importance of this plant, the available genomic and genetic information is still limited. De novo transcriptome construction is recognized as an essential method for transcriptome research when reference genome information is incomplete. In this study, we constructed three individual transcriptome sets (unigene sets) for detailed analysis of the phenylpropanoid biosynthesis pathway, a major metabolite of A. membranaceus. Set-1 was a circular consensus sequence (CCS) generated using PacBio sequencing (PacBio-seq). Set-2 consisted of hybridized assembled unigenes with Illumina sequencing (Illumina-seq) reads and PacBio CCS using rnaSPAdes. Set-3 unigenes were assembled from Illumina-seq reads using the Trinity software. Construction of multiple unigene sets provides several advantages for transcriptome analysis. First, it provides an appropriate expression filtering threshold for assembly-based unigenes: a threshold transcripts per million (TPM) ≥ 5 removed more than 88% of assembly-based unigenes, which were mostly short and low-expressing unigenes. Second, assembly-based unigenes compensated for the incomplete length of PacBio CCSs: the ends of the 5`/3` untranslated regions of phenylpropanoid-related unigenes derived from set-1 were incomplete, which suggests that PacBio CCSs are unlikely to be full-length transcripts. Third, more isoform unigenes could be obtained from multiple unigene sets; isoform unigenes missing in Set-1 were detected in set-2 and set-3. Finally, gene ontology and Kyoto Encyclopedia of Genes and Genomes analyses showed that phenylpropanoid biosynthesis and carbohydrate metabolism were highly activated in A. membranaceus roots. Various sequencing technologies and assemblers have been developed for de novo transcriptome analysis. However, no technique is perfect for de novo transcriptome analysis, suggesting the need to construct multiple unigene sets. This method enables efficient transcript filtering and detection of longer and more diverse transcripts.

Transcriptome Profiles Reveals ScDREB10 from Syntrichia caninervis Regulated Phenylpropanoid Biosynthesis and Starch/Sucrose Metabolism to Enhance Plant Stress Tolerance

Desiccation is a kind of extreme form of drought stress and desiccation tolerance (DT) is an ancient trait of plants that allows them to survive tissue water potentials reaching -100 MPa or lower. ScDREB10 is a DREB A-5 transcription factor gene from a DT moss named Syntrichia caninervis, which has strong comprehensive tolerance to osmotic and salt stresses. This study delves further into the molecular mechanism of ScDREB10 stress tolerance based on the transcriptome data of the overexpression of ScDREB10 in Arabidopsis under control, osmotic and salt treatments. The transcriptional analysis of weight gene co-expression network analysis (WGCNA) showed that "phenylpropanoid biosynthesis" and "starch and sucrose metabolism" were key pathways in the network of cyan and yellow modules. Meanwhile, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of differentially expressed genes (DEGs) also showed that "phenylpropanoid biosynthesis" and "starch and sucrose metabolism" pathways demonstrate the highest enrichment in response to osmotic and salt stress, respectively. Quantitative real-time PCR (qRT-PCR) results confirmed that most genes related to phenylpropanoid biosynthesis" and "starch and sucrose metabolism" pathways in overexpressing ScDREB10 Arabidopsis were up-regulated in response to osmotic and salt stresses, respectively. In line with the results, the corresponding lignin, sucrose, and trehalose contents and sucrose phosphate synthase activities were also increased in overexpressing ScDREB10 Arabidopsis under osmotic and salt stress treatments. Additionally, cis-acting promoter element analyses and yeast one-hybrid experiments showed that ScDREB10 was not only able to bind with classical cis-elements, such as DRE and TATCCC (MYBST1), but also bind with unknown element CGTCCA. All of these findings suggest that ScDREB10 may regulate plant stress tolerance by effecting phenylpropanoid biosynthesis, and starch and sucrose metabolism pathways. This research provides insights into the molecular mechanisms underpinning ScDREB10-mediated stress tolerance and contributes to deeply understanding the A-5 DREB regulatory mechanism.

Integrated Transcriptomics and Metabolomics Analysis Promotes the Understanding of Adventitious Root Formation in Eucommia ulmoides Oliver

As a primary approach to nutrient propagation for many woody plants, cutting roots is essential for the breeding and production of Eucommia ulmoides Oliver. In this study, hormone level, transcriptomics, and metabolomics analyses were performed on two E. ulmoides varieties with different adventitious root (AR) formation abilities. The higher JA level on the 0th day and the lower JA level on the 18th day promoted superior AR development. Several hub genes executed crucial roles in the crosstalk regulation of JA and other hormones, including F-box protein (EU012075), SAUR-like protein (EU0125382), LOB protein (EU0124232), AP2/ERF transcription factor (EU0128499), and CYP450 protein (EU0127354). Differentially expressed genes (DEGs) and metabolites of AR formation were enriched in phenylpropanoid biosynthesis, flavonoid biosynthesis, and isoflavonoid biosynthesis pathways. The up-regulated expression of PAL, CCR, CAD, DFR, and HIDH genes on the 18th day could contribute to AR formation. The 130 cis-acting lncRNAs had potential regulatory functions on hub genes in the module that significantly correlated with JA and DEGs in three metabolism pathways. These revealed key molecules, and vital pathways provided more comprehensive insight for the AR formation mechanism of E. ulmoides and other plants.

Dynamic Changes of Phenolic Composition, Antioxidant Capacity, and Gene Expression in 'Snow White' Loquat (Eriobotrya japonica Lindl.) Fruit throughout Development and Ripening

The newly released 'Snow White' (SW), a white-fleshed loquat (Eriobotrya japonica Lindl.) cultivar, holds promise for commercial production. However, the specifics of the phenolic composition in white-fleshed loquats, along with the antioxidant substances and their regulatory mechanisms, are not yet fully understood. In this study, we examined the dynamic changes in the phenolic compounds, enzyme activities, antioxidant capacity, and gene expression patterns of SW during the key stages of fruit development and ripening. A total of 18 phenolic compounds were identified in SW, with chlorogenic acid, neochlorogenic acid, and coniferyl alcohol being the most predominant. SW demonstrated a stronger antioxidant capacity in the early stages of development, largely due to total phenolics and flavonoids. Neochlorogenic acid may be the most significant antioxidant contributor in loquat. A decline in enzyme activities corresponded with fruit softening. Different genes within a multigene family played distinct roles in the synthesis of phenolics. C4H1, 4CL2, 4CL9, HCT, CCoAOMT5, F5H, COMT1, CAD6, and POD42 were implicated in the regulation of neochlorogenic acid synthesis and accumulation. Consequently, these findings enhance our understanding of phenolic metabolism and offer fresh perspectives on the development of germplasm resources for white-fleshed loquats.

Integrated comparative transcriptome and physiological analysis reveals the metabolic responses underlying genotype variations in NH4+ tolerance

Several mechanisms have been proposed to explain NH4 + toxicity. However, the core information about the biochemical regulation of plants in response to NH4 + toxicity is still lacking. In this study, the tissue NH4 + concentration is an important factor contributing to variations in plant growth even under nitrate nutrition and NH4 + tolerance under ammonium nutrition. Furthermore, NH4 + led to the reprogramming of the transcriptional profile, as genes related to trehalose-6-phosphate and zeatin biosynthesis were downregulated, whereas genes related to nitrogen metabolism, camalexin, stilbenoid and phenylpropanoid biosynthesis were upregulated. Further analysis revealed that a large number of genes, which enriched in phenylpropanoid and stilbenoid biosynthesis, were uniquely upregulated in the NH4 +- tolerant ecotype Or-1. These results suggested that the NH4 +-tolerant ecotype showed a more intense response to NH4 + by activating defense processes and pathways. Importantly, the tolerant ecotype had a higher 15NH4 + uptake and nitrogen utilization efficiency, but lower NH4 +, indicating the tolerant ecotype maintained a low NH4 + level, mainly by promoting NH4 + assimilation rather than inhibiting NH4 + uptake. The carbon and nitrogen metabolism analysis revealed that the tolerant ecotype had a stronger carbon skeleton production capacity with higher levels of hexokinase, pyruvate kinase, and glutamate dehydrogenase activity to assimilate free NH4 +, Taken together, the results revealed the core mechanisms utilized by plants in response to NH4 +, which are consequently of ecological and agricultural importance.

Physiological and transcriptomic analyses of response of walnuts (Juglans regia) to Pantoea agglomerans infection

Introduction

Walnut blight is a serious bacterial disease that affects the yield and quality of walnuts. Pantoea agglomerans is one of the main causative agents of walnut blight. However, there have been few studies on the response of walnuts to P. agglomerans infection.

Methods

In this study, the soluble sugar, photosynthesis, antioxidant enzyme activities, and secondary metabolites were measured, and the transcriptomic analysis was performed to determine the response of walnut tissue cultures to P. agglomerans infection.

Results

After pathogen inoculation, the soluble sugar content decreased, and photosynthesis was inhibited. Antioxidant enzyme (superoxide dismutase and peroxidase) activities and secondary metabolites (phenol and flavonoid) contents increased, especially in the early stages of inoculation. Transcriptomic analysis revealed that the phenylpropanoid biosynthesis pathway is induced after infection, and pathogen infection promotes ABA and ethylene signal transduction and inhibits auxin signaling. In addition, SA and JA-related gene expression was altered after inoculation with P. agglomerans, and the FLS- and calcium-mediated disease resistance signaling pathways were activated. Furthermore, our results suggested an involvement of the R-protein RPM-mediated disease resistance pathway in the response of walnuts to bacterial infections.

Discussion

Our findings indicated that phenylpropanoid biosynthesis, hormone signal transduction, and plant-pathogen interaction have key roles in pathogenic inoculation, which provide insights into the molecular mechanisms in the response of walnuts to P. agglomerans infection.

Effects of Altering Three Steps of Monolignol Biosynthesis on Sorghum Responses to Stalk Pathogens and Water Deficit

The drought-resilient crop sorghum (Sorghum bicolor [L.] Moench) is grown worldwide for multiple uses, including forage or potential lignocellulosic bioenergy feedstock. A major impediment to biomass yield and quality are the pathogens Fusarium thapsinum and Macrophomina phaseolina, which cause Fusarium stalk rot and charcoal rot, respectively. These fungi are more virulent with abiotic stresses such as drought. Monolignol biosynthesis plays a critical role in plant defense. The genes Brown midrib (Bmr)6, Bmr12, and Bmr2 encode the monolignol biosynthesis enzymes cinnamyl alcohol dehydrogenase, caffeic acid O-methyltransferase, and 4-coumarate:CoA ligase, respectively. Plant stalks from lines overexpressing these genes and containing bmr mutations were screened for pathogen responses with controlled adequate or deficit watering. Additionally, near-isogenic bmr12 and wild-type lines in five backgrounds were screened for response to F. thapsinum with adequate and deficit watering. All mutant and overexpression lines were no more susceptible than corresponding wild-type under both watering conditions. The bmr2 and bmr12 lines, near-isogenic to wild-type, had significantly shorter mean lesion lengths (were more resistant) than RTx430 wild-type when inoculated with F. thapsinum under water deficit. Additionally, bmr2 plants grown under water deficit had significantly smaller mean lesions when inoculated with M. phaseolina than under adequate-water conditions. When well-watered, bmr12 in cultivar Wheatland and one of two Bmr2 overexpression lines in RTx430 had shorter mean lesion lengths than corresponding wild-type lines. This research demonstrates that modifying monolignol biosynthesis for increased usability may not impair plant defenses but can even enhance resistance to stalk pathogens under drought conditions.

How Antioxidants, Osmoregulation, Genes and Metabolites Regulate the Late Seeding Tolerance of Rapeseeds (Brassica napus L.) during Wintering

Rapeseed seeding dates are largely delayed under the rice-rape rotation system, but how rapeseeds adapt to the delayed environment remains unclear. Here, five seeding dates (20 October, 30 October, 10 November, 20 November and 30 November, T1 to T5) were set and the dynamic differences between two late-seeding-tolerant (LST) and two late-seeding-sensitive (LSS) rapeseed cultivars were investigated in a field experiment. The growth was significantly repressed and the foldchange (LST/LSS) of yield increased from 1.50-T1 to 2.64-T5 with the delay in seeding. Both LST cultivars showed higher plant coverage than the LSS cultivars according to visible/hyperspectral imaging and the vegetation index acquired from an unmanned aerial vehicle. Fluorescence imaging, DAB and NBT staining showed that the LSS cultivars suffered more stress damage than the LST cultivars. Antioxidant enzymes (SOD, POD, CAT, APX) and osmoregulation substances (proline, soluble sugar, soluble protein) were decreased with the delay in seeding, while the LST cultivar levels were higher than those of the LSS cultivars. A comparative analysis of transcriptomes and metabolomes showed that 55 pathways involving 123 differentially expressed genes (DEGs) and 107 differentially accumulated metabolites (DAMs) participated in late seeding tolerance regulation, while 39 pathways involving 60 DEGs and 68 DAMs were related to sensitivity. Levanbiose, α-isopropylmalate, s-ribosyl-L-homocysteine, lauroyl-CoA and argino-succinate were differentially accumulated in both cultivars, while genes including isocitrate dehydrogenase, pyruvate kinase, phosphoenolpyruvate carboxykinase and newgene_7532 were also largely regulated. This study revealed the dynamic regulation mechanisms of rapeseeds on late seeding conditions, which showed considerable potential for the genetic improvement of rapeseed.

Phenolic, Polysaccharides Composition, and Texture Properties during Ripening and Storage Time of New Table Grape Cultivars in Chile

The aim of this study is to determine the phenolic and polysaccharidic composition, texture properties, and gene expression of new seedless table grape cultivars Timco™ and Krissy™ and compare them to the traditional table grape variety Crimson Seedless (Vitis vinifera L.), during ripening and in commercial postharvest conditions. According to the results, phenolic compounds were present in very different proportions. The total anthocyanins responsible for skin color increased during maturation and the majority anthocyanin in the three cultivars was peonidin-3-glucoside, followed by malvidin-3-glucoside. The phenolic compounds presented a different behavior (decreasing or increasing) during postharvest. The total skin soluble polysaccharides decreased during ripening and postharvest in Crimson Seedless and Krissy™ and remained constant from technological maturity to postharvest storage in Timco™. In all cultivars, the majority soluble polysaccharide fraction was that with a molecular mass between 500 and 35 KDa. The skin mechanical properties of table grapes were good parameters for differentiating varieties, with better results for the new cultivars, compared to the traditional Crimson Seedless, especially in postharvest. Genes involved in the flavonoid pathway and cell wall metabolism in skins exhibited an increase in expression from veraison to remaining constant at the end of the berry ripening.

Transcriptome and Metabolome Analysis of Selenium Treated Alfalfa Reveals Influence on Phenylpropanoid Biosynthesis to Enhance Growth

Selenium (Se) plays an important role in the growth of plants. Alfalfa (Medicago sativa L.) is a perennial legume forage crop with high nutritional value and Se-rich functions. Many studies have shown that selenium can promote alfalfa growth, but few have explored the molecular biology mechanisms behind this effect. In this study, alfalfa was divided into two groups. One group was sprayed with sodium selenite (Na2SeO3) and the other group was sprayed with distilled water as a control. This study determined the growth, reproductive traits, physiological changes, transcriptome and metabolome of both groups of alfalfa. We found that foliar spraying of 100 mg/L Na2SeO3 could significantly increase the growth rate, dry weight, total Se content, amount of pollen per flower, pollen viability, pod spirals, and seed number per pod of alfalfa plants. The level of chlorophyll, soluble protein, proline, and glutathione also increased dramatically in Na2SeO3-sprayed alfalfa seedlings. After transcriptome and metabolome analysis, a total of 614 differentially expressed genes (DEGs) and 1500 differentially expressed metabolites (DEMs), including 26 secondary differentially metabolites were identified. The DEGs were mainly enriched in MAPK signaling pathway, phenylpropanoid biosynthesis, isoflavonoid biosynthesis, cutin, suberine, and wax biosynthesis, and glycerolipid metabolism. The DEMs were mainly enriched in flavone and flavonol biosynthesis, carbon metabolism, glyoxylate and dicarboxylate metabolism, nitrogen metabolism, and phenylpropanoid biosynthesis. Integrative analysis of transcriptome and metabolome showed that the foliar spraying of Na2SeO3 mainly affects phenylpropanoid biosynthesis to promote alfalfa growth.

Potential Toxic Mechanisms of Neonicotinoid Insecticides in Rice: Inhibiting Auxin-Mediated Signal Transduction

Inappropriate application of pesticides not only causes sub-lethal effects on ecosystem service providers but also reduces crop yield and quality. As a xenobiotic signal molecule, pesticides may interact with signal transduction receptors in crops, resulting in oxidative damage and even metabolic perturbations. We discovered that three neonicotinoid insecticides (NIs), namely, imidacloprid, thiamethoxam, and clothianidin, at 0.06-0.12 kg ai/ha significantly inhibited the auxin signal pathway in rice leaves, thereby reducing the intracellular auxin (IAA) content. Molecular simulation further confirmed that NIs occupied the binding site where auxin transporter-like proteins 1 (LAX11) and 2 (LAX12), in which Thr253 and Asn66 of LAX11, as well as Thr244 and Asn57 of LAX12, were bound to the nitroguanidine of NIs via H-bonds. Meanwhile, Asn66 of LAX11 and Asn57 of LAX12 interacted with nitroguanidine via aromatic H-bonds. Moreover, phenylpropanoid biosynthesis was significantly disturbed because of the inhibited auxin signal pathway. Notably, peroxidase-coding genes were downregulated with a maximum value greater than 10-fold, resulting in decreased antioxidant metabolites flavone (37.82%) and lignin content (20.15%). Ultimately, rice biomass was reduced by up to 25.41% due to the decline in IAA content and antioxidant capacity. This study deeply explored the molecular mechanism of metabolic perturbations in crops stressed by pesticides, thus providing a scientific basis for pesticide environmental risk assessment and agricultural product safety.

Benzothiadiazole Affects Grape Polyphenol Metabolism and Wine Quality in Two Greek Cultivars: Effects during Ripening Period over Two Years

Grape berries are one of the most important sources of phenolic compounds, either consumed fresh or as wine. A pioneer practice aiming to enrich grape phenolic content has been developed based on the application of biostimulants such as agrochemicals initially designed to induce resistance against plant pathogens. A field experiment was conducted in two growing seasons (2019-2020) to investigate the effect of benzothiadiazole on polyphenol biosynthesis during grape ripening in Mouhtaro (red-colored) and Savvatiano (white-colored) varieties. Grapevines were treated at the stage of veraison with 0.3 mM and 0.6 mM benzothiadiazole. The phenolic content of grapes, as well as the expression level of genes involved in the phenylpropanoid pathway were evaluated and showed an induction of genes specifically engaged in anthocyanins and stilbenoids biosynthesis. Experimental wines deriving from benzothiadiazole-treated grapes exhibited increased amounts of phenolic compounds in both varietal wines, as well as an enhancement in anthocyanin content of Mouhtaro wines. Taken together, benzothiadiazole can be utilized to induce the biosynthesis of secondary metabolites with oenological interest and to improve the quality characteristics of grapes produced under organic conditions.

Transcriptional profiling analysis providing insights into desiccation tolerance mechanisms of the desert moss Syntrichia caninervis

Syntrichia caninervis is a desiccation tolerant moss and is the dominant bryophyte found in biological soil crusts in the Gurbantunggut desert. In this study, we assessed the transcriptome profiles of S. caninervis gametophytes during the dehydration-rehydration (D-R) process (across 9 time points) using Illumina sequencing. In total, 22489 transcripts were identified, including 5337 novel transcripts, that mapped to the reference genome. A total of 12548 transcripts exhibited significant alterations in the D-R samples compared with the control samples. The differentially expressed transcripts (DETs) possessed several enriched Gene Ontology terms, such as "water stress response", "oxidation-reduction process", "membrane metabolism", "photosynthesis", and "transcription factor activity". Moreover, during early dehydration stress, the DETs were significantly enriched in stress-related pathways from the Kyoto Encyclopedia of Genes and Genomes, such as "phenylpropanoid biosynthesis", "alpha-linolenic acid metabolism", and "fructose and mannose metabolism". Photosynthesis-related transcripts (e.g., ScPsa H, ScRubisco, and ScLhcb1) were inhibited during the dehydration treatment and significantly accumulated during the late rehydration period. Most transcripts from the late embryogenesis abundant proteins (LEA) and early light-inducible protein (ELIP) families strongly accumulated at the late dehydration stage. These pathways were positively correlated with the content changes of absolute water content and Fv/Fm values, alongside peroxidase and superoxide dismutase activities. Seven transcription factor families, including AP2-ERF, bHLH, G2-like, MYB, NAC, WRKY, and bZIP, were enriched in DETs during D-R treatment. This study is the first transcriptome analysis using the S. caninervis genome for gene annotation and multigroup D-R treatment points. Our results demonstrated the detailed dynamic changes in the transcriptome of S. caninervis during the D-R process. These results also improve understanding of desiccation tolerant plants' adaptations to desiccation stress at the transcription level and provide promising gene resources for transgenic crop breeding.

Identification of drought-responsive phenolic compounds and their biosynthetic regulation under drought stress in Ligularia fischeri

Ligularia fischeri, a leafy edible plant found in damp shady regions, has been used as an herbal medicine and is also consumed as a horticultural crop. In this study, we investigated the physiological and transcriptomic changes, especially those involved in phenylpropanoid biosynthesis, induced by severe drought stress in L. fischeri plants. A distinguishing characteristic of L. fischeri is a color change from green to purple due to anthocyanin biosynthesis. We chromatographically isolated and identified two anthocyanins and two flavones upregulated by drought stress using liquid chromatography-mass spectrometry and nuclear magnetic resonance analyses in this plant for the first time. In contrast, all types of caffeoylquinic acids (CQAs) and flavonol contents were decreased under drought stress. Further, we performed RNA sequencing to examine the molecular changes in these phenolic compounds at the transcriptome level. In an overview of drought-inducible responses, we identified 2,105 hits for 516 distinct transcripts as drought-responsive genes. Moreover, differentially expressed genes (DEGs) associated with phenylpropanoid biosynthesis accounted for the greatest number of both up- and downregulated DEGs by Kyoto Encyclopedia of Genes and Genomes enrichment analysis. We identified 24 meaningful DEGs based on the regulation of phenylpropanoid biosynthetic genes. Potential drought-responsive genes included upregulated flavone synthase (LfFNS, TRINITY DN31661 c0 g1 i1) and anthocyanin 5-O-glucosyltransferase (LfA5GT1, TRINITY DN782 c0 g1 i1), which could contribute to the high levels of flavones and anthocyanins under drought stress in L. fischeri. In addition, the downregulated shikimate O-hydroxycinnamolytransferase (LfHCT, TRINITY DN31661 c0 g1 i1) and hydroxycinnamoyl-CoA quinate/shikimate transferase (LfHQT4, TRINITY DN15180 c0 g1 i1) genes led to a reduction in CQAs. Only one or two BLASTP hits for LfHCT were obtained for six different Asteraceae species. It is possible that the HCT gene plays a crucial role in CQAs biosynthesis in these species. These findings expand our knowledge of the response mechanisms to drought stress, particularly regarding the regulation of key phenylpropanoid biosynthetic genes in L. fischeri.

The ScAPD1-like gene from the desert moss Syntrichia caninervis enhances resistance to Verticillium dahliae via phenylpropanoid gene regulation

Soloist is a member of a distinct and small subfamily within the AP2/ERF transcriptional factor family that play important roles in plant biotic and abiotic stress responses. There are limited studies of Soloist genes and their functions are poorly understood. We characterized the abiotic and biotic stress tolerance function of the ScSoloist gene (designated as ScAPD1-like) from the desert moss Syntrichia caninervis. ScAPD1-like responded to multiple abiotic, biotic stresses and plant hormone treatments. ScAPD1-like protein located to the nucleus and bound to several DNA elements. Overexpression of ScAPD1-like in Arabidopsis did not alter abiotic stress resistance or inhibit Pseudomonas syringae pv. tomato (Pst) DC3000 infection. However, overexpression of ScAPD1-like significantly increased the resistance of transgenic Arabidopsis and S. caninervis to Verticillium dahliae infection, decreased reactive oxygen species accumulation and improved reactive oxygen species scavenging activity. ScAPD1-like overexpression plants altered the abundance of transcripts for lignin synthesis and promoted lignin accumulation in Arabidopsis. ScAPD1-like directly bind to RAV1, AC elements, and TATA-box in the promoters of AtPAL1 and AtC4H genes, respectively, in vitro. Chromatin immunoprecipitation-quantitative polymerase chain reaction assays demonstrated ScAPD1-like directly bound to PAL and C4H genes promoters in Arabidopsis and their homologs in S. caninervis. In S. caninervis, ScAPD1-like overexpression and RNAi directly regulated the abundance of ScPAL and ScC4H transcripts and modified the metabolites of phenylpropanoid pathway. We provide insight into the function of Soloist in plant defense mechanisms that likely occurs through activation of the phenylpropanoid biosynthesis pathway. ScAPD1-like is a promising candidate gene for breeding strategies to improve resistance to Verticillium wilt.

Integrated metabolomic and transcriptomic analysis reveals the role of phenylpropanoid biosynthesis pathway in tomato roots during salt stress

As global soil salinization continues to intensify, there is a need to enhance salt tolerance in crops. Understanding the molecular mechanisms of tomato (Solanum lycopersicum) roots' adaptation to salt stress is of great significance to enhance its salt tolerance and promote its planting in saline soils. A combined analysis of the metabolome and transcriptome of S. lycopersicum roots under different periods of salt stress according to changes in phenotypic and root physiological indices revealed that different accumulated metabolites and differentially expressed genes (DEGs) associated with phenylpropanoid biosynthesis were significantly altered. The levels of phenylpropanoids increased and showed a dynamic trend with the duration of salt stress. Ferulic acid (FA) and spermidine (Spd) levels were substantially up-regulated at the initial and mid-late stages of salt stress, respectively, and were significantly correlated with the expression of the corresponding synthetic genes. The results of canonical correlation analysis screening of highly correlated DEGs and construction of regulatory relationship networks with transcription factors (TFs) for FA and Spd, respectively, showed that the obtained target genes were regulated by most of the TFs, and TFs such as MYB, Dof, BPC, GRAS, and AP2/ERF might contribute to the regulation of FA and Spd content levels. Ultimately, FA and Spd attenuated the harm caused by salt stress in S. lycopersicum, and they may be key regulators of its salt tolerance. These findings uncover the dynamics and possible molecular mechanisms of phenylpropanoids during different salt stress periods, providing a basis for future studies and crop improvement.

RNA-Seq Analysis Demonstrates Different Strategies Employed by Tiger Nuts (Cyperus esculentus L.) in Response to Drought Stress

Drought stress, an important abiotic stress, has affected global agricultural production by limiting the yield and the quality of crops. Tiger nuts (Cyperus esculentus L.) are C4 crops in the Cyperaceae family, which have high-quality wholesome ingredients. However, data on mechanisms underlying the response of tiger nuts to drought stress are few. Here, the variety of Jisha 1 and 15% polyethylene glycol (PEG; a drought stress simulator) were used to study the mechanisms of stress response in tiger nuts. Our evaluation of the changes in physiological indicators such as electrolyte leakage (El), malondialdehyde (MDA), hydrogen peroxide (H₂O₂), superoxide anion (O₂⁻) and activities of reactive oxygen species (ROS) showed that 12 h was the most suitable time point to harvest and analyze the response to drought stress. Thereafter, we performed transcriptome (RNA-Seq) analysis in the control (CK) and stress treatment groups and showed that there was a total of 1760 differentially expressed genes (DEGs). Gene Ontology (GO) analysis showed that the DEGs were enriched in abscisic acid (ABA) terms, and pathways such as starch and sucrose metabolism (ko00500), phenylpropanoid biosynthesis (ko00940) and plant hormone signal transduction (ko04075) were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. In addition, quantitative real-time PCR (qRT-PCR) analysis of the DEGs demonstrated an upregulation of ABA and lignin content, as well as enzyme activities in enriched pathways, which validated the RNA-Seq data. These results revealed the pathways and mechanisms adopted by the tiger nuts in response to drought stress.

Halotolerant Bacillus altitudinis WR10 improves salt tolerance in wheat via a multi-level mechanism

Soil salinity is an important abiotic stress factor that seriously affects the crop growth and yield. Use of plant-derived microorganisms is a promising strategy to alleviate salt stress. In a previous study, the endophytic strain Bacillus altitudinis WR10 isolated from wheat roots showed high salt resistance. In this study, we investigated the efficacy of WR10 in improving the salt tolerance of wheat and its potential mechanisms using a hydroponic test. Under salt stress, WR10 inoculation significantly increased the lengths and dry weights of the roots and shoots, indicating that WR10 improves wheat salt tolerance at the seedling stage. WR10 inoculation significantly reduced Na⁺ accumulation and enhanced K⁺, P, and Ca²⁺ uptake in salt-stressed plants, which can be attributed to the upregulated gene expression of H⁺-ATPase as well as the P-solubilizing and biofilm-producing characteristics of WR10. At the transcriptional level, L-ascorbate peroxidase (APX), glutathione (GSH) synthetase related to GSH biosynthesis, and phenylpropanoid biosynthesis genes (CYP73A, 4CL, and CAD) were significantly upregulated, whereas those of GSH metabolism genes (glutathione S-transferase and gamma-glutamyltranspeptidase) were significantly downregulated in WR10-applied wheat roots under salt stress. These changes increased the APX activity and GSH levels and resulted in a decrease in hydrogen peroxide levels. Additionally, a decrease in proline content was observed in WR10-inoculated plants under salt stress because of WR10-induced upregulation of proline dehydrogenase gene expression. These results provide supporting evidence that WR10 improves wheat salt tolerance via more than one mechanism and open a window of opportunity for WR10 application in salinized soil.

OsGRP3 Enhances Drought Resistance by Altering Phenylpropanoid Biosynthesis Pathway in Rice (Oryza sativa L.)

As a sessile organism, rice often faces various kinds of abiotic stresses, such as drought stress. Drought stress seriously harms plant growth and damages crop yield every year. Therefore, it is urgent to elucidate the mechanisms of drought resistance in rice. In this study, we identified a glycine-rich RNA-binding protein, OsGRP3, in rice. Evolutionary analysis showed that it was closely related to OsGR-RBP4, which was involved in various abiotic stresses. The expression of OsGRP3 was shown to be induced by several abiotic stress treatments and phytohormone treatments. Then, the drought tolerance tests of transgenic plants confirmed that OsGRP3 enhanced drought resistance in rice. Meanwhile, the yeast two-hybrid assay, bimolecular luminescence complementation assay and bimolecular fluorescence complementation assay demonstrated that OsGRP3 bound with itself may affect the RNA chaperone function. Subsequently, the RNA-seq analysis, physiological experiments and histochemical staining showed that OsGRP3 influenced the phenylpropanoid biosynthesis pathway and further modulated lignin accumulation. Herein, our findings suggested that OsGRP3 enhanced drought resistance in rice by altering the phenylpropanoid biosynthesis pathway and further increasing lignin accumulation.

Transcriptomic and Non-Targeted Metabolomic Analyses Reveal the Flavonoid Biosynthesis Pathway in Auricularia cornea

Flavonoids, which are abundant in plants, are recognized for their antioxidant and anticancer roles in clinical applications. However, little is known about the molecular basis of flavonoid biosynthesis in fungi. In this study, we found that inclusion of leachate of Korshinsk peashrub (Caragana korshinskii) in the fermentation medium increased the total flavonoid content of the edible fungus Auricularia cornea by 23.6% relative to that grown in a control medium. Combined transcriptomic and non-targeted metabolomic analysis of the flavonoid biosynthesis pathway in A. cornea illustrated that there are important metabolites in the phenylpropanoid, coumarin and isoflavonoid biosynthesis pathways. In addition, we found that certain homologous genes encode phenylalanine ammonia lyase (PAL), polyphenol oxidase (PPO) and chalcone isomerase (CHI) in these biosynthesis pathways. These results, in this study, provide a new line for studying the regulation of flavonoid production in edible fungi.

Growth, Stoichiometry, and Palatability of Suaeda salsa From Different Habitats Are Demonstrated by Differentially Expressed Proteins and Their Enriched Pathways

Suaeda salsa (L.) Pall., a medicinal and edible plant, has green and red-violet ecotypes that exhibit different phenotypes, tastes, and growth characteristics. However, few studies have focused on these differences from the aspect of differentially expressed proteins under the conditions of different habitats in the field. In this study, two ecotypes of S. salsa from the intertidal (control) and supratidal (treatment) habitats of the Yellow River Delta were selected. A total of 30 individual leaves were mixed into six samples (three biological replicates for each) and subjected to protein extraction by using tandem mass tag-labeled quantitative proteomic technology. A total of 4771 proteins were quantitated. They included 317 differentially expressed proteins (2.0-fold change, p < 0.05), among which 143 were upregulated and the remaining 174 were downregulated. These differentially expressed proteins mainly participated in biological processes, such as response to stimulus, stress, and biotic stimulus; in molecular functions, such as methyltransferase activity, transferase activity, one-C group transfer, and tetrapyrrole binding; and in cell components, such as non-membrane-bound organelles, intracellular non-membrane-bound organelles, chromosomes, and photosystems. The differentially expressed proteins were mainly enriched in eight pathways, among which the ribosome, phenylpropanoid biosynthesis, and photosynthesis pathways had higher protein numbers than the other pathways. The upregulation of differentially expressed proteins related to the ribosome and photosynthesis increased the relative growth rate and reduced the N:P ratio of S. salsa from the supratidal habitat, thereby improving its palatability. By contrast, most of the differentially expressed proteins involved in phenylpropanoid biosynthesis were downregulated in S. salsa from the intertidal habitat. This result indicated that S. salsa from the intertidal habitat might accumulate flavonoids, lignin, and other secondary metabolites in its leaves that confer a bitter taste. However, these secondary metabolites might increase the medicinal value of S. salsa from the intertidal habitat. This work could provide a theoretical basis and data support for the sustainable and high-value utilization of medicinal and edible plants from coastal wetlands.

The evolution of the phenylpropanoid pathway entailed pronounced radiations and divergences of enzyme families

Land plants constantly respond to fluctuations in their environment. Part of their response is the production of a diverse repertoire of specialized metabolites. One of the foremost sources for metabolites relevant to environmental responses is the phenylpropanoid pathway, which was long thought to be a land-plant-specific adaptation shaped by selective forces in the terrestrial habitat. Recent data have, however, revealed that streptophyte algae, the algal relatives of land plants, have candidates for the genetic toolkit for phenylpropanoid biosynthesis and produce phenylpropanoid-derived metabolites. Using phylogenetic and sequence analyses, we here show that the enzyme families that orchestrate pivotal steps in phenylpropanoid biosynthesis have independently undergone pronounced radiations and divergence in multiple lineages of major groups of land plants; sister to many of these radiated gene families are streptophyte algal candidates for these enzymes. These radiations suggest a high evolutionary versatility in the enzyme families involved in the phenylpropanoid-derived metabolism across embryophytes. We suggest that this versatility likely translates into functional divergence, and may explain the key to one of the defining traits of embryophytes: a rich specialized metabolism.

Transcriptomics and Metabolomics Reveal Purine and Phenylpropanoid Metabolism Response to Drought Stress in Dendrobium sinense, an Endemic Orchid Species in Hainan Island

Drought stress is a bottleneck factor for plant growth and development, especially in epiphytic orchids that absorb moisture mainly from the air. Recent studies have suggested that there are complex transcriptional regulatory networks related to drought stress in Dendrobium sinense. In this study, the transcription and metabolite alterations involved in drought stress response in D. sinense were investigated through RNA-seq and metabolomics. A total of 856 metabolites were identified from stressed and control samples, with 391 metabolites showing significant differences. With PacBio and Illumina RNA sequencing, 72,969 genes were obtained with a mean length of 2,486 bp, and 622 differentially expressed genes (DEGs) were identified. Correlation analysis showed 7 differential genes, and 39 differential metabolites were involved in interaction networks. The network analysis of differential genes and metabolites suggested that the pathways of purine metabolism and phenylpropanoid biosynthesis may play an important role in drought response in D. sinense. These results provide new insights and reference data for culturally important medicinal plants and the protection of endangered orchids.

Pichia galeiformis Induces Resistance in Postharvest Citrus by Activating the Phenylpropanoid Biosynthesis Pathway

This study aimed to investigate the effect of Pichia galeiformis on disease resistance and elucidate the changes in phenylpropane biosynthesis treated by P. galeiformis in postharvest citrus. The results showed that P. galeiformis reduced the disease incidence and lesion diameters without direct contact with the pathogen Penicillium digitatum. Transcriptome analysis revealed that phenylpropanoid biosynthesis was triggered by P. galeiformis. Genes encoding phenylpropanoid biosynthesis were upregulated, including phenylalanine ammonia-lyase (PAL), 4-coumaroyl-CoA ligase (4CL), cinnamate-4-hydroxylase (C4H), peroxidase (POD), cinnamyl alcohol dehydrogenase (CAD), O-methyltransferase, and hydroxyl cinnamoyl transferase. Moreover, P. galeiformis increased the activity of PAL, 4CL, C4H, POD, polyphenol oxidase, and CAD in citrus pericarp. In addition, P. galeiformis treated citrus displayed higher levels of total phenolic compounds, flavonoid, and lignin and higher amounts of ferulic and sinapic acid. In conclusion, these results suggested that P. galeiformis could induce resistance through modulating the pathway of phenylpropanoid biosynthesis in postharvest citrus.

Combined Transcriptomic and Metabolomic Analysis Reveals the Role of Phenylpropanoid Biosynthesis Pathway in the Salt Tolerance Process of Sophora alopecuroides

Salt stress is the main abiotic stress that limits crop yield and agricultural development. Therefore, it is imperative to study the effects of salt stress on plants and the mechanisms through which plants respond to salt stress. In this study, we used transcriptomics and metabolomics to explore the effects of salt stress on Sophora alopecuroides. We found that salt stress incurred significant gene expression and metabolite changes at 0, 4, 24, 48, and 72 h. The integrated transcriptomic and metabolomic analysis revealed that the differentially expressed genes (DEGs) and differential metabolites (DMs) obtained in the phenylpropanoid biosynthesis pathway were significantly correlated under salt stress. Of these, 28 DEGs and seven DMs were involved in lignin synthesis and 23 DEGs and seven DMs were involved in flavonoid synthesis. Under salt stress, the expression of genes and metabolites related to lignin and flavonoid synthesis changed significantly. Lignin and flavonoids may participate in the removal of reactive oxygen species (ROS) in the root tissue of S. alopecuroides and reduced the damage caused under salt stress. Our research provides new ideas and genetic resources to study the mechanism of plant responses to salt stress and further improve the salt tolerance of plants.

Overexpression of SlGATA17 Promotes Drought Tolerance in Transgenic Tomato Plants by Enhancing Activation of the Phenylpropanoid Biosynthetic Pathway

GATA transcription factors (TFs) are widely distributed in eukaryotes. Some GATA TFs have been shown to be related to photosynthesis, germination, circadian rhythm, and other functions in plants. Our previous study found that some members of this family have obvious responses when tomato plants are subjected to drought stress, in which the SlGATA17 gene is significantly upregulated. To further verify the function of this gene under drought stress, we constructed tomato lines with this gene overexpressed. Phenotypic and physiological indicators indicated that the SlGATA17-overexpressing plants were more drought tolerant than the wild-type plants. Transcriptomic sequencing results showed that the overexpression of the SlGATA17 gene improved the activity of the phenylpropanoid biosynthesis pathway. The PAL enzyme activity assay results confirmed that the initial activity of this pathway was enhanced in transgenic plants, especially in the initial response stage, indicating that the SlGATA17 gene regulates the drought resistance of tomato plants by regulating the activity of the phenylpropanoid biosynthesis pathway.

Bacillus pumilus LZP02 Promotes Rice Root Growth by Improving Carbohydrate Metabolism and Phenylpropanoid Biosynthesis

Elucidation of the underlying mechanisms of plant growth promotion of rhizobacteria is very important. This study explored the mechanism by which Bacillus pumilus LZP02 promotes growth in rice roots through proteomic, transcriptomic, and metabolomic techniques. The results showed that B. pumilus LZP02 promoted the absorption of phosphorous, calcium, and magnesium ions by colonization of rice roots and enhanced peroxidase, catalase, superoxide dismutase, and Ca²⁺Mg²⁺ adenosine triphosphatase activities and chlorophyll contents in rice. The proteomic results showed that most of the differentially expressed proteins were involved in carbohydrate metabolism and that the biosynthesis of other secondary metabolites was also increased. According to RNA-seq and reverse transcription-quantitative PCR analyses, expression of some genes involved in carbohydrate metabolism and phenylpropanoid biosynthesis was upregulated in rice roots. Regarding metabolomics, phenylpropanoid biosynthesis, starch and sucrose metabolism, the pentose phosphate pathway, and glyoxylate and dicarboxylate metabolism were increased. The results indicated that B. pumilus LZP02 promoted the growth of rice roots by enhancing carbohydrate metabolism and phenylpropanoid biosynthesis.

Barley Root Proteome and Metabolome in Response to Cytokinin and Abiotic Stimuli

Cytokinin is a phytohormone involved in the regulation of diverse developmental and physiological processes in plants. Its potential in biotechnology and for development of higher-yield and more resilient plants has been recognized, yet the molecular mechanisms behind its action are far from understood. In this report, the roots of barley seedlings were explored as a new source to reveal as yet unknown cytokinin-responsive proteins for crop improvement. Here we found significant differences reproducibly observed for 178 proteins, for which some of the revealed cytokinin-responsive pathways were confirmed in metabolome analysis, including alterations phenylpropanoid pathway, amino acid biosynthesis and ROS metabolism. Bioinformatics analysis indicated a significant overlap between cytokinin response and response to abiotic stress. This was confirmed by comparing proteome and metabolome profiles in response to drought, salinity or a period of temperature stress. The results illustrate complex abiotic stress response in the early development of model crop plant and confirm an extensive crosstalk between plant hormone cytokinin and response to temperature stimuli, water availability or salinity stress.

Maize phenylalanine ammonia-lyases contribute to resistance to Sugarcane mosaic virus infection, most likely through positive regulation of salicylic acid accumulation

Sugarcane mosaic virus (SCMV) is a pathogen of worldwide importance that causes dwarf mosaic disease on maize (Zea mays). Until now, few maize genes/proteins have been shown to be involved in resistance to SCMV. In this study, we characterized the role of maize phenylalanine ammonia-lyases (ZmPALs) in accumulation of the defence signal salicylic acid (SA) and in resistance to virus infection. SCMV infection significantly increased SA accumulation and expression of SA-responsive pathogenesis-related protein genes (PRs). Interestingly, exogenous SA treatment decreased SCMV accumulation and enhanced resistance. Both reverse transcription-coupled quantitative PCR and RNA-Seq data confirmed that expression levels of at least four ZmPAL genes were significantly up-regulated upon SCMV infection. Knockdown of ZmPAL expression led to enhanced SCMV infection symptom severity and virus multiplication, and simultaneously resulted in decreased SA accumulation and PR gene expression. Intriguingly, application of exogenous SA to SCMV-infected ZmPAL-silenced maize plants decreased SCMV accumulation, showing that ZmPALs are required for SA-mediated resistance to SCMV infection. In addition, lignin measurements and metabolomic analysis showed that ZmPALs are also involved in SCMV-induced lignin accumulation and synthesis of other secondary metabolites via the phenylpropanoid pathway. In summary, our results indicate that ZmPALs are required for SA accumulation in maize and are involved in resistance to virus infection by limiting virus accumulation and moderating symptom severity.

Proteomics Analysis to Identify Proteins and Pathways Associated with the Novel Lesion Mimic Mutant E40 in Rice Using iTRAQ-Based Strategy

A novel rice lesion mimic mutant (LMM) was isolated from the mutant population of Japonica rice cultivar Hitomebore generated by ethyl methane sulfonate (EMS) treatment. Compared with the wild-type (WT), the mutant, tentatively designated E40, developed necrotic lesions over the whole growth period along with detectable changes in several important agronomic traits including lower height, fewer tillers, lower yield, and premature death. To understand the molecular mechanism of mutation-induced phenotypic differences in E40, a proteomics-based approach was used to identify differentially accumulated proteins between E40 and WT. Proteomic data from isobaric tags for relative and absolute quantitation (iTRAQ) showed that 233 proteins were significantly up- or down-regulated in E40 compared with WT. These proteins are involved in diverse biological processes, but phenylpropanoid biosynthesis was the only up-regulated pathway. Differential expression of the genes encoding some candidate proteins with significant up- or down-regulation in E40 were further verified by qPCR. Consistent with the proteomic results, substance and energy flow in E40 shifted from basic metabolism to secondary metabolism, mainly phenylpropanoid biosynthesis, which is likely involved in the formation of leaf spots.

Regulation of resveratrol biosynthesis in grapevine: new approaches for disease resistance?

This article comments on: Jiang J, Xi H, Dai Z, Lecourieux F, Yuan L, Liu X, Patra B, Wei Y, Li S, Wang L. 2019. VvWRKY8 negatively regulates VvSTS through direct interaction with VvMYB14 to balance resveratrol biosynthesis in grapevine. Journal of Experimental Botany 70, 715–729.

Identification of Candidate Genes and Biosynthesis Pathways Related to Fertility Conversion by Wheat KTM3315A Transcriptome Profiling

The Aegilops kotschyi thermo-sensitive cytoplasmic male sterility (K-TCMS) system may facilitate hybrid wheat (Triticum aestivum L.) seed multiplication and production. The K-TCMS line is completely male sterile during the normal wheat-growing season, whereas its fertility can be restored in a high-temperature environment. To elucidate the molecular mechanisms responsible for male sterility/fertility conversion and candidate genes involved with pollen development in K-TCMS, we employed RNA-seq to sequence the transcriptomes of anthers from K-TCMS line KTM3315A during development under sterile and fertile conditions. We identified 16840 differentially expressed genes (DEGs) in different stages including15157 known genes (15135 nuclear genes and 22 plasmagenes) and 1683 novel genes. Bioinformatics analysis identified possible metabolic pathways involved with fertility based on KEGG pathway enrichment of the DEGs expressed in fertile and sterile plants. We found that most of the genes encoding key enzyme in the phenylpropanoid biosynthesis and jasmonate biosynthesis pathways were significant upregulated in uninucleate, binuclate or trinucleate stage, which both interact with MYB transcription factors, and that link between all play essential roles in fertility conversion. The relevant DEGs were verified by quantitative RT-PCR. Thus, we suggested that phenylpropanoid biosynthesis and jasmonate biosynthesis pathways were involved in fertility conversion of K-TCMS wheat. This will provide a new perspective and an effective foundation for the research of molecular mechanisms of fertility conversion of CMS wheat. Fertility conversion mechanism in thermo-sensitive cytoplasmic male sterile/fertile wheat involves the phenylpropanoid biosynthesis pathway, jasmonate biosynthesis pathway, and MYB transcription factors.