Maize is an important cereal crop worldwide and is sensitive to abiotic stresses in fluctuant environments that seriously affect its growth, yield, and quality. The small heat shock protein (HSP20) plays a crucial role in protecting plants from abiotic stress. However, little is known about HSP20 in maize (ZmHSP20). In this study, 44 ZmHSP20s were identified, which were unequally distributed over 10 chromosomes, and 6 pairs of ZmHSP20s were tandemly presented. The gene structure of ZmHSP20s was highly conserved, with 95% (42) of the genes having no more than one intron. The analysis of the cis-element in ZmHSP20s promoter demonstrated large amounts of elements related to hormonal and abiotic stress responses, including abscisic acid (ABA), high temperature, and hypoxia. The ZmHSP20s protein had more than two conserved motifs that were predictably localized in the cytoplasm, nucleus, endoplasmic reticulum, peroxisome, mitochondria, and plasma. Phylogenetic analysis using HSP20s in Arabidopsis, rice, maize, and Solanum tuberosum indicated that ZmHSP20s were classified into 11 categories, of which each category had unique subcellular localization. Approximately 80% (35) of ZmHSP20 were upregulated under heat stress at the maize seedling stage, whereas the opposite expression profiling of 10 genes under 37 and 48 °C was detected. A total of 20 genes were randomly selected to investigate their expression under treatments of ABA, gibberellin (GA), ethylene, low temperature, drought, and waterlogging, and the results displayed that more than half of these genes were downregulated while ZmHSP20-3, ZmHSP20-7, ZmHSP20-24, and ZmHSP20-44 were upregulated under 1 h treatment of ethylene. A yeast-one-hybrid experiment was conducted to analyze the binding of four heat stress transcription factors (ZmHSFs) with eight of the ZmHSP20s promoter sequences, in which ZmHSF3, ZmHSF13, and ZmHSF17 can bind to most of these selected ZmHSP20s promoters. Our results provided a valuable resource for studying HSP20s function and offering candidates for genetic improvement under abiotic stress.

Transcriptome analysis based on high-throughput sequencing of a cDNA library has been widely applied to functional genomic studies. However, the cDNA dependence of most RNA sequencing techniques constrains their ability to detect base modifications on RNA, which is an important element for the post-transcriptional regulation of gene expression. To comprehensively profile the N6-methyladenosine (m6A) and N5-methylcytosine (m5C) modifications on RNA, direct RNA sequencing (DRS) using the latest Oxford Nanopore Technology was applied to analyze the transcriptome of six tissues in rice. Approximately 94 million reads were generated, with an average length ranging from 619 nt to 1013 nt, and a total of 45,707 transcripts across 34,763 genes were detected. Expression profiles of transcripts at the isoform level were quantified among tissues. Transcriptome-wide mapping of m6A and m5C demonstrated that both modifications exhibited tissue-specific characteristics. The transcripts with m6A modifications tended to be modified by m5C, and the transcripts with modifications presented higher expression levels along with shorter poly(A) tails than transcripts without modifications, suggesting the complexity of gene expression regulation. Gene Ontology analysis demonstrated that m6A- and m5C-modified transcripts were involved in central metabolic pathways related to the life cycle, with modifications on the target genes selected in a tissue-specific manner. Furthermore, most modified sites were located within quantitative trait loci that control important agronomic traits, highlighting the value of cloning functional loci. The results provide new insights into the expression regulation complexity and data resource of the transcriptome and epitranscriptome, improving our understanding of the rice genome.