Quantitative trait loci-dependent analysis of a gene co-expression network associated with Fusarium head blight resistance in bread wheat (Triticum aestivum L.)

Molecular Investigations to Improve Fusarium Head Blight Resistance in Wheat: An Update Focusing on Multi-Omics Approaches

Fusarium head blight (FHB) is mainly caused by Fusarium graminearum (Fg) and is a very widespread disease throughout the world, leading to severe damage to wheat with losses in both grain yield and quality. FHB also leads to mycotoxin contamination in the infected grains, being toxic to humans and animals. In spite of the continuous advancements to elucidate more and more aspects of FHB host resistance, to date, our knowledge about the molecular mechanisms underlying wheat defense response to this pathogen is not comprehensive, most likely due to the complex wheat-Fg interaction. Recently, due to climate changes, such as high temperature and heavy rainfall, FHB has become more frequent and severe worldwide, making it even more urgent to completely understand wheat defense mechanisms. In this review, after a brief description of the first wheat immune response to Fg, we discuss, for each FHB resistance type, from Type I to Type V resistances, the main molecular mechanisms involved, the major quantitative trait loci (QTLs) and candidate genes found. The focus is on multi-omics research helping discover crucial molecular pathways for each resistance type. Finally, according to the emerging examined studies and results, a wheat response model to Fg attack, showing the major interactions in the different FHB resistance types, is proposed. The aim is to establish a useful reference point for the researchers in the field interested to adopt an interdisciplinary omics approach.

Genome-wide analysis of salicylic acid and jasmonic acid signalling marker gene families in wheat

Jasmonic acid (JA) and salicylic acid (SA) phytohormone pathways are important regulators of stress tolerance. Knowledge regarding the diversity, phylogeny and functionality of wheat genes involved in JA and SA response is limited. Using Arabidopsis, rice and wheat genomic and wheat disease transcriptomic data, we deduced the size, phylogenetic diversity and pathogen-responsiveness of seven hormone-responsive gene families, and thus selected 14 candidates as potential hormone responsive gene markers. Gene-specific expression studies assessed the impact of exogenous JA and SA on their transcriptional activation in leaves of two distinct wheat cultivars. RNAseq data were interrogated to assess their disease responsiveness and tissue-specific expression. This study elucidated the number, phylogeny and pathogen-responsiveness of wheat genes from seven families, including 12 TaAOS, 6 TaJAMyb, 256 TaWRKY group III, 85 TaPR1, 205 TaPR2, 76 TaPR3 and 124 TaPR5. This included the first description of the wheat AOS, JAMyb, PR2, PR3 and PR5 gene families. Gene expression studies delineated TaAOS1-5B and TaJAMyb-4A as JA-responsive in leaves, but not significantly responsive to SA treatment, while TaWRKY45-B was a SA- but not a JA-responsive marker. Other candidate genes were either unresponsive or non-specific to SA or JA. Our findings highlight that all seven gene families are greatly expanded in wheat as compared to other plants (up to 7.6-fold expansion), and demonstrate disparity in the response to biotic stress between some homoeologous and paralogous sequences within these families. The SA- and JA-responsive marker genes identified herein will prove useful tools to monitor these signalling pathways in wheat.

Exploring the genetic architecture of powdery mildew resistance in wheat through QTL meta-analysis

Powdery mildew (PM), caused by Blumeria graminis f. sp. tritici, poses a significant threat to wheat production, necessitating the development of genetically resistant varieties for long-term control. Therefore, exploring genetic architecture of PM in wheat to uncover important genomic regions is an important area of wheat research. In recent years, the utilization of meta-QTL (MQTL) analysis has gained prominence as an essential tool for unraveling the complex genetic architecture underlying complex quantitative traits. The aim of this research was to conduct a QTL meta-analysis to pinpoint the specific genomic regions in wheat responsible for governing PM resistance. This study integrated 222 QTLs from 33 linkage-based studies using a consensus map with 54,672 markers. The analysis revealed 39 MQTLs, refined to 9 high-confidence MQTLs (hcMQTLs) with confidence intervals of 0.49 to 12.94 cM. The MQTLs had an average physical interval of 41.00 Mb, ranging from 0.000048 Mb to 380.71 Mb per MQTL. Importantly, 18 MQTLs co-localized with known resistance genes like Pm2, Pm3, Pm8, Pm21, Pm38, and Pm41. The study identified 256 gene models within hcMQTLs, providing potential targets for marker-assisted breeding and genomic prediction programs to enhance PM resistance. These MQTLs would serve as a foundation for fine mapping, gene isolation, and functional genomics studies, facilitating a deeper understanding of molecular mechanisms. The identification of candidate genes opens up exciting possibilities for the development of PM-resistant wheat varieties after validation.

Identification of candidate genes for Fusarium head blight resistance from QTLs using RIL population in wheat

Fusarium head blight (FHB) stands out as one of the most devastating wheat diseases and leads to significantly grain yield losses and quality reductions in epidemic years. Exploring quantitative trait loci (QTL) for FHB resistance is a critical step for developing new FHB-resistant varieties. We previously constructed a genetic map of unigenes (UG-Map) according to the physical positions using a set of recombinant-inbred lines (RILs) derived from the cross of 'TN18 × LM6' (TL-RILs). Here, the number of diseased spikelets (NDS) and relative disease index (RDI) for FHB resistance were investigated under four environments using TL-RILs, which were distributed across 13 chromosomes. A number of 36 candidate genes for NDS and RDI from of 19 stable QTLs were identified. The average number of candidate genes per QTL was 1.89, with 14 (73.7%), two (10.5%), and three (15.8%) QTLs including one, two, and 3-10 candidate genes, respectively. Among the 24 candidate genes annotated in the reference genome RefSeq v1.1, the homologous genes of seven candidate genes, including TraesCS4B02G227300 for QNds/Rdi-4BL-4553, TraesCS5B02G303200, TraesCS5B02G303300, TraesCS5B02G303700, TraesCS5B02G303800 and TraesCS5B02G304000 for QNds/Rdi-5BL-9509, and TraesCS7A02G568400 for QNds/Rdi-7AL-14499, were previously reported to be related to FHB resistance in wheat, barely or Brachypodium distachyon. These genes should be closely associated with FHB resistance in wheat. In addition, the homologous genes of five genes, including TraesCS1A02G037600LC for QNds-1AS-2225, TraesCS1D02G017800 and TraesCS1D02G017900 for QNds-1DS-527, TraesCS1D02G018000 for QRdi-1DS-575, and TraesCS4B02G227400 for QNds/Rdi-4BL-4553, were involved in plant defense responses against pathogens. These genes should be likely associated with FHB resistance in wheat.

Update on the Basic Understanding of Fusarium graminearum Virulence Factors in Common Wheat Research

Wheat is one of the most important food crops, both in China and worldwide. Wheat production is facing extreme stresses posed by different diseases, including Fusarium head blight (FHB), which has recently become an increasingly serious concerns. FHB is one of the most significant and destructive diseases affecting wheat crops all over the world. Recent advancements in genomic tools provide a new avenue for the study of virulence factors in relation to the host plants. The current review focuses on recent progress in the study of different strains of Fusarium infection. The presence of genome-wide repeat-induced point (RIP) mutations causes genomic mutations, eventually leading to host plant susceptibility against Fusarium invasion. Furthermore, effector proteins disrupt the host plant resistance mechanism. In this study, we proposed systematic modification of the host genome using modern biological tools to facilitate plant resistance against foreign invasion. We also suggested a number of scientific strategies, such as gene cloning, developing more powerful functional markers, and using haplotype marker-assisted selection, to further improve FHB resistance and associated breeding methods.

Genome-wide association analysis of type II resistance to Fusarium head blight in common wheat

Background

Fusarium head blight (FHB) is a disease affecting wheat spikes caused by some Fusarium species and leads to cases of severe yield reduction and seed contamination. Identifying resistance genes/QTLs from wheat germplasm may help to improve FHB resistance in wheat production.

Methods

Our study evaluated 205 elite winter wheat cultivars for FHB resistance. A high-density 90K SNP array was used for genotyping the panel. A genome-wide association study (GWAS) from cultivars from three different environments was performed using a mixed linear model (MLM).

Results

Sixty-six significant marker-trait associations (MTAs) were identified (P < 0.001) on fifteen chromosomes that explained the phenotypic variation ranging from 5.4 to 11.2%. Some important new MTAs in genomic regions involving FHB resistance were found on chromosomes 2A, 3B, 5B, 6A, and 7B. Six MTAs at 92 cM on chromosome 7B were found in cultivars from two different environments. Moreover, there were 11 MTAs consistently associated with diseased spikelet rate and diseased rachis rate as pleiotropic effect loci and D_contig74317_533 on chromosome 5D was novel for FHB resistance. Eight new candidate genes of FHB resistance were predicated in wheat in this study. Three candidate genes, TraesCS5D02G006700, TraesCS6A02G013600, and TraesCS7B02G370700 on chromosome 5DS, 6AS, and 7BL, respectively, were perhaps important in defending against FHB by regulating intramolecular transferase activity, GTP binding, or chitinase activity in wheat, but further validation in needed. In addition, a total of five favorable alleles associated with wheat FHB resistance were discovered. These results provide important genes/loci for enhancing FHB resistance in wheat breeding by marker-assisted selection.

Transcriptional Profiling Reveals the Wheat Defences against Fusarium Head Blight Disease Regulated by a NAC Transcription Factor

The wheat NAC transcription factor TaNACL-D1 enhances resistance to the economically devastating Fusarium head blight (FHB) disease. The objective of this study was to decipher the alterations in gene expression, pathways and biological processes that led to enhanced resistance as a result of the constitutive expression of TaNACL-D1 in wheat. Transcriptomic analysis was used to determine the genes and processes enhanced in wheat due to TaNACL-D1 overexpression, both in the presence and absence of the causal agent of FHB, Fusarium graminearum (0- and 1-day post-treatment). The overexpression of TaNACL-D1 resulted in more pronounced transcriptional reprogramming as a response to fungal infection, leading to the enhanced expression of genes involved in detoxification, immune responses, secondary metabolism, hormone biosynthesis, and signalling. The regulation and response to JA and ABA were differentially regulated between the OE and the WT. Furthermore, the results suggest that the OE may more efficiently: (i) regulate the oxidative burst; (ii) modulate cell death; and (iii) induce both the phenylpropanoid pathway and lignin synthesis. Thus, this study provides insights into the mode of action and downstream target pathways for this novel NAC transcription factor, further validating its potential as a gene to enhance FHB resistance in wheat.

Genome-Wide Identification of Candidate Genes Associated with Heat Stress in Mulberry (Morus alba L.)

Mulberry (Morus alba L.) is an economically important plant for the silk industry and has the possibility of contributing immensely to Chinese pharmacopeia because of its health benefits. Domesticated silkworms feed only on mulberry leaves, meaning that the worms' survival depends on the mulberry tree. Mulberry production is threatened by climate change and global warming. However, the regulatory mechanisms of mulberry responses to heat are poorly understood. We performed transcriptome analysis of high-temperature-stressed (42 °C) M. alba seedlings using RNA-Seq technologies. A total of 703 differentially expressed genes (DEGs) were discovered from 18,989 unigenes. Among these, 356 were up-regulated, and 347 were down-regulated. KEGG analysis revealed that most DEGs were enriched in valine, leucine and isoleucine degradation, and in starch and sucrose metabolism, alpha-linolenic acid metabolism, carotenoid biosynthesis and galactose metabolism, among others. In addition, TFs such as the NAC, HSF, IAA1, MYB, AP2, GATA, WRKY, HLH and TCP families were actively involved in response to high temperatures. Moreover, we used RT-qPCR to confirm the expression changes of eight genes under heat stress observed in the RNA-Seq analysis. This study provides M. alba transcriptome profiles under heat stress and provides theoretical bases to researchers for better understanding mulberry heat response mechanisms and breeding heat-tolerant mulberry plants.

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.

Comprehensive analysis of pathogen-responsive wheat NAC transcription factors: new candidates for crop improvement

Wheat NAC (TaNAC) transcription factors are important regulators of stress responses and developmental processes. This study proposes a new TaNAC nomenclature and identified defense-associated TaNACs based on the analysis of RNA-sequencing datasets of wheat tissue infected with major fungal pathogens. A total of 146 TaNACs were pathogen-responsive, of which 52 were orthologous with functionally characterized defense-associated NACs from barley, rice, and Arabidopsis, as deduced via phylogenetic analysis. Next, we focused on the phylogenetic relationship of the pathogen-responsive TaNACs and their expression profiles in healthy and diseased tissues. Three subfamilies ("a," "e," and "f") were significantly enriched in pathogen-responsive TaNACs, of which the majority were responsive to at least 2 pathogens (universal pathogen response). Uncharacterized TaNACs from subfamily "a" enriched with defense-associated NACs are promising candidates for functional characterization in pathogen defense. In general, pathogen-responsive TaNACs were expressed in at least 2 healthy organs. Lastly, we showed that the wheat NAM domain is significantly divergent in sequence in subfamilies "f," "g," and "h" based on HMMER and motif analysis. New protein motifs were identified in both the N- and C-terminal parts of TaNACs. Three of those identified in the C-terminal part were linked to pathogen responsiveness of the TaNACs and 2 were linked to expression in grain tissue. Future studies should benefit from this comprehensive in silico analysis of pathogen-responsive TaNACs as a basis for selecting the most promising candidates for functional validation and crop improvement.

Meta-analysis reveals consensus genomic regions associated with multiple disease resistance in wheat (Triticum aestivum L.)

In wheat, meta-QTLs (MQTLs) and candidate genes (CGs) were identified for multiple disease resistance (MDR). For this purpose, information was collected from 58 studies for mapping QTLs for resistance to one or more of the five diseases. As many as 493 QTLs were available from these studies, which were distributed in five diseases as follows: septoria tritici blotch (STB) 126 QTLs; septoria nodorum blotch (SNB), 103 QTLs; fusarium head blight (FHB), 184 QTLs; karnal bunt (KB), 66 QTLs; and loose smut (LS), 14 QTLs. Of these 493 QTLs, only 291 QTLs could be projected onto a consensus genetic map, giving 63 MQTLs. The CI of the MQTLs ranged from 0.04 to 15.31 cM with an average of 3.09 cM per MQTL. This is a ~ 4.39 fold reduction from the CI of QTLs, which ranged from 0 to 197.6 cM, with a mean of 13.57 cM. Of 63 MQTLs, 60 were anchored to the reference physical map of wheat (the physical interval of these MQTLs ranged from 0.30 to 726.01 Mb with an average of 74.09 Mb). Thirty-eight (38) of these MQTLs were verified using marker-trait associations (MTAs) derived from genome-wide association studies. As many as 874 CGs were also identified which were further investigated for differential expression using data from five transcriptome studies, resulting in 194 differentially expressed candidate genes (DECGs). Among the DECGs, 85 genes had functions previously reported to be associated with disease resistance. These results should prove useful for fine mapping and cloning of MDR genes and marker-assisted breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01282-z.

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.

Response of Wheat Genotypes to Drought Stress Stimulated by PEG

Wheat is a cereal grain crop that is commonly cultivated and is a good source of nutrients that are beneficial to human health. In recent years, the productivity of wheat has been steadily declining, with abiotic pressures accounting for almost half of all yield losses. Drought stress is a significant limiting factor for plant development and production around the planet. The influence of polyethylene glycol (PEG) (at concentrations of 5, 10, and 15%)-induced drought stress on the morphological, physiological, and biochemical characteristics of fifteen wheat genotypes was investigated in this work. Overall, it was discovered that morphological and physiological indicators such as germination % and shoot-root lengths during the seedling stage had reduced significantly. The proline content, on the other hand, was shown to be positively correlated with the concentration of PEG treatments. There was a significant difference between the genotypes HD2733, HD2888, and RAJ3765 regarding tolerance to abiotic stress caused by drought. A further finding was that under stressful settings, the first three main components explained 56.65 percent, 65.06 percent, and 72.47 percent of the total variability in PEG treatment levels of five, ten, and fifteen percent, respectively. These collective morphological and physiological parameters, and analyses of their diverse responses, could be used for screening of drought tolerance among the 15 wheat genotypes to select for significant drought tolerance and diverse molecular responses during breeding of stress resistant forms.

Interpreting Functional Impact of Genetic Variations by Network QTL for Genotype–Phenotype Association Study

An enormous challenge in the post-genome era is to annotate and resolve the consequences of genetic variation on diverse phenotypes. The genome-wide association study (GWAS) is a well-known method to identify potential genetic loci for complex traits from huge genetic variations, following which it is crucial to identify expression quantitative trait loci (eQTL). However, the conventional eQTL methods usually disregard the systematical role of single-nucleotide polymorphisms (SNPs) or genes, thereby overlooking many network-associated phenotypic determinates. Such a problem motivates us to recognize the network-based quantitative trait loci (QTL), i.e., network QTL (nQTL), which is to detect the cascade association as genotype → network → phenotype rather than conventional genotype → expression → phenotype in eQTL. Specifically, we develop the nQTL framework on the theory and approach of single-sample networks, which can identify not only network traits (e.g., the gene subnetwork associated with genotype) for analyzing complex biological processes but also network signatures (e.g., the interactive gene biomarker candidates screened from network traits) for characterizing targeted phenotype and corresponding subtypes. Our results show that the nQTL framework can efficiently capture associations between SNPs and network traits (i.e., edge traits) in various simulated data scenarios, compared with traditional eQTL methods. Furthermore, we have carried out nQTL analysis on diverse biological and biomedical datasets. Our analysis is effective in detecting network traits for various biological problems and can discover many network signatures for discriminating phenotypes, which can help interpret the influence of nQTL on disease subtyping, disease prognosis, drug response, and pathogen factor association. Particularly, in contrast to the conventional approaches, the nQTL framework could also identify many network traits from human bulk expression data, validated by matched single-cell RNA-seq data in an independent or unsupervised manner. All these results strongly support that nQTL and its detection framework can simultaneously explore the global genotype–network–phenotype associations and the underlying network traits or network signatures with functional impact and importance.

Bioprospecting Phenols as Inhibitors of Trichothecene-Producing Fusarium: Sustainable Approaches to the Management of Wheat Pathogens

Fusarium spp. are ubiquitous fungi able to cause Fusarium head blight and Fusarium foot and root rot on wheat. Among relevant pathogenic species, Fusarium graminearum and Fusarium culmorum cause significant yield and quality loss and result in contamination of the grain with mycotoxins, mainly type B trichothecenes, which are a major health concern for humans and animals. Phenolic compounds of natural origin are being increasingly explored as fungicides on those pathogens. This review summarizes recent research activities related to the antifungal and anti-mycotoxigenic activity of natural phenolic compounds against Fusarium, including studies into the mechanisms of action of major exogenous phenolic inhibitors, their structure-activity interaction, and the combined effect of these compounds with other natural products or with conventional fungicides in mycotoxin modulation. The role of high-throughput analysis tools to decipher key signaling molecules able to modulate the production of mycotoxins and the development of sustainable formulations enhancing potential inhibitors' efficacy are also discussed.

Genome-Wide Diversity of MADS-Box Genes in Bread Wheat is Associated with its Rapid Global Adaptability

MADS-box gene family members play multifarious roles in regulating the growth and development of crop plants and hold enormous promise for bolstering grain yield potential under changing global environments. Bread wheat (Triticum aestivum L.) is a key stable food crop around the globe. Until now, the available information concerning MADS-box genes in the wheat genome has been insufficient. Here, a comprehensive genome-wide analysis identified 300 high confidence MADS-box genes from the publicly available reference genome of wheat. Comparative phylogenetic analyses with Arabidopsis and rice MADS-box genes classified the wheat genes into 16 distinct subfamilies. Gene duplications were mainly identified in subfamilies containing unbalanced homeologs, pointing towards a potential mechanism for gene family expansion. Moreover, a more rapid evolution was inferred for M-type genes, as compared with MIKC-type genes, indicating their significance in understanding the evolutionary history of the wheat genome. We speculate that subfamily-specific distal telomeric duplications in unbalanced homeologs facilitate the rapid adaptation of wheat to changing environments. Furthermore, our in-silico expression data strongly proposed MADS-box genes as active guardians of plants against pathogen insurgency and harsh environmental conditions. In conclusion, we provide an entire complement of MADS-box genes identified in the wheat genome that could accelerate functional genomics efforts and possibly facilitate bridging gaps between genotype-to-phenotype relationships through fine-tuning of agronomically important traits.

Network organization of the plant immune system: from pathogen perception to robust defense induction

The plant immune system has been explored essentially through the study of qualitative resistance, a simple form of immunity, and from a reductionist point of view. The recent identification of genes conferring quantitative disease resistance revealed a large array of functions, suggesting more complex mechanisms. In addition, thanks to the advent of high-throughput analyses and system approaches, our view of the immune system has become more integrative, revealing that plant immunity should rather be seen as a distributed and highly connected molecular network including diverse functions to optimize expression of plant defenses to pathogens. Here, we review the recent progress made to understand the network complexity of regulatory pathways leading to plant immunity, from pathogen perception, through signaling pathways and finally to immune responses. We also analyze the topological organization of these networks and their emergent properties, crucial to predict novel immune functions and test them experimentally. Finally, we report how these networks might be regulated by environmental clues. Although system approaches remain extremely scarce in this area of research, a growing body of evidence indicates that the plant response to combined biotic and abiotic stresses cannot be inferred from responses to individual stresses. A view of possible research avenues in this nascent biology domain is finally proposed.

Plant PhysioSpace: a robust tool to compare stress response across plant species

Generalization of transcriptomics results can be achieved by comparison across experiments. This generalization is based on integration of interrelated transcriptomics studies into a compendium. Such a focus on the bigger picture enables both characterizations of the fate of an organism and distinction between generic and specific responses. Numerous methods for analyzing transcriptomics datasets exist. Yet, most of these methods focus on gene-wise dimension reduction to obtain marker genes and gene sets for, for example, pathway analysis. Relying only on isolated biological modules might result in missing important confounders and relevant contexts. We developed a method called Plant PhysioSpace, which enables researchers to compute experimental conditions across species and platforms without a priori reducing the reference information to specific gene sets. Plant PhysioSpace extracts physiologically relevant signatures from a reference dataset (i.e. a collection of public datasets) by integrating and transforming heterogeneous reference gene expression data into a set of physiology-specific patterns. New experimental data can be mapped to these patterns, resulting in similarity scores between the acquired data and the extracted compendium. Because of its robustness against platform bias and noise, Plant PhysioSpace can function as an inter-species or cross-platform similarity measure. We have demonstrated its success in translating stress responses between different species and platforms, including single-cell technologies. We have also implemented two R packages, one software and one data package, and a Shiny web application to facilitate access to our method and precomputed models.

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.

Extensive structural variation in the Bowman-Birk inhibitor family in common wheat (Triticum aestivum L.)

BACKGROUND

Bowman-Birk inhibitors (BBI) are a family of serine-type protease inhibitors that modulate endogenous plant proteolytic activities during different phases of development. They also inhibit exogenous proteases as a component of plant defense mechanisms, and their overexpression can confer resistance to phytophagous herbivores and multiple fungal and bacterial pathogens. Dicot BBIs are multifunctional, with a "double-headed" structure containing two separate inhibitory loops that can bind and inhibit trypsin and chymotrypsin proteases simultaneously. By contrast, monocot BBIs have a non-functional chymotrypsin inhibitory loop, although they have undergone internal duplication events giving rise to proteins with multiple BBI domains.

RESULTS

We used a Hidden Markov Model (HMM) profile-based search to identify 57 BBI genes in the common wheat (Triticum aestivum L.) genome. The BBI genes are unevenly distributed, with large gene clusters in the telomeric regions of homoeologous group 1 and 3 chromosomes that likely arose through a series of tandem gene duplication events. The genomes of wheat progenitors also contain contiguous clusters of BBI genes, suggesting this family underwent expansion before the domestication of common wheat. However, the BBI gene family varied in size among different cultivars, showing this family remains dynamic. Because of these expansions, the BBI gene family is larger in wheat than other monocots such as maize, rice and Brachypodium. We found BBI proteins in common wheat with intragenic homologous duplications of cysteine-rich functional domains, including one protein with four functional BBI domains. This diversification may expand the spectrum of target substrates. Expression profiling suggests that some wheat BBI proteins may be involved in regulating endogenous proteases during grain development, while others were induced in response to biotic and abiotic stresses, suggesting a role in plant defense.

CONCLUSIONS

Genome-wide characterization reveals that the BBI gene family in wheat is subject to a high rate of homologous tandem duplication and deletion events, giving rise to a diverse set of encoded proteins. This information will facilitate the functional characterization of individual wheat BBI genes to determine their role in wheat development and stress responses, and their potential application in breeding.

Merging Genomics and Transcriptomics for Predicting Fusarium Head Blight Resistance in Wheat

Genomic selection with genome-wide distributed molecular markers has evolved into a well-implemented tool in many breeding programs during the last decade. The resistance against Fusarium head blight (FHB) in wheat is probably one of the most thoroughly studied systems within this framework. Aside from the genome, other biological strata like the transcriptome have likewise shown some potential in predictive breeding strategies but have not yet been investigated for the FHB-wheat pathosystem. The aims of this study were thus to compare the potential of genomic with transcriptomic prediction, and to assess the merit of blending incomplete transcriptomic with complete genomic data by the single-step method. A substantial advantage of gene expression data over molecular markers has been observed for the prediction of FHB resistance in the studied diversity panel of breeding lines and released cultivars. An increase in prediction ability was likewise found for the single-step predictions, although this can mostly be attributed to an increased accuracy among the RNA-sequenced genotypes. The usage of transcriptomics can thus be seen as a complement to already established predictive breeding pipelines with pedigree and genomic data, particularly when more cost-efficient multiplexing techniques for RNA-sequencing will become more accessible in the future.

Genome-wide identification and analysis of NPR family genes in Brassica juncea var. tumida

Nonexpressor of pathogenesis-related (NPR) genes are bona fide transcription cofactors in the signal transduction pathway of salicylic acid (SA) and play critical regulatory roles in plant immunity. However, the NPR family genes in Brassica juncea var. tumida have not yet been comprehensively identified and analyzed as of yet. In the present study, NPR genes in B. juncea var. tumida seedlings were identified, and the tissue-specific expression patterns of NPR genes in the seedling were analyzed under salt stress (200 mM) treatment and infection by Plasmodiophora brassicae. A total of 19 NPR family genes clustering into six separate groups were identified in the genome of B. juncea var. tumida. These BjuNPR family genes were located in 11 of 18 chromosomes of B. juncea var. tumida and each possessed 1-5 exons. The BjuNPR family members had similar protein structures and conserved motifs. The BjuNPR genes exhibited tissue-specific expression patterns in the root, stem, leaf, flower and pod. Some BjuNPR genes were sensitive to salt stress and showed up-regulated or down-regulated expression patterns and most BjuNPR genes were up-regulated upon infection by P. brassicae. This study provides a foundation for further research into BjuNPR genes regulation in plant growth, development, and abiotic stress tolerance.

Applying genomic resources to accelerate wheat biofortification

Wheat has low levels of the micronutrients iron and zinc in the grain, which contributes to 2 billion people suffering from micronutrient deficiency globally. While wheat flour is commonly fortified during processing, an attractive and more sustainable solution is biofortification, which could improve micronutrient content in the human diet, without the sustainability issues and costs associated with conventional fortification. Although many studies have used quantitative trait loci mapping and genome-wide association to identify genetic loci to improve micronutrient contents, recent developments in genomics offer an opportunity to accelerate marker discovery and use gene-focussed approaches to engineer improved micronutrient content in wheat. The recent publication of a high-quality wheat genome sequence, alongside gene expression atlases, variation datasets and sequenced mutant populations, provides a foundation to identify genetic loci and genes controlling micronutrient content in wheat. We discuss how novel genomic resources can identify candidate genes for biofortification, integrating knowledge from other cereal crops, and how these genes can be tested using gene editing, transgenic and TILLING approaches. Finally, we highlight key challenges remaining to develop wheat cultivars with high levels of iron and zinc.

Mapping and characterization of major QTL for spike traits in common wheat

The spike traits of wheat can directly affect yield. F₂ and F_2:3 lines derived from the cross of the multi-spikelet female 10-A and the uni-spikelet male BE89 were used to detect QTLs for spike length (SL), total spikelet number per spike (TSS), kernel number per spike (KNS) and thousand-kernel weight (TKW) in four different environments. A total of 1098 SNP and 5 SSR were used to construct genetic map of 2398.1 cM with the average distance of 2.2 cM between markers. A total of 11 QTLs were identified for spike traits, including three QTLs for SL, five QTLs for TSS, two QTLs for KNS and one QTL for TKW. The QTLs mapped to chromosomes 2D, 4A, 6A, 7A and 7B explained 8.2-37.8% of the phenotypic variation in single environment. The major QTL confidence interval with distance of 0.5 cM was located on chromosome 4A and detected in multiple environments, which can explain more than 30% of the phenotypic variation for SL, TSS and KNS. Combining IWGSC RefSeq v1.0 and RNA-seq data for 10-A and BE89, we identified 16 genes expressed on spike or grain in four QTL regions. These findings provide insights into improving wheat yield through increasing spikletes in wheat, particularly through the use of the multi-spikelet female 10-A for breeding.

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.

Genome-wide association analyses of quantitative disease resistance in diverse sets of soybean [Glycine max (L.) Merr.] plant introductions

Phytophthora sojae is one of the costliest soybean pathogens in the US. Quantitative disease resistance (QDR) is a vital part of Phytophthora disease management. In this study, QDR was measured in 478 and 495 plant introductions (PIs) towards P. sojae isolates OH.121 and C2.S1, respectively, in genome-wide association (GWA) analyses to identify genetic markers linked to QDR loci (QDRL). Populations were generated by sampling PIs from the US, the Republic of Korea, and the full collection of PIs maintained by the USDA. Additionally, a meta-analysis of QDRL reported from bi-parental studies was done to compare past and present findings. Twenty-four significant marker-trait associations were identified from the 478 PIs phenotyped with OH.121, and an additional 24 marker-trait associations were identified from the 495 PIs phenotyped with C2.S1. In total, 48 significant markers were distributed across 16 chromosomes and based on linkage analysis, represent a total of 44 QDRL. The majority of QDRL were identified with only one of the two isolates, and only a region on chromosome 13 was consistently identified. Regions on chromosomes 3, 13, and 17 were identified in previous GWA-analyses and were re-identified in this study. Five QDRL co-localized with P. sojae meta-QDRL identified from QDRL reported in previous biparental mapping studies. The remaining regions represent novel QDRL, in the soybean-P. sojae pathosystem and were primarily identified in germplasm from the Republic of Korea. Overall, the number of loci identified in this study highlights the complexity of QDR to P. sojae.

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.

Taxonomically Restricted Wheat Genes Interact With Small Secreted Fungal Proteins and Enhance Resistance to Septoria Tritici Blotch Disease

Understanding the nuances of host/pathogen interactions are paramount if we wish to effectively control cereal diseases. In the case of the wheat/Zymoseptoria tritici interaction that leads to Septoria tritici blotch (STB) disease, a 10,000-year-old conflict has led to considerable armaments being developed on both sides which are not reflected in conventional model systems. Taxonomically restricted genes (TRGs) have evolved in wheat to better allow it to cope with stress caused by fungal pathogens, and Z. tritici has evolved specialized effectors which allow it to manipulate its' host. A microarray focused on the latent phase response of a resistant wheat cultivar (cv. Stigg) and susceptible wheat cultivar (cv. Gallant) to Z. tritici infection was mined for TRGs within the Poaceae. From this analysis, we identified two TRGs that were significantly upregulated in response to Z. tritici infection, Septoria-responsive TRG6 and 7 (TaSRTRG6 and TaSRTRG7). Virus induced silencing of these genes resulted in an increased susceptibility to STB disease in cvs. Gallant and Stigg, and significantly so in the latter (2.5-fold increase in STB disease). In silico and localization studies categorized TaSRTRG6 as a secreted protein and TaSRTRG7 as an intracellular protein. Yeast two-hybrid analysis and biofluorescent complementation studies demonstrated that both TaSRTRG6 and TaSRTRG7 can interact with small proteins secreted by Z. tritici (potential effector candidates). Thus we conclude that TRGs are an important part of the wheat-Z. tritici co-evolution story and potential candidates for modulating STB resistance.

Characterization of QTL and eQTL controlling early Fusarium graminearum infection and deoxynivalenol levels in a Wuhan 1 x Nyubai doubled haploid wheat population

BACKGROUND

Fusarium head blight (FHB) is a major disease of cereal crops, caused by the fungal pathogen Fusarium graminearum and related species. Breeding wheat for FHB resistance contributes to increase yields and grain quality and to reduce the use of fungicides. The identification of genes and markers for FHB resistance in different wheat genotypes has nevertheless proven challenging.

RESULTS

In this study, early infection by F. graminearum was analyzed in a doubled haploid population derived from the cross of the moderately resistant wheat genotypes Wuhan 1 and Nyubai. Three quantitative trait loci (QTL) were identified: 1AL was associated with lower deoxynivalenol content, and 4BS and 5A were associated with reduced F. graminearum infection at 2 days post inoculation. Early resistance alleles were inherited from Wuhan 1 for QTL 1AL and 4BS and inherited from Nyubai for the 5A QTL. Cis and trans expression QTL (eQTL) were identified using RNA-seq data from infected head samples. Hotspots for trans eQTL were identified in the vicinity of the 1AL and 4BS QTL peaks. Among differentially expressed genes with cis eQTL within the QTL support intervals, nine genes had higher expression associated with FHB early resistance, and four genes had higher expression associated with FHB early susceptibility.

CONCLUSIONS

Our analysis of genotype and gene expression data of wheat infected by F. graminearum identified three QTL associated with FHB early resistance, and linked genes with eQTL and differential expression patterns to those QTL. These findings may have applications in breeding wheat for early resistance to FHB.

Genome-Wide Identification and Analysis of the NPR1-Like Gene Family in Bread Wheat and Its Relatives

NONEXPRESSOR OF PATHOGENESIS-RELATED GENES 1 (NPR1), and its paralogues NPR3 and NPR4, are bona fide salicylic acid (SA) receptors and play critical regulatory roles in plant immunity. However, comprehensive identification and analysis of the NPR1-like gene family had not been conducted so far in bread wheat and its relatives. Here, a total of 17 NPR genes in Triticum aestivum, five NPR genes in Triticum urartu, 12 NPR genes in Triticum dicoccoides, and six NPR genes in Aegilops tauschii were identified using bioinformatics approaches. Protein properties of these putative NPR1-like genes were also described. Phylogenetic analysis showed that the 40 NPR1-like proteins, together with 40 NPR1-related proteins from other plant species, were clustered into three major clades. The TaNPR1-like genes belonging to the same Arabidopsis subfamilies shared similar exon-intron patterns and protein domain compositions, as well as conserved motifs and amino acid residues. The cis-regulatory elements related to SA were identified in the promoter regions of TaNPR1-like genes. The TaNPR1-like genes were intensively mapped on the chromosomes of homoeologous groups 3, 4, and 5, except TaNPR2-D. Chromosomal distribution and collinearity analysis of NPR1-like genes among bread wheat and its relatives revealed that the evolution of this gene family was more conservative following formation of hexaploid wheat. Transcriptome data analysis indicated that TaNPR1-like genes exhibited tissue/organ-specific expression patterns and some members were induced under biotic stress. These findings lay the foundation for further functional characterization of NPR1-like proteins in bread wheat and its relatives.

A wheat NAC interacts with an orphan protein and enhances resistance to Fusarium head blight disease

Taxonomically-restricted orphan genes play an important role in environmental adaptation, as recently demonstrated by the fact that the Pooideae-specific orphan TaFROG (Triticum aestivum Fusarium Resistance Orphan Gene) enhanced wheat resistance to the economically devastating Fusarium head blight (FHB) disease. Like most orphan genes, little is known about the cellular function of the encoded protein TaFROG, other than it interacts with the central stress regulator TaSnRK1α. Here, we functionally characterized a wheat (T. aestivum) NAC-like transcription factor TaNACL-D1 that interacts with TaFROG and investigated its' role in FHB using studies to assess motif analyses, yeast transactivation, protein-protein interaction, gene expression and the disease response of wheat lines overexpressing TaNACL-D1. TaNACL-D1 is a Poaceae-divergent NAC transcription factor that encodes a Triticeae-specific protein C-terminal region with transcriptional activity and a nuclear localisation signal. The TaNACL-D1/TaFROG interaction was detected in yeast and confirmed in planta, within the nucleus. Analysis of multi-protein interactions indicated that TaFROG could form simultaneously distinct protein complexes with TaNACL-D1 and TaSnRK1α in planta. TaNACL-D1 and TaFROG are co-expressed as an early response to both the causal fungal agent of FHB, Fusarium graminearum and its virulence factor deoxynivalenol (DON). Wheat lines overexpressing TaNACL-D1 were more resistant to FHB disease than wild type plants. Thus, we conclude that the orphan protein TaFROG interacts with TaNACL-D1, a NAC transcription factor that forms part of the disease response evolved within the Triticeae.

Overlapping transcriptional expression response of wheat zinc-induced facilitator-like transporters emphasize important role during Fe and Zn stress

BACKGROUND

Hexaploid wheat is an important cereal crop that has been targeted to enhance grain micronutrient content including zinc (Zn) and iron (Fe). In this direction, modulating the expression of plant transporters involved in Fe and Zn homeostasis has proven to be one of the promising approaches. The present work was undertaken to identify wheat zinc-induced facilitator-like (ZIFL) family of transporters. The wheat ZIFL genes were characterized for their transcriptional expression response during micronutrient fluctuations and exposure to multiple heavy metals.

RESULTS

The genome-wide analyses resulted in identification of fifteen putative TaZIFL-like genes, which were distributed only on Chromosome 3, 4 and 5. Wheat ZIFL proteins subjected to the phylogenetic analysis showed the uniform distribution along with rice, Arabidopsis and maize. In-silico analysis of the promoters of the wheat ZIFL genes demonstrated the presence of multiple metal binding sites including those which are involved in Fe and heavy metal homeostasis. Quantitative real-time PCR analysis of wheat ZIFL genes suggested the differential regulation of the transcripts in both roots and shoots under Zn surplus and also during Fe deficiency. Specifically, in roots, TaZIFL2.3, TaZIFL4.1, TaZIFL4.2, TaZIFL5, TaZIFL6.1 and TaZIFL6.2 were significantly up-regulated by both Zn and Fe. This suggested that ZIFL could possibly be regulated by both the nutrient stress in a tissue specific manner. When exposed to heavy metals, TaZIFL4.2 and TaZIFL7.1 show significant up-regulation, whereas TaZIFL5 and TaZIFL6.2 remained almost unaffected.

CONCLUSION

This is the first report for detailed analysis of wheat ZIFL genes. ZIFL genes also encode for transporter of mugineic acid (TOM) proteins, that are involved in the release of phytosiderophores to enhance Fe/Zn uptake. The detailed expression analysis suggests the varying expression patterns during development of wheat seedlings and also against abiotic/biotic stresses. Overall, this study will lay foundation to prioritize functional assessment of the candidate ZIFL as a putative TOM protein in wheat.

Serpins: Genome-Wide Characterisation and Expression Analysis of the Serine Protease Inhibitor Family in Triticum aestivum

The serine protease inhibitor (serpin) gene family is the largest family of protease inhibitors. Serine protease inhibitors have an active, but under-characterized, role in grain development and defense against pathogen attack in cereal crops. By exploiting publicly available genomic, transcriptomic and proteomic data for wheat (Triticum aestivum), we have identified and annotated the entire ’serpinome’ of wheat and constructed a high-quality and robust phylogenetic tree of the gene family, identifying paralogous and homeologous clades from the hexaploid wheat genome, including the Serpin-Z group that have been well characterized in barley. Using publicly available RNAseq data (http://www.wheat-expression.com/), expression profiles of the wheat serpins were explored across a variety of tissues from the developing grain, spikelet and spike. We show that the SERPIN-Z clade, among others, are highly expressed during grain development, and that there is homeologous and paralogous functional redundancy in this gene family. Further to their role in grain development, serpins play an important but under-explored role in response to fungal pathogens. Using 13 RNAseq datasets of wheat tissues infected by fungal pathogens, we identified 37 serpins with a significant disease response. The majority of the disease-responsive serpins were upregulated by Fusarium graminearum, a destructive fungal pathogen that attacks the spike and developing grain of wheat. As serpins are ubiquitous in wheat grain, the genes encoding serpins may be linked to grain development, with their disease response a result of pleiotropy.

Isobaric tags for relative and absolute quantification-based proteomic analysis of defense responses triggered by the fungal pathogen Fusarium graminearum in wheat

Fusarium head blight (FHB) is a devastating disease worldwide that is predominantly caused by the fungal pathogen Fusarium graminearum. The aim of this work was to study differentially abundant protein species of near-isogenic lines A061-3 and A061-4 with the final goal of elucidating the molecular mechanisms of their differential resistance to F. graminearum. The objectives were accomplished using isobaric tags for relative and absolute quantification (iTRAQ) with mass spectrometry (MS). Lines A061-3 and A061-4 were resistant and susceptible to F. graminearum, respectively. At four post-inoculation points, 11,070 protein species were identified, of which 762 were differentially abundant. Gene Ontology enrichment analysis showed that most differentially abundant protein species participated in 18 biological processes after inoculation. Further analysis demonstrated that crucial metabolic pathways like plant-pathogen interaction had increased abundance. Real-time quantitative PCR (qRT-PCR) analysis revealed increased gene products of eight selected genes in plant-pathogen interaction. This investigation provides a basic bioinformatics-based characterization of differentially abundant protein species during early stages against F. graminearum. SIGNIFICANCE: FHB leads to severe yield loss and reduction in grain quality in wheat and other small grain cereals. Although extensive studies have focused on wheat resistance against F. graminearum, the molecular mechanism of FHB resistance in wheat remains to be further elucidated. In the present study, Kyoto Encyclopedia of Genes and Genomes analysis indicated that ten pathways were putatively associated with FHB resistance. Principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (OPLS-DA) showed that a valuable set of differentially abundant protein species including pathogenesis-related protein species were identified for further discovery of candidate genes for FHB resistance. This investigation provides new insights into the molecular mechanisms associated with FHB resistance and as well as a foundation for future studies.

Identification of Fusarium graminearum-responsive miRNAs and their targets in wheat by sRNA sequencing and degradome analysis

Fusarium head blight (FHB), a prevalent disease of bread wheat (Triticum aestivum L.) caused by Fusarium graminearum, leads to considerable losses of yield and quality in wheat production. MicroRNAs (miRNAs) are important regulators of plant defense responses. Here, to better understand the F. graminearum-responsive miRNAs, we constructed sRNA libraries for wheat cultivar Sumai 3 challenged with F. graminearum and sterile water, respectively. As a result, a total of 203 known miRNAs from 46 families and 68 novel miRNAs were identified. Among them, 18 known and six novel miRNAs were found to be differentially expressed between the F. graminearum-infected samples and the controls and thus were considered to be responsive to F. graminearum. The expression patterns of eight miRNAs were further validated by stem-loop qRT-PCR. Meanwhile, target genes were validated by degradome sequencing. Integrative analysis of the differentially expressed miRNAs and their targets revealed complex miRNA-mediated regulatory networks involved in the response of wheat to F. graminearum infection. Our findings are expected to facilitate a better understanding of the miRNA regulation in wheat-F. graminearum interaction.

A deletion mutation in TaHRC confers Fhb1 resistance to Fusarium head blight in wheat

Fusarium head blight (FHB), which is mainly caused by Fusarium graminearum, is a destructive wheat disease that threatens global wheat production. Fhb1, a quantitative trait locus discovered in Chinese germplasm, provides the most stable and the largest effect on FHB resistance in wheat. Here we show that TaHRC, a gene that encodes a putative histidine-rich calcium-binding protein, is the key determinant of Fhb1-mediated resistance to FHB. We demonstrate that TaHRC encodes a nuclear protein conferring FHB susceptibility and that a deletion spanning the start codon of this gene results in FHB resistance. Identical sequences of the TaHRC-R allele in diverse accessions indicate that Fhb1 had a single origin, and phylogenetic and haplotype analyses suggest that the TaHRC-R allele most likely originated from a line carrying the Dahongpao haplotype. This discovery opens a new avenue to improve FHB resistance in wheat, and possibly in other cereal crops, by manipulating TaHRC sequence through bioengineering approaches.

Identification and mapping of expressed genes associated with the 2DL QTL for fusarium head blight resistance in the wheat line Wuhan 1

Background

Fusarium head blight (FHB) is a problem of great concern in small grain cereals, especially wheat. A quantitative trait locus (QTL) for FHB resistance (FHB_SFI) located on the long arm of chromosome 2D in the spring wheat genotype Wuhan 1 is a resistance locus which has potential to improve the FHB resistance of bread wheat since it confers effective resistance to wheat breeding lines. Recently, differentially expressed genes (DEG) have been identified by comparing near isogenic lines (NIL) carrying the susceptible and resistant alleles for the 2DL QTL, using RNA-Seq. In the present study, we aimed to identify candidate genes located within the genetic interval for the 2DL QTL for FHB resistance, as assessed by single floret inoculation (FHB_SFI), and possibly contributing to it.

Results

Combining previous and additional bioinformatics analyses, 26 DEG that were located on chromosome arm 2DL were selected for further characterization of their expression profile by RT-qPCR. Seven of those DEG showed a consistent differential expression profile between either three pairs of near isogenic lines or other genotypes carrying the R and S alleles for the 2DL QTL for FHB resistance. UN25696, which was identified in previous expression work using microarray was also confirmed to have a differential expression pattern. Those eight candidate genes were further characterized in 85 lines of a double haploid mapping population derived from the cross Wuhan 1/Nyubai, the population where the 2DL QTL was originally identified. The expression QTL for gene Traes_2DL_179570792 overlapped completely with the mapping interval for the 2DL QTL for FHB_SFI while the expression QTL for UN25696 mapped near the QTL, but did not overlap with it. None of the other genes had a significant eQTL on chromosome 2DL. Higher expression of Traes_2DL_179570792 and UN25696 was associated with the resistant allele at that locus.

Conclusions

Of the 26 DEG from the 2DL chromosome further characterized in this study, only two had an expression QTL located in or near the interval for the 2DL QTL. Traes_2DL_179570792 is the first expression marker identified as associated with the 2DL QTL.

Electronic supplementary material

The online version of this article (10.1186/s12863-019-0748-6) contains supplementary material, which is available to authorized users.

Applying the latest advances in genomics and phenomics for trait discovery in polyploid wheat

Improving traits in wheat has historically been challenging due to its large and polyploid genome, limited genetic diversity and in-field phenotyping constraints. However, within recent years many of these barriers have been lowered. The availability of a chromosome-level assembly of the wheat genome now facilitates a step-change in wheat genetics and provides a common platform for resources, including variation data, gene expression data and genetic markers. The development of sequenced mutant populations and gene-editing techniques now enables the rapid assessment of gene function in wheat directly. The ability to alter gene function in a targeted manner will unmask the effects of homoeolog redundancy and allow the hidden potential of this polyploid genome to be discovered. New techniques to identify and exploit the genetic diversity within wheat wild relatives now enable wheat breeders to take advantage of these additional sources of variation to address challenges facing food production. Finally, advances in phenomics have unlocked rapid screening of populations for many traits of interest both in greenhouses and in the field. Looking forwards, integrating diverse data types, including genomic, epigenetic and phenomics data, will take advantage of big data approaches including machine learning to understand trait biology in wheat in unprecedented detail.

Meta-Analysis of the QTLome of Fusarium Head Blight Resistance in Bread Wheat: Refining the Current Puzzle

Background: Fusarium Head Blight (FHB) is a worldwide devastating disease of bread wheat (Triticum aestivum L.). Genetic resistance is the most effective way to control FHB and many QTL related to this trait have been mapped on the wheat genetic map. This information, however, must be refined to be more efficiently used in breeding programs and for the advance of the basic research. The objective of the present study was to in-depth analyze the QTLome of FHB resistance in bread wheat, further integrating genetic, genomic, and transcriptomic data, aiming to find candidate genes. Methods: An exhaustive bibliographic review on 76 scientific papers was carried out collecting information about QTL related to FHB resistance mapped on bread wheat. A dense genetic consensus map with 572,862 loci was generated for QTL projection. Meta-analysis could be performed on 323 QTL. Candidate gene mining was carried out within the most refined loci, containing genes that were cross-validated with publicly available transcriptional expression data of wheat under Fusarium infection. Most highlighted genes were investigated for protein evidence. Results: A total of 556 QTL were found in the literature, distributed on all sub-genomes and chromosomes of wheat. Meta-analysis generated 65 meta-QTL, and this refinement allows one to find markers more tightly linked to these regions. Candidate gene mining within the most refined meta-QTL, meta-QTL 1/chr. 3B, harvested 324 genes and transcriptional data cross-validated 10 of these genes, as responsive to FHB. One is of these genes encodes a Glycosiltransferase and the other encodes for a Cytochrome P450, and these such proteins have already been verified as being responsible for FHB resistance, but the remaining eight genes still have to be further studied, as promising loci for breeding. Conclusions: The QTLome of FHB resistance in wheat was successfully assembled and a refinement in terms of number and length of loci was obtained. The integration of the QTLome with genomic and transcriptomic data has allowed for the discovery of promising candidate genes for use in breeding programs.

Molecular characterization of the TaWTG1 in bread wheat (Triticum aestivum L.)

OsWTG1 (LOC_Os08g42540.1) functions as an important factor determining grain size and shape in rice. Our understanding on its ortholog in wheat, TaWTG1, is limited. Here, we identified and mapped TaWTG1 in wheat, characterized its gene and protein structures, predicted transcription factor binding sites of its promoter, and the expression patterns was also analysed bases on real-time quantitative PCR and public available microarray data. The WTG1 orthologs in barley (HvWTG1), rice (OsWTG1), Aegilops tauschii (AtWTG1), Triticum urartu (TuWTG1), Triticum turgidum (TtWTG1) and Brachypodium distachyon (BdWTG1) were also identified for comparative analyses. TaWTG1 was mapped onto the short arms of group 7 chromosomes (7AS, 7BS, and 7DS). Multiple alignments indicated that WTG1 possesses eight exons and seven introns in all of the orthologs, except for the orthologs on 7A of wild emmer and on 7D of A. tauschii (seven exons and six introns). An exon-intron junction composed of intron 2 to intron 3 and exon 2 to exon 4 was highly conserved. The protein of WTG1 exists a conserved domain (Peptidase_C65). WTG1 was mainly expressed in wheat roots, spikes and grains, in barley caryopsis and roots, and in rice anthers. Drought and heat stresses significantly regulated the expression of TaWTG1 in wheat. In barley, WTG1 was significantly down-regulated under Fusarium at late stage. In addition, significant correlations between the expression patterns of predicted transcription factors and WTG1 were also detected. Overall, the results presented here broaden our knowledge on WTG1 and will be helpful for its manipulation aiming at dissecting its function in plants.

A wheat cytochrome P450 enhances both resistance to deoxynivalenol and grain yield

The mycotoxin deoxynivalenol (DON) serves as a plant disease virulence factor for the fungi Fusarium graminearum and F. culmorum during the development of Fusarium head blight (FHB) disease on wheat. A wheat cytochrome P450 gene from the subfamily CYP72A, TaCYP72A, was cloned from wheat cultivar CM82036. TaCYP72A was located on chromosome 3A with homeologs present on 3B and 3D of the wheat genome. Using gene expression studies, we showed that TaCYP72A variants were activated in wheat spikelets as an early response to F. graminearum, and this activation was in response to the mycotoxic Fusarium virulence factor deoxynivalenol (DON). Virus induced gene silencing (VIGS) studies in wheat heads revealed that this gene family contributes to DON resistance. VIGS resulted in more DON-induced discoloration of spikelets, as compared to mock VIGS treatment. In addition to positively affecting DON resistance, TaCYP72A also had a positive effect on grain number. VIGS of TaCYP72A genes reduced grain number by more than 59%. Thus, we provide evidence that TaCYP72A contributes to host resistance to DON and conclude that this gene family warrants further assessment as positive contributors to both biotic stress resistance and grain development in wheat.

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.

A multiple near isogenic line (multi-NIL) RNA-seq approach to identify candidate genes underpinning QTL

This study demonstrates how identification of genes underpinning disease-resistance QTL based on differential expression and SNPs can be improved by performing transcriptomic analysis on multiple near isogenic lines. Transcriptomic analysis has been widely used to understand the genetic basis of a trait of interest by comparing genotypes with contrasting phenotypes. However, these approaches identify such large sets of differentially expressed genes that it proves difficult to isolate which genes underpin the phenotype of interest. This study tests whether using multiple near isogenic lines (NILs) can improve the resolution of RNA-seq-based approaches to identify genes underpinning disease-resistance QTL. A set of NILs for a major effect Fusarium crown rot-resistance QTL in barley on the 4HL chromosome arm were analysed under Fusarium crown rot using RNA-seq. Differential gene expression and single nucleotide polymorphism detection analyses reduced the number of putative candidates from thousands within individual NIL pairs to only one hundred and two genes, which were differentially expressed or contained SNPs in common across NIL pairs and occurred on 4HL. Our findings support the value of performing RNA-seq analysis using multiple NILs to remove genetic background effects. The enrichment analyses indicated conserved differences in the response to infection between resistant and sensitive isolines suggesting that sensitive isolines are impaired in systemic defence response to Fusarium pseudograminearum.

A Pathogen-Responsive Leucine Rich Receptor Like Kinase Contributes to Fusarium Resistance in Cereals

Receptor-like kinases form the largest family of receptors in plants and play an important role in recognizing pathogen-associated molecular patterns and modulating the plant immune responses to invasive fungi, including cereal defenses against fungal diseases. But hitherto, none have been shown to modulate the wheat response to the economically important Fusarium head blight (FHB) disease of small-grain cereals. Homologous genes were identified on barley chromosome 6H (HvLRRK-6H) and wheat chromosome 6DL (TaLRRK-6D), which encode the characteristic domains of surface-localized receptor like kinases. Gene expression studies validated that the wheat TaLRRK-6D is highly induced in heads as an early response to both the causal pathogen of FHB disease, Fusarium graminearum, and its' mycotoxic virulence factor deoxynivalenol. The transcription of other wheat homeologs of this gene, located on chromosomes 6A and 6B, was also up-regulated in response to F. graminearum. Virus-induced gene silencing (VIGS) of the barley HvLRRK-6H compromised leaf defense against F. graminearum. VIGS of TaLRRK-6D in two wheat cultivars, CM82036 (resistant to FHB disease) and cv. Remus (susceptible to FHB), confirmed that TaLRRK-6D contributes to basal resistance to FHB disease in both genotypes. Although the effect of VIGS did not generally reduce grain losses due to FHB, this experiment did reveal that TaLRRK-6D positively contributes to grain development. Further gene expression studies in wheat cv. Remus indicated that VIGS of TaLRRK-6D suppressed the expression of genes involved in salicylic acid signaling, which is a key hormonal pathway involved in defense. Thus, this study provides the first evidence of receptor like kinases as an important component of cereal defense against Fusarium and highlights this gene as a target for enhancing cereal resistance to FHB disease.

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.

Time-course expression QTL-atlas of the global transcriptional response of wheat to Fusarium graminearum

Fusarium head blight is a devastating disease of small grain cereals such as bread wheat (Triticum aestivum). The pathogen switches from a biotrophic to a nectrotrophic lifestyle in course of disease development forcing its host to adapt its defence strategies. Using a genetical genomics approach, we illustrate genome-wide reconfigurations of genetic control over transcript abundances between two decisive time points after inoculation with the causative pathogen Fusarium graminearum. Whole transcriptome measurements have been recorded for 163 lines of a wheat doubled haploid population segregating for several resistance genes yielding 15 552 at 30 h and 15 888 eQTL at 50 h after inoculation. The genetic map saturated with transcript abundance-derived markers identified of a novel QTL on chromosome 6A, besides the previously reported QTL Fhb1 and Qfhs.ifa-5A. We find a highly different distribution of eQTL between time points with about 40% of eQTL being unique for the respective assessed time points. But also for more than 20% of genes governed by eQTL at either time point, genetic control changes in time. These changes are reflected in the dynamic compositions of three major regulatory hotspots on chromosomes 2B, 4A and 5A. In particular, control of defence-related biological mechanisms concentrated in the hotspot at 4A shift to hotspot 2B as the disease progresses. Hotspots do not colocalize with phenotypic QTL, and within their intervals no higher than expected number of eQTL was detected. Thus, resistance conferred by either QTL is mediated by few or single genes.