Starch is a major component of crop grains, and its content affects food quality and taste. Tartary buckwheat is a traditional pseudo-cereal used in food as well as medicine. Starch content, granule morphology, and physicochemical properties have been extensively studied in Tartary buckwheat. However, the complex regulatory network related to its starch biosynthesis needs to be elucidated. Here, we performed RNA-seq analyses using seven Tartary buckwheat varieties differing in starch content and combined the RNA-seq data with starch content by weighted correlation network analysis (WGCNA). As a result, 10,873 differentially expressed genes (DEGs) were identified and were functionally clustered to six hierarchical clusters. Fifteen starch biosynthesis genes had higher expression level in seeds. Four trait-specific modules and 3131 hub genes were identified by WGCNA, with the lightcyan and brown modules positively correlated with starch-related traits. Furthermore, two potential gene regulatory networks were proposed, including the co-expression of FtNAC70, FtPUL, and FtGBSS1-3 in the lightcyan module and FtbHLH5, C3H, FtBE2, FtISA3, FtSS3-5, and FtSS1 in the brown. All the above genes were preferentially expressed in seeds, further suggesting their role in seed starch biosynthesis. These results provide crucial guidance for further research on starch biosynthesis and its regulatory network in Tartary buckwheat.
Keywords: DEGs; RNA-seq; Tartary buckwheat; amylopectin; amylose; gene regulatory network; starch biosynthesis; transcription factor.