Genome-wide transcriptional profiling provides clues to molecular mechanisms underlying cold tolerance in chickpea

Chickpea is an important food legume cultivated in several countries. A sudden drop in autumn temperature, freezing winter temperature, and late spring cold events result in significant losses in chickpea production. The current study used RNA sequencing of two cold tolerant (Saral) and sensitive (ILC533) Kabuli chickpea genotypes to identify cold tolerance-associated genes/pathways. A total of 200.85 million raw reads were acquired from the leaf samples by Illumina sequencing, and around 86% of the clean reads (199 million) were mapped to the chickpea reference genome. The results indicated that 3710 (1980 up- and 1730 down-regulated) and 3473 (1972 up- and 1501 down-regulated) genes were expressed differentially under cold stress in the tolerant and sensitive genotypes, respectively. According to the GO enrichment analysis of uniquely down-regulated genes under cold stress in ILC533, photosynthetic membrane, photosystem II, chloroplast part, and photosystem processes were enriched, revealing that the photosynthesis is severely sensitive to cold stress in this sensitive genotype. Many remarkable transcription factors (CaDREB1E, CaMYB4, CaNAC47, CaTCP4, and CaWRKY33), signaling/regulatory genes (CaCDPK4, CaPP2C6, CaMKK2, and CaHSFA3), and protective genes (CaCOR47, CaLEA3, and CaGST) were identified among the cold-responsive genes of the tolerant genotype. These findings would help improve cold tolerance across chickpea genotypes by molecular breeding or genetic engineering.