Abstract
Maize (Zea mays), a major food crop worldwide, is susceptible to infection by the saprophytic fungus Aspergillus flavus that can produce the carcinogenic metabolite aflatoxin (AF) especially under climate change induced abiotic stressors that favor mold growth. Several studies have used "-omics" approaches to identify genetic elements with potential roles in AF resistance, but there is a lack of research identifying the involvement of small RNAs such as microRNAs (miRNAs) in maize-A. flavus interaction. In this study, we compared the miRNA profiles of three maize lines (resistant TZAR102, moderately resistant MI82, and susceptible Va35) at 8 h, 3 d, and 7 d after A. flavus infection to investigate possible regulatory antifungal role of miRNAs. A total of 316 miRNAs (275 known and 41 putative novel) belonging to 115 miRNA families were identified in response to the fungal infection across all three maize lines. Eighty-two unique miRNAs were significantly differentially expressed with 39 miRNAs exhibiting temporal differential regulation irrespective of the maize genotype, which targeted 544 genes (mRNAs) involved in diverse molecular functions. The two most notable biological processes involved in plant immunity, namely cellular responses to oxidative stress (GO:00345990) and reactive oxygen species (GO:0034614) were significantly enriched in the resistant line TZAR102. Coexpression network analysis identified 34 hubs of miRNA-mRNA pairs where nine hubs had a node in the module connected to their target gene with potentially important roles in resistance/susceptible response of maize to A. flavus. The miRNA hubs in resistance modules (TZAR102 and MI82) were mostly connected to transcription factors and protein kinases. Specifically, the module of miRNA zma-miR156b-nb - squamosa promoter binding protein (SBP), zma-miR398a-3p - SKIP5, and zma-miR394a-5p - F-box protein 6 combinations in the resistance-associated modules were considered important candidates for future functional studies.