Abstract
Aluminum toxicity is a major abiotic stress on acidic soils, leading to restricted root growth and reduced plant yield. Long non-coding RNAs are crucial signaling molecules regulating the expression of downstream genes, particularly under abiotic stress conditions. However, the extent to which lncRNAs participate in the response to aluminum (Al) stress in barley remains largely unknown. Here, we conducted RNA sequencing of root samples under aluminum stress and compared the lncRNA transcriptomes of two Tibetan wild barley genotypes, XZ16 (Al-tolerant) and XZ61 (Al-sensitive), as well as the aluminum-tolerant cultivar Dayton. In total, 268 lncRNAs were identified as aluminum-responsive genes on the basis of their differential expression profiles under aluminum treatment. Through target gene prediction analysis, we identified 938 candidate lncRNA-messenger RNA (mRNA) pairs that function in a cis-acting manner. Subsequently, enrichment analysis showed that the genes targeted by aluminum-responsive lncRNAs were involved in diterpenoid biosynthesis, peroxisome function, and starch/sucrose metabolism. Further analysis of genotype differences in the transcriptome led to the identification of 15 aluminum-responsive lncRNAs specifically altered by aluminum stress in XZ16. The RNA sequencing data were further validated by RT-qPCR. The functional roles of lncRNA-mRNA interactions demonstrated that these lncRNAs are involved in the signal transduction of secondary messengers, and a disease resistance protein, such as RPP13-like protein 4, is probably involved in aluminum tolerance in XZ16. The current findings significantly contribute to our understanding of the regulatory roles of lncRNAs in aluminum tolerance and extend our knowledge of their importance in plant responses to aluminum stress.