Global identification and integrated analysis of heat-responsive long non-coding RNAs in contrasting rice cultivars.

key messageIdentified 2743 rice lncRNAs LncRNAs in response to heat stress Function prediction of HRLs Network among HRLs, genes and miRNAs co-localization of HRLs with QTLs Significant motifs in HRL sequences Long non-coding RNAs (lncRNAs) play vital roles in plant responses to environmental challenges. A better understanding of the gene regulation mediated by lncRNAs and their systematic identification would provide great benefits for modern agriculture. In this study, we performed strand-specific RNA sequencing for two rice varieties, heat-tolerant ZS97B and heat-susceptible SYD2 under heat stress. In total, 2743 putative lncRNAs were identified, and their expression profiles in response to heat treatments were established. We identified 231 differentially expressed lncRNAs (DELs) under heat stress, including 31 DELs common to both varieties and 103 and 97 specific to ZS97B and SYD2, respectively, all defined as heat-responsive lncRNAs (HRLs). The target-coding genes of HRLs were predicted, and GO and KEGG annotations of HRL targets revealed functions in which HRLs might be involved. The interaction network between HRLs, target genes and relevant miRNAs was constructed. The HRLs and their targets were compared with publicly available QTLs for rice seedling growth under heat stimulus. Ten HRLs and twelve target genes were linked with five heat stress-relevant QTLs. Sequence analysis revealed several motifs significantly enriched within the 231 HRL sequences. Our findings provide a valuable resource for further characterization of lncRNAs in terms of heat response and plant heat tolerance improvement.