Comparative analysis of developing grain transcriptome reveals candidate genes and pathways improving GPC in wheat lines derived from wild emmer

The grain protein content (GPC) in modern wheat is inherently low. Wild emmer wheat (Triticum turgidumssp.dicoccoides, 2n = 4x = 28, AABB) gene pool harbors wide genotypic variations in GPC. However, the characterization of candidate genes associated with high GPC is a challenge due to the complex characteristic of this trait. In the current study, we performed RNA-seq analysis on developing grains of wild emmer genotype D1, common wheat CN16, and their hexaploid wide hybrid BAd107-4 with contrasting GPC. We have found a total of 39,795 expressed genes on chromosomes A and B, of which 24,152 were shared between D1, CN16, and BAd107-4. From 1744 differentially expressed genes (DEGs), 1203 were downregulated and 541 were upregulated in the high GPC (D1+BAd107-4) compared with low GPC (CN16) groups. The majority of DEGs were associated with protein processing in endoplasmic reticulum, starch and sucrose metabolism, galactose metabolism, and protein export pathways. Expression levels of nine randomly selected genes were verified by qRT-PCR, which was consistent with the transcriptome data. The present database will help us to understand the potential regulation networks underlying wheat grain protein accumulation and provide the foundation for simultaneous improvement of grain protein content and yield in wheat breeding programs.