Comparative transcriptome and weighted correlation network analyses reveal candidate genes involved in chlorogenic acid biosynthesis in sweet potato.

Jing Xu, Jiahong Zhu, Yanhui Lin, Honglin Zhu, Liqiong Tang, Xinhua Wang, Xiaoning Wang
Author Information
  1. Jing Xu: Institute of Cereal Crops, Hainan Academy of Agricultural Sciences, Key Laboratory of Crop Genetics and Breeding of Hainan Province, Haikou, 571100, China.
  2. Jiahong Zhu: Institute of Tropical Biosciences and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou, 571101, China. zhujiahong@itbb.org.cn.
  3. Yanhui Lin: Institute of Cereal Crops, Hainan Academy of Agricultural Sciences, Key Laboratory of Crop Genetics and Breeding of Hainan Province, Haikou, 571100, China.
  4. Honglin Zhu: Institute of Cereal Crops, Hainan Academy of Agricultural Sciences, Key Laboratory of Crop Genetics and Breeding of Hainan Province, Haikou, 571100, China.
  5. Liqiong Tang: Institute of Cereal Crops, Hainan Academy of Agricultural Sciences, Key Laboratory of Crop Genetics and Breeding of Hainan Province, Haikou, 571100, China.
  6. Xinhua Wang: Institute of Cereal Crops, Hainan Academy of Agricultural Sciences, Key Laboratory of Crop Genetics and Breeding of Hainan Province, Haikou, 571100, China.
  7. Xiaoning Wang: Institute of Cereal Crops, Hainan Academy of Agricultural Sciences, Key Laboratory of Crop Genetics and Breeding of Hainan Province, Haikou, 571100, China. wxning2599@163.com.

Abstract

Chlorogenic acids (CGAs) are important secondary metabolites produced in sweet potato. However, the mechanisms of their biosynthesis and regulation remain unclear. To identify potential genes involved in CGA biosynthesis, analysis of the dynamic changes in CGA components and RNA sequencing were performed on young leaves (YL), mature leaves (ML), young stems (YS), mature stems (MS) and storage roots (SR). Accordingly, we found that the accumulation of six CGA components varied among the different tissues and developmental stages, with YS and YL recording the highest levels, while SR exhibited low levels. Moreover, the transcriptome analysis yielded 59,287 unigenes, 3,767 of which were related to secondary-metabolite pathways. The differentially expressed genes (DEGs) were identified based on CGA content levels by comparing the different samples, including ML vs. YL, MS vs. YS, SR vs. YL and SR vs. YS. A total of 501 common DEGs were identified, and these were mainly implicated in the secondary metabolites biosynthesis. Additionally, eight co-expressed gene modules were identified following weighted gene co-expression network analysis, while genes in darkgrey module were highly associated with CGA accumulation. Darkgrey module analysis revealed that 12 unigenes encoding crucial enzymes (PAL, 4CL, C4H, C3H and HCT/HQT) and 42 unigenes encoding transcription factors (MYB, bHLH, WD40, WRKY, ERF, MADS, GARS, bZIP and zinc finger protein) had similar expression patterns with change trends of CGAs, suggesting their potential roles in CGA metabolism. Our findings provide new insights into the biosynthesis and regulatory mechanisms of CGA pathway, and will inform future efforts to build a genetically improve sweet potato through the breeding of high CGA content varieties.

References

  1. Payyavula, R. S. et al. Synthesis and regulation of chlorogenic acid in potato: Rerouting phenylpropanoid flux in HQT-silenced lines. Plant Biotechnol. J. 13, 551–564 (2015). [PMID: 25421386]
  2. Clifford, M. N., Jaganath, B., Ludwig, I. A. & Crozier, A. Chlorogenic acids and the acyl-quinic acids: Discovery, biosynthesis, bioavailability and bioactivity. Nat. Prod. Rep. 34, 1391–1421 (2017). [PMID: 29160894]
  3. Clifford, M. N., Kerimi, A. & Williamson, G. Bioavailability and metabolism of chlorogenic acids (acyl-quinic acids) in humans. Compr. Rev. Food. Sci. Food Saf. 19, 1299–1352 (2020). [PMID: 33337099]
  4. He, L. et al. Transcriptome analysis of buds and leaves using 454 pyrosequencing to discover genes associated with the biosynthesis of active ingredients in Lonicera japonica Thunb. PLoS ONE 8, e62922 (2013). [PMID: 23638167]
  5. Niggeweg, R., Michael, A. J. & Martin, C. Engineering plants with increased levels of the antioxidant chlorogenic acid. Nat. Biotechnol. 22, 746–754 (2004). [PMID: 15107863]
  6. Villegas, R. J. A. & Kojima, M. Purification and characterization of hydroxycinnamoyl D-glucose quinate hydroxycinnamoyl transferase in the root of sweet potato, Ipomoea batatas Lam. J. Biol. Chem. 261, 8729–8733 (1986). [PMID: 3722170]
  7. Hoffmann, L., Maury, S., Martz, F., Geoffroy, P. & Legrand, M. Purification, cloning, and properties of an acyltransferase controlling shikimate and quinate ester intermediates in phenylpropanoid metabolism. J. Biol. Chem. 278, 95–103 (2003). [PMID: 12381722]
  8. Tuominen, L. K., Johnson, V. E. & Tsai, C. J. Differential phylogenetic expansions in BAHD acyltransferases across five angiosperm taxa and evidence of divergent expression among Populus paralogues. BMC Genom. 12, 236 (2011). [DOI: 10.1186/1471-2164-12-236]
  9. Cheevarungnapakul, K., Khaksar, G., Panpetch, P., Boonjing, P. & Sirikantaramas, S. Identification and functional characterization of genes involved in the biosynthesis of caffeoylquinic acids in Sunflower (Helianthus annuus L.). Front. Plant Sci. 10, 968 (2019). [PMID: 31417585]
  10. Zhang, J. et al. Genome-wide association studies and expression-based quantitative trait loci analyses reveal roles of HCT2 in caffeoylquinic acid biosynthesis and its regulation by defense-responsive transcription factors in Populus. New Phytol. 220, 502–516 (2018). [PMID: 29992670]
  11. Moglia, A. et al. Dual catalytic activity of hydroxycinnamoyl-coenzyme A quinate transferase from tomato allows it to moonlight in the synthesis of both mono- and dicaffeoylquinic acids. Plant Physiol. 166, 1777–1787 (2014). [PMID: 25301886]
  12. Moglia, A. et al. Genome-wide identification of BAHD acyltransferases and in vivo characterization of HQT-like enzymes involved in caffeoylquinic acid synthesis in Globe Artichoke. Front. Plant Sci. 7, 1424 (2016). [PMID: 27721818]
  13. Chen, Z. X., Liu, G. H., Liu, Y. Q., Xian, Z. Q. & Tang, N. Overexpression of the LmHQT1 gene increases chlorogenic acid production in Lonicera macranthoides Hand-Mazz. Acta. Physiol. Plant 39, 27 (2016). [DOI: 10.1007/s11738-016-2310-8]
  14. Mudau, S. P. et al. Metabolomics-guided investigations of unintended effects of the expression of the hydroxycinnamoyl quinate hydroxycinnamoyltransferase (hqt1) gene from Cynara cardunculus var scolymus in Nicotiana tabacum cell cultures. Plant Physiol. Biochem. 127, 287–298 (2018). [PMID: 29649745]
  15. Li, Y. et al. Correlation of the temporal and spatial expression patterns of HQT with the biosynthesis and accumulation of chlorogenic acid in Lonicera japonica flowers. Hortic. Res. 6, 73 (2019). [PMID: 31231531]
  16. Sheikha, A. F. E. I. & Ray, R. C. Potential impacts of bioprocessing of sweet potato: review. Crit. Rev. Food Sci. Nutr. 57, 455–471 (2017). [PMID: 25975980]
  17. Mohanraj, R. & Sivasankar, S. Sweet potato (Ipomoea batatas [L.] Lam)–a valuable medicinal food: A review. J. Med. Food. 17, 733–741 (2014). [PMID: 24921903]
  18. Zheng, W. & Clifford, M. N. Profiling the chlorogenic acids of sweet potato (Ipomoea batatas) from China. Food Chem. 106, 147–152 (2008). [DOI: 10.1016/j.foodchem.2007.05.053]
  19. Tanaka, M., Ishiguro, K., Oki, T. & Okuno, S. Functional components in sweet potato and their genetic improvement. Breed. Sci. 67, 52–61 (2017). [PMID: 28465668]
  20. Liao, Y., Zeng, L., Rao, S., Gu, D. & Yang, Z. Induced biosynthesis of chlorogenic acid in sweetpotato leaves confers the resistance against sweet potato weevil attack. J. Adv. Res. 24, 513–522 (2020). [PMID: 32612857]
  21. Yu, Y. et al. Overexpression of IbPAL1 promotes chlorogenic acid biosynthesis in sweet potato. Crop J. 9, 204–215 (2021). [DOI: 10.1016/j.cj.2020.06.003]
  22. Grienenberger, E. et al. A BAHD acyltransferase is expressed in the tapetum of Arabidopsis anthers and is involved in the synthesis of hydroxycinnamoyl spermidines. Plant J. 58, 246–259 (2009). [PMID: 19077165]
  23. Santana-Gálvez, J., Cisneros-Zevallos, L. & Jacobo-Velázquez, D. A. Chlorogenic acid: Recent advances on its dual role as a food additive and a nutraceutical against metabolic syndrome. Molecules 22, 358 (2017). [>PMCID: ]
  24. Bate, N. J. et al. Quantitative relationship between phenylalanine ammonia-lyase levels and phenylpropanoid accumulation in transgenic tobacco identifies a rate-determining step in natural product synthesis. Proc. Natl. Acad. Sci. USA 91, 7608–7612 (1994). [PMID: 8052628]
  25. Shadle, G. L. et al. Phenylpropanoid compounds and disease resistance in transgenic tobacco with altered expression of L-phenylalanine ammonia-lyase. Phytochemistry 64, 153–161 (2003). [PMID: 12946414]
  26. Zhang, X. & Liu, C. J. Multifaceted regulations of gateway enzyme phenylalanine ammonia-lyase in the biosynthesis of phenylpropanoids. Mol. Plant 8, 17–27 (2015). [PMID: 25578269]
  27. Masuo, S., Kobayashi, Y., Oinuma, K. & Takaya, N. Alternative fermentation pathway of cinnamic acid production via phenyllactic acid. Appl. Microbiol. Biotechnol. 100, 8701–8709 (2016). [PMID: 27225472]
  28. Chen, Z. et al. Transcriptome analysis reveals the mechanism underlying the production of a high quantity of chlorogenic acid in young leaves of Lonicera macranthoides Hand.-Mazz. PLoS ONE 10, 7212 (2015).
  29. Kim, Y. B. et al. Characterization of genes for a putative hydroxycinnamoyl-coenzyme A quinate transferase and p-coumarate 3’-hydroxylase and chlorogenic acid accumulation in tartary buckwheat. J. Agric. Food Chem. 61, 4120–4126 (2013). [PMID: 23550515]
  30. Knollenberg, B. J., Liu, J., Yu, S., Lin, H. & Tian, L. Cloning and functional characterization of a p-coumaroyl quinate/shikimate 3’-hydroxylase from potato (Solanum tuberosum). Biochem. Biophys. Res. Commun. 496, 462–467 (2018). [PMID: 29337064]
  31. Legrand, G. et al. Identification and characterization of five BAHD acyltransferases involved in hydroxycinnamoyl ester metabolism in Chicory. Front. Plant Sci. 7, 741 (2016). [PMID: 27375627]
  32. Strack, D. & Gross, W. Properties and activity changes of chlorogenic acid: Glucaric acid caffeoyl-transferase from tomato (Lycopersicon esculentum). Plant Physiol. 92, 41–47 (1990). [PMID: 16667263]
  33. Dong, N. Q. & Lin, H. X. Contribution of phenylpropanoid metabolism to plant development and plant-environment interactions. J. Integr. Plant Biol. 63, 180–209 (2021). [PMID: 33325112]
  34. Tuan, P. A. et al. Enhancement of chlorogenic acid production in hairy roots of Platycodon grandiflorum by over-expression of an Arabidopsis thaliana transcription factor AtPAP1. Int. J. Mol. Sci. 15, 14743–14752 (2014). [PMID: 25153629]
  35. Sitarek, P. et al. Over-expression of AtPAP1 transcriptional factor enhances phenolic acid production in transgenic roots of Leonurus sibiricus L. and their biological activities. Mol. Biotechnol. 60, 74–82 (2018). [PMID: 29196986]
  36. Mehrtens, F., Kranz, H., Bednarek, P. & Weisshaar, B. The Arabidopsis transcription factor MYB12 is a flavonol-specific regulator of phenylpropanoid biosynthesis. Plant Physiol. 138, 1083–1096 (2005). [PMID: 15923334]
  37. Luo, J. et al. AtMYB12 regulates caffeoyl quinic acid and flavonol synthesis in tomato: Expression in fruit results in very high levels of both types of polyphenol. Plant J. 56, 316–326 (2008). [PMID: 18643978]
  38. Pandey, A. et al. AtMYB12 expression in tomato leads to large scale differential modulation in transcriptome and flavonoid content in leaf and fruit tissues. Sci. Rep. 5, 12412 (2015). [PMID: 26206248]
  39. Feyissa, D. N., Løvdal, T., Olsen, K. M., Slimestad, R. & Lillo, C. The endogenous GL3, but not EGL3, gene is necessary for anthocyanin accumulation as induced by nitrogen depletion in Arabidopsis rosette stage leaves. Planta 230, 747–754 (2009). [PMID: 19621239]
  40. Liu, C. C. et al. The bZip transcription factor HY5 mediates CRY1a-induced anthocyanin biosynthesis in tomato. Plant Cell Environ. 41, 1762–1775 (2018). [PMID: 29566255]
  41. Zha, L. et al. DNA methylation influences chlorogenic acid biosynthesis in Lonicera japonica by mediating LjbZIP8 to regulate phenylalanine ammonia-lyase 2 expression. Front. Plant Sci. 8, 1178 (2017). [PMID: 28740500]
  42. Habib, S., Lwin, Y. Y. & Li, N. Down-regulation of SlGRAS10 in tomato confers abiotic stress tolerance. Genes (Basel) 12, 623 (2021). [DOI: 10.3390/genes12050623]
  43. Zhao, C. et al. Three AP2/ERF family members modulate flavonoid synthesis by regulating type IV chalcone isomerase in citrus. Plant Biotechnol. J. 19, 671–688 (2021). [PMID: 33089636]
  44. Ma, R. et al. AP2/ERF Transcription factor, Ii049, positively regulates lignan biosynthesis in Isatis indigotica through activating salicylic acid signaling and lignan/lignin pathway genes. Front. Plant Sci. 8, 1361 (2017). [PMID: 28824690]
  45. Garceau, D. C., Batson, M. K. & Pan, I. L. Variations on a theme in fruit development: The PLE lineage of MADS-box genes in tomato (TAGL1) and other species. Planta 246, 313–321 (2017). [PMID: 28660293]
  46. Tao, Z. et al. BnLATE, a Cys2/His2-type zinc-finger protein, enhances silique shattering resistance by negatively regulating lignin accumulation in the silique walls of Brassica napus. PLoS ONE 12, e0168046 (2017). [PMID: 28081140]
  47. Lian, J. et al. Silencing of BnTT1 family genes affects seed flavonoid biosynthesis and alters seed fatty acid composition in Brassica napus. Plant Sci. 254, 32–47 (2017). [PMID: 27964783]
  48. Tian, H. & Wang, S. TRANSPARENT TESTA GLABRA1, a Key Regulator in plants with multiple roles and multiple function mechanisms. Int. J. Mol. Sci. 21, 4881 (2020). [>PMCID: ]
  49. Grabherr, M. G. et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 29, 644–652 (2011). [PMID: 21572440]
  50. Pertea, G. et al. TIGR Gene Indices clustering tools (TGICL): A software system for fast clustering of large EST datasets. Bioinformatics 19, 651–652 (2003). [PMID: 12651724]
  51. Deng, Y. et al. Integrated nr database in protein annotation system and its localization. Comput. Eng. 32, 71–74 (2006).
  52. Ashburner, M. et al. Gene ontology: Tool for the unification of biology. The gene ontology consortium. Nat. Genet. 25, 25–29 (2000). [PMID: 10802651]
  53. Koonin, E. V. et al. A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes. Genome Biol. 5, R7 (2004). [PMID: 14759257]
  54. Huerta-Cepas, J. et al. eggNOG 45: A hierarchical orthology framework with improved functional annotations for eukaryotic, prokaryotic and viral sequences. Nucl. Acids Res. 44, D286–D293 (2016). [PMID: 26582926]
  55. Finn, R. D. et al. Pfam: The protein families database. Nucleic Acids Res. 42, D222–D230 (2014). [PMID: 24288371]
  56. Apweiler, R. et al. UniProt: the Universal Protein knowledgebase. Nucl. Acids Res. 32, D115–D119 (2004). [PMID: 14681372]
  57. Kanehisa, M., Goto, S., Kawashima, S., Okuno, Y. & Hattori, M. The KEGG resource for deciphering the genome. Nucl. Acids Res. 32, D277–D280 (2004). [PMID: 14681412]
  58. Tatusov, R. L. et al. The COG database: an updated version includes eukaryotes. BMC Bioinf. 4, 41 (2003). [DOI: 10.1186/1471-2105-4-41]
  59. Langfelder, P. & Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinf. 9, 559 (2008). [DOI: 10.1186/1471-2105-9-559]
  60. Schmittgen, T. D. & Livak, K. J. Analyzing real-time PCR data by the comparative CT method. Nat. Methods 3, 1101–1108 (2008).

MeSH Term

Chlorogenic Acid
Gene Expression Regulation, Plant
Ipomoea batatas
Plant Leaves
Plant Proteins
Plant Roots
RNA-Seq
Transcriptome

Chemicals

Plant Proteins
Chlorogenic Acid

Word Cloud

Created with Highcharts 10.0.0CGAbiosynthesisgenesanalysisYLYSSRvssweetpotatolevelsunigenesidentifiedCGAssecondarymetabolitesmechanismspotentialinvolvedcomponentsyoungleavesmatureMLstemsMSaccumulationdifferenttranscriptomegeneweightednetworkmoduleencodingChlorogenicacidsimportantproducedHoweverregulationremainunclearidentifydynamicchangesRNAsequencingperformedstoragerootsAccordinglyfoundsixvariedamongtissuesdevelopmentalstagesrecordinghighestexhibitedlowMoreoveryielded592873767relatedsecondary-metabolitepathwaysdifferentiallyexpressedDEGsbasedcontentcomparingsamplesincludingtotalof 501common DEGsthese were mainlyimplicatedAdditionallyeightco-expressedmodulesfollowingco-expressiondarkgreyhighlyassociatedwith CGADarkgreyrevealed12crucialenzymesPAL4CLC4HC3HHCT/HQT42transcriptionfactorsMYBbHLHWD40WRKYERFMADSGARSbZIPzincfingerproteinsimilarexpressionpatternschangetrendssuggestingrolesmetabolismfindingsprovidenewinsightsregulatorypathwaywillinformfutureeffortsbuildgeneticallyimprovebreeding of high CGA contentvarietiesComparativecorrelationanalysesrevealcandidatechlorogenicacid

Similar Articles

Cited By