Integration of mRNA and miRNA analysis reveals the molecular mechanisms of sugar beet (Beta vulgaris L.) response to salt stress

The continuous increase of saline-alkali areas worldwide has led to the emergence of saline-alkali conditions, which are the primary abiotic stress or hindering the growth of plants. Beet is among the main sources of sugar, and its yield and sugar content are notably affected by saline-alkali stress. Despite sugar beet being known as a salt-tolerant crop, there are few studies on the mechanisms underlying its salt tolerance, and previous studies have mainly delineated the crop’s response to stress induced by NaCl. Recently, advancements in miRNA-mRNA network analysis have led to an increased understanding of how plants, including sugar beet, respond to stress. In this study, seedlings of beet variety "N98122" were grown in the laboratory using hydroponics culture and were exposed to salt stress at 40 days of growth. According to the phenotypic adaptation of the seedlings' leaves from a state of turgidity to wilting and then back to turgidity before and after exposure, 18 different time points were selected to collect samples for analysis. Subsequently, based on the data of real-time quantitative PCR (qRT-PCR) of salt-responsive genes, the samples collected at the 0, 2.5, 7.5, and 16 h time points were subjected to further analysis with experimental materials. Next, mRNA-seq data led to the identification of 8455 differentially expressed mRNAs (DEMs) under exposure to salt stress. In addition, miRNA-seq based investigation retrieved 3558 miRNAs under exposure to salt stress, encompassing 887 known miRNAs belonging to 783 families and 2,671 novel miRNAs. With the integrated analysis of miRNA-mRNA network, 57 miRNA-target gene pairs were obtained, consisting of 55 DEMIs and 57 DEMs. Afterwards, we determined the pivotal involvement of aldh2b7, thic, and δ-oat genes in the response of sugar beet to the effect of salt stress. Subsequently, we identified the miRNAs novel-m035-5p and novel-m0365-5p regulating the aldh gene and miRNA novel-m0979-3p regulating the thic gene. The findings of miRNA and mRNA expression were validated by qRT-PCR.