Investigation into the underlying regulatory mechanisms shaping inflorescence architecture in Chenopodium quinoa

Abstract

Background: Inflorescence architecture is denoted by the spatial arrangement of various lateral branches and florets formed on them, which is shaped by a complex of regulators. Unveiling of the regulatory mechanisms underlying inflorescence architecture is pivotal for improving crop yield potential. Quinoa (Chenopodium quinoa Willd), a pseudo cereal originated from Andean region of South America, has been widely recognized as a functional super food due to its excellent nutritional elements. Increasing worldwide consumption of this crop urgently calls for its yield improvement. However, dissection of the regulatory networks underlying quinoa inflorescence patterning is lacking. Results: In this study, we performed RNA-seq analysis on quinoa inflorescence samples collected from six developmental stages, yielding a total of 138.8 GB data. We screened 21,610 differentially expressed genes (DEGs) among all the stages through comparative analysis. Weighted Gene Co-Expression Network Analysis (WGCNA) was performed to categorize the DEGs into ten different modules. Subsequently, we placed emphasis on investigating the modules associated with none branched and branched inflorescence samples. We manually refined the coexpression networks with stringent edge weight cutoffs, and generated core networks using transcription factors and key inflorescence architecture related genes as seed nodes. The core networks were visualized and analyzed by Cytoscape to obtain hub genes in each network. Our finding indicates that the specific occurrence of B3, TALE, WOX, LSH, LFY, GRAS, bHLH, EIL, DOF, G2-like and YABBY family members in early reproductive stage modules, and of TFL, ERF, bZIP, HD-ZIP, C2H2, LBD, NAC, C3H, Nin-like and FAR1 family members in late reproductive stage modules, as well as the several different MADS subfamily members identified in both stages may account for shaping quinoa inflorescence architecture. Conclusion: In this study we carried out comparative transcriptome analysis of six different stages quinoa inflorescences, and using WGCNA we obtained the most highly potential central hubs for shaping inflorescence. The data obtained from this study will enhance our understanding of the gene network regulating quinoa inflorescence architecture, as well will supply with valuable genetic resources for high-yield elite breeding in the future.