Comparative physiological and transcriptome analysis reveals the potential mechanism of selenium accumulation and tolerance to selenate toxicity of Broussonetia papyrifera

Broussonetia papyrifera is an important fodder tree that is widely distributed in China. Enhancing the selenium (Se) content in B. papyrifera may help to improve the nutritional value of the feed. In this study, sodium selenite and selenate were foliar applied to investigate the mechanisms of Se tolerance and accumulation in B. papyrifera. The results showed that both Se forms significantly increased the total Se content, and the proportion of organic Se was significantly higher in the sodium selenite treatment than in the control. In addition, the soluble sugar, phenolic acid and flavonoid contents and antioxidant enzyme activities were increased by exogenous Se. The de novo RNA sequencing results showed that 644 and 1804 differentially expressed genes were identified in the selenite and selenate comparison groups, respectively. Pathway enrichment analysis demonstrated that 24 of the 108 pathways were significantly enriched, of which sulfur assimilation genes in the sodium selenite-treated groups were upregulated, whereas Se conjugation and transporter genes, such as SBP1, PCS, GSTs, ABCs and GPX, were significantly induced under selenate treatment. The hub genes identified by weighted-gene co-expression network analysis further confirmed that sulfur assimilation, conjugation and transporter genes might play a vital role in Se assimilation and tolerance. From this, a model of Se metabolism in B. papyrifera was proposed based on the above physiological and RNA sequencing data. This study is the first study to report that B. papyrifera has a strong ability to accumulate and tolerate exogenous Se, thereby providing a foundation for further characterization of the accumulation and tolerance mechanism of B. papyrifera. Our findings can provide technical support for producing Se-enriched fodder.