Weighted gene co-expression network analysis unveils gene networks regulating folate biosynthesis in maize endosperm

Folates are essential elements for human growth and development, and their deficiency can lead to serious disorders. Waxy maize is a rich source of folates; however, the regulatory mechanism underlying folate biosynthesis in the endosperm remains unclear. Here, we examined changes in the folate content of maize endosperm collected at 15, 18, 21, 24, and 27 days after pollination (DAP) using liquid chromatograph-mass spectrometry and identified genes related to folate biosynthesis using transcriptome sequencing data. The results showed that 5-methyl-tetrahydrofolate and 5,10-methylene tetrahydrofolate were the main storage forms of folates in the endosperm, and their contents were relatively high at 21–24 days. We also identified 569, 3183, 4365, and 5513 differentially expressed genes (DEGs) in different days around milk stage. Functional annotation revealed 518 transcription factors (TFs) belonging to 33 families exhibiting specific expression in at least one sampling time. The key hub genes involved in folate biosynthesis were identified by weighted gene co-expression network analysis. In total, 24,976 genes were used to construct a co-expression network with 29 co-expression modules, among which the brown and purple modules were highly related to folate biosynthesis. Further, 187 transcription factors in the brown and purple modules were considered potential transcription factors related to endosperm folate biosynthesis. These results may improve the understanding of the molecular mechanism underlying folate biosynthesis in waxy maize and lead to the development of nutritionally fortified varieties.