Transcriptome reveals the gene expression patterns of sulforaphane metabolism in broccoli florets

Sulforaphane is a new and effective anti-cancer component that is abundant in broccoli. In the past few years, the patterns of variability in glucosinolate content and its regulation in A. thaliana have been described in detail. However, the diversity of glucosinolate and sulforaphane contents in different organs during vegetative and reproductive stages has not been clearly explained. In this paper, we firstly investigated the transcriptome profiles of the developing buds and leaves at bolting stage of broccoli (B52) to further assess the gene expression patterns involved in sulforaphane synthesis. The CYP79F1 gene, as well as nine other genes related to glucorahpanin biosynthesis, MAM1, MAM3, St5b-2, FMO GSOX1, MY, AOP2, AOP3, ESP and ESM1 were selected by digital gene expression analysis and were validated by quantitative real-time PCR (qRT-PCR). Meanwhile, the compositions of glucosinolates and sulforaphane were detected for correlation analysis with related genes. Finally the RNA sequencing libraries generated 147 957 344 clean reads, and 8 539 unigene assemblies were produced. In digital result, only CYP79F1, in the glucoraphanin pathway, was up-regulated in young buds but absent from the other organs, which was consistent with the highest level of sulforaphane content being in this organ compared to mature buds, buds one day before flowering, flowers and leaves. The sequencing results also presented that auxin and cytokinin might affect glucoraphanin accumulation. The study revealed that up-regulated expression of CYP79F1 plays a fundamental and direct role in sulforaphane production in inflorescences. Two genes of MAM1 and St5b-2 could up-regulated glucoraphanin generation. Synergistic expression of MAM1, MAM3, St5b-2, FMO GSOX1, MY, ESP and ESM1 was found in sulforaphane metabolism. This study will be beneficial for understanding the diversity of sulforaphane in broccoli organs.