Transcriptomic Analysis Reveals the Role of Long Non-Coding RNAs in Response to Drought Stress in Tibetan Hulless Barley

LncRNAs, a type of RNAs exceeding 200 nucleotides (nt) and lacking representative open reading frames (ORFs), have emerged as crucial regulatory molecules that modulate numerous growth and development processes in plants. While substantial progress has been made in interpreting the functions and regulatory mechanisms of coding RNAs, the study of lncRNAs in Tibetan hulless barley remains incomplete. To elucidate the coordination of drought stress responses in Tibetan hulless barely by lncRNAs, we analyzed the previously published RNA-seq data from two cultivars of hulless barley, drought-tolerant (Z772) and drought-sensitive (Z013), subjected to varying durations of drought treatment (0, 1, and 5 h). Initially, we identified a total of 2877 lncRNAs through a strict pipeline, of which 2179 were co-expressed in both cultivars. Additionally, 331 and 367 lncRNAs showed cultivar-specific expression patterns in Z772 and Z013, respectively. Given the trans-regulatory functions of lncRNAs, we utilized WGCNA and uncovered 11 modules that were enriched in drought-responsive pathways. Within these modules, lncRNAs and neighboring PCGs were co-clustered in key control modules. The GO enrichment analysis of potential lncRNA-PCG pairs primarily involved processes related to the response to water deprivation, regulation of abiotic stress, and RNA metabolic processes. Notably, 12 high-confidence lncRNA-PCG pairs displayed concordant expression profiles, with some annotated as TFs. Two of these pairs were validated by qRT-PCR in the Tibetan hulless barley cultivar Kunlun 14. These findings suggested that lncRNAs may participate in regulatory networks involved in drought adaptation in Tibetan hulless barley, offering novel insights into the drought resistance mechanisms of Poaceae crops and potential targets for breeding drought-resistant varieties.