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.

A common bean truncated CRINKLY4 kinase controls gene-for-gene resistance to the fungus Colletotrichum lindemuthianum.

Identifying the molecular basis of resistance to pathogens is critical to promote a chemical-free cropping system. In plants, nucleotide-binding leucine-rich repeat constitute the largest family of disease resistance (R) genes, but this resistance can be rapidly overcome by the pathogen, prompting research into alternative sources of resistance. Anthracnose, caused by the fungus Colletotrichum lindemuthianum, is one of the most important diseases of common bean. This study aimed to identify the molecular basis of Co-x, an anthracnose R gene conferring total resistance to the extremely virulent C. lindemuthianum strain 100. To that end, we sequenced the Co-x 58 kb target region in the resistant JaloEEP558 (Co-x) common bean and identified KTR2/3, an additional gene encoding a truncated and chimeric CRINKLY4 kinase, located within a CRINKLY4 kinase cluster. The presence of KTR2/3 is strictly correlated with resistance to strain 100 in a diversity panel of common beans. Furthermore, KTR2/3 expression is up-regulated 24 hours post-inoculation and its transient expression in a susceptible genotype increases resistance to strain 100. Our results provide evidence that Co-x encodes a truncated and chimeric CRINKLY4 kinase probably resulting from an unequal recombination event that occurred recently in the Andean domesticated gene pool. This atypical R gene may act as a decoy involved in indirect recognition of a fungal effector.

Genome-Wide Identification of Key Components of RNA Silencing in Two Phaseolus vulgaris Genotypes of Contrasting Origin and Their Expression Analyses in Response to Fungal Infection.

RNA silencing serves key roles in a multitude of cellular processes, including development, stress responses, metabolism, and maintenance of genome integrity. Dicer, Argonaute (AGO), double-stranded RNA binding (DRB) proteins, RNA-dependent RNA polymerase (RDR), and DNA-dependent RNA polymerases known as Pol IV and Pol V form core components to trigger RNA silencing. Common bean (Phaseolus vulgaris) is an important staple crop worldwide. In this study, we aimed to unravel the components of the RNA-guided silencing pathway in this non-model plant, taking advantage of the availability of two genome assemblies of Andean and Meso-American origin. We identified six PvDCLs, thirteen PvAGOs, 10 PvDRBs, 5 PvRDRs, in both genotypes, suggesting no recent gene amplification or deletion after the gene pool separation. In addition, we identified one PvNRPD1 and one PvNRPE1 encoding the largest subunits of Pol IV and Pol V, respectively. These genes were categorized into subgroups based on phylogenetic analyses. Comprehensive analyses of gene structure, genomic localization, and similarity among these genes were performed. Their expression patterns were investigated by means of expression models in different organs using online data and quantitative RT-PCR after pathogen infection. Several of the candidate genes were up-regulated after infection with the fungus Colletotrichum lindemuthianum.

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.

Genomic and epigenomic immunity in common bean: the unusual features of NB-LRR gene family.

In plants, a key class of genes comprising most of disease resistance (R) genes encodes Nucleotide-binding leucine-rich repeat (NL) proteins. Access to common bean (Phaseolus vulgaris) genome sequence provides unparalleled insight into the organization and evolution of this large gene family (∼400 NL) in this important crop. As observed in other plant species, most common bean NL are organized in cluster of genes. However, a particularity of common bean is that these clusters are often located in subtelomeric regions close to terminal knobs containing the satellite DNA khipu. Phylogenetically related NL are spread between different chromosome ends, suggesting frequent exchanges between non-homologous chromosomes. NL peculiar location, in proximity to heterochromatic regions, led us to study their DNA methylation status using a whole-genome cytosine methylation map. In common bean, NL genes displayed an unusual body methylation pattern since half of them are methylated in the three contexts, reminiscent of the DNA methylation pattern of repeated sequences. Moreover, 90 NL were also abundantly targeted by 24 nt siRNA, with 90% corresponding to methylated NL genes. This suggests the existence of a transcriptional gene silencing mechanism of NL through the RdDM (RNA-directed DNA methylation) pathway in common bean that has not been described in other plant species.

Evolutionary dynamics of satellite DNA repeats from Phaseolus beans.

Common bean (Phaseolus vulgaris) subtelomeres are highly enriched for khipu, the main satellite DNA identified so far in this genome. Here, we comparatively investigate khipu genomic organization in Phaseolus species from different clades. Additionally, we identified and characterized another satellite repeat, named jumper, associated to khipu. A mixture of P. vulgaris khipu clones hybridized in situ confirmed the presence of khipu-like sequences on subterminal chromosome regions in all Phaseolus species, with differences in the number and intensity of signals between species and when species-specific clones were used. Khipu is present as multimers of ∼500 bp and sequence analyses of cloned fragments revealed close relationship among khipu repeats. The new repeat, named jumper, is a 170-bp satellite sequence present in all Phaseolus species and inserted into the nontranscribed spacer (NTS) of the 5S rDNA in the P. vulgaris genome. Nevertheless, jumper was found as a high-copy repeat at subtelomeres and/or pericentromeres in the Phaseolus microcarpus lineage only. Our data argue for khipu as an important subtelomeric satellite DNA in the genus and for a complex satellite repeat composition of P. microcarpus subtelomeres, which also contain jumper. Furthermore, the differential amplification of these repeats in subtelomeres or pericentromeres reinforces the presence of a dynamic satellite DNA library in Phaseolus.

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.

The "one-step" Bean pod mottle virus (BPMV)-derived vector is a functional genomics tool for efficient overexpression of heterologous protein, virus-induced gene silencing and genetic mapping of BPMV R-gene in common bean (Phaseolus vulgaris L.).

BACKGROUND: Over the last two years, considerable advances have been made in common bean (Phaseolus vulgaris L.) genomics, especially with the completion of the genome sequence and the availability of RNAseq data. However, as common bean is recalcitrant to stable genetic transformation, much work remains to be done for the development of functional genomics tools adapted to large-scale studies. RESULTS: Here we report the successful implementation of an efficient viral vector system for foreign gene expression, virus-induced gene silencing (VIGS) and genetic mapping of a BPMV resistance gene in common bean, using a "one-step" BPMV vector originally developed in soybean. With the goal of developing this vector for high-throughput VIGS studies in common bean, we optimized the conditions for rub-inoculation of infectious BPMV-derived plasmids in common bean cv. Black Valentine. We then tested the susceptibility to BPMV of six cultivars, and found that only Black Valentine and JaloEEP558 were susceptible to BPMV. We used a BPMV-GFP construct to detect the spatial and temporal infection patterns of BPMV in vegetative and reproductive tissues. VIGS of the PHYTOENE DESATURASE (PvPDS) marker gene was successfully achieved with recombinant BPMV vectors carrying fragments ranging from 132 to 391 bp. Finally, we mapped a gene for resistance to BPMV (R-BPMV) at one end of linkage group 2, in the vicinity of a locus (I locus) previously shown to be involved in virus resistance. CONCLUSIONS: The "one-step" BPMV vector system therefore enables rapid and simple functional studies in common bean, and could be suitable for large-scale analyses. In the post-genomic era, these advances are timely for the common bean research community.

Differential gene expression and metabolomic analyses of Brachypodium distachyon infected by deoxynivalenol producing and non-producing strains of Fusarium graminearum.

BACKGROUND: Fusarium Head Blight (FHB) caused primarily by Fusarium graminearum (Fg) is one of the major diseases of small-grain cereals including bread wheat. This disease both reduces yields and causes quality losses due to the production of deoxynivalenol (DON), the major type B trichothecene mycotoxin. DON has been described as a virulence factor enabling efficient colonization of spikes by the fungus in wheat, but its precise role during the infection process is still elusive. Brachypodium distachyon (Bd) is a model cereal species which has been shown to be susceptible to FHB. Here, a functional genomics approach was performed in order to characterize the responses of Bd to Fg infection using a global transcriptional and metabolomic profiling of B. distachyon plants infected by two strains of F. graminearum: a wild-type strain producing DON (Fgdon+) and a mutant strain impaired in the production of the mycotoxin (Fgdon-). RESULTS: Histological analysis of the interaction of the Bd21 ecotype with both Fg strains showed extensive fungal tissue colonization with the Fgdon+ strain while the florets infected with the Fgdon- strain exhibited a reduced hyphal extension and cell death on palea and lemma tissues. Fungal biomass was reduced in spikes inoculated with the Fgdon- strain as compared with the wild-type strain. The transcriptional analysis showed that jasmonate and ethylene-signalling pathways are induced upon infection, together with genes encoding putative detoxification and transport proteins, antioxidant functions as well as secondary metabolite pathways. In particular, our metabolite profiling analysis showed that tryptophan-derived metabolites, tryptamine, serotonin, coumaroyl-serotonin and feruloyl-serotonin, are more induced upon infection by the Fgdon+ strain than by the Fgdon- strain. Serotonin was shown to exhibit a slight direct antimicrobial effect against Fg. CONCLUSION: Our results show that Bd exhibits defense hallmarks similar to those already identified in cereal crops. While the fungus uses DON as a virulence factor, the host plant preferentially induces detoxification and the phenylpropanoid and phenolamide pathways as resistance mechanisms. Together with its amenability in laboratory conditions, this makes Bd a very good model to study cereal resistance mechanisms towards the major disease FHB.

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.

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

KEY MESSAGE: 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.

The secondary metabolism glycosyltransferases UGT73B3 and UGT73B5 are components of redox status in resistance of Arabidopsis to Pseudomonas syringae pv. tomato.

Secondary metabolism plant glycosyltransferases (UGTs) ensure conjugation of sugar moieties to secondary metabolites (SMs) and glycosylation contributes to the great diversity, reactivity and regulation of SMs. UGT73B3 and UGT73B5, two UGTs of Arabidopsis thaliana (Arabidopsis), are involved in the hypersensitive response (HR) to the avirulent bacteria Pseudomonas syringae pv. tomato (Pst-AvrRpm1), but their function in planta is unknown. Here, we report that ugt73b3, ugt73b5 and ugt73b3 ugt73b5 T-DNA insertion mutants exhibited an accumulation of reactive oxygen species (ROS), an enhanced cell death during the HR to Pst-AvrRpm1, whereas glutathione levels increased in the single mutants. In silico analyses indicate that UGT73B3 and UGT73B5 belong to the early salicylic acid (SA)-induced genes whose pathogen-induced expression is co-regulated with genes related to cellular redox homeostasis and general detoxification. Analyses of metabolic alterations in ugt mutants reveal modification of SA and scopoletin contents which correlate with redox perturbation, and indicate quantitative modifications in the pattern of tryptophan-derived SM accumulation after Pst-AvrRpm1 inoculation. Our data suggest that UGT73B3 and UGT73B5 participate in regulation of redox status and general detoxification of ROS-reactive SMs during the HR to Pst-AvrRpm1, and that decreased resistance to Pst-AvrRpm1 in ugt mutants is tightly linked to redox perturbation.

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.

Identification and characterization of functional centromeres of the common bean.

In higher eukaryotes, centromeres are typically composed of megabase-sized arrays of satellite repeats that evolve rapidly and homogenize within a species' genome. Despite the importance of centromeres, our knowledge is limited to a few model species. We conducted a comprehensive analysis of common bean (Phaseolus vulgaris) centromeric satellite DNA using genomic data, fluorescence in situ hybridization (FISH), immunofluorescence and chromatin immunoprecipitation (ChIP). Two unrelated centromere-specific satellite repeats, CentPv1 and CentPv2, and the common bean centromere-specific histone H3 (PvCENH3) were identified. FISH showed that CentPv1 and CentPv2 are predominantly located at subsets of eight and three centromeres, respectively. Immunofluorescence- and ChIP-based assays demonstrated the functional significance of CentPv1 and CentPv2 at centromeres. Genomic analysis revealed several interesting features of CentPv1 and CentPv2: (i) CentPv1 is organized into an higher-order repeat structure, named Nazca, of 528 bp, whereas CentPv2 is composed of tandemly organized monomers; (ii) CentPv1 and CentPv2 have undergone chromosome-specific homogenization; and (iii) CentPv1 and CentPv2 are not likely to be commingled in the genome. These findings suggest that two distinct sets of centromere sequences have evolved independently within the common bean genome, and provide insight into centromere satellite evolution.

Gene duplication within the Green Lineage: the case of TEL genes.

Recent years have witnessed a breathtaking increase in the availability of genome sequence data, providing evidence of the highly duplicate nature of eukaryotic genomes. Plants are exceptional among eukaryotic organisms in that duplicate loci compose a large fraction of their genomes, partly because of the frequent occurrence of polyploidy (or whole-genome duplication) events. Tandem gene duplication and transposition have also contributed to the large number of duplicated genes in plant genomes. Evolutionary analyses allowed the dynamics of duplicate gene evolution to be studied and several models were proposed. It seems that, over time, many duplicated genes were lost and some of those that were retained gained new functions and/or expression patterns (neofunctionalization) or subdivided their functions and/or expression patterns between them (subfunctionalization). Recent studies have provided examples of genes that originated by duplication with successive diversification within plants. In this review, we focused on the TEL (TERMINAL EAR1-like) genes to illustrate such mechanisms. Emerged from the mei2 gene family, these TEL genes are likely to be land plant-specific. Phylogenetic analyses revealed one or two TEL copies per diploid genome. TEL gene degeneration and loss in several Angiosperm species such as in poplar and maize seem to have occurred. In Arabidopsis thaliana, whose genome experienced at least three polyploidy events followed by massive gene loss and genomic reorganization, two TEL genes were retained and two new shorter TEL-like (MCT) genes emerged. Molecular and expression analyses suggest for these genes sub- and neofunctionalization events, but confirmation will come from their functional characterization.

The function of the RNA-binding protein TEL1 in moss reveals ancient regulatory mechanisms of shoot development.

The shoot represents the basic body plan in land plants. It consists of a repeated structure composed of stems and leaves. Whereas vascular plants generate a shoot in their diploid phase, non-vascular plants such as mosses form a shoot (called the gametophore) in their haploid generation. The evolution of regulatory mechanisms or genetic networks used in the development of these two kinds of shoots is unclear. TERMINAL EAR1-like genes have been involved in diploid shoot development in vascular plants. Here, we show that disruption of PpTEL1 from the moss Physcomitrella patens, causes reduced protonema growth and gametophore initiation, as well as defects in gametophore development. Leafy shoots formed on ΔTEL1 mutants exhibit shorter stems with more leaves per shoot, suggesting an accelerated leaf initiation (shortened plastochron), a phenotype shared with the Poaceae vascular plants TE1 and PLA2/LHD2 mutants. Moreover, the positive correlation between plastochron length and leaf size observed in ΔTEL1 mutants suggests a conserved compensatory mechanism correlating leaf growth and leaf initiation rate that would minimize overall changes in plant biomass. The RNA-binding protein encoded by PpTEL1 contains two N-terminus RNA-recognition motifs, and a third C-terminus non-canonical RRM, specific to TEL proteins. Removal of the PpTEL1 C-terminus (including this third RRM) or only 16-18 amino acids within it seriously impairs PpTEL1 function, suggesting a critical role for this third RRM. These results show a conserved function of the RNA-binding PpTEL1 protein in the regulation of shoot development, from early ancestors to vascular plants, that depends on the third TEL-specific RRM.

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

*In plants, the evolution of specific resistance is poorly understood. Pseudomonas syringae effectors AvrB and AvrRpm1 are recognized by phylogenetically distinct resistance (R) proteins in Arabidopsis thaliana (Brassicaceae) and soybean (Glycine max, Fabaceae). In soybean, these resistances are encoded by two tightly linked R genes, Rpg1-b and Rpg1-r. To study the evolution of these specific resistances, we investigated AvrB- and AvrRpm1-induced responses in common bean (Phaseolus vulgaris, Fabaceae). *Common bean genotypes of various geographical origins were inoculated with P. syringae strains expressing AvrB or AvrRpm1. A common bean recombinant inbred line (RIL) population was used to map R genes to AvrRpm1. *No common bean genotypes recognized AvrB. By contrast, multiple genotypes responded to AvrRpm1, and two independent R genes conferring AvrRpm1-specific resistance were mapped to the ends of linkage group B11 (Rpsar-1, for resistance to Pseudomonas syringae effector AvrRpm1 number 1) and B8 (Rpsar-2). Rpsar-1 is located in a region syntenic with the soybean Rpg1 cluster. However, mapping of specific Rpg1 homologous genes suggests that AvrRpm1 recognition evolved independently in common bean and soybean. *The conservation of the genomic position of AvrRpm1-specific genes between soybean and common bean suggests a model whereby specific clusters of R genes are predisposed to evolve recognition of the same effector molecules.

Three highly similar formate dehydrogenase genes located in the vicinity of the B4 resistance gene cluster are differentially expressed under biotic and abiotic stresses in Phaseolus vulgaris.

In higher plants, formate dehydrogenase (FDH, EC1.2.1.2.) catalyzes the NAD-linked oxidation of formate to CO(2), and FDH transcript accumulation has been reported after various abiotic stresses. By sequencing a Phaseolus vulgaris BAC clone encompassing a CC-NBS-LRR gene rich region of the B4 resistance gene cluster, we identified three FDH-encoding genes. FDH is present as a single copy gene in the Arabidopsis thaliana genome, and public database searches confirm that FDH is a low copy gene in plant genomes, since only 33 FDH homologs were identified from 27 plant species. Three independent prediction programs (Predotar, TargetP and Mitoprot) used on this large subset of 33 plant FDHs, revealed that mitochondrial localization of FDH might be the rule in higher plants. A phylogenetic analysis suggests a scenario of local FDH gene duplication in an ancestor of the Phaseoleae followed by another more recent duplication event after bean/soybean divergence. The expression levels of two common bean FDH genes under different treatments were investigated by quantitative RT-PCR analysis. FDH genes are differentially up-regulated after biotic and abiotic stresses (infection with the fungus Colletotrichum lindemuthianum, and dark treatment, respectively). The present study provides the first report of FDH transcript accumulation after biotic stress, suggesting the involvement of FDH in the pathogen resistance process.

A nomadic subtelomeric disease resistance gene cluster in common bean.

The B4 resistance (R) gene cluster is one of the largest clusters known in common bean (Phaseolus vulgaris [Pv]). It is located in a peculiar genomic environment in the subtelomeric region of the short arm of chromosome 4, adjacent to two heterochromatic blocks (knobs). We sequenced 650 kb spanning this locus and annotated 97 genes, 26 of which correspond to Coiled-Coil-Nucleotide-Binding-Site-Leucine-Rich-Repeat (CNL). Conserved microsynteny was observed between the Pv B4 locus and corresponding regions of Medicago truncatula and Lotus japonicus in chromosomes Mt6 and Lj2, respectively. The notable exception was the CNL sequences, which were completely absent in these regions. The origin of the Pv B4-CNL sequences was investigated through phylogenetic analysis, which reveals that, in the Pv genome, paralogous CNL genes are shared among nonhomologous chromosomes (4 and 11). Together, our results suggest that Pv B4-CNL was derived from CNL sequences from another cluster, the Co-2 cluster, through an ectopic recombination event. Integration of the soybean (Glycine max) genome data enables us to date more precisely this event and also to infer that a single CNL moved from the Co-2 to the B4 cluster. Moreover, we identified a new 528-bp satellite repeat, referred to as khipu, specific to the Phaseolus genus, present both between B4-CNL sequences and in the two knobs identified at the B4 R gene cluster. The khipu repeat is present on most chromosomal termini, indicating the existence of frequent ectopic recombination events in Pv subtelomeric regions. Our results highlight the importance of ectopic recombination in R gene evolution.

Differential accumulation of retroelements and diversification of NB-LRR disease resistance genes in duplicated regions following polyploidy in the ancestor of soybean.

The genomes of most, if not all, flowering plants have undergone whole genome duplication events during their evolution. The impact of such polyploidy events is poorly understood, as is the fate of most duplicated genes. We sequenced an approximately 1 million-bp region in soybean (Glycine max) centered on the Rpg1-b disease resistance gene and compared this region with a region duplicated 10 to 14 million years ago. These two regions were also compared with homologous regions in several related legume species (a second soybean genotype, Glycine tomentella, Phaseolus vulgaris, and Medicago truncatula), which enabled us to determine how each of the duplicated regions (homoeologues) in soybean has changed following polyploidy. The biggest change was in retroelement content, with homoeologue 2 having expanded to 3-fold the size of homoeologue 1. Despite this accumulation of retroelements, over 77% of the duplicated low-copy genes have been retained in the same order and appear to be functional. This finding contrasts with recent analyses of the maize (Zea mays) genome, in which only about one-third of duplicated genes appear to have been retained over a similar time period. Fluorescent in situ hybridization revealed that the homoeologue 2 region is located very near a centromere. Thus, pericentromeric localization, per se, does not result in a high rate of gene inactivation, despite greatly accelerated retrotransposon accumulation. In contrast to low-copy genes, nucleotide-binding-leucine-rich repeat disease resistance gene clusters have undergone dramatic species/homoeologue-specific duplications and losses, with some evidence for partitioning of subfamilies between homoeologues.