Skip to main content
  • Research article
  • Open access
  • Published:

MicroRNA-mediated responses to long-term magnesium-deficiency in Citrus sinensis roots revealed by Illumina sequencing

Abstract

Background

Magnesium (Mg)-deficiency occurs most frequently in strongly acidic, sandy soils. Citrus are grown mainly on acidic and strong acidic soils. Mg-deficiency causes poor fruit quality and low fruit yield in some Citrus orchards. For the first time, we investigated Mg-deficiency-responsive miRNAs in ‘Xuegan’ (Citrus sinensis) roots using Illumina sequencing in order to obtain some miRNAs presumably responsible for Citrus Mg-deficiency tolerance.

Results

We obtained 101 (69) miRNAs with increased (decreased) expression from Mg-starved roots. Our results suggested that the adaptation of Citrus roots to Mg-deficiency was related to the several aspects: (a) inhibiting root respiration and related gene expression via inducing miR158 and miR2919; (b) enhancing antioxidant system by down-regulating related miRNAs (miR780, miR6190, miR1044, miR5261 and miR1151) and the adaptation to low-phosphorus (miR6190); (c) activating transport-related genes by altering the expression of miR6190, miR6485, miR1044, miR5029 and miR3437; (d) elevating protein ubiquitination due to decreased expression levels of miR1044, miR5261, miR1151 and miR5029; (e) maintaining root growth by regulating miR5261, miR6485 and miR158 expression; and (f) triggering DNA repair (transcription regulation) by regulating miR5176 and miR6485 (miR6028, miR6190, miR6485, miR5621, miR160 and miR7708) expression. Mg-deficiency-responsive miRNAs involved in root signal transduction also had functions in Citrus Mg-deficiency tolerance.

Conclusions

We obtained several novel Mg-deficiency-responsive miRNAs (i.e., miR5261, miR158, miR6190, miR6485, miR1151 and miR1044) possibly contributing to Mg-deficiency tolerance. These results revealed some novel clues on the miRNA-mediated adaptation to nutrient deficiencies in higher plants.

Background

Magnesium (Mg)-deficiency, a common problem in many agricultural crops, occurs most frequently in strongly acidic, sandy soils, where Mg is very prone to leaching [1]. Citrus are grown mainly on acidic and strong acidic soils and Mg-deficiency is responsible for the poor fruit quality and the reduction in fruit yield in some Citrus orchards [2]. According to our investigation in 2011, over 90% and 77% of Citrus grandis orchard soils from Pinghe county, Fujian province had a pH less than 5.0 and a soil exchange Mg content less than the optimum range, respectively [3]. What’s worse, crop Mg-deficiency, which is becoming more and more popular due to soil acidification and improper farmer practices such as intensive crop production systems and highly fortified rotation, has been considered to be an urgent agricultural problem [3, 4]. Although Mg is one of the most important nutrients in higher plants and plays essential roles in numerous cellular processes such as chlorophyll biosynthesis, gas exchanges [2, 5,6,7], conformational stabilization of proteins, nucleic acids, cell walls and membranes [8], partitioning and utilization of photoassimilates [7, 9], activation of enzymes [9, 10] and reactive oxygen species (ROS) generation [9]. Despite the important roles of Mg in higher plants, Mg has been less paid attention by agronomists and botanists relative to the other nutrients and is considered to be “the forgotten element” [4, 11]. Therefore, it is very important to elucidate the molecular mechanisms on Mg-deficiency impairments and tolerance in higher plants. To our knowledge, such data are rare [10, 12, 13].

Evidence demonstrates that microRNA (miRNA)-mediated posttranscriptional regulation of gene expression plays a role in plant adaptive responses to deficiencies of phosphorus (P), potassium (K), nitrogen (N), sulfur (S), manganese (Mn), boron (B), zinc (Zn) and iron (Fe) [14,15,16,17,18,19,20]. Numerous differentially expressed miRNAs have been isolated from P-starved Arabidopsis, white lupin, Medicago truncatula, common bean, rice, barley, tomato and soybean [21,22,23,24,25]. The roles of P-deficiency-induced up-regulation of plant miR399 and miR827 in the maintenance of P homeostasis via inhibiting their targets ubiquitin-conjugating enzyme E2 24 (UBC24) and N limitation adaptation (NLA), respectively have been well characterized [14, 23, 26, 27].

Nitrogen-deficiency-induced alterations of miRNA profiles have been reported on several higher plants including maize, Arabidopsis, soybean, common bean [20, 24, 25, 28, 29]. Several N-deficiency-responsive miRNAs have been characterized in some details. For example, root modulation under N-deficiency was coordinated by miR160, miR167 and miR171 and root growth was promoted by down-regulating miR167 expression and up-regulating miR160 and miR171 expression [24, 28]. N-deficiency-induced down-regulation of miR169 has been demonstrated to be an adaptive strategy of plants to N-starvation via N-uptake and remobilization [24, 30].

Despite the vital roles of K in higher plants, little is known about K-deficiency-responsive miRNAs. In a study, Yan et al. [31] examined K-deficiency-induced alterations in expression of miR444a and its targets (i.e., MADS-57, MADS-27b, MADS-27a and MADS-23) in rice roots, and found that miR444a was slightly down-regulated and MADS-23 was greatly up-regulated.

In addition, many differentially expressed miRNAs have been identified in B-deprived C. sinensis roots and leaves [17, 18], Cu-starved [32] and Fe-deficient [33] Arabidopsis, S-deprived Brassica rapus [34], Mn-limited Phaseolus vulgaris [25] and Zn-deficient Sorghum bicolor [35].

Although the effects of nutrient deficiencies on miRNA expression in higher plants have been explored by some workers, most of these studies have been paid to herbaceous plants. Little is known about Mg-deficiency-induced alterations of miRNA expression in woody plants. Previously, we examined Mg-deficiency-responsive miRNAs in C. sinensis leaves revealed by Illumina sequencing and identified 71 down- and 75 up-regulated miRNAs, implying the potential roles of miRNAs in Citrus Mg-deficiency tolerance [36]. On this basis, we used Illumina sequencing to sequence two small RNA libraries from Mg-sufficient (control) and -deficient C. sinensis roots in order to distinguish the differences in Mg-deficiency-induced alterations of miRNA profiles between C. sinensis roots and leaves and to obtain some miRNAs presumably responsible for Citrus Mg-deficiency tolerance.

Results

Root dry weight (DW) and root and leaf Mg

Root DW and root and leaf Mg levels were lower in 0 mM Mg-treated seedlings than in 1 mM Mg-treated ones, and Mg level in leaves from 0 mM Mg-treated seedlings was much less than the sufficient range (Fig. 1) [37]. Based on these data and our previous reports [6, 12], these seedlings submitted to 0 and 1 mM Mg were regarded as Mg-deficient and -sufficient (control), respectively.

Fig. 1
figure 1

Root DW (a), root (b) and leaf (c) Mg concentrations in response to Mg-deficiency. Bars represent mean ± SD (n = 5 for root and leaf Mg and 9 for root DW). Different letters above the bars indicate a significant difference at P < 0.05

Illumina sequencing and miRNA annotation

Using high-throughput sequencing, we got 20,726,716 (22,139,574) raw reads from sRNA library constructed from control (Mg-deficient) roots. After the adaptors, low quality tags and contaminants being removed, the control and Mg-deficient root sRNA libraries generated 20,325,777 (5,561,214) and 21,783,568 (6,124,980) clear reads (unique reads), respectively (Table 1). As shown in Fig. 2, the majority of the clear reads fell within the range of 18–25 nt. The most abundant clear reads were 24 nt length, followed by 21, 22, 23 and 20 nt length. This agrees with the previous data obtained on leaves, roots [17, 18] and fruits [38] of C. sinensis, and fruits and flowers of Citrus trifoliata [39]. Therefore, these data obtained via high-throughput sequencing of sRNA libraries are reliable. Mg-deficiency increased and decreased the abundances of 24 and 21 nt reads, respectively.

Table 1 Summary of sRNA sequencing data from Mg-sufficient and -deficient Citrus sinensis roots
Fig. 2
figure 2

Small RNA length distribution from Mg-deficient and -sufficient C. sinensis roots

Here, 13,624,836 clean reads (3,077,845 unique reads) from Mg-sufficient roots and 14,510,776 clear reads (3,378,231 unique reads) from Mg-deficient roots were mapped to C. sinensis genome (JGIversion1.1, http://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Csinensis) using SOAP [40]. Thereafter, we used the unannotated 5,050,734 and 5,550,820 unique reads from Mg-sufficient and -deficient roots, respectively to predict novel miRNAs (Table 1).

Identification and prediction of root miRNAs

As shown in Additional file 1, we identified 733 known miRNAs in C. sinensis roots. To avert false results due to the use of low abundant miRNAs, these known miRNAs with a transcript per million (TPM) value <10 in both Mg-sufficient and -deficient roots were not utilized for further analysis [17, 41]. The remained 300 miRNAs with a TPM value ≥10 in Mg-sufficient and/or -deficient roots were utilized for Mg-deficiency-responsive miRNA analysis (Additional file 2). As shown in Additional file 3, we obtained 71 up- and 54 down-regulated known miRNAs from Mg-deficient roots.

As shown in Additional files 4, 5, and 6, we identified 181 novel miRNAs in both Mg-sufficient and -deficient roots, and 30 up- and 15 down-regulated novel miRNAs in Mg-deficient roots with a TPM value ≥ ten in Mg-deficient and/or -sufficient roots.

Validation of sequencing data by stem-loop qRT-PCR

The expression levels of 27 Mg-deficiency-responsive miRNAs were assayed by stem-loop qRT-PCR. Except for miR1222, the expression patterns of all miRNAs obtained by stem-loop qRT-PCR and Illumiona sequencing were similar (Fig. 3 and Table 2). Thus, the results produced by Illumiona sequencing were reliable.

Fig. 3
figure 3

Relative expression levels of selected Mg-deficiency-responsive known miRNAs in Mg-deficient and control roots revealed by qRT-PCR. Bars represent mean ± SD (n = 3). For the same miRNA, different letters above the bars indicate a significant difference at P < 0.05. All the values were expressed relative to the control roots

Table 2 qRT-PCR analysis of predicted target genes for selected Mg-deficiency-responsive known miRNAs in C. sinensis roots

Prediction and GO analysis of targets for Mg-deficiency-responsive miRNAs

Here, we predicted 239 and 130 target genes from the 46 known and 15 novel Mg-deficiency-responsive-miRNAs, respectively (Additional files 7 and 8). As shown in Fig. 4a, the targets for known (novel) Mg-deficiency-responsive miRNAs were associated with 12 (nine) biological processes. The most abundant three GO terms were response to stress, transport and protein process for known miRNA targets and response to stress, regulation of transcription and transport for novel miRNA targets, respectively. On the basis of the molecular function, the highest percentages of three groups for known and novel miRNA targets were nucleic acid binding, other activity and kinase activity, and other activity, metal ion binding and transporter activity, respectively (Fig. 4b). As shown in Fig. 4c, the targets for known (novel) Mg-deficiency-responsive miRNAs were related to 12 (eight) cellular components. The most abundant component for known and novel miRNAs was nucleus.

Fig. 4
figure 4

GO categories of the predicted target genes for 46 (15) Mg-deficiency-responsive known (novel) miRNAs in Citrus sinensis roots. MiRNAs target genes were grouped based on biological process (a), molecular function (b) and cellular component (c)

Validation of target genes by qRT-PCR

As shown in Table 2, 105 targets for 11 up- and 16 down-regulated miRNAs were validated by qRT-PCR. As expected, we found that 65 (61.9%) targets and their corresponding miRNAs displayed opposite trends in expression profiles in Mg-deprived roots, suggesting that miRNAs might play a role in gene regulation by cleaving mRNAs. However, 34 (32.4%) targets displayed the same expression profiles with their corresponding miRNAs in Mg-deficient roots or were not significantly affected by Mg-deficiency. The remaining 6 (0.06%) targets were not detected in roots. It is worth mentioning that 4 target genes (i.e., range1.1g005482m, orange1.1g004896m, orange1.1g005075m and orange1.1g008078m) belonging to auxin responsive factor (ARF) family have been validated by us in C. sinensis [42], suggesting that the target prediction was accurate.

Discussion

Little is known about the possible roles of miRNAs in plant Mg homeostasis [36, 43]. Here, we first investigated the Mg-deficiency-induced alterations of miRNA profiles in Citrus roots and obtained 101 up- and 69 down-regulated miRNAs (Additional files 3 and 6), demonstrating that miRNAs might be involved in Mg-deficiency responses. We obtained similar amount of miRNAs (71 miRNAs) with decreased expression, but less amount of miRNAs (75 miRNAs) with increased expression from Mg-deficient C. sinensis leaves compared with Mg-deficient C. sinensis roots [36]. Moreover, most of these miRNAs were isolated only from Mg-deprived roots or leaves, only 30 Mg-deficiency-responsive miRNAs were shared by the two. Among the 30 overlapping miRNAs, only 15 miRNAs displayed similar expression trends in Mg-deprived roots and leaves (Table 3). Thus, great differences existed in Mg-deficiency-induced alterations of miRNA profiles between roots and leaves. This agrees with our report that the physiological and biochemical responses to long-term Mg-deficiency greatly differed between C. sinensis roots and leaves [7].

Table 3 List of Mg-deficiency-responsive known miRNAs shared by both C. sinensis roots and leaves

We observed that miR158 was induced in Mg-deprived roots (Table 2). Similar results have been obtained on Mg-deficient C. sinensis leaves [36], P-deficient tomato roots and leaves [22], and B-starved C. sinensis roots and leaves [17, 18]. As expected, its target gene SPFH (stomatins, prohibitins, flotillins and HflK/C)/Band 7/PHB domain-containing membrane-associated protein family (AT5G62740) was repressed in Mg-deprived roots. Wang et al. [44] found that Arabidopsis phb3–3 mutants were less sensitive to salt-stress-induced inhibition of primary root growth. Thus, the down-regulation of AT5G62740 might contribute to Citrus Mg-deficiency tolerance via alleviating Mg-deficiency-induced inhibition of root growth (Fig. 1a). Gehl et al. [45] observed that the basal tissue respiration rate in stomatin-like protein 1 (slp1) knockout Arabidopsis roots was reduced by 30% compared with wild-type. In addition, miR2919 expression was induced and its target: phosphoenolpyruvate carboxylase 3 (PEPC3) was inhibited in Mg-deprived C. sinensis roots (Table 2). Therefore, root respiration might be decreased in Mg-starved C. sinensis roots. This agrees with our reports that the abundances of pyruvate decarboxylase (gi|255,579,310) and phosphoglycerate kinase (gi|332,195,235) in glycolysis and the activities of key enzymes in glycolysis and tricarboxylic acid (TCA) cycle were reduced in Mg-deprived C. sinensis roots accompanied by decreased accumulation of carbohydrates and lower respiration [7, 12].

Both root miR6278 and its targets: NB-ARC domain-containing disease resistance protein involved in disease resistance and DnaJ/Hsp40 cysteine-rich domain superfamily protein (AT3G14470.1) were induced by Mg-deficiency. In addition, NB-ARC domain-containing disease resistance protein (AT4G27220.1) targeted by miR780 was induced in Mg-depprived roots (Table 2). Similarly, the abundances of Grp94 (HSP; gi|23,477,636) and disease resistance protein (gi|227,438,123) was increased in Mg-starved C. sinensis roots [12]. Thus, disease resistance might be elevated in these roots with increased levels of Ca and K [46], which contribute to plant disease resistance [47, 48].

MiR414 mainly targets transcriptional regulators including MYB, bZIP family transcription factors, WRKY and scarecrow and might have key roles in plant growth and development [49]. As expected, miR414 was up-regulated and its target gene: poly(A) polymerase 1 was inhibited in Mg-deprived roots (Table 2). This was also supported by our reports that the abundances of transcription factor homolog BTF3-like protein (gi|33,945,882), spliceosome RNA helicase BAT1 (gi|226,528,292) and RNA polymerase β chain (gi|90,403,817) were lowered in Mg-deficient C. sinensis roots [12].

We found that miR1847 was inhibited in Mg-deprived roots (Table 2). This agrees with the results obtained on B-deficient roots [17] and disagrees with the data obtained on B-deficient leaves [18]. As expected, its target genes: ribosomal protein S3 family proteins were up-regulated in these roots. In addition, VALRS targeted by miR6485 were induced or was little affected in Mg-starved roots (Table 2). Thus, protein biosynthesis might not be lowered in Mg-starved roots, as shown by unchanged concentration of total soluble proteins in Mg-deprived C. sinensis roots [7, 12]. Also, the reduced dilution due to the decrease in root DW (Fig. 1a) might account for the unchanged protein level.

MiR5176 was induced in Mg-deprived roots (Table 2), as found on B-starved C. sinensis roots [17]. DNA mismatch repair (MMR) system is required for the correction of DNA biosynthetic errors [50]. MUTL-homologue 1 (MLH1) participates in DNA MMR, correcting DNA damage and insertion-deletion loops arising from DNA replication [51]. MLH1 targeted by miR5176 were induced rather than inhibited in Mg-deprived roots (Table 2). Thus, MMR system might be up-regulated in these roots, thus enhancing Citrus Mg-deficiency tolerance via correcting DNA biosynthetic errors. Similarly, DNA repair and meiosis protein (Mre11) targeted by miR5261 was induced in Mg-deprived roots (Table 2).

Plant leucine-rich repeat receptor-like kinase proteins play crucial roles in abiotic stresses [52]. MiR5198 and its target gene: leucine-rich receptor-like protein kinase (LR-RLK) family protein were repressed and induced in Mg-starved roots, respectively (Table 2). Similarly, leucine-rich repeat receptor-like protein kinase (ACN59310.1) was up-regulated in Mg-starved C. reticulata roots [13]. Thus, miR5198 might be involved in Citrus Mg-deficiency responses.

Root miR780 was repressed by Mg-deficiency (Table 2), as found on B-starved C. sinensis roots [17]. As expected, its targets: NB-ARC domain-containing disease resistance protein and phytoene synthase (PSY) were up-regulated in Mg-starved roots (Table 2). Various transgenic plants over-expressing bacterial or plant gene encoding PSY, a major rate-limiting carotenoid (Car) enzyme, displayed increased Car level [53, 54]. Therefore, Car biosynthesis might be enhanced in Mg-starved C. sinensis roots, thus increasing their antioxidant ability.

We observed that miR7121 and its target gene: phosphatidic acid phosphatase (PAP2) family protein was up- and down-regulated in Mg-starved roots, respectively (Table 2). Nakano et al. [55] demonstrated that the inhibition of PAP2 expression or function conferred resistance to Ralstonia solanacearum via rapidly triggering plant defenses in Nicotiana benthamiana. Thus, the down-regulation of PAP2 family protein might contribute to plant disease-resistance.

As shown in Table 2, miR6190 was down-regulated and its target genes [i.e., purple acid phosphatase 26 (PPAP26), Sec23/Sec24 protein transport family protein and ortholog of human splicing factor SC35, also known as serine/arginine-rich (SR) splicing factor SC35] were up-regulated in Mg-deficient roots. SR proteins are required for regulating alternative splicing. In higher plants, great alterations in alternative splicing due to various abiotic stresses demonstrate the roles of SR proteins in the adaptation to environmental stress [56]. Induction of acid phosphatases (APases) by P-starvation is a well-documented mechanism of plant P-deficiency tolerance. Hurley et al. [57] demonstrated that AtPPAP26 was the major contributor to P-deficiency-inducible APase activity. In addition, AtPPAP26 also showed alkaline peroxidase (POD) activity. Mg-deficiency-induced up-regulation of root PPAP26 (Table 2) agrees with the report that AtPPAP26 was induced in P-deficient Arabidopsis roots, shoots and suspension cells [58] because C. sinensis leaf, stem and root P levels were reduced by Mg-deficiency [46]. Coat protein complex II (COPII) vesicles play an essential role for the export of secretory cargo from the endoplasmic reticulum (ER) to the Golgi complex in all eukaryotes [59]. Mg-deficiency-induced up-regulation of root gene encoding Sec23/Sec24 protein transport family protein (Table 2), a subset of the COPII components, agrees with our report that the abundance of Sec23/Sec24 protein transport family protein was elevated in B-deficient roots [60].

In Arabidopsis, miRNA160 negatively regulates the repressor auxin response factor (ARF) family: ARF17 [61], ARF16 [62] and ARF10 [63]. The repression of these genes by miR160 is required for seed germination and the normal development of roots, stems and leaves. Li et al. [64] demonstrated that soybean miR160a negatively regulated the progress of leaf senescence via repressing its targets: ARF10, ARF16 and ARF17. We found that miR160 was induced in Mg-deficient roots (Table 2), as obtained on P-starved Lupinus albus roots [65] and N-deficient maize roots [66]. Therefore, the induction of root miR160 by Mg-deprivation might be an adaptive response. Unexpectedly, its targets: ARF10, ARF16 and ARF17 were up-regulated in Mg-deprived C. sinensis roots (Table 2). Endogenous target mimics (eTMs) can impede the interaction between miRNAs and their authentic targets via binding to miRNAs [67, 68]. Lin et al. [69] demonstrated that eTMs repressed miR160-mediated cleavage of ARF10, ARF16 and ARF17 during longan somatic embryogenesis. No negative correlations were observed among the levels of miR160 and ARF10, ARF16 and ARF17 transcripts in longan vegetative and generative tissues. Thus, the correlations between miR160 and its targets in C. sinensis roots can be explained in this way.

Root miR6485 was repressed and its several target genes were up-regulated by Mg-deprivation (Table 2). Li et al. [70] reported that ARF7 (AT5G20730) is necessary for both auxin signaling and ethylene responses in Arabidopsis roots. Okushima et al. [71] observed that lateral root formation was badly damaged in Arabidopsis arf7 arf19 double knockout mutant, concluding that ARFs directly activated LATERAL ORGAN BOUNDARIES DOMAIN/ASYMMETRIC LEAVES2-LIKE (LBD/ASL) genes, thus regulating lateral root formation. Thus, the induction of root transcriptional factor B3 family protein/auxin-responsive factor AUX/IAA-related by Mg-deficiency might play a part in Mg-deficiency tolerance via maintaining lateral root formation. Similarly, SnRNA activating complex family protein (SDR2), which is associated with auxin-activated signaling pathway, was induced in Mg-deprived roots (Table 2). Ohtani et al. [72] reported that srd2 mutation repressed the expression of PIN-FORMED proteins, which might account for the failure to generate an auxin gradient, thus leading to different abnormalities in root morphogenesis in Arabidopsis mutant. Nuclear pore complex protein NUP88 is necessary for systemic acquired resistance and R protein-mediated defense [73]. The induction of root NUP88 by Mg-deficiency (Table 2) agrees with the above inference that that disease-resistance was improved in Mg-deficient roots.

Root miR1044 was repressed and its target genes [i.e., NF-X-like 1 (NFXL1), POD superfamily protein and ATPase, AAA-type, CDC48 protein] were induced by Mg-deprivation except for integrase-type DNA-binding superfamily protein (Table 2). Lisso et al. [74] observed that AtNFXL1 was induced in roots under salt and osmotic stress, and that both AtNFXL1-antisense plants and atnfxl1–1 knock-out mutants had lower growth and survival rates than wild-type plants when exposed to salt or osmotic stress. CDC48, a member of AAA-ATPase family proteins that provides energy for plant development via regulating ATPase, is required for plant cell division, expansion and differentiation [75]. Wang et al. [76] suggested that the induction of PpCDC48II by low temperature played a key role in cold-induced freezing tolerance of Physcomitrella patens cells.

MiR5261 and its target genes were repressed and induced in Mg-starved roots, respectively (Table 2). The induction of root RNA binding (RRM/RBD/RNP motifs) family protein by Mg-deficiency agrees with our report that the abundance of RNA binding (RRM/RBD/RNP motifs) family protein was elevated in B-deficient roots [60]. A typical mitogen-activated protein kinase (MAPK) cascade is composed of three sequentially activated protein kinases, namely MAPK, MAPK kinase (MAPKK) and MAPKK kinase (MAPKKK). Stress-tolerance of some plants such as Arabidopsis, tobacco and cereals has been enhanced by genetically altering the abundances and/or the activities of some MAPK components [77, 78]. The induction of root protein phosphatase 2A (PP2A) regulatory B subunit family protein by Mg-deficiency (Table 2) agrees with the report that wheat root PP2AbB"-α was up-regulated when exposed to various abiotic stresses. Transgenic wheat lines over-expressing TaPP2AbB"-α displayed better lateral root development, especially under NaCl and mannitol stresses [79].

Ubiquitination-proteasomal pathway has been shown to function in plant senescence and in stress response by facilitating the degradation of bulk proteins for N recycling [80]. Transgenic tobacco lines over-expressing a maize gene encoding E3 ubiquitin ligase (UPL) displayed increased drought tolerance accompanied by higher activities of superoxide dismutase (SOD) and catalase, more accumulation of proline and less accumulation of malondialdehyde (MDA) and ROS when exposed to drought stress [81]. Over-expression of TaFBA1 encoding F-box protein conferred drought and oxidative stress tolerance in tobacco plants via up-regulating the activities of SOD, catalase, ascorbate peroxidase (APX) and POD, and lowering the levels of ROS and MDA [82, 83]. Thus, up-regulation of UPL1 and RING/U-box superfamily protein targeted by miR5261, NFXL1 targeted by miR1044, F-box family protein targeted by miR1151 and proteasome subunit beta type-6 (PBA1) targeted by miR5029 in Mg-starved roots (Table 2) might confer stress-tolerance, thus contributing to Mg-deficiency tolerance in Citrus plants. Similarly, the expression levels of UPL5 (XP 003594229.1) and F-box family protein (XP 003612153.1) in C. reticulata roots [13] and the abundances of putative proteasome subunit alpha type (gi|255,584,432) in C. sinensis roots [12] were elevated by Mg-deficiency.

Conclusions

For the first time, we used Illumina sequencing to identify 71 known and 30 novel miRNAs with increased expressed, and 54 known and 15 novel miRNAs with decreased expression in Mg-deficient C. sinensis roots, demonstrating that miRNAs might be involved in Citrus Mg-deficiency tolerance. Through integrating our findings with the previous data, we put forward a potential scheme for the responses of miRNAs to Mg-deficiency in Citrus roots (Fig. 5). Here, we obtained several novel Mg-deficiency-responsive miRNAs (i.e., miR5261, miR158, miR6190, miR6485, miR1151 and miR1044) possibly responsible for Citrus Mg-deficiency tolerance. Our findings results not only increased our knowledge on the functions of plant miRNAs under nutrient deficiencies, but also established foundation to improve Mg-deficiency tolerance via manipulating the actions of miRNAs.

Fig. 5
figure 5

A potential scheme for responses of C. sinensis roots miRNAs to Mg-deficiency. LRRTPK: Leucine-rich repeat transmembrane protein kinase; LYK5: LysM-containing receptor-like kinase 5; PEX14: Peroxin 14; RBFP: RNA binding (RRM/RBD/RNP motifs) family protein; RLP: Receptor like protein; RPS3: Ribosomal protein S3 family protein; WRKY72: WRKY DNA-binding protein 72

Methods

Citrus sinensis Seedling culture and long-term Mg-deficient treatments

Seedling culture and long-term Mg-deficient treatments were carried out as described previously [12]. In short, 15-week-old ‘Xuegan’ [Citrus sinensis (L.) Osbeck] seedlings, which were grown in 6 L pots (two seedlings per pot) filled with clean river sand in a greenhouse under natural photoperiod at Fujian Agriculture and Forestry University, Fuzhou, were supplied every other day until dripping with nutrient solution at a Mg concentration of 0 mM (Mg-deficiency) or 1 mM (Mg-sufficiency, control) from MgSO4. S at the nutrient solution was kept at a constant level by adding equivalent moles of Na2SO4 in replace of MgSO4. After 16 weeks, ~ 5-mm-long root apices from new white fibrous roots were harvested and immediately frozen in liquid N2, then stored at −80 °C until extraction. The seedlings not being sampled were used for the measurements of root DW, leaf and root Mg.

Root DW and root and leaf Mg

For each treatment, roots from nine seedlings (one seedling per pot) were taken. Root DW was measured after being dried at 70 °C to a constant weight (~ 48 h).

Fibrous roots and ~7-week-old leaves (midribs and petioles removed) were harvested and then dried at 70 °C to a constant weight. Dried roots and leaves were ground to pass a 40 mesh sieve, finally digested with 1 N HCl [84]. Mg concentration in the solution was measured by atomic absorption spectroscopy.

Root sRNAs library construction, high-throughput sequencing, annotation and miRNA identification

Equal amounts of frozen root apices from five seedlings (one seedling per pot) were pooled as a biological replicate. There was one biological replicate for each treatment. Approximately 0.1 g mixed frozen Mg-deficient or control root apices were used to extract total RNA with TRIzol reagent (Invitrogen, Carlsbad, CA). Construction of sRNA libraries was performed as described by Lu et al. [17]. Illumina sequencing was carried out with a Solexa sequencer at the Beijing Genomics Institute (BGI), Shenzhen, China.

Both sRNA annotation and miRNA identification were made as described previously [17, 18]. After raw data being analyzed with a software developed by BGI, clean reads were then utilized to assay length distribution. Finally, the clear reads were mapped to C. sinensis genome (JGIversion 1.1, http://phytozome.jgi.doe.gov/pz/portal.html#!info?alias=Org_Csinensis) using SOAP, only perfectly mapped sequences were retained and analyzed further. rRNAs, tRNAs, snRNAs and snoRNAs were removed from the sRNAs sequences through BLASTn search using NCBI Genebank database (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi/) and Rfam (12.0) database (http://www.sanger.ac.uk/resources/databases/rfam.html) (e = 0.01). The remaining sequences were aligned with known plant miRNAs from miRBase 21 (http://www.mirbase.org/). Only the perfectly matched sequences were considered to be conserved miRNAs. Reads not being annotated were used for the prediction of novel miRNAs using Mireap (http://sourceforge.net/projects/mireap/), a software developed by BGI. Also, both DNAMAN 8 (http://www.lynnon.com/pc/framepc.html) and MTide (http://bis.zju.edu.cn/MTide) [85] were used for the prediction of novel miRNAs. Only these miRNA candidates being simultaneously predicted by the three softwares were regarded to be novel miRNAs.

Differentially expressed miRNAs and target prediction

Both the fold change between Mg-deficient and -sufficient roots and the P-value were calculated from the normalized expression of TPM [86]. A miRNA was regarded to be differentially expressed when it had both a P-value <0.01 and a log2-fold change >1.5 [17]. Target prediction of miRNAs was carried out by RNAhybrid according to the rules proposed by Schwab et al. [87] and Allen et al. [88].

Categories of the potential targets predicted from Mg-deficiency-responsive miRNAs

All target genes predicted from Mg-deficiency-responsive miRNAs were mapped to GO terms in the database (http://www.geneontology.org/), and gene numbers for each term was calculated. All these targets were grouped into three categories: biological process, molecular function, cellular component [17].

Validation of Mg-deficiency-responsive miRNAs by stem-loop qRT-PCR and of target genes by qRT-PCR

Stem-loop qRT-PCR analysis of miRNAs was carried out as described previously [18]. Stem-loop primers for reverse transcription and primers for qRT-PCR were summarized in Additional file 9. qRT-PCR analysis of target genes was carried out with an ABI 7500 Real Time System as described by Lu et al. [17]. The sequences of the F and R primers used were given in Additional file 10. Equal amounts of frozen root apices from five seedlings (one seedling per pot) were pooled as a biological replicate. For each treatment, there were three biological replicates and two technical replicates. Relative miRNA expression was calculated using ddCt algorithm. Actin (AEK97331.1) was used as an internal standard and the roots from Mg-sufficient seedlings were used as reference sample, which was set to 1.

Experimental design and data analysis

For each treatment, there were 20 pot seedlings in a completely randomized design. Experiments were carried out with 3 replicates except for high-throughput sequencing (n = 1), root and leaf Mg (n = 5), and root DW (n = 9). Unpaired t-test was performed for the significant test between two means (Mg-sufficiency and -deficiency).

Abbreviations

AHL9:

AT-hook motif nuclear-localized protein 9

APase:

Acid phosphatase

APX:

Ascorbate peroxidase

ARF:

Auxin responsive factor

Car:

Carotenoid

CDC48:

A member of AAA-ATPase family proteins

COPII:

Coat protein complex II

ER:

endoplasmic reticulum

eTMs:

Endogenous target mimics

LR-RLK:

Leucine-rich receptor-like protein kinase

MAPK:

Mitogen-activated protein kinase

MAPKK:

MAPK kinase

MAPKKK:

MAPKK kinase

MDA:

Malondialdehyde

Mg:

Magnesium

MLH1:

MUTL-homologue 1

MMR:

DNA mismatch repair

NUP88:

Nuclear pore complex protein

NFXL1:

NF-X-like 1

NLA:

Nitrogen limitation adaptation

PAP:

Phosphatidic acid phosphatase

PBA1:

Proteasome subunit beta type-6

PEPC:

Phosphoenolpyruvate carboxylase

PHB:

Prohibitin

PLA2A:

Phospholipase A 2A

POD:

Peroxidase

PP2A:

Protein phosphatase 2A

PPAP:

Purple acid phosphatase

PSY:

Phytoene synthase

ROS:

Reactive oxygen species

SC35:

Ortholog of human splicing factor

SDR2:

SnRNA activating complex family protein

SLP1:

Stomatin-like protein 1

SOD:

Superoxide dismutase

SPFH:

Stomatins, prohibitins, flotillins and HflK/C

TPM:

Transcript per million

UCB24:

Ubiquitin-conjugating enzyme E2 24

UPL1:

Ubiquitin-protein ligase 1 (E3 ubiquitin-protein ligase UPL1)

VALRS:

Valyl-tRNA synthetase / valine-tRNA ligase.

References

  1. Hopkins WG, Hüner NPA. Introduction to plant physiology. 4th ed. New York: John Wiley & Sons, Inc.; 2008.

    Google Scholar 

  2. Ling LL, Peng LZ, Cao L, Jiang CL, Chun CP, Zhang GY, Wang ZX. Effect of magnesium deficiency on photosynthesis characteristic of Beibei 447 Jinchen orange. J Fruit Sci. 2009;26:275–80.

    CAS  Google Scholar 

  3. Li Y, Han MQ, Lin F, Ten Y, Lin J, Zhu DH, Guo P, Weng YB, Chen LS. Soil chemical properties, ‘Guanximiyou’ pummelo leaf mineral nutrient status and fruit quality in the southern region of Fujian province. China J Soil Sci Plant Nutr. 2015;15:615–28.

    CAS  Google Scholar 

  4. Guo W, Hussain H, Liang Z, Yang D. Magnesium deficiency in plants: An urgent problem. Crop J. 2016;4:83–91.

    Article  Google Scholar 

  5. Tang N, Li Y, Chen LS. Magnesium deficiency-induced impairment of photosynthesis in leaves of fruiting Citrus reticulata trees accompanied by up-regulation of antioxidant metabolism to avoid photooxidative damage. J Plant Nutr Soil Sci. 2012;175:784–93.

    Article  CAS  Google Scholar 

  6. Yang GH, Yang LT, Jiang HX, Wang P, Chen LS. Physiological impacts of magnesium-deficiency in Citrus seedlings: photosynthesis, antioxidant system and carbohydrates. Trees Struct Funct. 2012;26:1237–50.

    Article  CAS  Google Scholar 

  7. Li CP, Qi YP, Zhang J, Yang LT, Wang DH, Ye X, Lai NW, Tan LL, Lin D, Chen LS. Magnesium-deficiency-induced alterations of gas exchange, major metabolites and key enzymes differ among roots, lower and upper leaves of Citrus sinensis seedlings. Tree Physiol. 2017;DOI: 10.1093/treephys/tpx067

  8. Cowan JA. Structural and catalytic chemistry of magnesium-dependent enzymes. Biometals. 2002;15:225–35.

    Article  CAS  PubMed  Google Scholar 

  9. Cakmak I, Kirkby EA. Role of magnesium in carbon partitioning and alleviating photooxidative damage. Physiol Plant. 2008;133:692–704.

    Article  CAS  PubMed  Google Scholar 

  10. Verbruggen N, Hermans C. Physiological and molecular responses to magnesium nutritional imbalance in plants. Plant Soil. 2013;368:87–99.

    Article  CAS  Google Scholar 

  11. Cakmak I, Yazici AM. Magnesium: a forgotten element in crop production. Better Crops. 2010;94:23–5.

    Google Scholar 

  12. Peng HY, Qi YP, Lee J, Yang LT, Guo P, Jiang HX, Chen LS. Proteomic analysis of Citrus sinensis roots and leaves in response to long-term magnesium-deficiency. BMC Genomics. 2015;16:253.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Jin XL, Ma CL, Yang LT, Chen LS. Alterations of physiology and gene expression due to long-term magnesium-deficiency differ between leaves and roots of Citrus reticulata. J Plant Physiol. 2016;198:103–15.

    Article  CAS  PubMed  Google Scholar 

  14. Shukla LI, Chinnusamy V, Sunkar R. The role of microRNAs and other endogenous small RNAs in plant stress responses. Biochim Biophys Acta. 2008;1779:743–8.

    Article  CAS  PubMed  Google Scholar 

  15. Khraiwesh B, Zhu JK, Zhu J. Role of miRNAs and siRNAs in biotic and abiotic stress responses of plants. Biochim Biophys Acta. 1819;2012:137–48.

    Google Scholar 

  16. Kulcheski FR, Côrrea R, Gomes IA, de Lima JC, Margis R. NPK macronutrients and microRNA homeostasis. Front Plant Sci. 2015;6:451.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Lu YB, Yang LT, Qi YP, Li Y, Li Z, Chen YB, Huang ZR, Chen LS. Identification of boron-deficiency-responsive microRNAs in Citrus sinensis roots by Illumina sequencing. BMC Plant Biol. 2014;14:123.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Lu YB, Qi YP, Yang LT, Guo P, Li Y, Chen LS. Boron-deficiency-responsive microRNAs and their targets in Citrus sinensis leaves. BMC Plant Biol. 2015;15:271.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Paul S, Datta SK, Datta K. MiRNA regulation of nutrient homeostasis in plants. Front Plant Sci. 2015;6:232.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Zeng H, Wang G, Hu X, Wang H, Du L, Zhu Y. Role of microRNAs in plant responses to nutrient stress. Plant Soil. 2014;374:1005–21.

    Article  CAS  Google Scholar 

  21. Branscheid A, Sieh D, Pant BD, May P, Devers EA, Elkrog A, Schauser L, Scheible WR, Krajinski F. Expression pattern suggests a role of miR399 in the regulation of the cellular response to local pi increase during arbuscular mycorrhizal symbiosis. Mol Plant-Microbe Interact. 2010;23:915–26.

    Article  CAS  PubMed  Google Scholar 

  22. Gu M, Xu K, Chen A, Zhu Y, Tang G, Xu G. Expression analysis suggests potential roles of microRNAs for phosphate and arbuscular mycorrhizal signaling in Solanum lycopersicum. Physiol Plant. 2010;138:226–37.

    Article  CAS  PubMed  Google Scholar 

  23. Hsieh LC, Lin SI, Shih ACC, Chen JW, Lin WY, Tseng CY, Li WH, Chiou TJ. Uncovering small RNA-mediated responses to phosphate deficiency in Arabidopsis by deep sequencing. Plant Physiol. 2009;151:2120–32.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Nguyen GN, Rothstein SJ, Spangenberg G, Kant S. Role of microRNAs involved in plant response to nitrogen and phosphorous limiting conditions. Front Plant Sci. 2015;6:629.

    PubMed  PubMed Central  Google Scholar 

  25. Valdés-López O, Yang SS, Aparicio-Fabre R, Graham PH, Reyes JL, Vance CP, Hernández G. MicroRNA expression profile in common bean (Phaseolus vulgaris) under nutrient deficiency stresses and manganese toxicity. New Phytol. 2010;187:805–18.

    Article  PubMed  Google Scholar 

  26. Chiou TJ. The role of microRNAs in sensing nutrient stress. Plant Cell Environ. 2007;30:323–32.

    Article  CAS  PubMed  Google Scholar 

  27. Liu TY, Lin WY, Huang TK, Chiou TJ. MicroRNA-mediated surveillance of phosphate transporters on the move. Trends Plant Sci. 2014;19:647–55.

    Article  CAS  PubMed  Google Scholar 

  28. Liang G, He H, Yu D. Identification of nitrogen starvation-responsive microRNAs in Arabidopsis thaliana. PLoS One. 2012;7:e48951.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Zhao M, Tai H, Sun S, Zhang F, Xu Y, Li WX. Cloning and characterization of maize miRNAs involved in responses to nitrogen deficiency. PLoS One. 2012;7:e29669.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Zhao M, Ding H, Zhu JK, Zhang F, Li WX. Involvement of miR169 in the nitrogen-starvation responses in Arabidopsis. New Phytol. 2011;190:906–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Yan Y, Wang H, Hamera S, Chen X, Fang R. MiR444a has multiple functions in the rice nitrate-signaling pathway. Plant J. 2014;78:44–55.

    Article  CAS  PubMed  Google Scholar 

  32. Abdel-Ghany SE, Pilon M. MicroRNA-mdiated systemic down-regulation of copper protein expression in response to low copper availability in Arabidopsis. J Biol Chem. 2008;283:15932–45.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Waters BM, Mclnturf SA, Stein RJ. Rosette iron deficiency transcript and microRNA profiling reveals links between copper and iron homeostasis in Arabidopsis thaliana. J Exp Bot. 2012;63:5903–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Liang G, Yang F, Yu D. MicroRNA395 mediates regulation of sulfate accumulation and allocation in Arabidopsis thaliana. Plant J. 2010;62:1046–57.

    CAS  PubMed  Google Scholar 

  35. Li Y, Zhang Y, Shi D, Liu X, Qin J, Ge Q, Xu L, Pan X, Li W, Zhu Y, Xu J. Spatial-temporal analysis of zinc homeostasis reveals the response mechanisms to acute zinc deficiency in Sorghum bicolor. New Phytol. 2013;200:1102–15.

    Article  CAS  PubMed  Google Scholar 

  36. Ma CL, Qi YP, Yang LT, Lu YB, Guo P, Ye X, Chen LS. MicroRNA regulatory mechanisms on Citrus sinensis leaves to magnesium-deficiency. Front Plant Sci. 2016;7:201.

    PubMed  PubMed Central  Google Scholar 

  37. Chapman HD. The mineral nutrition of Citrus. In: Reuther W, Webber HJ, Batchelor LD, editors. The Citrus industry, vol. 2. CA: Division of Agricultural Sciences, University of California; 1968. p. 127–89.

    Google Scholar 

  38. Xu Q, Liu Y, Zhu A, Wu X, Ye J, Yu K, Guo W, Deng X. Discovery and comparative and comparative profiling of microRNAs in a sweet orange red-flesh mutant and its wild type. BMC Genomics. 2010;11:246.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Song C, Wang C, Zhang C, Korir NK, Yu H, Ma Z, Fang J. Deep sequencing discovery of novel and conserved microRNAs in trifoliate orange (Citrus trifoliata). BMC Genomics. 2010;11:431.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Li R, Li Y, Kristiansen K, Wang J. SOAP: short oligonucleotide alignment program. Bioinformatics. 2008;24:713–4.

    Article  CAS  PubMed  Google Scholar 

  41. Chen L, Wang T, Zhao M, Tian Q, Zhang WH. Identification of aluminum-responsive microRNAs in Medicago truncatula by genome-wide high-throughput sequencing. Planta. 2012;235:375–86.

    Article  CAS  PubMed  Google Scholar 

  42. Huang JH, Qi YP, Wen SX, Guo P, Chen XM, Chen LS. Illumina microRNA profiles reveal the involvement of miR397a in Citrus adaptation to long-term boron toxicity via modulating secondary cell-wall biosynthesis. Sci Rep. 2016;6:22900.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Hermans C, Conn SJ, Chen J, Xiao Q, Verbruggen N. An update on magnesium homeostasis mechanisms in plants. Metallomics. 2013;5:1170–83.

    Article  CAS  PubMed  Google Scholar 

  44. Wang Y, Ries A, Wu K, Yang A, Crawford NM. The Arabidopsis prohibitin gene PHB3 functions in nitric oxide-mediated responses and in hydrogen peroxide-induced nitric oxide accumulation. Plant Cell. 2010;22:249–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Gehl B, Lee CP, Bota P, Blatt MR, Sweetlove LJ. An Arabidopsis stomatin-like protein affects mitochondrial respiratory supercomplex organization. Plant Physiol. 2014;164:1389–400.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Xu J. Effects of magnesium deficiency on Citrus elements and root and leaf anatomical structures, Master thesis. Fuzhou: Fujian Agriculture and Forestry University, Fuzhou; 2015.

    Google Scholar 

  47. Amtmann A, Troufflard S, Armengaud P. The effect of potassium nutrition on pest and disease resistance in plants. Physiol Plant. 2008;133:682–91.

    Article  CAS  PubMed  Google Scholar 

  48. Huber DM, Jones JB. The role of magnesium in plant disease. Plant Soil. 2013;368:73–85.

    Article  CAS  Google Scholar 

  49. Guleria P, Yadav SK. Identification of miR414 and expression analysis of conserved miRNAs from Stevia rebaudiana. Genomics Proteomics Bioinformatics. 2011;9:211–7.

    Article  CAS  PubMed  Google Scholar 

  50. Galles C, Spampinato CP. Yeast mutator phenotype enforced by Arabidopsis PMS1 expression. Mol Biol Rep. 2013;40:2107–14.

    Article  CAS  PubMed  Google Scholar 

  51. Dion E, Li L, Jean M, Belzile F. An Arabidopsis MLH1 mutant exhibits reproductive defects and reveals a dual role for this gene in mitotic recombination. Plant J. 2007;51:431–40.

    Article  CAS  PubMed  Google Scholar 

  52. Dievart A, Perin C, Hirsch J, Bettembourg M, Lanau N, Artus F, Bureau C, Noel N, Droc G, Peyramard M, Pereira S, Courtois B, Morel JB, Guiderdoni E. The phenome analysis of mutant alleles in leucine-rich repeat receptor-like kinase genes in rice reveals new potential targets for stress tolerant cereals. Plant Sci. 2016;242:240–9.

    Article  CAS  PubMed  Google Scholar 

  53. Rodríguez-Villalón A, Gas E, Rodríguez-Concepción M. Phytoene synthase activity controls the biosynthesis of carotenoids and the supply of their metabolic precursors in dark-grown Arabidopsis seedlings. Plant J. 2009;60:424–35.

    Article  PubMed  Google Scholar 

  54. Zhang J, Tao N, Xu Q, Zhou W, Cao H, Xu J, Deng X. Functional characterization of Citrus PSY gene in Hongkong kumquat (Fortunella hindsii Swingle). Plant Cell Rep. 2009;28:1737–46.

    Article  CAS  PubMed  Google Scholar 

  55. Nakano M, Nishihara M, Yoshioka H, Takahashi H, Sawasaki T, Ohnishi K, Hikichi Y, Kiba A. Suppression of DS1 phosphatidic acid phosphatase confirms resistance to Ralstonia solanacearum in Nicotiana benthamiana. PLoS One. 2013;8:e75124.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Duque P. A role for SR proteins in plant stress responses. Plant Signal Behav. 2011;6:49–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Hurley BA, Tran HT, Marty NJ, Park J, Snedden WA, Mullen RT, Plaxton WC. The dual-targeted purple acid phosphatase isozyme AtPAP26 is essential for efficient acclimation of Arabidopsis to nutritional phosphate deprivation. Plant Physiol. 2010;153:1112–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Veljanovski V, Vanderbeld B, Knowles VL, Snedden WA, Plaxton WC. Biochemical and molecular characterization of AtPAP26, a vacuolar purple acid phosphatase up-regulated in phosphate-deprived Arabidopsis suspension cells and seedlings. Plant Physiol. 2006;142:1282–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Malhotra V, Erlmann P. Protein export at the ER: loading big collagens into COPII carriers. EMBO J. 2011;30:3475–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Yang LT, Qi YP, Lu YB, Guo P, Sang W, Feng H, Zhang HX, Chen LS. iTRAQ protein profile analysis of Citrus sinensis roots in response to long-term boron-deficiency. J Proteome. 2013;93:179–206.

    Article  CAS  Google Scholar 

  61. Mallory AC, Bartel DP, Bartel B. MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes. Plant Cell. 2005;17:1360–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Wang JW, Wang LJ, Mao YB, Cai WJ, Xue HW, Chen XY. Control of root cap formation by microRNA-targeted auxin response factors in Arabidopsis. Plant Cell. 2005;17:2204–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Liu PP, Montgomery TA, Fahlgren N, Kasschau KD, Nonogaki H, Carrington JC. Repression of AUXIN RESPONSE FACTOR10 by microRNA160 is critical for seed germination and post-germination stages. Plant J. 2007;52:133–46.

    Article  CAS  PubMed  Google Scholar 

  64. Li XP, Zeng QF, Zhang GS, Zhao J. GmMIR160A, a class of soybean microRNA gene, negatively regulates progress of leaf senescence. Guihaia. 2015;35:84–91.

    Google Scholar 

  65. Zhu YY, Zeng HQ, Dong CX, Yin XM, Shen QR, Yang ZM. MicroRNA expression profiles associated with phosphorus deficiency in white lupin (Lupinus albus L.). Plant Sci. 2010;178:23–9.

    Article  CAS  Google Scholar 

  66. Xu Z, Zhong S, Li X, Li W, Rothstein SJ, Zhang S, Bi Y, Xie C. Genome-wide identification of microRNAs in response to low nitrate availability in maize leaves and roots. PLoS One. 2011;6:e28009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Franco-Zorrilla JM1, Valli A, Todesco M, Mateos I, Puga MI, Rubio-Somoza I, Leyva A, Weigel D, García JA, Paz-Ares J. Target mimicry provides a new mechanism for regulation of microRNA activity. Nat Genet. 2007;39:1033–7.

    Article  CAS  PubMed  Google Scholar 

  68. Wu HJ, Wang ZM, Wang M, Wang XJ. Widespread long noncoding RNAs as endogenous target mimics for microRNAs in plants. Plant Physiol. 2013;161:1875–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Lin Y, Lai Z, Tian Q, Lin L, Lai R, Yang M, Zhang D, Chen Y, Zhang Z. Endogenous target mimics down-regulate miR160 mediation of ARF10, −16, and −17 cleavage during somatic embryogenesis in Dimocarpus longan lour. Front Plant Sci. 2015;6:956.

    PubMed  PubMed Central  Google Scholar 

  70. Li J, Dai X, Zhao Y. A role for auxin response factor 19 in auxin and ethylene signaling in Arabidopsis. Plant Physiol. 2006;140:899–908.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  71. Okushima Y, Fukaki H, Onoda M, Theologis A, Tasaka M. ARF7 and ARF19 regulate lateral root formation via direct activation of LBD/ASL genes in Arabidopsis. Plant Cell. 2007;19:118–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Ohtani M, Demura T, Sugiyama M. Particular significance of SRD2-dependent snRNA accumulation in polarized pattern generation during lateral root development of Arabidopsis. Plant Cell Physiol. 2010;51:2002–12.

    Article  CAS  PubMed  Google Scholar 

  73. Cheng YT, Germain H, Wiermer M, Bi D, Xu F, García AV, Wirthmueller L, Després C, Parker JE, Zhang Y, Li X. Nuclear pore complex component MOS7/Nup88 is required for innate immunity and nuclear accumulation of defense regulators in Arabidopsis. Plant Cell. 2009;21:2503–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Lisso J, Altmann T, Müssig C. The AtNFXL1 gene encodes a NF-X1 type zinc finger protein required for growth under salt stress. FEBS Lett. 2006;580:4851–6.

    Article  CAS  PubMed  Google Scholar 

  75. Park S, Rancour DM, Bednarek SY. In planta analysis of the cell cycle-dependent localization of AtCDC48A and its critical roles in cell division, expansion, and differentiation. Plant Physiol. 2008;148:246–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  76. Wang H, Jin S, Chen X, Gen XK, He YK. Target deletion of the AAA ATPase PpCDC48II in Physcomitrella patens results in freezing sensitivity after cold acclimation. Sci China Life Sci. 2012;55:150–7.

    Article  CAS  PubMed  Google Scholar 

  77. Šamajová O, Plíhal O, Al-Yousif M, Hirt H, Šamaj J. Improvement of stress tolerance in plants by genetic manipulation of mitogen-activated protein kinases. Biotechnol Adv. 2013;31:118–28.

    Article  PubMed  Google Scholar 

  78. Xing Y, Chen WH, Jia W, Zhang J. Mitogen-activated protein kinase kinase 5 (MKK5)-mediated signalling cascade regulates expression of iron superoxide dismutase gene in Arabidopsis under salinity stress. J Exp Bot. 2015;66:5971–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Liu D, Li A, Mao X, Jing R. Cloning and characterization of TaPP2AbB”-α, a member of the PP2A regulatory subunit in wheat. PLoS One 2014;9:e94430.

  80. Zhou GA, Chang RZ, Qiu LJ. Overexpression of soybean ubiquitin-conjugating enzyme gene GmUBC2 confers enhanced drought and salt tolerance through modulating abiotic stress-responsive gene expression in Arabidopsis. Plant Mol Biol. 2010;72:357–67.

    Article  CAS  PubMed  Google Scholar 

  81. Liu J, Xia Z, Wang M, Zhang X, Yang T, Wu J. Overexpression of a maize E3 ubiquitin ligase gene enhances drought tolerance through regulating stomatal aperture and antioxidant system in transgenic tobacco. Plant Physiol Biochem. 2013;73:114–20.

    Article  CAS  PubMed  Google Scholar 

  82. Zhou S, Sun X, Yin S, Kong X, Zhou S, Xu Y, Luo Y, Wang W. The role of the F-box gene TaFBA1 from wheat (Triticum aestivum L.) in drought tolerance. Plant Physiol Biochem. 2014;84:213–23.

    Article  CAS  PubMed  Google Scholar 

  83. Zhou SM, Kong XZ, Kang HH, Sun XD, Wang W. The involvement of wheat F-box protein gene TaFBA1 in the oxidative stress tolerance of plants. PLoS One. 2015;10:e012211.

    Google Scholar 

  84. Kushizaki M. An extraction procedure of plant materials for the rapid determination of Mn, cu, Zn and mg by the atomic absorption analysis. J Sci Soil Manure Japan. 1968;39:489–90.

    CAS  Google Scholar 

  85. Zhang Z, Jiang L, Wang J, Gu P, Chen M. MTide: an integrated tool for the identification of miRNA-target interaction in plants. Bioinformatics. 2015;31:290–1.

    Article  CAS  PubMed  Google Scholar 

  86. Wang T, Chen L, Zhao M, Tian Q, Zhang WH. Identification of drought-responsive microRNAs in Medicago truncatula by genome-wide high throughout sequencing. BMC Genomics. 2011;12:367.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M. Specific effects of microRNAs on the plant transcriptome. Dev Cell. 2005;8:517–27.

    Article  CAS  PubMed  Google Scholar 

  88. Allen E, Xie Z, Gustafson AM, Carrington JC. MicroRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005;121:207–21.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by the National Natural Science Foundation of China (No. 31572081), the Natural Science Foundation of Fujian Province, China (No. 2016 J01115) and the earmarked fund for China Agriculture Research System (No. CARS-27). The funding bodies had no role in the design, collection and interpretation of data or the manuscript writing.

Availability of data and materials

All data generated or analyzed during this study are included in this published article [and its supplementary information files].

Primary RNA-Seq data are submitted to Gene Expression Omnibus (GEO) under accession number GSE80713 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE80713).

Author information

Authors and Affiliations

Authors

Contributions

WWL carried out most of the experiment, analyzed the data and drafted the manuscript; JHH and LTY participated in the direction of this study; CPL and XY participated in the measurement of Mg and data analysis; DL participated in the cultivation of seedlings and qRT-PCR analysis; LSC designed and directed the study and revised the manuscript. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Li-Song Chen.

Ethics declarations

Ethics approval and consent to participate

Citrus sinensis seeds, which were public and available for non-commercial purpose, were collected from Minan village, Tingjiang town, Mawei district, Fuzhou city, China. Collection of seeds complied with the institutional, national and international guidelines. No specific permits were required.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

List of known miRNAs in C. sinensis roots. (DOC 995 kb)

Additional file 2:

List of known miRNAs in C. sinensis roots after removing these miRNAs with normalized read-count less than 10 TPM in the two miRNA libraries constructed from Mg-sufficient and -deficient roots. (DOC 447 kb)

Additional file 3:

List of Mg-deficiency-responsive known miRNAs in C. sinensis roots. (DOC 205 kb)

Additional file 4:

List of novel miRNAs in C. sinensis roots. (DOC 288 kb)

Additional file 5:

List of novel miRNAs in C. sinensis roots after removing these miRNAs with normalized read-count less than 10 TPM in the two miRNA libraries constructed from Mg-sufficient and -deficient roots. (DOC 160 kb)

Additional file 6:

List of Mg-deficiency-responsive novel miRNAs in C. sinensis roots. (DOC 105 kb)

Additional file 7:

List of target genes for parts of known miRNAs in C. sinensis roots. (DOC 262 kb)

Additional file 8:

List of target genes for parts of novel miRNAs in C. sinensis roots. (DOCX 25 kb)

Additional file 9:

Stem loop primer sequences for qRT-PCR analysis of miRNAs. (DOCX 17 kb)

Additional file 10:

Specific primer pairs used for qRT-PCR expression analysis of selected miRNA target genes. (DOCX 31 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liang, WW., Huang, JH., Li, CP. et al. MicroRNA-mediated responses to long-term magnesium-deficiency in Citrus sinensis roots revealed by Illumina sequencing. BMC Genomics 18, 657 (2017). https://doi.org/10.1186/s12864-017-3999-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-017-3999-5

Keywords