Castanea sativa is one of the main multipurpose tree species valued for its timber and nuts. This species is susceptible to two major diseases, ink disease and chestnut blight, caused by Phytophthora spp. and Cryphonectria parasitica, respectively. The loss-of-function mutations of genes required for the onset of pathogenesis, referred to as plant susceptibility (S) genes, are one mechanism of plant resistance against pathogens. On the basis of sequence homology, functional domain identification, and phylogenetic analyses, we report for the first time on the identification of S-genes (mlo1, dmr6, dnd1, and pmr4) in the Castanea genus. The expression dynamics of S-genes were assessed in C. sativa and C. crenata plants inoculated with P. cinnamomi and C. parasitica. Our results highlighted the upregulation of pmr4 and dmr6 in response to pathogen infection. Pmr4 was strongly expressed at early infection phases of both pathogens in C. sativa, whereas in C. crenata, no significant upregulation was observed. The infection of P. cinnamomi led to a higher increase in the transcript level of dmr6 in C. sativa compared to C. crenata-infected samples. For a better understanding of plant responses, the transcript levels of defense genes gluB and chi3 were also analyzed.
The
genus belongs to the Fagaceae family and includes four major species of commercial and ecosystemic interest:
Mill. (European chestnut),
Sieb. et Zucc. (Japanese chestnut),
Bl. (Chinese chestnut), and
Borkh (American chestnut).
is a woody species common in all Mediterranean countries and Asia Minor. It has been widely used since ancient times, not only for the consumption of its edible nuts, but also for wood and the products of its ecosystem, such as mushrooms and honey. It is a forest tree, relevant for landscape ecology and biodiversity of mountain and rural areas [1].
Over the last century, the number of chestnut trees decreased in growing areas in Europe due to the depopulation of mountains, climate change, and the spread of two severe diseases: ink disease and chestnut blight [2][3]. Ink disease is caused by the Oomycete
and
. Both species are pathogenic to
, although
generally displays greater virulence than
species, only
exhibits high tolerance to
[6]. The disease, which affects both young and old trees, leads to subcortical necrosis of the root system and the basal part of the stem; this is followed by the appearance of wasting symptoms in the foliage until the total desiccation and death of the plant occur [7][8][9][10]. These pathogens spread mainly through the movement of soil harboring inoculum and the dissemination of asexual flagellated spores (i.e., zoospores) that can actively travel short distances or passively travel long distances in flowing water [10][11]. The use of resistant rootstocks represents one possible solution to protect against these pathogens, although, at present, only tolerant selections obtained from hybridization between
and
are available [12].
Chestnut blight stands among the most destructive fungal tree diseases ever [10][13]. The causal agent,
, infects trees through dead plant tissue and wounds, including those caused by pruning, graft, and hail [13][14]. The symptoms involve bark cankers that can develop on suckers, young branches, and adult branches and trunks [15]. Chestnut blight was one of the causes of the abandonment of chestnut orchards in Europe until the end of the 1970s, when the natural spread of the hypovirulent form of the fungus favored a slow but progressive recovery of chestnut orchards and coppices. However, the fungus still represents a relevant problem in many areas of Europe. It is very harmful to young grafted trees in particular, hampering the establishment of new orchards in many areas [10][13].
forests in Eastern North America were wiped out by
in the early 20th century [16]. Extensive studies and breeding activities have been carried out to restore the American chestnut species introgressing resistance genes of
[17][18]. More recently, researchers discovered that the onset of the disease is associated with the release of oxalic acid by the pathogen during infection. Blight-resistant
trees were obtained by transferring a wheat gene that encodes oxalate oxidase [19].
Recently, a new interest and sensitivity towards the preservation of the local landscape generated a growing interest in silviculture and chestnut trees [20]. Moreover, the market demand for chestnuts in European countries has been strong in the last two decades and has often been supplied by importations. This has been due in part to the gall wasp (
Yasumatsu) infestation, which only recently has been controlled effectively [1][21], and to the general difficulty of developing a modern chestnut industry based on quality cultivars of
that are more tolerant to pathogens. The elucidation of the genetic mechanism behind host–pathogen interaction could thus be useful for the development of novel breeding strategies aimed at achieving resistance or higher tolerance to these pathogens.
Plants take advantage of different defense mechanisms during pathogen attack, and pathogens trigger counter-defense mechanisms. Plants carry pattern recognition receptors (PRRs) able to perceive pathogen-associated molecular patterns (PAMPs); this perception leads to intracellular signal transduction culminating in PAMP-triggered immunity (PTI). PTI is characterized by the production of reactive oxygen species (ROS), the secretion of antimicrobial compounds, and hydrolytic enzymes targeting the pathogen cell wall (chitinase and glucanase) and local cell wall fortifications (through callose deposition) [22].
To suppress PTI, pathogens developed effector molecules able to facilitate pathogen infection by manipulating the host response to support compatibility. Plant resistance (R) genes can detect effectors and trigger effector-triggered immunity (ETI) [23]. The recognition between R genes and effectors causes a cascade of responses involving jasmonic acid (JA) and salicylic acid (SA), culminating in a hypersensitive response (HR) [24].
Most pathogens require the cooperation of the host to establish a compatible interaction. Plant genes supporting compatibility and facilitating infection are called susceptibility (S) genes. S-genes can be divided into three main classes: (a) genes required for the early pathogen infection step (basic compatibility); (b) genes encoding negative regulators of plant immunity; (c) genes necessary for pathogen proliferation (sustained compatibility) [22].
The mutation or loss of an S-gene can thus limit the ability of the pathogen to infect the host and the spread of the disease. The resistance mediated by the S-gene mutation can be pathogen-specific or broad-spectrum. In the former case, the pathway can be implicated in the penetration phase; in the latter, one of the target genes can be involved in constitutive defense responses [22]. Resistance due to the loss of S-genes is generally recessive, differing from the generally dominant resistance mediated by R genes.
Among the S-genes,
(
),
(
),
(
), and
(
) have been characterized in many plant species. The
gene family, encoding seven transmembrane domain proteins, has been characterized in many plant species [25]. Some
homologs act as PM susceptibility factors, as their loss of function results in a distinguished type of resistance known as
resistance. Originally discovered in barley (
L.),
resistance was later shown to occur in several monocot and eudicot species, namely Arabidopsis, tomato (
L.), pea (
L.), pepper (
L.), tobacco (
L.), and wheat (
L.) [26][27] plants. The callose synthase encoded by
is responsible for the production of callose in response to biotic and abiotic stresses. In tomato and potato plants, the knockout and silencing of
led to
and
is involved in the conversion of salicylic acid (SA) to 2,3-dihydroxybenzoic acid (2,3-DHBA) and negatively regulates defense gene expression [30]. Its silencing caused resistance to hemi-biotrophic
,
, and
[31]. Mutants of
, encoding for a cyclic nucleotide-gated cation channel, showed
resistance [29].
At the moment, studies on S-genes in woody plant species have been carried out only for
genes in rubber trees [32], poplar trees [33], apple trees, and grapevines [34]. In our work, we report on the S-genes identification and characterization in
on the basis of sequence homology, functional domain detection and phylogenetic relationships. In addition, the expression dynamics of S-genes were assessed in
and
plants inoculated with the two pathogens,
and
, belonging to different kingdoms. Using the same plant material, the transcription levels of key genes involved in pathogen resistance,
(
) and
(
, were also determined. Our analysis revealed the strong activation of
and
in response to infection by both
and
.
is a European woody tree species commonly used across the globe in the food and timber industries. This chestnut species is susceptible to the two major pathogens,
P. cinnamomiand
C. parasitica [10,36]. In contrast, the Asian chestnut species[10][35]. In contrast, the Asian chestnut species
C. crenataand
C. mollissimaexhibit higher tolerance to
P. cinnamomiand
C. parasitica [6]. Achieving tolerance or resistance to pathogens is the major aim of rootstock breeding. Blight-resistant trees were obtained through backcross breeding of introgression genes from Asian species into American chestnut trees. [37]. However, this approach, although successful in developing blight-resistant American chestnut selections has been slowed by a lack of genetic tools. In Europe, ink disease tolerant hybrids were obtained through interspecific crosses between[6]. Achieving tolerance or resistance to pathogens is the major aim of rootstock breeding. Blight-resistant trees were obtained through backcross breeding of introgression genes from Asian species into American chestnut trees. [36]. However, this approach, although successful in developing blight-resistant American chestnut selections has been slowed by a lack of genetic tools. In Europe, ink disease tolerant hybrids were obtained through interspecific crosses between
C. sativaand
C. crenata, although the nut quality produced by these hybrids is below current market standards [38,39]., although the nut quality produced by these hybrids is below current market standards [37][38].
It has long been recognized that a deep understanding of a pathogen’s biology, host–pathogen interactions, and the resistance mechanisms are fundamental to improving breeding programs. Genomic and transcriptomic analyses have provided the first genetic insights into mechanisms underlying susceptible and resistant chestnut species responses to
P. cinnamomiand
C. parasitica [37,38,40,41,42]. Santos et al. [40] reported the upregulation of a set of candidate genes (e.g.,[36][37][39][40][41]. Santos et al. [39] reported the upregulation of a set of candidate genes (e.g.,
Cast_Gnk2-likeand
Calcium-dependent protein kinase) after
P. cinnamomiinfection, which may trigger HR-like cell death in
C. crenata cells. A significant number of genes involved in the defense against chestnut blight were identified [37].cells. A significant number of genes involved in the defense against chestnut blight were identified [36].
Traditionally, the introduction of resistance gene analogues into plants was the most promising approach to facilitate the acquisition of resistance. However, it did not prove to be durable enough because the widespread use of R genes caused the selection of pathogens capable of overcoming them [24]. Susceptibility (S) genes can be interesting candidates to be used in target breeding programs [22,23,24]. Furthermore, on the basis of previous studies, it was highlighted that the disabling of susceptibility genes may facilitate durable resistance since the pathogen needs to gain a new function to replace the lost host factor it was exploiting [43].Traditionally, the introduction of resistance gene analogues into plants was the most promising approach to facilitate the acquisition of resistance. However, it did not prove to be durable enough because the widespread use of R genes caused the selection of pathogens capable of overcoming them [24]. Susceptibility (S) genes can be interesting candidates to be used in target breeding programs [22][23][24]. Furthermore, on the basis of previous studies, it was highlighted that the disabling of susceptibility genes may facilitate durable resistance since the pathogen needs to gain a new function to replace the lost host factor it was exploiting [42].
In woody species, the investigation of S-genes has been performed only for MLO genes in rubber trees [32], poplar trees [33], apple trees, and grapevines [34,44]. Due to the absence of aIn woody species, the investigation of S-genes has been performed only for MLO genes in rubber trees [32], poplar trees [33], apple trees, and grapevines [34][43]. Due to the absence of a
C. sativagenome, highly similar S-genes were selected using the
C. mollissimav 1.1 genome. Based on the blastn survey, four loci with high similarity were identified in the
C. mollissima genome and attributed to different subclasses of S-genes [31,45,46,47] due to the presence of specific domains:genome and attributed to different subclasses of S-genes [31][44][45][46] due to the presence of specific domains:
mlo1,
dmr6,
dnd1, and
pmr4(
Figure 1,
Table 1). As previously observed [31], in the phylogenetic trees, monocot proteins formed a separate clade with respect to those of dicotyledonous species, supporting the hypothesis that an independent evolution occurred for these genes (
Figure 2). Quantitative PCR analysis has been carried out to identify the differential expression of candidate S-genes in response to
P. cinnamomiand
C. parasiticain the stems of a susceptible species,
C. sativa, and of a tolerant one,
C. crenata. Lesion analysis and DNA quantification of the pathogen (
Figure 4 and3 and
Figure 6) confirmed the higher tolerance level of4) confirmed the higher tolerance level of
C. crenatain response to both
P. cinnamomiand
C. parasiticainfection. Our qPCR results highlighted the main upregulation of
pmr4and
dmr6in response to infection by both
P. cinnamomiand
C. parasitica. As expected, a greater increase in the transcription of these susceptibility genes was observed in the susceptible species
C. sativa. Remarkably,
pmr4was strongly expressed at early infection phases of both pathogens in
C. sativa; in the tolerant
C. crenata, significant upregulation was observed (
Figure 5and
Figure 7).6).
Pmr4codifies for a callose synthase, which is necessary to create a physical barrier to avoid pathogen penetration and is also implicated in plant-triggered immunity suppression.
Pmr4 is thus not only involved in the synthesis of callose, but it also acts as a negative regulator of the salicylic acid pathway [28].is thus not only involved in the synthesis of callose, but it also acts as a negative regulator of the salicylic acid pathway [28].
Figure 1. Chestnut S-genes and their protein structures. The graphical representations of gene exon/intron structures were generated using the http://wormweb.org/exonintron tool (accessed on 31 March 2021) and are shown in the left panel. The exons are indicated with black boxes, whereas introns are shown with lines. In the right panel, the protein structural domains are displayed.
Table 1. S-genes detected in the C. mollissima v1.1 genome and protein domain annotations.
Gene Name | Scaffold | ORF Length (bp) | N° Exons | Size (aa) | Domains | PFAM DOMAINS |
---|---|---|---|---|---|---|
MLO1 | scaffold00101 | 1425 | 13 | 474 | Mlo | PF03094 |
DMR6 | scaffold02358 | 1128 | 4 | 375 | 2OG-FeII_Oxy; DIOX_N | PF03171;PF14226 |
DND1 | scaffold00410 | 1407 | 6 | 468 | cNMP_binding; Ion_trans | PF00027;PF00520 |
PMR4 | scaffold00300 | 5346 | 1 | 1781 | FKS1_dom1; Glucan_synthase | PF14288;PF02364 |
Figure 2. Phylogenetic analysis of the S-genes. The 4 phylogenetic trees of mlo1 (A), dnd1 (B), pmr4 (C), and dmr6 (D) were constructed using MEGAX software by aligning chestnut S-gene coding sequences with NCBI S-gene ortholog coding sequences. The colors indicate the main clades detected, and the arrows underline the location of C. mollissima.
Figure 3. qRT-PCR pathogen DNA quantification after P. cinnamomi inoculation. Data were quantified using the 2−ΔΔCt method based on the Ct values of pathogen genes (ypt and mf1) and actin-7 as a housekeeping gene. Data are the means of three biological replicates ± SE. C. sativa data are normalized with C. sativa 0 hpi control; C. crenata data are normalized with C. crenata 0 hpi control. Different letters associated with the set of means indicate a significant difference based on Tukey’s HSD test (p ≤ 0.05).
Figure 4. qRT-PCR pathogen DNA quantification after C. parasitica inoculation. Data were quantified using the 2−ΔΔCt method based on the Ct values of fungal genes (ypt and mf1) with actin-7 as a housekeeping gene. Data are the means of three biological replicates ± SE. C. sativa data are normalized with the C. sativa 0 hpi control; C. crenata data are normalized with C. crenata 0 hpi control. Different letters associated with the set of means indicate a significant difference based on Tukey’s HSD test (p ≤ 0.05).
Figure 5. qRT-PCR-based transcription profiling after P. cinnamomi inoculation. (A) The S-gene transcription profiles in C. sativa (blue) and C. crenata (green) chestnut species. (B) The expression analysis of genes coding for several pathogenesis-related (PR) proteins in C. sativa (blue) and C. crenata (green) species. In all analyses, Cm7-actin was used as a housekeeping gene. Data are the means of three biological replicates ± SE. C. sativa data are normalized with C. sativa 0 hpi control; C. crenata data are normalized with C. crenata 0 hpi control. Different letters associated with the set of means indicate a significant difference based on Tukey’s HSD test (p ≤ 0.05).
Figure 6. qRT-PCR-based transcription profiling after C. parasitica inoculation. (A) The S-gene transcription profile in C. sativa (blue) and C. crenata (green) chestnut species. (B) The expression analysis of genes coding for several pathogenesis-related (PR) proteins of C. sativa (blue) and C. crenata (green) species. In all the analyses, Cm7-actin was used as the housekeeping gene. The data are the means of three biological replicates ± SE. C. sativa data are normalized with C. sativa 0 hpi control; C. crenata data are normalized with C. crenata 0 hpi control. Different letters associated with the set of means indicate a significant difference based on Tukey’s HSD test (p ≤ 0.05).
Huibers et al. [47] demonstrated that resistance due to the silencing of
Pmr4 is associated with salicylic acid (SA) accumulation rather than with the callose deposition absence. Salicylic acid signaling plays a key role protecting against biotrophic pathogens through the establishment of a hypersensitive response (HR). Saiz-Fernandez et al. [49] revealed the increment of SA levels inis associated with salicylic acid (SA) accumulation rather than with the callose deposition absence. Salicylic acid signaling plays a key role protecting against biotrophic pathogens through the establishment of a hypersensitive response (HR). Saiz-Fernandez et al. [48] revealed the increment of SA levels in
P. cinnamomiinoculated stems, indicating that
P. cinnamomiactivates a defense response similar to that triggered by biotrophic pathogens. Inoculation with both virulent and hypovirulent strains of
C. parasitica led to SA accumulation in European chestnut plantlets that were grown in vitro [50]. Transcriptome analyses carried out in bothled to SA accumulation in European chestnut plantlets that were grown in vitro [49]. Transcriptome analyses carried out in both
C. dentataand
C. mollissima highlighted activation of salicylic-acid-related genes in canker tissues [37].highlighted activation of salicylic-acid-related genes in canker tissues [36].
In chestnut trees, callose deposition around
P. cinnamomi hyphae was detected early in the infection process; however, it does not seem to play a key role in the associated interactions since the pathogen can reach the vascular cylinder in both susceptible and resistant plant genotypes [51]. This result was validated by transcriptomes analyses ofhyphae was detected early in the infection process; however, it does not seem to play a key role in the associated interactions since the pathogen can reach the vascular cylinder in both susceptible and resistant plant genotypes [50]. This result was validated by transcriptomes analyses of
C. sativaand
C. crenata, in which no overexpression of
Callose synthasesafter
P. cinnamomi infection was observed [38].infection was observed [37].
Based on our results and the literature, we can hypothesize that callose accumulation due to the
pmr4upregulation in inoculated
C. sativalines may not be responsible for controlling
P. cinnamomicolonization. We suggest that the upregulation of
pmr4could lead to a negative regulation of the SA pathway that in turn provokes the susceptibility of
C. sativato both
P. cinnamomiand
C. parasitica. A clear link with SA pathway has emerged even with the other chestnut gene candidate
dmr6(downy mildew resistance 6). The mutation of
Arabidopsis dmr6gene, associated with salicylic acid (SA) homeostasis [31], results in the generation of plants that are resistant to bacteria and oomycetes.
Dmr6is involved in the conversion of salicylic acid (SA) to 2,3-dihydroxybenzoic acid (2,3-DHBA) and negatively regulates the expression of defense genes (PR-1, PR-2, and PR-5) [30].
The expression trend of the
Dmr6gene in response to
P. cinnamomiinfection turned out to be similar to the profile of
pmr4. Indeed,
dmr6was strongly expressed at early infection phases of
P. cinnamomiin
C. sativa; in
C. crenatano significant upregulation was detected (
Figure 5). No upregulation of
dmr6in response to
C. parasiticawas highlighted in both plant species (
Figure 7). We can thus hypothesize that6). We can thus hypothesize that
dmr6upregulation observed in
C. sativacould negatively regulate defense gene expression, leading to susceptibility to
P. cinnamomi.
Plants produce a variety of hydrolytic defense enzymes against pathogens, including chitinases, proteases, and also glucanases [52]. The genes coding for several pathogenesis-related (PR) proteins,Plants produce a variety of hydrolytic defense enzymes against pathogens, including chitinases, proteases, and also glucanases [51]. The genes coding for several pathogenesis-related (PR) proteins,
Acidic 26 kDa endochitinasegene (
chi3) and
Glucan endo-1,3-beta-glucosidase Bgene (
gluB), were selected in our analysis because they are considered as responsive to SA-dependent signaling [53,54].), were selected in our analysis because they are considered as responsive to SA-dependent signaling [52][53].
Chi3and
gluB are enzymes that cause the lysis of hyphae of various pathogens, resulting in growth inhibition [55,56,57].are enzymes that cause the lysis of hyphae of various pathogens, resulting in growth inhibition [54][55][56].
In both
C. sativaand
C. crenataplants inoculated with
C. parasitica,
chi3and
gluBwere significantly upregulated. The transcription of
chi3was higher in
C. crenatathan in
C. sativa, presumably as a consequence of the improved defense mechanism against
C. parasitica. Our results are in agreement with Shain et al. [58], who demonstrated the involvement of. Our results are in agreement with Shain et al. [57], who demonstrated the involvement of
b-1,3-glucanaseand
chitinasein chestnut species affected by
C. parasitica. Studies on the role of
chitinasein blight infection mostly involved
C. sativa as a model system [50,59]. In bothas a model system [49][58]. In both
C. dentataand
C. mollissima, transcripts of several compounds expressing
chitinase accumulated more in canker tissues than healthy stem tissues [37]. In order to obtain chestnut plants with potentially increased resistance/tolerance to chestnut blight, the endogenous Ch3gene encoding a chitinase-like protein was over-expressed in the European chestnut through Agrobacterium-mediated transformation [60].accumulated more in canker tissues than healthy stem tissues [36]. In order to obtain chestnut plants with potentially increased resistance/tolerance to chestnut blight, the endogenous Ch3gene encoding a chitinase-like protein was over-expressed in the European chestnut through Agrobacterium-mediated transformation [59].
The emergent CRISPR/Cas9 technology is expected to play a key role in future crop breeding as it allows highly efficient gene editing and generates genetic changes indistinguishable from those arising spontaneously in nature or through conventional breeding [61]. Several examples of edited plants resistant to fungal pathogens have been described [62,63]. For example, genome editing was successfully applied to knock outThe emergent CRISPR/Cas9 technology is expected to play a key role in future crop breeding as it allows highly efficient gene editing and generates genetic changes indistinguishable from those arising spontaneously in nature or through conventional breeding [60]. Several examples of edited plants resistant to fungal pathogens have been described [61][62]. For example, genome editing was successfully applied to knock out
mlo S-genes, leading to Powdery mildew (PM) resistance [44,64,65,66].S-genes, leading to Powdery mildew (PM) resistance [43][63][64][65].
Pmr4and
dmr6 loss-of-function through CRISPR/Cas reduced the susceptibility to PM in tomato plants [28,67]. In our laboratory, we are setting up a CRISPR/Cas9 transformation protocol inloss-of-function through CRISPR/Cas reduced the susceptibility to PM in tomato plants [28][66]. In our laboratory, we are setting up a CRISPR/Cas9 transformation protocol in
Castanea sativa. Our future goal will be to perform the functional characterization using the CRISPR/Cas9 approach of the two candidate genes (
dmr6and
pmr4), while checking if the two genes may also play a role in the interaction between
C. sativaand the emergent nut rot and canker agent
Gnomoniopsis castaneae [68].[67].