Drought-induced circular RNAs in maize roots: separating signal from noise.

Circular RNAs (CircRNAs) play an important role in diverse biological processes; however, their origin and functions, especially in plants, remain largely unclear. Here, we used two maize (Zea mays) inbred lines, as well as 14 of their derivative RILs with different drought sensitivity, to systematically characterize 8,790 circRNAs in maize roots under well-watered (WW) and water-stress (WS) conditions. We found that a diverse set of circRNAs expressed at significantly higher levels under WS. Enhanced expression of circRNAs was associated with longer flanking introns and an enrichment of long interspersed nuclear element (LINE) retrotransposable elements. The epigenetic marks found at the back-splicing junctions of circRNA-producing genes were markedly different from canonical splicing, characterized by increased levels of H3K36me3/H3K4me1, as well as decreased levels of H3K9Ac/H3K27Ac. We found that genes expressing circRNAs are subject to relaxed selection. The significant enrichment of trait-associated sites along their genic regions suggested that genes giving rise to circRNAs were associated with plant survival rate under drought stress, implying that circRNAs play roles in plant drought responses. Furthermore, we found that overexpression of circMED16, one of the drought-responsive circRNAs, enhances drought tolerance in Arabidopsis (Arabidopsis thaliana). Our results provide a framework for understanding the intricate interplay of epigenetic modifications and how they contribute to the fine-tuning of circRNA expression under drought stress.

Identification and Functional Analysis of Drought-Responsive Long Noncoding RNAs in Maize Roots.

Long noncoding RNAs (lncRNAs) are transcripts with lengths of more than 200 nt and limited protein-coding potential. They were found to play important roles in plant stress responses. In this study, the maize drought-tolerant inbred line AC7643 and drought-sensitive inbred line AC7729/TZSRW, as well as their recombinant inbred lines (RILs) were selected to identify drought-responsive lncRNAs in roots. Compared with non-responsive lncRNAs, drought-responsive lncRNAs had different sequence characteristics in length of genes and number of exons. The ratio of down-regulated lncRNAs induced by drought was significantly higher than that of coding genes; and lncRNAs were more widespread expressed in recombination sites in the RILs. Additionally, by integration of the modifications of DNA 5-methylcytidine (5mC), histones, and RNA N6-methyladenosine (m6A), it was found that the enrichment of histone modifications associated with transcriptional activation in the genes generated lncRNAs was lower that coding genes. The lncRNAs-mRNAs co-expression network, containing 15,340 coding genes and 953 lncRNAs, was constructed to investigate the molecular functions of lncRNAs. There are 13 modules found to be associated with survival rate under drought. We found nine SNPs located in lncRNAs among the modules associated with plant survival under drought. In conclusion, we revealed the characteristics of lncRNAs responding to drought in maize roots based on multiomics studies. These findings enrich our understanding of lncRNAs under drought and shed light on the complex regulatory networks that are orchestrated by the noncoding RNAs in response to drought stress.