Regulatory annotation identifies KAN1, MYB44 and PIF4 as regulators of Arabidopsis lincRNAs expressed in root

Thousands of long intergenic noncoding RNAs (lincRNAs) have been identified in plant genomes. While some lincRNAs have been characterized as important regulators in different biological processes, little is known about the transcriptional regulation for most plant lincRNAs. Through the integration of eight annotation resources, we defined 6,599 high-confidence lincRNA loci in Arabidopsis thaliana. For lincRNAs belonging to different evolutionary age categories, we identified major differences in sequence and chromatin features, as well as in the level of conservation and purifying selection acting during evolution. Spatiotemporal gene expression profiles combined with transcription factor (TF) chromatin immunoprecipitation data were used to construct a TF- lincRNA regulatory network containing 2,659 lincRNAs and 15,686 interactions. We found that properties characterizing lincRNA expression, conservation and regulation differ between plants and animals. Experimental validation confirmed the role of three TFs, KAN1, MYB44, and PIF4, as key regulators controlling root- specific lincRNA expression, demonstrating the predictive power of our network. Furthermore, we identified 58 lincRNAs, regulated by these TFs, showing strong root cell-type specific expression or chromatin accessibility, which are linked with GWAS genetic associations related to root system development and growth. The multi-level genome-wide characterization covering chromatin state information, promoter conservation, and ChIP-based TF binding, for all detectable lincRNAs across 769 expression samples, permits to rapidly define the biological context and relevance of lincRNAs in Arabidopsis regulatory networks. One-line sentence A multi-level Arabidopsis gene regulatory network identifies novel regulators controlling root-specific lincRNA expression, offering a promising strategy to identify lincRNAs involved in plant biology.