Comparative transcriptomes and WGCNA reveal hub genes for spike germination in different quinoa lines

Abstract

Background Quinoa, as a new food crop, has attracted extensive attention at home and abroad. However, the natural disaster of spike germination seriously threatens the quality and yield of quinoa. Currently, there are limited reports on the molecular mechanisms associated with spike germination in quinoa. Results In this study, we utilized transcriptome sequencing technology and successfully obtained 154.51 Gb of high-quality data with a comparison efficiency of more than 88%, which fully demonstrates the extremely high reliability of the sequencing results and lays a solid foundation for subsequent analysis. Using these data, we constructed a weighted gene co-expression network (WGCNA) related to starch, sucrose, α-amylase, and phenolic acid metabolites, and screened six co-expression modules closely related to spike germination traits. Two of the modules associated with physiological indicators were analyzed in depth, and nine core genes were finally predicted. Further functional annotation revealed four key transcription factors involved in the regulation of dormancy and germination processes: gene LOC110698065, gene LOC110696037, gene LOC110736224, and gene LOC110705759, belonging to the bHLH, NF-YA, MYB, and FAR1 gene families, respectively. Conclusions These results provide clues to identify the core genes involved in quinoa spike germination. This will ultimately provide a theoretical basis for breeding new quinoa varieties with resistance.