Exploring the genetic potential of Pakistani soybean cultivars through RNA-seq based transcriptome analysis.

BACKGROUND: Soybean is largely grown and considered among the top oilseed crops. Three Pakistani cultivars, NARC-II (N), Swat-84 (S), and Rawal-I (R) were employed for RNA-Seq based transcriptome analysis to explore their genetic potential and performance in our local environment. METHODS AND RESULTS: We grew the plants in glass house at same conditions and sampled leaves for RNA-Seq analysis in triplicate for each variety. We retrieved 2225 differentially expressed genes (DEGs) between S vs R, 2591 DEGs between S vs N, and 1221 DEGs between R vs N cultvars. These genes consist of transcription factors representing Basic Helix-loop Helix, myeloblastosis, ethylene response factors, and WRKY amino acid motif (WRKY) type major families that were up-regulated. KEGG pathway analysis revealed that MAPK, plant hormone signal transduction, and Phenylpropanoid biosynthesis pathways were the most dominant pathways involved in plant defense and growth. Comparative analysis showed that Swat-84 (S) cultivar had better gene expression among these varieties having higher number of DEGs, where mostly genes related to important phenotypic traits were up regulated. CONCLUSIONS: This is a pilot study to investigate and functionally characterise the DEG involved in the stress response in the cultivars studied.