Integrated Assays of Genome-Wide Association Study, Multi-Omics Co-Localization, and Machine Learning Associated Calcium Signaling Genes with Oilseed Rape Resistance to Sclerotinia sclerotiorum.

Sclerotinia sclerotiorum (Ss) is one of the most devastating fungal pathogens, causing huge yield loss in multiple economically important crops including oilseed rape. Plant resistance to Ss pertains to quantitative disease resistance (QDR) controlled by multiple minor genes. Genome-wide identification of genes involved in QDR to Ss is yet to be conducted. In this study, we integrated several assays including genome-wide association study (GWAS), multi-omics co-localization, and machine learning prediction to identify, on a genome-wide scale, genes involved in the oilseed rape QDR to Ss. Employing GWAS and multi-omics co-localization, we identified seven resistance-associated loci (RALs) associated with oilseed rape resistance to Ss. Furthermore, we developed a machine learning algorithm and named it Integrative Multi-Omics Analysis and Machine Learning for Target Gene Prediction (iMAP), which integrates multi-omics data to rapidly predict disease resistance-related genes within a broad chromosomal region. Through iMAP based on the identified RALs, we revealed multiple calcium signaling genes related to the QDR to Ss. Population-level analysis of selective sweeps and haplotypes of variants confirmed the positive selection of the predicted calcium signaling genes during evolution. Overall, this study has developed an algorithm that integrates multi-omics data and machine learning methods, providing a powerful tool for predicting target genes associated with specific traits. Furthermore, it makes a basis for further understanding the role and mechanisms of calcium signaling genes in the QDR to Ss.

Comparative Transcriptomics of Soybean Genotypes with Partial Resistance Toward Phytophthora sojae, Conrad, and M92-220 to Moderately Susceptible Fast Neutron Mutant Soybeans and Sloan.

The breeding of disease-resistant soybeans cultivars to manage Phytophthora root and stem rot caused by the pathogen Phytophthora sojae involves combining quantitative disease resistance (QDR) and Rps gene-mediated resistance. To identify and confirm potential mechanisms of QDR toward P. sojae, we conducted a time course study comparing changes in gene expression among Conrad and M92-220 with high QDR to susceptible genotypes, Sloan, and three mutants derived from fast neutron irradiation of M92-220. Differentially expressed genes from Conrad and M92-220 indicated several shared defense-related pathways at the transcriptomic level but also defense pathways unique to each cultivar, such as stilbenoid, diarylheptanoid, and gingerol biosynthesis and monobactam biosynthesis. Gene Ontology pathway analysis showed that the susceptible fast neutron mutants lacked enrichment of three terpenoid-related pathways and two cell wall-related pathways at either one or both time points, in contrast to M92-220. The susceptible mutants also lacked enrichment of potentially important Kyoto Encyclopedia of Genes and Genomes pathways at either one or both time points, including sesquiterpenoid and triterpenoid biosynthesis; thiamine metabolism; arachidonic acid; stilbenoid, diarylheptanoid, and gingerol biosynthesis; and monobactam biosynthesis. Additionally, 31 genes that were differentially expressed in M92-220 following P. sojae infection were not expressed in the mutants. These 31 genes have annotations related to unknown proteins; valine, leucine, and isoleucine biosynthesis; and protein and lipid metabolic processes. The results of this study confirm previously proposed mechanisms of QDR, provide evidence for potential novel QDR pathways in M92-220, and further our understanding of the complex network associated with QDR mechanisms in soybean toward P. sojae.

Molecular complexity of quantitative immunity in plants: from QTL mapping to functional and systems biology.

In nature, plants defend themselves against pathogen attack by activating an arsenal of defense mechanisms. During the last decades, work mainly focused on the understanding of qualitative disease resistance mediated by a few genes conferring an almost complete resistance, while quantitative disease resistance (QDR) remains poorly understood despite the fact that it represents the predominant and more durable form of resistance in natural populations and crops. Here, we review our past and present work on the dissection of the complex mechanisms underlying QDR in Arabidopsis thaliana. The strategies, main steps and challenges of our studies related to one atypical QDR gene, RKS1 (Resistance related KinaSe 1), are presented. First, from genetic analyses by QTL (Quantitative Trait Locus) mapping and GWAs (Genome Wide Association studies), the identification, cloning and functional analysis of this gene have been used as a starting point for the exploration of the multiple and coordinated pathways acting together to mount the QDR response dependent on RKS1. Identification of RKS1 protein interactors and complexes was a first step, systems biology and reconstruction of protein networks were then used to decipher the molecular roadmap to the immune responses controlled by RKS1. Finally, exploration of the potential impact of key components of the RKS1-dependent gene network on leaf microbiota offers interesting and challenging perspectives to decipher how the plant immune systems interact with the microbial communities' systems.

Dans la nature, les plantes se défendent contre les attaques pathogènes en activant tout un arsenal de mécanismes de défense. Au cours des décennies passées, la recherche s'est principalement focalisée sur la compréhension de la résistance qualitative médiée par quelques gènes majeurs conférant une résistance quasi complète, alors que la résistance quantitative (QDR) demeure peu comprise bien qu'elle représente la forme de résistance prédominante et la plus durable dans les populations naturelles ou les cultures. Nous donnons ici une revue de nos travaux passés et présents sur la dissection des mécanismes complexes qui sous-tendent la QDR chez Arabidopsis thaliana. Les stratégies, étapes clés et défis de nos études concernant un gène QDR atypique, RKS1 (Resistance related KinaSe 1), sont rapportés. En premier lieu, à partir d'analyses génétiques par cartographie de QTL et GWA, l'identification, le clonage et l'analyse fonctionnelle de ce gène ont été utilisés comme point de départ à l'exploration des voies multiples et coordonnées agissant ensemble pour le développement de la réponse QDR dépendante de RKS1. L'identification des interacteurs et complexes protéiques impliquant RKS1 a été une première étape, la biologie des systèmes et la reconstruction de réseaux d'interactions protéines-protéines ont ensuite été mises en œuvre pour décoder les voies moléculaires conduisant aux réponses immunitaires contrôlées par RKS1. Finalement, l'exploration de l'impact potentiel de composantes clés du réseau de gènes dépendant de RKS1 sur le microbiote, offre des perspectives intéressantes et ambitieuses pour comprendre comment le système immunitaire de la plante interagit avec le système des communautés microbiennes.

Inhibition of ethylene involved in resistance to E. turcicum in an exotic-derived double haploid maize population.

Northern corn leaf blight (NCLB) is an economically important disease of maize. While the genetic architecture of NCLB has been well characterized, the pathogen is known to overcome currently deployed resistance genes, and the role of hormones in resistance to NCLB is an area of active research. The objectives of the study were (i) to identify significant markers associated with resistance to NCLB, (ii) to identify metabolic pathways associated with NCLB resistance, and (iii) to examine role of ethylene in resistance to NCLB. We screened 252 lines from the exotic-derived double haploid BGEM maize population for resistance to NCLB in both field and greenhouse environments. We used a genome wide association study (GWAS) and stepwise regression to identify four markers associated with resistance, followed by a pathway association study tool (PAST) to identify important metabolic pathways associated with disease severity and incubation period. The ethylene synthesis pathway was significant for disease severity and incubation period. We conducted a greenhouse assay in which we inhibited ethylene to examine the role of ethylene in resistance to NCLB. We observed a significant increase in incubation period and a significant decrease in disease severity between plants treated with the ethylene inhibitor and mock-treated plants. Our study confirms the potential of the BGEM population as a source of novel alleles for resistance. We also confirm the role of ethylene in resistance to NCLB and contribute to the growing body of literature on ethylene and disease resistance in monocots.

Medicago truncatula quantitative resistance to a new strain of Verticillium alfalfae from Iran revealed by a genome-wide association study.

Verticillium wilt is a major threat to many crops, among them alfalfa (Medicago sativa). The model plant Medicago truncatula, a close relative of alfalfa was used to study the genetic control of resistance towards a new Verticillium alfalfae isolate. The accidental introduction of pathogen strains through global trade is a threat to crop production and such new strains might also be better adapted to global warming. Isolates of V. alfalfae were obtained from alfalfa fields in Iran and characterized. The Iranian isolate AF1 was used in a genome-wide association study (GWAS) involving 242 accessions from the Mediterranean region. Root inoculations were performed with conidia at 25°C and symptoms were scored regularly. Maximum Symptom Score and Area under Disease Progess Curve were computed as phenotypic traits to be used in GWAS and for comparison to a previous study with French isolate V31.2 at 20°C. This comparison showed high correlation with a shift to higher susceptibility, and similar geographical distribution of resistant and susceptible accessions to AF1 at 25°C, with resistant accessions mainly in the western part. GWAS revealed 30 significant SNPs linked to resistance towards isolate AF1. None of them were common to the previous study with isolate V31.2 at 20°C. To confirm these loci, the expression of nine underlying genes was studied. All genes were induced in roots following inoculation, in susceptible and resistant plants. However, in resistant plants induction was higher and lasted longer. Taken together, the use of a new pathogen strain and a shift in temperature revealed a completely different genetic control compared to a previous study that demonstrated the existence of two major QTLs. These results can be useful for Medicago breeding programs to obtain varieties better adapted to future conditions.

Association Mapping Combined with Whole Genome Sequencing Data Reveals Candidate Causal Variants for Sclerotinia Stem Rot Resistance in Brassica napus.

Canola (Brassica napus) yield can be significantly reduced by the disease sclerotinia stem rot (SSR), which is caused by Sclerotinia sclerotiorum, a necrotrophic fungal pathogen with an unusually large host range. Breeding cultivars that are physiologically resistant to SSR is desirable to enhance crop productivity. However, the development of resistant varieties has proved challenging due to the highly polygenic nature of S. sclerotiorum resistance. Here, we identified regions of the B. napus genome associated with SSR resistance using data from a previous study by association mapping. We then validated their contribution to resistance in a follow-up screen. This follow-up screen also confirmed high levels of SSR resistance in several genotypes from the previous study. Using publicly available whole genome sequencing data for a panel of 83 B. napus genotypes, we identified nonsynonymous polymorphisms linked to the SSR resistance loci. A qPCR analysis showed that two of the genes containing these polymorphisms were transcriptionally responsive to S. sclerotiorum infection. In addition, we provide evidence that homologues of three of the candidate genes contribute to resistance in the model Brassicaceae species Arabidopsis thaliana. The identification of resistant germplasm and candidate genomic loci associated with resistance are important findings that can be exploited by breeders to improve the genetic resistance of canola varieties.

Exploring the Relationship Between Pattern-Triggered Immunity and Quantitative Resistance to Xanthomonas vasicola pv. vasculorum in Maize.

Bacterial leaf streak (BLS) of maize is an emerging foliar disease of maize in the Americas. It is caused by the gram-negative nonvascular bacterium Xanthomonas vasicola pv. vasculorum. There are no chemical controls available for BLS, and thus, host resistance is crucial for managing X. vasicola pv. vasculorum. The objective of this study was to examine the genetic determinants of resistance to X. vasicola pv. vasculorum in maize, as well as the relationship between other defense-related traits and BLS resistance. Specifically, we examined the correlations among BLS severity, severity for three fungal diseases, flg-22 response, hypersensitive response, and auricle color. We conducted quantitative trait locus (QTL) mapping for X. vasicola pv. vasculorum resistance using the maize recombinant inbred line population Z003 (B73 × CML228). We detected three QTLs for BLS resistance. In addition to the disease resistance QTL, we detected a single QTL for auricle color. We observed significant, yet weak, correlations among BLS severity, levels of pattern-triggered immunity response and leaf flecking. These results will be useful for understanding resistance to X. vasicola pv. vasculorum and mitigating the impact of BLS on maize yields.

Overexpression of SQUALENE SYNTHASE Reduces Nicotiana benthamiana Resistance against Phytophthora infestans.

Plant triterpenoids play a critical role in plant resistance against Phytophthora infestans de Bary, the causal pathogen of potato and tomato late blight. However, different triterpenoids could have contrasting functions on plant resistance against P. infestans. In this study, we targeted the key biosynthetic gene of all plant triterpenoids, SQUALENE SYNTHASE (SQS), to examine the function of this gene in plant-P. infestans interactions. A post-inoculation, time-course gene expression analysis revealed that SQS expression was induced in Nicotiana benthamiana but was transiently suppressed in Solanum lycopersicum. Consistent with the host-specific changes in SQS expression, concentrations of major triterpenoid compounds were only induced in S. lycopersicum. A stable overexpression of SQS in N. benthamiana reduced plant resistance against P. infestans and induced the hyperaccumulation of stigmasterol. A comparative transcriptomics analysis of the transgenic lines showed that diverse plant physiological processes were influenced by SQS overexpression, suggesting that phytosterol content regulation may not be the sole mechanism through which SQS promotes plant susceptibility towards P. infestans. This study provides experimental evidence for the host-specific transcriptional regulation and function of SQS in plant interactions with P. infestans, offering a novel perspective in examining the quantitative disease resistance against late blight.

Tomato deploys defence and growth simultaneously to resist bacterial wilt disease.

Plant disease limits crop production, and host genetic resistance is a major means of control. Plant pathogenic Ralstonia causes bacterial wilt disease and is best controlled with resistant varieties. Tomato wilt resistance is multigenic, yet the mechanisms of resistance remain largely unknown. We combined metaRNAseq analysis and functional experiments to identify core Ralstonia-responsive genes and the corresponding biological mechanisms in wilt-resistant and wilt-susceptible tomatoes. While trade-offs between growth and defence are common in plants, wilt-resistant plants activated both defence responses and growth processes. Measurements of innate immunity and growth, including reactive oxygen species production and root system growth, respectively, validated that resistant plants executed defence-related processes at the same time they increased root growth. In contrast, in wilt-susceptible plants roots senesced and root surface area declined following Ralstonia inoculation. Wilt-resistant plants repressed genes predicted to negatively regulate water stress tolerance, while susceptible plants repressed genes predicted to promote water stress tolerance. Our results suggest that wilt-resistant plants can simultaneously promote growth and defence by investing in resources that act in both processes. Infected susceptible plants activate defences, but fail to grow and so succumb to Ralstonia, likely because they cannot tolerate the water stress induced by vascular wilt.

Quantitative disease resistance: Multifaceted players in plant defense.

In contrast to large-effect qualitative disease resistance, quantitative disease resistance (QDR) exhibits partial and generally durable resistance and has been extensively utilized in crop breeding. The molecular mechanisms underlying QDR remain largely unknown but considerable progress has been made in this area in recent years. In this review, we summarize the genes that have been associated with plant QDR and their biological functions. Many QDR genes belong to the canonical resistance gene categories with predicted functions in pathogen perception, signal transduction, phytohormone homeostasis, metabolite transport and biosynthesis, and epigenetic regulation. However, other "atypical" QDR genes are predicted to be involved in processes that are not commonly associated with disease resistance, such as vesicle trafficking, molecular chaperones, and others. This diversity of function for QDR genes contrasts with qualitative resistance, which is often based on the actions of nucleotide-binding leucine-rich repeat (NLR) resistance proteins. An understanding of the diversity of QDR mechanisms and of which mechanisms are effective against which classes of pathogens will enable the more effective deployment of QDR to produce more durably resistant, resilient crops.

Pathogen-derived mechanical cues potentiate the spatio-temporal implementation of plant defense.

BACKGROUND: The ongoing adaptation of plants to their environment is the basis for their survival. In this adaptation, mechanoperception of gravity and local curvature plays a role of prime importance in finely regulating growth and ensuring a dynamic balance preventing buckling. However, the abiotic environment is not the exclusive cause of mechanical stimuli. Biotic interactions between plants and microorganisms also involve physical forces and potentially mechanoperception. Whether pathogens trigger mechanoperception in plants and the impact of mechanotransduction on the regulation of plant defense remains however elusive. RESULTS: Here, we found that the perception of pathogen-derived mechanical cues by microtubules potentiates the spatio-temporal implementation of plant immunity to fungus. By combining biomechanics modeling and image analysis of the post-invasion stage, we reveal that fungal colonization releases plant cell wall-born tension locally, causing fluctuations of tensile stress in walls of healthy cells distant from the infection site. In healthy cells, the pathogen-derived mechanical cues guide the reorganization of mechanosensing cortical microtubules (CMT). The anisotropic patterning of CMTs is required for the regulation of immunity-related genes in distal cells. The CMT-mediated mechanotransduction of pathogen-derived cues increases Arabidopsis disease resistance by 40% when challenged with the fungus Sclerotinia sclerotiorum. CONCLUSIONS: CMT anisotropic patterning triggered by pathogen-derived mechanical cues activates the implementation of early plant defense in cells distant from the infection site. We propose that the mechano-signaling triggered immunity (MTI) complements the molecular signals involved in pattern and effector-triggered immunity.

Identification of Candidate Genes for a Major Quantitative Disease Resistance Locus From Soybean PI 427105B for Resistance to Phytophthora sojae.

Phytophthora root and stem rot is a yield-limiting soybean disease caused by the soil-borne oomycete Phytophthora sojae. Although multiple quantitative disease resistance loci (QDRL) have been identified, most explain <10% of the phenotypic variation (PV). The major QDRL explaining up to 45% of the PV were previously identified on chromosome 18 and represent a valuable source of resistance for soybean breeding programs. Resistance alleles from plant introductions 427105B and 427106 significantly increase yield in disease-prone fields and result in no significant yield difference in fields with less to no disease pressure. In this study, high-resolution mapping reduced the QDRL interval to 3.1 cm, and RNA-seq analysis of near-isogenic lines (NILs) varying at QDRL-18 pinpointed a single gene of interest which was downregulated in inoculated NILs carrying the resistant allele compared to inoculated NILs with the susceptible allele. This gene of interest putatively encodes a serine-threonine kinase (STK) related to the AtCR4 family and may be acting as a susceptibility factor, based on the specific increase of jasmonic acid concentration in inoculated NILs. This work facilitates further functional analyses and marker-assisted breeding efforts by prioritizing candidate genes and narrowing the targeted region for introgression.

The evolutionary and molecular features of the broad-host-range plant pathogen Sclerotinia sclerotiorum.

Sclerotinia sclerotiorum is a pathogenic fungus that infects hundreds of plant species, including many of the world's most important crops. Key features of S. sclerotiorum include its extraordinary host range, preference for dicotyledonous plants, relatively slow evolution, and production of protein effectors that are active in multiple host species. Plant resistance to this pathogen is highly complex, typically involving numerous polymorphisms with infinitesimally small effects, which makes resistance breeding a major challenge. Due to its economic significance, S. sclerotiorum has been subjected to a large amount of molecular and evolutionary research. In this updated pathogen profile, we review the evolutionary and molecular features of S. sclerotiorum and discuss avenues for future research into this important species.

Quantitative Resistance Loci to Southern Rust Mapped in a Temperate Maize Diversity Panel.

Southern rust is a severe foliar disease of maize (Zea mays) resulting from infection with the obligate biotrophic fungus Puccinia polysora. This disease reduces photosynthetic productivity, which in turn reduces yields, with the greatest yield losses (up to 50%) associated with earlier onset infections. P. polysora urediniospores overwinter only in tropical and subtropical regions but cause outbreaks when environmental conditions favor initial infection. Increased temperatures and humidity during the growing season combined with an increased frequency of moderate winters are likely to increase the frequency of severe southern rust outbreaks in the U.S. Corn Belt. In summer 2020, a severe outbreak of southern rust was observed in eastern Nebraska, United States. We scored a replicated maize association panel planted in Lincoln, NE for disease severity and found that disease incidence and severity showed significant variation among maize genotypes. Genome-wide association studies identified four loci associated with significant quantitative variation in disease severity. These loci were associated with candidate genes with plausible links to quantitative disease resistance. A transcriptome-wide association study identified additional genes associated with disease severity. Together, these results indicate that substantial diversity in resistance to southern rust exists among current temperate-adapted maize germplasm, including several candidate loci that may explain the observed variation in resistance to southern rust.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

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

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

Differential Regulation of Maize and Sorghum Orthologs in Response to the Fungal Pathogen Exserohilum turcicum.

Pathogens that infect more than one host offer an opportunity to study how resistance mechanisms have evolved across different species. Exserohilum turcicum infects both maize and sorghum and the isolates are host-specific, offering a unique system to examine both compatible and incompatible interactions. We conducted transcriptional analysis of maize and sorghum in response to maize-specific and sorghum-specific E. turcicum isolates and identified functionally related co-expressed modules. Maize had a more robust transcriptional response than sorghum. E. turcicum responsive genes were enriched in core orthologs in both crops, but only up to 16% of core orthologs showed conserved expression patterns. Most changes in gene expression for the core orthologs, including hub genes, were lineage specific, suggesting a role for regulatory divergent evolution. We identified several defense-related shared differentially expressed (DE) orthologs with conserved expression patterns between the two crops, suggesting a role for parallel evolution of those genes in both crops. Many of the differentially expressed genes (DEGs) during the incompatible interaction were related to quantitative disease resistance (QDR). This work offers insights into how different hosts with relatively recent divergence interact with a common pathogen. Our results are important for developing resistance to this critical pathogen and understanding the evolution of host-pathogen interactions.

Hormesis-Inducing Essential Oil Nanodelivery System Protects Plants against Broad Host-Range Necrotrophs.

Botrytis cinerea, a broad host-range necrotrophic (BHN) phytopathogen, establishes compatible interactions with hosts by deploying multigene infection strategies, rendering simply inherited resistance ineffective to fight off this pathogen. Since essential oils (EOs) serve as intermediators during phytobiome communication, we hypothesize that they have the potential to enhance the quantitative disease resistance against BHN by eliciting the adaptive stress response (hormesis) in plants. However, using EOs is challenging due to their poor solubility in water. Nanoemulsification of EOs enhances not only the solubility of EOs but also their potency and stability. Here, we demonstrate the potential use of essential oil nanoemulsions (EONEs) to control infections caused by BHN. Using basic engineering principles of nanocarrier design, we demonstrate the efficacy of a robust EONEs design for controlling B. cinerea infection in a model plant, Arabidopsis thaliana. Our nanoemulsion delivery system significantly enhanced the disease resistance of the host by reducing the necrotic area by up to 50% compared to untreated plants. RNA-seq analysis indicated that successful treatments upregulated autophagy, ROS scavenging, and activation of the jasmonic acid signaling pathway.

A gene-for-gene interaction involving a 'late' effector contributes to quantitative resistance to the stem canker disease in Brassica napus.

The control of stem canker disease of Brassica napus (rapeseed), caused by the fungus Leptosphaeria maculans is based largely on plant genetic resistance: single-gene specific resistance (Rlm genes) or quantitative, polygenic, adult-stage resistance. Our working hypothesis was that quantitative resistance partly obeys the gene-for-gene model, with resistance genes 'recognizing' fungal effectors expressed during late systemic colonization. Five LmSTEE (stem-expressed effector) genes were selected and placed under the control of the AvrLm4-7 promoter, an effector gene highly expressed at the cotyledon stage of infection, for miniaturized cotyledon inoculation test screening of a gene pool of 204 rapeseed genotypes. We identified a rapeseed genotype, 'Yudal', expressing hypersensitive response to LmSTEE98. The LmSTEE98-RlmSTEE98 interaction was further validated by inactivation of the LmSTEE98 gene with a CRISPR-Cas9 approach. Isolates with mutated versions of LmSTEE98 induced more severe stem symptoms than the wild-type isolate in 'Yudal'. This single-gene resistance was mapped in a 0.6 cM interval of the 'Darmor_bzh' × 'Yudal' genetic map. One typical gene-for-gene interaction contributes partly to quantitative resistance when L. maculans colonizes the stems of rapeseed. With numerous other effectors specific to stem colonization, our study provides a new route for resistance gene discovery, elucidation of quantitative resistance mechanisms and selection for durable resistance.

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.