Wheat transcriptome profiling reveals abscisic and gibberellic acid treatments regulate early-stage phytohormone defense signaling, cell wall fortification, and metabolic switches following Fusarium graminearum-challenge

A bird's-eye view: exploration of the flavin-containing monooxygenase superfamily in common wheat

The Flavin Monooxygenase (FMO) gene superfamily in plants is involved in various processes most widely documented for its involvement in auxin biosynthesis, specialized metabolite biosynthesis, and plant microbial defense signaling. The roles of FMOs in defense signaling and disease resistance have recently come into focus as they may present opportunities to increase immune responses in plants including leading to systemic acquired resistance, but are not well characterized. We present a comprehensive catalogue of FMOs found in genomes across vascular plants and explore, in depth, 170 wheat TaFMO genes for sequence architecture, cis-acting regulatory elements, and changes due to Transposable Element insertions. A molecular phylogeny separates TaFMOs into three clades (A, B, and C) for which we further report gene duplication patterns, and differential rates of homoeologue expansion and retention among TaFMO subclades. We discuss Clade B TaFMOs where gene expansion is similarly seen in other cereal genomes. Transcriptome data from various studies point towards involvement of subclade B2 TaFMOs in disease responses against both biotrophic and necrotrophic pathogens, substantiated by promoter element analysis. We hypothesize that certain TaFMOs are responsive to both abiotic and biotic stresses, providing potential targets for enhancing disease resistance, plant yield and other important agronomic traits. Altogether, FMOs in wheat and other crop plants present an untapped resource to be exploited for improving the quality of crops.

Sesquiterpenes of the ectomycorrhizal fungus Pisolithus microcarpus alter root growth and promote host colonization

Trees form symbioses with ectomycorrhizal (ECM) fungi, maintained in part through mutual benefit to both organisms. Our understanding of the signaling events leading to the successful interaction between the two partners requires further study. This is especially true for understanding the role of volatile signals produced by ECM fungi. Terpenoids are a predominant class of volatiles produced by ECM fungi. While several ECM genomes are enriched in the enzymes responsible for the production of these volatiles (i.e., terpene synthases (TPSs)) when compared to other fungi, we have limited understanding of the biochemical products associated with each enzyme and the physiological impact of specific terpenes on plant growth. Using a combination of phylogenetic analyses, RNA sequencing, and functional characterization of five TPSs from two distantly related ECM fungi (Laccaria bicolor and Pisolithus microcarpus), we investigated the role of these secondary metabolites during the establishment of symbiosis. We found that despite phylogenetic divergence, these TPSs produced very similar terpene profiles. We focused on the role of P. microcarpus terpenes and found that the fungus expressed a diverse array of mono-, di-, and sesquiterpenes prior to contact with the host. However, these metabolites were repressed following physical contact with the host Eucalyptus grandis. Exposure of E. grandis to heterologously produced terpenes (enriched primarily in γ -cadinene) led to a reduction in the root growth rate and an increase in P. microcarpus-colonized root tips. These results support a very early putative role of fungal-produced terpenes in the establishment of symbiosis between mycorrhizal fungi and their hosts.

Comparative Transcriptome Sequencing and Endogenous Phytohormone Content of Annual Grafted Branches of Zelkova schneideriana and Its Dwarf Variety HenTianGao

Zelkova schneideriana is a fast-growing tree species endemic to China. Recent surveys and reports have highlighted a continued decline in its natural populations; therefore, it is included in the Red List of Threatened Species by The International Union for Conservation of Nature. A new variety "HenTianGao" (H) has been developed with smaller plant height, slow growth, and lower branching points. In this study, we attempted to understand the differences in plant height of Z. schneideriana (J) and its dwarf variety H. We determined the endogenous hormone content in the annual grafted branches of both J and H. J exhibited higher gibberellic acid (GA)-19 and trans-Zeatin (tZ) levels, whereas H had higher levels of indole-3-acetic acid (IAA) catabolite 2-oxindole-3-acetic acid (OxIAA), IAA-Glu conjugate, and jasmonic acid (JA) (and its conjugate JA-Ile). The transcriptome comparison showed differential regulation of 20,944 genes enriched in growth and development, signaling, and metabolism-related pathways. The results show that the differential phytohormone level (IAA, JA, tZ, and GA) was consistent with the expression of the genes associated with their biosynthesis. The differences in relative OxIAA, IAA-Glu, GA19, trans-Zeatin, JA, and JA-Ile levels were linked to changes in respective signaling-related genes. We also observed significant differences in the expression of cell size, number, proliferation, cell wall biosynthesis, and remodeling-related genes in J and H. The differences in relative endogenous hormone levels, expression of biosynthesis, and signaling genes provide a theoretical basis for understanding the plant height differences in Z. schneideriana.

Barley shows reduced Fusarium head blight under drought and modular expression of differentially expressed genes under combined stress

Plants often face simultaneous abiotic and biotic stress conditions; however, physiological and transcriptional responses under such combined stress conditions are still not fully understood. Spring barley (Hordeum vulgare) is susceptible to Fusarium head blight (FHB), which is strongly affected by weather conditions. We therefore studied the potential influence of drought on FHB severity and plant responses in three varieties of different susceptibility. We found strongly reduced FHB severity in susceptible varieties under drought. The number of differentially expressed genes (DEGs) and strength of transcriptomic regulation reflected the concentrations of physiological stress markers such as abscisic acid or fungal DNA contents. Infection-related gene expression was associated with susceptibility rather than resistance. Weighted gene co-expression network analysis revealed 18 modules of co-expressed genes that reflected the pathogen- or drought-response in the three varieties. A generally infection-related module contained co-expressed genes for defence, programmed cell death, and mycotoxin detoxification, indicating that the diverse genotypes used a similar defence strategy towards FHB, albeit with different degrees of success. Further, DEGs showed co-expression in drought- or genotype-associated modules that correlated with measured phytohormones or the osmolyte proline. The combination of drought stress with infection led to the highest numbers of DEGs and resulted in a modular composition of the single-stress responses rather than a specific transcriptional output.

Regulation of hormone pathways in wheat infested by Blumeria graminis f. sp. tritici

BACKGROUND

Wheat powdery mildew is an obligate biotrophic pathogen infecting wheat, which can pose a serious threat to wheat production. In this study, transcriptome sequencing was carried out on wheat leaves infected by Blumeria graminis f. sp. tritici from 0 h to 7 d.

RESULTS

KEGG and GO enrichment analysis revealed that the upstream biosynthetic pathways and downstream signal transduction pathways of salicylic acid, jasmonic acid, and ethylene were highly enriched at all infection periods. Trend analysis showed that the expressions of hormone-related genes were significantly expressed from 1 to 4 d, suggesting that 1 d-4 d is the main period in which hormones play a defensive role. During this period of time, the salicylic acid pathway was up-regulated, while the jasmonic acid and ethylene pathways were suppressed. Meanwhile, four key modules and 11 hub genes were identified, most of which were hormone related.

CONCLUSION

This study improves the understanding of the dynamical responses of wheat to Blumeria graminis f. sp. tritici infestation at the transcriptional level and provides a reference for screening core genes regulated by hormones.

Herbicide 2,4-dichlorophenoxyacetic acid interferes with MAP kinase signaling in Fusarium graminearum and is inhibitory to fungal growth and pathogenesis

Plant hormones are important for regulating growth, development, and plant-pathogen interactions. Some of them are inhibitory to growth of fungal pathogens but the underlying mechanism is not clear. In this study, we found that hyphal growth of Fusarium graminearum was significantly reduced by high concentrations of IAA and its metabolically stable analogue 2,4-dichlorophenoxyacetic acid (2,4-D). Besides inhibitory effects on growth rate, treatments with 2,4-D also caused significant reduction in conidiation, conidium germination, and germ tube growth. Treatments with 2,4-D had no obvious effect on sexual reproduction but significantly reduced TRI gene expression, toxisome formation, and DON production. More importantly, treatments with 2,4-D were inhibitory to infection structure formation and pathogenesis at concentrations higher than 100 µM. The presence of 1000 µM 2,4-D almost completely inhibited plant infection and invasive growth. In F. graminearum, 2,4-D induced ROS accumulation and FgHog1 activation but reduced the phosphorylation level of Gpmk1 MAP kinase. Metabolomics analysis showed that the accumulation of a number of metabolites such as glycerol and arabitol was increased by 2,4-D treatment in the wild type but not in the Fghog1 mutant. Transformants expressing the dominant active FgPBS2S451D T455D allele were less sensitive to 2,4-D and had elevated levels of intracellular glycerol and arabitol induced by 2,4-D in PH-1. Taken together, our results showed that treatments with 2,4-D interfere with two important MAP kinase pathways and are inhibitory to hyphal growth, DON biosynthesis, and plant infection in F. graminearum.

Differential Expression Feature Extraction (DEFE): A Case Study in Wheat FHB RNA-Seq Data Analysis

In differential gene expression data analysis, one objective is to identify groups of co-expressed genes from a large dataset in order to detect the association between such a group of genes and an experimental condition. This is often done through a clustering approach, such as k-means or bipartition hierarchical clustering, based on particular similarity measures in the grouping process. In such a dataset, the gene differential expression itself is an innate attribute that can be used in the feature extraction process. For example, in a dataset consisting of multiple treatments versus their controls, the expression of a gene in each treatment would have three possible behaviors, upregulated, downregulated, or unchanged. We present in this chapter, a differential expression feature extraction (DEFE) method by using a string consisting of three numerical values at each character to denote such behavior, i.e., 1 = up, 2 = down, and 0 = unchanged, which results in up to 3^B differential expression patterns across all B comparisons. This approach has been successfully applied in many research projects, and among these, we demonstrate the strength of DEFE in a case study on RNA-sequencing (RNA-seq) data analysis of wheat challenged with the phytopathogenic fungus, Fusarium graminearum. Combinations of multiple schemes of DEFE patterns revealed groups of genes putatively associated with resistance or susceptibility to FHB.

RNA-Seq Data Processing in Plant-Pathogen Interaction System: A Case Study

In RNA-seq data processing, short reads are usually aligned from one species against its own genome sequence; however, in plant-pathogen interaction systems, reads from both host and pathogen samples are blended together. In contrast with single-genome analyses, both pathogen and host reference genomes are involved in the alignment process. In such circumstances, the order in which the alignment is carried out, whether the host or pathogen is aligned first, or if both genomes are aligned simultaneously, influences the read counts of certain genes. This is a problem, especially at advanced infection stages. It is crucial to have an appropriate strategy for aligning the reads to their respective genomes, yet the existing strategies of either sequential or parallel alignment become problematic when mapping mixed reads to their corresponding reference genomes. The challenge lies in the determination of which reads belong to which species, especially when homology exists between the host and pathogen genomes. This chapter proposes a combo-genome alignment strategy, which was compared with existing alignment scenarios. Simulation results demonstrated that the degree of discrepancy in the results is correlated with phylogenetic distance of the two species in the mixture which was attributable to the extent of homology between the two genomes involved. This correlation was also found in the analysis using two real RNA-seq datasets of Fusarium-challenged wheat plants. Comparisons of the three RNA-seq processing strategies on three simulation datasets and two real Fusarium-infected wheat datasets showed that an alignment to a combo-genome, consisting of both host and pathogen genomes, improves mapping quality as compared to sequential alignment procedures.

Fusarium graminearum Ste3 G-Protein Coupled Receptor: A Mediator of Hyphal Chemotropism and Pathogenesis

Fungal hyphal chemotropism has been shown to be a major contributor to host-pathogen interactions. Previous studies on Fusarium species have highlighted the involvement of the Ste2 G-protein-coupled receptor (GPCR) in mediating polarized hyphal growth toward host-released peroxidase. Here, the role of the opposite mating type GPCR, Ste3, is characterized with respect to Fusarium graminearum chemotropism and pathogenicity. Fgste3Δ deletion strains were found to be compromised in the chemotropic response toward peroxidase, development of lesions on germinating wheat, and infection of Arabidopsis thaliana leaves. In the absence of FgSte3 or FgSte2, F. graminearum cells exposed to peroxidase showed no phosphorylation of the cell-wall integrity, mitogen-activated protein kinase pathway component Mgv1. In addition, transcriptomic gene expression profiling yielded a list of genes involved in cellular reorganization, cell wall remodeling, and infection-mediated responses that were differentially modulated by peroxidase when FgSte3 was present. Deletion of FgSte3 yielded the downregulation of genes associated with mycotoxin biosynthesis and appressorium development, compared to the wild-type strain, both in the presence of peroxidase. Together, these findings contribute to our understanding of the mechanism underlying fungal chemotropism and pathogenesis while raising the novel hypothesis that FgSte2 and FgSte3 are interdependent on each other for the mediation of the redirection of hyphal growth in response to host-derived peroxidase. IMPORTANCE Fusarium head blight of wheat, caused by the filamentous fungus Fusarium graminearum, leads to devastating global food shortages and economic losses. Fungal hyphal chemotropism has been shown to be a major contributor to host-pathogen interactions. Here, the role of the opposite mating type GPCR, Ste3, is characterized with respect to F. graminearum chemotropism and pathogenicity. These findings contribute to our understanding of the mechanisms underlying fungal chemotropism and pathogenesis.

Manipulating GA-Related Genes for Cereal Crop Improvement

The global population is projected to experience a rapid increase in the future, which poses a challenge to global food sustainability. The "Green Revolution" beginning in the 1960s allowed grain yield to reach two billion tons in 2000 due to the introduction of semi-dwarfing genes in cereal crops. Semi-dwarfing genes reduce the gibberellin (GA) signal, leading to short plant stature, which improves the lodging resistance and harvest index under modern fertilization practices. Here, we reviewed the literature on the function of GA in plant growth and development, and the role of GA-related genes in controlling key agronomic traits that contribute to grain yield in cereal crops. We showed that: (1) GA is a significant phytohormone in regulating plant development and reproduction; (2) GA metabolism and GA signalling pathways are two key components in GA-regulated plant growth; (3) GA interacts with other phytohormones manipulating plant development and reproduction; and (4) targeting GA signalling pathways is an effective genetic solution to improve agronomic traits in cereal crops. We suggest that the modification of GA-related genes and the identification of novel alleles without a negative impact on yield and adaptation are significant in cereal crop breeding for plant architecture improvement. We observed that an increasing number of GA-related genes and their mutants have been functionally validated, but only a limited number of GA-related genes have been genetically modified through conventional breeding tools and are widely used in crop breeding successfully. New genome editing technologies, such as the CRISPR/Cas9 system, hold the promise of validating the effectiveness of GA-related genes in crop development and opening a new venue for efficient and accelerated crop breeding.

Tissue-specific transcriptome responses to Fusarium head blight and Fusarium root rot

Fusarium head blight (FHB) and Fusarium root rot (FRR) are important diseases of small-grain cereals caused by Fusarium species. While host response to FHB has been subject to extensive study, very little is known about response to FRR and the transcriptome responses of FHB and FRR have not been thoroughly compared. Brachypodium distachyon (Bd) is an effective model for investigating host responses to both FHB and FRR. In this study the transcriptome response of Bd to F. graminearum (Fg) infection of heads and roots was investigated. An RNA-seq analysis was performed on both Bd FHB and FRR during the early infection. Additionally, an RNA-seq analysis was performed on in vitro samples of Fg for comparison with Fg gene expression in planta. Differential gene expression and gene-list enrichment analyses were used to compare FHB and FRR transcriptome responses in both Bd and Fg. Differential expression of selected genes was confirmed using RT-qPCR. Most genes associated with receptor signalling, cell-wall modification, oxidative stress metabolism, and cytokinin and auxin biosynthesis and signalling genes were generally upregulated in FHB or were downregulated in FRR. In contrast, Bd genes involved in jasmonic acid and ethylene biosynthesis and signalling, and antimicrobial production were similarly differentially expressed in both tissues in response to infection. A transcriptome analysis of predicted Fg effectors with the same infected material revealed elevated expression of core tissue-independent genes including cell-wall degradation enzymes and the gene cluster for DON production but also several tissue-dependent genes including those for aurofusarin production and cutin degradation. This evidence suggests that Fg modulates its transcriptome to different tissues of the same host.

Exploration of Mycotoxin Accumulation and Transcriptomes of Different Wheat Cultivars during Fusarium graminearum Infection

Fusarium graminearum is one of the most devastating diseases of wheat worldwide, and can cause Fusarium head blight (FHB). F. graminearum infection and mycotoxin production mainly present in wheat and can be influenced by environmental factors and wheat cultivars. The objectives of this study were to examine the effect of wheat cultivars and interacting conditions of temperature and water activity (aw) on mycotoxin production by two strains of F. graminearum and investigate the response mechanisms of different wheat cultivars to F. graminearum infection. In this regard, six cultivars of wheat spikes under field conditions and three cultivars of post-harvest wheat grains under three different temperature conditions combined with five water activity (aw) conditions were used for F. graminearum infection in our studies. Liquid chromatography tandem mass spectrometry (LC–MS/MS) analysis showed significant differences in the concentration of Fusarium mycotoxins deoxynivalenol (DON) and its derivative deoxynivalenol-3-glucoside (D3G) resulting from wheat cultivars and environmental factors. Transcriptome profiles of wheat infected with F. graminearum revealed the lower expression of disease defense-factor-related genes, such as mitogen-activated protein kinases (MAPK)-encoding genes and hypersensitivity response (HR)-related genes of infected Annong 0711 grains compared with infected Sumai 3 grains. These findings demonstrated the optimal temperature and air humidity resulting in mycotoxin accumulation, which will be beneficial in determining the conditions of the relative level of risk of contamination with FHB and mycotoxins. More importantly, our transcriptome profiling illustrated differences at the molecular level between wheat cultivars with different FHB resistances, which will lay the foundation for further research on mycotoxin biosynthesis of F. graminearum and regulatory mechanisms of wheat to F. graminearum.

Characterization and Expression of the Pirin Gene Family in Triticum aestivum

Pirins are nuclear bicupin proteins, encoded by genes that are one of several gene families that comprise the Cupin superfamily in plants. Pirin genes have been implicated in stress response pathways studied in Arabidopsis and At-Pirin1 has been shown to interact with the heterotrimeric G-protein alpha subunit (GPA1). The aim of this study was to identify the members of the Pirin gene family in Triticum aestivum, to correct their annotations in the whole genome and gain an insight into their tissue-specific expression as well as their response to abiotic and biotic stresses. The Pirin gene family in T. aestivum is comprised of 18 genes that represent six paralogous gene copies, each having an A, B and D homeolog. Expression analysis of the Pirin genes in T. aestivum Illumina RNA-seq libraries, which included sampling from differing tissue types as well as abiotic and biotic stresses, indicates that the members of the Pirin gene family have specialized expression and play a role in stress responses. Pirin gene families are also identified in other monocots including Aegilops tauschii, Hordeum vulgare, Brachypodium distachyon, Oryza sativa, Zea mays, Sorghum bicolor and the dicot Arabidopsis thaliana.

BdGUCD1 and Cyclic GMP Are Required for Responses of Brachypodium distachyon to Fusarium pseudograminearum in the Mechanism Involving Jasmonate

Guanosine 3′,5′-cyclic monophosphate (cGMP) is an important signaling molecule in plants. cGMP and guanylyl cyclases (GCs), enzymes that catalyze the synthesis of cGMP from GTP, are involved in several physiological processes and responses to environmental factors, including pathogen infections. Using in vitro analysis, we demonstrated that recombinant BdGUCD1 is a protein with high guanylyl cyclase activity and lower adenylyl cyclase activity. In Brachypodium distachyon, infection by Fusarium pseudograminearum leads to changes in BdGUCD1 mRNA levels, as well as differences in endogenous cGMP levels. These observed changes may be related to alarm reactions induced by pathogen infection. As fluctuations in stress phytohormones after infection have been previously described, we performed experiments to determine the relationship between cyclic nucleotides and phytohormones. The results revealed that inhibition of cellular cGMP changes disrupts stress phytohormone content and responses to pathogen. The observations made here allow us to conclude that cGMP is an important element involved in the processes triggered as a result of infection and changes in its levels affect jasmonic acid. Therefore, stimuli-induced transient elevation of cGMP in plants may play beneficial roles in priming an optimized response, likely by triggering the mechanisms of feedback control.

Linking Multi-Omics to Wheat Resistance Types to Fusarium Head Blight to Reveal the Underlying Mechanisms

Fusarium head blight (FHB) caused by Fusarium graminearum is a worldwide disease which has destructive effects on wheat production, resulting in severe yield reduction and quality deterioration, while FHB-infected wheat grains are toxic to people and animals due to accumulation of fungal toxins. Although impressive progress towards understanding host resistance has been achieved, our knowledge of the mechanism underlying host resistance is still quite limited due to the complexity of wheat-pathogen interactions. In recent years, disease epidemics, the resistance germplasms and components, the genetic mechanism of FHB, and disease management and control, etc., have been well reviewed. However, the resistance mechanism of FHB is quite complex with Type I, II to V resistances. In this review, we focus on the potential resistance mechanisms by linking different resistance types to multi-omics and emphasize the pathways or genes that may play significant roles in the different types of resistance. Deciphering the complicated mechanism of FHB resistance types in wheat at the integral levels based on multi-omics may help discover the genes or pathways that are critical for different FHB resistance, which could then be utilized and manipulated to improve FHB resistance in wheat breeding programs by using transgenic approaches, gene editing, or marker assisted selection strategies.

Molecular mechanisms underlying multi-level defense responses of horticultural crops to fungal pathogens

The horticultural industry helps to enrich and improve the human diet while contributing to growth of the agricultural economy. However, fungal diseases of horticultural crops frequently occur during pre- and postharvest periods, reducing yields and crop quality and causing huge economic losses and wasted food. Outcomes of fungal diseases depend on both horticultural plant defense responses and fungal pathogenicity. Plant defense responses are highly sophisticated and are generally divided into preformed and induced defense responses. Preformed defense responses include both physical barriers and phytochemicals, which are the first line of protection. Induced defense responses, which include innate immunity (pattern-triggered immunity and effector-triggered immunity), local defense responses, and systemic defense signaling, are triggered to counterstrike fungal pathogens. Therefore, to develop regulatory strategies for horticultural plant resistance, a comprehensive understanding of defense responses and their underlying mechanisms is critical. Recently, integrated multi-omics analyses, CRISPR-Cas9-based gene editing, high-throughput sequencing, and data mining have greatly contributed to identification and functional determination of novel phytochemicals, regulatory factors, and signaling molecules and their signaling pathways in plant resistance. In this review, research progress on defense responses of horticultural crops to fungal pathogens and novel regulatory strategies to regulate induction of plant resistance are summarized, and then the problems, challenges, and future research directions are examined.