Population-scale gene expression analysis reveals the contribution of expression diversity to the modern wheat improvement

Changes in gene expression are crucial for crop breeding, yet population genomics has primarily focused on sequence polymorphisms rather than gene expression diversity. The strategy of using single genome reference for RNA-seq analysis could not handle introgression bias, especially for hexaploidy wheat. Here, we conducted RNA-seq for 328 wheat lines, including representative diverse landraces and elite cultivars from China and the United States, to investigate the role of gene expression variation in shaping agronomic traits. Using pan-genome resources, we identified 20,615 more transcripts than using the ‘Chinese Spring’ reference genome alone. We constructed a pan-gene atlas regulatory map through eQTL analysis, demonstrating that genes introgressed from wild relatives were under tight genetic control. Genes responding to environmental stress show higher activity after introgressed into the wheat genome, demonstrating how long-term breeding selection impacted the gene expression regulation of targeted introgression. Multi-omics modeling identified 231 high-confidence candidate genes for 34 field agronomic traits and the seedling resistance phenotypes of 8 powdery mildew isolates. More than one fifth of those candidates have no homolog in ‘Chinese Spring’ reference genome. By utilizing the indexed KN9204 EMS library, 80% candidates showed significant trait difference between wild type and mutant lines. Furthermore, directional shifts in genes of which expression were changed by breeding improvement demonstrated distinct adaptations to local environments. Our study constructed a pan-gene atlas to correct the reference bias of reads mapping in RNA-seq studies and revealed the expression patterns of introgressed genes within the wheat genome and their regulatory mechanisms, which highlighted the impact of breeding selection on gene expression of the world’s most important crop.