Comparative Transcriptome Analysis Identifies Key Defense Genes and Mechanisms in Mulberry (Morus alba) Leaves against Silkworms (Bombyx mori)

As a consequence of long-term coevolution and natural selection, the leaves of mulberry (Morus alba) trees have become the best food source for silkworms (Bombyx mori). Nevertheless, the molecular and genomic basis of defense response remains largely unexplored. In the present study, we assessed changes in the transcriptome changes of mulberry in response to silkworm larval feeding at 0, 3, and 6 h. A total of 4709 (up = 2971, down = 1738) and 3009 (up = 1868, down = 1141) unigenes were identified after 3 and 6 h of silkworm infestation, respectively. MapMan enrichment analysis results show structural traits such as leaf surface wax, cell wall thickness and lignification form the first physical barrier to feeding by the silkworms. Cluster analysis revealed six unique temporal patterns of transcriptome changes. We predicted that mulberry promoted rapid changes in signaling and other regulatory processes to deal with mechanical damage, photosynthesis impairment, and other injury caused by herbivores within 3–6 h. LRR-RK coding genes (THE1, FER) was predicted participated in perception of cell wall perturbation in mulberry responding to silkworm feeding. Ca(2+) signal sensors (CMLs), ROS (OST1, SOS3), RBOHD/F, CDPKs, and ABA were part of the regulatory network after silkworm feeding. Jasmonic acid (JA) signal transduction was predicted to act in silkworm feeding response, 10 JA signaling genes (such as OPR3, JAR1, and JAZ1) and 21 JA synthesis genes (such as LOX2, AOS, and ACX1) were upregulated after silkworm feeding for 3 h. Besides, genes of “alpha-Linolenic acid metabolism” and “phenylpropanoid biosynthesis” were activated in 3 h to reprogram secondary metabolism. Collectively, these findings provided valuable insights into silkworm herbivory-induced regulatory and metabolic processes in mulberry, which might help improve the coevolution of silkworm and mulberry.