Integrated small RNA and mRNA expression profiles reveal miRNAs and their target genes in response to Aspergillus flavus growth in peanut seeds

Abstract Background MicroRNAs are important gene expression regulators in plants immune system. Aspergillus flavus is the most common causal agents of aflatoxin contamination in peanuts, but information on the function of miRNA in peanut-A. flavus interaction is lacking. In this study, the resistant cultivar (GT-C20) and susceptible cultivar (Tifrunner) were used to investigate regulatory roles of miRNAs in response to A. flavus growth. Results A total of 30 miRNAs, 447 genes and 21 potential miRNA/mRNA pairs were differentially expressed significantly when treated with A. flavus. A total of 62 miRNAs, 451 genes and 44 potential miRNA/mRNA pairs exhibited differential expression profiles between two peanut varieties. Gene Ontology (GO) analysis showed that metabolic-process related GO terms were enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses further supported the GO results, in which many enriched pathways were related with biosynthesis and metabolism, such as biosynthesis of secondary metabolites and metabolic pathways. Correlation analysis of small RNA, transcriptome and degradome indicated that miR156/SPL pairs might regulate the accumulation of flavonoids in resistant and susceptible genotypes. The miR482/2118 family might regulate NBS-LRR gene which had the higher expression level in resistant genotype. These results provided useful information for further understanding the roles of miR156/157/SPL and miR482/2118/NBS-LRR pairs. Conclusions Integration analysis of the transcriptome, miRNAome and degradome of resistant and susceptible peanut varieties were performed in this study. The knowledge gained will help to understand the roles of miRNAs of peanut in response to A. flavus.