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.

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.

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.

Molecular analysis of a large subtelomeric nucleotide-binding-site-leucine-rich-repeat family in two representative genotypes of the major gene pools of Phaseolus vulgaris.

In common bean, the B4 disease resistance gene cluster is a complex cluster localized at the end of linkage group (LG) B4, containing at least three R specificities to the fungus Colletotrichum lindemuthianum. To investigate the evolution of this R cluster since the divergence of Andean and Mesoamerican gene pools, DNA sequences were characterized from two representative genotypes of the two major gene pools of common bean (BAT93: Mesoamerican; JaloEEP558: Andean). Sequences encoding 29 B4-CC nucleotide-binding-site-leucine-rich-repeat (B4-CNL) genes were determined-12 from JaloEEP558 and 17 from BAT93. Although sequence exchange events were identified, phylogenetic analyses revealed that they were not frequent enough to lead to homogenization of B4-CNL sequences within a haplotype. Genetic mapping based on pulsed-field gel electrophoresis separation confirmed that the B4-CNL family is a large family specific to one end of LG B4 and is present at two distinct blocks separated by 26 cM. Fluorescent in situ hybridization on meiotic pachytene chromosomes revealed that two B4-CNL blocks are located in the subtelomeric region of the short arm of chromosome 4 on both sides of a heterochromatic block (knob), suggesting that this peculiar genomic environment may favor the proliferation of a large R gene cluster.

Replication of nonautonomous retroelements in soybean appears to be both recent and common.

Retrotransposons and their remnants often constitute more than 50% of higher plant genomes. Although extensively studied in monocot crops such as maize (Zea mays) and rice (Oryza sativa), the impact of retrotransposons on dicot crop genomes is not well documented. Here, we present an analysis of retrotransposons in soybean (Glycine max). Analysis of approximately 3.7 megabases (Mb) of genomic sequence, including 0.87 Mb of pericentromeric sequence, uncovered 45 intact long terminal repeat (LTR)-retrotransposons. The ratio of intact elements to solo LTRs was 8:1, one of the highest reported to date in plants, suggesting that removal of retrotransposons by homologous recombination between LTRs is occurring more slowly in soybean than in previously characterized plant species. Analysis of paired LTR sequences uncovered a low frequency of deletions relative to base substitutions, indicating that removal of retrotransposon sequences by illegitimate recombination is also operating more slowly. Significantly, we identified three subfamilies of nonautonomous elements that have replicated in the recent past, suggesting that retrotransposition can be catalyzed in trans by autonomous elements elsewhere in the genome. Analysis of 1.6 Mb of sequence from Glycine tomentella, a wild perennial relative of soybean, uncovered 23 intact retroelements, two of which had accumulated no mutations in their LTRs, indicating very recent insertion. A similar pattern was found in 0.94 Mb of sequence from Phaseolus vulgaris (common bean). Thus, autonomous and nonautonomous retrotransposons appear to be both abundant and active in Glycine and Phaseolus. The impact of nonautonomous retrotransposon replication on genome size appears to be much greater than previously appreciated.

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.

BAC end sequences corresponding to the B4 resistance gene cluster in common bean: a resource for markers and synteny analyses.

In common bean, a complex disease resistance (R) gene cluster, harboring many specific R genes against various pathogens, is located at the end of the linkage group B4. A BAC library of the Meso-american bean genotype BAT93 was screened with PRLJ1, a probe previously shown to be specific to the B4 R gene cluster, leading to the identification of 73 positive BAC clones. BAC-end sequencing (BES) of the 73 positive BACs generated 75 kb of sequence. These BACs were organized into 6 contigs, all mapped at the B4 R gene cluster. To evaluate the potential of BES for marker development, BES-derived specific primers were used to check for linkage with two allelic anthracnose R specificities Co-3 and Co-3 ( 2 ), through the analysis of pairs of Near Isogenic Lines (NILs). Out of 32 primer pairs tested, two revealed polymorphisms between the NILs, confirming the suspected location of Co-3 and Co-3 ( 2 ) at the B4 cluster. In order to identify the orthologous region of the B4 R gene cluster in the two model legume genomes, bean BESs were used as queries in TBLASTX searches of Medicago truncatula and Lotus japonicus BAC clones. Putative orthologous regions were identified on chromosome Mt6 and Lj2, in agreement with the colinearity observed between Mt and Lj for these regions.

Resistance to Colletotrichum lindemuthianum in Phaseolus vulgaris: a case study for mapping two independent genes.

Anthracnose, caused by the hemibiotrophic fungal pathogen Colletotrichum lindemuthianum is a devastating disease of common bean. Resistant cultivars are economical means for defense against this pathogen. In the present study, we mapped resistance specificities against 7 C. lindemuthianum strains of various geographical origins revealing differential reactions on BAT93 and JaloEEP558, two parents of a recombinant inbred lines (RILs) population, of Meso-american and Andean origin, respectively. Six strains revealed the segregation of two independent resistance genes. A specific numerical code calculating the LOD score in the case of two independent segregating genes (i.e. genes with duplicate effects) in a RILs population was developed in order to provide a recombination value (r) between each of the two resistance genes and the tested marker. We mapped two closely linked Andean resistance genes (Co-x, Co-w) at the end of linkage group (LG) B1 and mapped one Meso-american resistance genes (Co-u) at the end of LG B2. We also confirmed the complexity of the previously identified B4 resistance gene cluster, because four of the seven tested strains revealed a resistance specificity near Co-y from JaloEEP558 and two strains identified a resistance specificity near Co-9 from BAT93. Resistance genes found within the same cluster confer resistance to different strains of a single pathogen such as the two anthracnose specificities Co-x and Co-w clustered at the end of LG B1. Clustering of resistance specificities to multiple pathogens such as fungi (Co-u) and viruses (I) was also observed at the end of LG B2.

Distinct post-transcriptional modifications result into seven alternative transcripts of the CC-NBS-LRR gene JA1tr of Phaseolus vulgaris.

The generation of splice variants has been reported for various plant resistance (R) genes, suggesting that these variants play an important role in disease resistance. Most of the time these R genes belong to the Toll and mammalian IL-1 receptor-nucleotide-binding site-leucine-rich repeat (TIR-NBS-LRR) class of R genes. In Phaseolus vulgaris, a resistance gene cluster (referred to as the B4 R-gene cluster) has been identified at the end of linkage group B4. At this complex resistance cluster, three R specificities (Co-9, Co-y and Co-z) and two R QTLs effective against the fungal pathogen Colletotrichum lindemuthianum, the causal agent of anthracnose, have been identified. At the molecular level, four resistance gene candidates encoding putative full-length, coiled-coil (CC)-NBS-LRR R-like proteins, with LRR numbers ranging from 18 to 20, have been previously characterized. In the present study, seven cDNA corresponding to truncated R-like transcripts, belonging to the CC-NBS-LRR class of plant disease R genes, have been identified. These seven transcripts correspond to a single gene named JA1tr, which encodes, at most, only five LRRs. The seven JA1tr transcript variants result from distinct post-transcriptional modifications of JA1tr, corresponding to alternative splicing events of two introns, exon skipping and multiple 'aberrant splicing' events in the open reading frame (ORF). JA1tr was mapped at the B4 R-gene cluster identified in common bean. These post-transcriptional modifications of the single gene JA1tr could constitute an efficient source of diversity. The present results provide one of the few reports of transcript variants with truncated ORFs resulting from a CC-NBS-LRR gene.