Identification of Salt-Stress-Responding Genes by Weighted Gene Correlation Network Analysis and Association Analysis in Wheat Leaves

The leaf is not only the main site of photosynthesis, but also an important organ reflecting plant salt tolerance. Discovery of salt-stress-responding genes in the leaf is of great significance for the molecular improvement of salt tolerance in wheat varieties. In this study, transcriptome sequencing was conducted on the leaves of salt-tolerant wheat germplasm CH7034 seedlings at 0, 1, 6, 24, and 48 h after NaCl treatment. Based on weighted gene correlation network analysis of differentially expressed genes (DEGs) under salt stress, 12 co-expression modules were obtained, of which, 9 modules containing 4029 DEGs were related to the salt stress time-course. These DEGs were submitted to the Wheat Union database, and a total of 904,588 SNPs were retrieved from 114 wheat germplasms, distributed on 21 wheat chromosomes. Using the R language package and GAPIT program, association analysis was performed between 904,588 SNPs and leaf salt injury index of 114 wheat germplasms. The results showed that 30 single nucleotide polymorphisms (SNPs) from 15 DEGs were associated with salt tolerance. Then, nine candidate genes, including four genes (TaBAM, TaPGDH, TaGluTR, and TaAAP) encoding enzymes as well as five genes (TaB12D, TaS40, TaPPR, TaJAZ, and TaWRKY) encoding functional proteins, were identified by converting salt tolerance-related SNPs into Kompetitive Allele-Specifc PCR (KASP) markers for validation. Finally, interaction network prediction was performed on TaBAM and TaAAP, both belonging to the Turquoise module. Our results will contribute to a further understanding of the salt stress response mechanism in plant leaves and provide candidate genes and molecular markers for improving salt-tolerant wheat varieties.