Next Article in Journal
Evaluation of Methane Emission Reduction Potential of Water Management and Chinese Milk Vetch Planting in Hunan Paddy Rice Fields
Next Article in Special Issue
Boron Nutrition in Coffee Improves Drought Stress Resistance and, Together with Calcium, Improves Long-Term Productivity and Seed Composition
Previous Article in Journal
Growth and Physiological Characteristics of Sour Jujube Seedlings in Different Substrate Formulations
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agronomy
Volume 13
Issue 7
10.3390/agronomy13071798
agronomy-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Function of Soybean miR159 Family Members in Plant Responses to Low Phosphorus, High Salinity, and Abscisic Acid Treatment

by
Bodi Li
1,†,
Ping Tao
1,2,†,
Feng Xu
1,
Pingan He
1 and
Jinxiang Wang
1,3,*
1
College of Natural Resources and Environment & Root Biology Center, South China Agricultural University, Guangzhou 510642, China
2
Institute of Biological and Medical Engineering, Guangdong Academy of Sciences, Guangzhou 510316, China
3
Key Laboratory of Agricultural and Rural Pollution Control and Environmental Safety in Guangdong Province, Guangzhou 510642, China
*
Author to whom correspondence should be addressed.
†
These authors equally contributed to this work.
Agronomy 2023, 13(7), 1798; https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071798
Submission received: 1 June 2023 / Revised: 21 June 2023 / Accepted: 3 July 2023 / Published: 5 July 2023
(This article belongs to the Special Issue Effect of Nutrient Deficiencies on Stress Tolerance of Plants and Its Mechanism)
Browse Figures
Versions Notes

Abstract

:
MicroRNAs (miRNAs) regulate plant growth and development and plant responses to biotic and abiotic stresses. Although extensive studies show that miR159 family members regulate leaf and flower development in Arabidopsis thaliana, the roles of miRNAs in soybean (Glycine max) are poorly understood. Here, we identified six MIR159 genes in soybean, MIR159a–MIR159f, and investigate their expression patterns in plants under low-phosphorus (low-P), NaCl, or abscisic acid (ABA) treatments. In soybean leaves, MIR159e and MIR159f expression was induced by low-P treatment, while in roots, MIR159b, MIR159c, MIR159e, and MIR159f expression was upregulated. In flowers, low-P led to upregulation of MIR159a, MIR159b, MIR159c, and MIR159f but downregulation of MIR159d and MIR159e. In soybean nodules, MIR159b was upregulated but MIR159a, MIR159c, and MIR159d was downregulated under P deficiency. NaCl treatment induced MIR159a, MIR159b, MIR159c, and MIR159e expression in leaves and MIR159a–MIR159f expression in roots. ABA treatment upregulated MIR159a, MIR159b, and MIR159c but downregulated MIR159d, MIR159e, and MIR159f in leaves. These results suggest that miR159 family members function in plant abiotic stress responses. Moreover, total P content in leaves was significantly lower in plants overexpressing MIR159e than in the wild type, suggesting that miR159e may regulate P absorption and transport in soybean plants.

1. Introduction

MicroRNAs (miRNAs) are non-coding small (20~24 nucleotide) RNAs that are widespread in plants and animals as well as in single-celled algae [1]. MiRNAs inhibit gene expression and translation after transcription [2]. Several conserved miRNA families have been identified. In plants, for example, miR156, miR159, miR164, miR168, miR169, miR171, miR319, miR390, miR393, miR399, and miR827 play crucial roles in growth and development as well as in response to biotic and abiotic stress.
As one of conserved miRNA family, miR159s have been found in Arabidopsis [3], rice [4], soybean [5], cotton [6], tomato [7], potato [8], strawberry [9], rapeseed [10], wheat [11], apple [12], lily [13], and tobacco [14], etc. As has been well documented, the miR159 family plays important roles in plant growth and development [15]. The levels of Arabidopsis (Arabidopsis thaliana) miR159 are regulated by gibberellic acid (GA), and overexpressing miR159 reduces the transcripts of LEAFY, inhibits flowering in short-day conditions, and perturbs anther development [3]. Arabidopsis MYB101 and MYB33 are targets of miR159; these genes encode MYB transcription factors that function in seed germination [16]. Arabidopsis MiR159 is present in pollen where it has a crucial role in fertility [17,18]. In contrast to Arabidopsis wild type, the double mutant miR159a miR159b has a larger meristem; accordingly, the transcript levels of MYB33, MYB65, and MYB101 in miR159a miR159b are increased [19]. Rice miR159 positively regulates grain filling via downregulating the expression of miR167 that targets OsARF12 [4]. As previously documented, tomato miR159 (Sly-miR159) control fruit growth and development via promoting GA biosynthesis. The target gene of Sly-miR159 SlGAMYB2 directly inhibits the activities of SlGA3ox2 that catalyze the degradation of GA [7]. Overexpressing Brassica napus miR159 (bna-miR159) in Arabidopsis results in lower seed setting and shorter siliques [10].
miR159s are involved in the responses of plants to biotic stress. Soybean miR159-3p is down-regulated during cyst nematode infection; consequently, overexpressing miR159 in soybean hairy roots increases the resistance to Heterodera glycines [5]. Lily miR159a play positive roles to cope with grey mold via down-regulating its target gene GAMYB [13]. Tobacco miR159 inhibits constitutive pathogen defense [14].
miR159s play crucial roles in land plants to cope with abiotic stresses. For example, wheat miR159 is induced by osmotic stress [20]. miR159s, in particular, regulate plant responses to nutrient stress, such as low phosphate (low-P) stress. MIR159a, MIR319a, MIR396a, MIR389b, and MIR1507a are upregulated under low-P stress [21]. Nitrogen is a very important macronutrient for plant growth and development. Previous studies have revealed some low nitrate-responsive miRNAs such as miR169 and miR393 [22,23]. Moreover, ABA has been proved to regulate plants’ responses to low nitrogen [24]. On the other hand, interactions between phosphate and nitrogen have been reported [25]. However, the roles of soybean miR159s in phosphorus and nitrogen nutrition remain elusive. Eight miRNAs associated with salt tolerance have been identified in peanut (Arachis hypogea): miR159-1, miR159-2, miR159-3, miR164-2, miR167-3, miR319-1, miR319-2, and miR211-1 [26]. A study of the responses of salt tolerance related miRNAs to high-salt habitats in the mangrove companion plant Sesuvium portulacastrum showed that miR159 is involved in the responses of different tissues to high-salt stress [27], whereas the responses of soybean (Glycine max) MIR159 genes to abiotic stress and their roles in abiotic stress responses are unclear. It is well known that soil salinity and drought are common threats to crops. As global warming becomes more and more severe, the salinity and drought stress will be more common, which could lead to decreased crop yields worldwide. Moreover, the involvement of stress hormone ABA in salinity and drought has been well documented. Hence, in this study, we focus on the responses of soybean miR159 family members to low phosphate, salinity, and ABA.
The Arabidopsis hyponastic leaves 1 (hyl1) mutant, which fails to process primary miRNA transcripts (pri-miRNAs), is sensitive to abscisic acid (ABA) during germination [28], suggesting that the ABA signaling pathway might be regulated by miRNAs. Indeed, miR159 responds to ABA and drought treatment [16,29]. Analysis of the upstream regions of three MIR159 genes revealed the presence of ABA-responsive elements (ABREs) and associated stress factors such as AtMYC2 binding sites [30]. In line with this finding, ABRE-like elements were identified in the upstream region of the MIR159a promoter [16].
Here, we identified six MIR159 genes in soybean: MIR159a–MIR159f. We then explored the responses of the soybean MIR159 gene family to low-P, NaCl, and ABA treatment. Finally, we overexpressed MIR159e in soybean. Under LP conditions, the total P content of MIR159e-overexpressing transgenic plants was significantly lower than the wild type in leaves but not roots. These results suggest that miR159e modulates the absorption and transport of P in soybean.

2. Materials and Methods

2.1. Plant Materials and Growth Conditions

Seeds of soybean (Glycine max) genotype YC03-3 were sterilized with 10% NaClO, germinated in sand, and transferred to nutrient solution when the cotyledons were open and the apical bud had developed. To explore the responses of soybean miR159 family members to low-P stress, NaCl stress, and ABA treatment, seedlings with the first developed trifoliate leaf were transferred to high phosphate (Pi), high N (HP, 250 µM KH2PO4; HN, 5.3 mM nitrogen), N-deficient (HP, 250 µM KH2PO4; LN, 530 µM nitrogen), or Pi-deficient (LP, 5 µM KH2PO4; HN, 5.3 mM nitrogen; LN, 530 µM nitrogen) medium for 0 days (0 D), 7 days (7 D), and 40 days (40 D), respectively.
The plants were grown in growth chambers under a 16 hour (h) light/8 h dark cycle. The nutrient solutions were aerated for 15 min every 3 h and replaced with fresh nutrient solutions every 2 days. The roots, leaves, flowers and nodules were sampled, frozen in liquid nitrogen, and stored at −80 °C. After 0 or 6 h of ABA treatment, including −ABA (0 μM) or +ABA (300 μM), the treated seedlings were sampled, quickly frozen in liquid nitrogen, and stored at −80 °C. After 0 or 6 h of NaCl treatment, including −NaCl (0 mM) or +NaCl (200 mM), the seedlings were sampled, frozen in liquid nitrogen, and stored at −80 °C.

2.2. Extraction of Total RNA, Reverse-Transcription, and Extraction of Genomic DNA

Total RNA was extracted from soybean seedlings using the TRIzol method. The RNA was reverse transcribed to cDNA using M-MLV reverse transcriptase. DNA was extracted from the samples using the CTAB method. Other molecular experiments were performed using standard methods as described [31].

2.3. qRT-PCR Analysis

To quantify transcript levels in soybean genotype YC03-3 under different P conditions, quantitative real-time PCR (qRT-PCR) was performed. The roots and leaves of YC03-3 seedlings that were subjected to low-P stress were quickly frozen in liquid N and stored at −80 °C. Total RNA extraction and cDNA synthesis were performed using standard methods as described [31]. Reverse transcription reaction was performed using PrimeScriptTM RT regent Kit (Takara) in a final volume of 10 µL, including 2 µL of 5× PrimeScript™ Buffer, 0.5 µL of PrimeScript™ RT Enzyme Mix I, and 1 µM of Universal RT-primer, and it was incubated at 37 °C for 15 min followed by enzyme inactivation at 85 °C for 5 s. Soybean EF1a was used to normalize the PCR data. The expression levels of the six MIR159 genes (MIR159a, MIR159b, MIR159c, MIR159d, MIR159e, MIR159f) were measured in leaves, roots, and nodules. Forward and reverse primers for miRNA qRT-PCR amplification were designed according to previous studies [12,32] (Supplemental Table S2).

2.4. Measurement of Fresh Weight, Soluble Phosphate (Pi), and Total P Contents in Soybean Seedlings

Samples of soybean plants that were subjected to low-P stress for 40 days, as described, were weighed to quantify the fresh weight, including the fresh root and leaf. The lengths of roots and shoots were measured, and these values were used to calculate the root-to-shoot ratio. Soybean seedlings were dried at 105 °C for 30 min, oven-dried at 75 °C and then weighed. Soluble Pi and total P concentrations were measured exactly as described [32,33].

2.5. Vector Construction and Soybean Transformation

In short, the overexpressing soybean miR159e vector construction and whole soybean plant transformation were made based on the methods described [34]. The sequence of miR159e precursor was PCR amplified from the cDNA of soybean YC03-3 genotype using specific primers (forward: AAcccgggTAGCAAGGGTTTAGGTGGTG, reverse primer: AAtctagaAGAGCAAGAACGAGATTATGG) that contain SmaI and XbaI site, respectively. Then, the PCR product was cloned into the binary vector pTF101.1 that harbor a Bar gene under the control of 2× 35S CaMV promoter.
The cotyledonary-node was transformed. Primary transformants were grown to maturity in the greenhouse. Overexpressing miR159e lines derived from double 35S promoters were determined by PCR and Bar resistance [34]. Single copy T-DNA insertion homozygous transformant lines were used in this study.

2.6. Phytogenetic Tree Reconstruction and Data Analysis

MEGA 6.0 software [35] was used to construct the phylogenetic tree using the neighbor-joining method. The promoters were analyzed using the TSSP program of Phytozome (https://phytozome-next.jgi.doe.gov/) (accessed on 4 July 2023) and softberry (http://www.softberry.com/) (accessed on 4 July 2023); target genes were predicted using psRNA Target (http://plantgrn.noble.org/psRNATarget/) (accessed on 4 July 2023) [36].
All data were analyzed with Excel 2010. Student’s t-test was employed to identify significant differences between treatment groups. GraphPad Prism 8 was used to draw the figures.

3. Results

3.1. Identification of MIR159 Gene Family Members in Soybean

To explore the MIR159 gene family, their target genes, and the evolution of this family in soybean, we obtained sequences of miR159 family members in soybean, rice (Oryza sativa), Arabidopsis, alfalfa (Medicago truncatula), maize (Zea mays), and Brachypodium distachyon from the miRBase website (http://www.mirbase.org) (accessed on 4 July 2023). We found that the soybean MIR159 family consists of MIR159a, MIR159b, MIR159c, MIR159d, MIR159e, and MIR159f. MIR159a and MIR159d are located on chromosome 9, MIR159b and MIR159e are on chromosome 7, and MIR159c and MIR159f are on chromosome 16 (Table 1).
To decipher the evolutionary history of MIR159 in plants, we compared the MIR159 sequences of soybean with those from the plant species listed above. After aligning them using the CLUSTALW program, we analyzed their sequences using MrBayes3.2, finding that plant MIR159s belong to three subgroups (Figure 1). The soybean MIR159 genes fell into two subgroups: MIR159a, MIR159d, and MIR159e belong to subgroup I, and MIR159b, MIR159c, and MIR159f belong to subgroup II. Subgroup III contains only the MIR159 gene from monocotyledons (Figure 1). These results indicate that MIR159 genes are conserved across plant species, implying that they might play important roles in plant growth and development.
ABREs can be found in the promoter regions of MIR159 genes in Arabidopsis [16], and MIR159e responds to low-P stress in soybean [32]. To identify the cis-elements in these gene promoters, we analyzed the 2000-nucleotide sequence upstream of the start codon (ATG) in all six MIR159 genes (Table 2). As expected, we found phosphorus-responsive elements, NaCl-responsive elements, and ABREs in the promoters of all six soybean MIR159 genes (Table 2). Among the phosphorus-responsive elements, we identified three or four TATA-box and W-box binding elements in each of the six MIR159 gene promoters, but no TATA-box-like binding elements, and we found that only MIR159c and MIR159f contain P1BS binding elements. Among the NaCl-responsive elements, all six MIR159 gene promoters contain one ABRE-like and one ACGT sequence binding element, but MIR159b and MIR159c lack RD22 binding elements and MIR159f lacks AtMYB2 binding elements. Finally, among ABREs, all six MIR159 promoters contain a DPBF binding element, while MIR159c and MIR159f lack ABRE, and only MIR159e contains an RY binding element (Table 2). These results suggest that the transcription of MIR159 genes in soybean seem to be regulated by PHR1, WRKY, and ABF transcriptional regulators.

3.2. Responses of the Soybean MIR159 Gene Family to Low-P Stress

Various miRNAs are involved in plant responses to low-P (LP) stress. Although MIR159a is upregulated in soybean under LP stress [21], the responses of other MIR159 gene family members to LP stress were not known. To examine this issue, we analyzed the expression of MIR159 genes in soybean roots, leaves, flowers, and nodules after 7 and 40 days of LP stress via quantitative reverse-transcription PCR (qRT-PCR). Compared to those in the high-phosphorus (HP) control, after 7 days of LP stress, MIR159a and MIR159e transcript levels were not significantly altered in leaves (Figure 2A,I) but were inhibited in roots (Figure 2B,J), whereas MIR159b expression was not significantly changed in leaves but was induced in roots (Figure 2C,D), MIR159c expression was inhibited in both leaves and roots (Figure 2E,F), and MIR159d and MIR159f expression was unaltered. After 40 days of LP stress, however, MIR159a expression was increased in leaves (Figure 2A), MIR159c expression was induced in roots (Figure 2F), and MIR159e and MIR159f expression was induced in both leaves and roots (Figure 2I,L). These results suggest that miR159 family members may regulate the responses of soybean to low-P stress in various ways.
Given that MiR159 family members regulate flower development in Arabidopsis, and overexpressing MIR159 delays flowering in short-day conditions [3]. Hence, in this study, we analyzed the expression levels of soybean MIR159 gene family members in flowers after 25 days of LP stress. Compared with that under HP stress, the expression of MIR159a, MIR159b, MIR159c and MIR159f was significantly induced in flowers under LP stress, whereas MIR159d and MIR159e expression was inhibited by this treatment (Figure 3A).
Soybean is a legume crop capable of biological nitrogen fixation; some miR159 family members appear to regulate nodule growth and development in legumes [37]. Given that low-nitrogen (N) promotes nodulation, we cultured soybean plants under high-phosphorus and low-N (HPLN, control) and low-phosphorus and low-N (LPLN) conditions for 7 days and examined the samples 33 days after inoculation with rhizobia. We then examined the expression of MIR159 gene family members in nodules. Compared to HPLN, under LPLN conditions, the expression of MIR159a, MIR159c, and MIR159d in nodules was significantly inhibited, MIR159b expression was induced, and MIR159e and MIR159f expression did not significantly change (Figure 3B). These results indicate that MIR159 genes respond to low P stress in nodules.

3.3. Responses of the Soybean MIR159 Gene Family to Salt Stress

MiR159 responds to high-salt stress in S. portulacastrum [6,27]. We therefore used qRT-PCR to measure the transcript levels of MIR159 family genes in leaves and roots of soybean plants subjected to short-term (6 h) treatment with 200 mM NaCl. MIR159a (Figure 4A,B), MIR159b (Figure 4C,D), MIR159c (Figure 4E,F), and MIR159e (Figure 4I,J) were significantly upregulated in leaves and roots under this treatment; MIR159d was upregulated only in roots (Figure 4G,H); MIR159f was down-regulated in leaves but upregulated in roots (Figure 4K,L). These results suggest that the soybean MIR159 gene family plays important roles in plant responses to salt stress.

3.4. Responses of the Soybean MIR159 Gene Family to ABA

To further assess the effects of MIR159 family genes in stress response, we investigated their expression upon exposure to the hormone ABA, a downstream factor in stress response pathways. MIR159a responds to ABA treatment in Arabidopsis [16], but the responses of other soybean MIR159 family genes to ABA remained unclear. In this study, qRT-PCR demonstrated that MIR159 gene expression was induced or inhibited in response to ABA. Compared to the control (−ABA), MIR159a was upregulated in leaves (Figure 5A), but not in roots, after 6 h of ABA treatment (Figure 5B). In response to ABA treatment, MIR159b was upregulated in leaves and roots (Figure 5C,D), MIR159C was upregulated in leaves and downregulated in roots (Figure 5E,F), and MIR159d was downregulated in leaves (Figure 5G) and upregulated in roots (Figure 5H). MIR159e was down-regulated in leaves and roots (Figure 5I,J), and MIR159f was down-regulated in leaves and up-regulated in roots (Figure 5K,L). These results suggest that MIR159 family genes play important roles in regulating ABA responses in soybean.

3.5. Overexpressing MIR159e Decrease Total Phosphorus Content in Soybean Leaves under LP Conditions

To examine its functions, we overexpressed MIR159e driven by the constitutive CaMV 35S promoter in soybean genotype YC03-3. By screening for Basta resistance conferred by the Bar gene, we identified dozens of transgenic lines [34]. Following confirmation using a chi-squared test, we selected two single-copy T-DNA insertion lines, #2 and #4, for subsequent experiments. After 40 days of growth under HP (250 μM KH2PO4) or LP (5 μM KH2PO4) conditions, the growth performance of lines #2 and #4 was not significantly different from that of wild-type YC03-3 (Supplemental Figure S1). Under HP conditions, the P (Figure 6A) and total nitrogen contents (Figure 6C) in #2 and #4 plants also were not significantly different from YC03-3 in leaves or roots. However, under LP conditions, the total P contents were significantly reduced in leaves of #2 and #4 plants compared to YC03-3 (Figure 6B), although the total nitrogen contents were not significantly different from YC03-3 (Figure 6D). Finally, under both HP and LP conditions, the soluble P concentrations of #2 and #4 plants in old leaves, new leaves, and roots were not significantly different from those of wild-type plants (Figure 6F). These results indicate that overexpressing MIR159e alters total phosphorus contents in soybean.

4. Discussion

Increasing evidence indicates that miRNAs play crucial roles in plant adaptation to nutrient stress [22,23]. In the past decades, most studies on the roles of miRNAs in nutrient stress have been performed in model plants, such as Arabidopsis and rice, but little is known about their roles in legumes. Soybean is an important leguminous crop with great ability to fix nitrogen from the atmosphere, providing protein and edible oil for human consumption. MiRNAs are a class of non-coding small RNAs that function in the post-transcriptional regulation of their target genes by forming RNA-induced silencing complexes to shear the transcripts of their target genes or to inhibit the translation of these transcripts [38]. MiR159 is a highly conserved miRNA family whose members are important in regulating vegetative growth, flowering, anther development, and seed germination in various plants [3,15]. Of note, a recent study verified that rice miR159 regulates grain filling via negatively affecting miR167, and the transcription factor OsARF12 is the target of miR167, thus overexpressing OsARF12 enhance field yield [4]. In future, it is necessary to explore whether soybean miR159s modulate flowering development, seed set or grain filling.
Here, we identified six MIR159 genes in the soybean genome (Table 1), predicted cis-acting elements in their promoters, and identified candidate miR159 target genes (Supplementary Table S1). We also examined the expression patterns of MIR159 genes in response to P (Figure 2 and Figure 3), NaCl (Figure 4), and ABA (Figure 5) in soybean roots, leaves, flowers, and nodules. The highly conserved miRNA159 family is found in vascular plants and bryophytes [15]. Evolutionary analysis suggested that miR159 shares a common ancestor with miR319, another highly conserved miRNA family [39]. Soybean MIR159 genes are arranged in clusters (Table 1), similar to soybean MIR399 genes [32]. There are 3, 6, 11, 3, 3, and 1 MIR159 genes in alfalfa, maize, rice, Arabidopsis, Brachypodium distachyon, and common bean, respectively (https://www.mirbase.org/) (accessed on 4 July 2023). Arabidopsis contains three MIR159 genes [3,39], but the soybean genome contains six, suggesting that the roles of soybean miR159 family members are more complex. In addition, the larger size of the soybean MIR159 gene family might be related to the two genome duplication events that occurred during the evolution of the soybean genome [40].
Like structural genes, MIRNA genes are also regulated at the transcriptional level [2]. In this study, three or four TATA-box and W-box binding elements were found in the promoters of all six soybean MIR159 genes. Thus, it is interesting to screen WRKYs to regulate the transcription of soybean MIR159 in the near future. However, MIR159e lacks TATA-box elements, and only MIR159c and MIR159f contain P1BS binding elements (Table 2). Hence, MIR159c and MIR159f transcription might be regulated by PHR-type transcriptional regulators in soybean [41]. MIR159a was upregulated in leaves in response to 40 days of low-P stress but was downregulated in roots after 7 days of this treatment (Figure 2A), which is consistent with the results of previous studies [21]. MIR159c expression was inhibited in leaves after 7 days of low-P stress (Figure 2E) but induced in roots after 40 days (Figure 2F). MIR159e was upregulated in both leaves and roots under low-P stress (Figure 2I,J), which is consistent with our previous reports [32]. MIR159 was upregulated in leaves and roots after 40 days of low-P stress (Figure 2K,L). In addition, in flowers, MIR159a, MIR159b, MIR159c, and MIR159f expression was induced, while MIR159d and MIR159e expression was inhibited after 25 days of low-P stress (Figure 3A). Therefore, the transcription of these MIR159 genes is regulated by low-P stress and may be related to the presence of low-P-responsive elements (Table 2).
Overexpressing MIR159a in gloxinia (Sinningia speciosa) delays flowering and downregulates the expression of LEAFY (LFY), AGAMOUS, APETALA1 (AP1), and AP3 genes in flower buds [42]. Overexpressing Arabidopsis MIR159 delays flowering and anther development [3]. MIR159 and MYB33 are co-transcribed in the aleurone layer and embryo in germinating Arabidopsis seeds, and the two genes show the same spatiotemporal expression pattern [18]. The Arabidopsis mir159a mir159b double mutant shows plant dwarfing, reduced apical dominance, reduced fertility, and an irregular seed shape [18]. However, whether soybean miR159 is functionally conserved with Arabidopsis miR159 requires further study. Soybean MIR159 is expressed in floral organs and is induced by low-P stress (Figure 3A), but the role of soybean miR159 in regulating flower development remains unclear. Unlike some other stresses, low-P stress delays flowering [43]. In this study, we determined that the expression of MIR159a, MIR159b, MIR159c, and MIR159f increases under low-P conditions. This likely leads to the inhibition of flowering-related genes, such as LFY and AP, thereby delaying flowering under low-P conditions. Of course, further research is needed to confirm or disprove this theory.
Nodules are the sites of symbiosis between leguminous crops and rhizobia. Since nitrogen fixation in nodules requires ATP, maintaining the balance between P and other nutrients is crucial [44,45]. MIR159 is expressed at a relatively high level in alfalfa nodules [46]. Here, we showed that MIR159a, MIR159c, and MIR159d were downregulated, MIR159b was upregulated, and MIR159e and MIR159f expression was unchanged under low-P stress (Figure 3B). We also noted that MIR159 genes in soybean showed different response patterns to low-P stress in roots (Figure 2) and nodules (Figure 3B), possibly due to different nutritional conditions.
MIR159 expression is induced in peanut by salt stress, and miR159 regulates the expression of its target genes [26]. In addition, miR159 in different tissues of S. portulacastrum is involved in plant responses to high-salt stress [27]. In the current study, qRT-PCR results showed that MIR159 responds to high-NaCl concentrations in soybean (Figure 4). One ABRE-like and one ACGT sequence binding element were identified in each of the six MIR159 gene promoters in soybean (Table 2). However, MIR159b and MIR159c lack RD22 binding elements, which are found in the other four soybean MIR159 genes. MIR159a−MIR159e contain two, one, two, three, and four AtMYB2 binding elements, respectively, whereas MIR159f did not lack AtMYB2 binding elements (Table 2). Therefore, to understand the functions of soybean miR159s, it will be important to identify the transcriptional regulators of the response of this plant to salt stress and the binding of these regulators to the MIR159 promoters. Consistently, miR159 in potato (Solanum tuberosum subsp. andigena) is induced by salinity [8]. We noticed that the transcripts of wheat miR159 and its target gene TaGRAS178 are increased by osmotic stress [20], and tomato miR159 is repressed by drought, thus its target gene is induced [29]. Taken together, these studies indicates that plant miR159 play important roles in responses to drought, salinity, and osmotic stress.
Analysis of ABREs showed that all six MIR159 promoters in soybean contain a DPBF-binding element, whereas MIR159c and MIR159f lack ABREs and only MIR159e contains RY-binding elements (Table 2), suggesting that soybean miR159s may respond to ABA signals. In line with this notion, both MIR159a and MIR159b were upregulated by ABA treatment in Arabidopsis [16]. ABA plays an important role in regulating root development, seed maturation and germination, and drought and salt stress tolerance. Several ABA-responsive miRNAs have been identified, such as miR399f in Arabidopsis [47]. Previous studies demonstrated that the Arabidopsis hyl1 mutant is sensitive to ABA during seed germination [28], suggesting that the ABA signaling pathway might be regulated by miRNAs.
ABA signaling directly regulates stomatal opening and transpiration and is the prime candidate responsible for deciding the cellular fate under stress conditions. Moreover, ABA sensitivity is tightly associated with plant drought tolerance and seed dormancy [48,49]. Analysis of the upstream regions of the three Arabidopsis MIR159 genes revealed ABREs and binding sites for stress-related transcription factors such as AtMYC2 [30]. Jiang et al. (2022) established a signaling cascade involving miR159, MYB33, and ABA-responsive basic leucine zipper transcription factor ABI5 regulating Arabidopsis seed germination under drought stress [50]. A loss of MIR159 increases Arabidopsis drought tolerance and ABA sensitivity [51]. Sly-miR159 as an important regulator of fruit morphology in tomato, a model species of fleshy-fruit development [7]. The expression of MIR159a and MIR159b is induced by ABA in Arabidopsis at the seed germination stage, but not at the seedling stage [16]. In the current study, the results in Figure 5 suggest that miR159 family members play important roles in regulating plant responses to ABA in soybean. All six MIR159 gene promoters contain DPNF binding elements, whereas MIR159e and MIR159f lack ABREs, and only MIR159e contains RY elements (Table 2), which might help explain the different responses of individual soybean MIIR159 family members to exogenous ABA. On the other hand, strawberry MIR159a is induced by exogenous GA, but not miR159b [9], GA application enhanced the abundance of miR159, lowered the expression of its target gene GmMYB33, and delayed the development of Heterodera glycines [5]. Thus, further studies are needed to reveal the responses of soybean miR159s to other phytohormones such as GA, brassnosteroids, auxin, ethylene, and cytokinins.
In transgenic soybean plants overexpressing MIR159e and wild-type soybean plants under different P treatments for 40 days, the fresh weights of both transgenic and wild-type plants were lower under low-P vs. high-P treatment (Supplemental Figure S1A), which is consistent with previous findings [32,43]. Under HP and LP conditions, there was no significant difference in fresh weight, primary root length, root area, total P, soluble P, or total nitrogen contents between MIR159e-overexpressing transgenic and wild-type YC03-3 plants (Supplemental Figure S1). However, the primary roots of the Arabidopsis miR159ab double mutant are longer than those of the wild type, and its meristem is enlarged [19]. Under HPLN conditions, MIR159e-overexpressing soybean plants showed reduced total P contents in leaves and roots (Supplemental Figure S2A) and reduced soluble P concentrations in old leaves and nodules compared to the wild type (Supplemental Figure S2C). These results indicate that overexpressing MIR159e affects the uptake and reuse of P in soybean roots under low-N stress. MIR159e might be involved in N and P nutrition, but the physiological and molecular mechanisms are unclear. Interestingly, overexpressing TamiR159 delays rice heading, decreases male sterility, and increases the sensitivity to heat [11].
Whether overexpressing MIR159e affected the expression of a phosphate transporter gene and the abundance of its encoded protein in soybean; this should be further studied under low P or low-nitrogen conditions. Under low-N conditions, the total nitrogen content was higher in nodules than in leaves and roots (Supplemental Figure S2B), indicating that nodules play a role in nitrogen accumulation and provide nitrogen for plant growth. Overexpressing MIR159e did not affect the total nitrogen content in leaves, roots, or nodules under HPHN, LPHN, or HPLN conditions (Figure 6, Supplemental Figure S2). These results suggest that miR159e might not be involved in regulating nitrogen balance in soybean. Overexpressing MIR159 delayed flowering in ornamental gloxinia (Sinningia speciosa) [42]. However, in this study, overexpressing MIR159e in soybean did not affect flowering. It will be important to design experiments to determine whether overexpressing MIR159e alters the sensitivity of soybean to ABA and NaCl. In the future, CRISPR-Cas9 or miRNA-STTM technology is needed to knock out or knock down MIR159 expression to further analyze their functions.

5. Conclusions

In this study, we demonstrated that the soybean genome contains six MIR159 family members, which are divided into two subgroups. The promoter regions of soybean MIR159 genes contain phosphorus signal response and NaCl and ABA response elements, and the transcription of various soybean MIR159 genes is regulated by low-P stress, NaCl, and ABA, suggesting that the miR159 family plays important roles in abiotic stress in soybean. Overexpressing soybean MIR159e alters P nutrition in soybean roots under low-N conditions.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/agronomy13071798/s1, Supplemental Table S1. Predicted miR159 target genes in soybean; Supplemental Figure S1. Effects of MIR159e over expression on growth and development. WT, wild type (YC03-3); OE-2, miR159eOE-2; OE-4, miR159eOE-4; HP, plants treated with a high phosphorus concentration of 250 μM; LP, plants treated with a low-P concentration of 5 μM; treatment time was 40 days. Results are means ± SE from three independent experiments. Student’s t-test was used to compare the differences between Col-0 and transgenic MIR159eOE plants (* p < 0.05); Supplemental Figure S2. Total phosphorus, nitrogen, and soluble phosphorus contents in soybean overexpressing MIR159e under low-nitrogen conditions. WT, wild type (YC03-3); HPLN, high phosphorus and low nitrogen; treatment was performed for 40 days. Results are means ± SE from three independent experiments. Student’s t-test was used to compare the differences between Col-0 and the transgenic line MIR159eOE in the same plant part sampled at the same time point (* p < 0.05); Supplemental Table S2. List of primer pairs used in this study.

Author Contributions

J.W., B.L. and P.T. conceived and designed the experiments. B.L., P.T., F.X. and P.H. performed the experiments. J.W. and B.L. analyzed the data. J.W. and B.L. wrote the original manuscript. J.W. revised and approved the final version of the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Natural Science Foundation of China (Grant number: 31572184).

Data Availability Statement

All data are provided as figures and table which are included in this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jonesrhoades, M.W.; Bartel, D.P.; Bartel, B. MicroRNAs and their regulatory roles in plants. Annu. Rev. Plant Biol. 2006, 57, 19–53. [Google Scholar] [CrossRef] [PubMed]
  2. Rogers, K.; Chen, X. Biogenesis, turnover, and mode of action of plant microRNAs. Plant Cell 2013, 25, 2383–2399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Achard, P.; Herr, A.; Baulcombe, D.C.; Harberd, N.P. Modulation of floral development by a gibberellin-regulated microRNA. Development 2004, 131, 3357–3365. [Google Scholar] [CrossRef] [Green Version]
  4. Zhao, Y.; Zhang, X.; Cheng, Y.; Du, X.; Teotia, S.; Miao, C.; Sun, H.; Fan, G.; Tang, G.; Xue, H.; et al. The miR167-OsARF12 module regulates rice grain filling and grain size downstream of miR159. Plant Commun. 2023, 4, 100604. [Google Scholar] [CrossRef]
  5. Lei, P.; Qi, N.; Zhou, Y.; Wang, Y.; Zhu, X.; Xuan, Y.; Liu, X.; Fan, H.; Chen, L.; Duan, Y. Soybean miR159-GmMYB33 Regulatory network involved in gibberellin-modulated resistance to Heterodera glycines. Int J Mol Sci. 2021, 22, 13172. [Google Scholar] [CrossRef]
  6. Wang, M.; Wang, Q.; Zhang, B. Response of miRNAs and their targets to salt and drought stresses in cotton (Gossypium hirsutum L.). Gene 2013, 530, 26–32. [Google Scholar] [CrossRef]
  7. Zhao, P.; Wang, F.; Deng, Y.; Zhong, F.; Tian, P.; Lin, D.; Deng, J.; Zhang, Y.; Huang, T. Sly-miR159 regulates fruit morphology by modulating GA biosynthesis in tomato. Plant Biotechnol. J. 2022, 20, 833–845. [Google Scholar] [CrossRef] [PubMed]
  8. Kitazumi, A.; Kawahara, Y.; Onda, T.S.; De Koeyer, D.; de los Reyes, B.G. Implications of miR166 and miR159 induction to the basal response mechanisms of an andigena potato (Solanum tuberosum subsp. andigena). Genome 2015, 58, 13–24. [Google Scholar] [CrossRef]
  9. Csukasi, F.; Donaire, L.; Casañal, A.; Martínez-Priego, L.; Botella, M.A.; Medina-Escobar, N.; Llave, C.; Valpuesta, V. Two strawberry miR159 family members display developmental-specific expression patterns in the fruit receptacle and cooperatively regulate Fa-GAMYB. New Phytol. 2012, 195, 47–57. [Google Scholar] [CrossRef]
  10. Jiang, J.; Xu, P.; Li, Y.; Li, Y.; Zhou, X.; Jiang, M.; Zhang, J.; Zhu, J.; Wang, W.; Yang, L. Identification of miRNAs and their target genes in genic male sterility lines in Brassica napus by small RNA sequencing. BMC Plant Biol. 2021, 9, 520. [Google Scholar] [CrossRef]
  11. Wang, Y.; Sun, F.; Cao, H.; Peng, H.; Ni, Z.; Sun, Q.; Yao, Y. TamiR159 directed wheat TaGAMYB cleavage and its involvement in anther development and heat response. PLoS ONE 2012, 7, e48445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Garighan, J.; Dvorak, E.; Estevan, J.; Loridon, K.; Huettel, B.; Sarah, G.; Farrera, I.; Leclercq, J.; Grynberg, P.; Coiti Togawa, R.; et al. The identification of small RNAs differentially expressed in apple buds reveals a potential role of the Mir159-MYB regulatory module during dormancy. Plants 2021, 3, 2665. [Google Scholar] [CrossRef] [PubMed]
  13. Gao, X.; Zhang, Q.; Zhao, Y.Q.; Yang, J.; He, H.B.; Jia, G.X. The lre-miR159a-LrGAMYB pathway mediates resistance to grey mould infection in Lilium regale. Mol. Plant Pathol. 2020, 21, 749–760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zheng, Z.; Wang, N.; Jalajakumari, M.; Blackman, L.; Shen, E.; Verma, S.; Wang, M.B.; Millar, A.A. miR159 represses a constitutive pathogen defense response in tobacco. Plant Physiol. 2020, 182, 2182–2198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Millar, A.A.; Lohe, A.; Wong, G. Biology and function of miR159 in plants. Plants 2019, 8, 255. [Google Scholar] [CrossRef] [Green Version]
  16. Reyes, J.L.; Chua, N. ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination. Plant J. 2007, 49, 592–606. [Google Scholar] [CrossRef]
  17. Millar, A.A.; Gubler, F. The Arabidopsis GAMYB-Like genes, MYB33 and MYB65, are microRNA-regulated genes that redundantly facilitate anther development. Plant Cell 2005, 17, 705–721. [Google Scholar] [CrossRef] [Green Version]
  18. Allen, R.S.; Li, J.; Stahle, M.I. Genetic analysis reveals functional redundancy and the major target genes of the Arabidopsis miR159 family. Proc. Natl. Acad. Sci. USA 2007, 104, 16371–16376. [Google Scholar] [CrossRef]
  19. Xue, T.; Liu, Z.; Dai, X.; Xiang, F. Primary root growth in Arabidopsis thaliana is inhibited by the miR159 mediated repression of MYB33, MYB65 and MYB101. Plant Sci. 2017, 262, 182–189. [Google Scholar] [CrossRef]
  20. Mishra, S.; Chaudhary, R.; Sharma, P. Temporal expression analysis of microRNAs and their target GRAS genes induced by osmotic stress in two contrasting wheat genotypes. Mol Biol Rep. 2023, 13, 1–13. [Google Scholar] [CrossRef]
  21. Zeng, H.Q.; Zhu, Y.Y.; Huang, S.Q.; Yang, Z.M. Analysis of phosphorus-deficient responsive miRNAs and cis-elements from soybean (Glycine max L.). J. Plant Physiol. 2010, 167, 1289–1297. [Google Scholar] [CrossRef]
  22. Zhao, M.; Ding, H.; Zhu, J.K.; Zhang, F.; Li, W.X. Involvement of miR169 in the nitrogen-starvation responses in Arabidopsis. New Phytol. 2011, 190, 906–915. [Google Scholar] [CrossRef] [Green Version]
  23. Vidal, E.A.; Araus, V.; Lu, C.; Parry, G.; Green, P.J.; Coruzzi, G.M.; Gutiérrez, R.A. Nitrate-responsive miR393/AFB3 regulatory module controls root system architecture in Arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 2010, 107, 4477–4482. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, M.; Zhang, P.; Liu, Q.; Li, G.; Di, D.; Xia, G.; Kronzucker, H.J.; Fang, S.; Chu, J.; Shi, W. TaANR1-TaBG1 and TaWabi5-TaNRT2s/NARs link ABA metabolism and nitrate acquisition in wheat roots. Plant Physiol. 2020, 182, 1440–1453. [Google Scholar] [CrossRef] [Green Version]
  25. Hu, B.; Jiang, Z.; Wang, W.; Qiu, Y.; Zhang, Z.; Liu, Y.; Li, A.; Gao, X.; Liu, L.; Qian, Y.; et al. Nitrate-NRT1.1B-SPX4 cascade integrates nitrogen and phosphorus signalling networks in plants. Nat. Plants. 2019, 5, 401–413. [Google Scholar] [CrossRef]
  26. Zhang, T.T.; Yan, C.X.; Zhao, X.B.; Li, C.J.; Shi, C.R.; Yu, J.J.; Chan, S.H.; Li, R.G. Identification and functional characterization of salt tolerance related microRNAs in roots of peanut (Arachis hypogaea L.). J. Peanut Sci. 2016, 45, 36–42, (In Chinese with English Abstract). [Google Scholar]
  27. Tang, M.N.; Liu, Q.; Zhang, Y. Response of salt-tolerant miRNAs to high-salt habitats in mangrove companion plant hippocampal teeth. MPB 2017, 15, 1137–1142, (In Chinese with English Abstract). [Google Scholar]
  28. Lu, C.; Fedoroff, N. A Mutation in the Arabidopsis HYL1 gene encoding a dsRNA binding protein affects responses to abscisic acid, auxin, and cytokinin. Plant Cell 2000, 12, 2351–2366. [Google Scholar] [CrossRef] [Green Version]
  29. López-Galiano, M.J.; García-Robles, I.; González-Hernández, A.I.; Camañes, G.; Vicedo, B.; Real, M.D.; Rausell, C. Expression of mir159 is altered in tomato plants undergoing drought stress. Plants 2019, 8, 201. [Google Scholar] [CrossRef] [Green Version]
  30. AAbe, H.; Urao, T.; Ito, T.; Seki, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) Function as Transcriptional Activators in Abscisic Acid Signaling. Plant Cell 2003, 15, 63–78. [Google Scholar] [CrossRef] [Green Version]
  31. Huo, W.G.; Li, B.D.; Kuang, J.B. Functional characterization of the steroid reductase genes GmDET2 and GmDET2b from Glycine max. Int J Mol Sci. 2018, 19, 726. [Google Scholar] [CrossRef] [Green Version]
  32. Xu, F.; Liu, Q.; Chen, L.; Kuang, J.B.; Liao, H. Genome-wide identification of soybean microRNAs and their targets reveals their organ-specificity and responses to phosphate starvation. BMC Genom. 2013, 14, 66. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Zhou, J.; Jiao, F.; Wu, Z.; Li, Y.; Wang, X.; He, X.; Zhong, W.; Wu, P. OsPHR2 is involved in phosphate-starvation signaling and excessive phosphate accumulation in shoots of plants. Plant Physiol. 2008, 146, 1673–1686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wang, X.; Wang, Y.; Tian, J.; Lim, B.L.; Yan, X.; Liao, H. Overexpressing AtPAP15 enhances phosphorus efficiency in soybean. Plant Physiol. 2009, 151, 233–240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Kumar, S.; Stecher, G.; Peterson, D.; Tamura, K. MEGA-CC: Computing core of molecular evolutionary genetics analysis program for automated and iterative data analysis. Bioinformatics 2012, 28, 2685–2686. [Google Scholar] [CrossRef] [Green Version]
  36. Dai, X.; Zhao, P.X. psRNATarget: A plant small RNA target analysis server. Nucleic Acids Res. 2011, 39, W155–W159. [Google Scholar] [CrossRef] [Green Version]
  37. Figueredo, M.S.; Formey, D.; Rodriguez, J.; Ibanez, F.; Hernandez, G.; Fabra, A. Identification of miRNAs linked to peanut nodule functional processes. J Biosci. 2020, 45, 62. [Google Scholar] [CrossRef]
  38. Khraiwesh, B.; Zhu, J.K.; Zhu, J. Role of miRNAs and siRNAs in biotic and abiotic stress responses of plants. Biochim. Biophys. Acta (BBA)-Gene Regul. Mech. 2012, 1819, 137–148. [Google Scholar] [CrossRef] [Green Version]
  39. Li, Y.; Jin, Y.X. Evolution of MIR159/319 microRNA genes and their post-transcriptional regulatory link to siRNA pathways. BMC Evol. Biol. 2011, 11, 122. [Google Scholar] [CrossRef] [Green Version]
  40. Schmutz, J.; Cannon, S.B.; Schlueter, J.; Ma, J.X.; Mitros, T.; Nelson, W.; Hyten, D. Genome sequence of the paleopolyploid soybean (Glycine max (L.) Merr.). Nature 2010, 463, 178–183. [Google Scholar] [CrossRef] [Green Version]
  41. Bustos, R.; Castrillo, G.; Linhares, F.; Puga, M.I.; Rubio, V.; P’erez-Perez, J.; Solano, R.; Leyva, A.; Paz-Ares, J. A central regulatory system largely controls transcriptional activation and repression responses to phosphate starvation in Arabidopsis. PLoS Genet 2010, 6, e1001102. [Google Scholar] [CrossRef] [Green Version]
  42. Bian, X.Y.; Song, D.F.; Han, S.Y.; Wang, J.H.; Zhu, M.Y. Flowering time control in ornamental gloxinia (Sinningia speciosa) by manipulation of miR159 expression. Ann. Bot. 2013, 111, 791–799. [Google Scholar]
  43. Nord, E.A.; Lynch, J.P. Delayed reproduction in Arabidopsis thaliana improves fitness in soil with suboptimal phosphorus availability. Plant Cell Environ. 2008, 31, 1432–1441. [Google Scholar] [CrossRef]
  44. Olivera, M.; Tejera, N.; Iribarne, C.; Ocana, A.; Lluch, C. Growth, nitrogen fixation and ammonium assimilation in common bean (Phaseolus vulgaris): Effect of phosphorus. Physiol Plant. 2004, 121, 498–505. [Google Scholar] [CrossRef]
  45. Schulze, J.; Temple, G.; Temple, S.J.; Beschow, H.; Vance, C.P. Nitrogen fixation by white lupin under phosphorus deficiency. Ann. Bot. 2006, 98, 731–740. [Google Scholar] [CrossRef] [Green Version]
  46. Pokoo, R.; Ren, S.; Wang, Q.; Motes, C.M.; Hernandez, T.D.; Ahmadi, S.; Monteros, M.J.; Zheng, Y.; Sunkar, R. Genotype- and tissue-specific miRNA profiles and their targets in three alfalfa (Medicago sativa L.) genotypes. BMC Genom. 2018, 19, 913. [Google Scholar] [CrossRef] [PubMed]
  47. Panda, S.K.; Sunkar, R. Nutrient- and other stress-responsive microRNAs in plants: Role for thiol-based redox signaling. Plant Signal. Behav. 2015, 10, e1010916. [Google Scholar] [PubMed]
  48. Zucchero, J.C.; Caspi, M.; Dunn, K. ngl9: A third MADS box gene expressed in alfalfa root nodules. Mol. Plant-Microbe Interact. 2001, 14, 1463–1467. [Google Scholar] [CrossRef] [Green Version]
  49. Mega, R.; Abe, F.; Kim, J.S.; Tsuboi, Y.; Tanaka, K.; Kobayashi, H.; Sakata, Y.; Hanada, K.; Tsujimoto, H.; Kikuchi, J.; et al. Tuning water-use efficiency and drought tolerance in wheat using abscisic acid receptors. Nat. Plants 2019, 5, 153–159. [Google Scholar] [CrossRef]
  50. Jiang, Y.; Wu, X.; Shi, M.; Yu, J.; Guo, C. The miR159-MYB33-ABI5 module regulates seed germination in Arabidopsis. Physiol. Plantarum. 2022, 174, e13659. [Google Scholar] [CrossRef]
  51. Manish, T. Blessing in disguise: A loss of miR159 makes plant drought tolerant and ABA sensitive. Physiol. Plantarum. 2022, 174, e13763. [Google Scholar]
Agronomy 13 01798 g001 550
Figure 1. Phylogenetic analysis of the MIR159 family in soybean and other plant species. Note: Abbreviations: gma, Glycine max; ath, Arabidopsis thaliana; osa, Oryza sativa; mtr, Medicago truncatula; zma, Zea mays; bdi, Brachypodium distachyon; pvu, Phaseolus vulgaris.
Figure 1. Phylogenetic analysis of the MIR159 family in soybean and other plant species. Note: Abbreviations: gma, Glycine max; ath, Arabidopsis thaliana; osa, Oryza sativa; mtr, Medicago truncatula; zma, Zea mays; bdi, Brachypodium distachyon; pvu, Phaseolus vulgaris.
Agronomy 13 01798 g001
Agronomy 13 01798 g002 550
Figure 2. Responses of the soybean MIR159 gene family to low-P stress in leaves and roots. MIR159 transcript levels were quantified in samples from soybean plants of genotype YC03-3 that were transplanted when the first ternately compound leaf was fully unfolded and the second ternately compound leaf was not fully unfolded and treated with high and low P (HP, 250 μM KH2PO4; LP, 5 μM KH2PO4) for 7 days (7 D) or 40 days (40 D). (A) MIR159a expression in leaf, (B) MIR159a expression in root, (C) MIR159b expression in leaf, (D) MIR159b expression in root, (E) MIR159c expression in leaf, (F) MIR159c expression in root, (G) MIR159d expression in leaf, (H) MIR159d expression in root, (I) MIR159e expression in leaf, (J) MIR159e expression in root, (K) MIR159f expression in leaf, (L) MIR159f expression in root. Results are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient-deficiency conditions (ns, no significance; * p < 0.05; *** p < 0.001, **** p < 0.0001).
Figure 2. Responses of the soybean MIR159 gene family to low-P stress in leaves and roots. MIR159 transcript levels were quantified in samples from soybean plants of genotype YC03-3 that were transplanted when the first ternately compound leaf was fully unfolded and the second ternately compound leaf was not fully unfolded and treated with high and low P (HP, 250 μM KH2PO4; LP, 5 μM KH2PO4) for 7 days (7 D) or 40 days (40 D). (A) MIR159a expression in leaf, (B) MIR159a expression in root, (C) MIR159b expression in leaf, (D) MIR159b expression in root, (E) MIR159c expression in leaf, (F) MIR159c expression in root, (G) MIR159d expression in leaf, (H) MIR159d expression in root, (I) MIR159e expression in leaf, (J) MIR159e expression in root, (K) MIR159f expression in leaf, (L) MIR159f expression in root. Results are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient-deficiency conditions (ns, no significance; * p < 0.05; *** p < 0.001, **** p < 0.0001).
Agronomy 13 01798 g002
Agronomy 13 01798 g003 550
Figure 3. Responses of the soybean MIR159 gene family to low-P stress in flowers and nodules. (A) MIR159 gene expression in flowers; (B) MIR159 gene expression in nodules. Soybean genotype YC03-3 was transplanted when the first ternately compound leaf leaf was fully unfolded and the second ternately compound leaf was not fully unfolded, and treated with high and low P (HP, 250 μM KH2PO4; LP, 5 μM KH2PO4) for 25 days. Results are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient deficiency conditions (ns, no significance; * p < 0.05; *** p < 0.001).
Figure 3. Responses of the soybean MIR159 gene family to low-P stress in flowers and nodules. (A) MIR159 gene expression in flowers; (B) MIR159 gene expression in nodules. Soybean genotype YC03-3 was transplanted when the first ternately compound leaf leaf was fully unfolded and the second ternately compound leaf was not fully unfolded, and treated with high and low P (HP, 250 μM KH2PO4; LP, 5 μM KH2PO4) for 25 days. Results are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient deficiency conditions (ns, no significance; * p < 0.05; *** p < 0.001).
Agronomy 13 01798 g003
Agronomy 13 01798 g004 550
Figure 4. Responses of the soybean MIR159 gene family to NaCl stress in leaves and roots. MIR159 transcript levels were quantified in samples from soybean plants of genotype YC03-3 that were transplanted when the first ternately compound leaf compound leaf was fully unfolded and the second ternately compound leaf was not fully unfolded, and treated NaCl with for 0 or 6 h. −NaCl, 0 mM; +NaCl, 200 mM; h, hours. (A) MIR159a expression in leaf, (B) MIR159a expression in root, (C) MIR159b expression in leaf, (D) MIR159b expression in root, (E) MIR159c expression in leaf, (F) MIR159c expression in root, (G) MIR159d expression in leaf, (H) MIR159d expression in root, (I) MIR159e expression in leaf, (J) MIR159e expression in root, (K) MIR159f expression in leaf, (L) MIR159f expression in root. Results are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient deficiency conditions (ns, no significance; * p < 0.05; *** p < 0.001).
Figure 4. Responses of the soybean MIR159 gene family to NaCl stress in leaves and roots. MIR159 transcript levels were quantified in samples from soybean plants of genotype YC03-3 that were transplanted when the first ternately compound leaf compound leaf was fully unfolded and the second ternately compound leaf was not fully unfolded, and treated NaCl with for 0 or 6 h. −NaCl, 0 mM; +NaCl, 200 mM; h, hours. (A) MIR159a expression in leaf, (B) MIR159a expression in root, (C) MIR159b expression in leaf, (D) MIR159b expression in root, (E) MIR159c expression in leaf, (F) MIR159c expression in root, (G) MIR159d expression in leaf, (H) MIR159d expression in root, (I) MIR159e expression in leaf, (J) MIR159e expression in root, (K) MIR159f expression in leaf, (L) MIR159f expression in root. Results are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient deficiency conditions (ns, no significance; * p < 0.05; *** p < 0.001).
Agronomy 13 01798 g004
Agronomy 13 01798 g005 550
Figure 5. Responses of the soybean MIR159 gene family to ABA in leaves and roots. MiR159 transcript levels were quantified in samples from soybean plants of genotype YC03-3 that were transplanted when the first ternately compound leaf was fully unfolded and the ternately compound leaf was not fully unfolded and treated with ABA for 0 h or 6 h. −ABA, 0 μM; +ABA: 300 μM; h, hour. (A) MIR159a expression in leaf, (B) MIR159a expression in root, (C) MIR159b expression in leaf, (D) MIR159b expression in root, (E) MIR159c expression in leaf, (F) MIR159c expression in root, (G) MIR159d expression in leaf, (H) MIR159d expression in root, (I) MIR159e expression in leaf, (J) MIR159e expression in root, (K) MIR159f expression in leaf, (L) MIR159f expression in rootResults are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient-deficiency conditions (ns, no significance; * p < 0.05; ** p < 0.01; *** p < 0.001).
Figure 5. Responses of the soybean MIR159 gene family to ABA in leaves and roots. MiR159 transcript levels were quantified in samples from soybean plants of genotype YC03-3 that were transplanted when the first ternately compound leaf was fully unfolded and the ternately compound leaf was not fully unfolded and treated with ABA for 0 h or 6 h. −ABA, 0 μM; +ABA: 300 μM; h, hour. (A) MIR159a expression in leaf, (B) MIR159a expression in root, (C) MIR159b expression in leaf, (D) MIR159b expression in root, (E) MIR159c expression in leaf, (F) MIR159c expression in root, (G) MIR159d expression in leaf, (H) MIR159d expression in root, (I) MIR159e expression in leaf, (J) MIR159e expression in root, (K) MIR159f expression in leaf, (L) MIR159f expression in rootResults are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient-deficiency conditions (ns, no significance; * p < 0.05; ** p < 0.01; *** p < 0.001).
Agronomy 13 01798 g005
Agronomy 13 01798 g006 550
Figure 6. Overexpressing MIR159e affects total phosphorus, total nitrogen, and soluble phosphorus contents in soybean. WT, wild type (YC03-3); OL, old leaves; YL, young leaves. Soybean genotype YC03-3 was transplanted when the first ternately compound leaf was fully unfolded and the second ternately compound leaf was not fully unfolded, and treated with high and low P (HP, 250 μM KH2PO4; LP, 5 μM KH2PO4) for 40 days (40 D). (A) Total phosphorus contents in HP conditions, (B) Total phosphorus contents in LP conditions, (C) total nitrogen content in HP conditions, (D) total nitrogen content in LP conditions, (E) soluble phosphorus concentration in HP conditions, (F) soluble phosphorus concentration in LP conditions. Results are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient deficiency conditions (ns, no significance; ** p < 0.01; *** p < 0.001).
Figure 6. Overexpressing MIR159e affects total phosphorus, total nitrogen, and soluble phosphorus contents in soybean. WT, wild type (YC03-3); OL, old leaves; YL, young leaves. Soybean genotype YC03-3 was transplanted when the first ternately compound leaf was fully unfolded and the second ternately compound leaf was not fully unfolded, and treated with high and low P (HP, 250 μM KH2PO4; LP, 5 μM KH2PO4) for 40 days (40 D). (A) Total phosphorus contents in HP conditions, (B) Total phosphorus contents in LP conditions, (C) total nitrogen content in HP conditions, (D) total nitrogen content in LP conditions, (E) soluble phosphorus concentration in HP conditions, (F) soluble phosphorus concentration in LP conditions. Results are means ± SE from 3 independent experiments. Student’s t-test was used to determine the differences between control and nutrient deficiency conditions (ns, no significance; ** p < 0.01; *** p < 0.001).
Agronomy 13 01798 g006
Table
Table 1. The MIR159 family in soybean.
Table 1. The MIR159 family in soybean.
NameMature miRNA SequencePosition on Chromosome
MIR159aGAGCUCCUUGAAGUCCAAUUGGm09: 40,266,722–40,266,935 +
MIR159bGAGUUCCCUGCACUCCAAGUCGm07: 5,424,789–5,424,974 −
MIR159cAUUGGAGUGAAGGGAGCUCCGGm16: 2,830,034–2,830,218 −
MIR159dAGCUGCUUAGCUAUGGAUCCCGm09: 40,267,077–40,267,097 +
MIR159eGAGCUCCUUGAAGUCCAAUUGm07: 9,561,934–9,562,144 −
MIR159fGAGUUCCCUGCACUCCAAGUCGm16: 2,819,636–2,819,815 −
Note: Gm, Glycine max; +, sense strand; −, antisense strand.
Table
Table 2. Analysis of elements in the MIR159 promoters in soybean.
Table 2. Analysis of elements in the MIR159 promoters in soybean.
Low-P Responsive ElementsNaCl-Responsive ElementsABA-Responsive Elements
TATA -BoxTATA -Box LikeW-boxPHR1 ElementABRE-LikeACGT Sequencerd22AtMYB2MYC2ABREDPBFRY
Elements
MIR159a414 1112111
MIR159b413 11 1111
MIR159c414111 22 1
MIR159d414 1113111
MIR159e4 4 11143111
MIR159f4131111 1 1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, B.; Tao, P.; Xu, F.; He, P.; Wang, J. Function of Soybean miR159 Family Members in Plant Responses to Low Phosphorus, High Salinity, and Abscisic Acid Treatment. Agronomy 2023, 13, 1798. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071798

AMA Style

Li B, Tao P, Xu F, He P, Wang J. Function of Soybean miR159 Family Members in Plant Responses to Low Phosphorus, High Salinity, and Abscisic Acid Treatment. Agronomy. 2023; 13(7):1798. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071798

Chicago/Turabian Style

Li, Bodi, Ping Tao, Feng Xu, Pingan He, and Jinxiang Wang. 2023. "Function of Soybean miR159 Family Members in Plant Responses to Low Phosphorus, High Salinity, and Abscisic Acid Treatment" Agronomy 13, no. 7: 1798. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy13071798

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop