In total, 519 (486 conserved and 33 novel) miRNAs were identified, of which, 41 miRNAs had significant differential expression among the comparisons (p < 0.05). They also identified 55,366 unigenes by RNA-Seq, where 12,535 unigenes were differently expressed. Finally,degradome analysis revealed that 5950 transcripts were targeted by 487 miRNAs. A correlation analysis identified that miRNA ata-miR164c-3p and its target heat shock protein family A (HSP70) member 5 gene comp59407_c0 (BIPE3) may be essential in organ-specific plant drought stress response and/or adaptation in orchardgrass. Additionally, Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analyses found that “antigen processing and presentation” was the most enriched downregulated pathway in adaptation to drought conditions.
Taken together, these results serve as a valuable genetic resource for future studies focusing on how plants adapted to drought conditions.
Differentially expressed genes (DEGs) in orchardgrass treated with drought stress. (A) Venn diagram of DEGs from root and leaf after drought treatment; pink color represent to comparison of drought treat root to controlled root, light-yellow color represent to comparison of drought treat leaf to controlled leaf; (B) column graph of DEGs in different treatment groups showing both upregulated and downregulated genes; (C) gene ontology annotation for all 5272 DEGs; (D) Kyoto encyclopedia of genes and genomes pathway annotation for all 5272 DEGs.