Integrated Analysis of the Transcriptome and Metabolome Reveals Genes Involved in Terpenoid and Flavonoid Biosynthesis in the Loblolly Pine (Pinus taeda L.)

Multi-omics explores the potential regulatory role of acetylation modification in flavonoid biosynthesis of Ginkgo biloba

Flavonoids are crucial medicinal active ingredients in Ginkgo biloba L. However, the effect of protein post-translational modifications on flavonoid biosynthesis remains poorly explored. Lysine acetylation, a reversible post-translational modification, plays a crucial role in metabolic regulation. This study aims to investigate the potential role of acetylation in G. biloba flavonoid biosynthesis. Through comprehensive analysis of transcriptomes, metabolomes, proteomes and acetylated proteins in different tissues, a total of 11,788 lysine acetylation sites were identified on 4324 acetylated proteins, including 89 acetylation sites on 23 proteins. Additionally, 128 types of differentially accumulated flavonoids were identified among tissues, and a dataset of differentially expressed genes related to the flavonoid biosynthesis pathway was constructed. Twelve (CHI, C3H1, ANR, DFR, CCoAOMT1, F3H1, F3H2, CCoAOMT2, C3H2, HCT, F3'5'H and FG2) acetylated proteins that might be involved in flavonoid biosynthesis were identified. Specifically, we found that the modification levels of CCoAOMT1 and F3'5'H sites correlated with the catalytic production of homoeriodictyol and dihydromyricetin, respectively. Inhibitors of lysine deacetylase (trichostatin A) impacted total flavonoid content in different tissues and increased flavonoid levels in G. biloba roots. Treatment with trichostatin A revealed that expression levels of GbF3'5'H and GbCCoAOMT1 in stems and leaves aligned with total flavonoid content variations, while in roots, expression levels of GbC3H2 and GbFG2 corresponded to total flavonoid content changes. Collectively, these findings reveal for the first time the important role of acetylation in flavonoid biosynthesis.

Genome-wide identification and characterization of the AP2/ERF gene family in loblolly pine (Pinus taeda L.)

The loblolly pine (Pinus taeda L.) is one of the most profitable forest species worldwide owing to its quick growth, high wood yields, and strong adaptability. The AP2/ERF gene family plays a widespread role in the physiological processes of plant defense responses and the biosynthesis of metabolites. Nevertheless, there are no reports on this gene family in loblolly pine (P. taeda). In this study, a total of 303 members of the AP2/ERF gene family were identified. Through multiple sequence alignment and phylogenetic analysis, they were classified into four subfamilies, including AP2 (34), RAV (17), ERF (251), and Soloist (1). An analysis of the conservation domains, conserved motifs, and gene structure revealed that every PtAP2/ERF transcription factor (TF) had at least one AP2 domain. While evolutionary conservation was displayed within the same subfamilies, the distribution of conserved domains, conserved motifs, and gene architectures varied between subfamilies. Cis-element analysis revealed abundant light-responsive elements, phytohormone-responsive elements, and stress-responsive elements in the promoter of the PtAP2/ERF genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of potential target genes showed that the AP2/ERF gene family might play a critical role in plant growth and development, the response to environmental stresses, and metabolite biosynthesis. Utilizing quantitative real-time PCR (qRT-PCR), we examined the expression patterns of 10 randomly selected genes from Group IX after 6 h of treatments with mechanical injury, ethephon (Eth), and methyl jasmonate (MeJA). The AP2/ERF gene family in the loblolly pine was systematically analyzed for the first time in this study, offering a theoretical basis for exploring the functions and applications of AP2/ERF genes.

The Integration of the Metabolome and Transcriptome for Dendrobium nobile Lindl. in Response to Methyl Jasmonate

Dendrobium nobile Lindl., as an endangered medicinal plant within the genus Dendrobium, is widely distributed in southwestern China and has important ecological and economic value. There are a variety of metabolites with pharmacological activity in D. nobile. The alkaloids and polysaccharides contained within D. nobile are very important active components, which mainly have antiviral, anti-tumor, and immunity improvement effects. However, the changes in the compounds and functional genes of D. nobile induced by methyl jasmonate (MeJA) are not clearly understood. In this study, the metabolome and transcriptome of D. nobile were analyzed after exposure to MeJA. A total of 377 differential metabolites were obtained through data analysis, of which 15 were related to polysaccharide pathways and 35 were related to terpenoids and alkaloids pathways. Additionally, the transcriptome sequencing results identified 3256 differentially expressed genes that were discovered in 11 groups. Compared with the control group, 1346 unigenes were differentially expressed in the samples treated with MeJA for 14 days (TF14). Moreover, the expression levels of differentially expressed genes were also significant at different growth and development stages. According to GO and KEGG annotations, 189 and 99 candidate genes were identified as being involved in terpenoid biosynthesis and polysaccharide biosynthesis, respectively. In addition, the co-expression analysis indicated that 238 and 313 transcription factors (TFs) may contribute to the regulation of terpenoid and polysaccharide biosynthesis, respectively. Through a heat map analysis, fourteen terpenoid synthetase genes, twenty-three cytochrome P450 oxidase genes, eight methyltransferase genes, and six aminotransferase genes were identified that may be related to dendrobine biosynthesis. Among them, one sesquiterpene synthase gene was found to be highly expressed after the treatment with MeJA and was positively correlated with the content of dendrobine. This study provides important and valuable metabolomics and transcriptomic information for the further understanding of D. nobile at the metabolic and molecular levels and provides candidate genes and possible intermediate compounds for the dendrobine biosynthesis pathway, which lays a certain foundation for further research on and application of Dendrobium.

Phenolic Acids and Flavonoids Play Important Roles in Flower Bud Differentiation in Mikania micrantha: Transcriptomics and Metabolomics

Mikania micrantha is a highly invasive vine, and its ability to sexually reproduce is a major obstacle to its eradication. The long-distance dissemination of M. micrantha depends on the distribution of seeds; therefore, inhibiting M. micrantha flowering and seed production is an effective control strategy. The number of blooms of M. micrantha differs at different altitudes (200, 900, and 1300 m). In this study, we used a combination of metabolomics and transcriptomics methods to study the patterns of metabolite accumulation in the flower buds of M. micrantha. Using LC-MS/MS, 658 metabolites were found in the flower buds of M. micrantha at three different altitudes (200, 900, and 1300 m). Flavonoids and phenolic acids were found to be the main differential metabolites, and their concentrations were lower at 900 m than at 200 m and 1300 m, with the concentrations of benzoic acid, ferulic acid, and caffeic acid being the lowest. The biosynthesis pathways for flavonoids and phenolic compounds were significantly enriched for differentially expressed genes (DEGs), according to the results of transcriptome analysis. The production of flavonoid and phenolic acids was strongly linked with the expressions of phenylalanine ammonia-lyase (PAL), caffeoyl-CoA O-methyltransferase (COMT), and 4-coumarate-CoA ligase (4CL), according to the results of the combined transcriptome and metabolome analysis. These genes' roles in the regulation of distinct phenolic acids and flavonoids during M. micrantha bud differentiation are still unknown. This study adds to our understanding of how phenolic acids and flavonoids are regulated in M. micrantha flower buds at various altitudes and identifies regulatory networks that may be involved in this phenomenon, offering a new approach for the prevention and management of M. micrantha.

Temporospatial pattern of flavonoid metabolites and potential regulatory pathway of PbMYB211-coordinated kaempferol-3-O-rhamnoside biosynthesis in Phoebe bournei

Phoebe is a well-known timber tree species that contains abundant metabolites characterized by flavonoids that are widely used in the pharmaceutical industry. Nevertheless, temporospatial flavonoid metabolism variations substantially impact the Phoebe industry. Thus, a metabolomics analysis was carried out and identified 465 metabolites (102 flavonoids) in P. bournei, revealing distinct distribution patterns among five studied organs, and most of the flavonoids were dominant in the leaves. Furthermore, three kaempferol glycoside derivatives were significantly accumulated in the leaves and showed higher contents in young leaves than in mature leaves and differences between spring and autumn. For instance, greater accumulation of kaempferol-7-O-rhamnoside was detected in spring, whereas higher contents of kaempferol-3-O-arabinofuranoside and kaempferol-3-O-rhamnoside were found in autumn. Integrated metabolomics and transcriptomics identified 20 transcription factors (TFs) and 12 structural genes that participate in kaempferol derivative synthesis and elucidated a potential regulatory mechanism in P. bournei. Of the identified genes, PbMYB211 might contribute significantly to the kaempferol-3-O-rhamnoside content by regulating the target structural gene PbUGT139, as revealed by transient overexpression analysis. Overall, this study illuminated the temporospatial accumulation of flavonoids among different organs, seasons, and developmental stages in P. bournei and elucidated a potential regulatory pathway of kaempferol-3-O-rhamnoside. The results provide important insights into harvest techniques and a theoretical basis for the comprehensive utilization of P. bournei.

Nervonic acid improves liver inflammation in a mouse model of Parkinson's disease by inhibiting proinflammatory signaling pathways and regulating metabolic pathways

BACKGROUND

Nervonic acid (NA) - a type of bioactive fatty acid that is found in natural sources - can inhibit inflammatory reactions and regulate immune system balance. Therefore, the use of NA for the treatment of neurodegenerative diseases has received considerable attention. Our previous study found that NA inhibited inflammatory responses in the brain of Parkinson's disease (PD) mouse models. In addition to the brain, PD is also associated with visceral organ dysfunction, especially impaired liver function. Thus, studying the role of NA in PD-mediated inflammation of the liver is particularly important.

METHODS

A combined transcriptome and metabolomic approach was utilized to investigate the anti-inflammatory effects of NA on the liver of PD mice. Inflammatory signaling molecules and metabolic pathway-related genes were examined in the liver using real-time PCR and western blotting.

RESULTS

Liver transcriptome analysis revealed that NA exerted anti-inflammatory effects by controlling several pro-inflammatory signaling pathways, such as the down-regulation of the tumor necrosis factor and nuclear factor kappa B signaling pathways, both of which were essential in the development of inflammatory disease. In addition, liver metabolomic results revealed that metabolites related to steroid hormone biosynthesis, arachidonic acid metabolism, and linoleic acid metabolism were up-regulated and those related to valine, leucine, and isoleucine degradation pathways were down-regulated in NA treatment groups compared with the PD model. The integration of metabolomic and transcriptomic results showed NA significantly exerted its anti-inflammatory function by regulating the transcription and metabolic pathways of multiple genes. Particularly, linoleic acid metabolism, arachidonic acid metabolism, and steroid hormone biosynthesis were the crucial pathways of the anti-inflammatory action of NA. Key genes in these metabolic pathways and key molecules in inflammatory signaling pathways were also verified, which were consistent with transcriptomic results.

CONCLUSION

These findings provide novel insights into the liver protective effects of NA against PD mice. This study also showed that NA could be a useful dietary element for improving and treating PD-induced liver inflammation.

Changes in the content of pollen total lipid and TAG in Arabidopsis thaliana DGAT1 mutant as11

In mature pollen grains, lipids are primarily stored in the form of lipid droplets that provide energy and act as a carbon source for normal pollen development and germination. Triacylglycerol (TAG) is the major form of stored plant lipids. Diacylglycerol transferase, which is encoded by DGAT1 in Arabidopsis thaliana, is an important enzyme regulating triacylglycerol synthesis. Within the seeds of the DGAT1 mutant as11, the content of TAG is significantly decreased and the fatty acid composition also differs from the wild type. Transcriptome data of mature anthers showed that the genes involved in the TAG synthesis pathway were downregulated in as11. Analysis of gene expression patterns via transcriptome data also revealed that the expression of PDAT1, which functions in a manner complementary to the DGAT1 gene, was significantly decreased in as11, whereas the amylopectin synthase genes SS1 and SS2 were upregulated in mutant as11. We also detected lower total lipid, TAG and fatty acid contents in mature as11 pollen, with palmitic acid (C16:0) and linolenic acid (C18:3) being the major fatty acids in mature pollen. The cytological results showed that the lipid droplet content was reduced in mature as11 pollen. In the binuclear pollen grain II stage, WT pollen contained lipid droplets that were primarily accumulated around the generative nucleus, whereas the pollen in the mutant as11 was rich in starch grains that were primarily distributed around the vegetative nucleus. Ultrastructural analysis indicated that during pollen development in as11, the amount of endoplasmic reticulum in tapetal cells and pollen grains decreased, whereas the Golgi body content increased, which directly or indirectly led to a decrease in the levels of lipidosomes and an increase in the starch content in as11. Changes in the lipid content and fatty acid composition of the mutant as11 differ from those in the wild type during pollen development.

Integrated transcriptomic and metabolomic analyses revealed the molecular mechanism of terpenoid formation for salicylic acid resistance in Pulsatilla chinensis callus

As a kind of traditional Chinese medicine, Pulsatilla chinensis (Bunge) Regel is well known for its anti-inflammation and anti-cancer activities, which are attributed to its active components including total saponins and monomers. To clarify the synthesis and metabolism mechanisms of class components in callus terpenes of P. chinensis, a certain concentration of salicylic acid (SA) hormone elicitor was added to the callus before being analysed by transcriptomic and metabolomic techniques. Results showed that the content of Pulsatilla saponin B4 in the callus suspension culture was significantly increased up to 1.99% with the addition of SA. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the differentially expressed genes were mainly enriched in 122 metabolic pathways, such as terpenoid metabolism-related pathways: terpenoid skeleton synthesis pathway, monoterpenoid biosynthesis pathways, diterpenoid biosynthesis pathways, and ubiquinone and other terpenoid-quinone biosynthesis pathways. A total of 31 differentially accumulated metabolites were obtained from four differential groups. Amongst 21 kinds of known chemical components in P. chinensis, deoxyloganic acid was the only monoterpenoid; the others are triterpenoids. In summary, this study found that SA elicitors can affect the metabolic changes of terpenoids in P. chinensis callus, which provided a basis for analysing the genetic regulation of terpenoid components of leucons.

Transcriptome reveals insights into biosynthesis of ginseng polysaccharides

BACKGROUND

Ginseng polysaccharides, have been used to treat various diseases as an important active ingredient. Nevertheless, the biosynthesis of ginseng polysaccharides is poorly understood. To elucidate the biosynthesis mechanism of ginseng polysaccharides, combined the transcriptome analysis and polysaccharides content determination were performed on the roots, stems, and leaves collected from four cultivars of ginseng.

RESULTS

The results indicated that the total contents of nine monosaccharides were highest in the roots. Moreover, the total content of nine monosaccharides in the roots of the four cultivars were different but similar in stems and leaves. Glucose (Glc) was the most component of all monosaccharides. In total, 19 potential enzymes synthesizing of ginseng polysaccharides were identified, and 17 enzymes were significantly associated with polysaccharides content. Among these genes, the expression of phosphoglucomutase (PGM), glucose-6-phosphate isomerase (GPI), UTP-glucose-1-phosphate uridylyltransferase (UGP2), fructokinase (scrK), mannose-1-phosphate guanylyltransferase (GMPP), phosphomannomutase (PMM), UDP-glucose 4-epimerase (GALE), beta-fructofuranosidase (sacA), and sucrose synthase (SUS) were correlated with that of MYB, AP2/ERF, bZIP, and NAC transcription factors (TFs). These TFs may regulate the expression of genes involved in ginseng polysaccharides synthesis.

CONCLUSION

Our findings could provide insight into a better understanding of the regulatory mechanism of polysaccharides biosynthesis, and would drive progress in genetic improvement and plantation development of ginseng.

Comparative Analysis of Physiological, Hormonal and Transcriptomic Responses Reveal Mechanisms of Saline-Alkali Tolerance in Autotetraploid Rice (Oryza sativa L.)

Saline-alkali soil has posed challenges to the growth of agricultural crops, while polyploidy often show greater adaptability in diverse and extreme environments including saline-alkali stress, but its defense mechanisms in rice remain elusive. Herein, we explored the mechanisms of enhanced saline-alkali tolerance of autotetraploid rice 93-11T relative to diploid rice 93-11D, based on physiological, hormonal and transcriptomic profilings. Physiologically, the enhanced saline-alkali tolerance in 93-11T was manifested in higher soluble sugar accumulation and stronger superoxide dismutase (SOD) and peroxidase (POD) activities in leaves during 24 h after saline-alkali shock. Furthermore, various hormone levels in leaves of 93-11T altered greatly, such as the negative correlation between salicylic acid (SA) and the other four hormones changed to positive correlation due to polyploidy. Global transcriptome profiling revealed that the upregulated differentially expressed genes (DEGs) in leaves and roots of 93-11T were more abundant than that in 93-11D, and there were more DEGs in roots than in leaves under saline-alkali stress. Genes related to phytohormone signal transduction of auxin (AUX) and SA in roots, lignin biosynthesis in leaves or roots, and wax biosynthesis in leaves were obviously upregulated in 93-11T compared with 93-11D under saline-alkali condition. Collectively, 93-11T subjected to saline-alkali stress possibly possesses higher osmotic regulation ability due to cuticular wax synthesis, stronger negative regulation of reactive oxygen species (ROS) production by increasing the SA levels and maintaining relative lower levels of IAA, and higher antioxidant capacity by increasing activities of SOD and POD, as well as lignin biosynthesis. Our research provides new insights for exploring the mechanisms of saline-alkali tolerance in polyploid rice and discovering new gene targets for rice genetic improvement.

Integrated transcriptomic and metabolic analyses provide insights into the maintenance of embryogenic potential and the biosynthesis of phenolic acids and flavonoids involving transcription factors in Larix kaempferi (Lamb.) Carr

Somatic embryogenesis (SE) techniques have been established for micropropagation or basic research related to plant development in many conifer species. The frequent occurrence of non-embryogenic callus (NEC) during SE has impose constraints on the application of somatic embryogenesis SE in Larix kaempferi (Lamb.) Carr, but the potential regulatory mechanisms are poorly understood. In this study, integrated transcriptomic and metabolomic analyses were performed in embryogenic callus (EC) and NEC originating from a single immature zygotic embryo to better decipher the key molecular and metabolic mechanisms required for embryogenic potential maintenance. The results showed that a total of 13,842 differentially expressed genes (DEGs) were found in EC and NEC, among which many were enriched in plant hormone signal transduction, starch and sucrose metabolism, phenylpropanoid biosynthesis, flavonoid biosynthesis, and the biosynthesis of amino acids pathways. Metabolite profiling showed that 441 differentially accumulated metabolites (DAMs) were identified in EC and NEC. Both EC and NEC had vigorous primary metabolic activities, while most secondary metabolites were upregulated in NEC. Many totipotency-related transcription factor (TF) genes such as BBMs, WUSs, and LEC1 showed higher expression levels in EC compared with NEC, which may result in the higher accumulation of indole 3-acetic acid (IAA) in EC. NEC was characterized by upregulation of genes and metabolites associated with stress responses, such as DEGs involved in jasmonic acid (JA) and ethylene (ETH) biosynthesis and signal transduction pathways, and DEGs and DAMs related to phenylpropanoid and flavonoid biosynthesis. We predicted and analyzed TFs that could target several key co-expressed structural DEGs including two C4H genes, two CcoAOMT genes and three HCT genes involved in phenylpropanoid and flavonoid biosynthesis. Based on the targeted relationship and the co-expression network, two ERFs (Lk23436 and Lk458687), one MYB (Lk34626) and one C2C2-dof (Lk37167) may play an important role in regulating phenolic acid and flavonoid biosynthesis by transcriptionally regulating the expression of these structural genes. This study shows an approach involving integrated transcriptomic and metabolic analyses to obtain insights into molecular events underlying embryogenic potential maintenance and the biosynthesis mechanisms of key metabolites involving TF regulation, which provides valuable information for the improvement of SE efficiency in L. kaempferi.

Integrated analyses reveal the response of peanut to phosphorus deficiency on phenotype, transcriptome and metabolome

BACKGROUND

Phosphorus (P) is one of the most essential macronutrients for crops. The growth and yield of peanut (Arachis hypogaea L.) are always limited by P deficiency. However, the transcriptional and metabolic regulatory mechanisms were less studied. In this study, valuable phenotype, transcriptome and metabolome data were analyzed to illustrate the regulatory mechanisms of peanut under P deficiency stress.

RESULT

In present study, two treatments of P level in deficiency with no P application (-P) and in sufficiency with 0.6 mM P application (+ P) were used to investigate the response of peanut on morphology, physiology, transcriptome, microRNAs (miRNAs), and metabolome characterizations. The growth and development of plants were significantly inhibited under -P treatment. A total of 6088 differentially expressed genes (DEGs) were identified including several transcription factor family genes, phosphate transporter genes, hormone metabolism related genes and antioxidant enzyme related genes that highly related to P deficiency stress. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that 117 genes were annotated in the phenylpropanoid biosynthesis pathway under P deficiency stress. A total of 6 miRNAs have been identified significantly differential expression between + P and -P group by high-throughput sequencing of miRNAs, including two up-regulated miRNAs (ahy-miR160-5p and ahy-miR3518) and four down-regulated miRNAs (ahy-miR408-5p, ahy-miR408-3p, ahy-miR398, and ahy-miR3515). Further, the predicted 22 target genes for 6 miRNAs and cis-elements in 2000 bp promoter region of miRNA genes were analyzed. A total of 439 differentially accumulated metabolites (DAMs) showed obviously differences in two experimental conditions.

CONCLUSIONS

According to the result of transcripome and metabolome analyses, we can draw a conclusion that by increasing the content of lignin, amino acids, and levan combining with decreasing the content of LPC, cell reduced permeability, maintained stability, raised the antioxidant capacity, and increased the P uptake in struggling for survival under P deficiency stress.

Comparative transcriptome analysis of Lupinus polyphyllus Lindl. provides a rich molecular resource for research on coloration mechanism

Lupinus polyphyllus is rich in color, making it a well-known horticultural ornamental plant. However, little is known about the genes related to anthocyanin and carotenoid biosynthesis in L. polyphyllus. In this study, transcriptome sequencing was performed on eight different colors of L. polyphyllus. A total of 1.13 billion clean reads were obtained and assembled into 89,124 unigenes, which were then aligned with six databases, resulting in the identification of 54,823 annotated unigenes. Among these unigenes, 76 and 101 were involved in the biosynthetic pathway of carotenoids and anthocyanins, respectively. In addition, 505 transcription factors were revealed, which belonged to the MYB, R2R3-MYB, NAC, bHLH, and WD40 families. A total of 6,700 differentially expressed genes (DEGs) were obtained by comparative transcriptome analysis. Among them, 17 candidate unigenes (four carotenoid genes, seven anthocyanin genes, and six TFs) were specifically up-regulated for one or more colors of L. polyphyllus. Eight representative candidate unigenes were analyzed by qRT-PCR. The findings enrich the transcriptome database of lupine, and provide a rich molecular resource for research on the coloration mechanism of L. polyphyllus.

Transcriptomic and Metabolomic Analyses Provide Insights into the Formation of the Peach-like Aroma of Fragaria nilgerrensis Schlecht. Fruits

Fragaria nilgerrensis Schlecht. is a wild diploid strawberry species. The intense peach-like aroma of its fruits makes F. nilgerrensis an excellent resource for strawberry breeding programs aimed at enhancing flavors. However, the formation of the peach-like aroma of strawberry fruits has not been comprehensively characterized. In this study, fruit metabolome and transcriptome datasets for F. nilgerrensis (HA; peach-like aroma) and its interspecific hybrids PA (peach-like aroma) and NA (no peach-like aroma; control) were compared. In total, 150 differentially accumulated metabolites were detected. The K-means analysis revealed that esters/lactones, including acetic acid, octyl ester, δ-octalactone, and δ-decalactone, were more abundant in HA and PA than in NA. These metabolites may be important for the formation of the peach-like aroma of F. nilgerrensis fruits. The significantly enriched gene ontology terms assigned to the differentially expressed genes (DEGs) were fatty acid metabolic process and fatty acid biosynthetic process. Twenty-seven DEGs were predicted to be associated with ester and lactone biosynthesis, including AAT, LOX, AOS, FAD, AIM1, EH, FAH, ADH, and cytochrome P450 subfamily genes. Thirty-five transcription factor genes were predicted to be associated with aroma formation, including bHLH, MYB, bZIP, NAC, AP2, GATA, and TCPfamily members. Moreover, we identified differentially expressed FAD, AOS, and cytochrome P450 family genes and NAC, MYB, and AP2 transcription factor genes that were correlated with δ-octalactone and δ-decalactone. These findings provide key insights into the formation of the peach-like aroma of F. nilgerrensis fruits, with implications for the increased use of wild strawberry resources.