Integrated transcriptomic and metabolomic analyses reveal key genes controlling flavonoid biosynthesis in Citrus grandis ‘Tomentosa’ fruits

As a well-recognized traditional Chinese medicine (TCM), immature fruits of Citrus grandis 'Tomentosa' (CGT) serve to cure chronic cough in humans. Specialized metabolites including flavonoids may have contribute to this curing effect. Knowledge about the molecular mechanisms underlying flavonoid biosynthesis in 'Tomentosa' fruits will, therefore, support the breeding of varieties with improved medicinal properties. Hence, we profiled the transcriptomes and metabolites of the fruits of two contrasting C. grandis varieties, namely 'Zheng-Mao' ('ZM') used in TCM production, and a locally cultivated pomelo, namely 'Guang-Qing' ('GQ'), at four developmental stages. A total of 39 flavonoids, including 14 flavanone/flavone, 5 isoflavonoids, 12 flavonols, and 6 anthocyanins, were identified, and 16 of which were quantitatively determined in the fruits of the two varieties. We found that 'ZM' fruits contain more flavonoids than 'GQ'. Specifically, rhoifolin levels were significantly higher in 'ZM' than in 'GQ'. We annotated 31,510 genes, including 1,387 previously unknown ones, via transcriptome sequencing of 'ZM' and 'GQ.' A total of 646 genes were found to be differentially expressed between 'ZM' and 'GQ' throughout at all four fruit developmental stages, indicating that they are robust expression markers for future breeding programs. Weighted gene co-expression network analysis identified 18 modules. Combined transcriptional and metabolic analysis revealed 25 genes related to flavonoid biosynthesis and 16 transcriptional regulators (MYBs, bHLHs, WD40) that may be involved in the flavonoids biosynthesis in C. grandis 'Tomentosa' fruits.