An atypical kinase under balancing selection confers broad-spectrum disease resistance in Arabidopsis

A landscape of resistance gene analogs in sour cherry (Prunus cerasus L.)

OBJECTIVE

This research aims to analyze the presence and distribution of resistance genes in the avium and fruticosa subgenomes of Prunus cerasus through computational methods and bioinformatics tools.

RESULTS

Analysis of genome and transcriptome sequencing data revealed a total of 19,570 transcripts with at least one resistance gene domain in Prunus cerasus subgenome avium and 19,142 in Prunus cerasus subgenome fruticosa. Key findings include the identification of 804 "complete" resistance gene transcripts in Prunus cerasus subgenome avium and 817 in Prunus cerasus subgenome fruticosa, with distinct distributions of resistance gene classes observed between the subgenomes. Phylogenetic analysis showed clustering of resistance genes, and unique resistance proteins were identified in each subgenome. Functional annotation comparisons with Arabidopsis thaliana highlighted shared and unique resistance genes, emphasizing the complexity of disease resistance in cherry species. Additionally, a higher diversity of RLKs and RLPs was observed, with 504 transcripts identified and 18 showing similarity to known reference genes.

Natural genetic variation underlying the negative effect of elevated CO2 on ionome composition in Arabidopsis thaliana

The elevation of atmospheric CO2 leads to a decline in plant mineral content, which might pose a significant threat to food security in coming decades. Although few genes have been identified for the negative effect of elevated CO2 on plant mineral composition, several studies suggest the existence of genetic factors. Here, we performed a large-scale study to explore genetic diversity of plant ionome responses to elevated CO2, using six hundred Arabidopsis thaliana accessions, representing geographical distributions ranging from worldwide to regional and local environments. We show that growth under elevated CO2 leads to a global decrease of ionome content, whatever the geographic distribution of the population. We observed a high range of genetic diversity, ranging from the most negative effect to resilience or even to a benefit in response to elevated CO2. Using genome-wide association mapping, we identified a large set of genes associated with this response, and we demonstrated that the function of one of these genes is involved in the negative effect of elevated CO2 on plant mineral composition. This resource will contribute to understand the mechanisms underlying the effect of elevated CO2 on plant mineral nutrition, and could help towards the development of crops adapted to a high-CO2 world.

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.

Two adjacent NLR genes conferring quantitative resistance to clubroot disease in Arabidopsis are regulated by a stably inherited epiallelic variation

Clubroot caused by the protist Plasmodiophora brassicae is a major disease affecting cultivated Brassicaceae. Using a combination of quantitative trait locus (QTL) fine mapping, CRISPR-Cas9 validation, and extensive analyses of DNA sequence and methylation patterns, we revealed that the two adjacent neighboring NLR (nucleotide-binding and leucine-rich repeat) genes AT5G47260 and AT5G47280 cooperate in controlling broad-spectrum quantitative partial resistance to the root pathogen P. brassicae in Arabidopsis and that they are epigenetically regulated. The variation in DNA methylation is not associated with any nucleotide variation or any transposable element presence/absence variants and is stably inherited. Variations in DNA methylation at the Pb-At5.2 QTL are widespread across Arabidopsis accessions and correlate negatively with variations in expression of the two genes. Our study demonstrates that natural, stable, and transgenerationally inherited epigenetic variations can play an important role in shaping resistance to plant pathogens by modulating the expression of immune receptors.

Repeated loss of the ability of a wild pepper disease resistance gene to function at high temperatures suggests that thermoresistance is a costly trait

Specificity in plant-pathogen gene-for-gene (GFG) interactions is determined by the recognition of pathogen proteins by the products of plant resistance (R) genes. The evolutionary dynamics of R genes in plant-virus systems is poorly understood. We analyse the evolution of the L resistance locus to tobamoviruses in the wild pepper Capsicum annuum var. glabriusculum (chiltepin), a crop relative undergoing incipient domestication. The frequency, and the genetic and phenotypic diversity, of the L locus was analysed in 41 chiltepin populations under different levels of human management over its distribution range in Mexico. The frequency of resistance was lower in Cultivated than in Wild populations. L-locus genetic diversity showed a strong spatial structure with no isolation-by-distance pattern, suggesting environment-specific selection, possibly associated with infection by the highly virulent tobamoviruses found in the surveyed regions. L alleles differed in recognition specificity and in the expression of resistance at different temperatures, broad-spectrum recognition of P0  + P1 pathotypes and expression above 32°C being ancestral traits that were repeatedly lost along L-locus evolution. Overall, loss of resistance co-occurs with incipient domestication and broad-spectrum resistance expressed at high temperatures has apparent fitness costs. These findings contribute to understand the role of fitness trade-offs in plant-virus coevolution.

Genome scans capture key adaptation and historical hybridization signatures in tetraploid wheat

Tetraploid wheats (Triticum turgidum L.), including durum wheat (T. turgidum ssp. durum (Desf.) Husn.), are important crops with high nutritional and cultural values. However, their production is constrained by sensitivity to environmental conditions. In search of adaptive genetic signatures tracing historical selection and hybridization events, we performed genome scans on two datasets: (1) Durum Global Diversity Panel comprising a total of 442 tetraploid wheat and wild progenitor accessions including durum landraces (n = 286), domesticated emmer (T. turgidum ssp. dicoccum (Schrank) Thell.; n = 103) and wild emmer (T. turgidum ssp. dicoccoides (Korn. ex Asch. & Graebn.) Thell.; n = 53) wheats genotyped using the 90K single nucleotide polymorphism (SNP) array, and (2) a second dataset comprising a total 121 accessions of nine T. turgidum subspecies including wild emmer genotyped with >100 M SNPs from whole-genome resequencing. The genome scan on the first dataset detected six outlier loci on chromosomes 1A, 1B, 3A (n = 2), 6A, and 7A. These loci harbored important genes for adaptation to abiotic stresses, phenological responses, such as seed dormancy, circadian clock, flowering time, and key yield-related traits, including pleiotropic genes, such as HAT1, KUODA1, CBL1, and ZFN1. The scan on the second dataset captured a highly differentiated region on chromosome 2B that shows significant differentiation between two groups: one group consists of Georgian (T. turgidum ssp. paleocolchicum A. Love & D. Love) and Persian (T. turgidum ssp. carthlicum (Nevski) A. Love & D. Love) wheat accessions, while the other group comprises all the remaining tetraploids including wild emmer. This is consistent with a previously reported introgression in this genomic region from T. timopheevii Zhuk. which naturally cohabit in the Georgian and neighboring areas. This region harbored several adaptive genes, including the thermomorphogenesis gene PIF4, which confers temperature-resilient disease resistance and regulates other biological processes. Genome scans can be used to fast-track germplasm housed in gene banks and in situ; which helps to identify environmentally resilient accessions for breeding and/or to prioritize them for conservation.

Genetic architecture of the response of Arabidopsis thaliana to a native plant-growth-promoting bacterial strain

By improving plant nutrition and alleviating abiotic and biotic stresses, plant growth-promoting bacteria (PGPB) can help to develop eco-friendly and sustainable agricultural practices. Besides climatic conditions, soil conditions, and microbe-microbe interactions, the host genotype influences the effectiveness of PGPB. Yet, most GWAS conducted to characterize the genetic architecture of response to PGPB are based on non-native interactions between a host plant and PGPB strains isolated from the belowground compartment of other plants. In this study, a GWAS was set up under in vitro conditions to describe the genetic architecture of the response of Arabidopsis thaliana to the PGPB Pseudomonas siliginis, by inoculating seeds of 162 natural accessions from the southwest of France with one strain isolated from the leaf compartment in the same geographical region. Strong genetic variation of plant growth response to this native PGPB was observed at a regional scale, with the strain having a positive effect on the vegetative growth of small plants and a negative effect on the vegetative growth of large plants. The polygenic genetic architecture underlying this negative trade-off showed suggestive signatures of local adaptation. The main eco-evolutionary relevant candidate genes are involved in seed and root development.

An atypical NLR gene confers bacterial wilt susceptibility in Arabidopsis

Quantitative disease resistance (QDR) remains the most prevalent form of plant resistance in crop fields and wild habitats. Genome-wide association studies (GWAS) have proved to be successful in deciphering the quantitative genetic basis of complex traits such as QDR. To unravel the genetics of QDR to the devastating worldwide bacterial pathogen Ralstonia solanacearum, we performed a GWAS by challenging a highly polymorphic local mapping population of Arabidopsis thaliana with four R. solanacearum type III effector (T3E) mutants, identified as key pathogenicity determinants after a first screen on an A. thaliana core collection of 25 accessions. Although most quantitative trait loci (QTLs) were highly specific to the identity of the T3E mutant (ripAC, ripAG, ripAQ, and ripU), we finely mapped a common QTL located on a cluster of nucleotide-binding domain and leucine-rich repeat (NLR) genes that exhibited structural variation. We functionally validated one of these NLRs as a susceptibility factor in response to R. solanacearum, named it Bacterial Wilt Susceptibility 1 (BWS1), and cloned two alleles that conferred contrasting levels of QDR. Further characterization indicated that expression of BWS1 leads to suppression of immunity triggered by different R. solanacearum effectors. In addition, we showed a direct interaction between BWS1 and RipAC T3E, and BWS1 and SUPPRESSOR OF G2 ALLELE OF skp1 (SGT1b), the latter interaction being suppressed by RipAC. Together, our results highlight a putative role for BWS1 as a quantitative susceptibility factor directly targeted by the T3E RipAC, mediating negative regulation of the SGT1-dependent immune response.

Genome-wide association studies in plant pathosystems: success or failure?

Harnessing natural genetic variation is an established alternative to artificial genetic variation for investigating the molecular dialog between partners in plant pathosystems. Herein, we review the successes of genome-wide association studies (GWAS) in both plants and pathogens. While GWAS in plants confirmed that the genetic architecture of disease resistance is polygenic, dynamic during the infection kinetics, and dependent on the environment, GWAS shortened the time of identification of quantitative trait loci (QTLs) and revealed both complex epistatic networks and a genetic architecture dependent upon the geographical scale. A similar picture emerges from the few GWAS in pathogens. In addition, the ever-increasing number of functionally validated QTLs has revealed new molecular plant defense mechanisms and pathogenicity determinants. Finally, we propose recommendations to better decode the disease triangle.

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.

A Genome-Wide Association study in Arabidopsis thaliana to decipher the adaptive genetics of quantitative disease resistance in a native heterogeneous environment

Pathogens are often the main selective agents acting in plant communities, thereby influencing the distribution of polymorphism at loci affecting resistance within and among natural plant populations. In addition, the outcome of plant-pathogen interactions can be drastically affected by abiotic and biotic factors at different spatial and temporal grains. The characterization of the adaptive genetic architecture of disease resistance in native heterogeneous environments is however still missing. In this study, we conducted an in situ Genome-Wide Association study in the spatially heterogeneous native habitat of a highly genetically polymorphic local mapping population of Arabidopsis thaliana, to unravel the adaptive genetic architecture of quantitative disease resistance. Disease resistance largely differed among three native soils and was affected by the presence of the grass Poa annua. The observation of strong crossing reactions norms among the 195 A. thaliana genotypes for disease resistance among micro-habitats, combined with a negative fecundity-disease resistance relationship in each micro-habitat, suggest that alternative local genotypes of A. thaliana are favored under contrasting environmental conditions at the scale of few meters. A complex genetic architecture was detected for disease resistance and fecundity. However, only few QTLs were common between these two traits. Heterogeneous selection in this local population should therefore promote the maintenance of polymorphism at only few candidate resistance genes.

ZAR1: Guardian of plant kinases

A key facet of innate immunity in plants entails the recognition of pathogen "effector" virulence proteins by host Nucleotide-Binding Leucine-Rich Repeat Receptors (NLRs). Among characterized NLRs, the broadly conserved ZAR1 NLR is particularly remarkable due to its capacity to recognize at least six distinct families of effectors from at least two bacterial genera. This expanded recognition spectrum is conferred through interactions between ZAR1 and a dynamic network of two families of Receptor-Like Cytoplasmic Kinases (RLCKs): ZED1-Related Kinases (ZRKs) and PBS1-Like Kinases (PBLs). In this review, we survey the history of functional studies on ZAR1, with an emphasis on how the ZAR1-RLCK network functions to trap diverse effectors. We discuss 1) the dynamics of the ZAR1-associated RLCK network; 2) the specificity between ZRKs and PBLs; and 3) the specificity between effectors and the RLCK network. We posit that the shared protein fold of kinases and the switch-like properties of their interactions make them ideal effector sensors, enabling ZAR1 to act as a broad spectrum guardian of host kinases.

Study of natural diversity in response to a key pathogenicity regulator of Ralstonia solanacearum reveals new susceptibility genes in Arabidopsis thaliana

Ralstonia solanacearum gram-negative phytopathogenic bacterium exerts its virulence through a type III secretion system (T3SS) that translocates type III effectors (T3Es) directly into the host cells. T3E secretion is finely controlled at the posttranslational level by helper proteins, T3SS control proteins, and type III chaperones. The HpaP protein, one of the type III secretion substrate specificity switch (T3S4) proteins, was previously highlighted as a virulence factor on Arabidopsis thaliana Col-0 accession. In this study, we set up a genome-wide association analysis to explore the natural diversity of response to the hpaP mutant of two A. thaliana mapping populations: a worldwide collection and a local population. Quantitative genetic variation revealed different genetic architectures in both mapping populations, with a global delayed response to the hpaP mutant compared to the GMI1000 wild-type strain. We have identified several quantitative trait loci (QTLs) associated with the hpaP mutant inoculation. The genes underlying these QTLs are involved in different and specific biological processes, some of which were demonstrated important for R. solanacearum virulence. We focused our study on four candidate genes, RKL1, IRE3, RACK1B, and PEX3, identified using the worldwide collection, and validated three of them as susceptibility factors. Our findings demonstrate that the study of the natural diversity of plant response to a R. solanacearum mutant in a key regulator of virulence is an original and powerful strategy to identify genes directly or indirectly targeted by the pathogen.

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.

Advances in Multi-Omics Approaches for Molecular Breeding of Black Rot Resistance in Brassica oleracea L

Brassica oleracea is one of the most important species of the Brassicaceae family encompassing several economically important vegetables produced and consumed worldwide. But its sustainability is challenged by a range of pathogens, among which black rot, caused by Xanthomonas campestris pv. campestris (Xcc), is the most serious and destructive seed borne bacterial disease, causing huge yield losses. Host-plant resistance could act as the most effective and efficient solution to curb black rot disease for sustainable production of B. oleracea. Recently, 'omics' technologies have emerged as promising tools to understand the host-pathogen interactions, thereby gaining a deeper insight into the resistance mechanisms. In this review, we have summarized the recent achievements made in the emerging omics technologies to tackle the black rot challenge in B. oleracea. With an integrated approach of the omics technologies such as genomics, proteomics, transcriptomics, and metabolomics, it would allow better understanding of the complex molecular mechanisms underlying black rot resistance. Due to the availability of sequencing data, genomics and transcriptomics have progressed as expected for black rot resistance, however, other omics approaches like proteomics and metabolomics are lagging behind, necessitating a holistic and targeted approach to address the complex questions of Xcc-Brassica interactions. Genomic studies revealed that the black rot resistance is a complex trait and is mostly controlled by quantitative trait locus (QTL) with minor effects. Transcriptomic analysis divulged the genes related to photosynthesis, glucosinolate biosynthesis and catabolism, phenylpropanoid biosynthesis pathway, ROS scavenging, calcium signalling, hormonal synthesis and signalling pathway are being differentially expressed upon Xcc infection. Comparative proteomic analysis in relation to susceptible and/or resistance interactions with Xcc identified the involvement of proteins related to photosynthesis, protein biosynthesis, processing and degradation, energy metabolism, innate immunity, redox homeostasis, and defence response and signalling pathways in Xcc-Brassica interaction. Specifically, most of the studies focused on the regulation of the photosynthesis-related proteins as a resistance response in both early and later stages of infection. Metabolomic studies suggested that glucosinolates (GSLs), especially aliphatic and indolic GSLs, its subsequent hydrolysis products, and defensive metabolites synthesized by jasmonic acid (JA)-mediated phenylpropanoid biosynthesis pathway are involved in disease resistance mechanisms against Xcc in Brassica species. Multi-omics analysis showed that JA signalling pathway is regulating resistance against hemibiotrophic pathogen like Xcc. So, the bonhomie between omics technologies and plant breeding is going to trigger major breakthroughs in the field of crop improvement by developing superior cultivars with broad-spectrum resistance. If multi-omics tools are implemented at the right scale, we may be able to achieve the maximum benefits from the minimum. In this review, we have also discussed the challenges, future prospects, and the way forward in the application of omics technologies to accelerate the breeding of B. oleracea for disease resistance. A deeper insight about the current knowledge on omics can offer promising results in the breeding of high-quality disease-resistant crops.

The Genomic Architecture of Competitive Response of Arabidopsis thaliana Is Highly Flexible Among Plurispecific Neighborhoods

Plants are daily challenged by multiple abiotic and biotic stresses. A major biotic constraint corresponds to competition with other plant species. Although plants simultaneously interact with multiple neighboring species throughout their life cycle, there is still very limited information about the genetics of the competitive response in the context of plurispecific interactions. Using a local mapping population of Arabidopsis thaliana, we set up a genome wide association study (GWAS) to estimate the extent of genetic variation of competitive response in 12 plant species assemblages, based on three competitor species (Poa annua, Stellaria media, and Veronica arvensis). Based on five phenotypic traits, we detected strong crossing reaction norms not only between the three bispecific neighborhoods but also among the plurispecific neighborhoods. The genetic architecture of competitive response was highly dependent on the identity and the relative abundance of the neighboring species. In addition, most of the enriched biological processes underlying competitive responses largely differ among neighborhoods. While the RNA related processes might confer a broad range response toolkit for multiple traits in diverse neighborhoods, some processes, such as signaling and transport, might play a specific role in particular assemblages. Altogether, our results suggest that plants can integrate and respond to different species assemblages depending on the identity and number of each neighboring species, through a large range of candidate genes associated with diverse and unexpected processes leading to developmental and stress responses.

A tomato receptor-like cytoplasmic kinase, SlZRK1, acts as a negative regulator in wound-induced jasmonic acid accumulation and insect resistance

Jasmonates accumulate rapidly and act as key regulators in response to mechanical wounding, but few studies have linked receptor-like cytoplasmic kinases (RLCKs) to wound-induced jasmonic acid (JA) signaling cascades. Here, we identified a novel wounding-induced RLCK-XII-2 subfamily member (SlZRK1) in tomato (Solanum lycopersicum) that was closely related to Arabidopsis HOPZ-ETI-DEFICIENT 1 (ZED1)-related kinases 1 based on phylogenetic analysis. SlZRK1 was targeted to the plasma membrane of tobacco mesophyll protoplasts as determined by transient co-expression with the plasma membrane marker mCherry-H+-ATPase. Catalytic residue sequence analysis and an in vitro kinase assay indicated that SlZRK1 may act as a pseudokinase. To further analyse the function of SlZRK1, we developed two stable knock-out mutants by CRISPR/Cas9. Loss of SlZRK1 significantly altered the expression of genes involved in JA biosynthesis, salicylic acid biosynthesis, and ethylene response. Furthermore, after mechanical wounding treatment, slzrk1 mutants increased transcription of early wound-inducible genes involved in JA biosynthesis and signaling. In addition, JA accumulation after wounding and plant resistance to herbivorous insects also were enhanced. Our findings expand plant regulatory networks in the wound-induced JA production by adding RLCKs as a new component in the wound signal transduction pathway.

One Hundred Years of Hybrid Necrosis: Hybrid Autoimmunity as a Window into the Mechanisms and Evolution of Plant-Pathogen Interactions

Hybrid necrosis in plants refers to a genetic autoimmunity syndrome in the progeny of interspecific or intraspecific crosses. Although the phenomenon was first documented in 1920, it has been unequivocally linked to autoimmunity only recently, with the discovery of the underlying genetic and biochemical mechanisms. The most common causal loci encode immune receptors, which are known to differ within and between species. One mechanism can be explained by the guard hypothesis, in which a guard protein, often a nucleotide-binding site-leucine-rich repeat protein, is activated by interaction with a plant protein that mimics standard guardees modified by pathogen effector proteins. Another surprising mechanism is the formation of inappropriately active immune receptor complexes. In this review, we summarize our current knowledge of hybrid necrosis and discuss how its study is not only informing the understanding of immune gene evolution but also revealing new aspects of plant immune signaling.

Genome-wide identification and evolutionary analysis of RLKs involved in the response to aluminium stress in peanut

BACKGROUND

As an important cash crop, the yield of peanut is influenced by soil acidification and pathogen infection. Receptor-like protein kinases play important roles in plant growth, development and stress responses. However, little is known about the number, location, structure, molecular phylogeny, and expression of RLKs in peanut, and no comprehensive analysis of RLKs in the Al stress response in peanuts have been reported.

RESULTS

A total of 1311 AhRLKs were identified from the peanut genome. The AhLRR-RLKs and AhLecRLKs were further divided into 24 and 35 subfamilies, respectively. The AhRLKs were randomly distributed across all 20 chromosomes in the peanut. Among these AhRLKs, 9.53% and 61.78% originated from tandem duplications and segmental duplications, respectively. The ka/ks ratios of 96.97% (96/99) of tandem duplication gene pairs and 98.78% (646/654) of segmental duplication gene pairs were less than 1. Among the tested tandem duplication clusters, there were 28 gene conversion events. Moreover, all total of 90 Al-responsive AhRLKs were identified by mining transcriptome data, and they were divided into 7 groups. Most of the Al-responsive AhRLKs that clustered together had similar motifs and evolutionarily conserved structures. The gene expression patterns of these genes in different tissues were further analysed, and tissue-specifically expressed genes, including 14 root-specific Al-responsive AhRLKs were found. In addition, all 90 Al-responsive AhRLKs which were distributed unevenly in the subfamilies of AhRLKs, showed different expression patterns between the two peanut varieties (Al-sensitive and Al-tolerant) under Al stress.

CONCLUSIONS

In this study, we analysed the RLK gene family in the peanut genome. Segmental duplication events were the main driving force for AhRLK evolution, and most AhRLKs subject to purifying selection. A total of 90 genes were identified as Al-responsive AhRLKs, and the classification, conserved motifs, structures, tissue expression patterns and predicted functions of Al-responsive AhRLKs were further analysed and discussed, revealing their putative roles. This study provides a better understanding of the structures and functions of AhRLKs and Al-responsive AhRLKs.

A New Catalog of Structural Variants in 1,301 A. thaliana Lines from Africa, Eurasia, and North America Reveals a Signature of Balancing Selection at Defense Response Genes

Genomic variation in the model plant Arabidopsis thaliana has been extensively used to understand evolutionary processes in natural populations, mainly focusing on single-nucleotide polymorphisms. Conversely, structural variation has been largely ignored in spite of its potential to dramatically affect phenotype. Here, we identify 155,440 indels and structural variants ranging in size from 1 bp to 10 kb, including presence/absence variants (PAVs), inversions, and tandem duplications in 1,301 A. thaliana natural accessions from Morocco, Madeira, Europe, Asia, and North America. We show evidence for strong purifying selection on PAVs in genes, in particular for housekeeping genes and homeobox genes, and we find that PAVs are concentrated in defense-related genes (R-genes, secondary metabolites) and F-box genes. This implies the presence of a "core" genome underlying basic cellular processes and a "flexible" genome that includes genes that may be important in spatially or temporally varying selection. Further, we find an excess of intermediate frequency PAVs in defense response genes in nearly all populations studied, consistent with a history of balancing selection on this class of genes. Finally, we find that PAVs in genes involved in the cold requirement for flowering (vernalization) and drought response are strongly associated with temperature at the sites of origin.

A genome-wide association study reveals cytokinin as a major component in the root defense responses against Ralstonia solanacearum

Bacterial wilt caused by the soil-borne pathogen Ralstonia solancearum is economically devastating, with no effective methods to fight the disease. This pathogen invades plants through their roots and colonizes their xylem, clogging the vasculature and causing rapid wilting. Key to preventing colonization are the early defense responses triggered in the host's root upon infection, which remain mostly unknown. Here, we have taken advantage of a high-throughput in vitro infection system to screen natural variability associated with the root growth inhibition phenotype caused by R. solanacearum in Arabidopsis during the first hours of infection. To analyze the genetic determinants of this trait, we have performed a genome-wide association study, identifying allelic variation at several loci related to cytokinin metabolism, including genes responsible for biosynthesis and degradation of cytokinin. Further, our data clearly demonstrate that cytokinin signaling is induced early during the infection process and cytokinin contributes to immunity against R. solanacearum. This study highlights a new role for cytokinin in root immunity, paving the way for future research that will help in understanding the mechanisms underpinning root defenses.

Arabidopsis thaliana Genes Associated with Cucumber mosaic virus Virulence and Their Link to Virus Seed Transmission

Virulence, the effect of pathogen infection on progeny production, is a major determinant of host and pathogen fitness as it affects host fecundity and pathogen transmission. In plant-virus interactions, ample evidence indicates that virulence is genetically controlled by both partners. However, the host genetic determinants are poorly understood. Through a genome-wide association study (GWAS) of 154 Arabidopsis thaliana genotypes infected by Cucumber mosaic virus (CMV), we identified eight host genes associated with virulence, most of them involved in response to biotic stresses and in cell wall biogenesis in plant reproductive structures. Given that virulence is a main determinant of the efficiency of plant virus seed transmission, we explored the link between this trait and the genetic regulation of virulence. Our results suggest that the same functions that control virulence are also important for CMV transmission through seeds. In sum, this work provides evidence of a novel role for some previously known plant defense genes and for the cell wall metabolism in plant virus interactions.

Arabidopsis natural variation in insect egg-induced cell death reveals a role for LECTIN RECEPTOR KINASE-I.1

In Arabidopsis (Arabidopsis thaliana), a hypersensitive-like response (HR-like response) is triggered underneath the eggs of the large white butterfly Pieris brassicae (P. brassicae), and this response is dependent on salicylic acid (SA) accumulation and signaling. Previous reports indicate that the clade I L-type LECTIN RECEPTOR KINASE-I.8 (LecRK-I.8) is involved in early steps of egg recognition. A genome-wide association study was used to better characterize the genetic structure of the HR-like response and discover loci that contribute to this response. We report here the identification of LecRK-I.1, a close homolog of LecRK-I.8, and show that two main haplotypes that explain part of the variation in HR-like response segregate among natural Arabidopsis accessions. Besides, signatures of balancing selection at this locus suggest that it may be ecologically important. Disruption of LecRK-I.1 results in decreased HR-like response and SA signaling, indicating that this protein is important for the observed responses. Furthermore, we provide evidence that LecRK-I.1 functions in the same signaling pathway as LecRK-I.8. Altogether, our results show that the response to eggs of P. brassicae is controlled by multiple LecRKs.

A natural diversity screen in Arabidopsis thaliana reveals determinants for HopZ1a recognition in the ZAR1-ZED1 immune complex

Pathogen pressure on hosts can lead to the evolution of genes regulating the innate immune response. By characterizing naturally occurring polymorphisms in immune receptors, we can better understand the molecular determinants of pathogen recognition. ZAR1 is an ancient Arabidopsis thaliana NLR (Nucleotide-binding [NB] Leucine-rich-repeat [LRR] Receptor) that recognizes multiple secreted effector proteins from the pathogenic bacteria Pseudomonas syringae and Xanthomonas campestris through its interaction with receptor-like cytoplasmic kinases (RLCKs). ZAR1 was first identified for its role in recognizing P. syringae effector HopZ1a, through its interaction with the RLCK ZED1. To identify additional determinants of HopZ1a recognition, we performed a computational screen for ecotypes from the 1001 Genomes project that were likely to lack HopZ1a recognition, and tested ~300 ecotypes. We identified ecotypes containing polymorphisms in ZAR1 and ZED1. Using our previously established Nicotiana benthamiana transient assay and Arabidopsis ecotypes, we tested for the effect of naturally occurring polymorphisms on ZAR1 interactions and the immune response. We identified key residues in the NB or LRR domain of ZAR1 that impact the interaction with ZED1. We demonstrate that natural diversity combined with functional assays can help define the molecular determinants and interactions necessary to regulate immune induction in response to pathogens.

A complex network of additive and epistatic quantitative trait loci underlies natural variation of Arabidopsis thaliana quantitative disease resistance to Ralstonia solanacearum under heat stress

Plant immunity is often negatively impacted by heat stress. However, the underlying molecular mechanisms remain poorly characterized. Based on a genome-wide association mapping approach, this study aims to identify in Arabidopsis thaliana the genetic bases of robust resistance mechanisms to the devastating pathogen Ralstonia solanacearum under heat stress. A local mapping population was phenotyped against the R. solanacearum GMI1000 strain at 27 and 30 °C. To obtain a precise description of the genetic architecture underlying natural variation of quantitative disease resistance (QDR), we applied a genome-wide local score analysis. Alongside an extensive genetic variation found in this local population at both temperatures, we observed a playful dynamics of quantitative trait loci along the infection stages. In addition, a complex genetic network of interacting loci could be detected at 30 °C. As a first step to investigate the underlying molecular mechanisms, the atypical meiotic cyclin SOLO DANCERS gene was validated by a reverse genetic approach as involved in QDR to R. solanacearum at 30 °C. In the context of climate change, the complex genetic architecture underlying QDR under heat stress in a local mapping population revealed candidate genes with diverse molecular functions.

What natural variation can teach us about resistance durability

Breeding a crop variety to be resistant to a pathogen usually takes years. This is problematic because pathogens, with short generation times and fluid genomes, adapt quickly to overcome resistance. The triumph of the pathogen is not inevitable, however, as there are numerous examples of durable resistance, particularly in wild plants. Which factors then contribute to such resistance stability over millennia? We review current knowledge of wild and agricultural pathosystems, detailing the importance of genetic, species and spatial heterogeneity in the prevention of pathogen outbreaks. We also highlight challenges associated with increasing resistance diversity in crops, both in light of pathogen (co-)evolution and breeding practices. Historically it has been difficult to incorporate heterogeneity into agriculture due to reduced efficiency in harvesting. Recent advances implementing computer vision and automation in agricultural production may improve our ability to harvest mixed genotype and mixed species plantings, thereby increasing resistance durability.

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.

Rapid identification of an Arabidopsis NLR gene as a candidate conferring susceptibility to Sclerotinia sclerotiorum using time-resolved automated phenotyping

The broad host range necrotrophic fungus Sclerotinia sclerotiorum is a devastating pathogen of many oil and vegetable crops. Plant genes conferring complete resistance against S. sclerotiorum have not been reported. Instead, plant populations challenged by S. sclerotiorum exhibit a continuum of partial resistance designated as quantitative disease resistance (QDR). Because of their complex interplay and their small phenotypic effect, the functional characterization of QDR genes remains limited. How broad host range necrotrophic fungi manipulate plant programmed cell death is for instance largely unknown. Here, we designed a time-resolved automated disease phenotyping pipeline enabling high-throughput disease lesion measurement with high resolution, low footprint at low cost. We could accurately recover contrasted disease responses in several pathosystems using this system. We used our phenotyping pipeline to assess the kinetics of disease symptoms caused by seven S. sclerotiorum isolates on six A. thaliana natural accessions with unprecedented resolution. Large effect polymorphisms common to the most resistant A. thaliana accessions identified highly divergent alleles of the nucleotide-binding site leucine-rich repeat gene LAZ5 in the resistant accessions Rubezhnoe and Lip-0. We show that impaired LAZ5 expression in laz5.1 mutant lines and in A. thaliana Rub natural accession correlate with enhanced QDR to S. sclerotiorum. These findings illustrate the value of time-resolved image-based phenotyping for unravelling the genetic bases of complex traits such as QDR. Our results suggest that S. sclerotiorum manipulates plant sphingolipid pathways guarded by LAZ5 to trigger programmed cell death and cause disease.

A High-Throughput, Seedling Screen for Plant Immunity

An understanding of how biological diversity affects plant-microbe interactions is becoming increasingly important, particularly with respect to components of the pathogen effector arsenal and the plant immune system. Although technological improvements have greatly advanced our ability to examine molecular sequences and interactions, relatively few advances have been made that facilitate high-throughput, in vivo pathology screens. Here, we present a high-throughput, microplate-based, nondestructive seedling pathology assay, and apply it to identify Arabidopsis thaliana effector-triggered immunity (ETI) responses against Pseudomonas syringae type III secreted effectors. The assay was carried out in a 48-well microplate format with spray inoculation, and disease symptoms were quantitatively recorded in a semiautomated manner, thereby greatly reducing both time and costs. The assay requires only slight modifications of common labware and uses no proprietary software. We validated the assay by recapitulating known ETI responses induced by P. syringae in Arabidopsis. We also demonstrated that we can quantitatively differentiate responses from a diversity of plant genotypes grown in the same microplate. Finally, we showed that the results obtained from our assay can be used to perform genome-wide association studies to identify host immunity genes, recapitulating results that have been independently obtained with mature plants.

HpaP Sequesters HrpJ, an Essential Component of Ralstonia solanacearum Virulence That Triggers Necrosis in Arabidopsis

The Gram-negative bacterium Ralstonia solanacearum, the causal agent of bacterial wilt, is a worldwide major crop pathogen whose virulence strongly relies on a type III secretion system (T3SS). This extracellular apparatus allows the translocation of proteins, called type III effectors (T3Es), directly into the host cells. To date, very few data are available in plant-pathogenic bacteria concerning the role played by type III secretion (T3S) regulators at the posttranslational level. We have demonstrated that HpaP, a putative T3S substrate specificity switch protein of R. solanacearum, controls T3E secretion. To better understand the role of HpaP on T3S control, we analyzed the secretomes of the GMI1000 wild-type strain as well as the hpaP mutant using a mass spectrometry experiment (liquid chromatography tandem mass spectrometry). The secretomes of both strains appeared to be very similar and highlighted the modulation of the secretion of few type III substrates. Interestingly, only one type III-associated protein, HrpJ, was identified as specifically secreted by the hpaP mutant. HrpJ appeared to be an essential component of the T3SS, essential for T3S and pathogenicity. We further showed that HrpJ is specifically translocated in planta by the hpaP mutant and that HrpJ can physically interact with HpaP. Moreover, confocal microscopy experiments demonstrated a cytoplasmic localization for HrpJ once in planta. When injected into Arabidopsis thaliana leaves, HrpJ is able to trigger a necrosis on 16 natural accessions. A genome-wide association mapping revealed a major association peak with 12 highly significant single-nucleotide polymorphisms located on a plant acyl-transferase.

Immunodiversity of the Arabidopsis ZAR1 NLR Is Conveyed by Receptor-Like Cytoplasmic Kinase Sensors

The Arabidopsis nucleotide-binding leucine-rich repeat protein ZAR1 can recognize at least six distinct families of pathogenic effector proteins to mount an effector-triggered immune response. This remarkable immunodiversity appears to be conveyed by receptor-like cytoplasmic kinase (RLCK) complexes, which associate with ZAR1 to sense several effector-induced kinase perturbations. Here we show that the recently identified ZAR1-mediated immune responses against the HopX1, HopO1, and HopBA1 effector families of Pseudomonas syringae rely on an expanded diversity of RLCK sensors. We show that individual sensors can recognize distinct effector families, thereby contributing to the expanded surveillance potential of ZAR1 and supporting its role as a guardian of the plant kinome.

Introgression Breeding in Cowpea [Vigna unguiculata (L.) Walp.]

The narrow base of genetic diversity characteristic of cowpea can be attributed to it being self-pollinating, evolving from narrow wild germplasm and exhibiting very limited gene flow between wild and cultivated types. Backcrossing to introduce simply inherited desirable traits and utilization of improved breeding lines and varieties as parents in crossing programs further narrowed the genetic base of cowpea varieties. In most cowpea breeding programs, genes for resistance and market traits were pyramided into lines characterized by high levels of acceptance to farmers and consumers. Besides predisposing widely distributed improved varieties to genetic vulnerability, a narrow base of genetic variation may be contributing to the plateauing in cowpea grain yield, which compromises genetic gains. Cross compatible wild relatives have not been used in variety development because breeders shy away from them due to their tiny seed size, unattractive seed coat color and texture, pod shattering, and susceptibility to viruses. A number of wild cowpea relatives, both within and outside section Catiang of Vigna species, have been evaluated for their reaction to cowpea insect pests and diseases. Vigna vexillata lines were resistant to the legume pod borer (Maruca vitrata), the cowpea weevil (Callosobruchus maculatus), and Striga gesnerioides but are cross incompatible with cultivated cowpea. Some lines among the cross compatible wild relative V. unguiculata ssp. dekindtiana were found to be resistant to aphid in the seedling stage, while others showed good levels of drought and heat tolerance. Molecular markers are being generated to identify quantitative trait loci (QTL) with effects on some desirable attributes in cowpea. Modern breeding tools, including transgenics, can be applied for the improvement of cowpea, bypassing the natural barriers of traditional breeding. Transgenic cowpea with Bt gene cry1Ab showing resistance to M. vitrata has been released in Nigeria. Genome editing, a powerful emerging tool, can also be used for developing improved cowpea varieties with durable resistance to pests and diseases.

Small RNAs from the plant pathogenic fungus Sclerotinia sclerotiorum highlight host candidate genes associated with quantitative disease resistance

Fungal plant pathogens secrete effector proteins and metabolites to cause disease. Additionally, some species transfer small RNAs (sRNAs) into plant cells to silence host mRNAs through complementary base pairing and suppress plant immunity. The fungus Sclerotinia sclerotiorum infects over 600 plant species, but little is known about the molecular processes that govern interactions with its many hosts. In particular, evidence for the production of sRNAs by S. sclerotiorum during infection is lacking. We sequenced sRNAs produced by S. sclerotiorum in vitro and during infection of two host species, Arabidopsis thaliana and Phaseolus vulgaris. We found that S. sclerotiorum produces at least 374 distinct highly abundant sRNAs during infection, mostly originating from repeat-rich plastic genomic regions. We predicted the targets of these sRNAs in A. thaliana and found that these genes were significantly more down-regulated during infection than the rest of the genome. Predicted targets of S. sclerotiorum sRNAs in A. thaliana were enriched for functional domains associated with plant immunity and were more strongly associated with quantitative disease resistance in a genome-wide association study (GWAS) than the rest of the genome. Mutants in A. thaliana predicted sRNA target genes SERK2 and SNAK2 were more susceptible to S. sclerotiorum than wild-type, suggesting that S. sclerotiorum sRNAs may contribute to the silencing of immune components in plants. The prediction of fungal sRNA targets in plant genomes can be combined with other global approaches, such as GWAS, to assist in the identification of plant genes involved in quantitative disease resistance.

Novel organ-specific genetic factors for quantitative resistance to late blight in potato

Potato, Solanum tuberosum, is one of the major consumed food in the world, being the basis of the diet of millions of people. The main limiting and destructive disease of potato is late blight, caused by Phytophtora infestans. Here, we present a multi-environmental analysis of the response to P. infestans using an association panel of 150 accessions of S. tuberosum Group Phureja, evaluated in two localities in Colombia. Disease resistance data were merged with a genotyping matrix of 83,862 SNPs obtained by 2b-restriction site–associated DNA and Genotyping by sequencing approaches into a Genome-wide association study. We are reporting 16 organ-specific QTL conferring resistance to late blight. These QTL explain between 13.7% and 50.9% of the phenotypic variance. Six and ten QTL were detected for resistance response in leaves and stem, respectively. In silico analysis revealed 15 candidate genes for resistance to late blight. Four of them have no functional genome annotation, while eleven candidate genes code for diverse proteins, including a leucine-rich repeat kinase.

WhoGEM: an admixture-based prediction machine accurately predicts quantitative functional traits in plants

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis. The method's prediction reliability equals or outperforms all existing algorithms for quantitative phenotype prediction. WhoGEM analysis produces evidence that variation in genome admixture proportions explains most of the phenotypic variation for quantitative phenotypes.

Cucumber mosaic virus infection as a potential selective pressure on Arabidopsis thaliana populations

It has been proposed that in wild ecosystems viruses are often plant mutualists, whereas agroecosystems favour pathogenicity. We seek evidence for virus pathogenicity in wild ecosystems through the analysis of plant-virus coevolution, which requires a negative effect of infection on the host fitness. We focus on the interaction between Arabidopsis thaliana and Cucumber mosaic virus (CMV), which is significant in nature. We studied the genetic diversity of A. thaliana for two defence traits, resistance and tolerance, to CMV. A set of 185 individuals collected in 76 A. thaliana Iberian wild populations were inoculated with different CMV strains. Resistance was estimated from the level of virus multiplication in infected plants, and tolerance from the effect of infection on host progeny production. Resistance and tolerance to CMV showed substantial genetic variation within and between host populations, and depended on the virus x host genotype interaction, two conditions for coevolution. Resistance and tolerance were co-occurring independent traits that have evolved independently from related life-history traits involved in adaptation to climate. The comparison of the genetic structure for resistance and tolerance with that for neutral traits (QST/FST analyses) indicated that both defence traits are likely under uniform selection. These results strongly suggest that CMV infection selects for defence on A. thaliana populations, and support plant-virus coevolution. Thus, we propose that CMV infection reduces host fitness under the field conditions of the wild A. thaliana populations studied.

Tracing the origin and evolution of pseudokinases across the tree of life

Protein phosphorylation by eukaryotic protein kinases (ePKs) is a fundamental mechanism of cell signaling in all organisms. In model vertebrates, ~10% of ePKs are classified as pseudokinases, which have amino acid changes within the catalytic machinery of the kinase domain that distinguish them from their canonical kinase counterparts. However, pseudokinases still regulate various signaling pathways, usually doing so in the absence of their own catalytic output. To investigate the prevalence, evolutionary relationships, and biological diversity of these pseudoenzymes, we performed a comprehensive analysis of putative pseudokinase sequences in available eukaryotic, bacterial, and archaeal proteomes. We found that pseudokinases are present across all domains of life, and we classified nearly 30,000 eukaryotic, 1500 bacterial, and 20 archaeal pseudokinase sequences into 86 pseudokinase families, including ~30 families that were previously unknown. We uncovered a rich variety of pseudokinases with notable expansions not only in animals but also in plants, fungi, and bacteria, where pseudokinases have previously received cursory attention. These expansions are accompanied by domain shuffling, which suggests roles for pseudokinases in plant innate immunity, plant-fungal interactions, and bacterial signaling. Mechanistically, the ancestral kinase fold has diverged in many distinct ways through the enrichment of unique sequence motifs to generate new families of pseudokinases in which the kinase domain is repurposed for noncanonical nucleotide binding or to stabilize unique, inactive kinase conformations. We further provide a collection of annotated pseudokinase sequences in the Protein Kinase Ontology (ProKinO) as a new mineable resource for the signaling community.

Expression polymorphism at the ARPC4 locus links the actin cytoskeleton with quantitative disease resistance to Sclerotinia sclerotiorum in Arabidopsis thaliana

Quantitative disease resistance (QDR) is a form of plant immunity widespread in nature, and the only one active against broad host range fungal pathogens. The genetic determinants of QDR are complex and largely unknown, and are thought to rely partly on genes controlling plant morphology and development. We used genome-wide association mapping in Arabidopsis thaliana to identify ARPC4 as associated with QDR against the necrotrophic fungal pathogen Sclerotinia sclerotiorum. Mutants impaired in ARPC4 showed enhanced susceptibility to S. sclerotiorum, defects in the structure of the actin filaments and in their responsiveness to S. sclerotiorum. Disruption of ARPC4 also alters callose deposition and the expression of defense-related genes upon S. sclerotiorum infection. Analysis of ARPC4 diversity in A. thaliana identified one haplotype (ARPC4^R ) showing a c. 1 kbp insertion in ARPC4 regulatory region and associated with higher level of QDR. Accessions from the ARPC4^R haplotype showed enhanced ARPC4 expression upon S. sclerotiorum challenge, indicating that polymorphisms in ARPC4 regulatory region are associated with enhanced QDR. This work identifies a novel actor of plant QDR against a fungal pathogen and provides a prime example of genetic mechanisms leading to the recruitment of cell morphology processes in plant immunity.

Contribution of recent technological advances to future resistance breeding

The development of durable host resistance strategies to control crop diseases is a primary need for sustainable agricultural production in the future. This article highlights the potential of recent progress in the understanding of host resistance for future cereal breeding. Much of the novel work is based on advancements in large-scale sequencing and genomics, rapid gene isolation techniques and high-throughput molecular marker technologies. Moreover, emerging applications on the pathogen side like effector identification or field pathogenomics are discussed. The combination of knowledge from both sides of cereal pathosystems will result in new approaches for resistance breeding. We describe future applications and innovative strategies to implement effective and durable strategies to combat diseases of major cereal crops while reducing pesticide dependency.

Genome-wide association study reveals new loci involved in Arabidopsis thaliana and Turnip mosaic virus (TuMV) interactions in the field

The genetic architecture of plant response to viruses has often been studied in model nonnatural pathosystems under controlled conditions. There is an urgent need to elucidate the genetic architecture of the response to viruses in a natural setting. A field experiment was performed in each of two years. In total, 317 Arabidopsis thaliana accessions were inoculated with its natural Turnip mosaic virus (TuMV). The accessions were phenotyped for viral accumulation, frequency of infected plants, stem length and symptoms. Genome-wide association mapping was performed. Arabidopsis thaliana exhibits extensive natural variation in its response to TuMV in the field. The underlying genetic architecture reveals a more quantitative picture than in controlled conditions. Ten genomic regions were consistently identified across the two years. RTM3 (Restricted TEV Movement 3) is a major candidate for the response to TuMV in the field. New candidate genes include Dead box helicase 1, a Tim Barrel domain protein and the eukaryotic translation initiation factor eIF3b. To our knowledge, this study is the first to report the genetic architecture of quantitative response of A. thaliana to a naturally occurring virus in a field environment, thereby highlighting relevant candidate genes involved in plant virus interactions in nature.

Long-term balancing selection drives evolution of immunity genes in Capsella

Genetic drift is expected to remove polymorphism from populations over long periods of time, with the rate of polymorphism loss being accelerated when species experience strong reductions in population size. Adaptive forces that maintain genetic variation in populations, or balancing selection, might counteract this process. To understand the extent to which natural selection can drive the retention of genetic diversity, we document genomic variability after two parallel species-wide bottlenecks in the genus Capsella. We find that ancestral variation preferentially persists at immunity related loci, and that the same collection of alleles has been maintained in different lineages that have been separated for several million years. By reconstructing the evolution of the disease-related locus MLO2b, we find that divergence between ancient haplotypes can be obscured by referenced based re-sequencing methods, and that trans-specific alleles can encode substantially diverged protein sequences. Our data point to long-term balancing selection as an important factor shaping the genetics of immune systems in plants and as the predominant driver of genomic variability after a population bottleneck.

Novel organ-specific genetic factors for quantitative resistance to late blight in potato

Potato, Solanum tuberosum, is one of the major consumed food in the world, being the basis of the diet of millions of people. The main limiting and destructive disease of potato is late blight, caused by Phytophtora infestans. Here, we present a multi-environmental analysis of the response to P. infestans using an association panel of 150 accessions of S. tuberosum Group Phureja, evaluated in two localities in Colombia. Disease resistance data were merged with a genotyping matrix of 83,862 SNPs obtained by 2b-restriction site-associated DNA and Genotyping by sequencing approaches into a Genome-wide association study. We are reporting 16 organ-specific QTL conferring resistance to late blight. These QTL explain between 13.7% and 50.9% of the phenotypic variance. Six and ten QTL were detected for resistance response in leaves and stem, respectively. In silico analysis revealed 15 candidate genes for resistance to late blight. Four of them have no functional genome annotation, while eleven candidate genes code for diverse proteins, including a leucine-rich repeat kinase.

Infection Assay for Xanthomonas campestris pv. campestris in Arabidopsis thaliana Mimicking Natural Entry via Hydathodes

Xanthomonas campestris pv. campestris (Xcc) causes the devastating disease Black rot in Brassicaceae. Typically Xcc enters the plant through specialized organs on the leaf margin, called hydathodes, and spreads from there through the vasculature. In order to mimic natural entry as closely as possible, we here describe a "hydathode guttation"-based entry assay for Xcc in Arabidopsis. This disease assay combines spray inoculation with the induction of guttation and allows reabsorption of guttation droplets by the plant. Moreover, our assay relies on a bioluminescent reporter strain of Xcc to allow direct visualization of both entry and subsequent spreading of Xcc in its host. The assay allows the routine infection from one to two hydathodes per Arabidopsis leaf. Infections are scored 14 days post inoculation, just before the infection goes systemic.

Recent Advances in Mechanisms of Plant Defense to Sclerotinia sclerotiorum

Sclerotinia sclerotiorum (Lib.) de Bary is an unusual pathogen which has the broad host range, diverse infection modes, and potential double feeding lifestyles of both biotroph and necrotroph. It is capable of infecting over 400 plant species found worldwide and more than 60 names have agriculturally been used to refer to diseases caused by this pathogen. Plant defense to S. sclerotiorum is a complex biological process and exhibits a typical quantitative disease resistance (QDR) response. Recent studies using Arabidopsis thaliana and crop plants have obtained new advances in mechanisms used by plants to cope with S. sclerotiorum infection. In this review, we focused on our current understanding on plant defense mechanisms against this pathogen, and set up a model for the defense process including three stages: recognition of this pathogen, signal transduction and defense response. We also have a particular interest in defense signaling mediated by diverse signaling molecules. We highlight the current challenges and unanswered questions in both the defense process and defense signaling. Essentially, we discussed candidate resistance genes newly mapped by using high-throughput experiments in important crops, and classified these potential gene targets into different stages of the defense process, which will broaden our understanding of the genetic architecture underlying quantitative resistance to S. sclerotiorum. We proposed that more powerful mapping population(s) will be required for accurate and reliable QDR gene identification.

Genome-wide association studies on the phyllosphere microbiome: Embracing complexity in host-microbe interactions

Environmental sequencing shows that plants harbor complex communities of microbes that vary across environments. However, many approaches for mapping plant genetic variation to microbe-related traits were developed in the relatively simple context of binary host-microbe interactions under controlled conditions. Recent advances in sequencing and statistics make genome-wide association studies (GWAS) an increasingly promising approach for identifying the plant genetic variation associated with microbes in a community context. This review discusses early efforts on GWAS of the plant phyllosphere microbiome and the outlook for future studies based on human microbiome GWAS. A workflow for GWAS of the phyllosphere microbiome is then presented, with particular attention to how perspectives on the mechanisms, evolution and environmental dependence of plant-microbe interactions will influence the choice of traits to be mapped.

Assessment of heterosis in two Arabidopsis thaliana common-reference mapping populations

Hybrid vigour, or heterosis, has been of tremendous importance in agriculture for the improvement of both crops and livestock. Notwithstanding large efforts to study the phenomenon of heterosis in the last decades, the identification of common molecular mechanisms underlying hybrid vigour remain rare. Here, we conducted a systematic survey of the degree of heterosis in Arabidopsis thaliana hybrids. For this purpose, two overlapping Arabidopsis hybrid populations were generated by crossing a large collection of naturally occurring accessions to two common reference lines. In these Arabidopsis hybrid populations the range of heterosis for several developmental and yield related traits was examined, and the relationship between them was studied. The traits under study were projected leaf area at 17 days after sowing, flowering time, height of the main inflorescence, number of side branches from the main stem or from the rosette base, total seed yield, seed weight, seed size and the estimated number of seeds per plant. Predominantly positive heterosis was observed for leaf area and height of the main inflorescence, whereas mainly negative heterosis was observed for rosette branching. For the other traits both positive and negative heterosis was observed in roughly equal amounts. For flowering time and seed size only low levels of heterosis were detected. In general the observed heterosis levels were highly trait specific. Furthermore, no correlation was observed between heterosis levels and the genetic distance between the parental lines. Since all selected lines were a part of the Arabidopsis genome wide association (GWA) mapping panel, a genetic mapping approach was applied to identify possible regions harbouring genetic factors causal for heterosis, with separate calculations for additive and dominance effects. Our study showed that the genetic mechanisms underlying heterosis were highly trait specific in our hybrid populations and greatly depended on the genetic background, confirming the elusive character of heterosis.

Natural variation among Arabidopsis thaliana accessions in tolerance to high magnesium supply

High magnesium (Mg2+) in some extreme serpentine soils or semi-arid regions is an important factor affecting crop growth and development. Specific loci that form the genetic framework underlying high Mg2+ homeostasis, however, are not well understood. By using GWA mapping on 388 accessions of Arabidopsis thaliana selected from a worldwide collection and genotyped at approximately 250,00 SNPs, we successfully identified 109 and 74 putative genetic regions associated in nutrient traits under normal (1,000 µM) and high Mg2+ (10,000 µM), respectively. Above 90% SNPs associated with nutrients including Mg2+ and only two SNPs shared between normal and high Mg2+. A single strong peak of SNPs associated with Ca concentration corresponding to candidate gene At1g60420 ARABIDOPSIS NUCLEOREDOXIN (AtNRX1) under high Mg2+ was further determined. Compared with wildtype, mutants of Atnrx1-1 and Atnrx1-2 supplied with high Mg2+ had higher Ca concentrations in the plant, and higher cytosolic Ca2+ concentrations during root elongation, as well as higher fresh weight and lateral-root number. This suggests that AtNRX1 was a critical gene negatively regulating Ca uptake under high Mg2+ conditions. The discovery could help to breed/select crops that can adapt to high-Mg2+ soils such as serpentine soils (high ratio of Mg2+: Ca2+) or Mars soil with high levels of magnesium sulfate.

Genome-wide association mapping in Arabidopsis identifies novel genes underlying quantitative disease resistance to Alternaria brassicae

Quantitative disease resistance (QDR) is the predominant form of resistance against necrotrophic pathogens. The genes and mechanisms underlying QDR are not well known. In the current study, the Arabidopsis-Alternaria brassicae pathosystem was used to uncover the genetic architecture underlying resistance to A. brassicae in a set of geographically diverse Arabidopsis accessions. Arabidopsis accessions revealed a rich variation in the host responses to the pathogen, varying from complete resistance to high susceptibility. Genome-wide association (GWA) mapping revealed multiple regions to be associated with disease resistance. A subset of genes prioritized on the basis of gene annotations and evidence of transcriptional regulation in other biotic stresses was analysed using a reverse genetics approach employing T-DNA insertion mutants. The mutants of three genes, namely At1g06990 (GDSL-motif lipase), At3g25180 (CYP82G1) and At5g37500 (GORK), displayed an enhanced susceptibility relative to the wild-type. These genes are involved in the development of morphological phenotypes (stomatal aperture) and secondary metabolite synthesis, thus defining some of the diverse facets of quantitative resistance against A. brassicae.

Two-way mixed-effects methods for joint association analysis using both host and pathogen genomes

Infectious diseases are often affected by specific pairings of hosts and pathogens and therefore by both of their genomes. The integration of a pair of genomes into genome-wide association mapping can provide an exquisitely detailed view of the genetic landscape of complex traits. We present a statistical method, ATOMM (Analysis with a Two-Organism Mixed Model), that maps a trait of interest to a pair of genomes simultaneously; this method makes use of whole-genome sequence data for both host and pathogen organisms. ATOMM uses a two-way mixed-effect model to test for genetic associations and cross-species genetic interactions while accounting for sample structure including interactions between the genetic backgrounds of the two organisms. We demonstrate the applicability of ATOMM to a joint association study of quantitative disease resistance (QDR) in the Arabidopsis thaliana-Xanthomonas arboricola pathosystem. Our method uncovers a clear host-strain specificity in QDR and provides a powerful approach to identify genetic variants on both genomes that contribute to phenotypic variation.

In situ relationships between microbiota and potential pathobiota in Arabidopsis thaliana

A current challenge in microbial pathogenesis is to identify biological control agents that may prevent and/or limit host invasion by microbial pathogens. In natura, hosts are often infected by multiple pathogens. However, most of the current studies have been performed under laboratory controlled conditions and by taking into account the interaction between a single commensal species and a single pathogenic species. The next step is therefore to explore the relationships between host-microbial communities (microbiota) and microbial members with potential pathogenic behavior (pathobiota) in a realistic ecological context. In the present study, we investigated such relationships within root-associated and leaf-associated bacterial communities of 163 ecologically contrasted Arabidopsis thaliana populations sampled across two seasons in southwest of France. In agreement with the theory of the invasion paradox, we observed a significant humped-back relationship between microbiota and pathobiota α-diversity that was robust between both seasons and plant organs. In most populations, we also observed a strong dynamics of microbiota composition between seasons. Accordingly, the potential pathobiota composition was explained by combinations of season-specific microbiota operational taxonomic units. This result suggests that the potential biomarkers controlling pathogen's invasion are highly dynamic.