Transcriptome Analysis Reveals the Heat Stress Response Genes by Fire Stimulation in Michelia macclurei Dandy

Abstract

Heat stress due to external heat sources such as fire is an ecological problem for plants. When forest plants suffer from fire, high temperatures cause an array of morphological, physiological, and biochemical changes, which affect growth and development. Michelia macclurei Dandy is an evergreen broad-leaved tree species with the characteristics of fast growth, strong adaptability, and good fire-resistance. Some studies have improved the understanding of how fire behavior affects physiology, function and mortality, but the extreme heat response genes and mechanisms need improved understanding. In this study, we conducted a fire experiment (slight and severe) and RNA-Seq in M. macclure. The de novo assembly obtained 104,052 unigenes, and 48.46% were annotated in at least one public database. Specifically, 4458 and 4810 differentially expressed genes (DEGs) were identified in slight and severe fire treatment groups, respectively. In two treatment groups, 612 unigenes were differentially expressed, which were enriched in ‘oxidoreductase activity’ in the molecular function (MF) category of Gene Ontology (GO) enrichment analysis, suggesting the core role of oxidoreductase activity in response to extremely high temperatures in M. macclurei. In KEGG enrichment analysis of DEGs, the ‘plant hormone signal transduction’ is overrepresented, suggesting that this process plays an important role during heat response in M. macclurei. In the pathways of cytokinine and salicylic acid, some vital DEGs were enriched, which were related to cell division, shoot initiation, and disease resistance, and the potential interactions during heat stress were discussed. Moreover, the DEGs linked to heat stress response were identified, including heat shock factors, stress enhanced protein, signal transduction, photosystem, and major transcription factors. The qRT-PCR examination of various tissues, expression dynamics, and treatments revealed that the genes coding for the heat shock protein HSF30, stress enhanced protein, and photosystem I reaction center subunit II exhibited particularities in leaf tissue. Genes coding for heat shock proteins displayed a distinct expression pattern between fire treatment and conventional heat stress, which could signify the distinctive function of HSPs and the mechanism of heat responses. Altogether, these may interact to respond to fire stress through alterations in cellular processes, signaling transduction, and the synthesis and degradation of response proteins in M. macclurei. The results of this study provide a crucial transcriptional profile influenced by heat stress in M. macclurei, and could be of great use to explore the fire prevention mechanisms of fire-resistant tree species.