Skip to main content

Comparative analysis of Fusarium crown rot resistance in synthetic hexaploid wheats and their parental genotypes

Abstract

Background

Fusarium crown rot (FCR) is a chronic disease of cereals worldwide. Compared with tetraploid wheat, hexaploid wheat is more resistant to FCR infection. The underlying reasons for the differences are still not clear. In this study, we compared FCR responses of 10 synthetic hexaploid wheats (SHWs) and their tetraploid and diploid parents. We then performed transcriptome analysis to uncover the molecular mechanism of FCR on these SHWs and their parents.

Results

We observed higher levels of FCR resistance in the SHWs compared with their tetraploid parents. The transcriptome analysis suggested that multiple defense pathways responsive to FCR infection were upregulated in the SHWs. Notably, phenylalanine ammonia lyase (PAL) genes, involved in lignin and salicylic acid (SA) biosynthesis, exhibited a higher level of expression to FCR infection in the SHWs. Physiological and biochemical analysis validated that PAL activity and SA and lignin contents of the stem bases were higher in SHWs than in their tetraploid parents.

Conclusion

Overall, these findings imply that improved FCR resistance in SHWs compared with their tetraploid parents is probably related to higher levels of response on PAL-mediated lignin and SA biosynthesis pathways.

Peer Review reports

Background

Hexaploid wheat, or common wheat (Triticum aestivum, AABBDD), is the youngest allohexaploid formed nearly 10,000 years ago as a hybrid between an early-cultivated tetraploid wheat (Triticum turgidum, BBAA) and diploid wild goatgrass (Aegilops tauschii, DD) [1, 2]. Compared with its tetraploid wheat progenitor, hexaploid wheat has greater physiological and ecological plasticity because of its higher ploidy level and complex genome composition [3]. Hexaploid wheat has a broader adaptability to different photoperiod and vernalization requirements than tetraploid wheat. It also has better stress tolerance to biotic and abiotic stresses [3]. Owing to this broader adaptability, hexaploid wheat has spread rapidly and is now one of the most important food crops worldwide [4, 5].

Fusarium crown rot (FCR), primarily caused by Fusarium pseudograminearum, is a soil-borne disease that seriously impacts wheat and barley production. This disease was first reported in Australia and now, it has become a prevalent disease, threating cereal production worldwide [6, 7]. Significant yield losses due to FCR have been documented globally. In Australia, estimated losses of wheat and barley due to FCR are approximately 97 million AUD, with total crop production losses of 10–20% reported [8, 9]. According to data from the Pacific Northwest in the United States, yield losses may be as high as 35% in winter wheat varieties under natural inoculum levels [10]. FCR has recently become one of the most important wheat diseases in the Huanghuai wheat-growing region since 2010 and affects different provinces with varying severity. For example, the highest frequency of FCR occurrence was observed in Henan province, where Xu et al. (2016) [11] reported that the diseased field ratio and white head ratio caused by FCR were 65.1% and 0.1%-31.5% in 2016, respectively. The yield loss in some affected fields reached as high as 51.6%. Although significant progress has been made on genetic resistance [12], currently, no wheat cultivars are completely immune to FCR.

A better understanding of the mechanisms of the host defenses against F. pseudograminearum can provide new strategies to develop varieties with better resistance to FCR. In wheat, several transcriptome analyses have been performed to investigate defense responses in wheat after FCR infection [13,14,15]. For example, comparative transcriptome analysis identified that genes encoding anti-microbial proteins, oxidative stress-related proteins, signaling molecules, and proteins involved in both primary and secondary metabolism were associated with FCR resistance [13]. In a combined transcriptome and metabolite analysis, Powell et al. (2017) found that genes associated with pathogen sensing and signaling, including transcription factor, cellular transport, and detoxification genes, were differentially expressed in the Australian bread wheat cultivar Chara after inoculation with FCR [14]. More recently, integrated transcriptome and metabolite analysis of Chinese wheat cultivars, highlighted that the benzoxazinoid biosynthesis pathway contributes to FCR resistance in wheat [15]. In fact, many responses observed in wheat during F. pseudograminearum infection are very similar to those induced by F. graminearum and F. culmorum [7].

An interesting phenomenon is that hexaploid wheat is generally more resistant to FCR infection than tetraploid wheat. Differences in FCR resistance between tetraploid and hexaploid wheat have been uncovered by several studies. For instance, Daniel et al. (2008) reported that FCR inoculation reduced the yield of durum wheat by 58%, whereas the average yield loss in bread wheat was only 25% across different soil types [16]. Liu et al. (2012) reported that durum wheat has a similar FCR severity as that of barley, whereas bread wheat exhibits less severe FCR stem-base browning symptoms [17]. Ma et al. (2012) assayed the FCR resistance of 2,500 wheat genotypes and found highly resistant genotypes in hexaploid wheat but not in tetraploid wheat [18]. The underlying reasons for the differences between tetraploid and hexaploid wheat are still not clear.

Recent studies have shown that reprogramming of transcriptomes, which occurs in newly formed hexaploid wheat, has distinct effects on growth vigor and growth adaptation following allohexaploidization [4, 19,20,21]. In this study, we compared the FCR responses of 10 synthetic hexaploid wheats (SHWs) (AABBDD) and their parental genotypes, T. turgidum (AABB) and Ae. tauschii (DD). We performed transcriptome analyses of two SHWs and their tetraploid and diploid parents following FCR infection to identify genes and pathways associated with improved FCR resistance of SHWs. We found evidence that defense pathways regulated by phenylalanine ammonia lyase (PAL) genes are altered in SHWs compared with their tetraploid and diploid parents. We also measured six physiological parameters related to FCR resistance. Our findings revealed the genetic basis of improved FCR resistance in SHWs, which may help to develop highly FCR-resistant varieties in future.

Materials and methods

Plant materials

In the study, 10 sets of SHWs and their tetraploid parents (including T. turgidum ssp. dicoccoides, T. turgidum ssp. turgidum and T. turgidum ssp. dicoccon) and diploid parents (including Ae. tauschii ssp. tauschii and Ae. tauschii ssp. strangulata) were used (Table 1). The tetraploid wheat T. turgidum (AABB) and diploid goat grass Ae. tauschii (DD) were used as parental lines to generate SHWs [22]. T. turgidum was used as a maternal parent and pollinated with Ae.tauschii anthers. S1 seeds were obtained from triploid F1 hybrids (ABD) after spontaneous chromosome doubling and were grown as the first generation of allohexaploid plants (AABBDD) [22]. Only plants with complete sets of chromosomes (28 from T. turgidum and 14 from Ae. tauschii) were used for further experimentation. Spontaneous chromosome doubling resulted from the union of un-reduced gametes [23]. No embryo rescue procedure or hormone treatment was applied for the production of the triploid F1 hybrid. All materials were preserved and provided by Triticeae Research Institute of Sichuan Agricultural University.

Table 1 FCR severity assessment of 10 SHWs and their tetraploid and diploid parents

FCR inoculation and assessment

A Fusarium pseudograminearum isolate (NLYY), collected from Henan Province of China, was provided by Henan Institute of Science and Technology. The inoculum preparation, inoculation and FCR assessment were performed as described previously [24]. Briefy, inoculum was prepared on ½ strength potato dextrose agar. Inoculated plates were kept for 12 days at room temperature, and then mycelium was harvested. The mycelia were oscillated in Carboxymethylcellulose sodium Fluid Medium for 5–7 days before being filtered out with sterile gauze and spores were harvested. The concentration of spore suspension was adjusted to 1 × 106 spores/mL and then used directly for inoculation or stored at − 20 °C until needed. Tween 20 was added (0.1% v/v) to the spore suspension prior to use.

Seedlings at 4 days post-germination, were inoculated with Fusarium pseudograminearum isolate (Fp inoculation) or distilled water (mock inoculation) following the protocol described previously [25]. Ten seedlings were used per biological replication with three replications. Briefly, Seeds were treated with 10% available hypochlorite solution for 20 min and then were thoroughly rinsed with distilled water. Surface sterilised seeds were then placed on three layers of filter paper saturated with water and left to germinate. Seedlings (4 days post-germination) were immersed in either spore suspension (Fp inoculation) or in distilled water (mock inoculation) for 1 min and seedlings were planted into the 50-well plastic trays containing sterilised nutrient soil. The plastic trays were arranged in the incubator and the settings for the incubator were: 25/16 (± 1) °C day/night temperature and 65%/85% day/night relative humidity.

FCR severity was assessed on a 0–5 scale at fourth weeks after inoculation, where “0” represents no symptom and “5” represents whole seedling completely necrotic [24]. A disease index (DI) was then calculated for each line following the formula: DI=(∑nX/5 N)×100, where X is the scale value of each plant, n is the number of plants in the category, and N is the total number of plants assessed for each material.

RNA sequencing and expression analysis

Two SHWs (SHW-1 from the cross of T. turgidum ssp. turgidum and Ae. tauschii ssp. tauschii, and SHW-2 from the cross of T. turgidum ssp. dicoccoides and Ae. tauschii ssp. strangulata) and their tetraploid and diploid parents were selected for RNA-seq. The experimental design for each material contained two treatments (mock and Fp inoculation), one time point (3 days post inoculation) and three biological replicates. And a total of ten samples for each biological replicate were used for RNA extraction. Samples were harvested by cutting the stem bases (1 cm) at 3 days post inoculation and snap-frozen in liquid nitrogen and kept at − 80 °C until processed. Total RNA was isolated using a HiPure HP Plant RNA Mini Kit (Magen, R4165–02) according to the manufacturer’s instructions.

A total of 36 libraries were sequenced in the 150 bp paired-end mode using Illumina NovaSeq6000 platform at Berry Genomics Corporation (Beijing, China). Sequence files were deposited at the National Genomics Data Center (NGDC) Genome Sequence Archive under BioProject ID PRJCA011781 (https://ngdc.cncb.ac.cn/search/). Sequence quality assessment, trimming and RNA-seq analysis was performed using CLC Genomics 12.0 software with default parameters. The processed RNA-seq data in SHWs were aligned to complete mRNA (cDNA) sequences of Chinese Spring Reference genome IWGSCv1.1 (a total of 105,200 genes). While RNA-seq data in tetraploid parents and diploid parents were respectively aligned to mRNA sequences of AABB sub-genomes (70,988 genes) and DD sub-genome (34,212 genes) from Chinese Spring Reference genome IWGSCv1.1. The expression level of each gene was calculated using TPM (Transcripts Per Kilobase of exon model per Million mapped reads ). Differentially Expressed Genes (DEGs) were restricted with the absolute value of log2 (fold Change) ≥ 1 and False Discovery Rate (FDR) < 0.05 as the threshold by performing pairwise comparisons of Fp treated and corresponding mock samples. DEGs were separated into up-regulated and down-regulated genes and used as individual test sets for GO enrichment analysis. GO term enrichment analysis was done using Triticeae-GeneTribe (http://wheat.cau.edu.cn/TGT/) [26].

RT-qPCR analysis

Thirteen DEGs were randomly selected and assessed using real-time quantitative PCR (RT-qPCR). Primers were designed using Primer-BLAST tool (http://www.ncbi.nlm.nih.gov/tools/primer-blast/) and listed in Table S1 (Additional file 1). EF-1α, elongation factor of wheat, was used as the internal reference gene. RNA isolation, cDNA synthesis and RT-qPCR analysis were performed as described in Wang et al. (2010) [27]. The relative fold changes were calculated using the comparative CT method (2−∆∆CT). Genes with Ct values > 40 cycles were regarded as having no expression value.

Physiological and biochemical assessment

PAL activity, chitinase activity, salicylic acid (SA) content, jasmonic acid (JA) content and abscisic acid (ABA) content of 5 SHWs (SHW-1, SHW-2, SHW-3, SHW-4 and SHW-5) and their tetraploid and diploid parents were determined using the ELISA method at FANKEL Industrial Corporation (Shanghai, China). The experimental design for each line contained one treatment (Fp inoculation), one time point (3 days post inoculation) and three biological replicates. And a total of ten samples for each biological replicate were used. Samples were harvested by cutting the stem bases and snap-frozen in liquid nitrogen and kept at − 80 °C until processed.

Lignin content of 10 SHWs and their tetraploid parents was measured by ultraviolet spectrophotometry at FANKEL Industrial Corporation (Shanghai, China). To measure basal lignin levels in stem bases, the experimental design for each line contained two time points (seedling and heading stage) and three biological replicates. And a total of ten samples for each biological replicate were used. Samples were harvested by cutting the stem bases and dried to constant weight.

Statistical analysis

Paired samples T-test was performed in Excel 2019 and significant differences was claimed at p < 0.05.

Results

Assessment of FCR severity in SHWs and their parents

In this study, we performed an assessment of the FCR severity of 10 SHWs and their parents at the seedling stage in a greenhouse. Eight SHWs had lower disease index (DI) values than their diploid parents, whereas all SHWs had lower DI values than their tetraploid parents (Fig. 1A–B). Average DI values of SHWs, tetraploid parents, and diploid parents were 41.8 (26.3–50.0), 61.6 (49.0–74.1), and 52.3 (38.5–62.6), respectively (Table 1; Fig. 1B–C). Compared with those of their tetraploid parents, DI values of SHWs were reduced by 16.0–29.0%. Clearly, SHWs had significant higher levels of resistance to FCR (p < 0.01) infection than their parents (Fig. 1C).

Fig. 1
figure 1

– The difference in FCR severity among SHWs and their tetraploid and diploid parents. (A) A representative figure of the FCR responses in 2 SHWs and their tetraploid and diploid parents. (B) DI values of 10 SHWs and their tetraploid and diploid parents. (C) Box plot displays average DI values of 10 SHWs and their tetraploid and diploid parents. **p < 0.01. “SHW” represents synthetic hexaploid wheat, “TW” represents tetraploid wheat, “DW” represents diploid wheat

Transcriptomic response of two SHWs and their parents to FCR infection

To gain insight of the molecular mechanism of the improved FCR resistance in hexaploid wheat, we performed RNA-seq analysis on two SHWs and their corresponding tetraploid and diploid donors. An average of 12.9, 13.1, and 12.3 Gb of sequencing data were generated from the SHWs, tetraploid parents, and diploid parents, respectively (Additional file 1: Table S2). The filtered clean reads were mapped to the Chinese Spring reference genome (IWGSC 1.1). Average mapped ratios of SHWs and their tetraploid and diploid parents were 98.6%, 87.2%, and 87.4%, respectively (Additional file 1: Table S2). Most reads from all samples were successfully mapped, thus indicating that the selected reference genome (IWGSC 1.1) sufficiently reflected the genomic expression of the SHWs and their parents.

To verify the RNA-seq results, we carried out a RT-qPCR analysis using 13 genes (10 up-regulated and 3 down-regulated) randomly selected from the DEGs of two SHWs (Additional file 1: Table S1). Results of RT-PCR analysis were in agreement with the RNA-seq analysis (Additional file 2: Figure S1).

The number of genes differentially expressed in response to FCR were much higher in SHWs than in their tetraploid and diploid parents, and most of these DEGs were up-regulated (Fig. 2A, Additional file 1: Table S3). A total of 431 DEGs (408 up-regulated and 23 down-regulated) were shared between the two SHWs, whereas 122 (90 up-regulated and 32 down-regulated) and 131 (101 up-regulated and 30 down-regulated) were identical between the tetraploid and diploid parents, respectively (Fig. 2B). These shared DEGs may be important signatures of basal resistance of wheat lines with different ploidy levels. In particular, 25 identical DEGs (22 up-regulated and 3 down-regulated) encoding proteins such as dolabradiene monooxygenase, WRKY transcription factors, flavanone 3-dioxygenase, receptor-like protein kinase HSL, and crocetin glucosyltransferase were highly detected in the AB subgenomes of the two SHWs (Fig. 2B–C). 22 identical DEGs from the D subgenome (17 up-regulated and 5 down-regulated) were highly detected in the two SHWs (Fig. 2B-C); these genes encode proteins such as IAA-amino acid hydrolase ILR1-like protein, WRKY transcription factors, sucrose 1-fructosyltransferase, and zinc finger CCCH domain-containing protein.

Fig. 2
figure 2

– Transcriptome differences among 2 SHWs and their tetraploid and diploid parents after FCR infection. (A) The number of up-regulated and down-regulated differentially expressed genes (DEGs). (B) Venn diagrams showing the number of DEGs among 2 SHWs and their tetraploid and diploid parents. Red and blue indicated the number of up-regulated and down-regulated DEGs, respectively. (C) Heat maps showing the identical and specific DEGs in SHWs. Log2 (fold change) in two SHWs ≥ 1 and all parents < 1, or two SHWs ≤ -1 and all parents > -1. (D) Heat maps showing the identical up-regulated GO terms in SHWs. Purple highlighted the differential GO terms among 2 SHWs and their parents. (E) Heat maps showing the identical down-regulated GO terms in SHWs.

Gene ontology (GO) enrichment analysis was performed to classify the functions of DEGs. This analysis revealed that 90, 48, and 26 terms related to biological processes were enriched (FDR < 0.05) in up-regulated DEGs of SHW-1 and its tetraploid and diploid parents, respectively (Additional file 1: Table S4), whereas 24, 36, and 17 functional biological process categories were enriched in down-regulated DEGs (Additional file 1: Table S5). In regard to SHW-2 and its tetraploid and diploid parents, 62, 45, and 20 functional categories under biological process were respectively enriched in up-regulated DEGs (Additional file 1: Table S4), whereas 33, 31, and 14 categories were enriched in down-regulated DEGs (Additional file 1: Table S5). GO terms enriched in up-regulated DEGs of SHWs were more abundant than those of their parents. Several defense-related pathways were commonly enriched in up-regulated DEGs of SHWs and their tetraploid and diploid parents, including cinnamic acid biosynthetic process (GO:0009800), chitin catabolic process (GO:0006032), L-phenylalanine catabolic process (GO:0006559), regulation of defense response (GO:0031347), glutathione metabolic process (GO:0006749), response to oxidative stress (GO:0006979) (Fig. 2D), while pathways related to response to heat (GO:0009408), protein folding (GO:0006457), response to ethanol (GO:0045471), protein complex oligomerization (GO:0051259), and response to salt stress (GO:0009651) were enriched in down-regulated DEGs (Fig. 2E). Some pathways were significantly enriched only in up-regulated DEGs of SHWs, namely, induced systemic resistance (GO:0009864), regulation of response to water deprivation (GO:2,000,070), negative regulation of leaf senescence (GO:1,900,056), and response to salicylic acid (GO:0009751) (Fig. 2D). In addition, the categories of L-phenylalanine catabolic process (GO:0006559) and cinnamic acid biosynthetic process (GO:0009800) were more significantly enriched in DEGs of SHWs than in those of their tetraploid and diploid parents (Fig. 2D).

Differences in expression patterns of defense pathways regulated by PAL genes among SHWs and their tetraploid and diploid parents

PAL, which plays an important role in lignin and salicylic acid biosynthesis, catalyzes the deamination of phenylalanine to form trans-cinnamic acid [28, 29]. L-phenylalanine catabolic process (GO:0006559) and cinnamic acid biosynthetic process (GO:0009800) are synonymous GO terms for functions that are mainly regulated by PAL genes. In this study, 38 PAL genes were found to be up-regulated DEGs in two SHWs (SHW-1 and SHW-2) following FCR infection; among these genes, 19 were significantly up-regulated in both SHWs. Almost all of these PAL genes had a stronger response to FCR infection in the SHWs than in their tetraploid and diploid parents (Fig. 3A–B). At the same time, many genes involved in SA response, such as those encoding WRKY transcription factors, MYB transcription factors, ethylene-responsive transcription factors, wall-associated receptor kinase, and NDR1/HIN1-like protein, had higher expression levels in the two SHWs than in their tetraploid and diploid parents following FCR infection (Fig. 3A). Several genes involved in the lignin biosynthesis pathway were also up-regulated in SHWs and their tetraploid and diploid parents following FCR infection (Fig. 3A). Specifically, part of genes involved in the lignin biosynthesis pathway were likely to be more active in both two SHWs following FCR infection (Fig. 3A, C), including genes encoding cinnamoyl-CoA reductase (TraesCS5A02G061600, TraesCS5A02G061700, TraesCS5B02G069600, TraesCS5B02G160700, TraesCS5D02G073300, TraesCS5D02G168400 and TraesCS7D02G150500), probable 4-coumarate-CoA ligase (TraesCS2A02G145800 and TraesCS2B02G171200), and peroxidase (TraesCS3A02G180000 and TraesCS7A02G353400).

Fig. 3
figure 3

– The difference of expression pattern in defense pathways regulated by Phenylalanine ammonia-lyase genes (PAL) among SHWs and their tetraploid and diploid parents. (A) Expression of genes involved in L-phenylalanine catabolic process (GO:0006559), response to salicylic acid (GO:0009751) and lignin biosynthesis (GO:0009809) in two SHWs and their tetraploid and diploid parents following FCR infection. (B) Scatter plots of fold changes for PAL genes. n represents number of genes. (C) The difference of expression pattern in lignin and salicylic acid biosynthesis pathways between SHWs and their tetraploid parents. C4H, cinnamic acid 4-hydroxylase; C3H, coumarate 3-hydroxylase; 4CL, 4-coumarate coenzyme A ligase; CCoAOMT, caffeoyl CoA O-methyltransferase; CCoA3H: p-coumarate3-hydroxylase; CCR, cinnamoyl-CoA reductase; F5H, ferulate 5-hydroxylase; COMT, caffeic acid O-methyltransferase; CAD, cinnamyl alcohol dehydrogenase; AIM1, abnormal inflorescence meristem 1;BA2H, benzoic acid 2-hydroxylase

Physiological and biochemical parameter associated with resistance to FCR infection in SHWs and their parents

Six physiological parameters related to FCR resistance were also analyzed in the SHWs and their parents. First, we measured PAL and chitinase activities and SA, JA, and ABA contents of the stem base of five SHWs and their parents following FCR infection. As expected, average PAL activities (0.1285 U/g FW) were significantly higher in SHWs than in their tetraploid parents (0.0992 U/g FW) and diploid parents (0.0933 U/g FW) (Fig. 4A–B, Additional file 1: Table S6). Moreover, all five SHWs had higher SA contents than observed in their tetraploid and diploid parents; in particular, average SA contents of the five SHWs, their tetraploid parents, and diploid parents were 73.37, 51.38, and 63.39 ng/g FW, respectively. In contrast to the significantly higher (p < 0.01) SA contents of SHWs compared with their tetraploid parents (Fig. 4A–B, Additional file 1: Table S6), no significant differences in chitinase activity, JA content, and ABA content were found among SHWs and their tetraploid and diploid parents (Fig. 4A-B, Additional file 1: Table S6). We also measured the total lignin content of the stem base at seedling and heading stages. Most of the 10 SHWs had higher lignin contents than those of their tetraploid parents at seedling or heading stages (Fig. 4C, Additional file 1: Table S7), and the overall average lignin content of the 10 SHWs was significantly higher than that of their tetraploid parents both at seedling (p < 0.05) and heading (p < 0.01) stages (Fig. 4D, Additional file 1: Table S7).

Fig. 4
figure 4

– Measurements of six physiological parameters associated with resistance to FCR infection in SHWs and their parents. (A) The bar graphs display PAL activity, chitinase activity, SA content, JA content and ABA content of the stem base in 5 SHWs and their parents following FCR infection, respectively. (B) The box plots display average PAL activity, chitinase activity, SA content, JA content and ABA content of the stem base in 5 SHWs and their parents following FCR infection, respectively. (C) The bar graphs display lignin content of the stem base in 10 SHWs and their tetraploid parents at seedling and heading stages, respectively. (D) The box plots display average lignin content of the stem base in 10 SHWs and their tetraploid parents at seedling and heading stage. *p < 0.05, **p < 0.01. “SHW” represents synthetic hexaploid wheat, “TW” represents tetraploid wheat, “DW” represents diploid wheat

Discussion

FCR has dramatically increased worldwide as a result of drier seasonal conditions and the adoption of moisture-preserving cultural practices, such as minimum tillage and stubble retention [6, 7]. Tetraploid wheat is especially susceptible to FCR, with few highly resistant varieties and germplasms currently available for production and breeding. Similar to salt [4] and N-deficiency [5] tolerances, the FCR resistance of hexaploid wheat is stronger than that of tetraploid wheat [16,17,18]. Nevertheless, little is known about the contributions of polyploidization to FCR resistance in wheat. In this study, we confirmed that improved FCR resistance in the SHWs compared with their tetraploid and diploid parents. This resistance has seen consistently in all 8 out of 10 SHWs (except SHW-6 and SHW-10) derived from different tetraploid and diploid parents (including three T. turgidum and two Ae. tauschii subspecies), thereby exhibiting transgressive performance over both diploid and tetraploid parents.

Transcriptome reprogramming due to the introduction of the D genome occurred in SHWs, leading in turn to global changes in gene expression [19]. These gene expression changes may have led to phenotypic differences between SHW and its parental species and may thus have been important sources of the dominant phenotypes of SHW, such as its growth vigor and adaptation [19, 20]. In this study, we performed transcriptome analyses of two SHWs and their tetraploid and diploid parents following FCR infection. We found that more DEGs and up-regulated GO terms in the SHWs than in their tetraploid and diploid parents following FCR infection. In particular, induced systemic resistance and response to salicylic acid were only significantly enriched in up-regulated DEGs of SHWs. This suggests that SHWs are more likely to have a stronger induced systemic resistance to FCR compared to their parents. Additionally, the L-phenylalanine catabolic process and cinnamic acid biosynthetic process were the most significantly enriched plant defense pathways against FCR infection, and they were commonly enriched in SHWs and their tetraploid and diploid parents. However, these pathways were likely more active in the SHWs than in their parents.

As an entry-point enzyme in the phenylpropanoid biosynthesis pathway, PAL converts phenylalanine to trans-cinnamic acid and free ammonia, thereby playing an important role in lignin and salicylic acid biosynthesis. Many studies have shown that PAL acts as a positive regulator of resistance to various plant biotic stresses [30,31,32,33,34,35,36]. In the transcriptome analyses in the present study, cinnamic acid biosynthesis (GO:0009800) and L-phenylalanine catabolic (GO:0006559) processes were significantly enriched both in SHWs and their parents following FCR infection. A total of 38 PAL genes were detected as up-regulated DEGs in two SHWs following FCR infection, and most of these PAL genes had higher levels of expression in the SHWs than in their tetraploid and diploid parents. Physiological measurements later confirmed that PAL activity was significantly higher in the SHWs than in their tetraploid and diploid parents following FCR infection, which suggests that the SHWs possess a stronger ability to induce PAL activity against FCR infection than their parents.

Lignin, an important component of the plant cell wall, provides an effective barrier for plants to resist pathogen infections and plays a crucial role in plant disease resistance [37]. A recent study has documented that loss of function of the dirigent gene TaDIR-B1 improves resistance to FCR in wheat by increasing the accumulation of lignin [38], thus indicating that lignin plays an important role in FCR resistance. In our study, most up-regulated cinnamoyl-CoA reductase genes had stronger responses in SHWs than in their tetraploid and diploid parents following FCR infection, whereas higher lignin contents were found in the stem bases of the SHWs, both at seedling and heading stages. Alteration of PAL activity in SHWs could lead to changes in lignin contents and thicknesses of cell walls, making SHWs more resistant to FCR infection compared with their tetraploid parents.

The importance of SA in plant defense against pathogen attack has been thoroughly documented [35, 39,40,41,42]. Compared with that in tetraploid parents, higher PAL activity in SHWs is likely to increase levels of SA accumulation and SA-responsive gene expression following FCR infection. As expected, response to SA (GO:0009751) and induced systemic resistance (GO:0009864) were only significantly enriched in up-regulated DEGs of SHWs. Moreover, SA contents of the stem bases of the SHWs were significantly higher than those of their tetraploid parents following FCR infection, while there was no significant difference from their diploid parents. In addition, the SA content of diploid parents was intermediate between that of SHWs and tetraploid parents, consistent with the performance seen in our FCR severity assessment.

Higher levels of response on PAL-mediated lignin and SA biosynthesis pathways may contribute to improved FCR resistance in SHWs compared with their tetraploid parents (Fig. 5). However, the altered defense pathways in SHWs may be diverse and include broader basal resistance. The changes of defense pathways in SHWs may involve two main reasons. One is the resistance effect of D genome genes from diploid parents. Our study has shown that diploid parents were more resistant and had higher SA content than tetraploid parents, indicating that stronger induced systemic resistance in SHWs may be mainly due to enhanced functionality of the D genome. This is consistent with a similar finding in SHWs under salt stress [4]. The other reason is an additive effect of genes from two parents. In fact, most genes showed additive expression in SHWs [43], providing an additive performance for many morphophysiological traits in SHWs. Similarly, SHWs exhibited stronger PAL activity than their tetraploid and diploid parents after FCR infection, which may be due to the presence of more PAL genes and their higher levels of expression in SHWs. Anyway, improved FCR resistance in SHWs which may be accomplished via multiple mechanisms that need further investigation.

Fig. 5
figure 5

 A potential mechanism for improving Fusarium crown rot resistance in hexaploid wheat. Red indicated stronger response pathways in hexaploid wheat than in tetraploid wheat following F. pesudagraminearum infection

Previous transcriptome analyses have revealed various genes and pathways associated with FCR resistance in wheat, such as anti-microbial proteins, oxidative stress-related proteins, signaling molecules, primary and secondary metabolism, pathogen sensing and signaling, cellular transport and detoxification, and benzoxazinoid biosynthesis [13,14,15]. However, these studies mainly focused on comparing different wheat genotypes or cultivars with different FCR resistance levels. Unlike previous findings, our study highlighted the advantage of hexaploid wheat on FCR resistance. The differential defense pathways identified between hexaploid and tetraploid wheat not only help elucidate the genetic basis of FCR resistance in hexaploid wheat but also shed light on the reasons for higher susceptibility in tetraploid wheat. These findings have potential implications for breeding improvement for FCR resistance, especially for highly susceptible tetraploid wheat.

Conclusion

In this study, we compared the FCR responses of 10 SHWs and their parental genotypes, T. turgidum and Ae. tauschii. We confirmed improved FCR resistance in the SHWs compared with their parental genotypes. The transcriptome analysis of two SHWs and their tetraploid and diploid parents suggested that multiple defense pathways responsive to FCR infection were stronger in the SHWs than in the parents. Notably, dozens of PAL genes, involved in lignin and SA biosynthesis, exhibited a higher level of expression to FCR infection in the SHWs. Physiological and biochemical analysis validated that PAL activity and SA and lignin contents of the stem bases were higher in SHWs than their tetraploid parents. In summary, this study suggested that improved FCR resistance in SHWs compared with their tetraploid parents is probably related to higher levels of response on PAL-mediated lignin and SA biosynthesis pathways.

Data Availability

Sequence files were deposited at the National Genomics Data Center (NGDC) Genome Sequence Archive under BioProject ID PRJCA011781.

References

  1. Kihara H. Discovery of the DD-analyser, one of the ancestors of Triticum vulgare. Agric Hortic. 1944;19:889–90.

    Google Scholar 

  2. Mcfadden ES, Sears ER. The artificial synthesis of Triticum spleta. Rec Genet Soc Am. 1944;13:26–7.

    Google Scholar 

  3. Dubcovsky J, Dvorak J. Genome plasticity a key factor in the success of polyploid wheat under domestication. Science. 2007;316:1862–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Yang CW, Zhao L, Zhang HK, Yang ZZ, Wang H, Wen SS, Zhang CY, Rustgi SC, Von WD, Liu B. Evolution of physiological responses to salt stress in hexaploid wheat. Proc Natl Acad Sci U S A. 2014;111:11882–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Yang CW, Yang ZZ, Zhao L, Sun FS, Liu B. A newly formed hexaploid wheat exhibits immediate higher tolerance to nitrogen-deficiency than its parental lines. BMC Plant Biol. 2018;18:1–12.

    Article  CAS  Google Scholar 

  6. Chakraborty S, Liu CJ, Mitter V, Scott JB, Akinsanmi OA, Ali S, Dill MR, Nicol J, Backhouse D, Simpfendorfer S. Pathogen population structure and epidemiology are keys to wheat crown rot and Fusarium head blight management. Australas Plant Pathol. 2006;35:643–55.

    Article  Google Scholar 

  7. Kazan K, Gardiner DM. Fusarium crown rot caused by Fusarium pseudograminearum in cereal crops: recent progress and future prospects. Mol Plant Pathol. 2018;19:1547–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Murray GM, Brennan JP. Estimating disease losses to the australian wheat industry. Australas Plant Pathol. 2009;38:558–70.

    Article  Google Scholar 

  9. Murray GM, Brennan JP. Estimating disease losses to the australian barley industry. Australas Plant Pathol. 2010;39:85–96.

    Article  Google Scholar 

  10. Smiley RW, Gourlie JA, Easley SA, Patterson LM, Whittaker RG. Crop damage estimates for crown rot of wheat and barley in the Pacific Northwest. Plant Dis. 2005;89:595–604.

    Article  PubMed  Google Scholar 

  11. Xu F, Song YL, Zhou YL, Zhang H, Wang JM, Li YH, Han ZH. Occurrence dynamics and characteristics of Fusarium root and crown rot of wheat in Henan Province during 2013–2016. Plant Prot. 2016;42:126–32.

    CAS  Google Scholar 

  12. Liu CJ, Ogbonnaya FC. Resistance to Fusarium crown rot in wheat and barley: a review. Plant Breed. 2015;134:365–72.

    Article  Google Scholar 

  13. Desmond OJ, Manners JM, Schenk PM, Maclean DJ, Kazan K. Gene expression analysis of the wheat response to infection by Fusarium pseudograminearum, Physiol Mol. Plant Pathol. 2008;73:40–7.

    CAS  Google Scholar 

  14. Powell JJ, Carere J, Fitzgerald TL, Stiller J, Covarelli L, Xu Q, Gubler F, Colgrave ML, Gardiner DM, Manners JM, et al. The Fusarium crown rot pathogen Fusarium pseudograminearum triggers a suite of transcriptional and metabolic changes in bread wheat (Triticum aestivum L). Ann Bot. 2017;119:853–67.

    CAS  PubMed  Google Scholar 

  15. Duan SN, Jin JJ, Gao YT, Jin CL, Mu JY, Zhen WC, Sun QX, Xie CJ, Ma J. Integrated transcriptome and metabolite profiling highlights the role of benzoxazinoids in wheat resistance against Fusarium crown rot. Crop J. 2021;10:407–17.

    Article  Google Scholar 

  16. Daniel R, Alliance NG, Simpfendorfer S, Tamworth ND. The impact of crown rot on winter cereal yields. Grains Research and Development Corporation; 2008.

  17. Liu YX, Ma J, Yan W, Yan GJ, Zhou MX, Wei YM, Zheng YL, Liu CJ. Different tolerance in bread wheat, durum wheat and barley to Fusarium crown rot disease caused by Fusarium pseudograminearum. J Phytopathol. 2012;160:412–7.

    Article  Google Scholar 

  18. Ma J, Zhang CY, Liu YX, Yan GJ, Liu CJ. Enhancing Fusarium crown rot resistance of durum wheat by introgressing chromosome segments from hexaploid wheat. Euphytica. 2012;186:67–73.

    Article  Google Scholar 

  19. Li AL, Liu DC, Wu J, Zhao XB, Hao M, Geng SF, Yan J, Jiang XX, Zhang LQ, Wu JY, et al. mRNA and small RNA transcriptomes reveal insights into dynamic homoeolog regulation of allopolyploid heterosis in nascent hexaploid wheat. Plant Cell. 2014;26:1878–900.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Li AL, Geng SF, Zhang LQ, Liu DC, Mao L. Making the bread: insights from newly synthesized allohexaploid wheat. Mol Plant. 2015;8:847–59.

    Article  CAS  PubMed  Google Scholar 

  21. Hao M, Li AL, Shi TW, Luo JT, Zhang LQ, Zhang XC, Ning SZ, Yuan ZW, Zeng DY, Kong XC, et al. The abundance of homoeologue transcripts is disrupted by hybridization and is partially restored by genome doubling in synthetic hexaploid wheat. BMC Genomics. 2017;18:149.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Zhang LQ, Liu DC, Zheng YL, Yan ZH, Dai SF, Li YF, Jiang Q, Ye YQ, Yen Y. Frequent occurrence of unreduced gametes in Triticum turgidum-Aegilops tauschii hybrids. Euphytica. 2010;172:285–94.

    Article  Google Scholar 

  23. Zhang LQ, Yen Y, Zheng YL, Liu DC. Meiotic restriction in emmer wheat is controlled by one or more nuclear genes that continue to function in derived lines. Sex. Plant Reprod. 2007;20:159–66.

    Article  Google Scholar 

  24. Li XM, Liu CJ, Chakraborty S, Manners JM, Kazan K. A simple method for the assessment of crown rot disease severity in wheat seedlings inoculated with Fusarium pseudograminearum. J Phytopathol. 2008;156:751–4.

    Article  Google Scholar 

  25. Chen GD, Habib A, Wei YM, Zheng YL, Shabala S, Zhou MX, Liu CJ. Enhancing Fusarium crown rot resistance by pyramiding large-effect QTL in barley. Mol Breed. 2015;35:1–8.

    Article  Google Scholar 

  26. Chen YM, Song WJ, Xie XM, Wang ZH, Guan PF, Peng HR, Jiao YN, Ni ZF, Sun QX, Guo WL. A collinearity-incorporating homology inference strategy for connecting emerging assemblies in the triticeae tribe as a pilot practice in the plant pangenomic era. Mol Plant. 2020;13:1964–1708.

    Article  Google Scholar 

  27. Wang JR, Wang L, Gulden S, Rocheleau H, Balcerzak M, Hattori J, Cao W, Han F, Zheng YL, Fedak G, et al. RNA profiling of fusarium head blight-resistant wheat addition lines containing the Thinopyrum elongatum chromosome 7E. Can J Plant Pathol. 2010;32:188–214.

    Article  CAS  Google Scholar 

  28. Zhao Q, Dixon RA. Transcriptional networks for lignin biosynthesis: more complex than we thought. Trends Plant Sci. 2011;16:227–33.

    Article  CAS  PubMed  Google Scholar 

  29. Lefevere H, Bauters L, Gheysen G. Salicylic acid biosynthesis in plants. Front Plant Sci. 2020;11:338–44.

    Article  PubMed  PubMed Central  Google Scholar 

  30. He J, Liu YQ, Yuan DY, Duan MJ, Liu YL, Shen ZJ, Yang CY, Qiu ZY, Liu DM, Wen PZ, et al. An R2R3 MYB transcription factor confers brown planthopper resistance by regulating the phenylalanine ammonia-lyase pathway in rice. Proc Natl Acad Sci U S A. 2020;117:271–7.

    Article  CAS  PubMed  Google Scholar 

  31. Zhang CZ, Wang X, Zhang F, Dong LD, Wu JJ, Cheng Q, Qi DY, Yan XF, Jiang LY, Fan SJ, et al. Phenylalanine ammonia-lyase2.1 contributes to the soybean response towards Phytophthora sojae infection. Sci Rep. 2017;7:7242.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Shadle GL, Wesley SV, Korth KL, Chen F, Lamb C, Dixon RA. Phenylpropanoid compounds and disease resistance in transgenic tobacco with altered expression of l-phenylalanine ammonia-lyase. Phytochemistry. 2003;64:153–61.

    Article  CAS  PubMed  Google Scholar 

  33. Chen YP, Li FJ, Tian L, Huang MC, Deng RF, Li XL, Chen W, Wu PZ, Li MR, Jiang HW, et al. The phenylalanine ammonia lyase gene LjPAL1 is involved in plant defense responses to pathogens and plays diverse roles in Lotus japonicus-rhizobium symbioses. Mol Plant Microbe Interact. 2017;30:739–53.

    Article  CAS  PubMed  Google Scholar 

  34. Kim DS, Hwang BK. An important role of the pepper phenylalanine ammonia-lyase gene (PAL1) in salicylic acid-dependent signalling of the defence response to microbial pathogens. J Exp Bot. 2014;65:2295–306.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Zhang HL, Huang QL, Yi L, Song XN, Li L, Deng GB, Liang JJ, Chen F, Yu MQ, Long H. PAL-mediated SA biosynthesis pathway contributes to nematode resistance in wheat. Plant J. 2021;107:698–712.

    Article  CAS  PubMed  Google Scholar 

  36. Zhan C, Li YT, Li H, Wang MR, Gong SJ, Ma DF, Li Y. Phylogenomic analysis of phenylalanine ammonia-lyase (PAL) multigene family and their differential expression analysis in wheat (Triticum aestivum L.) suggested their roles during different stress responses. Front Plant Sci. 2022;13:982457.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Ma QH, Zhu HH, Qiao MY. Contribution of both lignin content and sinapyl monomer to disease resistance in tobacco. Plant Pathol. 2018;67:642–50.

    Article  CAS  Google Scholar 

  38. Yang X, Zhong SB, Zhang QJ, Ren Y, Sun CW, Chen F. A loss-of-function of the dirigent gene TaDIR-B1 improves resistance to Fusarium crown rot in wheat. Plant Biotechnol J. 2021;19:866–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Qi PF, Zhang YZ, Liu CH, Chen Q, Guo ZR, Wang Y, Xu BJ, Jiang YF, Zheng T, Gong X, et al. Functional analysis of FgNahG clarifies the contribution of salicylic acid to wheat (Triticum aestivum) resistance against Fusarium head blight. Toxins. 2019;11:59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Duan L, Liu HB, Li XH, Xiao JH, Wang SP. Multiple phytohormones and phytoalexins are involved in disease resistance to Magnaporthe oryzae invaded from roots in rice. Physiol Plant. 2014;152:486–500.

    Article  CAS  PubMed  Google Scholar 

  41. Pant SR, Irigoyen S, Liu JX, Bedre R, Christensen SA, Schmelz EA, Sedbrook JC, Scholthof KBG, Mandadi KK. Brachypodium phenylalanine ammonia lyase (PAL) promotes antiviral defenses against panicum mosaic virus and its satellites. mBio. 2021;12:e03518–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Yuan W, Jiang T, Du KT, Chen H, Cao YY, Xie JP, Li MF, Carr JP, Wu BM, Fan ZF, et al. Maize phenylalanine ammonia-lyases contribute to resistance to sugarcane mosaic virus infection, most likely through positive regulation of salicylic acid accumulation. Mol Plant Pathol. 2019;20:1365–78.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Zhang HK, Zhu B, Qi B, Gou XW, Dong YZ, Xu CM, Zhang BJ, Huang W, Liu C, Wang XT, Yang CW, Zhou H, et al. Evolution of the BBAA component of bread wheat during its history at the allohexaploid level. Plant Cell. 2014;26:2761–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by the National Natural Science Foundation of China (grant number 32001494) and the Science and Technology Department of Sichuan Province (grant numbers 2022ZDZX0014, 2021YFYZ0002 and 2021YJ0297).

Author information

Authors and Affiliations

Authors

Contributions

YC and YPW carried out the experiment, analyzed the data, and drafted the manuscript; FNG, LL, YQW, HL and MD carried out the phenotypic evaluation; YZZ, ZEP, WL, QTJ, JRW, YMW, JM, QX, HYK, PFQ, ZWY, LQZ, DCL and YLZ provided resources and technique guidance; YFJ and GYC designed and carried out the experiment, formulated the questions, analyzed the data and revised the manuscript. All authors have reviewed and approved the final manuscript.

Corresponding authors

Correspondence to Guoyue Chen or Yunfeng Jiang.

Ethics declarations

Ethics approval and consent to participate

All experimental research and field studies on plants complied with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Additional file 1: Table S1. Primers for real-time quantitative PCR.

Table S2. RNA-seq read information and mapping statistics of 2 SHWs and their tetraploid and diploid parents.

Table S3. Differentially expressed gene (DEGs) of 2 SHWs and their tetraploid and diploid parents.

Table S4. GO terms enriched by up-regulated DEGs of 2 SHWs and their tetraploid and diploid parents.

Table S5. GO terms enriched by down-regulated DEGs of 2 SHWs and their tetraploid and diploid parents.

Table S6. PAL activity, chitinase activity, SA content, JA content and ABA content of the stem base in 5 SHWs and their parents following FCR infection.

Table S7. Total lignin content of the stem base of 10 SHWs and their tetraploid parents at the seedling and heading stage.

Additional file 2: Figure S1. Validation of data of RNA-seq by RT-qPCR

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, Y., Wang, Y., Guan, F. et al. Comparative analysis of Fusarium crown rot resistance in synthetic hexaploid wheats and their parental genotypes. BMC Genomics 24, 178 (2023). https://doi.org/10.1186/s12864-023-09268-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09268-7

Keywords