Transcriptome Analysis and miRNA Target Profiling at Various Stages of Root-Knot Nematode Meloidogyne incognita Development for Identification of Potential Regulatory Networks

. The egg stage showed the maximum number of transcripts, and 12,803 gene transcripts were expressed in all stages. Functional Gene Ontology (GO) analysis resulted in three main GO classes: biological process, cellular components, and molecular function; the detected sequences were longer than sequences in the reference genome. Stage-wise selected fragments per kilobase of transcript per million mapped reads (FPKM) values of the top 10 stage-specific and common mRNAs were used to construct expression profiles, and 20 mRNAs were validated through quantitative real-time PCR (qRT-PCR). Next, we used three target prediction programs (miRanda, RNAhybrid, and PITA) to obtain 2431 potential target miRNA genes in RKNs, which regulate 8331 mRNAs. The predicted potential targets of miRNA were generally involved in cellular and metabolic processes, binding of molecules in the cell, membranes, and organelles. Stage-wise miRNA target analysis revealed that the egg stage contains heat shock proteins, transcriptional factors, and DNA repair proteins, whereas J2 includes DNA replication, heat shock, and ubiquitin-conjugating pathway-related proteins; the J3 and J4 stages are represented by the major sperm protein domain and translation-related proteins, respectively. In the female stage, we found proteins related to the maintenance of molybdopterin-binding domain-containing proteins and ubiquitin-mediated protein degradation. In total, 29 highly expressed stage-specific mRNA-targeting miRNAs were analyzed using qRT-PCR to validate the sequence analysis data. Overall, our findings provide new insights into the molecular mechanisms occurring at various developmental stages of the RKN life cycle, thus aiding in the identification of potential control strategies.