Genome-wide transcriptome profiling of mulberry (Morus alba) response to boron deficiency and toxicity reveal candidate genes associated with boron tolerance in leaves

Mulberry (Morus alba) is an essential plant with countless economic benefits; however, its growth and metabolic processes are hampered by boron (B) stresses. Very little research has been performed to elucidate boron tolerance and detoxification mechanisms in this species. The M. alba cultivar, Yu-711, was exposed to five different concentrations of boric acid (H3BO3), including deficient (T1; 0 mM) moderate B deficiency (T2; 0.02 mM), sufficient (CK; 0.1 mM) and toxic (T3 and T4; 0.5 and 1 mM) levels for 18 days of growth in pots experiment. Transcriptome analysis of B deficiency and toxicity treatments was performed on mulberry leaves. The transcriptome data reveal that a total of 6114 genes were differentially expressed (DEGs), of which 3830 were up-regulated and 2284 were down-regulated. A comparative analysis between treatment groups CK-vs-T1 (deficiency) and CK-vs-T4 (toxicity) indicates that 590 and 1383 genes were down-regulated in both deficiency and B toxicity, respectively. The results show that 206 genes were differentially expressed in all treatments. B deficiency and toxicity significantly altered the expression of the key aquaporins (PIP2-1, PIP2-7, PIP2-4 and NIP3-1) and high-affinity boron transporter genes (BOR1 and BOR7). In addition, boron stress also altered the expression of antioxidants and photosynthesis-related genes. B stresses were found to alter several transcription factors including ERF1B, which is associated with the regulation of boron uptake and the synthesis and signaling of phytohormones. Unravelling the mechanisms of B tolerance and detoxification is important and would give us further insight into how B stresses affect mulberry plants.