Transcriptome profiling of two maize inbreds with distinct responses to Gibberella ear rot disease to identify candidate resistance genes

Evaluation of Resistance Resources and Analysis of Resistance Mechanisms of Maize to Stalk Rot Caused by Fusarium graminearum

Stalk rot is one of the most destructive and widely distributed diseases in maize plants worldwide. Research on the performance and resistance mechanisms of maize against stem rot is constantly improving. In this study, among 120 inbred maize lines infected by Fusarium graminearum using the injection method, 4 lines (3.33%) were highly resistant to stalk rot, 28 lines (23.33%) were resistant, 57 lines (47.50%) were susceptible, and 31 lines (25.84%) were highly susceptible. The inbred lines 18N10118 and 18N10370 were the most resistant and susceptible with disease indices of 7.5 and 75.6, respectively. Treatment of resistant and susceptible maize inbred seedlings with F. graminearum showed that root hair growth of the susceptible inbred lines was significantly inhibited, and a large number of hyphae attached and adsorbed multiple conidia near the root system. However, the resistant inbred lines were delayed and inconspicuous, with only a few hyphae and spores appearing near the root system. Compared with susceptible inbred lines, resistant maize inbred line seedlings treated with F. graminearum exhibited elevated activities of catalase, phenylalanine ammonia-lyase, polyphenol oxidase, and superoxide dismutase. We identified 153 genes related to disease resistance by transcriptome analysis. The mitogen-activated protein kinase signaling and peroxisome pathways mainly regulated the resistance mechanism of maize inbred lines to F. graminearum infection. These two pathways might play an important role in the disease resistance mechanism, and the function of genes in the two pathways must be further studied, which might provide a theoretical basis for further understanding the molecular resistance mechanism of stalk rot and resistance gene mining.

Effectiveness of introgression of resistance loci for Gibberella ear rot from two European flint landraces into adapted elite maize (Zea mays L.)

European flint landraces are a major class of maize possessing favorable alleles for improving host resistance to Gibberella ear rot (GER) disease which reduces yield and contaminates the grains with mycotoxins. However, the incorporation of these landraces into breeding programs requires a clear understanding of the effectiveness of their introgression into elite materials. We evaluated 15 pre-selected doubled haploid (DH) lines from two European flint landraces, "Kemater Landmais Gelb" (KE) and "Petkuser Ferdinand Rot" (PE), together with two adapted elite flint lines and seven standard lines for GER severity as the main trait, and several adaptation traits (plant height, days to silking, seed-set, plant vigor) across four environments. From this evaluation, three KE DH lines and one PE DH line, with the lowest GER severity, were selected and used as donor parents that were crossed with the two adapted and GER susceptible flint lines (Flint1 and Flint2) to develop six bi-parental DH populations with 34-145 DH lines each. Each DH population was evaluated across two locations. Correlations between GER severity, which was the target trait, and adaptation traits were weak (-0.02 to 0.19). GER severity of lines from PE landrace was on average 2-fold higher than lines from KE landrace, indicating a clear superiority of the KE landrace lines. Mean GER severity of the DH populations was 39.4-61.0% lower than the adapted elite flint lines. All KE-derived DH populations were on average more resistant (27.0-36.7%) than the PE-derived population (51.0%). Highly resistant lines (1.3-5.2%) were found in all of the populations, suggesting that the DH populations can be successfully integrated into elite breeding programs. The findings demonstrate that selected KE landrace lines used as donors were effective in improving GER resistance of the adapted elite inbreds.

Comprehensive meta-analysis of QTL and gene expression studies identify candidate genes associated with Aspergillus flavus resistance in maize

Aflatoxin (AF) contamination, caused by Aspergillus flavus, compromises the food safety and marketability of commodities, such as maize, cotton, peanuts, and tree nuts. Multigenic inheritance of AF resistance impedes conventional introgression of resistance traits into high-yielding commercial maize varieties. Several AF resistance-associated quantitative trait loci (QTLs) and markers have been reported from multiple biparental mapping and genome-wide association studies (GWAS) in maize. However, QTLs with large confidence intervals (CI) explaining inconsistent phenotypic variance limit their use in marker-assisted selection. Meta-analysis of published QTLs can identify significant meta-QTLs (MQTLs) with a narrower CI for reliable identification of genes and linked markers for AF resistance. Using 276 out of 356 reported QTLs controlling resistance to A. flavus infection and AF contamination in maize, we identified 58 MQTLs on all 10 chromosomes with a 66.5% reduction in the average CI. Similarly, a meta-analysis of maize genes differentially expressed in response to (a)biotic stresses from the to-date published literature identified 591 genes putatively responding to only A. flavus infection, of which 14 were significantly differentially expressed (-1.0 ≤ Log2Fc ≥ 1.0; p ≤ 0.05). Eight MQTLs were validated by their colocalization with 14 A. flavus resistance-associated SNPs identified from GWAS in maize. A total of 15 genes were physically close between the MQTL intervals and SNPs. Assessment of 12 MQTL-linked SSR markers identified three markers that could discriminate 14 and eight cultivars with resistance and susceptible responses, respectively. A comprehensive meta-analysis of QTLs and differentially expressed genes led to the identification of genes and makers for their potential application in marker-assisted breeding of A. flavus-resistant maize varieties.

Cell wall-related genes and lignin accumulation contribute to the root resistance in different maize (Zea mays L.) genotypes to Fusarium verticillioides (Sacc.) Nirenberg infection

Introduction

The fungal pathogen Fusarium verticillioides (Sacc.) Nirenberg (Fv) causes considerable agricultural and economic losses and is harmful to animal and human health. Fv can infect maize throughout its long agricultural cycle, and root infection drastically affects maize growth and yield.

Methods

The root cell wall is the first physical and defensive barrier against soilborne pathogens such as Fv. This study compares two contrasting genotypes of maize (Zea mays L.) roots that are resistant (RES) or susceptible (SUS) to Fv infection by using transcriptomics, fluorescence, scanning electron microscopy analyses, and ddPCR.

Results

Seeds were infected with a highly virulent local Fv isolate. Although Fv infected both the RES and SUS genotypes, infection occurred faster in SUS, notably showing a difference of three to four days. In addition, root infections in RES were less severe in comparison to SUS infections. Comparative transcriptomics (rate +Fv/control) were performed seven days after inoculation (DAI). The analysis of differentially expressed genes (DEGs) in each rate revealed 733 and 559 unique transcripts that were significantly (P ≤0.05) up and downregulated in RES (+Fv/C) and SUS (+Fv/C), respectively. KEGG pathway enrichment analysis identified coumarin and furanocoumarin biosynthesis, phenylpropanoid biosynthesis, and plant-pathogen interaction pathways as being highly enriched with specific genes involved in cell wall modifications in the RES genotype, whereas the SUS genotype mainly displayed a repressed plant-pathogen interaction pathway and did not show any enriched cell wall genes. In particular, cell wall-related gene expression showed a higher level in RES than in SUS under Fv infection. Analysis of DEG abundance made it possible to identify transcripts involved in response to abiotic and biotic stresses, biosynthetic and catabolic processes, pectin biosynthesis, phenylpropanoid metabolism, and cell wall biosynthesis and organization. Root histological analysis in RES showed an increase in lignified cells in the sclerenchymatous hypodermis zone during Fv infection.

Discussion

These differences in the cell wall and lignification could be related to an enhanced degradation of the root hairs and the epidermis cell wall in SUS, as was visualized by SEM. These findings reveal that components of the root cell wall are important against Fv infection and possibly other soilborne phytopathogens.

Meta-analysis and co-expression analysis revealed stable QTL and candidate genes conferring resistances to Fusarium and Gibberella ear rots while reducing mycotoxin contamination in maize

Fusarium (FER) and Gibberella ear rots (GER) are the two most devastating diseases of maize (Zea mays L.) which reduce yield and affect grain quality worldwide, especially by contamination with mycotoxins. Genetic improvement of host resistance to effectively tackle FER and GER diseases requires the identification of stable quantitative trait loci (QTL) to facilitate the application of genomics-assisted breeding for improving selection efficiency in breeding programs. We applied improved meta-analysis algorithms to re-analyze 224 QTL identified in 15 studies based on dense genome-wide single nucleotide polymorphisms (SNP) in order to identify meta-QTL (MQTL) and colocalized genomic loci for fumonisin (FUM) and deoxynivalenol (DON) accumulation, silk (SR) and kernel (KR) resistances of both FER and GER, kernel dry-down rate (KDD) and husk coverage (HC). A high-resolution genetic consensus map with 36,243 loci was constructed and enabled the projection of 164 of the 224 collected QTL. Candidate genes (CG) mining was performed within the most refined MQTL, and identified CG were cross-validated using publicly available transcriptomic data of maize under Fusarium graminearum infection. The meta-analysis revealed 40 MQTL, of which 29 were associated each with 2-5 FER- and/or GER-related traits. Twenty-eight of the 40 MQTL were common to both FER and GER resistances and 19 MQTL were common to silk and kernel resistances. Fourteen most refined MQTL on chromosomes 1, 2, 3, 4, 7 and 9 harbored a total of 2,272 CG. Cross-validation identified 59 of these CG as responsive to FER and/or GER diseases. MQTL ZmMQTL2.2, ZmMQTL9.2 and ZmMQTL9.4 harbored promising resistance genes, of which GRMZM2G011151 and GRMZM2G093092 were specific to the resistant line for both diseases and encoded "terpene synthase21 (tps21)" and "flavonoid O-methyltransferase2 (fomt2)", respectively. Our findings revealed stable refined MQTL harboring promising candidate genes for use in breeding programs for improving FER and GER resistances with reduced mycotoxin accumulation. These candidate genes can be transferred into elite cultivars by integrating refined MQTL into genomics-assisted backcross breeding strategies.

Autotoxin affects the rhizosphere microbial community structure by influencing the secretory characteristics of grapevine roots

Autotoxins secreted by roots into the soil can trigger rhizosphere microecological imbalances and affect root secretory properties resulting in conditions such as replanting disease. However, information on the effect of autotoxins on root secretion characteristics and regulation of the composition of rhizosphere microorganisms by altered root exudates is limited. In this study, autotoxin ρ-hydroxybenzoic acid (4-HBA) was added to the soil of potted grapevine seedlings, CO₂ pulse-labeling, and DNA stable isotope probing were used to track the rhizosphere microbiome that assimilates root exudates. Bacterial and fungal microbiomes that assimilated plant-derived carbon were identified by high-throughput sequencing. Results showed that 4-HBA treatment altered bacterial and fungal communities in ¹³C-labeled organisms, with a lower abundance of beneficial bacteria (e.g., Gemmatimonas, Streptomyces, and Bacillus) and a higher abundance of potential pathogen fungi (e.g., Fusarium, Neocosmospora, Gibberella, and Fusicolla) by changing the composition of root exudates. The exogenous addition of upregulated compound mixtures of root exudates reduced the abundance of beneficial bacterial Bacillus and increased the abundance of potential pathogen fungi Gibberella. These results suggest that 4-HBA can alter root secretion properties and altered root exudates may enrich certain potential pathogens and reduce certain beneficial bacteria, thereby unbalancing the structure of the rhizosphere microbial community.

Network and Evolutionary Analysis Reveals Candidate Genes of Membrane Trafficking Involved in Maize Seed Development and Immune Response

The plant membrane-trafficking system plays a crucial role in maintaining proper cellular functions and responding to various developmental and environmental cues. Thus far, our knowledge of the maize membrane-trafficking system is still limited. In this study, we systematically identified 479 membrane-trafficking genes from the maize genome using orthology search and studied their functions by integrating transcriptome and evolution analyses. These genes encode the components of coated vesicles, AP complexes, autophagy, ESCRTs, retromers, Rab GTPases, tethering factors, and SNAREs. The maize genes exhibited diverse but coordinated expression patterns, with 249 genes showing elevated expression in reproductive tissues. Further WGCNA analysis revealed that five COPII components and four Rab GTPases had high connectivity with protein biosynthesis during endosperm development and that eight components of autophagy, ESCRT, Rab, and SNARE were strongly co-upregulated with defense-related genes and/or with secondary metabolic processes to confer basal resistance to Fusarium graminearum. In addition, we identified 39 membrane-trafficking genes with strong selection signals during maize domestication and/or improvement. Among them, ZmSec23a and ZmVPS37A were selected for kernel oil production during improvement and pathogen resistance during domestication, respectively. In summary, these findings will provide important hints for future appreciation of the functions of membrane-trafficking genes in maize.

Transcriptomic Analysis Reveals Candidate Genes Responding Maize Gray Leaf Spot Caused by Cercospora zeina

Gray leaf spot (GLS), caused by the fungal pathogen Cercospora zeina (C. zeina), is one of the most destructive soil-borne diseases in maize (Zea mays L.), and severely reduces maize production in Southwest China. However, the mechanism of resistance to GLS is not clear and few resistant alleles have been identified. Two maize inbred lines, which were shown to be resistant (R6) and susceptible (S8) to GLS, were injected by C. zeina spore suspensions. Transcriptome analysis was carried out with leaf tissue at 0, 6, 24, 144, and 240 h after inoculation. Compared with 0 h of inoculation, a total of 667 and 419 stable common differentially expressed genes (DEGs) were found in the resistant and susceptible lines across the four timepoints, respectively. The DEGs were usually enriched in 'response to stimulus' and 'response to stress' in GO term analysis, and 'plant-pathogen interaction', 'MAPK signaling pathways', and 'plant hormone signal transduction' pathways, which were related to maize's response to GLS, were enriched in KEGG analysis. Weighted-Genes Co-expression Network Analysis (WGCNA) identified two modules, while twenty hub genes identified from these indicated that plant hormone signaling, calcium signaling pathways, and transcription factors played a central role in GLS sensing and response. Combing DEGs and QTL mapping, five genes were identified as the consensus genes for the resistance of GLS. Two genes, were both putative Leucine-rich repeat protein kinase family proteins, specifically expressed in R6. In summary, our results can provide resources for gene mining and exploring the mechanism of resistance to GLS in maize.

Omics for the Improvement of Abiotic, Biotic, and Agronomic Traits in Major Cereal Crops: Applications, Challenges, and Prospects

Omics technologies, namely genomics, transcriptomics, proteomics, metabolomics, and phenomics, are becoming an integral part of virtually every commercial cereal crop breeding program, as they provide substantial dividends per unit time in both pre-breeding and breeding phases. Continuous advances in omics assure time efficiency and cost benefits to improve cereal crops. This review provides a comprehensive overview of the established omics methods in five major cereals, namely rice, sorghum, maize, barley, and bread wheat. We cover the evolution of technologies in each omics section independently and concentrate on their use to improve economically important agronomic as well as biotic and abiotic stress-related traits. Advancements in the (1) identification, mapping, and sequencing of molecular/structural variants; (2) high-density transcriptomics data to study gene expression patterns; (3) global and targeted proteome profiling to study protein structure and interaction; (4) metabolomic profiling to quantify organ-level, small-density metabolites, and their composition; and (5) high-resolution, high-throughput, image-based phenomics approaches are surveyed in this review.

Identification and Cloning of a CC-NBS-NBS-LRR Gene as a Candidate of Pm40 by Integrated Analysis of Both the Available Transcriptional Data and Published Linkage Mapping

Wheat powdery mildew, caused by the obligate parasite Blumeria graminis f. sp. tritici, severely reduces wheat yields. Identifying durable and effective genes against wheat powdery mildew and further transferring them into wheat cultivars is important for finally controlling this disease in wheat production. Pm40 has been widely used in wheat breeding programs in Southwest China due to the spectrum and potentially durable resistance to powdery mildew. In the present study, a resistance test demonstrated that Pm40 is still effective against the Bgt race E20. We identified and cloned the TraesCS7B01G164000 with a total length of 4883 bp, including three exons and two introns, and encoded a protein carrying the CC-NBS-NBS-LRR domain in the Pm40-linked region flanked by two EST markers, BF478514 and BF291338, by integrating analysis of gene annotation in wheat reference genome and both sequence and expression difference in available transcriptome data. Two missense mutations were detected at positions 68 and 83 in the CC domain. The results of both cosegregation linkage analysis and qRT-PCR also suggested that TraesCS7B01G164000 was a potential candidate gene of Pm40. This study allowed us to move toward the final successfully clone and apply Pm40 in wheat resistance improvement by gene engineering.

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

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

Comparative Transcriptome Analysis Reveals Regulatory Networks during the Maize Ear Shank Elongation Process

In maize, the ear shank is a short branch that connects the ear to the stalk. The length of the ear shank mainly affects the transportation of photosynthetic products to the ear, and also influences the dehydration of the grain by adjusting the tightness of the husks. However, the molecular mechanisms of maize shank elongation have rarely been described. It has been reported that the maize ear shank length is a quantitative trait, but its genetic basis is still unclear. In this study, RNA-seq was performed to explore the transcriptional dynamics and determine the key genes involved in maize shank elongation at four different developmental stages. A total of 8145 differentially expressed genes (DEGs) were identified, including 729 transcription factors (TFs). Some important genes which participate in shank elongation were detected via function annotation and temporal expression pattern analyses, including genes related to signal transduction hormones (auxin, brassinosteroids, gibberellin, etc.), xyloglucan and xyloglucan xyloglucosyl transferase, and transcription factor families. The results provide insights into the genetic architecture of maize ear shanks and developing new varieties with ideal ear shank lengths, enabling adjustments for mechanized harvesting in the future.

Commonly and Specifically Activated Defense Responses in Maize Disease Lesion Mimic Mutants Revealed by Integrated Transcriptomics and Metabolomics Analysis

Disease lesion mimic (Les/les) mutants display disease-like spontaneous lesions in the absence of pathogen infection, implying the constitutive activation of defense responses. However, the genetic and biochemical bases underlying the activated defense responses in those mutants remain largely unknown. Here, we performed integrated transcriptomics and metabolomics analysis on three typical maize Les mutants Les4, Les10, and Les17 with large, medium, and small lesion size, respectively, thereby dissecting the activated defense responses at the transcriptional and metabolomic level. A total of 1,714, 4,887, and 1,625 differentially expressed genes (DEGs) were identified in Les4, Les10, and Les17, respectively. Among them, 570, 3,299, and 447 specific differentially expressed genes (SGs) were identified, implying a specific function of each LES gene. In addition, 480 common differentially expressed genes (CGs) and 42 common differentially accumulated metabolites (CMs) were identified in all Les mutants, suggesting the robust activation of shared signaling pathways. Intriguingly, substantial analysis of the CGs indicated that genes involved in the programmed cell death, defense responses, and phenylpropanoid and terpenoid biosynthesis were most commonly activated. Genes involved in photosynthetic biosynthesis, however, were generally repressed. Consistently, the dominant CMs identified were phenylpropanoids and flavonoids. In particular, lignin, the phenylpropanoid-based polymer, was significantly increased in all three mutants. These data collectively imply that transcriptional activation of defense-related gene expression; increase of phenylpropanoid, lignin, flavonoid, and terpenoid biosynthesis; and inhibition of photosynthesis are generalnatures associated with the lesion formation and constitutively activated defense responses in those mutants. Further studies on the identified SGs and CGs will shed new light on the function of each LES gene as well as the regulatory network of defense responses in maize.

Comparative transcriptome profiling identifies maize line specificity of fungal effectors in the maize-Ustilago maydis interaction

The biotrophic pathogen Ustilago maydis causes smut disease on maize (Zea mays) and induces the formation of tumours on all aerial parts of the plant. Unlike in other biotrophic interactions, no gene-for-gene interactions have been identified in the maize-U. maydis pathosystem. Thus, maize resistance to U. maydis is considered a polygenic, quantitative trait. Here, we study the molecular mechanisms of quantitative disease resistance (QDR) in maize, and how U. maydis interferes with its components. Based on quantitative scoring of disease symptoms in 26 maize lines, we performed an RNA sequencing (RNA-Seq) analysis of six U. maydis-infected maize lines of highly distinct resistance levels. The different maize lines showed specific responses of diverse cellular processes to U. maydis infection. For U. maydis, our analysis identified 406 genes being differentially expressed between maize lines, of which 102 encode predicted effector proteins. Based on this analysis, we generated U. maydis CRISPR/Cas9 knock-out mutants for selected candidate effector sets. After infections of different maize lines with the fungal mutants, RNA-Seq analysis identified effectors with quantitative, maize line-specific virulence functions, and revealed auxin-related processes as a possible target for one of them. Thus, we show that both transcriptional activity and virulence function of fungal effector genes are modified according to the infected maize line, providing insights into the molecular mechanisms underlying QDR in the maize-U. maydis interaction.

Molecular evidence of the avocado defense response to Fusarium kuroshium infection: a deep transcriptome analysis using RNA-Seq

Fusarium kuroshium is a novel member of the Ambrosia Fusarium Clade (AFC) that has been recognized as one of the symbionts of the invasive Kuroshio shot hole borer, an Asian ambrosia beetle. This complex is considered the causal agent of Fusarium dieback, a disease that has severely threatened natural forests, landscape trees, and avocado orchards in the last 8 years. Despite the interest in this species, the molecular responses of both the host and F. kuroshium during the infection process and disease establishment remain unknown. In this work, we established an in vitro pathosystem using Hass avocado stems inoculated with F. kuroshium to investigate differential gene expression at 1, 4, 7 and 14 days post-inoculation. RNA-seq technology allowed us to obtain data from both the plant and the fungus, and the sequences obtained from both organisms were analyzed independently. The pathosystem established was able to mimic Fusarium dieback symptoms, such as carbohydrate exudation, necrosis, and vascular tissue discoloration. The results provide interesting evidence regarding the genes that may play roles in the avocado defense response to Fusarium dieback disease. The avocado data set comprised a coding sequence collection of 51,379 UniGenes, from which 2,403 (4.67%) were identified as differentially expressed. The global expression analysis showed that F. kuroshium responsive UniGenes can be clustered into six groups according to their expression profiles. The biologically relevant functional categories that were identified included photosynthesis as well as responses to stress, hormones, abscisic acid, and water deprivation. Additionally, processes such as oxidation-reduction, organization and biogenesis of the cell wall and polysaccharide metabolism were detected. Moreover, we identified orthologues of nucleotide-binding leucine-rich receptors, and their possible action mode was analyzed. In F. kuroshium, we identified 57 differentially expressed genes. Interestingly, the alcohol metabolic process biological category had the highest number of upregulated genes, and the enzyme group in this category may play an important role in the mechanisms of secondary metabolite detoxification. Hydrolytic enzymes, such as endoglucanases and a pectate lyase, were also identified, as well as some proteases. In conclusion, our research was conducted mainly to explain how the vascular tissue of a recognized host of the ambrosia complex responds during F. kuroshium infection since Fusarium dieback is an ambrosia beetle-vectored disease and many variables facilitate its establishment.

The eukaryotic MEP-pathway genes are evolutionarily conserved and originated from Chlaymidia and cyanobacteria

Background

Isoprenoids are the most ancient and essential class of metabolites produced in all organisms, either via mevalonate (MVA)-and/or methylerythritol phosphate (MEP)-pathways. The MEP-pathway is present in all plastid-bearing organisms and most eubacteria. However, no comprehensive study reveals the origination and evolutionary characteristics of MEP-pathway genes in eukaryotes.

Results

Here, detailed bioinformatics analyses of the MEP-pathway provide an in-depth understanding the evolutionary history of this indispensable biochemical route, and offer a basis for the co-existence of the cytosolic MVA- and plastidial MEP-pathway in plants given the established exchange of the end products between the two isoprenoid-biosynthesis pathways. Here, phylogenetic analyses establish the contributions of both cyanobacteria and Chlamydiae sequences to the plant’s MEP-pathway genes. Moreover, Phylogenetic and inter-species syntenic block analyses demonstrate that six of the seven MEP-pathway genes have predominantly remained as single-copy in land plants in spite of multiple whole-genome duplication events (WGDs). Substitution rate and domain studies display the evolutionary conservation of these genes, reinforced by their high expression levels. Distinct phenotypic variation among plants with reduced expression levels of individual MEP-pathway genes confirm the indispensable function of each nuclear-encoded plastid-targeted MEP-pathway enzyme in plant growth and development.

Conclusion

Collectively, these findings reveal the polyphyletic origin and restrict conservation of MEP-pathway genes, and reinforce the potential function of the individual enzymes beyond production of the isoprenoids intermediates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07448-x.

A temporal hierarchy underpins the transcription factor-DNA interactome of the maize UPR

Adverse environmental conditions reduce crop productivity and often increase the load of unfolded or misfolded proteins in the endoplasmic reticulum (ER). This potentially lethal condition, known as ER stress, is buffered by the unfolded protein response (UPR), a set of signaling pathways designed to either recover ER functionality or ignite programmed cell death. Despite the biological significance of the UPR to the life of the organism, the regulatory transcriptional landscape underpinning ER stress management is largely unmapped, especially in crops. To fill this significant knowledge gap, we performed a large-scale systems-level analysis of the protein-DNA interaction (PDI) network in maize (Zea mays). Using 23 promoter fragments of six UPR marker genes in a high-throughput enhanced yeast one-hybrid assay, we identified a highly interconnected network of 262 transcription factors (TFs) associated with significant biological traits and 831 PDIs underlying the UPR. We established a temporal hierarchy of TF binding to gene promoters within the same family as well as across different families of TFs. Cistrome analysis revealed the dynamic activities of a variety of cis-regulatory elements (CREs) in ER stress-responsive gene promoters. By integrating the cistrome results into a TF network analysis, we mapped a subnetwork of TFs associated with a CRE that may contribute to UPR management. Finally, we validated the role of a predicted network hub gene using the Arabidopsis system. The PDIs, TF networks, and CREs identified in our work are foundational resources for understanding transcription-regulatory mechanisms in the stress responses and crop improvement.

Genomics-Assisted Breeding for Quantitative Disease Resistances in Small-Grain Cereals and Maize

Generating genomics-driven knowledge opens a way to accelerate the resistance breeding process by family or population mapping and genomic selection. Important prerequisites are large populations that are genomically analyzed by medium- to high-density marker arrays and extensive phenotyping across locations and years of the same populations. The latter is important to train a genomic model that is used to predict genomic estimated breeding values of phenotypically untested genotypes. After reviewing the specific features of quantitative resistances and the basic genomic techniques, the possibilities for genomics-assisted breeding are evaluated for six pathosystems with hemi-biotrophic fungi: Small-grain cereals/Fusarium head blight (FHB), wheat/Septoria tritici blotch (STB) and Septoria nodorum blotch (SNB), maize/Gibberella ear rot (GER) and Fusarium ear rot (FER), maize/Northern corn leaf blight (NCLB). Typically, all quantitative disease resistances are caused by hundreds of QTL scattered across the whole genome, but often available in hotspots as exemplified for NCLB resistance in maize. Because all crops are suffering from many diseases, multi-disease resistance (MDR) is an attractive aim that can be selected by specific MDR QTL. Finally, the integration of genomic data in the breeding process for introgression of genetic resources and for the improvement within elite materials is discussed.

Genome-wide identification of the maize 2OGD superfamily genes and their response to Fusarium verticillioides and Fusarium graminearum

In maize, eat rot and stalk rot caused by Fusarium verticillioides and Fusarium graminearum lead to contamination of moldy grains to produce mycotoxins. Identification of resistance genes against these pathogens for maize breeding is an effective way for disease control. Several 2-oxoglutarate-dependent dioxygenase (2OGD) proteins have been found to confer resistance to different pathogens in diverse plant species. However, little is known about the 2OGD superfamily in maize. Here, we identified 103 putative 2OGD genes in maize from a genome-wide analysis, and divided them into three classes - DOXA, DOXB, and DOXC. We further comprehensively investigated their gene structure, chromosome distribution, phylogenetic tree, gene-function enrichment, and expression profiles among different tissues. The genes encoding three 2OGD proteins, ACO, F3H, and NCS involved in ethylene biosynthesis, flavonoids biosynthesis, and alkaloids biosynthesis pathways, respectively, were identified to be induced by F. verticillioides and F. graminearum. The promoters of the three genes contain the binding sites for the transcription factor ZmDOF and ZmHSF, which are also induced by the two pathogens. The results imply that the three 2OGDs and the two transcription factors might be involved in the resistance to the two pathogens. This study provided a comprehensive understanding of the 2OGD superfamily in maize and laid the foundation for the further functional analysis of their roles in maize resistance to eat rot and stalk rot.

Robustness of plant quantitative disease resistance is provided by a decentralized immune network

Quantitative disease resistance (QDR) represents the predominant form of resistance in natural populations and crops. Surprisingly, very limited information exists on the biomolecular network of the signaling machineries underlying this form of plant immunity. This lack of information may result from its complex and quantitative nature. Here, we used an integrative approach including genomics, network reconstruction, and mutational analysis to identify and validate molecular networks that control QDR in Arabidopsis thaliana in response to the bacterial pathogen Xanthomonas campestris To tackle this challenge, we first performed a transcriptomic analysis focused on the early stages of infection and using transgenic lines deregulated for the expression of RKS1, a gene underlying a QTL conferring quantitative and broad-spectrum resistance to XcampestrisRKS1-dependent gene expression was shown to involve multiple cellular activities (signaling, transport, and metabolism processes), mainly distinct from effector-triggered immunity (ETI) and pathogen-associated molecular pattern (PAMP)-triggered immunity (PTI) responses already characterized in Athaliana Protein-protein interaction network reconstitution then revealed a highly interconnected and distributed RKS1-dependent network, organized in five gene modules. Finally, knockout mutants for 41 genes belonging to the different functional modules of the network revealed that 76% of the genes and all gene modules participate partially in RKS1-mediated resistance. However, these functional modules exhibit differential robustness to genetic mutations, indicating that, within the decentralized structure of the QDR network, some modules are more resilient than others. In conclusion, our work sheds light on the complexity of QDR and provides comprehensive understanding of a QDR immune network.

MicroRNAs Are Involved in Maize Immunity Against Fusarium verticillioides Ear Rot

Fusarium ear rot (FER) caused by Fusarium verticillioides is one of the most common diseases affecting maize production worldwide. FER results in severe yield losses and grain contamination with health-threatening mycotoxins. Although most studies to date have focused on comprehensive analysis of gene regulation in maize during defense responses against F. verticillioides infection, less is known about the role of microRNAs (miRNAs) in this process. We used deep sequencing to compare small RNA libraries from the maize kernels of susceptible (N6) or resistant (BT-1) inbred lines from uninfected plants and upon F. verticillioides infection. We found that pathogen exposure was accompanied by dynamic alterations in expression levels of multiple miRNAs, including new members of previously annotated miRNA families. A combination of transcriptomic, degradomic, and bioinformatics analyses revealed that F. verticillioides-responsive miRNAs and their potential target genes displayed opposite expression patterns in the susceptible and resistant genotypes. Functional category analysis uncovered preferential enrichment of the pathogen-responsive miRNAs and their targets in the phenylpropanoid metabolic processes, plant-pathogen interactions, and plant phytohormone signal transduction pathways. Furthermore, transgenic maize plants overexpressing miR408b exhibited reduced resistance to F. verticillioides infection in a susceptible maize line. These findings provide new insights into the regulatory roles of miRNAs in maize immunity against FER and new resources for breeding disease resistance into maize.

Maize ZmFNSI Homologs Interact with an NLR Protein to Modulate Hypersensitive Response

Nucleotide binding, leucine-rich-repeat (NLR) proteins are the major class of resistance (R) proteins used by plants to defend against pathogen infection. The recognition between NLRs and their cognate pathogen effectors usually triggers a rapid localized cell death, termed the hypersensitive response (HR). Flavone synthase I (FNSI) is one of the key enzymes in the flavone biosynthesis pathway. It also displays salicylic acid (SA) 5-hydroxylase (S5H) activity. A close homolog of FNSI/S5H displays SA 3-hydroxylase (S3H) activity. Both FNSI/S5H and S3H play important roles in plant innate immunity. However, the underlying molecular mechanisms and the relationship between S5H and S3H with the NLR-mediated HR are not known in any plant species. In this study, we identified three genes encoding ZmFNSI-1, ZmFNSI-2 and ZmS3H that are significantly upregulated in a maize line carrying an autoactive NLR Rp1-D21 mutant. Functional analysis showed that ZmFNSI-1 and ZmFNSI-2, but not ZmS3H, suppressed HR conferred by Rp1-D21 and its signaling domain CCD21 when transiently expressed in N. benthamiana. ZmFNSI-1 and ZmFNSI-2 physically interacted with CCD21. Furthermore, ZmFNSI-1 and ZmFNSI-2 interacted with HCT, a key enzyme in lignin biosynthesis pathway, which can also suppress Rp1-D21-mediated HR. These results lay the foundation for the further functional analysis of the roles of FNSI in plant innate immunity.

Proteomic Changes during MCMV Infection Revealed by iTRAQ Quantitative Proteomic Analysis in Maize

Maize chlorotic mottle virus (MCMV) has been occurring frequently worldwide and causes severe yield losses in maize (Zea mays). To better investigate the destructive effects of MCMV infection on maize plants, isobaric tagging for relative and absolute quantitation (iTRAQ)-based comparative proteomic analysis was performed on MCMV infected maize cv. B73. A total of 972 differentially abundant proteins (DAPs), including 661 proteins with increased abundance and 311 proteins with reduced abundance, were identified in response to MCMV infection. Functional annotations of DAPs and measurement of photosynthetic activity revealed that photosynthesis was decreased, while the abundance of ribosomal proteins, proteins related to stress responses, oxidation-reduction and redox homeostasis was altered significantly during MCMV infection. Two DAPs, disulfide isomerases like protein ZmPDIL-1 and peroxiredoxin family protein ZmPrx5, were further analyzed for their roles during MCMV infection through cucumber mosaic virus-based virus-induced gene silencing (CMV-VIGS). The accumulation of MCMV was suppressed in ZmPDIL-1-silenced or ZmPrx5-silenced B73 maize, suggesting ZmPDIL-1 and ZmPrx5 might enhance host susceptibility to MCMV infection.

Genomics-assisted breeding for ear rot resistances and reduced mycotoxin contamination in maize: methods, advances and prospects

Genetic mapping, genomic profiling and bioinformatic approaches were used to identify putative resistance genes for ear rots and low mycotoxin contamination in maize. Genomic selection seems to have good perspectives. Maize is globally an indispensable crop for humans and livestock. About 30% of yield is lost by fungal diseases with Gibberella, Fusarium and Aspergillus ear rots (ERs) having a high economic impact in most maize-growing regions of the world. They reduce not only yield, but also contaminate grains with mycotoxins like deoxynivalenol, zearalenone, fumonisins and aflatoxins, respectively. These mycotoxins pose serious health problems to humans and animals. A number of studies have been conducted to dissect the genetic architecture of resistance to these three major ear rots over the past decade. The review concentrates on studies carried out to locate quantitative trait loci (QTL) and candidate genes (CG) on the maize genome as well as the application of genomic selection in maize for resistance against Fusarium graminearum, Fusarium verticillioides and Aspergillus flavus. QTL studies by linkage or genome-wide association mapping, omic technologies (genomics, proteomics, transcriptomics and metabolomics) and bioinformatics are the methods used in the current studies to propose resistance genes against ear rot pathogens. Though a number of QTL and CG are reported, only a few specific genes were found to directly confer ER resistance in maize. A combination of two or more gene identification methods would provide a more powerful and reliable tool. Genomic selection seems to be promising for ER resistance breeding, but there are only a limited number of studies in this area. A strategy that can accurately validate and predict genotypes with major effect QTL and CG for selection will be worthwhile for practical breeding against ERs and mycotoxin contamination in maize.

Transcriptome analysis reveals the molecular mechanisms of the defense response to gray leaf spot disease in maize

BACKGROUND

Gray leaf spot (GLS), which is caused by the necrotrophic fungi Cercospora zeae-maydis and Cercospora zeina, is one of the most impactful diseases in maize worldwide. The aim of the present study is to identify the resistance genes and understand the molecular mechanisms for GLS resistance.

RESULTS

Two cultivars, 'Yayu889' and 'Zhenghong532,' which are distinguished as resistant and susceptible cultivars, respectively, were challenged with the GLS disease and a RNA-seq experiment was conducted on infected plants at 81, 89, 91, and 93 days post planting (dap). Compared with the beginning stage at 81 dap, 4666, 1733, and 1166 differentially expressed genes (DEGs) were identified at 89, 91, and 93 dap, respectively, in 'Yayu889,' while relatively fewer, i.e., 4713, 881, and 722 DEGs, were identified in 'Zhenghong532.' Multiple pathways involved in the response of maize to GLS, including 'response to salicylic acid,' 'protein phosphorylation,' 'oxidation-reduction process,' and 'carotenoid biosynthetic process,' were enriched by combining differential expression analysis and Weighted Gene Co-expression Network Analysis (WGCNA). The expression of 12 candidate resistance proteins in these pathways were quantified by the multiple reaction monitoring (MRM) method. This approach identified two candidate resistance proteins, a calmodulin-like protein and a leucine-rich repeat receptor-like protein kinase with SNPs that were located in QTL regions for GLS resistance. Metabolic analysis showed that, compared with 'Zhenghong532,' the amount of salicylic acid (SA) and total carotenoids in 'Yayu889' increased, while peroxidase activity decreased during the early infection stages, suggesting that increased levels of SA, carotenoids, and reactive oxygen species (ROS) may enhance the defense response of 'Yayu889' to GLS.

CONCLUSION

By combining transcriptome and proteome analyses with comparisons of resistance QTL regions, calmodulin-like protein and leucine-rich repeat receptor-like protein kinase were identified as candidate GLS resistance proteins. Moreover, we found that the metabolic pathways for ROS, SA, and carotenoids are especially active in the resistant cultivar. These findings could lead to a better understanding of the GLS resistance mechanisms and facilitate the breeding of GLS-resistant maize cultivars.

Plant host habitat and root exudates shape fungal diversity

The rhizospheric microbiome is clearly affected by plant species and certain of their functional traits. These functional traits allow plants to adapt to their environmental conditions by acquiring or conserving nutrients, thus defining different ecological resource-use plant strategies. In the present study, we investigated whether plants with one of the two nutrient-use strategies (conservative versus exploitative) could influence fungal communities involved in soil organic matter degradation and root exudate assimilation, as well as those colonizing root tissues. We applied a DNA-based, stable-isotope probing (DNA-SIP) approach to four grass species distributed along a gradient of plant nutrient resource strategies, ranging from conservative to exploitative species, and analyzed their associated mycobiota composition using a fungal internal transcribed spacer (ITS) and Glomeromycotina 18S rRNA gene metabarcoding approach. Our results demonstrated that fungal taxa associated with exploitative and conservative plants could be separated into two general categories according to their location: generalists, which are broadly distributed among plants from each strategy and represent the core mycobiota of soil organic matter degraders, root exudate consumers in the root-adhering soil, and root colonizers; and specialists, which are locally abundant in one species and more specifically involved in soil organic matter degradation or root exudate assimilation on the root-adhering soil and the root tissues. Interestingly, for arbuscular mycorrhizal fungi analysis, all plant roots were mainly colonized by Glomus species, whereas an increased diversity of Glomeromycotina genera was observed for the exploitative plant species Dactylis glomerata.