Downregulation of caffeic acid 3-O-methyltransferase and caffeoyl CoA 3-O-methyltransferase in transgenic alfalfa. impacts on lignin structure and implications for the biosynthesis of G and S lignin

Functional analysis of two caffeoyl-coenzyme 3 a-o-methyltransferase involved in pear lignin metabolism

Caffeoyl-coenzyme 3 A-O-methyltransferase (CCoAOMT) plays a crucial role in the lignin synthesis in many higher plants. In this study, nine PbCCoAOMT genes in total were identified from pear, and classified into six categories. We treated pear fruits with hormones abscisic acid (ABA) and methyl jasmonate (MeJA) and salicylic acid (SA) and observed differential expression levels of these genes. Through qRT-PCR, we also preliminarily identified candidate PbCCoAOMT gene, potentially involved in lignin synthesis in pear fruits. Additionally, the overexpression of PbCCoAOMT1/2 in Arabidopsis and pear fruits increased in lignin content. Enzymatic assays showed that recombinant PbCCoAOMT1/2 proteins have similar enzymatic activity in vitro. The Y1H (Yeast one-hybrid) and dual luciferase (dual-LUC) experiments demonstrated that PbMYB25 can bind to the AC elements in the promoter region of the PbCCoAOMT1 gene. Our findings suggested that the PbCCoAOMT1 and PbCCoAOMT2 genes may contribute to the synthesis of lignin and provide insights into the mechanism of lignin biosynthesis and stone cell development in pear fruits.

Exploring the leaf regeneration cycles response of Zostera japonica to ocean acidification

Ocean acidification is one of the major global environmental problems facing humankind today, and it has far-reaching impacts on marine organisms and the entire marine ecosystem. Zostera japonica, an important supporting species of intertidal seagrass beds, exhibits high photosynthetic productivity and plays an important role in the carbon cycle of nearshore waters. However, little is known about the characteristics, processes, and mechanisms of its response to ocean acidification. In this study, we conducted a 120-day acidification experiment in Z. japonica; here, plants underwent four leaf regeneration cycles to reveal the response mechanism of Z. japonica to ocean acidification (OA). We found that acidification significantly affected the seedling stage of Z. japonica, impacting leaf regeneration cycles by altering physiological and molecular responses. In one leaf regeneration cycle, the short-term exposure to CO2 affected the seagrass parameters, such as the regulation of inorganic carbon uptake modes and the regulation of photosynthesis between the dark and light reactions, with the potential to affect the carbon sinks of the marine organisms. The long-term effects on the regulation of antioxidant enzymes and antioxidant metabolites, caused an improvement in the marine life adaptation to OA. In a comparison of the different leaf regeneration cycles, the response pattern of Z. japonica showed an offset of the acidification during the short cycles and an adaption to the acidification during the long cycles. This study revealed the response mechanism of Z. japonica to OA at different time scales and could provide a theoretical basis for accurately assessing the impact of OA on seagrass and the entire seagrass bed ecosystem.

Systematic Identification and Characterization of O-Methyltransferase Gene Family Members Involved in Flavonoid Biosynthesis in Chrysanthemum indicum L

Chrysanthemum indicum L. capitulum is an enriched source of flavonoids with broad-ranging biological activities, mainly due to their anti-inflammatory, anti-cancer, immune regulation, anti-microbial activity, hepatoprotective, and neuroprotective effects. The O-methylation of various secondary metabolites has previously been demonstrated to be mainly catalyzed by S-adenosyl-L-methionine-dependent O-methyltransferase (OMT) proteins encoded by the OMT gene family. However, limited comprehensive study was published on the OMT gene family, especially the CCoAOMT subfamily, involved in the O-methylation of flavonoids in Chrysanthemum. Here, we analyzed the spatiotemporal expression patterns of C. indicum OMT genes in leaf and flower at different developmental stages. Transcriptome sequencing and qRT-PCR analysis showed that COMTs were mainly highly expressed in capitulum, especially in full bloom, while CCoAOMTs were mainly highly expressed in leaves. Correlation analysis of OMT gene expression and flavonoids accumulation revealed that four OMTs (CHR00029120, CHR00029783, CHR00077404, and CHR00078333) were putatively involved in most methylated flavonoids biosynthesis in the capitulum. Furthermore, we identified a true CCoAOMT enzyme, CiCCoAOMT1, and found that it catalyzed O-methylation of quercetin and luteolin at the 3'-OH position. In summary, this work provides an important theoretical basis for further research on the biological functions of OMTs in C. indicum.

Exploring Lignin Biosynthesis Genes in Rice: Evolution, Function, and Expression

Lignin is nature's second most abundant vascular plant biopolymer, playing significant roles in mechanical support, water transport, and stress responses. This study identified 90 lignin biosynthesis genes in rice based on phylogeny and motif constitution, and they belong to PAL, C4H, 4CL, HCT, C3H, CCoAOMT, CCR, F5H, COMT, and CAD families. Duplication events contributed largely to the expansion of these gene families, such as PAL, CCoAOMT, CCR, and CAD families, mainly attributed to tandem and segmental duplication. Microarray data of 33 tissue samples covering the entire life cycle of rice suggested fairly high PAL, HCT, C3H, CCoAOMT, CCR, COMT, and CAD gene expressions and rather variable C4H, 4CL, and F5H expressions. Some members of lignin-related genes (OsCCRL11, OsHCT1/2/5, OsCCoAOMT1/3/5, OsCOMT, OsC3H, OsCAD2, and OsPAL1/6) were expressed in all tissues examined. The expression patterns of lignin-related genes can be divided into two major groups with eight subgroups, each showing a distinct co-expression in tissues representing typically primary and secondary cell wall constitutions. Some lignin-related genes were strongly co-expressed in tissues typical of secondary cell walls. Combined HPLC analysis showed increased lignin monomer (H, G, and S) contents from young to old growth stages in five genotypes. Based on 90 genes' microarray data, 27 genes were selected for qRT-PCR gene expression analysis. Four genes (OsPAL9, OsCAD8C, OsCCR8, and OsCOMTL4) were significantly negatively correlated with lignin monomers. Furthermore, eleven genes were co-expressed in certain genotypes during secondary growth stages. Among them, six genes (OsC3H, OsCAD2, OsCCR2, OsCOMT, OsPAL2, and OsPAL8) were overlapped with microarray gene expressions, highlighting their importance in lignin biosynthesis.

The MYB transcription factor PpMYB5 regulates Pp4CL1/Pp4CL2 expression to promote lignin biosynthesis of fruit russeting in the flat nectarine

KEY MESSAGE

Transcription factor PpMYB5 promotes lignin synthesis by directly binding to the Pp4CL1/Pp4CL2 promoter and affecting their expression, which may be related to nectarine russeting formation. Nectarine russeting is usually considered to be a non-invasive physiological disease that usually occurs on late-maturing cultivars and seriously affects their appearance quality and commercial value. The cause of nectarine fruit rust is currently unknown. In this study, we compared two flat nectarine cultivars, 'zhongyoupanweidi' (HD; russeting-free cultivar) and 'zhongyoupanweihou' (TH; russeting-prone cultivar), with respect to nectarine russeting by means of microscopy, transcriptomics, and hormone analysis. Compared to HD fruits, TH fruits had a broken cuticle, missing wax layer, and heavy lignin deposition. RNA sequencing (RNA-seq) revealed significant alternations in the expression of genes related to lignin synthesis. Moreover, structure genes Pp4CL1 and Pp4CL2, MYB transcription factor (TF) gene PpMYB5 were identified through weighted gene co-expression network analysis (WGCNA). Molecular experiments and transgenic evidence suggested that PpMYB5 regulates Pp4CL1/Pp4CL2 expression to promote lignin synthesis. Overall, in addition to providing new insights into the formation of mechanisms for nectarine russeting, our study also establishes a foundation for nectarine russeting prevention.

Transcription Factor VlbZIP14 Inhibits Postharvest Grape Berry Abscission by Directly Activating VlCOMT and Promoting Lignin Biosynthesis

Sulfur dioxide (SO2) is the most effective preservative for table grapes as it reduces the respiratory intensity of berries and inhibits mold growth. However, excessive SO2 causes berry abscission during storage, resulting in an economic loss postharvest. In this study, grapes were exogenously treated with SO2, SO2 + 1.5% chitosan, SO2 + 1.5% eugenol, and SO2 + eugenol-loaded chitosan nanoparticles (SN). In comparison to SO2 treatment, SN treatment reduced the berries' abscission rate by 74% while maintaining the quality of the berries. Among the treatments, SN treatment most effectively inhibited berry abscission and maintained berry quality. RNA-sequencing (RNA-seq) revealed that SN treatment promoted the expression of genes related to cell wall metabolism. Among these genes, VlCOMT was detected as the central gene, playing a key role in mediating the effects of SN. Dual luciferase and yeast one-hybrid (Y1H) assays demonstrated that VlbZIP14 directly activated VlCOMT by binding to the G-box motif in the latter's promoter, which then participated in lignin synthesis. Our results provide key insights into the molecular mechanisms underlying the SN-mediated inhibition of berry abscission and could be used to improve the commercial value of SO2-treated postharvest table grapes.

PagARGOS promotes low-lignin wood formation in poplar

Wood formation, which occurs mainly through secondary xylem development, is important not only for supplying raw material for the 'ligno-chemical' industry but also for driving the storage of carbon. However, the complex mechanisms underlying the promotion of xylem formation remain to be elucidated. Here, we found that overexpression of Auxin-Regulated Gene involved in Organ Size (ARGOS) in hybrid poplar 84 K (Populus alba × Populus tremula var. glandulosa) enlarged organ size. In particular, PagARGOS promoted secondary growth of stems with increased xylem formation. To gain further insight into how PagARGOS regulates xylem development, we further carried out yeast two-hybrid screening and identified that the auxin transporter WALLS ARE THIN1 (WAT1) interacts with PagARGOS. Overexpression of PagARGOS up-regulated WAT1, activating a downstream auxin response promoting cambial cell division and xylem differentiation for wood formation. Moreover, overexpressing PagARGOS caused not only higher wood yield but also lower lignin content compared with wild-type controls. PagARGOS is therefore a potential candidate gene for engineering fast-growing and low-lignin trees with improved biomass production.

Breeding for improved digestibility and processing of lignocellulosic biomass in Zea mays

Forage maize is a versatile crop extensively utilized for animal nutrition in agriculture and holds promise as a valuable resource for the production of fermentable sugars in the biorefinery sector. Within this context, the carbohydrate fraction of the lignocellulosic biomass undergoes deconstruction during ruminal digestion and the saccharification process. However, the cell wall's natural resistance towards enzymatic degradation poses a significant challenge during both processes. This so-called biomass recalcitrance is primarily attributed to the presence of lignin and ferulates in the cell walls. Consequently, maize varieties with a reduced lignin or ferulate content or an altered lignin composition can have important beneficial effects on cell wall digestibility. Considerable efforts in genetic improvement have been dedicated towards enhancing cell wall digestibility, benefiting agriculture, the biorefinery sector and the environment. In part I of this paper, we review conventional and advanced breeding methods used in the genetic improvement of maize germplasm. In part II, we zoom in on maize mutants with altered lignin for improved digestibility and biomass processing.

Genome-Wide Association Study Reveals the Genetic Basis of Crude Fiber Components in Brassica napus L. Shoots at Stem Elongation Stage

Brassica napus is currently the principal field crop for producing materials for primary, secondary and tertiary industries. B. napus shoots at stem elongation stage are rich in anthocyanins, vitamin C and mineral elements such as selenium, calcium and zinc, and represent a new type of green vegetable. However, the high crude fiber (CF) content of B. napus shoots affects their taste, and few studies have focused on the quality traits of these vegetables. In this study, we investigated five traits related to the CF components, including neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), hemicellulose (Hem) and cellulose (Cel), of B. napus shoots. Whole-genome resequencing at a depth of ∼20× was utilized to genotype an association panel of 202 diverse accessions, which resulted in the identification of 6,093,649 single nucleotide polymorphisms (SNPs) and 996,252 indels, respectively. A genome-wide association study (GWAS) was performed for the five CF-related traits based on the phenotypic data observed in four environments. A total of 1,285 significant SNPs were detected at the threshold of -log10 (p) = 5.16, and 97 significant association regions were obtained. In addition, seven candidate genes located on chromosomes A2 (one gene), A8 (three genes), A9 (two genes) and C9 (one gene) related to CF traits were identified, and ten lines containing low CF contents were selected as excellent germplasm resources for breeding. Our results contributed new insights into the genetic basis of CF traits and suggested germplasm resources for the quality improvement of B. napus shoots.

Biosynthesis and metabolic engineering of isoflavonoids in model plants and crops: a review

Isoflavonoids, the major secondary metabolites within the flavonoid biosynthetic pathway, play important roles in plant defense and exhibit free radical scavenging properties in mammals. Recent advancements in understanding the synthesis, transport, and regulation of isoflavonoids have identified their biosynthetic pathways as promising targets for metabolic engineering, offering potential benefits such as enhanced plant resistance, improved biomass, and restoration of soil fertility. This review provides an overview of recent breakthroughs in isoflavonoid biosynthesis, encompassing key enzymes in the biosynthetic pathway, transporters influencing their subcellular localization, molecular mechanisms regulating the metabolic pathway (including transcriptional and post-transcriptional regulation, as well as epigenetic modifications). Metabolic engineering strategies aimed at boosting isoflavonoid content in both leguminous and non-leguminous plants. Additionally, we discuss emerging technologies and resources for precise isoflavonoid regulation. This comprehensive review primarily focuses on model plants and crops, offering insights for more effective and sustainable metabolic engineering approaches to enhance nutritional quality and stress tolerance.

Evolution and Analysis of Caffeic Acid Transferase (COMT) in Seed Plants

Caffeic acid transferase (COMT) is a key enzyme in the lignin and melatonin synthesis pathways and plays an important role in plant growth and development. All seed plants have two characteristics: they have vascular tissues, phloem, and xylem, and they can produce and reproduce seeds. In order to understand the distribution and evolution of COMTs in seed plants, we performed physicochemical property analysis, subcellular localization, phylogenetic analysis, conserved motif analysis, and protein interaction network analysis of 44 COMT homologs from 26 seed plants through in silico. The results showed that in seed plants, the structure of COMT genes tends to be stable in different plant taxa, while the relationship between the chromosomal positions of different COMT genes in the same plant was more intricate. The conserved distribution of COMT in seed plants reflected its highly specialized function.

Integrative Dissection of Lignin Composition in Tartary Buckwheat Seed Hulls for Enhanced Dehulling Efficiency

The rigid hull encasing Tartary buckwheat seeds necessitates a laborious dehulling process before flour milling, resulting in considerable nutrient loss. Investigation of lignin composition is pivotal in understanding the structural properties of tartary buckwheat seeds hulls, as lignin is key determinant of rigidity in plant cell walls, thus directly impacting the dehulling process. Here, the lignin composition of seed hulls from 274 Tartary buckwheat accessions is analyzed, unveiling a unique lignin chemotype primarily consisting of G lignin, a common feature in gymnosperms. Furthermore, the hardness of the seed hull showed a strong negative correlation with the S lignin content. Genome-wide detection of selective sweeps uncovered that genes governing the biosynthesis of S lignin, specifically two caffeic acid O-methyltransferases (COMTs) and one ferulate 5-hydroxylases, are selected during domestication. This likely contributed to the increased S lignin content and decreased hardness of seed hulls from more domesticated varieties. Genome-wide association studies identified robust associations between FtCOMT1 and the accumulation of S lignin in seed hull. Transgenic Arabidopsis comt1 plants expressing FtCOMT1 successfully reinstated S lignin content, confirming its conserved function across plant species. These findings provide valuable metabolic and genetic insights for the potential redesign of Tartary buckwheat seed hulls.

Overexpression of MsDREB1C Modulates Growth and Improves Forage Quality in Tetraploid Alfalfa (Medicago sativa L.)

DREB has been reported to be involved in plant growth and response to environmental factors. However, the function of DREB in growth and development has not been elucidated in alfalfa (Medicago sativa L.), a perennial tetraploid forage cultivated worldwide. In this study, an ortholog of MtDREB1C was characterized from alfalfa and named MsDREB1C accordingly. MsDREB1C was significantly induced by abiotic stress. The transcription factor MsDREB1C resided in the nucleus and had self-transactivation activity. The MsDREB1C overexpression (OE) alfalfa displayed growth retardation under both long-day and short-day conditions, which was supported by decreased MsGA20ox and upregulated MsGA2ox in the OE lines. Consistently, a decrease in active gibberellin (GA) was detected, suggesting a negative effect of MsDREB1C on GA accumulation in alfalfa. Interestingly, the forage quality of the OE lines was better than that of WT lines, with higher crude protein and lower lignin content, which was supported by an increase in the leaf-stem ratio (LSR) and repression of several lignin-synthesis genes (MsNST, MsPAL1, MsC4H, and Ms4CL). Therefore, this study revealed the effects of MsDREB1C overexpression on growth and forage quality via modifying GA accumulation and lignin synthesis, respectively. Our findings provide a valuable candidate for improving the critical agronomic traits of alfalfa, such as overwintering and feeding value of the forage.

Multifunctional 5-hydroxyconiferaldehyde O-methyltransferases (CAldOMTs) in plant metabolism

Lignin, flavonoids, melatonin, and stilbenes are plant specialized metabolites with diverse physiological and biological functions, supporting plant growth and conferring stress resistance. Their biosynthesis requires O-methylations catalyzed by 5-hydroxyconiferaldehyde O-methyltransferase (CAldOMT; also called caffeic acid O-methyltransferase, COMT). CAldOMT was first known for its roles in syringyl (S) lignin biosynthesis in angiosperm cell walls and later found to be multifunctional. This enzyme also catalyzes O-methylations in flavonoid, melatonin, and stilbene biosynthetic pathways. Phylogenetic analysis indicated the convergent evolution of enzymes with OMT activities towards the monolignol biosynthetic pathway intermediates in some gymnosperm species that lack S-lignin and Selaginella moellendorffii, a lycophyte which produces S-lignin. Furthermore, neofunctionalization of CAldOMTs occurred repeatedly during evolution, generating unique O-methyltransferases (OMTs) with novel catalytic activities and/or accepting novel substrates, including lignans, 1,2,3-trihydroxybenzene, and phenylpropenes. This review summarizes multiple aspects of CAldOMTs and their related proteins in plant metabolism and discusses their evolution, molecular mechanism, and roles in biorefineries, agriculture, and synthetic biology.

Genome-wide identification, expression profiling, and protein interaction analysis of the CCoAOMT gene family in the tea plant (Camellia sinensis)

BACKGROUND

The caffeoyl-CoA-O methyltransferase (CCoAOMT) family plays a crucial role in the oxidative methylation of phenolic substances and is involved in various plant processes, including growth, development, and stress response. However, there is a limited understanding of the interactions among CCoAOMT protein members in tea plants.

RESULTS

In this study, we identified 10 members of the CsCCoAOMT family in the genome of Camellia sinensis (cultivar 'HuangDan'), characterized by conserved gene structures and motifs. These CsCCoAOMT members were located on six different chromosomes (1, 2, 3, 4, 6, and 14). Based on phylogenetic analysis, CsCCoAOMT can be divided into two groups: I and II. Notably, the CsCCoAOMT members of group Ia are likely to be candidate genes involved in lignin biosynthesis. Moreover, through the yeast two-hybrid (Y2H) assay, we established protein interaction networks for the CsCCoAOMT family, revealing 9 pairs of members with interaction relationships.

CONCLUSIONS

We identified the CCoAOMT gene family in Camellia sinensis and conducted a comprehensive analysis of their classifications, phylogenetic and synteny relationships, gene structures, protein interactions, tissue-specific expression patterns, and responses to various stresses. Our findings shed light on the evolution and composition of CsCCoAOMT. Notably, the observed interaction among CCoAOMT proteins suggests the potential formation of the O-methyltransferase (OMT) complex during the methylation modification process, expanding our understanding of the functional roles of this gene family in diverse biological processes.

Lignin reduction in sugarcane by performing CRISPR/Cas9 site-direct mutation of SoLIM transcription factor

Genetic engineering of plant cell walls is limited for reducing lignocellulose recalcitrance, so mild and/or green-like pretreatment is still required for sequential enzymatic saccharification. Here, we report a method to reduce lignin content in sugarcane stalks using the CRISPR/Cas 9 technique. Three target sequences of SoLIM were designed and fused to pRGEB32. The cassette constructs were introduced into sugarcane calli cv. KK3 through Agrobacterium-mediated transformation. We produced one base substitution and one insertion line for the 1st target site; two insertions, one deletion, and one base substitution for the 2nd target site; and one base substitution and insertion for the 3rd target site. qRT-PCR analysis of SoLIM, SoPAL, SoC4H, and SoCAD showeded that downregulation of SoLIM by single nucleotide insertions or deletions reduced the expression of SoPAL, SoC4H, and SoCAD. Consequently, the edited lines contained 9.74 to 51.46% less lignin content compared to that in the wild-type plants. The syringyl/guaiacyl (S/G) ratio of the edited lines ranged between 0.23 and 0.49, while the wild-type was 0.22. The histochemical evaluation and scanning electron microscopy of the cell walls supported this observation. A low lignin content sugarcane will provide a better feedstock for second-generation bioethanol production.

Transcriptional alterations of peanut root during interaction with growth-promoting Tsukamurella tyrosinosolvens strain P9

The plant growth-promoting rhizobacterium Tsukamurella tyrosinosolvens P9 can improve peanut growth. In this study, a co-culture system of strain P9 and peanut was established to analyze the transcriptome of peanut roots interacting with P9 for 24 and 72 h. During the early stage of co-culturing, genes related to mitogen-activated protein kinase (MAPK) and Ca2+ signal transduction, ethylene synthesis, and cell wall pectin degradation were induced, and the up-regulation of phenylpropanoid derivative, flavonoid, and isoflavone synthesis enhanced the defense response of peanut. The enhanced expression of genes associated with photosynthesis and carbon fixation, circadian rhythm regulation, indoleacetic acid (IAA) synthesis, and cytokinin decomposition promoted root growth and development. At the late stage of co-culturing, ethylene synthesis was reduced, whereas Ca2+ signal transduction, isoquinoline alkaloid synthesis, and ascorbate and aldarate metabolism were up-regulated, thereby maintaining root ROS homeostasis. Sugar decomposition and oxidative phosphorylation and nitrogen and fatty acid metabolism were induced, and peanut growth was significantly promoted. Finally, the gene expression of seedlings inoculated with strain P9 exhibited temporal differences. The results of our study, which explored transcriptional alterations of peanut root during interacting with P9, provide a basis for elucidating the growth-promoting mechanism of this bacterial strain in peanut.

Targeted mutagenesis of flavonoid biosynthesis pathway genes reveals functional divergence in seed coat colour, oil content and fatty acid composition in Brassica napus L

Yellow-seed is widely accepted as a good-quality trait in Brassica crops. Previous studies have shown that the flavonoid biosynthesis pathway is essential for the development of seed colour, but its function in Brassica napus, an important oil crop, is poorly understood. To systematically explore the gene functions of the flavonoid biosynthesis pathway in rapeseed, several representative TRANSPARENT TESTA (TT) genes, including three structural genes (BnaTT7, BnaTT18, BnaTT10), two regulatory genes (BnaTT1, BnaTT2) and a transporter (BnaTT12), were selected for targeted mutation by CRISPR/Cas9 in the present study. Seed coat colour, lignin content, seed quality and yield-related traits were investigated in these Bnatt mutants together with Bnatt8 generated previously. These Bnatt mutants produced seeds with an elevated seed oil content and decreased pigment and lignin accumulation in the seed coat without any serious defects in the yield-related traits. In addition, the fatty acid (FA) composition was also altered to different degrees, i.e., decreased oleic acid and increased linoleic acid and α-linolenic acid, in all Bnatt mutants except Bnatt18. Furthermore, gene expression analysis revealed that most of BnaTT mutations resulted in the down-regulation of key genes related to flavonoid and lignin synthesis, and the up-regulation of key genes related to lipid synthesis and oil body formation, which may contribute to the phenotype. Collectively, our study generated valuable resources for breeding programs, and more importantly demonstrated the functional divergence and overlap of flavonoid biosynthesis pathway genes in seed coat colour, oil content and FA composition of rapeseed.

Investigation of the response of Platycodongrandiflorus (Jacq.) A. DC to salt stress using combined transcriptomics and metabolomics

BACKGROUND

Platycodon grandiflorus (Jacq.) A. DC is a famous traditional Chinese medicine in China and an authentic medicine in Inner Mongolia. It has been traditionally used as an expectorant in cough and also has anti-inflammatory and other pharmacological effects. As a homologous plant of medicine and food, P. grandiflorus is widely planted in Northeast China. Soil salinity isa limiting factor for its cultivation. In this study, we comprehensively described the physiological characteristics of P. grandiflorus and combined transcriptomics and metabolomics to study the response of roots of P. grandiflorus to salt stress.

RESULTS

Overall, 8,988 differentially expressed genes were activated and significantly altered the metabolic processes. In total, 428 differentially abundant metabolites were affected by salt stress. After moderate and severe salt stress, most of the differentially abundant metabolites were enriched in the L-phenylalanine metabolic pathway. Through the comprehensive analysis of the interaction between key genes and metabolites, the main pathways such as lignin compound biosynthesis and triterpene saponin biosynthesis were completed. The relative content of compounds related to lignin biosynthesis, such as caffeic acid, coniferin, and syringing, increased under salt stress, and the related genes such as PAL, C4H, and the key enzyme gene UGT72E2 were activated to adapt to the salt stress. Platycodon saponin is one of the major triterpene saponins in P. grandiflorus, and Platycodin D is its most abundant major bioactive component. Under severe salt stress, Platycodin D level increased by nearly 1.77-fold compared with the control group. Most of the genes involved insynthetic pathway of Platycodin D, such as HMGCR, GGPS, SE, and LUP, were upregulated under salt stress.

CONCLUSION

Salt stress led to a decrease in the biomass and affected the activities of antioxidant enzymes and contents of osmotic regulators in the plant. These results provided not only novel insights into the underlying mechanisms of response of P. grandiflorus to salt stress but also a foundation for future studies on the function of genes related to salt tolerance in the triterpenoid saponin biosynthesis pathway.

Two QTL regions for spike length showing pleiotropic effects on Fusarium head blight resistance and thousand-grain weight in bread wheat (Triticum aestivum L.)

Spike length (SL) plays an important role in the yield improvement of wheat and is significantly associated with other traits. Here, we used a recombinant inbred line (RIL) population derived from a cross between Yangmai 12 (YM12) and Yanzhan 1 (YZ1) to construct a genetic linkage map and identify quantitative trait loci (QTL) for SL. A total of 5 QTL were identified for SL, among which QSl.yaas-3A and QSl.yaas-5B are two novel QTL for SL. The YZ1 alleles at QSl.yaas-2D and QSl.yaas-5A, and the YM12 alleles at QSl.yaas-2A, QSl.yaas-3A, and QSl.yaas-5B conferred increasing SL effects. Two major QTL QSl.yaas-5A and QSl.yaas-5B explained 9.11-15.85% and 9.01-12.85% of the phenotypic variations, respectively. Moreover, the positive alleles of QSl.yaas-5A and QSl.yaas-5B could significantly increase Fusarium head blight (FHB) resistance (soil surface inoculation and spray inoculation were used) and thousand-grain weight (TGW) in the RIL population. Kompetitive allele-specific PCR (KASP) markers for QSl.yaas-5A and QSl.yaas-5B were developed and validated in an additional panel of 180 wheat cultivars/lines. The cultivars/lines harboring both the positive alleles of QSl.yaas-5A and QSl.yaas-5B accounted for only 28.33% of the validation populations and had the longest SL, best FHB resistance (using spray inoculation), and highest TGW. A total of 358 and 200 high-confidence annotated genes in QSl.yaas-5A and QSl.yaas-5B were identified, respectively. Some of the genes in these two regions were involved in cell development, disease resistance, and so on. The results of this study will provide a basis for directional breeding of longer SL, higher TGW, and better FHB resistance varieties and a solid foundation for fine-mapping QSl.yaas-5A and QSl.yaas-5B in future.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01427-8.

A bacterial type III effector targets plant vesicle-associated membrane proteins

The soilborne bacterial pathogen Ralstonia solanacearum is one of the most destructive plant pathogens worldwide, and its infection process involves the manipulation of numerous plant cellular functions. In this work, we found that the R. solanacearum effector protein RipD partially suppressed different levels of plant immunity triggered by R. solanacearum elicitors, including specific responses triggered by pathogen-associated molecular patterns and secreted effectors. RipD localized in different subcellular compartments in plant cells, including vesicles, and its vesicular localization was enriched in cells undergoing R. solanacearum infection, suggesting that this specific localization may be particularly relevant during infection. Among RipD-interacting proteins, we identified plant vesicle-associated membrane proteins (VAMPs). We also found that overexpression of Arabidopsis thaliana VAMP721 and VAMP722 in Nicotiana benthamiana leaves promoted resistance to R. solanacearum, and this was abolished by the simultaneous expression of RipD, suggesting that RipD targets VAMPs to contribute to R. solanacearum virulence. Among proteins secreted in VAMP721/722-containing vesicles, CCOAOMT1 is an enzyme required for lignin biosynthesis, and mutation of CCOAOMT1 enhanced plant susceptibility to R. solanacearum. Altogether our results reveal the contribution of VAMPs to plant resistance against R. solanacearum and their targeting by a bacterial effector as a pathogen virulence strategy.

Osmotic stress-induced lignin synthesis is regulated at multiple levels in alfalfa (Medicago sativa L.)

Alfalfa is an important forage crop. Yield and quality are frequently threatened by extreme environments such as drought and salt stress. As a component of the cell wall, lignin plays an important role in the abiotic stress response, the mechanisms of which have not been well clarified. In this study, we combined physiological, transcriptional, and metabolic analyses to reveal the changes in lignin content in alfalfa under mannitol-induced osmotic stress. Osmotic stress enhanced lignin accumulation by increasing G and S units, which was associated with increases in enzyme activities and decreases in 8 intermediate metabolites. Upon combined analysis of the transcriptome and metabolome, we identified five key structural genes and several coexpressed transcription factors, such as MYB and WRKY, which may play a core role in regulating lignin content and composition under osmotic stress. In addition, lignin synthesis was positively regulated by ABA but negatively regulated by ethylene under osmotic stress. These results provide new insight into the regulatory mechanism of lignin synthesis under abiotic stress.

Integrated transcriptomic and metabolomic analysis reveals the effects of polyploidization on the lignin content and metabolic pathway in Eucalyptus

BACKGROUND

Lignin is a major restriction factor for the industrial production of biomass resources, such as pulp and bioenergy. Eucalyptus is one of the most important sources of pulp and bioenergy. After polyploidization, the lignin content of forest trees is generally reduced, which is considered a beneficial genetic improvement. However, the differences in the lignin content between triploid and diploid Eucalyptus and the underlying regulatory mechanism are still unclear.

RESULTS

We conducted a comprehensive analysis at the phenotypic, transcriptional and metabolite levels between Eucalyptus urophylla triploids and diploids to reveal the effects of polyploidization on the lignin content and lignin metabolic pathway. The results showed that the lignin content of Eucalyptus urophylla triploid stems was significantly lower than that of diploids. Lignin-related metabolites were differentially accumulated between triploids and diploids, among which coniferaldehyde, p-coumaryl alcohol, sinapaldehyde and coniferyl alcohol had significant positive correlations with lignin content, indicating that they might be primarily contributing metabolites. Most lignin biosynthetic genes were significantly downregulated, among which 11 genes were significantly positively correlated with the lignin content and above metabolites. Furthermore, we constructed a co-expression network between lignin biosynthetic genes and transcription factors based on weighted gene co-expression network analysis. The network identified some putative orthologues of secondary cell wall (SCW)-related transcription factors, among which MYB52, MYB42, NAC076, and LBD15 were significantly downregulated in Eucalyptus urophylla triploids. In addition, potential important transcription factors, including HSL1, BEE3, HHO3, and NAC046, also had high degrees of connectivity and high edge weights with lignin biosynthetic genes, indicating that they might also be involved in the variation of lignin accumulation between triploid and diploid Eucalyptus urophylla.

CONCLUSIONS

The results demonstrated that some lignin-related metabolites, lignin biosynthetic genes and transcription factors in Eucalyptus urophylla triploids may be relatively sensitive in response to the polyploidization effect, significantly changing their expression levels, which ultimately correlated with the varied lignin content. The analysis of the underlying formation mechanism could provide beneficial information for the development and utilization of polyploid biomass resources, which will be also valuable for genetic improvement in other bioenergy plants.

Exogenous Methyl Jasmonate Mediated MiRNA-mRNA Network Improves Heat Tolerance of Perennial Ryegrass

Heat stress can hinder the growth of perennial ryegrass (Lolium perenne L.). Methyl jasmonate (MeJA) applied exogenously can increase heat stress tolerance in plants; however, the regulatory mechanisms involved in heat tolerance mediated by MeJA are poorly understood in perennial ryegrass. Here, the microRNA (miRNA) expression profiles of perennial ryegrass were assessed to elucidate the regulatory pathways associated with heat tolerance induced by MeJA. Plants were subjected to four treatments, namely, control (CK), MeJA pre-treatment (T), heat stress treatment (H), and MeJA pre-treatment + heat stress (TH). According to the results, 102 miRNAs were up-regulated in all treatments, with 20, 27, and 33 miRNAs being up-regulated in the T, H, and TH treatment groups, respectively. The co-expression network analysis between the deferentially expressed miRNAs and their corresponding target genes showed that 20 miRNAs modulated 51 potential target genes. Notably, the miRNAs that targeted genes related to with regards to heat tolerance were driven by MeJA, and they were involved in four pathways: novel-m0258-5p mediated signal transduction, novel-m0350-5p mediated protein homeostasis, miR397-z, miR5658-z, and novel-m0008-5p involved in cell wall component, and miR1144-z and miR5185-z dominated chlorophyll degradation. Overall, the findings of this research paved the way for more research into the heat tolerance mechanism in perennial ryegrass and provided a theoretical foundation for developing cultivars with enhanced heat tolerance.

Comparative transcriptome profiling of susceptible and tolerant citrus species at early and late stage of infection by "Candidatus Liberibacter asiaticus"

Citrus Huanglongbing (HLB), caused by "Candidatus Liberibacter asiaticus" (CLas), is the most destructive disease threatening global citrus industry. Most commercial cultivars were susceptible to HLB, although some showed tolerant to HLB phenotypically. Identifying tolerant citrus genotypes and understanding the mechanism correlated with tolerance to HLB is essential for breeding citrus variety tolerance/resistance to HLB. In this study, the graft assay with CLas-infected bud were performed in four citrus genotypes, including Citrus reticulata Blanco, C. sinensis, C. limon, and C. maxima. HLB tolerance was observed in C. limon and C. maxima, while C. Blanco and C. sinensis were susceptible to HLB. The time-course transcriptomic analysis revealed a significant variation in differentially expressed genes (DEGs) related to HLB between susceptible and tolerant cultivar group at early and late infection stage. Functional analysis of DEGs indicated that the activation of genes involved in SA-mediated defense response, PTI, cell wall associated immunity, endochitinase, phenylpropanoid and alpha-linolenic/linoleic lipid metabolism played an important in the tolerance of C. limon and C. maxima to HLB at early infection stage. In addition, the overactive plant defense combined with the stronger antibacterial activity (antibacterial secondary and lipid metabolism) and the suppression of pectinesterase were contributed to the long-term tolerance to HLB in C. limon and C. maxima at late infection stage. Particularly, the activation of ROS scavenging genes (catalases and ascorbate peroxidases) could help to reduce HLB symptoms in tolerant cultivars. In contrast, the overexpression of genes involved in oxidative burst and ethylene metabolism, as well as the late inducing of defense related genes could lead to the early HLB symptom development in susceptible cultivars at early infection stage. The weak defense response and antibacterial secondary metabolism, and the induce of pectinesterase were responsible for sensitivity to HLB in C. reticulata Blanco and C. sinensis at late infection stage. This study provided new insights into the tolerance/sensitivity mechanism against HLB and valuable guidance for breeding of HLB-tolerant/resistant cultivars.

Systematic approaches to C-lignin engineering in Medicago truncatula

BACKGROUND

C-lignin is a homopolymer of caffeyl alcohol present in the seed coats of a variety of plant species including vanilla orchid, various cacti, and the ornamental plant Cleome hassleriana. Because of its unique chemical and physical properties, there is considerable interest in engineering C-lignin into the cell walls of bioenergy crops as a high-value co-product of bioprocessing. We have used information from a transcriptomic analysis of developing C. hassleriana seed coats to suggest strategies for engineering C-lignin in a heterologous system, using hairy roots of the model legume Medicago truncatula.

RESULTS

We systematically tested strategies for C-lignin engineering using a combination of gene overexpression and RNAi-mediated knockdown in the caffeic acid/5-hydroxy coniferaldehyde 3/5-O-methyltransferase (comt) mutant background, monitoring the outcomes by analysis of lignin composition and profiling of monolignol pathway metabolites. In all cases, C-lignin accumulation required strong down-regulation of caffeoyl CoA 3-O-methyltransferase (CCoAOMT) paired with loss of function of COMT. Overexpression of the Selaginella moellendorffii ferulate 5-hydroxylase (SmF5H) gene in comt mutant hairy roots resulted in lines that unexpectedly accumulated high levels of S-lignin.

CONCLUSION

C-Lignin accumulation of up to 15% of total lignin in lines with the greatest reduction in CCoAOMT expression required the strong down-regulation of both COMT and CCoAOMT, but did not require expression of a heterologous laccase, cinnamyl alcohol dehydrogenase (CAD) or cinnamoyl CoA reductase (CCR) with preference for 3,4-dihydroxy-substituted substrates in M. truncatula hairy roots. Cell wall fractionation studies suggested that the engineered C-units are not present in a heteropolymer with the bulk of the G-lignin.

Genome-Wide Identification and Expression Analysis of Dendrocalamus farinosus CCoAOMT Gene Family and the Role of DfCCoAOMT14 Involved in Lignin Synthesis

As the main component of plant cell walls, lignin can not only provide mechanical strength and physical defense for plants, but can also be an important indicator affecting the properties and quality of wood and bamboo. Dendrocalamus farinosus is an important economic bamboo species for both shoots and timber in southwest China, with the advantages of fast growth, high yield and slender fiber. Caffeoyl-coenzyme A-O-methyltransferase (CCoAOMT) is a key rate-limiting enzyme in the lignin biosynthesis pathway, but little is known about it in D. farinosus. Here, a total of 17 DfCCoAOMT genes were identified based on the D. farinosus whole genome. DfCCoAOMT1/14/15/16 were homologs of AtCCoAOMT1. DfCCoAOMT6/9/14/15/16 were highly expressed in stems of D. farinosus; this is consistent with the trend of lignin accumulation during bamboo shoot elongation, especially DfCCoAOMT14. The analysis of promoter cis-acting elements suggested that DfCCoAOMTs might be important for photosynthesis, ABA/MeJA responses, drought stress and lignin synthesis. We then confirmed that the expression levels of DfCCoAOMT2/5/6/8/9/14/15 were regulated by ABA/MeJA signaling. In addition, overexpression of DfCCoAOMT14 in transgenic plants significantly increased the lignin content, xylem thickness and drought resistance of plants. Our findings revealed that DfCCoAOMT14 can be a candidate gene that is involved in the drought response and lignin synthesis pathway in plants, which could contribute to the genetic improvement of many important traits in D. farinosus and other species.

In-silico identification and characterization of O-methyltransferase gene family in peanut (Arachis hypogaea L.) reveals their putative roles in development and stress tolerance

Cultivated peanut (Arachis hypogaea) is a leading protein and oil-providing crop and food source in many countries. At the same time, it is affected by a number of biotic and abiotic stresses. O-methyltransferases (OMTs) play important roles in secondary metabolism, biotic and abiotic stress tolerance. However, the OMT genes have not been comprehensively analyzed in peanut. In this study, we performed a genome-wide investigation of A. hypogaea OMT genes (AhOMTs). Gene structure, motifs distribution, phylogenetic history, genome collinearity and duplication of AhOMTs were studied in detail. Promoter cis-elements, protein-protein interactions, and micro-RNAs targeting AhOMTs were also predicted. We also comprehensively studied their expression in different tissues and under different stresses. We identified 116 OMT genes in the genome of cultivated peanut. Phylogenetically, AhOMTs were divided into three groups. Tandem and segmental duplication events played a role in the evolution of AhOMTs, and purifying selection pressure drove the duplication process. AhOMT promoters were enriched in several key cis-elements involved in growth and development, hormones, light, and defense-related activities. Micro-RNAs from 12 different families targeted 35 AhOMTs. GO enrichment analysis indicated that AhOMTs are highly enriched in transferase and catalytic activities, cellular metabolic and biosynthesis processes. Transcriptome datasets revealed that AhOMTs possessed varying expression levels in different tissues and under hormones, water, and temperature stress. Expression profiling based on qRT-PCR results also supported the transcriptome results. This study provides the theoretical basis for further work on the biological roles of AhOMT genes for developmental and stress responses.

Evaluation of Rumen Degradation Kinetics of Low-Lignin Alfalfa ‘Hi-Gest® 360’ in Saskatchewan Canada

The objective of this study was to determine rumen degradation kinetics of new low-lignin alfalfa (Medicago sativa L.) cv. Hi-Gest®360 (HiGest) in comparison with conventional alfalfa cv. AC Grazeland (Grazeland) in monoculture and binary mixtures at different maturity stages. Two cultivars of alfalfa (HiGest, and AC Grazeland) and their binary mixtures with hybrid bromegrass (HBG; cv. AC Success), grown in 2019 at two locations (Saskatoon and Lanigan), were cut at three maturity stages of alfalfa (1 = 10% bloom; 2 = 40% bloom; and 3 = 100% bloom). Rumen degradation characteristics, including rapidly degradable fraction (S), potentially degradable fraction (D), undegradable fraction (U), degradation rate (Kd), lag time (T0), and effective degradability (ED) of each component were determined using in situ technique and were analyzed by a first-order kinetic equation described by Ørskov and McDonald with lag time. Generally, in alfalfa monoculture, S or D were decreased and U was increased without affecting Kd and T0, resulting in decreased ED fraction with increasing stage of maturity. In binary mixtures, plant maturity stages have negligible effects on rumen degradation characteristics of CP. HiGest had higher effective degradability of DM (EDDM) as well as of NDF (EDNDF) than Grazeland. In conclusion, HiGest had greater DM and NDF rumen degradation potential relative to Grazeland. HiGest and Grazeland were different in DM and CP degradation patterns, with HiGest having higher EDDM and EDCP than Grazeland.

Transcriptome analysis identifies differentially expressed genes involved in lignin biosynthesis in barley

Lignin is an essential metabolite for plant growth but negatively affects the quality of forage barley. Genetic modification of quality traits to improve the forage digestibility requires an understanding of the molecular mechanism of lignin biosynthesis. RNA-Seq was used to quantify transcripts differentially expressed among leaf, stem and spike tissues from two barley genotypes. A total of 13,172 differentially expressed genes (DEGs) were identified, of which much more up-regulated DEGs were detected from the contrasting groups of leaf vs spike (L-S) and stem vs spike (S-S), and down-regulated DEGs were dominant in the group of stem vs leaf (S-L). 47 DEGs were successfully annotated to the monolignol pathway and six of them were candidate genes regulating the lignin biosynthesis. The qRT-PCR assay verified the expression profiles of the six candidate genes. Among them, four genes might positively regulate the lignin biosynthesis during forage barley development in terms of the consistency of their expression levels and changes of lignin content among the tissues, while the other two genes may have the reverse effects. These findings provide target genes for further investigations on molecular regulatory mechanisms of lignin biosynthesis and genetic resources for improvement of forage quality in barley molecular breeding programme.

Lignocellulose molecular assembly and deconstruction properties of lignin-altered rice mutants

Bioengineering approaches to modify lignin content and structure in plant cell walls have shown promise for facilitating biochemical conversions of lignocellulosic biomass into valuable chemicals. Despite numerous research efforts, however, the effect of altered lignin chemistry on the supramolecular assembly of lignocellulose and consequently its deconstruction in lignin-modified transgenic and mutant plants is not fully understood. In this study, we aimed to close this gap by analyzing lignin-modified rice (Oryza sativa L.) mutants deficient in 5-HYDROXYCONIFERALDEHYDE O-METHYLTRANSFERASE (CAldOMT) and CINNAMYL ALCOHOL DEHYDROGENASE (CAD). A set of rice mutants harboring knockout mutations in either or both OsCAldOMT1 and OsCAD2 was generated in part by genome editing and subjected to comparative cell wall chemical and supramolecular structure analyses. In line with the proposed functions of CAldOMT and CAD in grass lignin biosynthesis, OsCAldOMT1-deficient mutant lines produced altered lignins depleted of syringyl and tricin units and incorporating noncanonical 5-hydroxyguaiacyl units, whereas OsCAD2-deficient mutant lines produced lignins incorporating noncanonical hydroxycinnamaldehyde-derived units. All tested OsCAldOMT1- and OsCAD2-deficient mutants, especially OsCAldOMT1-deficient lines, displayed enhanced cell wall saccharification efficiency. Solid-state nuclear magnetic resonance (NMR) and X-ray diffraction analyses of rice cell walls revealed that both OsCAldOMT1- and OsCAD2 deficiencies contributed to the disruptions of the cellulose crystalline network. Further, OsCAldOMT1 deficiency contributed to the increase of the cellulose molecular mobility more prominently than OsCAD2 deficiency, resulting in apparently more loosened lignocellulose molecular assembly. Such alterations in cell wall chemical and supramolecular structures may in part account for the variations of saccharification performance of the OsCAldOMT1- and OsCAD2-deficient rice mutants.

Low nitrogen stress-induced transcriptome changes revealed the molecular response and tolerance characteristics in maintaining the C/N balance of sugar beet (Beta vulgaris L.)

Nitrogen (N) is an essential macronutrient for plants, acting as a common limiting factor for crop yield. The application of nitrogen fertilizer is related to the sustainable development of both crops and the environment. To further explore the molecular response of sugar beet under low nitrogen (LN) supply, transcriptome analysis was performed on the LN-tolerant germplasm '780016B/12 superior'. In total, 580 differentially expressed genes (DEGs) were identified in leaves, and 1,075 DEGs were identified in roots (log₂ ^|FC| ≥ 1; q value < 0.05). Gene Ontology (GO), protein-protein interaction (PPI), and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses clarified the role and relationship of DEGs under LN stress. Most of the downregulated DEGs were closely related to "photosynthesis" and the metabolism of "photosynthesis-antenna proteins", "carbon", "nitrogen", and "glutathione", while the upregulated DEGs were involved in flavonoid and phenylalanine biosynthesis. For example, GLUDB (glutamate dehydrogenase B) was identified as a key downregulated gene, linking carbon, nitrogen, and glutamate metabolism. Thus, low nitrogen-tolerant sugar beet reduced energy expenditure mainly by reducing the synthesis of energy-consuming amino acids, which in turn improved tolerance to low nitrogen stress. The glutathione metabolism biosynthesis pathway was promoted to quench reactive oxygen species (ROS) and protect cells from oxidative damage. The expression levels of nitrogen assimilation and amino acid transport genes, such as NRT2.5 (high-affinity nitrate transporter), NR (nitrate reductase [NADH]), NIR (ferredoxin-nitrite reductase), GS (glutamine synthetase leaf isozyme), GLUDB, GST (glutathione transferase) and GGT3 (glutathione hydrolase 3) at low nitrogen levels play a decisive role in nitrogen utilization and may affect the conversion of the carbon skeleton. DFRA (dihydroflavonol 4-reductase) in roots was negatively correlated with NIR in leaves (coefficient = -0.98, p < 0.05), suggesting that there may be corresponding remote regulation between "flavonoid biosynthesis" and "nitrogen metabolism" in roots and leaves. FBP (fructose 1,6-bisphosphatase) and PGK (phosphoglycerate kinase) were significantly positively correlated (p < 0.001) with Ci (intercellular CO₂ concentration). The reliability and reproducibility of the RNA-seq data were further confirmed by real-time fluorescence quantitative PCR (qRT-PCR) validation of 22 genes (R² = 0.98). This study reveals possible pivotal genes and metabolic pathways for sugar beet adaptation to nitrogen-deficient environments.

Caffeoyl-CoA 3-O-methyltransferase gene family in jute: Genome-wide identification, evolutionary progression and transcript profiling under different quandaries

Jute (Corchorus sp.), is a versatile, naturally occurring, biodegradable material that holds the promising possibility of diminishing the extensive use of plastic bags. One of the major components of the cell wall, lignin plays both positive and negative roles in fiber fineness and quality. Although it gives mechanical strength to plants, an excess amount of it is responsible for the diminution of fiber quality. Among various gene families involved in the lignin biosynthesis, Caffeoyl-CoA 3-O-methyltransferase (CCoAOMT) is the most significant and has remained mostly unexplored. In this study, an extensive in-silico characterization of the CCoAOMT gene family was carried out in two jute species (C. capsularis L. and C. olitoroius L.) by analyzing their structural, functional, molecular and evolutionary characteristics. A total of 6 CCoAOMT gene members were identified in each of the two species using published reference genomes. These two jute species showed high syntenic conservation and the identified CCoAOMT genes formed four clusters in the phylogenetic tree. Histochemical assay of lignin in both jute species could shed light on the deposition pattern in stems and how it changes in response to abiotic stresses. Furthermore, expression profiling using qPCR showed considerable alteration of CCoAOMT transcripts under various abiotic stresses and hormonal treatment. This study will lay a base for further analysis and exploration of target candidates for overexpression of gene silencing using modern biotechnological techniques to enhance the quality of this economically important fiber crop.

Comparative Transcriptome Analysis of MeJA Responsive Enzymes Involved in Phillyrin Biosynthesis of Forsythia suspensa

Forsythia suspensa (Thunb.) has been widely used in traditional medicines in Asia. According to the 2020 edition of Chinese Pharmacopoeia, phillyrin is the main active ingredient in F. suspensa, which is effective in clearing heat, reducing swelling, and dispersing nodules. F. suspensa leaf is a non-toxic substance and it can be used to make a health tea. Here, we combine elicitors and transcriptomics to investigate the inducible biosynthesis of the phillyrin from the F. suspensa. After the fruits and leaves of F. suspensa were treated with different concentrations of methyl jasmonate (MeJA), the content of phillyrin in the fruits reached a peak at 200 µM MeJA for 12 h, but which was decreased in leaves. To analyze the differences in key enzyme genes involved in the phillyrin biosynthesis, we sequenced the transcriptome of F. suspensa leaves and fruits treated with 200 µM MeJA for 12 h. We hypothesized that nine genes related to coniferin synthesis including: F. suspensa UDP-glycosyltransferase (FsUGT); F. suspensa 4-coumarate coenzyme CoA ligase (Fs4CL); and F. suspensa Caffeoyl-CoA O-methyltransferase (FsCCoAOMT) etc. The qRT-PCR analysis of genes related to phillyrin biosynthesis was consistent with RNA-seq analysis. We also investigated the dynamic changes of genes in F. suspensa leaves and fruits at different time points after 200 µM MeJA treatment, which laid the foundation for further study of the molecular mechanisms regulating the biosynthesis of phillyrin.

The physiological and transcriptomic study of secondary growth in Neolamarckia cadamba stimulated by the ethylene precursor ACC

Though many biological roles of ethylene have been investigated intensively, the molecular mechanism of ethylene's action in woody plants remains unclear. In this study, we investigated the effects of exogenous 1-aminocyclopropane-1-carboxylic acid (ACC), the precursor of ethylene, on the growth of Neolamarckia cadamba seedlings, a fast-growing tropical tree. After 14 days of ACC treatment, the plants showed a reduced physiological morphology while stem diameter increased; however, this did not occur after the addition of 1-MCP. Meanwhile, the lignin content of N. cadamba also increased. Transcriptome analysis revealed that the expression of the ethylene biosynthesis and signaling genes ACC oxidase (ACO) and ethylene insensitive 3 (EIN3) were up-regulated mainly at the 6th hour and the 3rd day of the ACC treatment, respectively. The transcription levels of transcription factors, mainly in the basic helix-loop-helix (bHLH), ethylene response factor (ERF), WRKY and v-myb avian myeloblastosis viral oncogene homolog (MYB) families, involved in the ethylene signaling and secondary growth also increased significantly. Furthermore, in accordance to the increased lignification of the stem, the transcriptional level of key enzymes in the phenylalanine pathway were elevated after the ACC treatment. Our results revealed the physiological and molecular mechanisms underlying the secondary growth stimulated by exogenous ACC treatment on N. cadamba seedlings.

Comparative Transcriptome and Anatomic Characteristics of Stems in Two Alfalfa Genotypes

Stems are more important to forage quality than leaves in alfalfa. To understand lignin formation at different stages in alfalfa, lignin distribution, anatomical characteristics and transcriptome profile were employed using two alfalfa cultivars. The results showed that the in vitro true digestibility (IVTD) of stems in WL168 was significantly higher than that of Zhungeer, along with the significantly lower neutral detergent fiber (NDF), acid detergent fiber (ADF) and lignin contents. In addition, Zhungeer exhibited increased staining of the xylem areas in the stems of different developmental stages compared to WL168. Interestingly, the stems of WL168 appeared intracellular space from the stage 3, while Zhungeer did not. The comparative transcriptome analysis showed that a total of 1993 genes were differentially expressed in the stem between the cultivars, with a higher number of expressed genes in the stage 4. Of the differentially expressed genes, starch and sucrose metabolism as well as phenylpropanoid biosynthesis pathways were the most significantly enriched pathways. Furthermore, expression of genes involved in lignin biosynthesis such as PAL, 4CL, HCT, CAD, COMT and POD coincides with the anatomic characteristics and lignin accumulation. These results may help elucidate the regulatory mechanisms of lignin biosynthesis and improve forage quality in alfalfa.

Expression of a Hydroxycinnamoyl-CoA Shikimate/Quinate Hydroxycinnamoyl Transferase 4 Gene from Zoysia japonica (ZjHCT4) Causes Excessive Elongation and Lignin Composition Changes in Agrostis stolonifera

Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT) is considered to be an essential enzyme for regulating the biosynthesis and composition of lignin. To investigate the properties and function of ZjHCT4, the ZjHCT4 gene was cloned from Zoysia japonica with a completed coding sequence of 1284-bp in length, encoding 428 amino acids. The ZjHCT4 gene promoter has several methyl jasmonate (MeJA) response elements. According to analysis of expression patterns, it was up-regulated by MeJA, GA₃ (Gibberellin), and SA (Salicylic acid), and down-regulated by ABA (Abscisic acid). Ectopic ZjHCT4 expression in creeping bentgrass causes excessive plant elongation. In addition, the content of G-lingnin and H-lingnin fell in transgenic plants, whereas the level of S-lingnin increased, resulting in a considerable rise in the S/G unit ratio. Analysis of the expression levels of lignin-related genes revealed that the ectopic expression of ZjHCT4 altered the expression levels of a number of genes involved in the lignin synthesis pathway. Simultaneously, MeJA, SA, GA₃, IAA, BR (Brassinosteroid), and other hormones were dramatically enhanced in transgenic plants relative to control plants, whereas ABA concentration was significantly decreased. Expression of ZjHCT4 impacted lignin composition and plant growth via altering the phenylpropionic acid metabolic pathway and hormone response, as revealed by transcriptome analysis. HCTs may influence plant lignin composition and plant development by altering hormone content. These findings contributed to a deeper comprehension of the lignin synthesis pathway and set the stage for further investigation and application of the HCTs gene.

Transcriptome of the coralline alga Calliarthron tuberculosum (Corallinales, Rhodophyta) reveals convergent evolution of a partial lignin biosynthesis pathway

The discovery of lignins in the coralline red alga Calliarthron tuberculosum raised new questions about the deep evolution of lignin biosynthesis. Here we present the transcriptome of C. tuberculosum supported with newly generated genomic data to identify gene candidates from the monolignol biosynthetic pathway using a combination of sequence similarity-based methods. We identified candidates in the monolignol biosynthesis pathway for the genes 4CL, CCR, CAD, CCoAOMT, and CSE but did not identify candidates for PAL, CYP450 (F5H, C3H, C4H), HCT, and COMT. In gene tree analysis, we present evidence that these gene candidates evolved independently from their land plant counterparts, suggesting convergent evolution of a complex multistep lignin biosynthetic pathway in this red algal lineage. Additionally, we provide tools to extract metabolic pathways and genes from the newly generated transcriptomic and genomic datasets. Using these methods, we extracted genes related to sucrose metabolism and calcification. Ultimately, this transcriptome will provide a foundation for further genetic and experimental studies of calcifying red algae.

The Larix kaempferi genome reveals new insights into wood properties

Here, through single-molecule real-time sequencing, we present a high-quality genome sequence of the Japanese larch (Larix kaempferi), a conifer species with great value for wood production and ecological afforestation. The assembled genome is 10.97 Gb in size, harboring 45,828 protein-coding genes. Of the genome, 66.8% consists of repeat sequences, of which long terminal repeat retrotransposons are dominant and make up 69.86%. We find that tandem duplications have been responsible for the expansion of genes involved in transcriptional regulation and stress responses, unveiling their crucial roles in adaptive evolution. Population transcriptome analysis reveals that lignin content in L. kaempferi is mainly determined by the process of monolignol polymerization. The expression values of six genes (LkCOMT7, LkCOMT8, LkLAC23, LkLAC102, LkPRX148, and LkPRX166) have significantly positive correlations with lignin content. These results indicated that the increased expression of these six genes might be responsible for the high lignin content of the larches' wood. Overall, this study provides new genome resources for investigating the evolution and biological function of conifer trees, and also offers new insights into wood properties of larches.

Transcriptome and metabolite profiling to identify genes associated with rhizome lignification and the function of ZoCSE in ginger (Zingiber officinale)

Ginger (Zingiber officinale Roscoe) is an important spice crop in China, and fresh ginger rhizomes are consumed as vegetable in Sichuan and Chongqing. However, tissue lignification accelerates with rhizome maturation, resulting in the loss of edible quality. To understand the molecular mechanisms of texture modification during rhizome development, we investigated lignin accumulation patterns and identified the key genes associated with lignin biosynthesis using gas chromatography-mass spectrometry (GC-MS), liquid chromatography-tandem mass spectrometry (LC-MS/MS) and RNA-sequencing (RNA-Seq). Results showed that the contents of total lignin and its precursors exhibited notable declines with tissue maturation. However, the lignin composition was remarkably modified and syringyl lignin was deposited in mature rhizomes, leading to ginger lignification. Transcriptome analysis displayed 32 lignin biosynthetic genes were dramatically downregulated with rhizome development, including caffeoylshikimate esterase (CSE ), 4-coumarate-CoA ligase , laccase , cinnamoyl-CoA reductase , cinnamyl-alcohol dehydrogenase , peroxidase and caffeic acid 3-O-methyltransferase , indicating that lignin reduction might be attributed to deficiency in intermediates or the downregulation of key biosynthetic enzymes. Furthermore, overexpressing ZoCSE in Nicotiana benthamiana L. enhanced the total lignin content, suggesting its fundamental role in lignin biosynthesis. RNA-Seq also identified candidate lignin production regulators, including hormone-related genes and NAC/MYB transcription factors (ZoNAC1 , ZoNAC4 , ZoMYB14 and ZoMYB17 ). This result provides a molecular basis for lignin accumulation in ginger.

Down-regulation of PvSAMS impairs S-adenosyl-L-methionine and lignin biosynthesis, and improves cell wall digestibility in switchgrass

S-adenosyl- l-methionine (SAM) is the methyl donor involved in the biosynthesis of guaiacyl (G) and syringyl (S) lignins in vascular plants. SAM is synthesized from methionine through the catalysis of the enzyme S-adenosylmethionine synthase (SAMS). However, the detailed function of SAMS in lignin biosynthesis has not been widely investigated in plants, particularly in monocot species. In this study, we identified PvSAMS genes from switchgrass (Panicum virgatum L.), an important dual-purpose fodder and biofuel crop, and generated numerous transgenic switchgrass lines through PvSAMS RNA interference technology. Down-regulation of PvSAMS reduced the contents of SAM, G-lignins, and S-lignins in the transgenic switchgrass. The methionine and glucoside derivatives of caffeoyl alcohol were found to accumulate in the transgenic plants. Moreover, down-regulation of PvSAMS in switchgrass resulted in brownish stems associated with reduced lignin content and improved cell wall digestibility. Furthermore, transcriptomic analysis revealed that most sulfur deficiency-responsive genes were differentially expressed in the transgenic switchgrass, leading to a significant increase in total sulfur content; thus implying an important role of SAMS in the methionine cycle, lignin biosynthesis, and sulfur assimilation. Taken together, our results suggest that SAMS is a valuable target in lignin manipulation, and that manipulation of PvSAMS can simultaneously regulate the biosynthesis of SAM and methylated monolignols in switchgrass.

Molecular Mechanisms of Phenylpropane-Synthesis-Related Genes Regulating the Shoot Blight Resistance of Bambusa pervariabilis × Dendrocalamopsis grandis

Bambusa pervariabilis × Dendrocalamopsis grandis shoot blight caused by Arthrinium phaeospermum is a fungal disease that has affected a large area in China in recent years. However, it is not clear which genes are responsible for the disease resistance of B. pervariabilis × D. grandis. Based on the analysis of transcriptome and proteome data, two genes, CCoAOMT2 and CAD5, which may be involved in disease resistance, were screened. Two gene expression-interfering varieties, COF RNAi and CAD RNAi were successfully obtained using RNAi technology. Quantitative real-time fluorescence (qRT-PCR) results showed that CCoAOMT2 gene, CAD5 gene and seven related genes expression was down-regulated in the transformed varieties. After inoculating pathogen spore suspension, the incidence and disease index of cof-RNAi and cad-RNAi transformed plants increased significantly. At the same time, it was found that the content of total lignin and flavonoids in the two transformed varieties were significantly lower than that of the wild-type. The subcellular localization results showed that both CCoAOMT2 and CAD5 were localized in the nucleus and cytoplasm. The above results confirm that the CCoAOMT2 and CAD5 genes are involved in the resistance of B. pervariabilis × D.grandis to shoot blight through regulating the synthesis of lignin and flavonoids.

The Genetic Architecture of Nitrogen Use Efficiency in Switchgrass (Panicum virgatum L.)

Switchgrass (Panicum virgatum L.) has immense potential as a bioenergy crop with the aim of producing biofuel as an end goal. Nitrogen (N)-related sustainability traits, such as nitrogen use efficiency (NUE) and nitrogen remobilization efficiency (NRE), are important factors affecting switchgrass quality and productivity. Hence, it is imperative to develop nitrogen use-efficient switchgrass accessions by exploring the genetic basis of NUE in switchgrass. For that, we used 331 diverse field-grown switchgrass accessions planted under low and moderate N fertility treatments. We performed a genome wide association study (GWAS) in a holistic manner where we not only considered NUE as a single trait but also used its related phenotypic traits, such as total dry biomass at low N and moderate N, and nitrogen use index, such as NRE. We have evaluated the phenotypic characterization of the NUE and the related traits, highlighted their relationship using correlation analysis, and identified the top ten nitrogen use-efficient switchgrass accessions. Our GWAS analysis identified 19 unique single nucleotide polymorphisms (SNPs) and 32 candidate genes. Two promising GWAS candidate genes, caffeoyl-CoA O-methyltransferase (CCoAOMT) and alfin-like 6 (AL6), were further supported by linkage disequilibrium (LD) analysis. Finally, we discussed the potential role of nitrogen in modulating the expression of these two genes. Our findings have opened avenues for the development of improved nitrogen use-efficient switchgrass lines.

Structure basis of the caffeic acid O-methyltransferase from Ligusiticum chuanxiong to understand its selective mechanism

Caffeic acid O-methyltransferase from Ligusticum chuanxiong (LcCOMT) showed strict regiospecificity despite a relative degree of preference. Compared with caffeic acid, methyl caffeate was the preferential substrate by its low Km and high Kcat. In this study, we obtained the SAM binary (1.80 Å) and SAH binary (1.95 Å) complex LcCOMT crystal structures, and established the ternary complex structure with methyl caffeate by molecular docking. The active site of LcCOMT included phenolic substrate pocket, SAM/SAH ligand pocket and conserved catalytic residues as well. The regiospecificity of LcCOMT that permitted only 3-hydroxyl group to be methylated arise from the interactions between the active site and the phenyl ring. However, the propanoid tail governed the relative preference of LcCOMT. The ester group in methyl caffeate stabilized the anionic intermediate caused by His268-Asp269 pair, whereas caffeic acid was unable to stabilize the anionic intermediate due to the adjacent carboxylate anion in the propanoid tail. Ser183 residue formed an additional hydrogen bond with SAH and its role was identified by S183A mutation. Ile318 residue might be a potential site for determination of substrate preference, and its mutation led to the change of tertiary conformation. The results supported the selective mechanism of LcCOMT.

Proteomic differences in apple spur buds from high and non-cropping trees during floral initiation

The cropping behavior of biennial apple (Malus ×domestica Borkh.) cultivars is irregular and often follows a biennial bearing pattern with 'On' years (high crop load and inhibited floral bud formation) followed by 'Off' years (little crop load and a promoted formation of floral buds). To study proteomic differences between floral and vegetative buds, trees of the strongly alternating cultivar 'Fuji' and the regular bearing cultivar 'Gala' were either completely thinned or not thinned at full bloom to establish two cropping treatments with no ('Off') or a high ('On') crop load, respectively. Student's t-Tests indicated significant differences of protein profiles in buds from 2-year old spurs from both treatments at each sampling date. Abundance patterns of protein clusters coincided with the onset of floral bud initiation and were most noticeable in buds from 'On' trees with a decreased abundance of key enzymes of the phenylpropanoid and flavonoid pathways and an increased abundance of histone deacetylase and ferritins. Furthermore, an increased abundance of proteins involved in histone and DNA methylation was found in the buds from 'Off' trees. This study presents the first large-scale, label-free proteomic profiling of floral and vegetative apple buds during the period of floral bud initiation. SIGNIFICANCE: Although several studies exist that address the complex developmental processes associated with the formation of floral buds in apple (Malus ×domestica Borkh.) at transcriptomic level, no data is available for explaining the difference between floral and vegetative buds or biennial and regular bearing cultivars on a proteomic level. This study presents the first large-scale, label-free proteomic profiling of floral and vegetative apple buds from the two cultivars 'Fuji' and 'Royal Gala' during the period of floral bud initiation and renders possible the development of suitable biomarkers for biennial bearing in apple.

Quantitative Trait Loci and Transcriptome Analysis Reveal Genetic Basis of Fiber Quality Traits in CCRI70 RIL Population of Gossypium hirsutum

Upland cotton (Gossypium hirsutum) is widely planted around the world for its natural fiber, and producing high-quality fiber is essential for the textile industry. CCRI70 is a hybrid cotton plant harboring superior yield and fiber quality, whose recombinant inbred line (RIL) population was developed from two upland cotton varieties (sGK156 and 901-001) and were used here to investigate the source of high-quality related alleles. Based on the material of the whole population, a high-density genetic map was constructed using specific locus-amplified fragment sequencing (SLAF-seq). It contained 24,425 single nucleotide polymorphism (SNP) markers, spanning a distance of 4,850.47 centimorgans (cM) over 26 chromosomes with an average marker interval of 0.20 cM. In evaluating three fiber quality traits in nine environments to detect multiple environments stable quantitative trait loci (QTLs), we found 289 QTLs, of which 36 of them were stable QTLs and 18 were novel. Based on the transcriptome analysis for two parents and two RILs, 24,941 unique differentially expressed genes (DEGs) were identified, 473 of which were promising genes. For the fiber strength (FS) QTLs, 320 DEGs were identified, suggesting that pectin synthesis, phenylpropanoid biosynthesis, and plant hormone signaling pathways could influence FS, and several transcription factors may regulate fiber development, such as GAE6, C4H, OMT1, AFR18, EIN3, bZIP44, and GAI. Notably, the marker D13_56413025 in qFS-chr18-4 provides a potential basis for enhancing fiber quality of upland cotton via marker-assisted breeding and gene cloning of important fiber quality traits.

Protoplast isolation and transcriptome analysis of developing xylem in Pinus massoniana (Pinaceae)

BACKGROUND

With active physiological and biochemical activities, tissue-specific protoplasts from cambial derivatives, could serve as a specific source for information on xylogenesis for softwood species resistant to stable genetic transformation and lacking available mutants.

METHODS AND RESULTS

In this study, protoplasts were isolated from developing xylem of the Chinese red pine, Pinus massoniana, by enzymolysis. High-quality RNAs were extracted from developing xylem and their protoplasts for constructing transcriptome libraries. Using Illumina HiSeq 2500 PE150 platform, a total of 362,328,426 clean paired-end reads (54.35G) were generated from multiple cDNA libraries and assembled into 146,422 unigenes. The transcriptome data were further analysed to identify 1567 differentially expressed genes (DEGs) between the isolated protoplasts and developing xylem of P. massoniana (Masson pine), 1126 DEGs were upregulated in protoplasts relative to developing xylem cells and 441 were downregulated. Most of the differentially expressed genes in biological process terms are related to plant response, which may be due to the response to cell wall removal. Further, the expression pattern of 71 unigenes involved in lignin biosynthesis was verified by RNA-seq.

CONCLUSIONS

This study is the first to report the transcriptome profiles of the developing xylem and its protoplasts of coniferous trees, which provide a new perspective and valuable resource for tracking transcriptional regulatory events in wood formation of Masson pine.

The regulation of cell wall lignification and lignin biosynthesis during pigmentation of winter jujube

Fruit lignification is due to lignin deposition in the cell wall during cell development. However, there are few studies on the regulation of cell wall lignification and lignin biosynthesis during fruit pigmentation. In this study, we investigated the regulation of cell wall lignification and lignin biosynthesis during pigmentation of winter jujube. The cellulose content decreased, while the lignin content increased in the winter jujube pericarp during pigmentation. Safranin O-fast green staining showed that the cellulose content was higher in the cell wall of winter jujube prior to pigmentation, whereas the lignin in the cell wall increased after pigmentation. The thickness of the epidermal cells decreased with pericarp pigmentation. A combined metabolomics and transcriptomics analysis showed that guaiacyl-syringyl (G-S) lignin was the main lignin type in the pericarp of winter jujube, and F5H (LOC107424406) and CCR (LOC107420974) were preliminarily identified as the key genes modulating lignin biosynthesis in winter jujube. Seventeen MYB and six NAC transcription factors (TFs) with potential regulation of lignin biosynthesis were screened out based on phylogenetic analysis. Three MYB and two NAC TFs were selected as candidate genes and further studied in detail. Arabidopsis ectopic expression and winter jujube pericarp injection of the candidate genes indicated that the MYB activator (LOC107425254) and the MYB repressor (LOC107415078) control lignin biosynthesis by regulating CCR and F5H, while the NAC (LOC107435239) TF promotes F5H expression and positively regulates lignin biosynthesis. These findings revealed the lignin biosynthetic pathway and associated genes during pigmentation of winter jujube pericarp and provide a basis for further research on lignin regulation.

Melatonin synthesis genes N-acetylserotonin methyltransferases evolved into caffeic acid O-methyltransferases and both assisted in plant terrestrialization

Terrestrialization is one of the most momentous events in the history of plant life, which leads to the subsequent evolution of plant diversity. The transition species, in this process, had to acquire a range of adaptive mechanisms to cope with the harsh features of terrestrial environments compared to that of aquatic habitat. As an ancient antioxidant, a leading regulator of ROS signaling or homeostasis, and a presumed plant master regulator, melatonin likely assisted plants transition to land and their adaption to terrestrial ecosystems. N-acetylserotonin methyltransferases (ASMT) and caffeic acid O-methyltransferases (COMT), both in the O-methyltransferase (OMT) family, catalyze the core O-methylation reaction in melatonin biosynthesis. How these two enzymes with close relevance evolved in plant evolutionary history and whether they participated in plant terrestrialization remains unknown. Using combined phylogenetic evidence and protein structure analysis, it is revealed that COMT likely evolved from ASMT by gene duplication and subsequent divergence. Newly emergent COMT gained a significantly higher ASMT activity to produce greater amounts of melatonin for immobile plants to acclimate to the stressful land environments after evolving from the more environmentally-stable aquatic conditions. The COMT genes possess more conserved substrate-binding sites at the amino acid level and more open protein conformation compared to ASMT, and getting a new function to catalyze the lignin biosynthesis. This development directly contributed to the dominance of vascular plants among the Earth's flora and prompted plant colonization of land. Thus, ASMT, together with its descendant COMT, might play key roles in plant transition to land. The current study provides new insights into plant terrestrialization with gene duplication contributing to this process along with well-known horizontal gene transfer.

Transcriptomic and Metabolomic Analyses Provide Insights Into an Aberrant Tissue of Tea Plant (Camellia sinensis)

Tea plant (Camellia sinensis (L.) O. Kuntze) is one of the most important economic crops with multiple mutants. Recently, we found a special tea germplasm that has an aberrant tissue on its branches. To figure out whether this aberrant tissue is associated with floral bud (FB) or dormant bud (DB), we performed tissue section, transcriptome sequencing, and metabolomic analysis of these tissues. Longitudinal sections indicated the aberrant tissue internal structure was more like a special bud (SB), but was similar to that of DB. Transcriptome data analysis showed that the number of heterozygous and homozygous SNPs was significantly different in the aberrant tissue compared with FB and DB. Further, by aligning the unmapped sequences of the aberrant tissue to the Non-Redundant Protein Sequences (NR) database, we observed that 36.13% of unmapped sequences were insect sequences, which suggested that the aberrant tissue might be a variation of dormant bud tissue influenced by the interaction of tea plants and insects or pathogens. Metabolomic analysis showed that the differentially expressed metabolites (DEMs) between the aberrant tissue and DB were significantly enriched in the metabolic pathways of biosynthesis of plant hormones and biosynthesis of phenylpropanoids. Subsequently, we analyzed the differentially expressed genes (DEGs) in the above mentioned two tissues, and the results indicated that photosynthetic capacity in the aberrant tissue was reduced, whereas the ethylene, salicylic acid and jasmonic acid signaling pathways were activated. We speculated that exogenous infection induced programmed cell death (PCD) and increased the lignin content in dormant buds of tea plants, leading to the formation of this aberrant tissue. This study advanced our understanding of the interaction between plants and insects or pathogens, providing important clues about biotic stress factors and key genes that lead to mutations and formation of the aberrant tissue.