Specific resistances against Pseudomonas syringae effectors AvrB and AvrRpm1 have evolved differently in common bean (Phaseolus vulgaris), soybean (Glycine max), and Arabidopsis thaliana

Fine-mapping and evolutionary history of R-BPMV, a dominant resistance gene to Bean pod mottle virus in Phaseolus vulgaris L

KEY MESSAGE

R-BPMV is located within a recently expanded TNL cluster in the Phaseolus genus with suppressed recombination and known for resistance to multiple pathogens including potyviruses controlled by the I gene. Bean pod mottle virus (BPMV) is a comovirus that infects common bean and legumes in general. BPMV is distributed throughout the world and is a major threat on soybean, a closely related species of common bean. In common bean, BAT93 was reported to carry the R-BPMV resistance gene conferring resistance to BPMV and linked with the I resistance gene. To fine map R-BPMV, 182 recombinant inbred lines (RILs) derived from the cross BAT93 × JaloEEP558 were genotyped with polymerase chain reaction (PCR)-based markers developed using genome assemblies from G19833 and BAT93, as well as BAT93 BAC clone sequences. Analysis of RILs carrying key recombination events positioned R-BPMV to a target region containing at least 16 TIR-NB-LRR (TNL) sequences in BAT93. Because the I cluster presents a suppression of recombination and a large number of repeated sequences, none of the 16 TNLs could be excluded as R-BPMV candidate gene. The evolutionary history of the TNLs for the I cluster were reconstructed using microsynteny and phylogenetic analyses within the legume family. A single I TNL was present in Medicago truncatula and lost in soybean, mirroring the absence of complete BPMV resistance in soybean. Amplification of TNLs in the I cluster predates the divergence of the Phaseolus species, in agreement with the emergence of R-BPMV before the separation of the common bean wild centers of diversity. This analysis provides PCR-based markers useful in marker-assisted selection (MAS) and laid the foundation for cloning of R-BPMV resistance gene in order to transfer the resistance into soybean.

Diversified host target families mediate convergently evolved effector recognition across plant species

Recognition of pathogen effectors is a crucial step for triggering plant immunity. Resistance (R) genes often encode for nucleotide-binding leucine-rich repeat receptors (NLRs), and NLRs detect effectors from pathogens to trigger effector-triggered immunity (ETI). NLR recognition of effectors is observed in diverse forms where NLRs directly interact with effectors or indirectly detect effectors by monitoring host guardees/decoys (HGDs). HGDs undergo different biochemical modifications by diverse effectors and expand the effector recognition spectrum of NLRs, contributing robustness to plant immunity. Interestingly, in many cases of the indirect recognition of effectors, HGD families targeted by effectors are conserved across the plant species while NLRs are not. Notably, a family of diversified HGDs can activate multiple non-orthologous NLRs across plant species. Further investigation on HGDs would reveal the mechanistic basis of how the diversification of HGDs confers novel effector recognition by NLRs.

A Puccinia striiformis f. sp. tritici effector inhibits high-temperature seedling-plant resistance in wheat

Resistance to Pseudomonas syringae pv. maculicola 1 (RPM1)-induced protein kinase (RIPK) in Arabidopsis belongs to the receptor-like cytoplasmic kinase (RLCK) family and plays a vital role in immunity. However, the role of RLCKs in the high-temperature seedling-plant (HTSP) resistance of wheat (Triticum aestivum) to Puccinia striiformis f. sp. tritici (Pst), the stripe rust pathogen, remains unclear. Here, we identified a homologous gene of RIPK in wheat, namely TaRIPK. Expression of TaRIPK was induced by Pst inoculation and high temperatures. Silencing of TaRIPK reduced the expression level of TaRPM1, resulting in weaker HTSP resistance. Moreover, TaRIPK interacts with and phosphorylates papain-like cysteine protease 1 (TaPLCP1). Meanwhile, we found that the Pst-secreted protein PSTG_01766 targets TaPLCP1. Transient expression of PSTG_01766 inhibited basal immunity in tobacco (Nicotiana benthamiana) and wheat. The role of PSTG_01766 as an effector involved in HTSP resistance was further supported by host-induced gene silencing and bacterial type three secretion system-mediated delivery into wheat. PSTG_01766 inhibited the TaRIPK-induced phosphorylation of TaPLCP1. Furthermore, PSTG_01766 has the potential to influence the subcellular localization of TaPLCP1. Overall, we suggest that the TaRIPK-TaPLCP1-TaRPM1 module fits the guard model for disease resistance, participating in HTSP resistance. PSTG_01766 decreases HTSP resistance via targeting TaPLCP1. Guarded by wheat and attacked by Pst, TaPLCP1 may serve as a central hub of the defense response. Our findings improve the understanding of the molecular mechanism of wheat HTSP resistance, which may be an important strategy for controlling stripe rust in the face of global warming.

Genetics of autoimmunity in plants: an evolutionary genetics perspective

Autoimmunity in plants has been found in numerous hybrids as a form of hybrid necrosis and mutant panels. Uncontrolled cell death is a main cellular outcome of autoimmunity, which negatively impacts growth. Its occurrence highlights the vulnerable nature of the plant immune system. Genetic investigation of autoimmunity in hybrid plants revealed that extreme variation in the immune receptor repertoire is a major contributor, reflecting an evolutionary conundrum that plants face in nature. In this review, we discuss natural variation in the plant immune system and its contribution to fitness. The value of autoimmunity genetics lies in its ability to identify combinations of a natural immune receptor and its partner that are predisposed to triggering autoimmunity. The network of immune components for autoimmunity becomes instrumental in revealing mechanistic details of how immune receptors recognize cellular invasion and activate signaling. The list of autoimmunity-risk variants also allows us to infer evolutionary processes contributing to their maintenance in the natural population. Our approach to autoimmunity, which integrates mechanistic understanding and evolutionary genetics, has the potential to serve as a prognosis tool to optimize immunity in crops.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

BACKGROUND

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

RESULTS

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

CONCLUSIONS

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

Common bean resistance to Xanthomonas is associated with upregulation of the salicylic acid pathway and downregulation of photosynthesis

Background

Common bacterial blight (CBB) caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans is one of the major threats to common bean crops (Phaseolus vulgaris L.). Resistance to CBB is particularly complex as 26 quantitative resistance loci to CBB have been described so far. To date, transcriptomic studies after CBB infection have been very scarce and the molecular mechanisms underlying susceptibility or resistance are largely unknown.

Results

We sequenced and annotated the genomes of two common bean genotypes being either resistant (BAT93) or susceptible (JaloEEP558) to CBB. Reciprocal BLASTp analysis led to a list of 20,787 homologs between these genotypes and the common bean reference genome (G19833), which provides a solid dataset for further comparative analyses. RNA-Seq after inoculation with X. phaseoli pv. phaseoli showed that the susceptible genotype initiated a more intense and diverse biological response than the resistant genotype. Resistance was linked to upregulation of the salicylic acid pathway and downregulation of photosynthesis and sugar metabolism, while susceptibility was linked to downregulation of resistance genes and upregulation of the ethylene pathway and of genes involved in cell wall modification.

Conclusions

This study helps better understanding the mechanisms occurring during the early colonization phase of common bean by Xanthomonas and unveils new actors potentially important for resistance and susceptibility to CBB. We discuss the potential link between the pathways induced during bean colonization and genes induced by transcription activator-like effectors (TALEs), as illustrated in other Xanthomonas pathovars.

Genome-Wide Mining of Disease Resistance Gene Analogs Using Conserved Domains

The production of legume crop species is severely affected by disease, imposing a significant yield loss annually worldwide. Plant resistance gene analogs (RGAs) play specific roles in plant resistance responses, and their identification and subsequent application in breeding programs help to reduce this yield loss. RGAs contain conserved domains and motifs, which can be used for their identification and classification. Nucleotide-binding site-leucine-rich repeat (NLR), receptor like kinase (RLK), and receptor like protein (RLP) genes are the main types of RGAs. Computational identification and characterization of RGAs has been performed successfully among different plant species. Here, we explain the computational workflow for genome-wide RGA identification in legumes.

PRRs and NB-LRRs: From Signal Perception to Activation of Plant Innate Immunity

To ward off pathogens and pests, plants use a sophisticated immune system. They use pattern-recognition receptors (PRRs), as well as nucleotide-binding and leucine-rich repeat (NB-LRR) domains, for detecting nonindigenous molecular signatures from pathogens. Plant PRRs induce local and systemic immunity. Plasma-membrane-localized PRRs are the main components of multiprotein complexes having additional transmembrane and cytosolic kinases. Topical research involving proteins and their interactive partners, along with transcriptional and posttranscriptional regulation, has extended our understanding of R-gene-mediated plant immunity. The unique LRR domain conformation helps in the best utilization of a surface area and essentially mediates protein–protein interactions. Genome-wide analyses of inter- and intraspecies PRRs and NB-LRRs offer innovative information about their working and evolution. We reviewed plant immune responses with relevance to PRRs and NB-LRRs. This article focuses on the significant functional diversity, pathogen-recognition mechanisms, and subcellular compartmentalization of plant PRRs and NB-LRRs. We highlight the potential biotechnological application of PRRs and NB-LRRs to enhance broad-spectrum disease resistance in crops.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Common Bean Subtelomeres Are Hot Spots of Recombination and Favor Resistance Gene Evolution

Subtelomeres of most eukaryotes contain fast-evolving genes usually involved in adaptive processes. In common bean (Phaseolus vulgaris), the Co-2 anthracnose resistance (R) locus corresponds to a cluster of nucleotide-binding-site leucine-rich-repeat (NL) encoding sequences, the prevalent class of plant R genes. To study the recent evolution of this R gene cluster, we used a combination of sequence, genetic and cytogenetic comparative analyses between common bean genotypes from two distinct gene pools (Andean and Mesoamerican) that diverged 0.165 million years ago. Co-2 is a large subtelomeric cluster on chromosome 11 comprising from 32 (Mesoamerican) to 52 (Andean) NL sequences embedded within khipu satellite repeats. Since the recent split between Andean and Mesoamerican gene pools, the Co-2 cluster has experienced numerous gene-pool specific NL losses, leading to distinct NL repertoires. The high proportion of solo-LTR retrotransposons indicates that the Co-2 cluster is located in a hot spot of unequal intra-strand homologous recombination. Furthermore, we observe large segmental duplications involving both Non-Homologous End Joining and Homologous Recombination double-strand break repair pathways. Finally, the identification of a Mesoamerican-specific subtelomeric sequence reveals frequent interchromosomal recombinations between common bean subtelomeres. Altogether, our results highlight that common bean subtelomeres are hot spots of recombination and favor the rapid evolution of R genes. We propose that chromosome ends could act as R gene incubators in many plant genomes.

Dissection of Resistance Genes to Pseudomonas syringae pv. phaseolicola in UI3 Common Bean Cultivar

Few quantitative trait loci have been mapped for resistance to Pseudomonas syringae pv. phaseolicola in common bean. Two F₂ populations were developed from the host differential UI3 cultivar. The objective of this study was to further characterize the resistance to races 1, 5, 7 and 9 of Psp included in UI3. Using a QTL mapping approach, 16 and 11 main-effect QTLs for pod and primary leaf resistance were located on LG10, explaining up to 90% and 26% of the phenotypic variation, respectively. The homologous genomic region corresponding to primary leaf resistance QTLs detected tested positive for the presence of resistance-associated gene cluster encoding nucleotide-binding and leucine-rich repeat (NL), Natural Resistance Associated Macrophage (NRAMP) and Pentatricopeptide Repeat family (PPR) proteins. It is worth noting that the main effect QTLs for resistance in pod were located inside a 3.5 Mb genomic region that included the Phvul.010G021200 gene, which encodes a protein that has the highest sequence similarity to the RIN4 gene of Arabidopsis, and can be considered an important candidate gene for the organ-specific QTLs identified here. These results support that resistance to Psp from UI3 might result from the immune response activated by combinations of R proteins, and suggest the guard model as an important mechanism in pod resistance to halo blight. The candidate genes identified here warrant functional studies that will help in characterizing the actual defense gene(s) in UI3 genotype.

Comparative analysis of NBS-LRR genes and their response to Aspergillus flavus in Arachis

Studies have demonstrated that nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes respond to pathogen attack in plants. Characterization of NBS-LRR genes in peanut is not well documented. The newly released whole genome sequences of Arachis duranensis and Arachis ipaënsis have allowed a global analysis of this important gene family in peanut to be conducted. In this study, we identified 393 (AdNBS) and 437 (AiNBS) NBS-LRR genes from A. duranensis and A. ipaënsis, respectively, using bioinformatics approaches. Full-length sequences of 278 AdNBS and 303 AiNBS were identified. Fifty-one orthologous, four AdNBS paralogous, and six AiNBS paralogous gene pairs were predicted. All paralogous gene pairs were located in the same chromosomes, indicating that tandem duplication was the most likely mechanism forming these paralogs. The paralogs mainly underwent purifying selection, but most LRR 8 domains underwent positive selection. More gene clusters were found in A. ipaënsis than in A. duranensis, possibly owing to tandem duplication events occurring more frequently in A. ipaënsis. The expression profile of NBS-LRR genes was different between A. duranensis and A. hypogaea after Aspergillus flavus infection. The up-regulated expression of NBS-LRR in A. duranensis was continuous, while these genes responded to the pathogen temporally in A. hypogaea.

Molecular Soybean-Pathogen Interactions

Soybean hosts a wide variety of pathogens that cause significant yield losses. The importance of soybean as a major oilseed crop has led to research focused on its interactions with pathogens, such as Soybean mosaic virus, Pseudomonas syringae, Phytophthora sojae, Phakopsora pachyrhizi, and Heterodera glycines. Pioneering work on soybean's interactions with these organisms, which represent the five major pathogen groups (viruses, bacteria, oomycetes, fungi, and nematodes), has contributed to our understanding of the molecular mechanisms underlying virulence and immunity. These mechanisms involve conserved and unique features that validate the need for research in both soybean and homologous model systems. In this review, we discuss identification of effectors and their functions as well as resistance gene-mediated recognition and signaling. We also point out areas in which model systems and recent advances in resources and tools have provided opportunities to gain deeper insights into soybean-pathogen interactions.

Development of molecular markers linked to disease resistance genes in common bean based on whole genome sequence

Common bean (Phaseolus vulgaris) is the most important grain legume for direct human consumption in the world, particularly in developing countries where it constitutes the main source of protein. Unfortunately, common bean yield stability is constrained by a number of pests and diseases. As use of resistant genotypes is the most economic and ecologically safe means for controlling plant diseases, efforts have been made to genetically characterize resistance genes (R genes) in common bean. Despite its agronomic importance, genomic resources available in common bean were limited until the recent sequencing of common bean genome (Andean genotype G19833). Besides allowing the annotation of Nucleotide Binding-Leucine Rich Repeat (NB-LRR) encoding gene family, which is the prevalent class of disease R genes in plants, access to the whole genome sequence of common bean can be of great help for intense selection to increase the overall efficiency of crop improvement programs using marker-assisted selection (MAS). This review presents the state of the art of common bean NB-LRR gene clusters, their peculiar location in subtelomeres and correlation with genetically characterized monogenic R genes, as well as how the availability of the whole genome sequence can boost the development of molecular markers for MAS.

Dissection of Two Complex Clusters of Resistance Genes in Lettuce (Lactuca sativa)

Of the over 50 phenotypic resistance genes mapped in lettuce, 25 colocalize to three major resistance clusters (MRC) on chromosomes 1, 2, and 4. Similarly, the majority of candidate resistance genes encoding nucleotide binding-leucine rich repeat (NLR) proteins genetically colocalize with phenotypic resistance loci. MRC1 and MRC4 span over 66 and 63 Mb containing 84 and 21 NLR-encoding genes, respectively, as well as 765 and 627 genes that are not related to NLR genes. Forward and reverse genetic approaches were applied to dissect MRC1 and MRC4. Transgenic lines exhibiting silencing were selected using silencing of β-glucuronidase as a reporter. Silencing of two of five NLR-encoding gene families resulted in abrogation of nine of 14 tested resistance phenotypes mapping to these two regions. At MRC1, members of the coiled coil-NLR-encoding RGC1 gene family were implicated in host and nonhost resistance through requirement for Dm5/8- and Dm45-mediated resistance to downy mildew caused by Bremia lactucae as well as the hypersensitive response to effectors AvrB, AvrRpm1, and AvrRpt2 of the nonpathogen Pseudomonas syringae. At MRC4, RGC12 family members, which encode toll interleukin receptor-NLR proteins, were implicated in Dm4-, Dm7-, Dm11-, and Dm44-mediated resistance to B. lactucae. Lesions were identified in the sequence of a candidate gene within dm7 loss-of-resistance mutant lines, confirming that RGC12G confers Dm7.

Understanding plant immunity as a surveillance system to detect invasion

Various conceptual models to describe the plant immune system have been presented. The most recent paradigm to gain wide acceptance in the field is often referred to as the zigzag model, which reconciles the previously formulated gene-for-gene hypothesis with the recognition of general elicitors in a single model. This review focuses on the limitations of the current paradigm of molecular plant-microbe interactions and how it too narrowly defines the plant immune system. As such, we discuss an alternative view of plant innate immunity as a system that evolves to detect invasion. This view accommodates the range from mutualistic to parasitic symbioses that plants form with diverse organisms, as well as the spectrum of ligands that the plant immune system perceives. Finally, how this view can contribute to the current practice of resistance breeding is discussed.

Determining the GmRIN4 requirements of the soybean disease resistance proteins Rpg1b and Rpg1r using a nicotiana glutinosa-based agroinfiltration system

Rpg1b and Rpg1r are soybean disease resistance (R) genes responsible for conferring resistance to Pseudomonas syringae strains expressing the effectors AvrB and AvrRpm1, respectively. The study of these cloned genes would be greatly facilitated by the availability of a suitable transient expression system. The commonly used Niciotiana benthamiana-based system is not suitable for studying Rpg1b and Rpg1r function, however, because expression of AvrB or AvrRpm1 alone induces a hypersensitive response (HR), indicating that N. benthamiana contains endogenous R genes that recognize these effectors. To identify a suitable alternative host for transient expression assays, we screened 13 species of Nicotiana along with 11 accessions of N. tabacum for lack of response to transient expression of AvrB and AvrRpm1. We found that N. glutinosa did not respond to either effector and was readily transformable as determined by transient expression of β-glucuronidase. Using this system, we determined that Rpg1b-mediated HR in N. glutinosa required co-expression of avrB and a soybean ortholog of the Arabidopsis RIN4 gene. All four soybean RIN4 orthologs tested worked in the assay. In contrast, Rpg1r did not require co-expression of a soybean RIN4 ortholog to recognize AvrRpm1, but recognition was suppressed by co-expression with AvrRpt2. These observations suggest that an endogenous RIN4 gene in N. glutinosa can substitute for the soybean RIN4 ortholog in the recognition of AvrRpm1 by Rpg1r.

Evolutionary relationship of disease resistance genes in soybean and Arabidopsis specific for the Pseudomonas syringae effectors AvrB and AvrRpm1

In Arabidopsis (Arabidopsis thaliana), the Pseudomonas syringae effector proteins AvrB and AvrRpm1 are both detected by the RESISTANCE TO PSEUDOMONAS MACULICOLA1 (RPM1) disease resistance (R) protein. By contrast, soybean (Glycine max) can distinguish between these effectors, with AvrB and AvrRpm1 being detected by the Resistance to Pseudomonas glycinea 1b (Rpg1b) and Rpg1r R proteins, respectively. We have been using these genes to investigate the evolution of R gene specificity and have previously identified RPM1 and Rpg1b. Here, we report the cloning of Rpg1r, which, like RPM1 and Rpg1b, encodes a coiled-coil (CC)-nucleotide-binding (NB)-leucine-rich repeat (LRR) protein. As previously found for Rpg1b, we determined that Rpg1r is not orthologous with RPM1, indicating that the ability to detect both AvrB and AvrRpm1 evolved independently in soybean and Arabidopsis. The tightly linked soybean Rpg1b and Rpg1r genes share a close evolutionary relationship, with Rpg1b containing a recombination event that combined a NB domain closely related to Rpg1r with CC and LRR domains from a more distantly related CC-NB-LRR gene. Using structural modeling, we mapped polymorphisms between Rpg1b and Rpg1r onto the predicted tertiary structure of Rpg1b, which revealed highly polymorphic surfaces within both the CC and LRR domains. Assessment of chimeras between Rpg1b and Rpg1r using a transient expression system revealed that AvrB versus AvrRpm1 specificity is determined by the C-terminal portion of the LRR domain. The P. syringae effector AvrRpt2, which targets RPM1 INTERACTOR4 (RIN4) proteins in both Arabidopsis and soybean, partially blocked recognition of both AvrB and AvrRpm1 in soybean, suggesting that both Rpg1b and Rpg1r may detect these effectors via modification of a RIN4 homolog.

Fine mapping of Co-x, an anthracnose resistance gene to a highly virulent strain of Colletotrichum lindemuthianum in common bean

The Co - x anthracnose R gene of common bean was fine-mapped into a 58 kb region at one end of chromosome 1, where no canonical NB-LRR-encoding genes are present in G19833 genome sequence. Anthracnose, caused by the phytopathogenic fungus Colletotrichum lindemuthianum, is one of the most damaging diseases of common bean, Phaseolus vulgaris. Various resistance (R) genes, named Co-, conferring race-specific resistance to different strains of C. lindemuthianum have been identified. The Andean cultivar JaloEEP558 was reported to carry Co-x on chromosome 1, conferring resistance to the highly virulent strain 100. To fine map Co-x, 181 recombinant inbred lines derived from the cross between JaloEEP558 and BAT93 were genotyped with polymerase chain reaction (PCR)-based markers developed using the genome sequence of the Andean genotype G19833. Analysis of RILs carrying key recombination events positioned Co-x at one end of chromosome 1 to a 58 kb region of the G19833 genome sequence. Annotation of this target region revealed eight genes: three phosphoinositide-specific phospholipases C (PI-PLC), one zinc finger protein and four kinases, suggesting that Co-x is not a classical nucleotide-binding leucine-rich encoding gene. In addition, we identified and characterized the seven members of common bean PI-PLC gene family distributed into two clusters located at the ends of chromosomes 1 and 8. Co-x is not a member of Co-1 allelic series since these two genes are separated by at least 190 kb. Comparative analysis between soybean and common bean revealed that the Co-x syntenic region, located at one end of Glycine max chromosome 18, carries Rhg1, a major QTL contributing to soybean cyst nematode resistance. The PCR-based markers generated in this study should be useful in marker-assisted selection for pyramiding Co-x with other R genes.

A reference genome for common bean and genome-wide analysis of dual domestications

Common bean (Phaseolus vulgaris L.) is the most important grain legume for human consumption and has a role in sustainable agriculture owing to its ability to fix atmospheric nitrogen. We assembled 473 Mb of the 587-Mb genome and genetically anchored 98% of this sequence in 11 chromosome-scale pseudomolecules. We compared the genome for the common bean against the soybean genome to find changes in soybean resulting from polyploidy. Using resequencing of 60 wild individuals and 100 landraces from the genetically differentiated Mesoamerican and Andean gene pools, we confirmed 2 independent domestications from genetic pools that diverged before human colonization. Less than 10% of the 74 Mb of sequence putatively involved in domestication was shared by the two domestication events. We identified a set of genes linked with increased leaf and seed size and combined these results with quantitative trait locus data from Mesoamerican cultivars. Genes affected by domestication may be useful for genomics-enabled crop improvement.

Genetic analysis of the response to eleven Colletotrichum lindemuthianum races in a RIL population of common bean (Phaseolus vulgaris L.)

BACKGROUND

Bean anthracnose is caused by the fungus Colletotrichum lindemuthianum (Sacc. & Magnus) Lams.- Scrib. Resistance to C. lindemuthianum in common bean (Phaseolus vulgaris L.) generally follows a qualitative mode of inheritance. The pathogen shows extensive pathogenic variation and up to 20 anthracnose resistance loci (named Co-), conferring resistance to specific races, have been described. Anthracnose resistance has generally been investigated by analyzing a limited number of isolates or races in segregating populations. In this work, we analyzed the response against eleven C. lindemuthianum races in a recombinant inbred line (RIL) common bean population derived from the cross Xana × Cornell 49242 in which a saturated linkage map was previously developed.

RESULTS

A systematic genetic analysis was carried out to dissect the complex resistance segregations observed, which included contingency analyses, subpopulations and genetic mapping. Twenty two resistance genes were identified, some with a complementary mode of action. The Cornell 49242 genotype carries a complex cluster of resistance genes at the end of linkage group (LG) Pv11 corresponding to the previously described anthracnose resistance cluster Co-2. In this position, specific resistance genes to races 3, 6, 7, 19, 38, 39, 65, 357, 449 and 453 were identified, with one of them showing a complementary mode of action. In addition, Cornell 49242 had an independent gene on LG Pv09 showing a complementary mode of action for resistance to race 453. Resistance genes in genotype Xana were located on three regions involving LGs Pv01, Pv02 and Pv04. All resistance genes identified in Xana showed a complementary mode of action, except for two controlling resistance to races 65 and 73 located on LG Pv01, in the position of the previously described anthracnose resistance cluster Co-1.

CONCLUSIONS

Results shown herein reveal a complex and specific interaction between bean and fungus genotypes leading to anthracnose resistance. Organization of specific resistance genes in clusters including resistance genes with different modes of action (dominant and complementary genes) was also confirmed. Finally, new locations for anthracnose resistance genes were identified in LG Pv09.

The Subtelomeric khipu Satellite Repeat from Phaseolus vulgaris: Lessons Learned from the Genome Analysis of the Andean Genotype G19833

Subtelomeric regions in eukaryotic organisms are known for harboring species-specific tandemly repeated satellite sequences. However, studies on the molecular organization and evolution of subtelomeric repeats are scarce, especially in plants. Khipu is a satellite DNA of 528-bp repeat unit, specific of the Phaseolus genus, with a subtelomeric distribution in common bean, P. vulgaris. To investigate the genomic organization and the evolution of khipu, we performed genome-wide analysis on the complete genome sequence of the common bean genotype G19833. We identified 2,460 khipu units located at most distal ends of the sequenced regions. Khipu units are arranged in discrete blocks of 2-55 copies and are heterogeneously distributed among the different chromosome ends of G19833 (from 0 to 555 khipus units per chromosome arm). Phylogenetically related khipu units are spread between numerous chromosome ends, suggesting frequent exchanges between non-homologous subtelomeres. However, most subclades contain numerous khipu units from only one or few chromosome ends indicating that local duplication is also driving khipu expansion. Unexpectedly, we also identified 81 khipu units located at centromeres. All the centromeric khipu units belong to a single divergent clade also comprised of a few units from several subtelomeres, suggesting that a few sequence exchanges between centromeres and subtelomeres took place in the common bean genome. The divergence and low copy number of these centromeric units from the subtelomeric units could explain why they were not detected by FISH (Fluorescence in situ Hybridization) although it can not be excluded that these centromeric units may have resulted from errors in the pseudomolecule assembly. Altogether our data highlight extensive sequence exchanges in subtelomeres between non-homologous chromosomes in common bean and confirm that subtelomeres represent one of the most dynamic and rapidly evolving regions in eukaryotic genomes.

Genetic mapping of two genes conferring resistance to powdery mildew in common bean (Phaseolus vulgaris L.)

Powdery mildew (PM) is a serious disease in many legume species, including the common bean (Phaseolus vulgaris L.). This study investigated the genetic control behind resistance reaction to PM in the bean genotype, Cornell 49242. The results revealed evidence supporting a qualitative mode of inheritance for resistance and the involvement of two independent genes in the resistance reaction. The location of these resistance genes was investigated in a linkage genetic map developed for the XC RIL population. Contingency tests revealed significant associations for 28 loci out of a total of 329 mapped loci. Fifteen were isolated or formed groups with less than two loci. The thirteen remaining loci were located at three regions in linkage groups Pv04, Pv09, and Pv11. The involvement of Pv09 was discarded due to the observed segregation in the subpopulation obtained from the Xana genotype for the loci located in this region. In contrast, the two subpopulations obtained from the Xana genotype for the BM161 locus, linked to the Co-3/9 anthracnose resistance gene (Pv04), and from the Xana genotype for the SCAReoli locus, linked to the Co-2 anthracnose resistance gene (Pv11), exhibited monogenic segregations, suggesting that both regions were involved in the genetic control of resistance. A genetic dissection was carried out to verify the involvement of both regions in the reaction to PM. Two resistant recombinant lines were selected, according to their genotypes, for the block of loci included in the Co-2 and Co-3/9 regions, and they were crossed with the susceptible parent, Xana. Linkage analysis in the respective F2 populations supported the hypothesis that a dominant gene (Pm1) was located in the linkage group Pv11 and another gene (Pm2) was located in the linkage group Pv04. This is the first report showing the localization of resistance genes against powdery mildew in Phaseolus vulgaris and the results offer the opportunity to increase the efficiency of breeding programs by means of marker-assisted selection.

Evolution of a complex disease resistance gene cluster in diploid Phaseolus and tetraploid Glycine

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.

Pseudomonas syringae pv. phaseolicola effector HopF1 inhibits pathogen-associated molecular pattern-triggered immunity in a RIN4-independent manner in common bean (Phaseolus vulgaris)

Plant pathogens usually promote pathogenesis by secreting effector proteins into host plant cells. One of the secreted effectors of Pseudomonas syringae pv. phaseolicola, the causative agent of halo-blight disease in common bean (Phaseolus vulgaris), HopF1, activates effector-triggered immunity (ETI) in a bean cultivar containing R1 resistance gene, but displays virulence function in a bean cultivar without the R1 gene. The virulence mechanism of the effector remained unknown, although it was identified more than a decade ago. Here we demonstrated that HopF1 can inhibit pathogen-associated molecular pattern-triggered immunity (PTI) in a susceptible bean cultivar Tendergreen. HopF1 directly interacted with two RPM1-interacting protein 4 (RIN4) orthologs of bean, PvRIN4a and PvRIN4b. Like RIN4 in Arabidopsis, both PvRIN4 orthologs negatively regulated the PTI responses in bean. However, the virulence function of HopF1 was enhanced in Tendergreen silencing PvRIN4. Furthermore, silencing PvRIN4a compromised the avrβ1-induced hypersensitive response (HR), which previously was reported to be suppressed by HopF1. Together, these results demonstrated that PvRIN4 orthologs were not the virulence target of HopF1 for inhibiting PTI, but probably for interfering with ETI.