Mulberry is an economically significant crop for the sericulture industry worldwide. Stresses such as drought exposure have a significant influence on plant survival. Because metabolome directly reflects plant physiological condition, performing a global metabolomic analysis is one technique to examine this influence. Using a liquid chromatography-mass spectrometry (LC-MS) technique based on an untargeted metabolomic approach, the effect of drought stress on mulberry Yu-711 metabolic balance was examined. For this objective, Yu-711 leaves were subjected to two weeks of drought stress treatment and control without drought stress. Numerous differentially accumulated metabolic components in response to drought stress treatment were revealed by multivariate and univariate statistical analysis. Drought stress treatment (EG) revealed a more differentiated metabolite response than the control (CK). We found that the levels of total lipids, galactolipids, and phospholipids (PC, PA, PE) were significantly altered, producing 48% of the total differentially expressed metabolites. Fatty acyls components were the most abundant lipids expressed and decreased considerably by 73.6%. On the other hand, the prenol lipids class of lipids increased in drought leaves. Other classes of metabolites, including polyphenols (flavonoids and cinnamic acid), organic acid (amino acids), carbohydrates, benzenoids, and organoheterocyclic, had a dynamic trend in response to the drought stress. However, their levels under drought stress decreased significantly compared to the control. These findings give an overview for the understanding of global plant metabolic changes in defense mechanisms by revealing the mulberry plant metabolic profile through differentially accumulated compounds.

Intersexual differentiation is crucial for the speciation and maintenance of dioecious plants, but the underlying mechanisms, including the genes involved, are still poorly understood. Here, we focused on a typical dioicous plant Morus alba, to explore the molecular footprints relevant to sex evolution by revealing the differentially expressed genes (DEGs) between two sexes and the testing signals of selection for these DEGs. From the results, we found a total of 1543 DEGs. Interestingly, 333 and 66 genes expression were detected only in male and female inflorescences, respectively. Using comparative transcriptomics, the expression of 841 genes were found to be significantly higher in male than in female inflorescences and were mainly enriched in defense-related pathways including the biosynthesis of phenylpropanoids, cutin, suberine and waxes. Meanwhile, the expression of 702 genes was female-biased and largely enriched in pathways related to growth and development, such as carbohydrate metabolism, auxin signaling and cellular responses. In addition, 16.7% and 17.6% signals of selection were significantly detected in female- and male-biased genes, respectively, suggesting their non-negligible role in evolution. Our findings expanded the understanding of the molecular basis of intersexual differentiation and contribute to further research on sex evolution in dioecious plants.