Reviews in Agricultural Science
Online ISSN : 2187-090X
Transcriptional Regulation and Gene Expression of Disease-Responsive Genes in Model Plants and Hevea Brasiliensis
Afdholiatus SyafaahOkechukwu S. EzehYoshiharu Y. Yamamoto
著者情報
ジャーナル フリー HTML

2026 年 14 巻 2 号 p. 74-94

詳細
ABSTRACT

Plant immunity relies on coordinated transcriptional reprogramming that integrates pathogen perception, hormone signaling, and redox homeostasis to balance effective defense with growth and development. While these mechanisms have been extensively characterized in model plants, their conservation and diversification in woody perennial crops remain less well understood. Hevea brasiliensis, the primary source of natural rubber, represents a biologically important tree crop in which immune regulation must be tightly coordinated with specialized laticifer metabolism and long life cycles. Here, we review current advances in understanding the transcriptional regulation of immune responses in Hevea, with a focus on major transcription factor (TF) families (WRKY, AP2/ERF, NAC, and MYB) and cis-regulatory elements conserved in defense regulation in model plants. We additionally discuss pathogenesis-related genes and the roles of reactive oxygen species (ROS) biosynthesis and detoxification genes in Hevea immunity. Transcriptomic evidence indicates that Hevea activates core immune modules shared with model plants, while also exhibiting additional layers of distinctive regulatory complexity, including tissue-specific gene expression in laticifers during wounding and pathogen infection. We conclude by discussing current limitations in functional validation and highlight future directions integrating RNA-seq, cis-regulatory elements analysis, and targeted functional genomics to elucidate biotic stress networks. Such advances will be critical for accelerating molecular breeding strategies aimed at improving disease resistance while preserving natural rubber productivity.

1. Introduction

Biotic stressors such as fungi, bacteria, and viruses can infect plants, inhibit plant growth, and reduce yield. To protect against pathogen invasion, plants have developed complex defense systems, controlled by many genes that switch on in response to biotic stress. In the first layer of defense, plant immune responses to pathogens depend on the pathogen-associated molecular patterns (PAMPs) perceived by specific pattern recognition receptors (PRRs). This detection initiates a defense mechanism known as PAMP-triggered immunity (PTI) [1] while the second layer is the effector-triggered immunity (ETI) [2, 3]. The host cells activate transcriptional reprogramming including dynamic regulation of signaling networks, transcription factors (TFs), and effector genes [4]. Advances in high-throughput studies have shown the molecular basis of plant immunity in model plants, such as Arabidopsis and rice [5, 6]. Several reports consistently show that infection triggers the activation of defense-related genes, such as pathogenesis-related (PR) proteins, reactive oxygen species (ROS)-related proteins, and hormone-responsive regulators [7, 8, 9].

Transcriptomic studies examining the responses of tropical tree species to pathogen infection are beginning to emerge, yet a substantial knowledge gap remains compared with the deep understanding obtained from model and non-woody plants. Hevea brasiliensis, a major economic tree crop valued for latex production, is similarly exposed to diverse biotic challenges, but its defense regulatory mechanisms remain poorly characterized. Consequently, mechanistic interpretations in Hevea are often inferred based on findings from related or better-studied species rather than from direct experimental evidence. Recently, transcriptomic studies in Hevea reveal strong activation of WRKY, NAC, ERF, and MYB transcription factors, early ROS accumulation, hormone crosstalk, and induction of PR proteins [10, 11]. As a perennial species, Hevea brasiliensis experiences constant environmental fluctuations, requiring flexible and stimulus-dependent modulation of its defense responses when confronted with pathogen attack. Understanding how transcriptional regulation in Hevea compares with that of model plants is critical for breeding disease-resistant rubber and for elucidating immune mechanisms in woody plants. This review provides current information on the transcriptional regulation and gene expression of disease-responsive genes in model plants and Hevea brasiliensis, and proposes research priorities for the molecular understanding of defense regulation in this important tropical tree crop.

2. Transcriptional regulation during pathogen infection

Transcription factors regulate gene expression by binding to specific cis-regulatory elements in promoter regions, thereby modulating the expression of downstream regulatory and defense-related genes. This coordinated interaction establishes an integrated transcriptional network essential for plant stress responses.

Comparative transcriptomic analyses of two Hevea clones following infection with Corynespora cassiicola revealed distinct gene expression and transcription factor regulatory patterns associated with contrasting levels of disease susceptibility [12]. Notably, genes encoding several transcription factor families, including WRKY, NAC, and GATA, were upregulated during disease development in both clones [12], suggesting their involvement in the core defense response. These transcription factor families are well characterized in model plant systems, where they play central roles in coordinating immune signaling and defense gene expression during pathogen attack. Their conserved induction in Hevea brasiliensis during C. cassiicola infection supports the notion that core transcriptional regulatory frameworks governing biotic stress responses are evolutionary conserved between herbaceous model plants and woody perennial species.

2.1 Transcriptional regulators of defense

Plant defense signaling triggers transcriptional reprogramming mediated by transcription factors that serve as key modulators of plant immune responses. During pathogen infection, defense responses are governed by complex regulatory networks involving multiple signaling cascades and dynamic metabolic reprogramming. Among the transcription factor families involved, WRKY, AP2/ERF, NAC, and MYB are frequently implicated in the regulation of immune responses to biotic stress [13, 14, 15, 16].

2.1.1 WRKY transcription factors

The WRKY transcription factors constitute a prominent family of plant-specific regulatory proteins marked by the conserved WRKY domains and the zinc-finger motif [17]. Members are classified into three groups according to the number of WRKY domains and the nature of their zinc finger motifs. Each WRKY domain consists of an N-terminal DNA-binding domain (DBD) harboring the conserved WRKYGQK heptapeptide and a C-terminal zinc finger motif that is essential for DNA binding and protein stability [18]. WRKY transcription factors can act as activators or repressors depending on the cellular and environmental context [19], and play pivotal roles in plant responses to abiotic and biotic stresses, as well as in development and senescence, primarily by binding to W-box (TTGACC/T) elements in the promoters of defense and stress-responsive genes [18, 20].

In Arabidopsis, AtWRKY18, AtWRKY40, and AtWRKY60 act cooperatively and partially redundantly as negative and positive regulators of resistance to both the hemibiotrophic bacterial pathogen Pseudomonas syringae and necrotrophic fungal pathogen Botrytis cinerea, with AtWRKY18 playing a more crucial role [21]. The contrasting effects of mutations of the three WRKY genes to the microbial pathogens reflect the antagonism between salicylic acid (SA)- and jasmonic acid (JA)-mediated defense signaling pathways [21]. Birkenbihl et al. [22] demonstrated that AtWRKY33 expression is strongly induced during B. cinerea infection. The atwrky33 mutant exhibits increased susceptibility to the pathogen, attributed to a reduction in the JA-regulated expression of PDF1.2 and impaired camalexin biosynthesis [22, 23]. Conversely, overexpression of AtWRKY33 enhanced the resistance against B. cinerea [23].

In maize (Zea mays), WRKY transcription factors comprise approximately 120 members, of which seven show significant transcriptional upregulation in response to infection by Ustilago maydis. These WRKY TFs specifically bind to ZmSWEET4b (the sugar transporter) promoter, either subjecting it to self-activation or enhancing its transcriptional activity [24]. Taken together, WRKY TF functions as either an activator or a repressor in the biotic stressors.

In Hevea, 111 genes have been identified to encode the WRKY proteins. HbWRKY40 is strongly induced by Colletotrichum gloeosporioides and salicylic acid and directly binds to the promoters of MYB and ethylene-responsive transcription factor genes [25], indicating a role in coordinating cross-regulatory transcriptional network during pathogen infection. The overexpression of HbWRKY40 in tobacco triggers ROS burst and enhances resistance to B. cinerea in Arabidopsis [25]. Other HbWRKY TFs, such as HbWRKY1, HbWRKY27, and HbWRKY83, are also involved in wound responses, direct gene regulation, hormone-mediated signaling pathways, drought, salt stress, and natural rubber biosynthesis [26, 27, 28].

2.1.2 AP2/ERF transcription factors

AP2/ERF TF family is classified into four major subfamilies: APETALA2 (AP2), Dehydration-Responsive Element Binding Proteins (DREBs), Related to Abscisic Acid Insensitive 3/ Viviparous 1 (RAV), and Ethylene Responsive Element Binding Factors (ERFs) [29]. The AP2/ERF TFs are involved in the regulation of various downstream target genes related to plant growth, development, response to plant hormones such as abscisic acid (ABA) and ethylene (ET) [30, 31] and mediated responses to both biotic and abiotic stress regulation [31, 32]. However, their downstream regulation does not always confer a beneficial effect in the plant and can also be associated with adverse developmental outcomes [31]. In the context of plant defense, the AP2/ERF TF family is crucial in regulating plant defense mechanisms against biotic stress. Recent studies indicate that AP2/ERF TFs mediate immune signaling via complex regulatory networks, primarily acting downstream of mitogen-activated protein kinase (MAPK) cascades [33]. These TFs coordinate the expression of defense-related genes that are important for systemic acquired resistance (SAR) and localized defense responses against biotrophic and necrotrophic pathogens [33, 34].

Arabidopsis ORA59, an AP2/ERF transcription factor, has been identified as a major regulator of pathogen-induced phytoalexin synthesis that is involved in ET and JA signaling pathways [30, 35] in defense systems against B. cinerea necrotrophic pathogen [30]. Similarly, AtERF014-overexpressing plants exhibited enhanced expression of SA-responsive genes (AtPR1 and AtPR5), and increased ROS accumulation, following infection with Pseudomonas syringae pv. tomato. In contrast, overexpression of AtERF014 suppressed JA-responsive gene AtPDF1.2 during B. cinerea infection [36]. Thus, AP2/ERF transcription factors have dual functions in plant immunity, coordinating antagonistic SA and JA/ET pathway based on the stressors.

The AP2/ERF TFs have also been investigated in transgenic rice for their role in response to biotic stress. Overexpression of OsEREBP1 enhances tolerance to abiotic and biotic stresses, functioning downstream of signaling pathways activated during Xanthomonas oryzae infection. This gene regulates the expression of genes involved in JA and ABA signaling, as well as lipid metabolism [37]. Similarly, OsBIERF3 overexpression enhances immunity against Pyricularia oryzae (syn. Magnaporthe oryzae) and Xanthomona oryzae by activating PRs, MAPK pathways, and cell wall biosynthesis genes, including cellulose synthase genes (OsCes7 and OsCesA9) and glucan endo-1,3-beta-glucosidase genes [38]. In addition, OsERF83 overexpression in plants leads to the upregulation of PR genes, including PR1, PR2, PR3, PR5, and PR10 during P. oryzae infection [39]. In contrast, OsERF922 acts as a negative regulator, its overexpression suppresses defense-related genes, including PRs, PAL, and genes involved in phytoalexin biosynthesis [40]. Collectively, these studies demonstrate that AP2/ERF TFs contribute to pathogen resistance in rice by modulating the expression of key defense-related genes.

Transcriptomic analyses have identified 142 AP2/ERF transcription factor contigs in Hevea, comprising 20 AP2, 115 ERF, 4 RAV, and 3 soloist members [29]. Hevea AP2/ERF transcription factors respond to pathogen attack and are involved in hormone signaling pathways (ET, JA), latex metabolism, and tapping panel dryness [29, 41, 42, 43, 44], although the functional evidence is still limited. Gene expression analyses of HbERF genes following phytohormone treatments revealed that HbERF-IXc4 and HbERF-IXc5 are capable of transactivating the GCC-box, and it is assumed that they may have regulatory roles in modulating defense genes in laticifers [44].

2.1.3 NAC transcription factors

The NAC transcription factor family comprises more than 100 genes in each of Arabidopsis [15], rice [45, 46], maize [47, 48], and wheat [49]. The TFs are classified on basis of a highly conserved DNA-binding domain and by their distinct and sometimes overlapping functional roles in regulating plant development and defense mechanisms [15].

Members of the NAC TF family play important roles in immune response, including plant hormone signaling, hypersensitive response (HR), regulation of defense genes, and interaction with pathogen-derived effectors [14]. AtATAF2, a NAC-domain transcription factor, is induced in response to tobacco mosaic virus (TMV) infection, and its overexpression leads to a reduction in virus accumulation and the upregulation of defense-related marker gene PR1 in infected tissues [50]. In tomato, LeJA2 and LeJA2L genes play distinct roles in the regulation of stomatal movement during pathogen infection. Expression of LeJA2 is induced by the gene LeNCED1, a key gene in the ABA biosynthesis pathway, whereas LeJA2L expression is induced by JA and toxin coronatine (COR) to regulate the transcriptional activation of genes related to metabolism of SA [51].

Genome-wide studies have identified 51 NAC genes in Hevea, of which 11 HbNAC genes are homologous to Arabidopsis NACs [52]. The involvement of HbNAC genes in response to cold and drought stress [52, 53, 54], as well as in infection by C. cassiicola [12] has been documented. Promoter analysis of HbNAC genes by Yang et al. [52] revealed enrichment of cis-regulatory elements associated with environmental stress responses (anaerobic induction, low temperature, drought, defense, and wound responsiveness), hormone responsiveness (ABA, auxin, gibberellin, and SA) and developmental regulation (meristem expression, cell cycle control). Both Arabidopsis and Hevea share conserved NAC subfamilies, although the size of the several subfamilies differs between the two species. Notably, the HbNAC subfamily VIII appears to be specific to Hevea, suggesting species-specific evolutionary divergence that may confer specialized functional role in woody plants.

2.1.4 MYB transcription factor

MYB transcription factors have been identified in eukaryotes, including Arabidopsis, rice, maize, and soybean. The MYB superfamily of TFs is divided into four groups based on the number and position of MYB repeats, namely 1R-MYB, 2R-MYB, 3R-MYB, and 4R-MYB [55]. Among these, the 1R-MYB and 2R-MYB groups represent the largest MYB families in higher plants. The MYB TFs are known to play important roles in growth, development, and stress responses, particularly in signal transduction in response to pathogen invasion [5]. In the context of biotic stress responses, MYB TFs play pivotal roles in plant immunity, functioning as either negative or positive regulators of host defense during fungal pathogenicity [56, 57]. Moreover, the involvement of MYB TFs in hormonal crosstalk, such as SA, JA, and ET pathways, has been investigated as central to the Arabidopsis defense mechanism [5, 58].

In Arabidopsis, numerous MYB TFs, particularly those belonging to the R2R3-MYB subgroup, have been implicated in regulating plant defense responses. For instance, BOTRYTIS SUSCEPTIBLE (BOS1), a R2R3-MYB TF regulate plant responses to biotic stress by mediating signaling from reactive oxygen molecules [59]. Another key member, AtMYB30, is known as one of the HR activators upon pathogen infection [60, 61], modulating the biosynthesis of very-long-chain fatty acids (VLCFAs) and reprogramming toward the other signaling molecules involved in programed cell death [61]. Furthermore, AtMYB15 is involved in reinforcing the plant cell wall through the transcriptional activation of lignin biosynthesis genes upon pathogen challenge [62, 63]. The expression of OsMYB4 in Arabidopsis has been reported to be involved in phenylpropanoid biosynthesis, SAR pathway, and ROS scavenging during P. syringae pv. tomato and B. cinerea infections [64].

In rice, MYB TFs have also been involved in both basal and inducible resistance to various pathogens, including fungal and bacterial infections. For example, OsMYB30, in cooperation with OsMYB55 and OsMYB110, activates the hydroxycinnamic-acid amide (HCAA) in the cinnamate/monolignol pathway, leading to increased HCAA accumulation, such as ferulic acid (FA), which plays important roles in cell-wall reinforcement and exhibits antimicrobial activity against the blast fungus P. oryzae infection [65]. Li et al. [66] revealed that P. oryzae induces the OsMYB30 binding and transactivation of 4-coumarate:coenzyme A ligase (Os4CL3 and Os4CL5) genes, which are responsible for lignin biosynthesis, leading to sclerenchyma thickening.

In Hevea, genome-wide analysis identified over 40 HbMYB proteins from rubber tree laticifer cells, which were grouped into 17 subgroups based on Arabidopsis orthologs. Transient overexpression of Hblmyb19 and Hblmyb44 in tobacco activated the promoter of HbFDPS1 and HbSRPP, which participate in natural rubber biosynthesis [67]. Hbmyb1 expression is reported to be responsive to latex tapping and ethephon stimulation [68]. The effects of wounding and hormonal treatment highlight the pivotal role that MYB factors play in defense regulation. In disease context, the overexpression of Hbmyb1 in transgenic tobacco suppresses stress-induced cell death caused by B. cinerea. In addition, RNA-seq work on wood and lignification in Hevea reported by Meng et al. [69] showed that TFs from the MYB, C2H2, and NAC family are involved in phenylpropanoid biosynthesis pathway, which may be linked to potential stress responses. As a result, the studies on MYB TFs in Hevea during pathogen infection are still limited. The direct role of TFs in the plant-pathogen interaction in disease resistance requires further study. The summary of the involvement of WRKY, AP2/ERF, NAC, and MYB TFs in regulating immune responses to biotic stresses are shown in Table 1.

Table 1: Transcription factors involved in plant responses to biotic stress

Family/
Subfamily
Transcription factor Plant Pathogen Function Reference
WRKY AtWRKY18, AtWRKY40, AtWRKY60 Arabidopsis Pseudomonas Syringae Negative regulators of resistance; modulate SA-dependent defense signaling. [21]
AtWRKY18, AtWRKY40, AtWRKY60 Arabidopsis Botrytis cinerea Positive regulators of resistance; associated with JA-mediated defense pathway. [21]
AtWRKY33 (mutant) Arabidopsis Botrytis cinerea Reduction in JA-mediated regulation of PDF1.2 in atwrky33 mutant. [22, 23]
WRKY TFs Maize Ustilago maydis The WRKY TFs bind to the ZmSWEET4b promoter and self-activate or enhance its transcriptional activity. [24]
AP2/ERF AtORA59 Arabidopsis Botrytis cinerea pathogen-induced phytoalexin synthesis and involvement in ET and JA-signaling pathways. [30, 35]
AtERF15 (overexpression) Arabidopsis Pseudomonas Syringae; B. cinerea Reduction in ROS accumulation in overexpressing plants. [6]
AtERF014 (overexpression) Arabidopsis Botrytis cinerea Positively regulates SA-responsive genes (PR1, PR5) while acting as a suppressor of JA gene (PDF1.2). [36]
OsEREBP1 (overexpression) Rice Xanthomonas oryzae OsEREBP1 overexpression regulates JA and ABA signaling and also involved in lipid metabolism. [38]
OsBIERF3 (overexpression) Rice P. oryzae; Xanthomonas oryzae Enhances plant immunity through the activation of PRs, MAPK pathways, and cell wall biosynthesis genes (OsCes7 and OsCeA9) in overexpressing plants. [38]
OsERF83 (overexpression) Rice P. oryzae Overexpressing plants upregulate PR gene transcripts (PR1, PR2, PR3, PR5, PR10). [39]
OsERF922 (overexpression) Rice P. oryzae OsERF922 overexpression suppresses the activation of defense-related genes (PRs, PAL), and genes involved in phytoalexin biosynthesis. [40]
NAC AtATAF2 (overexpression) Arabidopsis Fusarium oxysporum Enhances resistance to the pathogen. [50]
JA2 and JA2-like Tomato P. Syringae pv. tomato DC3000 JA2 regulates an ABA biosynthetic gene, whereas JA2-like regulates the activation of genes involved in SA metabolism. [51]
MYB BOS1 Arabidopsis B. cinerea;
A. brassicicola
Mediates signaling from reactive oxygen molecules. [59]
AtMYB30 Arabidopsis n.a It acts as an HR activator, an early regulator in MAMP signaling pathway, and a direct transcriptional target of bri1-EMS-suppressor 1 (BES1) which activates brassinosteroid genes. [60, 61]
AtMYB15 Arabidopsis n.a Involved in reinforcing the plant cell wall. [62, 63]
OsMYB30 Rice Pyricularia syringae Acts cooperatively with OsMYB55 and OsMYB110 to increase hydroxycinnamic acid amide (HCCA) accumulation needed for cell-wall reinforcement. [65]
OsMYB30 Rice P. oryzae Transactivates Os4CL and Os4CL5 which are responsible for lignin biosynthesis, thereby leading to sclerenchyma thickening. [66]

n.a – not applicable

2.2 Cis-elements in resistance regulation

Plants constantly encounter a diverse array of biotic stressors and have evolved sophisticated defense systems that rely heavily on rapid and large-scale transcriptional reprogramming [70, 71]. Upon perception of pathogen- or damage-associated molecular patterns, immune signaling cascades involving SA, JA, ET, and other hormones channel signals toward the nucleus, where transcription factors orchestrate the activation or repression of defense-related genes [72, 73]. Central to this regulatory rewiring are cis-regulatory elements (CREs) in the promoter region of genes which govern the transcription of genes. CREs represent the primary interface between upstream immune signaling and downstream transcriptional outputs [74], and thus lie at the core of plant immune specificity, speed, and robustness. Despite the differences in gene promoter architecture of genes, the modulation of defense-related transcription by DNA binding proteins has led to the identification of common defense-associated motifs such as the W-box (TTGACC/T, WRKY TF binding), GCC-box (AGCCGCC, AP2/ERF TF binding), TGA/as‑1 elements (TGACG or TGACGTCA, TGA/bZIP binding) and MYB recognition elements (MREs) [20, 75, 76]. Enrichment patterns across promoters induced by different attackers reflect both shared immune circuitry and pathogen-class–specific signaling, permitting plants to balance general defenses with specialized responses. Below we discuss two common CREs as motif signatures enriched in promoters of genes induced by pathogens and wounding.

2.2.1 W-box

The W-box motif (TTGACC/T) is recognized by the WRKY TFs family via a conserved ~60 amino acids WRKY domain and a C2H2 zinc finger [20, 71, 77, 78]. PTI initiated by PRRs rapidly activates MAPK cascades that in turn induce and activate WRKY TF [79] which bind W-boxes in promoters of early defense genes and amplify transcriptional outputs [71, 80]. Promoter analysis of Arabidopsis early flg22-induced genes reveals the overrepresentation of the W-box cis-acting element [80, 81, 82], and this serves as fast-response binding site for WRKY-mediated transcriptional activation during bacterial perception. Enrichment of W-box and few other cis-acting elements have been recorded in rice during pathogen attack [83]. Mutational analyses of W-boxes in promoter regions of PR-1 [84], NPR1 [85], RLK4 [86], AtWRKY33 [87], RPP8 [88] have demonstrated their relevance to immune regulation in Arabidopsis. Evidence suggests that W-box can direct wounding- and pathogen-inducible gene expression in plants [89, 90]. Resistance (R) genes can be induced during ETI, and the majority of the genes encoding R protein contain a nucleotide-binding site and leucine rich repeats (NLR proteins), and W-box is significantly enriched in their promoter [88]. The W-boxes are important regulatory elements in genes that encode redox components, transcriptional regulators of defense, and antimicrobial genes, thus establishing WRKY TFs role in coordinating defense responsiveness (see [78, 91]). The induction of a pathogenesis-related (PR) gene requires W-box [25, 92, 93], an element which also modulates resistance to pathogen attack [85]. The W-box is suggested to be involved in autoregulation as it is statistically enriched in the promoters of activated AtWRKY genes triggered by SA-dependent defense [91, 94].

2.2.2 GCC-box

The GCC-motif (AGCCGCC) is a CRE primarily identified as an ethylene-responsive element in promoters of PR genes and it is bound by members of the AP2/ERF TF family [76, 90, 95, 96, 97, 98, 99]. Over decades of work, the GCC-box has been shown to play a central role in defense gene regulation in response to necrotrophic pathogens, wounding, and combinations of JA and ET signaling [71, 90, 95, 99, 100, 101]. The motif is found to be over-represented in the promoters of many PR genes [83, 95, 101]. Promoter mutagenesis studies have demonstrated the importance of GCC-box element in methyl jasmonate (MeJA) responsiveness of PDF1.2 [100, 102] and ET regulation of PDF1.2 [103]. Overexpression of ERF96 in Arabidopsis induces several defense genes, among which the genes with GCC-box in their promoters were reported to have the strongest activation and are directly bound by ERF96 [95]. ERF96 is JA and ET responsive, triggers the expression of JA/ET-responsive defense genes enriched in GCC elements in their regulatory region [95]. A tetrameric construct of GCC-box has been shown to be responsive to JA and ET [103], wounding and incompatible biotrophic oomycete Peronospora parasitica pv. Cala2 in transgenic plants [90]. A genome-wide survey in rice identified a large number of promoters harboring the GCC-box, as well as other pathogen-inducible cis-elements, and many of those genes encode disease-resistance/susceptibility proteins and transcription factors [83]. OsERF922, a transcriptional activator, acts on GCC-box–containing genes and influences resistance to the rice blast fungus P. oryzae [40]. The GCC-box defines a core element of JA/ET-driven immune transcription, linking pathogen perception to downstream defense gene activation in plants.

3. Core disease-responsive genes and their functions

3.1 Pathogenesis-related (PR) gene families

Pathogenesis-related (PR) proteins are among the earliest-discovered plant defense proteins, playing a crucial role in the plant defense mechanism against pathogenic invasion. Recently, PR proteins have been grouped into nineteen families (PR-1 to PR-19) based on their characteristics and functions [7, 9]. Depending on the plant-pathogen interactions, various PR proteins are typically involved in SAR defense mechanism, which leads to modifications of the cell wall and phytoalexin production.

Several PR proteins exhibit enzymatic activities that are directly involved in the plant defense mechanism. For instance, β-1,3-glucanases and chitinases degrade structural polysaccharides in pathogen cell walls, peroxidases help generate ROS activity, and ribonucleases may participate in phytohormone-dependent signaling [9]. By contrast, some PR proteins function as non-enzymatic antimicrobials, most notably thionins and defensins, which exhibit toxicity toward invading pathogens [104]. PR-1, PR-2, and PR-5 are SA-responsive and associate with biotrophic pathogens, whereas PR-3, PR-4, and PR-5 are involved in JA and ET signaling pathways and contribute to resistance against necrotrophs [7]. In this chapter, we discuss common PR proteins induced by pathogens in model plants and Hevea.

β -1,3-glucanases, members of the PR-2 protein family, are endonuclease enzymes (E.C.3.2.1.39) found in bacteria, fungi, invertebrates, and plants. β-1,3-glucanases play multiple roles in plant physiological and developmental processes, including cell division and seed germination. These enzymes also act as defenses against pathogens [105, 106, 107]. In the context of plant-pathogen interaction, tobacco class I β-1,3-glucanases inhibit Fusarium solani hyphal germination, resulting in lysis because of hydrolysis of the fungal cell walls [108]. Similarly, in wheat kernels, β-1,3-glucanase plays differential inhibitory effects on the hyphal growth of Fusarium graminearum, spore formation of Penicillium sp., and mycelial morphology of Alternaria spp. [109]. Moreover, the interaction between rice and P. oryzae demonstrated the antifungal function of Gns6, a β-1,3-glucanase, which can inhibit the germination of germ tubes and appressoria of P. oryzae, a pathogen that causes blast disease in rice [110].

Chitinases are hydrolytic enzymes (E.C. 3.2.1.14) that catalyze the hydrolysis of chitin, a polymer of N-acetyl-D-glucosamine. Chitinases in plants function in the degradation of chitin in fungal cell walls. These enzymes hydrolyze pathogen elicitors, transforming them into eliciting oligosaccharides [9]. Researchers have conducted functional studies on PR-3 proteins and their corresponding genes. Filyushin et al. [111] revealed that chitinase genes in garlic (AsCHI2, AsCHI3, and AsCHI7) were induced during Fusarium proliferatum infection and were linked to the SAR pathway, as evidenced by the presence of SA-responsive elements in the promoters of AsCHI2, AsCHI3, and AsCHI7. In addition, He et al. [112] reported that overexpressing strawberry FvChi-14 in Arabidopsis modulated in SA and JA signaling pathways, thereby conferring resistance to C. gloeosporioides. Moreover, chitinase activity from bean expressed in transgenic tobacco delayed disease development following infection with Rhizoctonia solani [113].

In Arabidopsis thaliana, At3g49120 (AtPCb) and At3g49110 (AtPCa) encode class III peroxidases that are involved in H2O2 production in response to Fusarium oxysporum infection [114]. Similarly, the tomato LePrx09 gene is transcriptionally induced by diverse stressors, including wounding, oxidative stress (H2O2 exposure), and infection with Alternaria solani [115]. Collectively, these findings highlight the importance role of peroxidase in plants during pathogen infection, specifically the modulation of H2O2 homeostasis and signaling [8]. Pathogen-responsive genes in plants are summarized in Table 2.

Table 2: Pathogen-responsive genes involved in plant responses to biotic stress

Protein/gene Plant Pathogen Function Reference
β-1,3-glucanase Tobacco Fusarium solani Involved in the inhibition of fungal hyphal germination via hydrolysis of fungal cell walls. [108]
β-1,3-glucanase Wheat Fusarium graminearum Penicillium sp., Alternaria spp. Inhibition of hyphal growth. [109]
β-1,3-glucanase Rice P. oryzae Inhibits the germination of germ tubes and appressoria. [110]
Chitinase
AsCHI2, AsCH13, AsCHI7
Garlic Fusarium proliferatum Involve in SAR through SA-mediated pathway. [111]
Chitinase
FvChi-14
(overexpression)
Arabidopsis Colletotrichum gloeosporioides Linked to SA and JA signaling pathway. [112]
Peroxidase LePrx09
(overexpression)
Tomato Alternaria solani LePrx09 participate in the H2O2 signaling pathway. [115]

In recent years, numerous studies have examined the roles of PR genes in Hevea during pathogen infection. PR-1, PR-2 (β-1,3-glucanase), PR-3 chitinase, and PR-10 genes are induced upon pathogen infection by Microcyclus ulei, Corynespora cassiicola, Phytophthora palmivora, Phytophthora meadii and Rigidoporus microporus [10, 12, 116, 117, 118, 119, 120, 121]. Hurtado Páez et al. [10] reported that M. ulei suppressed ET/JA-associated pathway by downregulating AP2/ERF genes, while inducing SA accumulation and cell-wall reinforcement. This phenomenon is consistent with classical SA-ET/JA crosstalk in plant immunity. During necrotrophic infection by Corynespora cassiicola, both tolerant and susceptible Hevea clones activate defense response, including production of PR protein, and secondary metabolism. In the tolerant clone, transcriptional reprogramming occurs earlier, whereas in the susceptible clone, the transcriptional response is delayed but exhibits a much stronger HR-like burst of gene expression [120]. The summary of pathogen-responsive genes in Hevea is presented in Table 3.

Table 3: Pathogen-responsive genes in Hevea responses to biotic stress

Pathogen/Stress/treatment Analysis/method Finding Gene Reference
Microcyclus ulei
(SALB)
Biotroph
RNA-Seq
(Resistance clone FX 3864)
Induction of PR/chitinase, phenylpropanoid, cell-wall reinforcement, SA-related genes.
Suppression of ET/JA-associated pathway via down-regulation of AP2/ERF genes.
PR1/PR5, chitinase, β-1,3-glucanase, PAL/4CL [10]
Corynespora cassiicola
Cassicolin Cas1 toxin
Necrotroph (toxin-mediated)
RNA-seq
(time course transcriptome analysis of a tolerant and a susceptible clone)
Both clones activate defense genes and secondary metabolisms.
The transcriptional response in the susceptible clone is delayed but exhibits a much stronger HR-like burst of gene expression.
AP2/ERF TFs, chitinase, defensin-like, kinases [120]
Phytophthora palmivora
(hemibiotroph oomycetes)
PTI/priming by elicitor and ABA treatment ABA activated with PRs and cell-wall reinforcement. HbPR1, HbASI, HbCAT [117]
Rigidoporus microporus
(white root disease; necrotroph)
qRT-PCR analysis Differential expression of defense-related genes among clones. Defense-related genes with clones specific
HbPR1, HbPR3, HbPAL (high in tolerance clone)
[119]
Phytophthora meadii RNA-seq; qRT-PCR Early induction of WRKY and MYB TFs and a robust hormonal signaling response during mid-phase of pathogen response in the tolerant rubber clone. HbPR genes; R genes, ROS-related genes [116, 121]

3.2 Reactive oxygen species (ROS) biosynthesis and detoxification genes

ROS are produced naturally from redox reaction during aerobic respiration and photosynthesis. ROS production is required for SAR and its production peaks during PTI have been extensively reviewed in [122]. To mitigate the potentially damaging effect of ROS and to preserve redox homeostasis, plant cells coordinate networks of enzymatic and non-enzymatic scavenging systems. During stress responses, plants synthesize ROS, known as oxidative burst. Microbial elicitors and pathogen attacks trigger a cascade of host gene expressions, some of which play important roles in the biosynthesis and detoxification of ROS. Here, we review the functions and transcription of genes involved in ROS biosynthesis and detoxification as affected by biotic assaults; however discussion of ROS signaling pathways in response to stress is beyond the scope of our review and has been addressed elsewhere (see [122, 123, 124, 125]).

Superoxide (O2-) production is catalyzed by respiratory burst oxidase homologs (Rbohs), a NADPH oxidase (NOX) homolog in plant. Arabidopsis Atrboh transcripts are modulated by environmental cues, where AtrbohC to F are activated by a variety of biotic stresses [126, 127]. AtrbohD and AtrbohF function in ROS production but have distinctive functions in hypersensitive resistance regulation during interactions with avirulent P. syringae and Hyaloperonospora arabidopsis [127]. Enhanced resistance to necrotrophic fungus A. brassicola is conferred in mutants of atrbohD [128]. Repression of Rboh activity altered redox-related metabolism and induced multiple pleiotropic developmental effects, in addition to impairing systemic wound responses [126, 129]. In tomato, silencing of SlRbohB, a homolog of Arabidopsis AtRbohD, reduced resistance to B. cinerea, accompanied by diminished pathogen-induced ROS accumulation and attenuated activation of defense-related genes [130]. Together, these findings highlight Rboh proteins as key enzymatic hubs integrating environmental cues with ROS-mediated defense signaling and developmental regulation.

ROS production is an integral component of phytohormone-mediated signaling pathways. SA is synthesized via the isochorismate and phenylalanine ammonia-lyase pathways [131], processes that are promoted by ROS accumulation [132]. ROS generated by Rboh act in several hormonal signaling pathways [126] and plant immunity responses during pathogen-host interactions [133, 134]. Pogány et al. [128] revealed that SA reinforces ROS signaling by transcriptionally activating ROS-producing genes, thereby promoting localized cell death in infected tissues. Mechanistically, ROS and SA engage in a positive feedback loop in which ROS promotes SA biosynthesis, while SA enhances ROS accumulation through coordinated transcriptional regulation.

Bacterial, fungal, viral, and herbivore infections alter the transcripts of ROS-related genes. A comprehensive computational analysis of transcriptomic datasets by Bilgin et al. [135] revealed coordinated changes in the expression of both ROS biosynthesis and detoxification genes in response to diverse biotic stressors. Notably, among genes encoding ROS-detoxifying enzymes, 30 were upregulated while 35 were downregulated [135], highlighting the dynamic and context-dependent modulation of redox homeostasis during plant-pathogen interactions.

In addition to the genome-wide transcriptomic profiling, specific regulatory mechanisms linking ROS homeostasis to immune signaling have been identified. Zhang et al. [136] reported that the transcriptional repressor MAPK Alfin-like 7 (MAPK-AL7) interacts with nucleotide-binding leucine-rich repeat (NLR) proteins to suppress the expression of ROS-scavenging genes. This suppression promotes NLR-mediated immunity during ETI, leading to enhanced ROS accumulation and increased disease resistance. Transcriptional activation of ROS-scavenging genes has been reported in a study on soybean. Hydrogen peroxide (H2O2) from the oxidative burst triggered by the elicitor Phytophthora sojae induces glutathione S-transferase (GST) and glutathione peroxidase (GPX) transcripts in soybean cells [137].

Genome-wide analyses of redox-related gene families in Hevea have revealed that 30 gene families are involved in ROS production and scavenging, including superoxidase dismutases (SODs), glutathione (GSH)-associated enzymes, class III peroxidases, thioredoxins, and other antioxidant systems. Many of these genes are associated with physiological stress called tapping panel dryness (TPD) [138, 139], and subsets of ROS-scavenging components have also been implicated in pathogen defense. For example, Koop et al. [11] demonstrated that a rubber genotype resistance to Pseudocercospora ulei exhibits strong induction of ROS-related genes, including RBOHs, peroxidases, superoxide dismutase (SOD), catalase (CAT), and ascorbate peroxidase (APX), directly linking resistance to NADPH oxidase-dependent ROS production. These findings suggest that pathogen infection triggers an early oxidative burst mediated by RBOHs and peroxidases, followed by coordinated upregulation of ROS-scavenging enzymes to maintain a balance between defense signaling and oxidative damage. In addition, Wang et al. [140] reported that HbRbohD, a homolog of AtRbohD, is localized to the plasma membrane and is involved in biotic stress responses by modulating SA signaling pathways. Despite limited mechanistic understanding of transcription factor-ROS gene regulation in Hevea, accumulating evidence from RNA-seq and gene expression studies [10, 116, 117, 119, 120, 121] supports the notion that transcriptional control of ROS production and scavenging plays a pivotal role in Hevea immune responses to biotic stress.

4. Comparative overview: Model plants vs Hevea brasiliensis

In classical plant immunity, pathogen perception by PRRs or NLRs activates downstream signaling cascades, including MAPK pathways. These signaling events are mediated by conserved families of transcription factors, such as WRKY, NAC, MYB, and AP2/ERF, which bind to cis-regulatory elements in the promoters of target genes, thereby activating defense-related transcriptional programs. Notably, some TFs function as dual regulators, acting as transcriptional activators or repressors of downstream defense genes depending on pathogen identity, infection stage, and cellular context.

Transcriptional regulation in Hevea shows considerable evolutionary conservation with model plants in core signaling mechanisms, particularly involving the WRKY, AP2/ERF, NAC, and MYB TF families. In Arabidopsis, these TF families are well-characterized and are established as central modulators of immune signaling through specific interactions with cis-regulatory elements, including the W-box, GCC-box, and NAC/MYB-binding elements in promoter regions, thereby promoting immune responses [96, 141].

This conserved transcriptional regulation is evidenced in Hevea by the identification of HbERF-IX genes, one of which is an ortholog of AtERF1 and has been shown to transactivate the GCC-box [44], as well as by the recent characterization of HbRbohD [140]. RbohD is a plasma membrane-localized respiratory burst oxidase that translates regulatory signals into a physical ROS burst. Experimental validation demonstrated that HbRbohD is up-regulated by biotic stimuli and SA, and that HbRbohD overexpression enhances plant resistance to fungal pathogens [140]. Together, these findings suggest that key immune control mechanisms governing ROS burst and defense gene activation are structurally conserved across species, with conserved HbWRKY40 homologs initiating SA-mediated resistance [25].

Despite these conserved elements, distinct functional divergence characterizes the Hevea immune systems, driven in part by the metabolic complexity of its laticifer network as a woody perennial. Unlike Arabidopsis, for which detailed functional information is available on the dual regulatory role of AtWRKY40 [21, 141], Hevea transcription factors exhibit specialized, tissue-specific expression within laticifers, enabling the integration of targeted immune responses that protect latex yield under biotic stress. This notion is further supported by the identification of the HbNAC family VIII, a unique lineage that likely evolved in this particular woody plant physiology [52]. Similarly, the HbMYB TF family has diverged to primarily regulate secondary metabolism, wood formation, and lignification [67, 69], suggesting that Hevea must continuously regulate laticifer tissues and produce wood-forming compounds that are not prioritized in model plants. In addition, although genomic studies in Hevea have identified numerous differentially expressed genes (DEGs) during pathogen infection [12, 116, 120], these studies largely lack functional validation using knockout mutants in model system.

The expansion of these TF families in the rubber genome suggests the presence of specialized regulatory networks that adapt conserved hormonal signaling pathways to the unique physiology of laticifers, while ensuring effective pathogen defense. Finally, by integrating the conserved ability of HbERF-IX proteins to bind the GCC-box with tissue-specific regulatory patterns, Hevea maintains a balance between pathogen defense and the preservation of natural rubber biosynthesis [44, 138]. A schematic representation of hormone and ROS signaling in plant biotic stress responses is shown in Figure 1.

Figure 1: Schematic representation of hormone and ROS signaling in plant biotic stress response. Salicylic acid (SA) is involved in systemic acquired resistance (SAR) in plants through the inhibition of catalase (CAT) and ascorbate peroxidase (APX) activities, thereby enhancing the accumulation of H₂O₂ [142, 143]. CAT and APX function by maintaining ROS homeostasis in cells [144], however, the resulting H₂O₂ induces the expression of defense-related genes associated with SAR [142, 143]. SA-mediated inhibition of antioxidant enzymes elevates ROS levels, thereby activating a positive feedback loop that amplifies SA biosynthesis and the activation of downstream defense genes [145]. This establishes that ROS signaling functions both upstream and downstream of SA signaling. Respiratory burst oxidase homolog (RBOH) proteins are required for ROS production following successful pathogen recognition [125]. SA and ethylene (ET) can promote the spread of cell death [128, 146]. Mutation of RBOHD alters the transcript levels of defense-related genes during necrotrophic fungal pathogen attack. RBOHD does not appear to function as a direct inducer of defense gene expression but instead modulates SA and ET hormone levels to regulate plant cell death [128]. At the transcriptional level, immune signaling converges on conserved transcription factor families, including WRKY, AP2/ERF, NAC, and MYB. WRKY and ERF factors regulate canonical defense genes and hormone-responsive pathways, while NAC and MYB transcription factors integrate stress responses with development, secondary metabolism, and cell wall modification. SA exhibits antagonistic effects on JA/ET-mediated transcriptional responses [147]. Specifically, SA suppresses the expression of a subset of JA/ET-induced genes through repression of ORA59, which encodes a master AP2/ERF transcription factor [148,149].

5. Conclusion and future perspectives

Although gene expression studies in Hevea remain at an early stage compared to model systems, emerging transcriptomic and proteomic evidence already points to the conservation of core plant immune mechanisms. Continued advances in RNA-seq profiling, cis-regulatory analysis, and functional genomics will be essential to deepen our understanding of transcriptional regulation in Hevea and will accelerate the development of disease-resistant rubber cultivars through molecular breeding and targeted genetic improvement.

CRediT authorship contribution statement

Afdholiatus Syafaah: Writing–original draft, Writing–review & editing, Visualization, Conceptualization. Okechukwu S. Ezeh: Writing–original draft, Writing–review & editing, Visualization. Yoshiharu Y. Yamamoto: Writing–review & editing, Conceptualization, Supervision.

Funding

This work was in part supported by the Science and Technology Research Partnership for Sustainable Development (SATREPS).

Conflict of Interest

The authors declare no conflict of interest.

Data availability statement

No new data were generated or analyzed in this study.

References
 
© 2026 The Uniited Graduate Schools of Agricultural Sciences, Japan
feedback
Top