Transcriptome analysis of wheat inoculated with Fusarium graminearum

Oligogalacturonide application increases resistance to Fusarium head blight in durum wheat

Fusariosis causes substantial yield losses in the wheat crop worldwide and compromises food safety because of the presence of toxins associated with the fungal disease. Among the current approaches to crop protection, the use of elicitors able to activate natural defense mechanisms in plants is a strategy gaining increasing attention. Several studies indicate that applications of plant cell-wall-derived elicitors, such as oligogalacturonides (OGs) derived from partial degradation of pectin, induce local and systemic resistance against plant pathogens. The aim of this study was to establish the efficacy of OGs in protecting durum wheat (Triticum turgidum subsp. durum), which is characterized by an extreme susceptibility to Fusarium graminearum. To evaluate the functionality of OGs, spikes and seedlings of cv. Svevo were inoculated with OGs, F. graminearum spores, and a co-treatment of both. Results demonstrated that OGs are active elicitors of wheat defenses, triggering typical immune marker genes and determining regulation of fungal genes. Moreover, bioassays on spikes and transcriptomic analyses on seedlings showed that OGs can regulate relevant physiological processes in Svevo with dose-dependent specificity. Thus, the OG sensing system plays an important role in fine tuning immune signaling pathways in durum wheat.

Integrative systems biology of wheat susceptibility to Fusarium graminearum uncovers a conserved gene regulatory network and identifies master regulators targeted by fungal core effectors

BACKGROUND

Plant diseases are driven by an intricate set of defense mechanisms counterbalanced by the expression of host susceptibility factors promoted through the action of pathogen effectors. In spite of their central role in the establishment of the pathology, the primary components of plant susceptibility are still poorly understood and challenging to trace especially in plant-fungal interactions such as in Fusarium head blight (FHB) of bread wheat. Designing a system-level transcriptomics approach, we leveraged the analysis of wheat responses from a susceptible cultivar facing Fusarium graminearum strains of different aggressiveness and examined their constancy in four other wheat cultivars also developing FHB.

RESULTS

In this study, we describe unexpected differential expression of a conserved set of transcription factors and an original subset of master regulators were evidenced using a regulation network approach. The dual-integration with the expression data of pathogen effector genes combined with database mining, demonstrated robust connections with the plant molecular regulators and identified relevant candidate genes involved in plant susceptibility, mostly able to suppress plant defense mechanisms. Furthermore, taking advantage of wheat cultivars of contrasting susceptibility levels, a refined list of 142 conserved susceptibility gene candidates was proposed to be necessary host's determinants for the establishment of a compatible interaction.

CONCLUSIONS

Our findings emphasized major FHB determinants potentially controlling a set of conserved responses associated with susceptibility in bread wheat. They provide new clues for improving FHB control in wheat and also could conceivably leverage further original researches dealing with a broader spectrum of plant pathogens.

Identification and Characterization of a Plant Endophytic Fungus Paraphaosphaeria sp. JRF11 and Its Growth-Promoting Effects

Endophytic fungi establish mutualistic relationships with host plants and can promote the growth and development of plants. In this study, the endophytic fungus JRF11 was isolated from Carya illinoinensis. Sequence analysis of the internal transcribed spacer (ITS) region and 18S rRNA gene combined with colonial and conidial morphology identified JRF11 as a Paraphaosphaeria strain. Plant-fungus interaction assays revealed that JRF11 showed significant growth-promoting effects on plants. In particular, JRF11 significantly increased the root biomass and soluble sugar content of plants. Furthermore, transcriptome analysis demonstrated that JRF11 treatment reprogrammed a variety of genes involved in plant mitogen-activated protein kinase (MAPK) signaling and starch and sucrose metabolism pathways through Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Our research indicates that beneficial endophytic fungi are able to interact with plants and exhibit outstanding plant growth-promoting activities.

Genome-wide identification and characterization of Puccinia striiformis-responsive lncRNAs in Triticum aestivum

Wheat stripe rust (yellow rust) caused by Puccinia striiformis f. sp. tritici (Pst) is a serious biotic stress factor limiting wheat production worldwide. Emerging evidence demonstrates that long non-coding RNAs (lncRNAs) participate in various developmental processes in plants via post-transcription regulation. In this study, RNA sequencing (RNA-seq) was performed on a pair of near-isogenic lines-rust resistance line FLW29 and rust susceptible line PBW343-which differed only in the rust susceptibility trait. A total of 6,807 lncRNA transcripts were identified using bioinformatics analyses, among which 10 lncRNAs were found to be differentially expressed between resistance and susceptible lines. In order to find the target genes of the identified lncRNAs, their interactions with wheat microRNA (miRNAs) were predicted. A total of 199 lncRNAs showed interactions with 65 miRNAs, which further target 757 distinct mRNA transcripts. Moreover, detailed functional annotations of the target genes were used to identify the candidate genes, pathways, domains, families, and transcription factors that may be related to stripe rust resistance response in wheat plants. The NAC domain protein, disease resistance proteins RPP13 and RPM1, At1g58400, monodehydroascorbate reductase, NBS-LRR-like protein, rust resistance kinase Lr10-like, LRR receptor, serine/threonine-protein kinase, and cysteine proteinase are among the identified targets that are crucial for wheat stripe rust resistance. Semiquantitative PCR analysis of some of the differentially expressed lncRNAs revealed variations in expression profiles of two lncRNAs between the Pst-resistant and Pst-susceptible genotypes at least under one condition. Additionally, simple sequence repeats (SSRs) were also identified from wheat lncRNA sequences, which may be very useful for conducting targeted gene mapping studies of stripe rust resistance in wheat. These findings improved our understanding of the molecular mechanism responsible for the stripe rust disease that can be further utilized to develop wheat varieties with durable resistance to this disease.

Identification of Key Gene Network Modules and Hub Genes Associated with Wheat Response to Biotic Stress Using Combined Microarray Meta-analysis and WGCN Analysis

Wheat (Triticum aestivum) is one of the major crops worldwide and a primary source of calories for human food. Biotic stresses such as fungi, bacteria, and diseases limit wheat production. Although plant breeding and genetic engineering for biotic stress resistance have been suggested as promising solutions to handle losses caused by biotic stress factors, a comprehensive understanding of molecular mechanisms and identifying key genes is a critical step to obtaining success. Here, a network-based meta-analysis approach based on two main statistical methods was used to identify key genes and molecular mechanisms of the wheat response to biotic stress. A total of 163 samples (21,792 genes) from 10 datasets were analyzed. Fisher Z test based on the p-value and REM method based on effect size resulted in 533 differentially expressed genes (p < 0.001 and FDR < 0.001). WGCNA analysis using a dynamic tree-cutting algorithm was used to construct a co-expression network and three significant modules were detected. The modules were significantly enriched by 16 BP terms and 4 KEGG pathways (Benjamini-Hochberg FDR < 0.001). A total of nine hub genes (a top 1.5% of genes with the highest degree) were identified from the constructed network. The identification of DE genes, gene-gene co-expressing network, and hub genes may contribute to uncovering the molecular mechanisms of the wheat response to biotic stress.

Biocontrol Effect of Clonostachys rosea on Fusarium graminearum Infection and Mycotoxin Detoxification in Oat (Avena sativa)

Oat (Avena sativa) is susceptible to Fusarium head blight (FHB). The quality of oat grain is threatened by the accumulation of mycotoxins, particularly the trichothecene deoxynivalenol (DON), which also acts as a virulence factor for the main pathogen Fusarium graminearum. The plant can defend itself, e.g., by DON detoxification by UGT-glycosyltransferases (UTGs) and accumulation of PR-proteins, even though these mechanisms do not deliver effective levels of resistance. We studied the ability of the fungal biocontrol agent (BCA) Clonostachys rosea to reduce FHB and mycotoxin accumulation. Greenhouse trials showed that C. rosea-inoculation of oat spikelets at anthesis 3 days prior to F. graminearum inoculation reduced both the amount of Fusarium DNA (79%) and DON level (80%) in mature oat kernels substantially. DON applied to C. rosea-treated spikelets resulted in higher conversion of DON to DON-3-Glc than in mock treated plants. Moreover, there was a significant enhancement of expression of two oat UGT-glycosyltransferase genes in C. rosea-treated oat. In addition, C. rosea treatment activated expression of genes encoding four PR-proteins and a WRKY23-like transcription factor, suggesting that C. rosea may induce resistance in oat. Thus, C. rosea IK726 has strong potential to be used as a BCA against FHB in oat as it inhibits F. graminearum infection effectively, whilst detoxifying DON mycotoxin rapidly.

Knockout of MITOGEN-ACTIVATED PROTEIN KINASE 3 causes barley root resistance against Fusarium graminearum

The roles of mitogen-activated protein kinases (MAPKs) in plant-fungal pathogenic interactions are poorly understood in crops. Here, microscopic, phenotypic, proteomic, and biochemical analyses revealed that roots of independent transcription activator-like effector nuclease (TALEN)-based knockout lines of barley (Hordeum vulgare L.) MAPK 3 (HvMPK3 KO) were resistant against Fusarium graminearum infection. When co-cultured with roots of the HvMPK3 KO lines, F. graminearum hyphae were excluded to the extracellular space, the growth pattern of extracellular hyphae was considerably deregulated, mycelia development was less efficient, and number of appressoria-like structures and their penetration potential were substantially reduced. Intracellular penetration of hyphae was preceded by the massive production of reactive oxygen species (ROS) in attacked cells of the wild-type (WT), but ROS production was mitigated in the HvMPK3 KO lines. Suppression of ROS production in these lines coincided with elevated abundance of catalase (CAT) and ascorbate peroxidase (APX). Moreover, differential proteomic analysis revealed downregulation of several defense-related proteins in WT, and the upregulation of pathogenesis-related protein 1 (PR-1) and cysteine proteases in HvMPK3 KO lines. Proteins involved in suberin formation, such as peroxidases, lipid transfer proteins (LTPs), and the GDSL esterase/lipase (containing "GDSL" aminosequence motif) were differentially regulated in HvMPK3 KO lines after F. graminearum inoculation. Consistent with proteomic analysis, microscopic observations showed enhanced suberin accumulation in roots of HvMPK3 KO lines, most likely contributing to the arrested infection by F. graminearum. These results suggest that TALEN-based knockout of HvMPK3 leads to barley root resistance against Fusarium root rot.

FgMet3 and FgMet14 related to cysteine and methionine biosynthesis regulate vegetative growth, sexual reproduction, pathogenicity, and sensitivity to fungicides in Fusarium graminearum

Fusarium graminearum is a destructive filamentous fungus, which widely exists in wheat and other cereal crops. Cysteine and Methionine are unique sulfur-containing amino acids that play an essential role in protein synthesis and cell life, but their functions and regulation in F. graminearum remain largely unknown. Here we identified two proteins, FgMet3 and FgMet14 in F. graminearum, which are related to the synthesis of cysteine and methionine. We found FgMet3 and FgMet14 were localized to the cytoplasm and there was an interaction between them. FgMet3 or FgMet14 deletion mutants (ΔFgMet3 and ΔFgMet14) were deficient in vegetative growth, pigment formation, sexual development, penetrability and pathogenicity. With exogenous addition of cysteine and methionine, the vegetative growth and penetrability could be completely restored in ΔFgMet3 and ΔFgMet14, while sexual reproduction could be fully restored in ΔFgMet3 and partially restored in ΔFgMet14. ΔFgMet3 and ΔFgMet14 exhibited decreased sensitivity to Congo red stress and increased sensitivity to SDS, NaCl, KCl, Sorbitol, Menadione, and Zn ion stresses. Moreover, FgMet3 and FgMet14 nonspecifically regulate the sensitivity of F. graminearum to fungicides. In conclusion, FgMet3 and FgMet14 interacted to jointly regulate the development, pathogenicity, pigment formation, sensitivity to fungicides and stress factors in F. graminearum.

Exogenous expression of barley HvWRKY6 in wheat improves broad-spectrum resistance to leaf rust, Fusarium crown rot, and sharp eyespot

Systemic acquired resistance (SAR) is a broad-spectrum plant defense phenomena controlled by the salicylic acid receptor NPR1. Key regulators of the SAR signaling pathway showed great potentials to improve crop resistance to various diseases. In our previous investigation, a barley transcription factor gene HvWRKY6 was identified as downstream of NPR1 during SAR. However, the broad-spectrum resistance features and molecular mechanisms of HvWRKY6 remain to be explored. In this study, a transgenic wheat line exogenously expressing HvWRKY6 showed improved resistance to leaf rust, Fusarium crown rot (FCR), and sharp eyespot. The model pathogen Pseudomonas syringae pv. tomato DC3000 was employed to induce the SAR response in wheat plants' leaf region adjacent to the infiltration area. Transcriptome sequencing revealed activation of broad-spectrum defense responses by expressing HvWRKY6 in a pathogen-independent manner. Based on the differentially expressed genes in plant hormone signal transduction, we speculated that the enhanced resistance in HvWRKY6-OE wheat transgenic line was associated with activation of the salicylic acid pathway and suppression of the abscisic acid and jasmonic acid pathways. These findings suggest that the transgenic line HvWRKY6-OE might be applied for the genetic improvement of wheat to several fungal diseases; the underlying resistance mechanism was clarified.

Proteome-Wide Response of Dormant Caryopses of the Weed, Avena fatua, After Colonization by a Seed-Decay Isolate of Fusarium avenaceum

Promoting seed decay is an ecological approach to reducing weed persistence in the soil seedbank. Previous work demonstrated that Fusarium avenaceum F.a.1 decays dormant Avena fatua (wild oat) caryopses and induces several defense enzyme activities in vitro. The objectives of this study were to obtain a global perspective of proteins expressed after F.a.1-caryopsis colonization by conducting proteomic evaluations on (i) leachates, soluble extrinsic (seed-surface) proteins released upon washing caryopses in buffer and (ii) proteins extracted from whole caryopses; interactions with aluminum (Al) were also evaluated in the latter study because soil acidification and associated metal toxicity are growing problems. Of the 119 leachate proteins classified as defense/stress, 80 were induced or repressed. Defense/stress proteins were far more abundant in A. fatua (35%) than in F.a.1 (12%). Avena defense/stress proteins were also the most highly regulated category, with 30% induced and 35% repressed by F.a.1. Antifungal proteins represented 36% of Avena defense proteins and were the most highly regulated, with 36% induced and 37% repressed by F.a.1. These results implicate selective regulation of Avena defense proteins by F.a.1. Fusarium proteins were also highly abundant in the leachates, with 10% related to pathogenicity, 45% of which were associated with host cell wall degradation. In whole caryopsis extracts, fungal colonization generally resulted in induction of a similar set of Avena proteins in the presence and absence of Al. Results advance the hypothesis that seed decay pathogens elicit intricate and dynamic biochemical responses in dormant seeds.

Melatonin delays dark-induced leaf senescence by inducing miR171b expression in tomato

Melatonin functions in multiple aspects of plant growth, development, and stress response. Nonetheless, the mechanism of melatonin in plant carbon metabolism remains largely unknown. In this study, we investigated the influence of melatonin on the degradation of starch in tomato leaves. Results showed that exogenous melatonin attenuated carbon starvation-induced chlorophyll degradation and leaf senescence. In addition, melatonin delayed leaf starch degradation and inhibited the transcription of starch-degrading enzymes after sunset. Interestingly, melatonin-alleviated symptoms of leaf senescence and starch degradation were compromised when the first key gene for starch degradation, α-glucan water dikinase (GWD), was overexpressed. Furthermore, exogenous melatonin significantly upregulated the transcript levels of several microRNAs, including miR171b. Crucially, the GWD gene was identified as a target of miR171b, and the overexpression of miR171b ameliorated the carbon starvation-induced degradation of chlorophyll and starch, and inhibited the expression of the GWD gene. Taken together, these results demonstrate that melatonin promotes plant tolerance against carbon starvation by upregulating the expression of miR171b, which can directly inhibit GWD expression in tomato leaves.

The Fungal Endophyte Penicillium olsonii ML37 Reduces Fusarium Head Blight by Local Induced Resistance in Wheat Spikes

The fungal endophyte Penicillium olsonii ML37 is a biocontrol agent of Fusarium head blight in wheat (caused by Fusarium graminearum), which has shown a limited direct inhibition of fungal growth in vitro. We used RNA-seq and LC-MS/MS analyses to elucidate metabolic interactions of the three-way system Penicillium–wheat–Fusarium in greenhouse experiments. We demonstrated that P. olsonii ML37 colonises wheat spikes and transiently activates plant defence mechanisms, as pretreated spikes show a faster and stronger expression of the defence metabolism during the first 24 h after pathogen inoculation. This effect was transient and the expression of the same genes was lower in the pathogen-infected spikes than in those infected by P. olsonii alone. This response to the endophyte includes the transcriptional activation of several WRKY transcription factors. This early activation is associated with a reduction in FHB symptoms and significantly lower levels of the F. graminearum metabolites 15-acetyl-DON and culmorin. An increase in the Penicillium-associated metabolite asperphanamate confirms colonisation by the endophyte. Our results suggest that the mode of action used by P. olsonii ML37 is via a local defence activation in wheat spikes, and that this fungus has potential as a novel biological alternative in wheat disease control.

Transcriptomics Analysis of Wheat Tassel Response to Tilletia laevis Kühn, Which Causes Common Bunt of Wheat

Tilletia laevis Kühn [synonym T. foetida (Wallr.) Liro] can lead to a wheat common bunt, which is one of the most serious diseases affecting kernels, a serious reduction in grain yield, and losses can reach up to 80% in favorable environments. To understand how wheat tassels respond to T. laevis, based on an RNA-Seq technology, we analyzed a host transcript accumulation on healthy wheat tassels and on tassels infected by the pathogen. Our results showed that 7,767 out of 15,658 genes were upregulated and 7,891 out of 15,658 genes were downregulated in wheat tassels. Subsequent gene ontology (GO) showed that differentially expressed genes (DEGs) are predominantly involved in biological processes, cellular components, and molecular functions. Additionally, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that 20 pathways were expressed significantly during the infection of wheat with T. laevis, while biosynthesis of amino acids, carbon metabolism, and starch and sucrose metabolism pathways were more highly expressed. Our findings also demonstrated that genes involved in defense mechanisms and myeloblastosis (MYB) transcription factor families were mostly upregulated, and the RNA-seq results were validated by quantitative real-time polymerase chain reaction (qRT-PCR). This is the first report on transcriptomics analysis of wheat tassels in response to T. laevis, which will contribute to understanding the interaction of T. laevis and wheat, and may provide higher efficiency control strategies, including developing new methods to increase the resistance of wheat crops to T. laevis-caused wheat common bunt.

The decreased expression of GW2 homologous genes contributed to the increased grain width and thousand‑grain weight in wheat-Dasypyrum villosum 6VS·6DL translocation lines

This study demonstrated that the aberrant transcription of DvGW2 contributed to the increased grain width and thousand-grain weight in wheat-Dasypyrum villosum T6VS·6DL translocation lines. Due to the high immunity to powdery mildew, Dasypyrum villosum 6VS has been one of the most successful applications of the wild relatives in modern wheat breeding. Along with the desired traits, side-effects could be brought when large alien chromosome fragments are introduced into wheat, but little is known about effects of 6VS on agronomic traits. Here, we found that T6VS·6DL translocation had significantly positive effects on grain weight, plant heightand spike length, and small negative effects on total spikelet number and spikelet compactness using recipient and wheat-D. villosum T6VS·6DL allohexaploid wheats, Wan7107 and Pm97033. Further analysis showed that the 6VS segment might exert direct genetic effect on grain width, then driving the increase of thousand-grain weight. Furthermore, comparative transcriptome analysis identified 2549 and 1282 differentially expressed genes (DEGs) and 2220 and 1496 specifically expressed genes (SEGs) at 6 days after pollination (DAP) grains and 15 DAP endosperms, respectively. Enrichment analysis indicated that the process of cell proliferation category was over-represented in the DEGs. Notably, two homologous genes, TaGW2-D1 and DvGW2, were identified as putative candidate genes associated with grain weight and yield. The expression analysis showed that DvGW2 had an aberrant expression in Pm97033, resulting in significantly lower total expression level of GW2 than Wan7107, which drives the increase of grain weight and width in Pm97033. Collectively, our data indicated that the compromised expression of DvGW2 is critical for increased grain width and weight in T6VS·6DL translocation lines.

Transcriptome analysis of sesame-Macrophomina phaseolina interactions revealing the distinct genetic components for early defense responses

Sesame (Sesamum indicum L.) is an oilseed crop challenged by many biotic stresses. Charcoal rot caused by Macrophomina phaseolina (MP) is one of the most devastating diseases of sesame. Till date, molecular mechanisms of resistance to charcoal rot in sesame is not yet reported. In this study, two sesame variety GT-10 (resistant) and RT-373 (susceptible) were identified with contrasting disease incidence when infected with MP. To get the molecular insight, root samples were collected at 0, 24, 48- and 72-h post inoculation (hpi) with the pathogen and generated RNAseq data was analyzed. A total of 1153 and 1226 differentially expressed genes (DEGS) were identified in GT-10 and RT-373, respectively. During the inoculation with MP, resistant genotype showed high number DEGs at early time point of 24 hpi and when compared to late expression in susceptible genotype at 48 hpi. Distinct clusters were represented for each time period represented by cytochrome P450 83B1-like, single anchor, hypothetical protein C4D60, kirola like and heat shock proteins in the resistant genotype contributing for resistance. Analysis of differentially expressed genes, catalogued the genes involved in synthesis of pathogenesis-related (PR) proteins, MYB, WRKY, leucine zipper protein, bHLH, bZIP and NAC transcription factors, ABC transporters (B, C and G subfamily), glutathione metabolism, secondary metabolites, fatty acid biosynthesis and phytohormones like auxin, abscisic acid, ethylene and gibberellic acid. Additionally, in the resistant response we have found three unique GO terms including ATP binding, ribonucleotide binding and nucleic acid binding in molecular function category. The molecular clues generated through this work will provide an important resource of genes contributing for disease resistance and could prioritize genes for functional validation in the important oil crop.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01039-6.

Variation in Gene Expression between Two Sorghum bicolor Lines Differing in Innate Immunity Response

Microbe associated molecular pattern (MAMPs) triggered immunity (MTI) is a key component of the plant innate immunity response to microbial recognition. However, most of our current knowledge of MTI comes from model plants (i.e., Arabidopsis thaliana) with comparatively less work done using crop plants. In this work, we studied the MAMP triggered oxidative burst (ROS) and the transcriptional response in two Sorghum bicolor genotypes, BTx623 and SC155-14E. SC155-14E is a line that shows high anthracnose resistance and the line BTx623 is susceptible to anthracnose. Our results revealed a clear variation in gene expression and ROS in response to either flagellin (flg22) or chitin elicitation between the two lines. While the transcriptional response to each MAMP and in each line was unique there was a considerable degree of overlap, and we were able to define a core set of genes associated with the sorghum MAMP transcriptional response. The GO term and KEGG pathway enrichment analysis discovered more immunity and pathogen resistance related DEGs in MAMP treated SC155-14E samples than in BTx623 with the same treatment. The results provide a baseline for future studies to investigate innate immunity pathways in sorghum, including efforts to enhance disease resistance.

Comparative transcriptome analysis reveals distinct gene expression profiles in Brachypodium distachyon infected by two fungal pathogens

BACKGROUND

The production of cereal crops is frequently affected by diseases caused by Fusarium graminearum and Magnaporthe oryzae, two devastating fungal pathogens. To improve crop resistance, many studies have focused on understanding the mechanisms of host defense against these two fungi individually. However, our knowledge of the common and different host defenses against these pathogens is very limited.

RESULTS

In this study, we employed Brachypodium distachyon as a model for cereal crops and performed comparative transcriptomics to study the dynamics of host gene expression at different infection stages. We found that infection with either F. graminearum or M. oryzae triggered massive transcriptomic reprogramming in the diseased tissues. Numerous defense-related genes were induced with dynamic changes during the time course of infection, including genes that function in pattern detection, MAPK cascade, phytohormone signaling, transcription, protein degradation, and secondary metabolism. In particular, the expression of jasmonic acid signaling genes and proteasome component genes were likely specifically inhibited or manipulated upon infection by F. graminearum.

CONCLUSIONS

Our analysis showed that, although the affected host pathways are similar, their expression programs and regulations are distinct during infection by F. graminearum and M. oryzae. The results provide valuable insight into the interactions between B. distachyon and two important cereal pathogens.

Fusarium head blight resistance in European winter wheat: insights from genome-wide transcriptome analysis

Background

Fusarium head blight (FHB) is a devastating disease of wheat worldwide. Resistance to FHB is quantitatively controlled by the combined effects of many small to medium effect QTL. Flowering traits, especially the extent of extruded anthers, are strongly associated with FHB resistance.

Results

To characterize the genetic basis of FHB resistance, we generated and analyzed phenotypic and gene expression data on the response to Fusarium graminearum (Fg) infection in 96 European winter wheat genotypes, including several lines containing introgressions from the highly resistant Asian cultivar Sumai3. The 96 lines represented a broad range in FHB resistance and were assigned to sub-groups based on their phenotypic FHB severity score. Comparative analyses were conducted to connect sub-group-specific expression profiles in response to Fg infection with FHB resistance level. Collectively, over 12,300 wheat genes were Fusarium responsive. The core set of genes induced in response to Fg was common across different resistance groups, indicating that the activation of basal defense response mechanisms was largely independent of the resistance level of the wheat line. Fg-induced genes tended to have higher expression levels in more susceptible genotypes. Compared to the more susceptible non-Sumai3 lines, the Sumai3-derivatives demonstrated higher constitutive expression of genes associated with cell wall and plant-type secondary cell wall biogenesis and higher constitutive and Fg-induced expression of genes involved in terpene metabolism. Gene expression analysis of the FHB QTL Qfhs.ifa-5A identified a constitutively expressed gene encoding a stress response NST1-like protein (TraesCS5A01G211300LC) as a candidate gene for FHB resistance. NST1 genes are key regulators of secondary cell wall biosynthesis in anther endothecium cells. Whether the stress response NST1-like gene affects anther extrusion, thereby affecting FHB resistance, needs further investigation.

Conclusion

Induced and preexisting cell wall components and terpene metabolites contribute to resistance and limit fungal colonization early on. In contrast, excessive gene expression directs plant defense response towards programmed cell death which favors necrotrophic growth of the Fg pathogen and could thus lead to increased fungal colonization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07800-1.

Transcriptome-based analysis of resistance mechanism to black point caused by Bipolaris sorokiniana in wheat

Black point is a cereal disease caused by complex pathogens, of which the pathogenicity of Bipolaris sorokiniana is the most serious in wheat. Resistance to black point is quantitative in nature, and thus the mechanism is poorly understood. We conducted a comparative transcriptome analysis to identify differentially expressed genes (DEGs) in black point-slightly susceptible and -highly susceptible wheat lines at different timepoints following B. sorokiniana inoculation. DEGs associated with photosynthesis were upregulated in black point-slightly susceptible lines. The top Gene Ontology enrichment terms for biological processes were oxidation–reduction, response to cold, salt stress, oxidative stress, and cadmium ion; terms for cellular component genes were mainly involved in plasma membrane and cytoplasmic membrane-bounded vesicle, whereas those for molecular function were heme binding and peroxidase activity. Moreover, activities of antioxidant enzymes superoxide dismutase, catalase, and peroxidase were higher in slightly susceptible lines than those in highly susceptible lines (except peroxidase 12–24 days post-inoculation). Thus, resistance to B. sorokiniana-caused black point in wheat was mainly related to counteracting oxidative stress, although the specific metabolic pathways require further study. This study presents new insights for understanding resistance mechanisms of selected wheat lines to black point.

Proteomics-Based Data Integration of Wheat Cultivars Facing Fusarium graminearum Strains Revealed a Core-Responsive Pattern Controlling Fusarium Head Blight

Fusarium head blight (FHB), mainly occurring upon Fusarium graminearum infection in a wide variety of small-grain cereals, is supposed to be controlled by a range of processes diverted by the fungal pathogen, the so-called susceptibility factors. As a mean to provide relevant information about the molecular events involved in FHB susceptibility in bread wheat, we studied an extensive proteome of more than 7,900 identified wheat proteins in three cultivars of contrasting susceptibilities during their interaction with three F. graminearum strains of different aggressiveness. No cultivar-specific proteins discriminated the three wheat genotypes, demonstrating the establishment of a core proteome regardless of unequivocal FHB susceptibility differences. Quantitative protein analysis revealed that most of the FHB-induced molecular adjustments were shared by wheat cultivars and occurred independently of the F. graminearum strain aggressiveness. Although subtle abundance changes evidenced genotype-dependent responses to FHB, cultivar distinction was found to be mainly due to basal abundance differences, especially regarding the chloroplast functions. Integrating these data with previous proteome mapping of the three F. graminearum strains facing the three same wheat cultivars, we demonstrated strong correlations between the wheat protein abundance changes and the adjustments of fungal proteins supposed to interfere with host molecular functions. Together, these results provide a resourceful dataset that expands our understanding of the specific molecular events taking place during the wheat-F. graminearum interaction.

Screening and functional characterization of candidate resistance genes to powdery mildew from Dasypyrum villosum#4 in a wheat line Pm97033

KEY MESSAGE

Three genes designated DvLox, Pm21#4, and Pm21#4-H identified in a wheat-Dasypyrum villosum#4 T6V#4S·6DL translocation line Pm97033 conferred wheat for powdery mildew resistance. Powdery mildew (PM) caused by Blumeria graminis f. sp. tritici (Bgt) is one of the most devastating diseases in wheat. To date, only a few genes conferring resistance to wheat PM are cloned. Dasypyrum villosum is a wild relative of wheat, which provides Pm21 conferring wheat immunity to PM. In this study, we obtained many differentially expressed genes (DEGs) from a wheat-D. villosum#4 T6V#4S·6DL translocation line Pm97033 using RNA-sequencing. Among them, 7 DEGs associated with pathogen resistance were up-regulated in front of Bgt infection. Virus-induced gene silencing and transformation assays demonstrated that two of them, DvLox and Pm21#4 encoding a lipoxygenase and a encoding coiled-coil/nucleotide-binding site/leucine-rich repeat resistance protein, conferred wheat PM resistance. The transgenic wheat plants expressing DvLox enhanced PM resistance, and the transgenic wheat plants expressing Pm21#4 showed PM immunity. The Pm21#4-silenced Pm97033 plants by the cluster regularly interspaced short palindromic repeats-associated endonuclease (CRISPR/Cas9) system were susceptible to PM. Thus, Pm21#4 is a key gene contributing PM immune resistance in Pm97033. Constitutively expression of Pm21#4-H, which is silenced in Pm97033 and D. villosum#4, endowed a PM-susceptible wheat variety Fielder with PM immune resistance.

Comparative transcriptomic analysis reveals that multiple hormone signal transduction and carbohydrate metabolic pathways are affected by Bacillus cereus in Nicotiana tabacum

Bacillus cereus is thought to be a beneficial bacterium for plants in several aspects, such as promoting plant growth and inducing plant disease resistance. However, there is no detailed report on the effect of Bacillus cereus acting on Nicotiana tabacum. In the present study, RNA-based sequencing (RNA-seq) was used to identify the molecular mechanisms of the interaction between B. cereus CGMCC 5977 and N. tabacum. A total of 7345 and 5604 differentially expressed genes (DEGs) were identified from leaves inoculated with Bacillus cereus at 6 and 24 hpi, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses revealed that the most DEGs could be significantly enriched in hormone signal transduction, the MAPK signaling pathway, photosynthesis, oxidative stress, and amino sugar, and nucleotide sugar metabolism. Furthermore, glycolysis/gluconeogenesis was severely affected by inoculation with Bacillus cereus. In the hormone signal pathway, multiple DEGs were involved in plant defense-related major hormones, including activation of jasmonic acid (JA), salicylic acid (SA), and ethylene (Eth). Further analyses showed that other hormone-related genes involved in abscisic acid (ABA), gibberellin (GA), auxin (AUX), and cytokinin (CK) also showed changes. Notably, a large number of genes associated with glycolysis/gluconeogenesis, catabolism of starch and oxidative stress were induced. In addition, the majority of DEGs related to nucleic acid sugar metabolism were also significantly upregulated. Biochemical assays showed that the starch content of B. cereus-treated leaves was reduced to 2.51 mg/g and 2.38 mg/g at 6 and 24 hpi, respectively, while that of the control sample was 5.42 mg/g. Overall, our results demonstrated that multiple hormone signal transduction and carbohydrate metabolic pathways are involved in the interaction of tobacco and B. cereus.

Epigenetic regulation of gene expression improves Fusarium head blight resistance in durum wheat

Eight advanced durum-breeding lines were treated with 5-methyl-azacytidine to test the feasibility of generating sources of Fusarium head blight (FHB) resistance. Of the 800 treated seeds, 415 germinated and were advanced up to four (M4) generations by selfing. Thirty-two of the resulting 415 M4 lines were selected following preliminary screening and were further tested for FHB resistance for three years at two field locations, and in the greenhouse. Five of the 32 M4 lines showed less than 30% disease severity, as compared to the parental lines and susceptible checks. Fusarium-damaged kernels and deoxynivalenol analyses supported the findings of the field and greenhouse disease assessments. Two of the most resistant M4 lines were crossed to a susceptible parent, advanced to third generation (BC₁:F₃) and were tested for stability and inheritance of the resistance. About, one third of the BC₁:F₃ lines showed FHB resistance similar to their M4 parents. The overall methylation levels (%) were compared using FASTmC method, which did not show a significant difference between M4 and parental lines. However, transcriptome analysis of one M4 line revealed significant number of differentially expressed genes related to biosynthesis of secondary metabolites, MAPK signaling, photosynthesis, starch and sucrose metabolism, plant hormone signal transduction and plant-pathogen interaction pathways, which may have helped in improved FHB resistance.

A Comparative Transcriptome Analysis, Conserved Regulatory Elements and Associated Transcription Factors Related to Accumulation of Fusariotoxins in Grain of Rye (Secale cereale L.) Hybrids

Detoxification of fusariotoxin is a type V Fusarium head blight (FHB) resistance and is considered a component of type II resistance, which is related to the spread of infection within spikes. Understanding this type of resistance is vital for FHB resistance, but to date, nothing is known about candidate genes that confer this resistance in rye due to scarce genomic resources. In this study, we generated a transcriptomic resource. The molecular response was mined through a comprehensive transcriptomic analysis of two rye hybrids differing in the build-up of fusariotoxin contents in grain upon pathogen infection. Gene mining identified candidate genes and pathways contributing to the detoxification of fusariotoxins in rye. Moreover, we found cis regulatory elements in the promoters of identified genes and linked them to transcription factors. In the fusariotoxin analysis, we found that grain from the Nordic seed rye hybrid "Helltop" accumulated 4 times higher concentrations of deoxynivalenol (DON), 9 times higher nivalenol (NIV), and 28 times higher of zearalenone (ZEN) than that of the hybrid "DH372" after artificial inoculation under field conditions. In the transcriptome analysis, we identified 6675 and 5151 differentially expressed genes (DEGs) in DH372 and Helltop, respectively, compared to non-inoculated control plants. A Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that DEGs were associated with glycolysis and the mechanistic target of rapamycin (mTOR) signaling pathway in Helltop, whereas carbon fixation in photosynthesis organisms were represented in DH372. The gene ontology (GO) enrichment and gene set enrichment analysis (GSEA) of DEGs lead to identification of the metabolic and biosynthetic processes of peptides and amides in DH372, whereas photosynthesis, negative regulation of catalytic activity, and protein-chromophore linkage were the significant pathways in Helltop. In the process of gene mining, we found four genes that were known to be involved in FHB resistance in wheat and that were differentially expressed after infection only in DH372 but not in Helltop. Based on our results, we assume that DH372 employed a specific response to pathogen infection that led to detoxification of fusariotoxin and prevented their accumulation in grain. Our results indicate that DH372 might resist the accumulation of fusariotoxin through activation of the glycolysis and drug metabolism via cytochrome P450. The identified genes in DH372 might be regulated by the WRKY family transcription factors as associated cis regulatory elements found in the in silico analysis. The results of this study will help rye breeders to develop strategies against type V FHB.

Germplasms, genetics and genomics for better control of disastrous wheat Fusarium head blight

Fusarium head blight (FHB), or scab, for its devastating nature to wheat production and food security, has stimulated worldwide attention. Multidisciplinary efforts have been made to fight against FHB for a long time, but the great progress has been achieved only in the genomics era of the past 20 years, particularly in the areas of resistance gene/QTL discovery, resistance mechanism elucidation and molecular breeding for better resistance. This review includes the following nine main sections, (1) FHB incidence, epidemic and impact, (2) causal Fusarium species, distribution and virulence, (3) types of host resistance to FHB, (4) germplasm exploitation for FHB resistance, (5) genetic control of FHB resistance, (6) fine mapping of Fhb1, Fhb2, Fhb4 and Fhb5, (7) cloning of Fhb1, (8) omics-based gene discovery and resistance mechanism study and (9) breeding for better FHB resistance. The advancements that have been made are outstanding and exciting; however, judged by the complicated nature of resistance to hemi-biotrophic pathogens like Fusarium species and lack of immune germplasm, it is still a long way to go to overcome FHB.

Dissecting the Genetic Complexity of Fusarium Crown Rot Resistance in Wheat

Fusarium crown rot (FCR) is one of the most important diseases of wheat (Triticum aestivum L.). FCR is mainly caused by the fungal pathogens Fusarium culmorum and F. pseudograminearum. In order to identify new sources of resistance to FCR and to dissect the complexity of FCR resistance, a panel of 161 wheat accessions was phenotyped under growth room (GR) and greenhouse conditions (GH). Analysis of variance showed significant differences in crown rot development among wheat accessions and high heritability of genotype-environment interactions for GR (0.96) and GH (0.91). Mixed linear model analysis revealed seven novel quantitative trait loci (QTLs) linked to F. culmorum on chromosomes 2AL, 3AS, 4BS, 5BS, 5DS, 5DL and 6DS for GR and eight QTLs on chromosomes on 3AS, 3BS, 3DL, 4BS (2), 5BS, 6BS and 6BL for GH. Total phenotypic variances (R²) explained by the QTLs linked to GR and GH were 48% and 59%, respectively. In addition, five favorable epistasis interactions among the QTLs were detected for both GR and GH with and without main effects. Epistatic interaction contributed additional variation up to 21% under GR and 7% under GH indicating strong effects of environment on the expression of QTLs. Our results revealed FCR resistance responses in wheat to be complex and controlled by multiple QTLs.

Searching for FHB Resistances in Bread Wheat: Susceptibility at the Crossroad

Fusarium head blight (FHB), primarily caused by Fusarium graminearum, is one of the most devastating fungal wheat diseases. During the past decades, many efforts have been deployed to dissect FHB resistance, investigating both the wheat responses to infection and, more recently, the fungal determinants of pathogenicity. Although no total resistance has been identified so far, they demonstrated that some plant functions and the expression of specific genes are needed to promote FHB. Associated with the increasing list of F. graminearum effectors able to divert plant molecular processes, this fact strongly argues for a functional link between susceptibility-related factors and the fate of this disease in wheat. In this review, we gather more recent data concerning the involvement of plant and fungal genes and the functions and mechanisms in the development of FHB susceptibility, and we discuss the possibility to use them to diversify the current sources of FHB resistance.

A Comparative Transcriptomic and Proteomic Analysis of Hexaploid Wheat's Responses to Colonization by Bacillus velezensis and Gaeumannomyces graminis, Both Separately and Combined

Tritrophic interactions involving a biocontrol agent, a pathogen, and a plant have been analyzed predominantly from the perspective of the biocontrol agent. To explore the adaptive strategies of wheat in response to beneficial, pathogenic, and combined microorganisms, we performed the first comprehensive transcriptomic, proteomic, and biochemical analysis in wheat roots after exposure to Bacillus velezensis CC09, Gaeumannomyces graminis var. tritici, and their combined colonization, respectively. The transcriptional or translational programming of wheat roots inoculated with beneficial B. velezensis showed mild alterations compared with that of pathogenic G. graminis var. tritici. However, the combination of B. velezensis and G. graminis var. tritici activated a larger transcriptional or translational program than for each single microorganism, although the gene expression pattern was similar to that of individual infection by G. graminis var. tritici, suggesting a prioritization of defense against G. graminis var. tritici infection. Surprisingly, pathogen-associated molecular pattern-triggered immunity and effector-triggered immunity made wheat pretreated with B. velezensis more sensitive to subsequent G. graminis var. tritici infection. Additionally, B. velezensis triggered a salicylic acid (SA)-dependent mode of induced systemic resistance that resembles pathogen-induced systemic acquired resistance. Wheat plants mainly depend on SA-mediated resistance, and not that mediated by jasmonic acid (JA), against the necrotrophic pathogen G. graminis var. tritici. Moreover, SA-JA interactions resulted in antagonistic effects regardless of the type of microorganisms in wheat. Further enhancement of SA-dependent defense responses such as lignification to the combined infection was shown to reduce the level of induced JA-dependent defense against subsequent infection with G. graminis var. tritici. Altogether, our results demonstrate how the hexaploid monocot wheat responds to beneficial or pathogenic microorganisms and prolongs the onset of take-all disease through modulation of cell reprogramming and signaling events.

Time-resolved dissection of the molecular crosstalk driving Fusarium head blight in wheat provides new insights into host susceptibility determinism

Fungal plant diseases are controlled by a complex molecular dialogue that involves pathogen effectors able to manipulate plant susceptibility factors at the earliest stages of the interaction. By probing the wheat-Fusarium graminearum pathosystem, we profiled the coregulations of the fungal and plant proteins shaping the molecular responses of a 96-hr-long infection's dynamics. Although no symptoms were yet detectable, fungal biomass swiftly increased along with an extremely diverse set of secreted proteins and candidate effectors supposed to target key plant organelles. Some showed to be early accumulated during the interaction or already present in spores, otherwise stored in germinating spores and detectable in an in vitro F. graminearum exudate. Wheat responses were swiftly set up and were evidenced before any visible symptom. Significant wheat protein abundance changes co-occurred along with the accumulation of putative secreted fungal proteins and predicted effectors. Regulated wheat proteins were closely connected to basal cellular processes occurring during spikelet ontogeny, and particular coregulation patterns were evidenced between chloroplast proteins and fungal proteins harbouring a predicted chloroplast transit peptide. The described plant and fungal coordinated responses provide a resourceful set of data and expand our understanding of the wheat-F. graminearum interaction.

Transcriptome Sequencing Analysis Provides Insights Into the Response to Fusarium oxysporum in Lilium pumilum

Lily basal rot, caused by Fusarium oxysporum f. sp. lilii, is one of the most serious diseases of lily. Although the lily germplasm which is resistant to F. oxysporum has been used in disease-resistant breeding, few studies on its molecular mechanism of disease resistance have been reported. To comprehensively study the mechanism of resistance to F. oxysporum, transcriptome sequencings of root tissues from Lilium pumilum inoculated with F. oxysporum or sterile water for 6, 12, or 24 h were performed. A total of 50 GB of data were obtained from the transcriptome sequencings of the 6 L. pumilum samples, and 217 098 Unigenes were obtained after the de novo assembly, of which 38.36% Unigenes were annotated. The sequencing results showed that the numbers of differentially expressed genes at 6, 12, and 24 h after inoculation compared with the control were 111, 254, and 2500, respectively. The functional enrichment analysis of the differentially expressed genes showed that several pathways were involved in responses of L. pumilum, mainly including starch and sucrose metabolism, glycolysis/gluconeogenesis, phenylpropanoid biosynthesis, plant hormone signal transduction, flavonoid biosynthesis, vitamin B6 (VB6) biosynthesis, acid biosynthesis, proteasome, and ribosome. Transcription factor analysis revealed that the WRKY and ERF families played important roles in responses of L. pumilum to F. oxysporum. The results of this study elucidate the molecular responses to F. oxysporum in lily and lay a theoretical foundation for improving lily breeding and strategies for lily basal rot resistance.

Proteomics Analysis Reveals That Caspase-Like and Metacaspase-Like Activities Are Dispensable for Activation of Proteases Involved in Early Response to Biotic Stress in Triticum aestivum L

Plants, including Triticum aestivum L., are constantly attacked by various pathogens which induce immune responses. Immune processes in plants are tightly regulated by proteases from different families within their degradome. In this study, a wheat degradome was characterized. Using profile hidden Markov model (HMMer) algorithm and Pfam database, comprehensive analysis of the T. aestivum genome revealed a large number of proteases (1544 in total) belonging to the five major protease families: serine, cysteine, threonine, aspartic, and metallo-proteases. Mass-spectrometry analysis revealed a 30% difference between degradomes of distinct wheat cultivars (Khakasskaya and Darya), and infection by biotrophic (Puccinia recondita Rob. ex Desm f. sp. tritici) or necrotrophic (Stagonospora nodorum) pathogens induced drastic changes in the presence of proteolytic enzymes. This study shows that an early immune response to biotic stress is associated with the same core of proteases from the C1, C48, C65, M24, M41, S10, S9, S8, and A1 families. Further liquid chromatography-mass spectrometry (LC-MS) analysis of the detected protease-derived peptides revealed that infection by both pathogens enhances overall proteolytic activity in wheat cells and leads to activation of proteolytic cascades. Moreover, sites of proteolysis were identified within the proteases, which probably represent targets of autocatalytic activation, or hydrolysis by another protease within the proteolytic cascades. Although predicted substrates of metacaspase-like and caspase-like proteases were similar in biotrophic and necrotrophic infections, proteolytic activation of proteases was not found to be associated with metacaspase-like and caspase-like activities. These findings indicate that the response of T. aestivum to biotic stress is regulated by unique mechanisms.

Degradation of salicylic acid by Fusarium graminearum

Fusarium head blight (FHB) is a major cereal crop disease, caused most frequently by the fungus Fusarium graminearum. We have previously demonstrated that F. graminearum can utilize SA as sole source of carbon to grow. In this current study, we further characterized selected four fungal SA-responsive genes that are predicted to encode salicylic acid (SA)-degrading enzymes and we used a gene replacement approach to characterize them further. These included two genes predicted to encode a salicylate 1-monooxygenase, FGSG_03657 and FGSG_09063, a catechol 1, 2-dioxygenase gene, FGSG_03667, and a 2, 3-dihydroxybenzoic acid decarboxylase gene, FGSG_09061. For each gene, three independent gene replacement strains were assayed for their ability to degrade salicylic acid in liquid culture. Salicylate 1-monooxygenase FGSG_03657 and catechol 1, 2-dioxygenase FGSG_03667 were shown to be essential for SA degradation, while a loss of 2, 3-dihydroxybenzoic acid decarboxylase FGSG_09061 caused only a partial reduction of SA degradation and a loss of salicylate 1-monooxygenase FGSG_09063 had no effect when compared to wild type culture. Salicylate 1-monooxygenase FGSG_03657 and catechol 1, 2-dioxygenase FGSG_03667 were identified as the first two key enzyme steps of SA degradation via catechol in the β-ketoadipate pathway. Expression profiles for all four genes were also determined in liquid culture and in planta. Salicylate 1-monooxygenase FGSG_03657 and catechol 1, 2-dioxygenase FGSG_03667 were co-expressed and their expression was substrate dependent in liquid culture; however their expression was uncoupled in planta. Disruption of the gene for catechol 1, 2-dioxygenase FGSG_03667 was shown to have no effect on fungal virulence on wheat. Our results with 2, 3-dihydroxybenzoic acid decarboxylase FGSG_09061 raise the possibility of an alternate non-oxidative decarboxylation pathway for the conversion of SA to catechol via 2, 3-dihydrozybenzoic acid and for a connection between the oxidative and the non-oxidative decarboxylation pathways for SA conversion.

H3K4/K9 acetylation and Lr28-mediated expression of six leaf rust responsive genes in wheat (Triticum aestivum)

Development of leaf rust-resistant cultivars is a priority during wheat breeding, since leaf rust causes major losses in yield. Resistance against leaf rust due to Lr genes is partly controlled by epigenetic modifications including histone acetylation that is known to respond to biotic/abiotic stresses. In the present study, enrichment of H3K4ac and H3K9ac in promoters of six defense responsive genes (N-acetyltransferase, WRKY 40, WRKY 70, ASR1, Peroxidase 12 and Sarcosine oxidase) was compared with their expression in a pair of near-isogenic lines (NILs) for the gene Lr28 following inoculation with leaf rust pathotype '77-5'; ChIP-qPCR was used for this purpose. The proximal and distal promoters of these genes contained a number of motifs that are known to respond to biotic stresses. The enrichment of two acetylation marks changed with passage of time; changes in expression of two of the six genes (N-acetyltransferase and peroxidase12), largely matched with changes in H3K4/H3K9 acetylation patterns of the two promoter regions. For example, enrichment of both the marks matched with higher expression of N-acetyltransferase gene in susceptible NIL and the deacetylation (H3K4ac) largely matched with reduced gene expression in resistant NIL. In peroxidase12, enrichment of H3K4ac and H3K9ac largely matched with higher expression in both the NILs. In the remaining four genes, changes in H3 acetylation did not always match with gene expression levels. This indicated complexity in the regulation of the expression of these remaining four genes, which may be controlled by other epigenetic/genetic regulatory mechanisms that need further analysis.

Transcriptome dynamics associated with resistance and susceptibility against fusarium head blight in four wheat genotypes

Background

Fusarium head blight (FHB) of wheat in North America is caused mostly by the fungal pathogen Fusarium graminearum (Fg). Upon exposure to Fg, wheat initiates a series of cellular responses involving massive transcriptional reprogramming. In this study, we analyzed transcriptomics data of four wheat genotypes (Nyubai, Wuhan 1, HC374, and Shaw), at 2 and 4 days post inoculation (dpi) with Fg, using RNA-seq technology.

Results

A total of 37,772 differentially expressed genes (DEGs) were identified, 28,961 from wheat and 8811 from the pathogen. The susceptible genotype Shaw exhibited the highest number of host and pathogen DEGs, including 2270 DEGs associating with FHB susceptibility. Protein serine/threonine kinases and LRR-RK were associated with susceptibility at 2 dpi, while several ethylene-responsive, WRKY, Myb, bZIP and NAC-domain containing transcription factors were associated with susceptibility at 4 dpi. In the three resistant genotypes, 220 DEGs were associated with resistance. Glutathione S-transferase (GST), membrane proteins and distinct LRR-RKs were associated with FHB resistance across the three genotypes. Genes with unique, high up-regulation by Fg in Wuhan 1 were mostly transiently expressed at 2 dpi, while many defense-associated genes were up-regulated at both 2 and 4 dpi in Nyubai; the majority of unique genes up-regulated in HC374 were detected at 4 dpi only. In the pathogen, most genes showed increased expression between 2 and 4 dpi in all genotypes, with stronger levels in the susceptible host; however two pectate lyases and a hydrolase were expressed higher at 2 dpi, and acetyltransferase activity was highly enriched at 4 dpi.

Conclusions

There was an early up-regulation of LRR-RKs, different between susceptible and resistant genotypes; subsequently, distinct sets of genes associated with defense response were up-regulated. Differences in expression profiles among the resistant genotypes indicate genotype-specific defense mechanisms. This study also shows a greater resemblance in transcriptomics of HC374 to Nyubai, consistent with their sharing of two FHB resistance QTLs on 3BS and 5AS, compared to Wuhan 1 which carries one QTL on 2DL in common with HC374.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5012-3) contains supplementary material, which is available to authorized users.

Transcriptome dynamics associated with resistance and susceptibility against fusarium head blight in four wheat genotypes

BACKGROUND

Fusarium head blight (FHB) of wheat in North America is caused mostly by the fungal pathogen Fusarium graminearum (Fg). Upon exposure to Fg, wheat initiates a series of cellular responses involving massive transcriptional reprogramming. In this study, we analyzed transcriptomics data of four wheat genotypes (Nyubai, Wuhan 1, HC374, and Shaw), at 2 and 4 days post inoculation (dpi) with Fg, using RNA-seq technology.

RESULTS

A total of 37,772 differentially expressed genes (DEGs) were identified, 28,961 from wheat and 8811 from the pathogen. The susceptible genotype Shaw exhibited the highest number of host and pathogen DEGs, including 2270 DEGs associating with FHB susceptibility. Protein serine/threonine kinases and LRR-RK were associated with susceptibility at 2 dpi, while several ethylene-responsive, WRKY, Myb, bZIP and NAC-domain containing transcription factors were associated with susceptibility at 4 dpi. In the three resistant genotypes, 220 DEGs were associated with resistance. Glutathione S-transferase (GST), membrane proteins and distinct LRR-RKs were associated with FHB resistance across the three genotypes. Genes with unique, high up-regulation by Fg in Wuhan 1 were mostly transiently expressed at 2 dpi, while many defense-associated genes were up-regulated at both 2 and 4 dpi in Nyubai; the majority of unique genes up-regulated in HC374 were detected at 4 dpi only. In the pathogen, most genes showed increased expression between 2 and 4 dpi in all genotypes, with stronger levels in the susceptible host; however two pectate lyases and a hydrolase were expressed higher at 2 dpi, and acetyltransferase activity was highly enriched at 4 dpi.

CONCLUSIONS

There was an early up-regulation of LRR-RKs, different between susceptible and resistant genotypes; subsequently, distinct sets of genes associated with defense response were up-regulated. Differences in expression profiles among the resistant genotypes indicate genotype-specific defense mechanisms. This study also shows a greater resemblance in transcriptomics of HC374 to Nyubai, consistent with their sharing of two FHB resistance QTLs on 3BS and 5AS, compared to Wuhan 1 which carries one QTL on 2DL in common with HC374.

Transcriptomics of cereal-Fusarium graminearum interactions: what we have learned so far

The ascomycete fungal pathogen Fusarium graminearum causes the globally important Fusarium head blight (FHB) disease on cereal hosts, such as wheat and barley. In addition to reducing grain yield, infection by this pathogen causes major quality losses. In particular, the contamination of food and feed with the F. graminearum trichothecene toxin deoxynivalenol (DON) can have many adverse short- and long-term effects on human and animal health. During the last decade, the interaction between F. graminearum and both cereal and model hosts has been extensively studied through transcriptomic analyses. In this review, we present an overview of how such analyses have advanced our understanding of this economically important plant-microbe interaction. From a host point of view, the transcriptomes of FHB-resistant and FHB-susceptible cereal genotypes, including near-isogenic lines (NILs) that differ by the presence or absence of quantitative trait loci (QTLs), have been studied to understand the mechanisms of disease resistance afforded by such QTLs. Transcriptomic analyses employed to dissect host responses to DON have facilitated the identification of the genes involved in toxin detoxification and disease resistance. From the pathogen point of view, the transcriptome of F. graminearum during pathogenic vs. saprophytic growth, or when infecting different cereal hosts or different tissues of the same host, have been studied. In addition, comparative transcriptomic analyses of F. graminearum knock-out mutants with altered virulence have provided new insights into pathogenicity-related processes. The F. graminearum transcriptomic data generated over the years are now being exploited to build a systems level understanding of the biology of this pathogen, with an ultimate aim of developing effective and sustainable disease prevention strategies.

Comparative Transcriptome Profiles of Near-Isogenic Hexaploid Wheat Lines Differing for Effective Alleles at the 2DL FHB Resistance QTL

Fusarium head blight (FHB), caused by the fungus Fusarium graminearum, represents one of the major wheat diseases worldwide, determining severe yield losses and reduction of grain quality due to the accumulation of mycotoxins. The molecular response associated with the wheat 2DL FHB resistance QTL was mined through a comprehensive transcriptomic analysis of the early response to F. graminearum infection, at 3 days post-inoculation, in spikelets and rachis. The analyses were conducted on two near isogenic lines (NILs) differing for the presence of the 2DL QTL (2-2618, resistant 2DL+ and 2-2890, susceptible null). The general response to fungal infection in terms of mRNAs accumulation trend was similar in both NILs, even though involving an higher number of DEGs in the susceptible NIL, and included down-regulation of the primary and energy metabolism, up-regulation of enzymes implicated in lignin and phenylpropanoid biosynthesis, activation of hormons biosynthesis and signal transduction pathways and genes involved in redox homeostasis and transcriptional regulation. The search for candidate genes with expression profiles associated with the 2DL QTL for FHB resistance led to the discovery of processes differentially modulated in the R and S NILs related to cell wall metabolism, sugar and JA signaling, signal reception and transduction, regulation of the redox status and transcription factors. Wheat FHB response-related miRNAs differentially regulated were also identified as putatively implicated in the superoxide dismutase activities and affecting genes regulating responses to biotic/abiotic stresses and auxin signaling. Altered gene expression was also observed for fungal non-codingRNAs. The putative targets of two of these were represented by the wheat gene WIR1A, involved in resistance response, and a gene encoding a jacalin-related lectin protein, which participate in biotic and abiotic stress response, supporting the presence of a cross-talk between the plant and the fungus.

Development and comparative genomic mapping of Dasypyrum villosum 6V#4S-specific PCR markers using transcriptome data

Twenty-five Dasypyrum villosum 6V#4S-specific PCR markers were developed using transcriptome data and further assigned to comparative genomic maps of wheat chromosome 6A, 6B, and 6D and barley chromosome 6H contrasting their homologous genes in these genomes. Two Dasypyrum villosum accessions, D.v#2 and No. 1026 from England and Russia, respectively, contain Pm21 on chromosome 6V#2S and PmV on chromosome 6V#4S. Both genes confer high resistance to powdery mildew (PM) in wheat. Even though several molecular markers have been developed to detect Pm21 and PmV, only the MBH1 marker can simultaneously detect both Pm21 and PmV. In this study, we first used a high-throughput sequencing technique to obtain the transcriptome sequences of a wheat-D. villosum translocation line, Pm97033-which contains chromosome 6V#4S carrying the PmV locus, under wheat PM pathogen induction. Twenty-five 6V#4S chromosome-specific markers were developed. Three of them were able to clearly distinguish chromosomes 6V#4S and 6V#2S by product size, four amplified the product specific for chromosome 6V#4S only, and the remaining 18 markers identified chromosome 6VS in wheat backgrounds. Two different D. villosum accessions, their derived translocation lines and wheat varieties carrying different chromosome 6VS were identified using these specific markers. The 25 newly developed markers together with the known PM resistance gene Stpk-V were used to construct comparative genomic maps with the homoeologous chromosome arms of wheat and barley. The colinearity of the identified gene sequences amplified by the 25 markers among wheat chromosomes 6A, 6B, and 6D and barley chromosome 6H was not very conserved and interrupted frequently by inversion and insertion. Our markers have potential in marker assisted selection for PM resistance breeding, and for locating other potential important genes and cloning the PmV gene on chromosome 6V#4S.

The Fusarium crown rot pathogen Fusarium pseudograminearum triggers a suite of transcriptional and metabolic changes in bread wheat (Triticum aestivum L.)

Background and Aims

Fusarium crown rot caused by the fungal pathogen Fusarium pseudograminearum is a disease of wheat and barley, bearing significant economic cost. Efforts to develop effective resistance to this disease have been hampered by the quantitative nature of resistance and a lack of understanding of the factors associated with resistance and susceptibility. Here, we aimed to dissect transcriptional responses triggered in wheat by F. pseudograminearum infection.

Methods

We used an RNA-seq approach to analyse host responses during a compatible interaction and identified >2700 wheat genes differentially regulated after inoculation with F. pseudograminearum . The production of a few key metabolites and plant hormones in the host during the interaction was also analysed.

Key Results

Analysis of gene ontology enrichment showed that a disproportionate number of genes involved in primary and secondary metabolism, signalling and transport were differentially expressed in infected seedlings. A number of genes encoding pathogen-responsive uridine-diphosphate glycosyltransferases (UGTs) potentially involved in detoxification of the Fusarium mycotoxin deoxynivalenol (DON) were differentially expressed. Using a F. pseudograminearum DON-non-producing mutant, DON was shown to play an important role in virulence during Fusarium crown rot. An over-representation of genes involved in the phenylalanine, tryptophan and tyrosine biosynthesis pathways was observed. This was confirmed through metabolite analyses that demonstrated tryptamine and serotonin levels are induced after F. pseudograminearum inoculation.

Conclusions

Overall, the observed host response in bread wheat to F. pseudograminearum during early infection exhibited enrichment of processes related to pathogen perception, defence signalling, transport and metabolism and deployment of chemical and enzymatic defences. Additional functional analyses of candidate genes should reveal their roles in disease resistance or susceptibility. Better understanding of host responses contributing to resistance and/or susceptibility will aid the development of future disease improvement strategies against this important plant pathogen.

Enabling Molecular Technologies for Trait Improvement in Wheat

Wheat is the major staple food crop and a source of calories for humans worldwide. A steady increase in the wheat production is essential to meet the demands of an ever-increasing global population and to achieve food security. The large size and structurally intricate genome of polyploid wheat had hindered the genomic analysis. However, with the advent of new genomic technologies such as next generation sequencing has led to genome drafts for bread wheat and its progenitors and has paved the way to design new strategies for crop improvement. Here we provide an overview of the advancements made in wheat genomics together with the available "omics approaches" and bioinformatics resources developed for wheat research. Advances in genomic, transcriptomic, and metabolomic technologies are highlighted as options to circumvent existing bottlenecks in the phenotypic and genomic selection and gene transfer. The contemporary reverse genetics approaches, including the novel genome editing techniques to inform targeted manipulation of a single/multiple genes and strategies for generating marker-free transgenic wheat plants, emphasize potential to revolutionize wheat improvement shortly.

Comparison of Fusarium graminearum Transcriptomes on Living or Dead Wheat Differentiates Substrate-Responsive and Defense-Responsive Genes

Fusarium graminearum is an opportunistic pathogen of cereals where it causes severe yield losses and concomitant mycotoxin contamination of the grains. The pathogen has mixed biotrophic and necrotrophic (saprophytic) growth phases during infection and the regulatory networks associated with these phases have so far always been analyzed together. In this study we compared the transcriptomes of fungal cells infecting a living, actively defending plant representing the mixed live style (pathogenic growth on living flowering wheat heads) to the response of the fungus infecting identical, but dead plant tissues (cold-killed flowering wheat heads) representing strictly saprophytic conditions. We found that the living plant actively suppressed fungal growth and promoted much higher toxin production in comparison to the identical plant tissue without metabolism suggesting that molecules signaling secondary metabolite induction are not pre-existing or not stable in the plant in sufficient amounts before infection. Differential gene expression analysis was used to define gene sets responding to the active or the passive plant as main impact factor and driver for gene expression. We correlated our results to the published F. graminearum transcriptomes, proteomes, and secretomes and found that only a limited number of in planta- expressed genes require the living plant for induction but the majority uses simply the plant tissue as signal. Many secondary metabolite (SM) gene clusters show a heterogeneous expression pattern within the cluster indicating that different genetic or epigenetic signals govern the expression of individual genes within a physically linked cluster. Our bioinformatic approach also identified fungal genes which were actively repressed by signals derived from the active plant and may thus represent direct targets of the plant defense against the invading pathogen.

Transcriptome dynamics of a susceptible wheat upon Fusarium head blight reveals that molecular responses to Fusarium graminearum infection fit over the grain development processes

In many plant/pathogen interactions, host susceptibility factors are key determinants of disease development promoting pathogen growth and spreading in plant tissues. In the Fusarium head blight (FHB) disease, the molecular basis of wheat susceptibility is still poorly understood while it could provide new insights into the understanding of the wheat/Fusarium graminearum (Fg) interaction and guide future breeding programs to produce cultivars with sustainable resistance. To identify the wheat grain candidate genes, a genome-wide gene expression profiling was performed in the French susceptible wheat cultivar, Recital. Gene-specific two-way ANOVA of about 40 K transcripts at five grain developmental stages identified 1309 differentially expressed genes. Out of these, 536 were impacted by the Fg effect alone. Most of these Fg-responsive genes belonged to biological and molecular functions related to biotic and abiotic stresses indicating the activation of common stress pathways during susceptibility response of wheat grain to FHB. This analysis revealed also 773 other genes displaying either specific Fg-responsive profiles along with grain development stages or synergistic adjustments with the grain development effect. These genes were involved in various molecular pathways including primary metabolism, cell death, and gene expression reprogramming. An increasingly complex host response was revealed, as was the impact of both Fg infection and grain ontogeny on the transcription of wheat genes. This analysis provides a wealth of candidate genes and pathways involved in susceptibility responses to FHB and depicts new clues to the understanding of the susceptibility determinism in plant/pathogen interactions.

Transcriptome dynamics of a susceptible wheat upon Fusarium head blight reveals that molecular responses to Fusarium graminearum infection fit over the grain development processes

In many plant/pathogen interactions, host susceptibility factors are key determinants of disease development promoting pathogen growth and spreading in plant tissues. In the Fusarium head blight (FHB) disease, the molecular basis of wheat susceptibility is still poorly understood while it could provide new insights into the understanding of the wheat/Fusarium graminearum (Fg) interaction and guide future breeding programs to produce cultivars with sustainable resistance. To identify the wheat grain candidate genes, a genome-wide gene expression profiling was performed in the French susceptible wheat cultivar, Recital. Gene-specific two-way ANOVA of about 40 K transcripts at five grain developmental stages identified 1309 differentially expressed genes. Out of these, 536 were impacted by the Fg effect alone. Most of these Fg-responsive genes belonged to biological and molecular functions related to biotic and abiotic stresses indicating the activation of common stress pathways during susceptibility response of wheat grain to FHB. This analysis revealed also 773 other genes displaying either specific Fg-responsive profiles along with grain development stages or synergistic adjustments with the grain development effect. These genes were involved in various molecular pathways including primary metabolism, cell death, and gene expression reprogramming. An increasingly complex host response was revealed, as was the impact of both Fg infection and grain ontogeny on the transcription of wheat genes. This analysis provides a wealth of candidate genes and pathways involved in susceptibility responses to FHB and depicts new clues to the understanding of the susceptibility determinism in plant/pathogen interactions.

RNA-seq Transcriptome Response of Flax (Linum usitatissimum L.) to the Pathogenic Fungus Fusarium oxysporum f. sp. lini

Fusarium oxysporum f. sp. lini is a hemibiotrophic fungus that causes wilt in flax. Along with rust, fusarium wilt has become an important factor in flax production worldwide. Resistant flax cultivars have been used to manage the disease, but the resistance varies, depending on the interactions between specific cultivars and isolates of the pathogen. This interaction has a strong molecular basis, but no genomic information is available on how the plant responds to attempted infection, to inform breeding programs on potential candidate genes to evaluate or improve resistance across cultivars. In the current study, disease progression in two flax cultivars [Crop Development Center (CDC) Bethune and Lutea], showed earlier disease symptoms and higher susceptibility in the later cultivar. Chitinase gene expression was also divergent and demonstrated and earlier molecular response in Lutea. The most resistant cultivar (CDC Bethune) was used for a full RNA-seq transcriptome study through a time course at 2, 4, 8, and 18 days post-inoculation (DPI). While over 100 genes were significantly differentially expressed at both 4 and 8 DPI, the broadest deployment of plant defense responses was evident at 18 DPI with transcripts of more than 1,000 genes responding to the treatment. These genes evidenced a reception and transduction of pathogen signals, a large transcriptional reprogramming, induction of hormone signaling, activation of pathogenesis-related genes, and changes in secondary metabolism. Among these, several key genes that consistently appear in studies of plant-pathogen interactions, had increased transcript abundance in our study, and constitute suitable candidates for resistance breeding programs. These included: an induced RPMI-induced protein kinase; transcription factors WRKY3, WRKY70, WRKY75, MYB113, and MYB108; the ethylene response factors ERF1 and ERF14; two genes involved in auxin/glucosinolate precursor synthesis (CYP79B2 and CYP79B3); the flavonoid-related enzymes chalcone synthase, dihydroflavonol reductase and multiple anthocyanidin synthases; and a peroxidase implicated in lignin formation (PRX52). Additionally, regulation of some genes indicated potential pathogen manipulation to facilitate infection; these included four disease resistance proteins that were repressed, indole acetic acid amido/amino hydrolases which were upregulated, activated expansins and glucanases, amino acid transporters and aquaporins, and finally, repression of major latex proteins.