Combined RNA-seq and Phenotype Analysis Reveals a Potential Molecular Mechanism of the Difference in Grain Size of Naked Barley From the Qinghai-Tibetan Plateau

To understand the molecular mechanism controlling the size of barley grains, a number of traits were analyzed and RNA-seq was conducted on grains of two barley materials with a significant difference in thousand-grain weight (TGW) after flowering. The trait dataset delineates the dynamic changes in grain size after flowering, and it provides an understanding of the source of the difference in TGW. By comparing the transcripts of barley grains at several stages after flowering, we identified the gene expression characteristics and significantly enriched pathways in each stage. At the early stage of grain development, genes involved in fatty acid metabolism, plant hormone signal transduction, and pathways involved in cytoskeleton formation were significantly upregulated. At the later stage of grain development, genes involved in starch synthesis, glucose metabolism, and other pathways were significantly upregulated. Further, we used weighted gene coexpression network analysis (WGCNA) and correlation analysis of trait datasets to identify the coexpressed gene modules significantly associated with traits, such as grain length (GL), grain width (GW), and dry weight (DW). After comparing the modules with the differentially expressed gene (DEG) set, 12 candidate genes were selected, and among these, four genes were homologous to genes that regulate grain size in rice and other plants. The combined analysis identified many potential key regulatory factors that may control barley grain size and yield potential, thus providing new insights into the molecular mechanism of barley grain size.