Novel low-nitrogen stress-responsive long non-coding RNAs (lncRNA) in barley landrace B968 (Liuzhutouzidamai) at seedling stage

Reducing the dependence of crop production on chemical fertilizer with its associated costs, carbon footprint and other environmental problems is a challenge for agriculture. New solutions are required to solve this problem, and crop breeding for high nitrogen use efficiency or tolerance of low nitrogen availability has been widely considered to be a promising approach. However, the molecular mechanisms of high nitrogen use efficiency or low-nitrogen tolerance in crop plants are still to be elucidated, including the role of long non-coding RNAs (lncRNAs).

In this study, we identified 498 lncRNAs in barley (Hordeum vulgare) landrace B968 (Liuzhutouzidamai), of which 487 were novel, and characterised 56 that were responsive to low-nitrogen stress. For functional analysis of differentially-expressed lncRNAs, the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of co-expressed and co-located protein-coding genes were analyzed, and interactions with annotated co-expressed protein coding genes or micro RNAs (miRNAs) were further predicted. Target mimicry prediction between differentially-expressed lncRNAs and miRNAs identified 40 putative target mimics of lncRNAs and 58 target miRNAs. Six differentially-expressed lncRNAs were further validated by qPCR, and one in particular showed consistent differential expression using both techniques. Expression levels of most of the lncRNAs were found to be very low, and this may be the reason for the apparent inconsistency between RNA-seq and qPCR data.

The analysis of lncRNAs that are differentially-expressed under low-nitrogen stress, as well as their co-expressed or co-located protein coding genes and target mimics, could elucidate complex and hitherto uncharacterised mechanisms involved in the adaptation to low-nitrogen stress in barley and other crop plants.