Digital gene expression analysis of corky split vein caused by boron deficiency in 'Newhall' Navel Orange (Citrus sinensis Osbeck) for selecting differentially expressed genes related to vascular hypertrophy

Transcriptome Analysis in Pyrus betulaefolia Roots in Response to Short-Term Boron Deficiency

Boron (B) deficiency stress is frequently observed in pear orchards and causes a considerable loss of productivity and fruit quality. Pyrus betulaefolia is one of the most important rootstocks that has been widely used in pear production. The present study confirmed that the boron form of different tissues showed various changes, and the free boron content was significantly decreased under the short-term B deficiency condition. Moreover, the ABA and JA content also significantly accumulated in the root after short-term B deficiency treatment. A comprehensive transcriptome analysis of 24 h B deficiency treatment P. betulaefolia root was performed in this study. Transcriptome results revealed a total of 1230 up-regulated and 642 down-regulated differentially expressed genes (DEGs), respectively. B deficiency significantly increased the expression of the key aquaporin gene NIP5-1. In addition, B deficiency also increased the expression of ABA (ZEP and NCED) and JA (LOX, AOS and OPR) synthesis genes. Several MYB, WRKY, bHLH and ERF transcription factors were induced by B deficiency stress, which may relate to the regulation of B uptake and plant hormone synthesis. Overall, these findings suggested that P. betulaefolia root had adaptive responses to short-term B deficiency stress by improved boron absorption ability and hormone (JA and ABA) synthesis. The transcriptome analysis provided further information for understanding the mechanism of the pear rootstock responses to B deficiency stress.

The regulation of plant lignin biosynthesis under boron deficiency conditions

Boron (B) is a required micronutrient that is crucial for the growth and development of vascular plants. A deficiency in B is generally regarded as a limiting factor affecting agricultural production in many parts of the world. Boron is involved in the metabolism of plant lignin and additionally, B deficiency can lead to the excessive accumulation of lignin in plant leaves/roots, resulting in corking symptoms and inhibited growth. However, the effect of B on lignin biosynthesis is not as well characterized as the specific function of B in the cell wall. In this article, recent studies on the regulation of lignin biosynthesis in plants under low-B stress conditions are reviewed. Moreover, the following possible mechanisms underlying the lignin synthesis promoted by B deficiency are discussed: (1) the accumulation of phenolic substances during B deficiency directly enhances lignin synthesis; (2) excess H₂ O₂ has a dual function to the enhancement of lignin under boron deficiency conditions, serving as a substrate and a signaling molecule; and (3) B deficiency regulates lignin synthesis through the expression of genes encoding transcription factors such as MYBs. Finally, future studies regarding physiology, molecules, and transcriptional regulation may reveal the mechanism(s) mediating the relationship between lignin synthesis and B deficiency. This review provides new insights and important references for future research and the enhancement of plant B nutrition.

JASMONATE RESISTANT 1 negatively regulates root growth under boron deficiency in Arabidopsis

Boron (B) is an essential micronutrient for plant growth and development. Jasmonic acid (JA) plays pivotal roles in plant growth, but the underlying molecular mechanism of JA involvement in B-deficiency-induced root growth inhibition is yet to be explored. In this study, we investigated the response of JA to B deficiency and the mechanism of JAR1-dependent JA signaling in root growth inhibition under B deficiency in Arabidopsis. B deficiency enhanced JA signaling in roots, and root growth inhibition was partially restored by JA biosynthesis inhibition. The jar1-1 (jasmonate-resistant 1, JAR1) mutant, and mutants of coronatine-insensitive 1 (coi1-2) and myc2 defective in JA signaling showed insensitivity to B deficiency. The ethylene-overproduction mutant eto1 and ethylene-insensitive mutant etr1 showed sensitivity and insensitivity to B deficiency, respectively, suggesting that ethylene is involved in the inhibition of primary root growth under B deficiency. Furthermore, after a decline in levels of EIN3, which may contribute to root growth, ethylene signaling was weakened in the jar1-1 mutant root under B deficiency. Under B deficiency, B concentrations were increased in the roots and shoots of the jar1-1 mutant, owing to the large root system and its activity. Therefore, our findings revealed that JA, which is involved in the inhibition of root growth under B deficiency, is regulated by JAR1-activated JA and ethylene signaling pathways.

Molecular mechanisms for magnesium-deficiency-induced leaf vein lignification, enlargement and cracking in Citrus sinensis revealed by RNA-Seq

Citrus sinensis (L.) Osbeck seedlings were fertigated with nutrient solution containing 2 [magnesium (Mg)-sufficiency] or 0 mM (Mg-deficiency) Mg(NO3)2 for 16 weeks. Thereafter, RNA-Seq was used to investigate Mg-deficiency-responsive genes in the veins of upper and lower leaves in order to understand the molecular mechanisms for Mg-deficiency-induced vein lignification, enlargement and cracking, which appeared only in the lower leaves. In this study, 3065 upregulated and 1220 downregulated, and 1390 upregulated and 375 downregulated genes were identified in Mg-deficiency veins of lower leaves (MDVLL) vs Mg-sufficiency veins of lower leaves (MSVLL) and Mg-deficiency veins of upper leaves (MDVUL) vs Mg-sufficiency veins of upper leaves (MSVUL), respectively. There were 1473 common differentially expressed genes (DEGs) between MDVLL vs MSVLL and MDVUL vs MSVUL, 1463 of which displayed the same expression trend. Magnesium-deficiency-induced lignification, enlargement and cracking in veins of lower leaves might be related to the following factors: (i) numerous transciption factors and genes involved in lignin biosynthesis pathways, regulation of cell cycle and cell wall metabolism were upregulated; and (ii) reactive oxygen species, phytohormone and cell wall integrity signalings were activated. Conjoint analysis of proteome and transcriptome indicated that there were 287 and 56 common elements between DEGs and differentially abundant proteins (DAPs) identified in MDVLL vs MSVLL and MDVUL vs MSVUL, respectively, and that among these common elements, the abundances of 198 and 55 DAPs matched well with the transcript levels of the corresponding DEGs in MDVLL vs MSVLL and MDVUL vs MSVUL, respectively, indicating the existence of concordances between protein and transcript levels.

From element to development: the power of the essential micronutrient boron to shape morphological processes in plants

Deficiency of the essential nutrient boron (B) in the soil is one of the most widespread micronutrient deficiencies worldwide, leading to developmental defects in root and shoot tissues of plants, and severe yield reductions in many crops. Despite this agricultural importance, the underlying mechanisms of how B shapes plant developmental and morphological processes are still not unequivocally understood in detail. This review evaluates experimental approaches that address our current understanding of how B influences plant morphological processes by focusing on developmental defects observed under B deficiency. We assess what is known about mechanisms that control B homeostasis and specifically highlight: (i) limitations in the methodology that is used to induce B deficiency; (ii) differences between mutant phenotypes and normal plants grown under B deficiency; and (iii) recent research on analyzing interactions between B and phytohormones. Our analysis highlights the need for standardized methodology to evaluate the roles of B in the cell wall versus other parts of the cell.

Genome-Wide Identification of Long Non-coding RNA in Trifoliate Orange (Poncirus trifoliata (L.) Raf) Leaves in Response to Boron Deficiency

Long non-coding RNAs (lncRNAs) play important roles in plant growth and stress responses. As a dominant abiotic stress factor in soil, boron (B) deficiency stress has impacted the growth and development of citrus in the red soil region of southern China. In the present work, we performed a genome-wide identification and characterization of lncRNAs in response to B deficiency stress in the leaves of trifoliate orange (Poncirus trifoliata), an important rootstock of citrus. A total of 2101 unique lncRNAs and 24,534 mRNAs were predicted. Quantitative real-time polymerase chain reaction (qRT-PCR) experiments were performed for a total of 16 random mRNAs and lncRNAs to validate their existence and expression patterns. Expression profiling of the leaves of trifoliate orange under B deficiency stress identified 729 up-regulated and 721 down-regulated lncRNAs, and 8419 up-regulated and 8395 down-regulated mRNAs. Further analysis showed that a total of 84 differentially expressed lncRNAs (DELs) were up-regulated and 31 were down-regulated, where the number of up-regulated DELs was 2.71-fold that of down-regulated. A similar trend was also observed in differentially expressed mRNAs (DEMs, 4.21-fold). Functional annotation of these DEMs was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, and the results demonstrated an enrichment of the categories of the biosynthesis of secondary metabolites (including phenylpropanoid biosynthesis/lignin biosynthesis), plant hormone signal transduction and the calcium signaling pathway. LncRNA target gene enrichment identified several target genes that were involved in plant hormones, and the expression of lncRNAs and their target genes was significantly influenced. Therefore, our results suggest that lncRNAs can regulate the metabolism and signal transduction of plant hormones, which play an important role in the responses of citrus plants to B deficiency stress. Co-expression network analysis indicated that 468 significantly differentially expressed genes may be potential targets of 90 lncRNAs, and a total of 838 matched lncRNA-mRNA pairs were identified. In summary, our data provides a rich resource of candidate lncRNAs and mRNAs, as well as their related pathways, thereby improving our understanding of the role of lncRNAs in response to B deficiency stress, and in symptom formation caused by B deficiency in the leaves of trifoliate orange.

Cytokinin at the Crossroads of Abiotic Stress Signalling Pathways

Cytokinin is a multifaceted plant hormone that plays major roles not only in diverse plant growth and development processes, but also stress responses. We summarize knowledge of the roles of its metabolism, transport, and signalling in responses to changes in levels of both macronutrients (nitrogen, phosphorus, potassium, sulphur) and micronutrients (boron, iron, silicon, selenium). We comment on cytokinin's effects on plants' xenobiotic resistance, and its interactions with light, temperature, drought, and salinity signals. Further, we have compiled a list of abiotic stress-related genes and demonstrate that their expression patterns overlap with those of cytokinin metabolism and signalling genes.

Boron deficiency inhibits root growth by controlling meristem activity under cytokinin regulation

Significant advances have been made in the last years trying to identify regulatory pathways that control plant responses to boron (B) deficiency. Still, there is a lack of a deep understanding of how they act regulating growth and development under B limiting conditions. Here, we analyzed the impact of B deficit on cell division leading to root apical meristem (RAM) disorganization. Our results reveal that inhibition of cell proliferation under the regulatory control of cytokinins (CKs) is an early event contributing to root growth arrest under B deficiency. An early recovery of QC46:GUS expression after transferring B-deficient seedlings to control conditions revealed a role of B in the maintenance of QC identity whose loss under deficiency occurred at later stages of the stress. Additionally, the D-type cyclin CYCD3 overexpressor and triple mutant cycd3;1-3 were used to evaluate the effect on mitosis inhibition at the G1-S boundary. Overall, this study supports the hypothesis that meristem activity is inhibited by B deficiency at early stages of the stress as it does cell elongation. Likewise, distinct regulatory mechanisms seem to take place depending on the severity of the stress. The results presented here are key to better understand early signaling responses under B deficiency.

The changes of intracellular calcium concentration and distribution in the hard end pear (Pyrus pyrifolia cv. 'Whangkeumbae') fruit

Hard end is a physiological disorder of pear fruit that is frequently observed in the 'Whangkeumbae' (Pyrus pyrifolia) variety, however, the mechanisms that are involved in its development are poorly understood. In this study, we explored the causes of hard end disorder in pear fruit in relation to calcium deficiency. During fruit development, the ratio of Ca/N, Ca/K, Ca/Mg and the content of B were significantly lower in the hard end fruit as compared to normal fruit. However, no calcium deficiency was detected in the soil and leaves of the orchard where the hard end fruit were located. Additionally, the Ca²⁺ influx in the calyx of hard end fruit was lower than that of normal fruit at 90 d after anthesis. The free Ca²⁺ and storage Ca²⁺ in the flesh cells of hard end fruit were less than that of normal fruit during fruit development, while an opposite tendency was observed at 120 d after anthesis (harvest day). In hard end fruit, the Ca²⁺ transport-related gene, PpCNGC1 (Cyclic nucleotide-gated ion channel 1), was up-regulated; whereas the Ca²⁺ sensor-related genes of PpCIPKs, PpCDPK28 and PpCML41 were all down-regulated. Spraying with a 2% calcium chloride (CaCl₂) solution inhibited the incidence rate of hard end disorder and decreased fruit firmness and lignin content during storage. Additionally, the ratio of Ca/N, Ca/K, Ca/Mg and the content of B all increased on harvest day. Our study suggests that low Ca²⁺ influx leads to less Ca²⁺ into the pear fruit, which results in an intracellular imbalance of Ca²⁺ and consequently triggers the development of hard end disorder.

Effect of boron deficiency on anatomical structure and chemical composition of petioles and photosynthesis of leaves in cotton (Gossypium hirsutum L.)

The main symptom of boron (B) deficiency in cotton is the formation of brown rings on leaf petioles. The objective of the present study was to determine the changes in the anatomical structure and chemical composition of petioles and photosynthesis of leaves in cotton under B deficiency. Compared to the control, B deficiency treatment resulted in large increases in the number of petioles with brown rings per plant (160.0%) and the number of rings on the petiole per functional leaf (711.1%) in cotton seedlings. The relative absorbance intensity in the fingerprint region of polysaccharide structure was decreased in petiole rings under B deficiency, while lignin contents were increased. Cotton plants mitigated the impairment of transport function in cotton petioles by increasing the areas of vascular bundles, phloem, xylem, and phloem fiber. Moreover, the stomatal conductance, photosynthetic rate, and transpiration rate in leaves were significantly decreased under B deficiency, thus impeding photosynthesis in cotton plants. Therefore, B deficiency reduces transport function in petioles and photosynthesis in leaves, and leads to the formation of noticeable brown rings on petioles of cotton seedlings.

Transcriptome Changes Associated with Boron Deficiency in Leaves of Two Citrus Scion-Rootstock Combinations

Boron (B) deficiency stress is frequently observed in citrus orchards and causes considerable loss of productivity and fruit quality. Carrizo citrange (Cc) has been reported as a rootstock more tolerant to B deficiency than Trifoliate orange (To). The 'Newhall' navel orange (Ns) performed better when grafted onto Cc (Ns/Cc) than when grafted onto To (Ns/To) under long-term B deficiency. The present study confirmed that Ns/Cc had higher boron content, leaf fresh weight, lower leaf chlorosis and stronger photosynthesis ability than Ns/To. Moreover, B-deficiency significantly reduced the chlorophyll and carotenoid content in Ns/To. The content of total soluble sugar and lignin were dramatically increased and the expression levels of photosynthesis-related genes were substantially down-regulated in Ns/To by B-deficient treatment. B-deficiency also strongly induced expression levels of chlorophyll decomposition-related genes, glucose synthesis-related genes and lignin synthesis-related genes, and significantly inhibited the expression of carotenoid synthesis-related genes in Ns/To. Overall, these findings suggested that the influence of To on the scion of Ns was worse than that of Cc due to differently regulating these metabolic pathways under the long term of B-deficiency. The transcriptome analysis provided further information for understanding the mechanism of the different responses of scion-rootstock combinations to B-deficiency stress.

Digital Gene Expression Analysis Provides Insight into the Transcript Profile of the Genes Involved in Aporphine Alkaloid Biosynthesis in Lotus (Nelumbo nucifera)

The predominant alkaloids in lotus leaves are aporphine alkaloids. These are the most important active components and have many pharmacological properties, but little is known about their biosynthesis. We used digital gene expression (DGE) technology to identify differentially-expressed genes (DEGs) between two lotus cultivars with different alkaloid contents at four leaf development stages. We also predicted potential genes involved in aporphine alkaloid biosynthesis by weighted gene co-expression network analysis (WGCNA). Approximately 335 billion nucleotides were generated; and 94% of which were aligned against the reference genome. Of 22 thousand expressed genes, 19,000 were differentially expressed between the two cultivars at the four stages. Gene Ontology (GO) enrichment analysis revealed that catalytic activity and oxidoreductase activity were enriched significantly in most pairwise comparisons. In Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, dozens of DEGs were assigned to the categories of biosynthesis of secondary metabolites, isoquinoline alkaloid biosynthesis, and flavonoid biosynthesis. The genes encoding norcoclaurine synthase (NCS), norcoclaurine 6-O-methyltransferase (6OMT), coclaurine N-methyltransferase (CNMT), N-methylcoclaurine 3'-hydroxylase (NMCH), and 3'-hydroxy-N-methylcoclaurine 4'-O-methyltransferase (4'OMT) in the common pathways of benzylisoquinoline alkaloid biosynthesis and the ones encoding corytuberine synthase (CTS) in aporphine alkaloid biosynthetic pathway, which have been characterized in other plants, were identified in lotus. These genes had positive effects on alkaloid content, albeit with phenotypic lag. The WGCNA of DEGs revealed that one network module was associated with the dynamic change of alkaloid content. Eleven genes encoding proteins with methyltransferase, oxidoreductase and CYP450 activities were identified. These were surmised to be genes involved in aporphine alkaloid biosynthesis. This transcriptomic database provides new directions for future studies on clarifying the aporphine alkaloid pathway.

Omics studies of citrus, grape and rosaceae fruit trees

Recent advance of bioinformatics and analytical apparatuses such as next generation DNA sequencer (NGS) and mass spectrometer (MS) has brought a big wave of comprehensive study to biology. Comprehensive study targeting all genes, transcripts (RNAs), proteins, metabolites, hormones, ions or phenotypes is called genomics, transcriptomics, proteomics, metabolomics, hormonomics, ionomics or phenomics, respectively. These omics are powerful approaches to identify key genes for important traits, to clarify events of physiological mechanisms and to reveal unknown metabolic pathways in crops. Recently, the use of omics approach has increased dramatically in fruit tree research. Although the most reported omics studies on fruit trees are transcriptomics, proteomics and metabolomics, and a few is reported on hormonomics and ionomics. In this article, we reviewed recent omics studies of major fruit trees, i.e. citrus, grapevine and rosaceae fruit trees. The effectiveness and prospects of omics in fruit tree research will as well be highlighted.

Boron-deficiency-responsive microRNAs and their targets in Citrus sinensis leaves

BACKGROUND

MicroRNAs play important roles in the adaptive responses of plants to nutrient deficiencies. Most research, however, has focused on nitrogen (N), phosphorus (P), sulfur (S), copper (Cu) and iron (Fe) deficiencies, limited data are available on the differential expression of miRNAs and their target genes in response to deficiencies of other nutrient elements. In this study, we identified the known and novel miRNAs as well as the boron (B)-deficiency-responsive miRNAs from citrus leaves in order to obtain the potential miRNAs related to the tolerance of citrus to B-deficiency.

METHODS

Seedlings of 'Xuegan' [Citrus sinensis (L.) Osbeck] were supplied every other day with B-deficient (0 μM H3BO3) or -sufficient (10 μM H3BO3) nutrient solution for 15 weeks. Thereafter, we sequenced two small RNA libraries from B-deficient and -sufficient (control) citrus leaves, respectively, using Illumina sequencing.

RESULTS

Ninety one (83 known and 8 novel) up- and 81 (75 known and 6 novel) down-regulated miRNAs were isolated from B-deficient leaves. The great alteration of miRNA expression might contribute to the tolerance of citrus to B-deficiency. The adaptive responses of miRNAs to B-deficiency might related to several aspects: (a) attenuation of plant growth and development by repressing auxin signaling due to decreased TIR1 level and ARF-mediated gene expression by altering the expression of miR393, miR160 and miR3946; (b) maintaining leaf phenotype and enhancing the stress tolerance by up-regulating NACs targeted by miR159, miR782, miR3946 and miR7539; (c) activation of the stress responses and antioxidant system through down-regulating the expression of miR164, miR6260, miR5929, miR6214, miR3946 and miR3446; (d) decreasing the expression of major facilitator superfamily protein genes targeted by miR5037, thus lowering B export from plants. Also, B-deficiency-induced down-regulation of miR408 might play a role in plant tolerance to B-deficiency by regulating Cu homeostasis and enhancing superoxide dismutase activity.

CONCLUSIONS

Our study reveals some novel responses of citrus to B-deficiency, which increase our understanding of the adaptive mechanisms of citrus to B-deficiency at the miRNA (post-transcriptional) level.

Boron deficiency in woody plants: various responses and tolerance mechanisms

Boron (B) is an essential microelement for higher plants, and its deficiency is widespread around the world and constrains the productivity of both agriculture and forestry. In the last two decades, numerous studies on model or herbaceous plants have contributed greatly to our understanding of the complex network of B-deficiency responses and mechanisms for tolerance. In woody plants, however, fewer studies have been conducted and they have not well been recently synthesized or related to the findings on model species on B transporters. Trees have a larger body size, longer lifespan and more B reserves than do herbaceous plants, indicating that woody species might undergo long-term or mild B deficiency more commonly and that regulation of B reserves helps trees cope with B deficiency. In addition, the highly heterozygous genetic background of tree species suggests that they may have more complex mechanisms of response and tolerance to B deficiency than do model plants. Boron-deficient trees usually exhibit two key visible symptoms: depression of growing points (root tip, bud, flower, and young leaf) and deformity of organs (root, shoot, leaf, and fruit). These symptoms may be ascribed to B functioning in the cell wall and membrane, and particularly to damage to vascular tissues and the suppression of both B and water transport. Boron deficiency also affects metabolic processes such as decreased leaf photosynthesis, and increased lignin and phenol content in trees. These negative effects will influence the quality and quantity of wood, fruit and other agricultural products. Boron efficiency probably originates from a combined effect of three processes: B uptake, B translocation and retranslocation, and B utilization. Root morphology and mycorrhiza can affect the B uptake efficiency of trees. During B translocation from the root to shoot, differences in B concentration between root cell sap and xylem exudate, as well as water use efficiency, may play key roles in tolerance to B deficiency. In addition, B retranslocation efficiency primarily depends on the extent of xylem-to-phloem transfer and the variety and amount of cis-diol moieties in the phloem. The B requirement for cell wall construction also contribute to the B use efficiency in trees. The present review will provide an update on the physiological and molecular responses and tolerance mechanisms to B deficiency in woody plants. Emphasis is placed on the roles of B reserves that are more important for tolerance to B deficiency in trees than in herbaceous plants and the possible physiological and molecular mechanisms of differential B efficiency in trees. We propose that B may be used to study the relationship between the cell wall and the membrane via the B-bridge. Transgenic B-efficient tree cultivars have considerable potential for forestry or fruit rootstock production on low B soils in the future.

miRNAome analysis associated with anatomic and transcriptomic investigations reveal the polar exhibition of corky split vein in boron deficient Citrus sinensis

Corky split vein can develop under long-term boron deficient conditions in Citrus sinensis L. Osbeck cv. Newhall. This symptom only occurs in the upper rather than the lower epidermis of old leaves. Our previous study demonstrated that vascular hypertrophy was involved in the symptoms, and the 3rd developmental stage of corky split vein (BD3) was the critical stage for phenotype formation. Here, we performed an intensive study on the BD3 vein and its control sample (CK3 vein). A lignin test demonstrated that the lignin content in BD3 vein was approximately 1.7 times more than the CK3 vein. Anatomical investigation of the corky split vein indicated that the upper epidermis was destroyed by overgrowing vascular cells, and the increased lignin may contribute to vascular cell differentiation and wounding-induced lignification. In a subsequent small RNA sequencing of the BD3 and CK3 veins, 99 known miRNAs and 22 novel miRNAs were identified. Comparative profiling of these miRNAs demonstrated that the 57 known miRNAs and all novel miRNAs exhibited significant expression differences between the two small RNAs libraries of the BD3 and CK3 veins. Associated with our corresponding digital gene expression data, we propose that the decreased expression of two miRNAs, csi-miR156b and csi-miR164, which leads to the up-regulation of their target genes, SPLs (csi-miR156b-targeted) and CUC2 (csi-miR164-targeted), may promote vascular cell division and orderless stage transition in old leaves.

Long-term boron-deficiency-responsive genes revealed by cDNA-AFLP differ between Citrus sinensis roots and leaves

Seedlings of Citrus sinensis (L.) Osbeck were supplied with boron (B)-deficient (without H3BO3) or -sufficient (10 μM H3BO3) nutrient solution for 15 weeks. We identified 54 (38) and 38 (45) up (down)-regulated cDNA-AFLP bands (transcript-derived fragments, TDFs) from B-deficient leaves and roots, respectively. These TDFs were mainly involved in protein and amino acid metabolism, carbohydrate and energy metabolism, nucleic acid metabolism, cell transport, signal transduction, and stress response and defense. The majority of the differentially expressed TDFs were isolated only from B-deficient roots or leaves, only seven TDFs with the same GenBank ID were isolated from the both. In addition, ATP biosynthesis-related TDFs were induced in B-deficient roots, but unaffected in B-deficient leaves. Most of the differentially expressed TDFs associated with signal transduction and stress defense were down-regulated in roots, but up-regulated in leaves. TDFs related to protein ubiquitination and proteolysis were induced in B-deficient leaves except for one TDF, while only two down-regulated TDFs associated with ubiquitination were detected in B-deficient roots. Thus, many differences existed in long-term B-deficiency-responsive genes between roots and leaves. In conclusion, our findings provided a global picture of the differential responses occurring in B-deficient roots and leaves and revealed new insight into the different adaptive mechanisms of C. sinensis roots and leaves to B-deficiency at the transcriptional level.

Identification of boron-deficiency-responsive microRNAs in Citrus sinensis roots by Illumina sequencing

BACKGROUND

Boron (B)-deficiency is a widespread problem in many crops, including Citrus. MicroRNAs (miRNAs) play important roles in nutrient deficiencies. However, little is known on B-deficiency-responsive miRNAs in plants. In this study, we first identified miRNAs and their expression pattern in B-deficient Citrus sinensis roots by Illumina sequencing in order to identify miRNAs that might be involved in the tolerance of plants to B-deficiency.

RESULTS

We isolated 52 (40 known and 12 novel) up-regulated and 82 (72 known and 10 novel) down-regulated miRNAs from B-deficient roots, demonstrating remarkable metabolic flexibility of roots, which might contribute to the tolerance of plants to B-deficiency. A model for the possible roles of miRNAs in the tolerance of roots to B-deficiency was proposed. miRNAs might regulate the adaptations of roots to B-deficiency through following several aspects: (a) inactivating reactive oxygen species (ROS) signaling and scavenging through up-regulating miR474 and down-regulating miR782 and miR843; (b) increasing lateral root number by lowering miR5023 expression and maintaining a certain phenotype favorable for B-deficiency-tolerance by increasing miR394 expression; (c) enhancing cell transport by decreasing the transcripts of miR830, miR5266 and miR3465; (d) improving osmoprotection (miR474) and regulating other metabolic reactions (miR5023 and miR821). Other miRNAs such as miR472 and miR2118 in roots increased in response to B-deficiency, thus decreasing the expression of their target genes, which are involved in disease resistance, and hence, the disease resistance of roots.

CONCLUSIONS

Our work demonstrates the possible roles of miRNAs and related mechanisms in the response of plant roots to B-deficiency.