Genome-Wide Identification of Candidate Genes Associated with Heat Stress in Mulberry (Morus alba L.)

Abstract

Mulberry (Morus alba L.) is an economically important plant for the silk industry and has the possibility of contributing immensely to Chinese pharmacopeia because of its health benefits. Domesticated silkworms feed only on mulberry leaves, meaning that the worms’ survival depends on the mulberry tree. Mulberry production is threatened by climate change and global warming. However, the regulatory mechanisms of mulberry responses to heat are poorly understood. We performed transcriptome analysis of high-temperature-stressed (42 °C) M. alba seedlings using RNA-Seq technologies. A total of 703 differentially expressed genes (DEGs) were discovered from 18,989 unigenes. Among these, 356 were up-regulated, and 347 were down-regulated. KEGG analysis revealed that most DEGs were enriched in valine, leucine and isoleucine degradation, and in starch and sucrose metabolism, alpha-linolenic acid metabolism, carotenoid biosynthesis and galactose metabolism, among others. In addition, TFs such as the NAC, HSF, IAA1, MYB, AP2, GATA, WRKY, HLH and TCP families were actively involved in response to high temperatures. Moreover, we used RT-qPCR to confirm the expression changes of eight genes under heat stress observed in the RNA-Seq analysis. This study provides M. alba transcriptome profiles under heat stress and provides theoretical bases to researchers for better understanding mulberry heat response mechanisms and breeding heat-tolerant mulberry plants.