Patterns and impacting factors of gene evolutionary rate between wild and cultivated emmer wheat (Triticum turgidum)

Tetraploid emmer wheat (Triticum turgidum L., BBAA) is the founder progenitor of bread wheat, providing the valuable genetic resource and gene pool for wheat improvement. However, the evolutionary trajectory of tetraploid wheat, especially the evolutionary fate of different types of genes has not been well studied. In this study, the rate of non-synonymous substitution (d_N) and synonymous substitution (d_S) was calculated by comparing the orthologs between the wild emmer and cultivated durum wheat at the whole genome and subgenome levels to obtain the positively selected genes (PSGs) and negatively selected genes (NSGs). Then, mutation rate, gene length, exon number, GC content, codon bias, and expression level were comprehensively investigated and compared between the PSGs and NSGs. Within both wild emmer and cultivated durum wheat, PSGs between A and B subgenome displayed shorter gene and exon lengths as well as fewer exon numbers compared with NSGs, whereas from wild emmer to cultivated durum wheat, PSGs showed longer gene length and more exon numbers. Furthermore, PSGs displayed much higher expression levels and stronger codon usage bias, but lower genetic diversity compared with NSGs. Finally, two PSGs TdER1-6B, and TdLC7-2A, were found to play the crucial roles in regulating grain width and plant height of tetraploid wheat, respectively. This study systematically investigated the evolutionary, structural, and functional difference between PSGs and NSGs in tetraploid wheat, which will contribute to a better understanding of the selective mode and evolutionary trajectory during wheat domestication and evolution.