ABSTRACT: Waterlogging stress (WS) in a dynamic environment seriously limits plant growth, development, and yield. The regulatory mechanism underlying WS conditions at an early stage in maize seedlings is largely unknown. In the present study, the primary root tips of B73 seedlings were sampled before (0 h) and after (2 h, 4 h, 6 h, 8 h, 10 h, and 12 h) WS and then subjected to transcriptome sequencing, resulting in the identification of differentially expressed protein-coding genes (DEpcGs) and long non-coding RNAs (DElncRs) in response to WS. These DEpcGs were classified into nine clusters, which were significantly enriched in several metabolic pathways, such as glycolysis and methionine metabolism. Several Transcription factor families, including AP2-EREBP, bZIP, NAC, bHLH, and MYB, were also significantly enriched. In total, 6099 lncRNAs were identified, of which 3190 were DElncRs. A co-expression analysis revealed lncRNAs to be involved in 11 transcription modules, 10 of which were significantly associated with WS. The DEpcGs in the four modules were enriched in the hypoxia response pathways, including phenylpropanoid biosynthesis, MAPK signaling, and carotenoid biosynthesis, in which 137 DElncRs were also co-expressed. Most of the co-expressed DElncRs were co-localized with previously identified quantitative trait loci associated with waterlogging tolerance. A quantitative RT-PCR analysis of DEpcG and DElncR expression among the 32 maize genotypes after 4 h of WS verified significant expression correlations between them as well as the significant correlation with phenotype of waterlogging tolerance. Moreover, the high proportion of hypoxia response elements in the promoter region increased the reliability of the DElncRs identified in this study. These results provide a comprehensive transcriptome in response to WS at an early stage of maize seedlings and expand our understanding of the regulatory network involved in hypoxia in plants.
