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

Genome-wide identification and examination of the wheat glycosyltransferase family 43 regulation during Fusarium graminearum infection

In Arabidopsis and rice, the glycosyltransferase (GT) 43 family is involved in xylan synthesis. However, there have been limited reports on the study of the TaGT43 family in wheat. In this study, 28 TaGT43 family members were identified in wheat (Triticum aestivum L.) and clustered into three major groups based on the similarity of amino acid sequences. The results of the TaGT43 family's conserved motif and gene structure analyses agree with this result. Collinearity analysis revealed segmental duplications mainly promoted TaGT43 family expansion. cis-Acting element analysis revealed that the TaGT43 genes were involved in the light response, phytohormone response, abiotic/biotic stress response, and growth and development. The TaGT43 family showed a tissue-specific expression pattern, primarily expressed in roots and stems. Besides, the transcriptional and expression levels of multiple TaGT43 genes were upregulated during the infection of F. graminearum. According to metabolomics studies, F. graminearum infection affected the phenylpropanoid biosynthesis pathway in wheat, a critical factor in cell wall construction. Furthermore, GO enrichment analysis indicated that the TaGT43 genes play a significant role in cell wall organization. Based on these results, it may be concluded that the TaGT43 family mediates cell wall organization in response to F. graminearum infection.

Chemical-tag-based semi-annotated metabolomics facilitates gene identification and specialized metabolic pathway elucidation in wheat

The importance of metabolite modification and species-specific metabolic pathways has long been recognized. However, linking the chemical structure of metabolites to gene function in order to explore the genetic and biochemical basis of metabolism has not yet been reported in wheat (Triticum aestivum). Here, we profiled metabolic fragment enrichment in wheat leaves and consequently applied chemical-tag-based semi-annotated metabolomics in a genome-wide association study in accessions of wheat. The studies revealed that all 1,483 quantified metabolites have at least one known functional group whose modification is tailored in an enzyme-catalyzed manner and eventually allows efficient candidate gene mining. A Triticeae crop-specific flavonoid pathway and its underlying metabolic gene cluster were elucidated in further functional studies. Additionally, upon overexpressing the major effect gene of the cluster TraesCS2B01G460000 (TaOMT24), the pathway was reconstructed in rice (Oryza sativa), which lacks this pathway. The reported workflow represents an efficient and unbiased approach for gene mining using forward genetics in hexaploid wheat. The resultant candidate gene list contains vast molecular resources for decoding the genetic architecture of complex traits and identifying valuable breeding targets and will ultimately aid in achieving wheat crop improvement.

New insights into molecular features of the genome-wide AOX family and their responses to various stresses in common wheat (Triticum aestivum L.)

Alternative oxidase (AOX) is an important terminal oxidase involved in the alternative oxidation pathway in plants, which is closely related to various biotic and abiotic stress responses. However, a comprehensive research on AOX gene family of wheat is still lacking. In this study, the members of wheat AOX (TaAOX) family were identified, and their molecular characteristics and gene expression patterns were systematically investigated. Seventeen TaAOX genes were identified from Chinese Spring (CS) genome, which were mapped on 7 chromosomes and mainly clustered on the long arm's distal end of the second homologous groups. Phylogenetic analysis showed that TaAOX genes were classified into four subgroups (Ia, Ib, Ic, and Id), and the Ia subgroup possessed the most members. Tandem duplication and segmental duplication events were found during the evolution of TaAOX genes and they were affected by purifying selection demonstrated by Ka/Ks analysis. The exon numbers of this family gene varied greatly from 1 to 9. Except for Ta3BSAOX14, all the proteins encoded by the other 16 TaAOX genes contained the amino acid residues of the key active sites in the AOX domain (cd01053). The expression patterns of TaAOX genes in various tissues and under abiotic and biotic stresses were analyzed using public transcriptome data, furthermore, qRT-PCR analysis was performed for some selected TaAOX genes, and the results suggested that most members of this gene family play an important role in response to different stresses in common wheat. Our results provide basic information and valuable reference for further exploring the gene function of TaAOX family by using gene editing, RNAi, VIGS, and other technologies.

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.

Expression Profiling of the Slow Rusting Resistance Genes Lr34/Yr18 and Lr67/Yr46 in Common Wheat (Triticum aestivum L.) and Associated miRNAs Patterns

The main efforts in common wheat (Triticum aestivum L.) breeding focus on yield, grain quality, and resistance to biotic and abiotic stresses. One of the major threats affecting global wheat cultivation and causing significant crop production losses are rust diseases, including leaf rust caused by a biotrophic fungus Puccinia triticina Eriks. Genetically determined resistance to leaf rust has been characterized in young plants (seedling resistance) as well as in plants at the adult plant stage. At the seedling stage, resistance is controlled vertically by major R genes, conferring a race-specific response that is highly effective but usually short-lived due to the rapid evolution of potentially virulent fungi. In mature plants, horizontal adult plant resistance (APR) was described, which provides long-term protection against multiple races of pathogens. A better understanding of molecular mechanisms underlying the function of APR genes would enable the development of new strategies for resistance breeding in wheat. Therefore, in the present study we focused on early transcriptomic responses of two major wheat APR genes, Lr34 and Lr67, and three complementary miRNAs, tae-miR9653b, tae-miR9773 and tae-miR9677b, to inoculation with P. triticina. Plant material consisted of five wheat reference varieties, Artigas, NP846, Glenlea, Lerma Rojo and TX89D6435, containing the Lr34/Yr18 and Lr67/Yr46 resistance genes. Biotic stress was induced by inoculation with fungal spores under controlled conditions in a phytotron. Plant material consisted of leaf tissue sampled before inoculation as well as 6, 12, 24 and 48 h postinoculation (hpi). The APR gene expression was quantified using real-time PCR with two reference genes, whereas miRNA was quantified using droplet digital PCR. This paper describes the resistance response of APR genes to inoculation with races of leaf rust-causing fungi that occur in central Europe. The study revealed high variability of expression profiles between varieties and time-points, with the prevalence of downregulation for APR genes and upregulation for miRNAs during the development of an early defense response. Nevertheless, despite the downregulation initially observed, the expression of Lr34 and Lr67 genes in studied cultivars was significantly higher than in a control line carrying wild (susceptible) alleles.

Genome-wide characterization of the VQ genes in Triticeae and their functionalization driven by polyploidization and gene duplication events in wheat

Valine-glutamine motif-containing (VQ) proteins are transcriptional cofactors widely involved in plant growth, development, and response to various stresses. Although the VQ family has been genome-wide identified in some species, but the knowledge regarding duplication-driven functionalization of VQ genes among evolutionarily related species is still lacking. Here, 952 VQ genes have been identified from 16 species, emphasizing seven Triticeae species including the bread wheat. Comprehensive phylogenetic and syntenic analyses allow us to establish the orthologous relationship of VQ genes from rice (Oryza sativa) to bread wheat (Triticum aestivum). The evolutionary analysis revealed that whole-genome duplication (WGD) drives the expansion of OsVQs, while TaVQs expansion is associated with a recent burst of gene duplication (RBGD). We also analyzed the motif composition and molecular properties of TaVQ proteins, enriched biological functions, and expression patterns of TaVQs. We demonstrate that WGD-derived TaVQs have become divergent in both protein motif composition and expression pattern, while RBGD-derived TaVQs tend to adopt specific expression patterns, suggesting their functionalization in certain biological processes or in response to specific stresses. Furthermore, some RBGD-derived TaVQs are found to be associated with salt tolerance. Several of the identified salt-related TaVQ proteins were located in the cytoplasm and nucleus and their salt-responsive expression patterns were validated by qPCR analysis. Yeast-based functional experiments confirmed that TaVQ27 may be a new regulator to salt response and regulation. Overall, this study lays the foundation for further functional validation of VQ family members within the Triticeae species.

Comparative Transcriptomics of Fusarium graminearum and Magnaporthe oryzae Spore Germination Leading up To Infection

For fungal plant pathogens, the germinating spore provides the first interaction with the host. Spore germlings move across the plant surface and use diverse penetration strategies for ingress into plant surfaces. Penetration strategies include pressurized melanized appressoria, which facilitate physically punching through the plant cuticle, and nonmelanized appressoria, which penetrate with the help of enzymes or cuticular damage to breach the plant surface. Two well-studied plant pathogens, Fusarium graminearum and Magnaporthe oryzae, are typical of these two modes of penetration. We applied comparative transcriptomics to Fusarium graminearum and Magnaporthe oryzae to characterize the genetic programming of the early host-pathogen interface. Four sequential stages of development following spore localization on the plant surface, from spore swelling to appressorium formation, were sampled for each species on culture medium and on barley sheaths, and transcriptomic analyses were performed. Gene expression in the prepenetration stages in both species and under both conditions was similar. In contrast, gene expression in the final stage was strongly influenced by the environment. Appressorium formation involved the greatest number of differentially expressed genes. Laser-dissection microscopy was used to perform detailed transcriptomics of initial infection points by F. graminearum. These analyses revealed new and important aspects of early fungal ingress in this species. Expression of the trichothecene genes involved in biosynthesis of deoxynivalenol by F. graminearum implies that toxisomes are not fully functional until after penetration and indicates that deoxynivalenol is not essential for penetration under our conditions. The use of comparative gene expression of divergent fungi promises to advance highly effective targets for antifungal strategies. IMPORTANCE Fusarium graminearum and Magnaporthe oryzae are two of the most important pathogens of cereal grains worldwide. Despite years of research, strong host resistance has not been identified for F. graminearum, so other methods of control are essential. The pathogen takes advantage of multiple entry points to infect the host, including breaches in the florets due to senescence of flower parts and penetration of the weakened trichome bases to breach the epidermis. In contrast, M. oryzae directly punctures leaves that it infects, and resistant cultivars have been characterized. The threat of either pathogen causing a major disease outbreak is ever present. Comparative transcriptomics demonstrated its potential to reveal novel and effective disease prevention strategies that affect the initial stages of disease. Shedding light on the basis of this diversity of infection strategies will result in development of increasingly specific control strategies.

Transcriptome Dynamic Analysis Reveals New Candidate Genes Associated with Resistance to Fusarium Head Blight in Two Chinese Contrasting Wheat Genotypes

In recent years, Fusarium head blight (FHB) has developed into a global disease that seriously affects the yield and quality of wheat. Effective measures to solve this problem include exploring disease-resistant genes and breeding disease-resistant varieties. In this study, we conducted a comparative transcriptome analysis to identify the important genes that are differentially expressed in FHB medium-resistant (Nankang 1) and FHB medium-susceptible (Shannong 102) wheat varieties for various periods after Fusarium graminearum infection using RNA-seq technology. In total, 96,628 differentially expressed genes (DEGs) were identified, 42,767 from Shannong 102 and 53,861 from Nankang 1 (FDR < 0.05 and |log2FC| > 1). Of these, 5754 and 6841 genes were found to be shared among the three time points in Shannong 102 and Nankang 1, respectively. After inoculation for 48 h, the number of upregulated genes in Nankang 1 was significantly lower than that of Shannong 102, but at 96 h, the number of DEGs in Nankang 1 was higher than that in Shannong 102. This indicated that Shannong 102 and Nankang 1 had different defensive responses to F. graminearum in the early stages of infection. By comparing the DEGs, there were 2282 genes shared at the three time points between the two strains. GO and KEGG analyses of these DEGs showed that the following pathways were associated with disease resistance genes: response to stimulus pathway in GO, glutathione metabolism, phenylpropanoid biosynthesis, plant hormone signal transduction, and plant-pathogen interaction in KEGG. Among them, 16 upregulated genes were identified in the plant-pathogen interaction pathway. There were five upregulated genes, TraesCS5A02G439700, TraesCS5B02G442900, TraesCS5B02G443300, TraesCS5B02G443400, and TraesCS5D02G446900, with significantly higher expression levels in Nankang 1 than in Shannong 102, and these genes may have an important role in regulating the resistance of Nankang 1 to F. graminearum infection. The PR proteins they encode are PR protein 1-9, PR protein 1-6, PR protein 1-7, PR protein 1-7, and PR protein 1-like. In addition, the number of DEGs in Nankang 1 was higher than that in Shannong 102 on almost all chromosomes, except chromosomes 1A and 3D, but especially on chromosomes 6B, 4B, 3B, and 5A. These results indicate that gene expression and the genetic background must be considered for FHB resistance in wheat breeding.

'Garlic-lipo'4Plants: Liposome-Encapsulated Garlic Extract Stimulates ABA Pathway and PR Genes in Wheat (Triticum aestivum)

Recently, environmentally friendly crop improvements using next-generation plant biostimulants (PBs) come to the forefront in agriculture, regardless of whether they are used by scientists, farmers, or industries. Various organic and inorganic solutions have been investigated by researchers and producers, focusing on tolerance to abiotic and biotic stresses, crop quality, or nutritional deficiency. Garlic has been considered a universal remedy ever since antiquity. A supercritical carbon dioxide garlic extract encapsulated in nanoscale liposomes composed of plant-derived lipids was examined as a possible PB agent. The present study focused on the characterization of the genes associated with the pathways involved in defense response triggered by the liposome nanoparticles that were loaded with supercritical garlic extracts. This material was applied to Triticum aestivum in greenhouse experiments using foliar spraying. The effects were examined in a large-scale genome-wide transcriptional profiling experiment by collecting the samples four times (0 min, used as a control, and 15 min, 24 h, and 48 h after spraying). Based on a time-course expression analysis, the dynamics of the cellular response were determined by examining differentially expressed genes and applying a cluster analysis. The results suggested an enhanced expression of abscisic acid (ABA) pathway and pathogenesis-related (PR) genes, of which positive regulation was found for the AP2-, C2H2-, HD-ZIP-, and MYB-related transcription factor families.

Time Course Metabolite Profiling of Fusarium Head Blight-Infected Hard Red Spring Wheat Using Ultra-High-Performance Liquid Chromatography Coupled with Quadrupole Time of Flight/MS

Wheat is an important food crop, yet its value is reduced by fungal infections (ex. Fusarium graminearum). Metabolite profiling is a useful tool for explaining resistance mechanisms. By analyzing near-isogenic lines (NILs) with contrasting Fhb1 alleles and three wheat varieties, a time course resulting in 61 relevant metabolites was studied. The presence of one metabolite as resistant related constitutive late in the time course was detected. Results confirm the presence of hydroxycinnamic acid amides conjugated with polyamine derivatives (hydroxycinnamic acid amides, HCAAs), which have been shown to induce thickening of cell walls. These compounds are shared by resistant and susceptible genotypes with no difference in intensities but vary in time as early- or late-occurring, suggesting that for the NIL studied here, HCAAs were a normal part of the host reaction. Overall, metabolites synthesized as a result of infection were observed regardless of susceptibility but occurred at different times after infection.

Genomic sequencing of Thinopyrum elongatum chromosome arm 7EL, carrying fusarium head blight resistance, and characterization of its impact on the transcriptome of the introgressed line CS-7EL

Background

The tall wheatgrass species Thinopyrum elongatum carries a strong fusarium head blight (FHB) resistance locus located on the long arm of chromosome 7 (7EL) as well as resistance to leaf and stem rusts, all diseases with a significant impact on wheat production. Towards understanding the contribution of Th. elongatum 7EL to improvement of disease resistance in wheat, the genomic sequence of the 7EL fragment present in the wheat Chinese Spring (CS) telosomic addition line CS-7EL was determined and the contribution and impact of 7EL on the rachis transcriptome during FHB infection was compared between CS and CS-7EL.

Results

We assembled the Th. elongatum 7EL chromosome arm using a reference-guided approach. Combining this assembly with the available reference sequence for CS hexaploid wheat provided a reliable reference for interrogating the transcriptomic differences in response to infection conferred by the 7EL fragment. Comparison of the transcriptomes of rachis tissues from CS and CS-7EL showed expression of Th. elongatum transcripts as well as modulation of wheat transcript expression profiles in the CS-7EL line. Expression profiles at 4 days after infection with Fusarium graminearum, the causal agent of FHB, showed an increased in expression of genes associated with an effective defense response, in particular glucan endo-1,3-beta-glucosidases and chitinases, in the FHB-resistant line CS-7EL while there was a larger increase in differential expression for genes associated with the level of fungal infection in the FHB-susceptible line CS. One hundred and seven 7EL transcripts were expressed in the smallest 7EL region defined to carry FHB resistance.

Conclusion

7EL contributed to CS-7EL transcriptome by direct expression and through alteration of wheat transcript profiles. FHB resistance in CS-7EL was associated with transcriptome changes suggesting a more effective defense response. A list of candidate genes for the FHB resistance locus on 7EL has been established.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08433-8.

Identification of Wheat LACCASEs in Response to Fusarium graminearum as Potential Deoxynivalenol Trappers

Fusarium graminearum (F. graminearum) can cause huge yield reductions and contamination of grain with deoxynivalenol (DON), and thus is one of the most problematic pathogen of wheat worldwide. Although great efforts have been paid and great achievements have been made to control the pathogens, there is still a wide gap for understanding the mechanism underlying F. graminearum resistance. Plant LACCASEs (LACs) catalyze the oxidative polymerization of monolignols by reinforcing cell-wall of various cell types to provide mechanical support, xylem sap transportation, and defense against pest and pathogens. To date, little has been known about LAC genes in bread wheat and their potential roles in wheat-F. graminearum interaction. Through systematic analysis of the genome-wide homologs and transcriptomes of wheat, a total of 95 Triticum aestivum laccases (TaLACs) were identified, and 14 of them were responsive to F. graminearum challenge. 3D structure modelings of the 14 TaLAC proteins showed that only TaLAC78 contains the entire activity center for oxidation and the others lack the type 1 copper ion ligand (T1Cu). Both amino acid sequence alignment and three-dimensional reconstruction after amino acid mutation showed that the loss of T1Cu is not only related to variation of the key amino acid coordinating T1Cu, but also closely related to the flanking amino acids. Significantly differential temporal expression patterns of TaLACs suggested that their subfunctionalization might occur. Promoter array analysis indicated that the induction of TaLACs may be closely associated with salicylic acid signaling, dehydration, and low-oxygen stress under F. graminearum infection. Molecular docking simulation demonstrated that TaLACs can not only catalyze lignin as a substrate, but also interact with DON, which may be docked into the binding position of the monolignols, where the LACs recognize substrates. The current study provides clues for exploring the novel functions of TaLACs in wheat resistance to F. graminearum, and TaLACs maybe candidates for conferring a high level of resistance against F. graminearum in wheat.

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

Background

Treatment of wheat with the phytohormones abscisic acid (ABA) and gibberellic acid (GA) has been shown to affect Fusarium head blight (FHB) disease severity. However, the molecular mechanisms underlying the elicited phenotypes remain unclear. Toward addressing this gap in our knowledge, global transcriptomic profiling was applied to the FHB-susceptible wheat cultivar ‘Fielder’ to map the regulatory responses effected upon treatment with ABA, an ABA receptor antagonist (AS6), or GA in the presence or absence of Fusarium graminearum (Fg) challenge.

Results

Spike treatments resulted in a total of 30,876 differentially expressed genes (DEGs) identified in ‘Fielder’ (26,004) and the Fg (4872) pathogen. Topology overlap and correlation analyses defined 9689 wheat DEGs as Fg-related across the treatments. Further enrichment analyses demonstrated that these included expression changes within ‘Fielder’ defense responses, cell structural metabolism, molecular transport, and membrane/lipid metabolism. Dysregulation of ABA and GA crosstalk arising from repression of ‘Fielder’ FUS3 was noted. As well, expression of a putative Fg ABA-biosynthetic cytochrome P450 was detected. The co-applied condition of Fg + ABA elicited further up-regulation of phytohormone biosynthesis, as well as SA and ET signaling pathways and cell wall/polyphenolic metabolism. In contrast, co-applied Fg + GA mainly suppressed phytohormone biosynthesis and signaling, while modulating primary and secondary metabolism and flowering. Unexpectedly, co-applied Fg + AS6 did not affect ABA biosynthesis or signaling, but rather elicited antagonistic responses tied to stress, phytohormone transport, and FHB disease-related genes.

Conclusions

Observed exacerbation (misregulation) of classical defense mechanisms and cell wall fortifications upon ABA treatment are consistent with its ability to promote FHB severity and its proposed role as a fungal effector. In contrast, GA was found to modulate primary and secondary metabolism, suggesting a general metabolic shift underlying its reduction in FHB severity. While AS6 did not antagonize traditional ABA pathways, its impact on host defense and Fg responses imply potential for future investigation. Overall, by comparing these findings to those previously reported for four additional plant genotypes, an additive model of the wheat-Fg interaction is proposed in the context of phytohormone responses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08069-0.

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.

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

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

De novo transcriptome assembly, functional annotation, and expression profiling of rye (Secale cereale L.) hybrids inoculated with ergot (Claviceps purpurea)

Rye is used as food, feed, and for bioenergy production and remain an essential grain crop for cool temperate zones in marginal soils. Ergot is known to cause severe problems in cross-pollinated rye by contamination of harvested grains. The molecular response of the underlying mechanisms of this disease is still poorly understood due to the complex infection pattern. RNA sequencing can provide astonishing details about the transcriptional landscape, hence we employed a transcriptomic approach to identify genes in the underlying mechanism of ergot infection in rye. In this study, we generated de novo assemblies from twelve biological samples of two rye hybrids with identified contrasting phenotypic responses to ergot infection. The final transcriptome of ergot susceptible (DH372) and moderately ergot resistant (Helltop) hybrids contain 208,690 and 192,116 contigs, respectively. By applying the BUSCO pipeline, we confirmed that these transcriptome assemblies contain more than 90% of gene representation of the available orthologue groups at Virdiplantae odb10. We employed a de novo assembled and the draft reference genome of rye to count the differentially expressed genes (DEGs) between the two hybrids with and without inoculation. The gene expression comparisons revealed that 228 genes were linked to ergot infection in both hybrids. The genome ontology enrichment analysis of DEGs associated them with metabolic processes, hydrolase activity, pectinesterase activity, cell wall modification, pollen development and pollen wall assembly. In addition, gene set enrichment analysis of DEGs linked them to cell wall modification and pectinesterase activity. These results suggest that a combination of different pathways, particularly cell wall modification and pectinesterase activity contribute to the underlying mechanism that might lead to resistance against ergot in rye. Our results may pave the way to select genetic material to improve resistance against ergot through better understanding of the mechanism of ergot infection at molecular level. Furthermore, the sequence data and de novo assemblies are valuable as scientific resources for future studies in rye.

Infection cushions of Fusarium graminearum are fungal arsenals for wheat infection

Fusarium graminearum is one of the most destructive plant pathogens worldwide, causing fusarium head blight (FHB) on cereals. F. graminearum colonizes wheat plant surfaces with specialized unbranched hyphae called runner hyphae (RH), which develop multicelled complex appressoria called infection cushions (IC). IC generate multiple penetration sites, allowing the fungus to enter the plant cuticle. Complex infection structures are typical for several economically important plant pathogens, yet with unknown molecular basis. In this study, RH and IC formed on the surface of wheat paleae were isolated by laser capture microdissection. RNA-Seq-based transcriptomic analyses were performed on RH and IC and compared to mycelium grown in complete medium (MY). Both RH and IC displayed a high number of infection up-regulated genes (982), encoding, among others, carbohydrate-active enzymes (CAZymes: 140), putative effectors (PE: 88), or secondary metabolism gene clusters (SMC: 12 of 67 clusters). RH specifically up-regulated one SMC corresponding to aurofusarin biosynthesis, a broad activity antibiotic. IC specifically up-regulated 248 genes encoding mostly putative virulence factors such as 7 SMC, including the mycotoxin deoxynivalenol and the newly identified fusaoctaxin A, 33 PE, and 42 CAZymes. Furthermore, we studied selected candidate virulence factors using cellular biology and reverse genetics. Hence, our results demonstrate that IC accumulate an arsenal of proven and putative virulence factors to facilitate the invasion of epidermal cells.

Transcriptome analysis of near-isogenic lines for glume hairiness of wheat

Hairiness, which is a phenotypic trait common among land plants, primarily affects the stem, leaf, and floral organs. Plant hairiness is associated with complex functions. For example, glume hairiness in wheat is related to the resistance to biotic and abiotic stresses, and may also influence human health. In the present study, two pairs of near-isogenic lines (NILs) for glume hairiness, which were derived from a cross between a Tibetan semi-wild wheat accession (Triticum aestivum ssp. tibetanum Q1028) and a common wheat cultivar (T. aestivum 'Zhengmai 9023'), underwent a glume transcriptome analysis. We detected 27,935 novel genes, of which 18,027 were annotated. Additionally, 488 and 600 differentially expressed genes (DEGs) were detected in NIL1 and NIL2, respectively, with 37 DEGs detected in both NIL pairs. Moreover, 987 and 1584 single nucleotide polymorphisms (SNPs) were detected in NIL1 and NIL2, respectively, with 39 SNPs detected in both NIL pairs, of which most were located in the Hairy glume (Hg) gene region on chromosome arm 1AS. The annotation of the DEGs with gene ontology terms revealed that genes associated with hairiness in Arabidopsis and rice were similarly enriched. The possible functions of these genes related to glume hairiness were examined. The study results provide useful information for identifying candidate genes and the fine-mapping of Hg in the wheat genome.

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.

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

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

Transcription Factor FOXO3a Is a Negative Regulator of Cytotoxicity of Fusarium mycotoxin in GES-1 Cells

Molecular mechanism and key factors responsible for cytotoxicity against mycotoxin deoxynivalenol (DON) from Fusarium pathogens are rarely elucidated. In this study, rapid increases of ROS were first observed in human gastric epithelial (GES-1) cells under DON exposure. Mitochondrial DNA damage, impaired respiratory chain, and decreased oxygen consumption rate (OCR) values, as well as G2/M cell cycle arrest and apoptosis, were also detected. Via combinatorial approaches of a large-scale microarray of differentially expressed genes, high content and RNAi analysis, a transcription factor of Forkhead box O3 (FOXO3a) was found with crucial functionalities, regulated some apoptotic genes associated with mitochondrial toxicity and cell death after activation by nuclear translocation. Namely, knockdown of FOXO3a decreased the cytotoxicity of DON to GES-1 cells. Moreover, knockdown of the FOXO ortholog DAF16 in Caenorhabditis elegans increased the resistance to DON-induced cytotoxicity. Simultaneously, the signaling pathway of ROS/JNK/FOXO3a of DON-induced cytotoxicity was newly proposed. In total, FOXO3a via ROS/JNK/FOXO3a plays a critical role to function as negative regulator associating with DON-induced cytotoxicity, with the potential extending to other substances.

Trichothecenes in Cereal Grains - An Update

Trichothecenes are sesquiterpenoid mycotoxins produced by fungi from the order Hypocreales, including members of the Fusarium genus that infect cereal grain crops. Different trichothecene-producing Fusarium species and strains have different trichothecene chemotypes belonging to the Type A and B class. These fungi cause a disease of small grain cereals, called Fusarium head blight, and their toxins contaminate host tissues. As potent inhibitors of eukaryotic protein synthesis, trichothecenes pose a health risk to human and animal consumers of infected cereal grains. In 2009, Foroud and Eudes published a review of trichothecenes in cereal grains for human consumption. As an update to this review, the work herein provides a comprehensive and multi-disciplinary review of the Fusarium trichothecenes covering topics in chemistry and biochemistry, pathogen biology, trichothecene toxicity, molecular mechanisms of resistance or detoxification, genetics of resistance and breeding strategies to reduce their contamination of wheat and barley.

Stable Isotope-Assisted Plant Metabolomics: Investigation of Phenylalanine-Related Metabolic Response in Wheat Upon Treatment With the Fusarium Virulence Factor Deoxynivalenol

The major Fusarium mycotoxin deoxynivalenol (DON) is a virulence factor in wheat and has also been shown to induce defense responses in host plant tissue. In this study, global and tracer labeling with 13C were combined to annotate the overall metabolome of wheat spikes and to evaluate the response of phenylalanine-related pathways upon treatment with DON. At anthesis, spikes of resistant and susceptible cultivars as well as two related near isogenic wheat lines (NILs) differing in the presence/absence of the major resistance QTL Fhb1 were treated with 1 mg DON or water (control), and samples were collected at 0, 12, 24, 48, and 96 h after treatment (hat). A total of 172 Phe-derived wheat constituents were detected with our untargeted approach employing 13C-labeled phenylalanine and subsequently annotated as flavonoids, lignans, coumarins, benzoic acid derivatives, hydroxycinnamic acid amides (HCAAs), as well as peptides. Ninety-six hours after the DON treatment, up to 30% of the metabolites biosynthesized from Phe showed significantly increased levels compared to the control samples. Major metabolic changes included the formation of precursors of compounds implicated in cell wall reinforcement and presumed antifungal compounds. In addition, also dipeptides, which presumably are products of proteolytic degradation of truncated proteins generated in the presence of the toxin, were significantly more abundant upon DON treatment. An in-depth comparison of the two NILs with correlation clustering of time course profiles revealed some 70 DON-responsive Phe derivatives. While several flavonoids had constitutively different abundance levels between the two NILs differing in resistance, other Phe-derived metabolites such as HCAAs and hydroxycinnamoyl quinates were affected differently in the two NILs after treatment with DON. Our results suggest a strong activation of the general phenylpropanoid pathway and that coumaroyl-CoA is mainly diverted towards HCAAs in the presence of Fhb1, whereas the metabolic route to monolignol(-conjugates), lignans, and lignin seems to be favored in the absence of the Fhb1 resistance quantitative trait loci.

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.

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.

Comparative transcriptome profiling of the response to Pyrenochaeta lycopersici in resistant tomato cultivar Mogeor and its background genotype-susceptible Moneymaker

Breeding for resistance is the most effective tool for controlling the corky root disease of tomato caused by Pyrenochaeta lycopersici. A comparative RNA-Seq-based transcriptomic analysis was conducted at 96 hpi (hours post infection) on two tomato cultivars: resistant Mogeor and its genetic background, and susceptible Moneymaker to investigate the differences in their transcriptomic response and identify the molecular bases of this plant-pathogen interaction. The number of differentially expressed genes (DEGs) identified was much higher in the susceptible than in the resistant genotype; however, the proportion of upregulated genes was higher in Mogeor (70.81%) than in Moneymaker (52.95%). Gene Ontology (GO) analysis enabled identification of 24 terms shared by the two cultivars that were consistent with responses to external stimulus, such as fungal infection. On the other hand, as many as 54 GO were enriched solely in Moneymaker, including terms related to defense response and cell wall metabolism. Our results could support the previous observations in other pathosystems, that susceptibility and resistance have overlapping signaling pathways and responses, suggesting that the P. lycopersici resistance gene pyl might be a recessive allele at a susceptibility locus, for which different candidate genes were identified based on the differences in induction or expression levels, observed between the resistant and susceptible genotype. MapMan analysis highlighted a complex hormone and transcription factors interplay where SA- and JA-induced pathways are modulated in a similar way in both genotypes and thus take part in a common response while the ethylene signaling pathways, induced mainly in susceptible Moneymaker, seem putatively contribute to its susceptibility.

A Comparison of sun, ovate, fs8.1 and Auxin Application on Tomato Fruit Shape and Gene Expression

Elongated tomato fruit shape is the result of the action of the fruit shape genes possibly in coordination with the phytohormone auxin. To investigate the possible link between auxin and the fruit shape genes, a series of auxin (2,4-D) treatments were performed on the wild-type and the fruit shape near-isogenic lines (NILs) in Solanum pimpinellifolium accession LA1589 background. Morphological and histological analyses indicated that auxin application approximately 3 weeks before anthesis led to elongated pear-shaped ovaries and fruits, which was mainly attributed to the increase of ovary/fruit proximal end caused by the increase of both cell number and cell size. Fruit shape changes caused by SUN, OVATE and fs8.1 were primarily due to the alterations of cell number along different growth axes. Particularly, SUN caused elongation by extending cell number along the entire proximal-distal axis, whereas OVATE caused fruit elongation in the proximal area, which was most similar to the effect of auxin on ovary shape. Expression analysis of flower buds at different stages in fruit shape NILs indicated that SUN had a stronger impact on the transcriptome than OVATE and fs8.1. The sun NIL differentially expressed genes were enriched in several biological processes, such as lipid metabolism, ion transmembrane and actin cytoskeleton organization. Additionally, SUN also shifted the expression of the auxin-related genes, including those involved in auxin biosynthesis, homeostasis, signal transduction and polar transport, indicating that SUN may regulate ovary/fruit shape through modifying the expression of auxin-related genes very early during the formation of the ovary in the developing flower.

Proteomics of Bulked Rachides Combined with Documented QTL Uncovers Genotype Nonspecific Players of the Fusarium Head Blight Responses in Wheat

Fusarium head blight (FHB) is a destructive disease of wheat that reduces yield and grain quality. High-throughput proteomic techniques have been used to identify a wide range of candidate proteins involved in host resistance. The majority of the published works on the proteomics of the wheat response to Fusarium graminearum infection are case specific. In the current study, a high-throughput quantitative label-free strategy was employed on bulked rachides of F. graminearum-infected wheat collected from multiple genotypes. Differentially accumulated proteins among the following four pools were identified: mock-inoculated FHB-resistant accessions (RM), mock-inoculated FHB-susceptible accessions (SM), F. graminearum-inoculated FHB-resistant accessions (RFg), and F. graminearum-inoculated FHB-susceptible accessions (SFg). Four pairs of comparisons were made: RFg versus RM, SFg versus SM, RM versus SM, and RFg versus SFg. Proteins were projected onto the consensus intervals of previously reported quantitative trait loci in the FHB-resistant pool by blasting against the Chinese Spring reference sequences. In addition to proteins previously reported in the host response to Fusarium spp., new candidates have emerged in association with resistance or susceptibility, including a group 3 late embryogenesis abundant as a resistance-related protein and a purple acid phosphatase as a susceptibility protein. The protein atlas presented here provides new perspectives on the interaction between F. graminearum and wheat.

Tiny Yet Indispensable Plant MicroRNAs Are Worth to Explore as Key Components for Combating Genotoxic Stresses

Plants being sessile are always exposed to various stresses including biotic and abiotic stresses. Some of these stresses are genotoxic to cells causing DNA damage by forming lesions which include altered bases, cross-links, and breaking of DNA strands, which in turn hamper the genomic integrity. In order to survive through all these adverse conditions, plants have evolved different DNA repair mechanisms. As seen from the mammalian system and different human diseases, various microRNAs (miRNAs) can target the 3'-untranslated region of mRNAs that code for the proteins involved in DNA repair pathways. Since miRNAs play an important role in plant cells by regulating various metabolic pathways, it can also be possible that miRNAs play an important role in DNA repair pathways too. However, till date, only a handful of plant miRNAs have been identified to play important role in combating genotoxic stresses in plants. Limitation of information regarding involvement of miRNAs in DNA repair as well as in ROS scavenging prompted us to gather information about plant miRNAs specific for these tasks. This mini-review aims to present pertinent literature dealing with different genotoxic stresses that cause genome instability as well as plant specific responses to survive the damage. This is intertwined with the involvement of miRNAs in genotoxic stress in plants, challenges of applying miRNAs as a tool to combat DNA damage along with ways to overcome these challenges, and finally, the future prospective of these understudied aspects.

De novo assembly of wheat root transcriptomes and transcriptional signature of longitudinal differentiation

Hidden underground, root systems constitute an important part of the plant for its development, nourishment and sensing the soil environment around it, but we know very little about its genetic regulation in crop plants like wheat. In the present study, we de novo assembled the root transcriptomes in reference cultivar Chinese Spring from RNA-seq reads generated by the 454-GS-FLX and HiSeq platforms. The FLX reads were assembled into 24,986 transcripts with completeness of 54.84%, and the HiSeq reads were assembled into 91,543 high-confidence protein-coding transcripts, 2,404 low-confidence protein-coding transcripts, and 13,181 non-coding transcripts with the completeness of >90%. Combining the FLX and HiSeq assemblies, we assembled a root transcriptome of 92,335 ORF-containing transcripts. Approximately 7% of the coding transcripts and ~2% non-coding transcripts are not present in the current wheat genome assembly. Functional annotation of both assemblies showed similar gene ontology patterns and that ~7% coding and >5% non-coding transcripts are root-specific. Transcription quantification identified 1,728 differentially expressed transcripts between root tips and maturation zone, and functional annotation of these transcripts captured a transcriptional signature of longitudinal development of wheat root. With the transcriptomic resources developed, this study provided the first view of wheat root transcriptome under different developmental zones and laid a foundation for molecular studies of wheat root development and growth using a reverse genetic approach.

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.