Integrated Physiochemical, Hormonal, and Transcriptomic Analysis Revealed the Underlying Mechanisms for Granulation in Huyou (Citrus changshanensis) Fruit

Juice sac granulation is a common internal physiological disorder of citrus fruit. In the present study, we compared the physiochemical characteristics and transcriptome profiles of juice sacs in different granulation levels from Huyou fruit (Citrus changshanensis). The accumulation of cell wall components, including the water-soluble pectin, protopectin, cellulose, and lignin, were significantly correlated with the granulation process, resulting in the firmness increase of the juice sac. The in situ labeling of the cell wall components indicated the early accumulation of cellulose and high-methylesterified pectin in the outer layer cells, as well as the late accumulation of lignin in the inner layer cells of the juice sac. Several phytohormones, including auxins, abscisic acids, cytokinins, jasmonic acid, salicylic acid, and/or their metabolites, were positively correlated to the granulation level, indicating an active and complex phytohormones metabolism in the granulation process. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 2940 differentially expressed genes (DEGs) were found to be positively correlated with the granulation level. Gene Ontology (GO) enrichment indicated that the selected DEGs were mainly involved in the cell wall organization and biogenesis, cell wall macromolecule metabolic process, carbohydrate metabolic process, and polysaccharide metabolic process. Among these selected genes, those encoding β-1,4-xylosyltransferase IRX9, cellulose synthase, xyloglucan: xyloglucosyl transferase, xyloglucan galactosyltransferase MUR3, α-1,4-galacturonosyltransferase, expansin, polygalacturonase, pectinesterase, β-glucosidase, β-galactosidase, endo-1,3(4)-β-glucanase, endoglucanase and pectate lyase that required for the biosynthesis or structural modification of cell wall were identified. In addition, NAC, MYB, bHLH, and MADS were the top abundant transcription factors (TFs) families positively correlated with the granulation level, while the LOB was the top abundant TFs family negatively correlated with the granulation level.