Integration of Dual Stress Transcriptomes and Major QTLs from a Pair of Genotypes Contrasting for Drought and Chronic Nitrogen Starvation Identifies Key Stress Responsive Genes in Rice

We report here the genome-wide changes resulting from low N (N-W+), low water (N+W-)) and dual stresses (N-W-) in root and shoot tissues of two rice genotypes, namely, IR 64 (IR64) and Nagina 22 (N22), and their association with the QTLs for nitrogen use efficiency. For all the root parameters, except for root length under N-W+, N22 performed better than IR64. Chlorophyll a, b and carotenoid content were higher in IR64 under N+W+ treatment and N-W+ and N+W- stresses; however, under dual stress, N22 had higher chlorophyll b content. While nitrite reductase, glutamate synthase (GS) and citrate synthase assays showed better specific activity in IR64, glutamate dehydrogenase showed better specific activity in N22 under dual stress (N-W-); the other N and C assimilating enzymes showed similar but low specific activities in both the genotypes. A total of 8926 differentially expressed genes (DEGs) were identified compared to optimal (N+W+) condition from across all treatments. While 1174, 698 and 903 DEGs in IR64 roots and 1197, 187 and 781 in N22 roots were identified, nearly double the number of DEGs were found in the shoot tissues; 3357, 1006 and 4005 in IR64 and 4004, 990 and 2143 in N22, under N-W+, N+W- and N-W- treatments, respectively. IR64 and N22 showed differential expression in 15 and 11 N-transporter genes respectively, under one or more stress treatments, out of which four showed differential expression also in N+W- condition. The negative regulators of N- stress, e.g., NIGT1, OsACTPK1 and OsBT were downregulated in IR64 while in N22, OsBT was not downregulated. Overall, N22 performed better under dual stress conditions owing to its better root architecture, chlorophyll and porphyrin synthesis and oxidative stress management. We identified 12 QTLs for seed and straw N content using 253 recombinant inbred lines derived from IR64 and N22 and a 5K SNP array. The QTL hotspot region on chromosome 6 comprised of 61 genes, of which, five were DEGs encoding for UDP-glucuronosyltransferase, serine threonine kinase, anthocyanidin 3-O-glucosyltransferase, and nitrate induced proteins. The DEGs, QTLs and candidate genes reported in this study can serve as a major resource for both rice improvement and functional biology.