Transcriptome Analysis Reveals Candidate Lignin-Related Genes and Transcription Factors during Fruit Development in Pomelo (Citrus maxima)

Juice sac granulation (a physiological disorder) leads to large postharvest losses of pomelo (Citrus maxima). Previous studies have shown that juice sac granulation is closely related to lignin accumulation, while the molecular mechanisms underlying this disorder remain elusive in pomelo. Our results showed that the lignin content in NC (near the core) and FC (far away from the core) juice sacs overall increased from 157 DPA (days post anthesis) to 212 DPA and reached a maximum at 212 DPA. Additionally, the lignin content of NC juice sacs was higher than that of FC juice sacs. In this study, we used transcriptome-based weighted gene co-expression network analysis (WGCNA) to address how lignin formation in NC and FC juice sacs is generated during the development of pomelo. After data assembly and bioinformatic analysis, we found a most correlated module (black module) to the lignin content, then we used the 11 DEGs in this module as hub genes for lignin biosynthesis. Among these DEGs, PAL (phenylalanine ammonia lyase), HCT (hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase), 4CL2 (4-coumarate: CoA ligase), C4H (cinnamate 4-hydroxylase), C3′H (p-coumarate 3-hydroxylase), and CCoAOMT1 (caffeoyl CoA 3-Omethyltransferase) were the most distinct DEGs in granulated juice sacs. Co-expression analysis revealed that the expression patterns of several transcription factors such as MYB, NAC, OFP6, and bHLH130 are highly correlated with lignin formation. In addition, the expression patterns of the DEGs related to lignin biosynthesis and transcription factors were validated by qRT-PCR, and the results were highly concordant with the RNA-seq results. These results would be beneficial for further studies on the molecular mechanism of lignin accumulation in pomelo juice sacs and would help with citrus breeding.