Abstract
Heterosis is a complex biological phenomenon regulated by genetic variations and epigenetic changes. However, the roles of small RNAs (sRNAs), an important epigenetic regulatory element, on plant heterosis are still poorly understood. Here, an integrative analysis was performed with sequencing data from multi-omics layers of maize hybrids and their two homologous parental lines to explore the potential underlying mechanisms of sRNAs in plant height (PH) heterosis. sRNAome analysis revealed that 59 (18.61%) microRNAs (miRNAs) and 64,534 (54.00%) 24-nt small interfering RNAs (siRNAs) clusters were non-additively expressed in hybrids. Transcriptome profiles showed that these non-additively expressed miRNAs regulated PH heterosis through activating genes involved in vegetative growth-related pathways while suppressing those related to reproductive and stress response pathways. DNA methylome profiles showed that non-additive methylation events were more likely to be induced by non-additively expressed siRNA clusters. Genes associated with low-parental expression (LPE) siRNAs and trans-chromosomal demethylation (TCdM) events were enriched in developmental processes as well as nutrients and energy metabolism, whereas genes associated with high-parental expression (HPE) siRNAs and trans-chromosomal methylation (TCM) events were gathered in stress response and organelle organization pathways. Our results provide insights into the expression and regulation patterns of sRNAs in hybrids and help to elucidate their potential targeting pathways contributing to PH heterosis.