Perturbation of nuclear-cytosolic shuttling of Rx1 compromises extreme resistance and translational arrest of potato virus X transcripts.

Many plant intracellular immune receptors mount a hypersensitive response (HR) upon pathogen perception. The concomitant localized cell death is proposed to trap pathogens, such as viruses, inside infected cells, thereby preventing their spread. Notably, extreme resistance (ER) conferred by the potato immune receptor Rx1 to potato virus X (PVX) does not involve the death of infected cells. It is unknown what defines ER and how it differs from HR-based resistance. Interestingly, Rx1 can trigger an HR, but only upon artificial (over)expression of PVX or its avirulence coat protein (CP). Rx1 has a nucleocytoplasmic distribution and both pools are required for HR upon transient expression of a PVX-GFP amplicon. It is unknown whether mislocalized Rx1 variants can induce ER upon natural PVX infection. Here, we generated transgenic Nicotiana benthamiana producing nuclear- or cytosol-restricted Rx1 variants. We found that these variants can still mount an HR. However, nuclear- or cytosol-restricted Rx1 variants can no longer trigger ER or restricts viral infection. Interestingly, unlike the mislocalized Rx1 variants, wild-type Rx1 was found to compromise CP protein accumulation. We show that the lack of CP accumulation does not result from its degradation but is likely to be linked with translational arrest of its mRNA. Together, our findings suggest that translational arrest of viral genes is a major component of ER and, unlike the HR, is required for resistance to PVX.

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.

Bean pod mottle virus: a new powerful tool for functional genomics studies in Pisum sativum.

Pea (Pisum sativum L.) is an important legume worldwide. The importance of pea in arable rotations and nutritional value for both human and animal consumption have fostered sustained production and different studies to improve agronomic traits of interest. Moreover, complete sequencing of the pea genome is currently underway and will lead to the identification of a large number of genes potentially associated with important agronomic traits. Because stable genetic transformation is laborious for pea, virus-induced gene silencing (VIGS) appears as a powerful alternative technology for determining the function of unknown genes. In this work, we present a rapid and efficient viral inoculation method using DNA infectious plasmids of Bean pod mottle virus (BPMV)-derived VIGS vector. Six pea genotypes with important genes controlling biotic and/or abiotic stresses were found susceptible to BPMV carrying a GFP reporter gene and showed fluorescence in both shoots and roots. In a second step, we investigated 37 additional pea genotypes and found that 30 were susceptible to BPMV and only 7 were resistant. The capacity of BPMV to induce silencing of endogenes was investigated in the most susceptible genotype using two visual reporter genes: PsPDS and PsKORRIGAN1 (PsKOR1) encoding PHYTOENE DESATURASE and a 1,4-ß-D-glucanase, respectively. The features of the 'one-step' BPMV-derived VIGS vector include (i) the ease of rub-inoculation, without any need for biolistic or agro-inoculation procedures, (ii) simple cost-effective procedure and (iii) noninterference of viral symptoms with silencing. These features make BPMV the most adapted VIGS vector in pea to make low- to high-throughput VIGS studies.

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.

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.

Epigenetic analyses and the distribution of repetitive DNA and resistance genes reveal the complexity of common bean (Phaseolus vulgaris L., Fabaceae) heterochromatin.

The common bean (Phaseolus vulgaris L.) is the main representative of its genus and one of most important sources of proteins in African and Latin American countries. Although it is a species with a small genome, its pericentromeric and subtelomeric heterochromatin fractions are interspersed with single-copy sequences and active genes, suggesting a less compartmentalized genome organization. The present study characterized its chromatin fractions, associating the distribution of repetitive sequences and resistance genes with histone and DNA epigenetic modifications with and without biotic stress. Immunostaining with H3K4me3 and H4K5ac were generally associated with euchromatic regions, whereas H3K9me2, H3K27me1, and 5mC preferentially labeled the pericentromeric heterochromatin. The 45S rDNA and centromeric DNA sequences were hypomethylated as were most of the terminal heterochromatic blocks. The largest of them, which is associated with resistance genes, was also hypomethylated after the plants were infected with virulent and avirulent strains of the fungus Colletotrichum lindemuthianum, suggesting no correlation with control of resistance gene expression. The results highlighted the differences between subtelomeric and pericentromeric heterochromatin as well as variation within the pericentromeric heterochromatin.

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.

Hydathode immunity protects the Arabidopsis leaf vasculature against colonization by bacterial pathogens.

Plants prevent disease by passively and actively protecting potential entry routes against invading microbes. For example, the plant immune system actively guards roots, wounds, and stomata. How plants prevent vascular disease upon bacterial entry via guttation fluids excreted from specialized glands at the leaf margin remains largely unknown. These so-called hydathodes release xylem sap when root pressure is too high. By studying hydathode colonization by both hydathode-adapted (Xanthomonas campestris pv. campestris) and non-adapted pathogenic bacteria (Pseudomonas syringae pv. tomato) in immunocompromised Arabidopsis mutants, we show that the immune hubs BAK1 and EDS1-PAD4-ADR1 restrict bacterial multiplication in hydathodes. Both immune hubs effectively confine bacterial pathogens to hydathodes and lower the number of successful escape events of an hydathode-adapted pathogen toward the xylem. A second layer of defense, which is dependent on the plant hormones' pipecolic acid and to a lesser extent on salicylic acid, reduces the vascular spread of the pathogen. Thus, besides glands, hydathodes represent a potent first line of defense against leaf-invading microbes.

Members of the ribosomal protein S6 (RPS6) family act as pro-viral factor for tomato spotted wilt orthotospovirus infectivity in Nicotiana benthamiana.

To identify host factors for tomato spotted wilt orthotospovirus (TSWV), a virus-induced gene silencing (VIGS) screen using tobacco rattle virus (TRV) was performed on Nicotiana benthamiana for TSWV susceptibility. To rule out any negative effect on the plants' performance due to a double viral infection, the method was optimized to allow screening of hundreds of clones in a standardized fashion. To normalize the results obtained in and between experiments, a set of controls was developed to evaluate in a consist manner both VIGS efficacy and the level of TSWV resistance. Using this method, 4532 random clones of an N. benthamiana cDNA library were tested, resulting in five TRV clones that provided nearly complete resistance against TSWV. Here we report on one of these clones, of which the insert targets a small gene family coding for the ribosomal protein S6 (RPS6) that is part of the 40S ribosomal subunit. This RPS6 family is represented by three gene clades in the genome of Solanaceae family members, which were jointly important for TSWV susceptibility. Interestingly, RPS6 is a known host factor implicated in the replication of different plant RNA viruses, including the negative-stranded TSWV and the positive-stranded potato virus X.

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.

Unlike Many Disease Resistances, Rx1-Mediated Immunity to Potato Virus X Is Not Compromised at Elevated Temperatures.

Specificity in the plant immune system is mediated by Resistance (R) proteins. Most R genes encode intracellular NLR-type immune receptors and these pathogen sensors require helper NLRs to activate immune signaling upon pathogen perception. Resistance conferred by many R genes is temperature sensitive and compromised above 28°C. Many Solanaceae R genes, including the potato NLR Rx1 conferring resistance to Potato Virus X (PVX), have been reported to be temperature labile. Rx1 activity, like many Solanaceae NLRs, depends on helper-NLRs called NRC's. In this study, we investigated Rx1 resistance at elevated temperatures in potato and in Nicotiana benthamiana plants stably expressing Rx1 upon rub-inoculation with GFP-expressing PVX particles. In parallel, we used susceptible plants as a control to assess infectiousness of PVX at a range of different temperatures. Surprisingly, we found that Rx1 confers virus resistance in N. benthamiana up to 32°C, a temperature at which the PVX::GFP lost infectiousness. Furthermore, at 34°C, an Rx1-mediated hypersensitive response could still be triggered in N. benthamiana upon PVX Coat-Protein overexpression. As the Rx1-immune signaling pathway is not temperature compromised, this implies that at least one N. benthamiana helper NRC and its downstream signaling components are temperature tolerant. This finding suggests that the temperature sensitivity for Solanaceous resistances is likely attributable to the sensor NLR and not to its downstream signaling components.

Activation of immune receptor Rx1 triggers distinct immune responses culminating in cell death after 4 hours.

Intracellular nucleotide-binding leucine-rich repeat (NLR)-type immune receptors are a fundamental part of plant immune systems. As infection occurs at foci, activation of immune responses is typically non-uniform and non-synchronized, hampering the systematic dissection of their cellular effects and determining their phasing. We investigated the potato NLR Rx1 using the CESSNA (Controlled Expression of effectors for Synchronized and Systemic NLR Activation) platform. CESSNA-mediated Potato virus X coat protein (CP) expression allowed the monitoring of Rx1-mediated immune responses in a quantitative and reproducible manner. Rx1 was found to trigger a reactive oxygen species (ROS) burst and ion leakage within 1 h and a change in autofluorescence within 2 h after the induction of CP production. After 2 h, HIN1 expression was increased and single-stranded DNA (ssDNA) damage and loss of cellular integrity became apparent, followed by double-stranded DNA (dsDNA) damage after 3 h and increased PR-1a, LOX, ERF1 and AOX1B expression and cell death at 4 h. Nuclear exclusion of Rx1 resulted in increased basal levels of ROS and permitted Rx1 activation by an Rx1-breaking CP variant. In contrast, nuclear-targeted Rx1 showed diminished basal ROS levels, and only avirulent CP could trigger a compromised ROS production. Both nuclear-excluded and nuclear-targeted Rx1 triggered a delayed ion leakage compared with non-modified Rx1, suggesting that ion leakage and ROS production originate from distinct signalling pathways. This work offers novel insights into the influence of Rx1 localization on its activity, and the interplay between Rx1-triggered processes.

Molecular mechanisms that limit the costs of NLR-mediated resistance in plants.

Crop diseases cause significant yield losses, and the use of resistant cultivars can effectively mitigate these losses and control many plant diseases. Most plant resistance (R) genes encode immune receptors composed of nucleotide-binding and leucine-rich repeat (NLR) domains. These proteins mediate the specific recognition of pathogen avirulence effectors to induce defence responses. However, NLR-triggered immunity can be associated with a reduction in growth and yield, so-called 'fitness costs'. Recent data have shown that plants use an elaborate interplay of different mechanisms to control NLR gene transcript levels, as well as NLR protein abundance and activity, to avoid the associated cost of resistance in the absence of a pathogen. In this review, we discuss the different levels of NLR regulation (transcriptional, post-transcriptional and at the protein level). We address the apparent need for plants to maintain diverse modes of regulation. A recent model suggesting an equilibrium 'ON/OFF state' of NLR proteins, in the absence of a pathogen, provides the context for our discussion.

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.

Plant Autoimmunity: When Good Things Go Bad.

A recent study finds that the Arabidopsis DM1 and DM2d proteins physically interact and trigger autoimmunity in plants. The DM1-DM2d interaction pattern differs from that of known immune receptor pairs, portraying the versatility in NLR functioning.

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.

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.

VIGS technology: an attractive tool for functional genomics studies in legumes.

Legume species are among the most important crops worldwide. In recent years, six legume genomes have been completely sequenced, and there is now an urgent need for reverse-genetics tools to validate genes affecting yield and product quality. As most legumes are recalcitrant to stable genetic transformation, virus-induced gene silencing (VIGS) appears to be a powerful alternative technology for determining the function of unknown genes. VIGS technology is based on the property of plant viruses to trigger a defence mechanism related to post-transcriptional gene silencing (PTGS). Infection by a recombinant virus carrying a fragment of a plant target gene will induce homology-dependent silencing of the endogenous target gene. Several VIGS systems have been developed for legume species since 2004, including those based on Bean pod mottle virus, Pea early browning virus, and Apple latent spherical virus, and used in reverse-genetics studies of a wide variety of plant biological processes. In this work, we give an overview of the VIGS systems available for legumes, and present their successful applications in functional genomics studies. We also discuss the limitations of these VIGS systems and the future challenges to be faced in order to use VIGS to its full potential in legume species.

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.