Monica Rodriguez

16S rDNA is the historical gold standard for bacterial identification, particularly in metabarcoding approaches reliant on sequence similarity thresholds. We analyzed 6,660 Rhizobium 16S rRNA gene sequences from GenBank to examine the relationship between sequence identity and three metadata: species name, host plant, and geographic origin. Using an iterative BLAST-based pipeline, we detected 116,069 pairwise matches and assessed concordance among sequences (average length 1,328 bp) sharing 100% identity. For those in which the organism name, host plant and country of isolation were present in the record, surprisingly, 66.59% of identical sequence pairs showed full discordance across all three metadata, while only 1.40% shared the same name, host, and country. The most widespread sequence, detected 371 times, was associated with over 56 different host plants across 25 countries and bore multiple species name designations. These results highlight a striking mismatch between the 16S barcode and the taxonomic, ecological, and phenotypic variability it is assumed to reflect, likely arising from the slow evolution of rRNA genes contrasted with the mobility of ecologically relevant genes via horizontal transfer on plasmids, transposons, and phages. Our findings further challenge the limitations of relying on 16S rRNA alone for fine-scale taxonomic and metadata-based inference in capturing the true functional and ecological diversity of bacteria, endorsing the critical importance of polyphasic taxonomic approaches that integrate genomic, phenotypic, and ecological data. An interesting byproduct of the analysis was to realize the possibility of treating these data as if they were ‘citations.’ The more one finds the same query sequence, the more that sequence can be considered biologically ‘cited’, i.e., re-proposed elsewhere in the world. Thus, one can also analyze the h-index of such a ranking. In our Rhizobium dataset, we calculated an h-index = 201, meaning the sequence ranked 201 ^st had 202 identical homologues in GenBank. Although the research effort on given species is directly connected with it, this number provides a quantitative indicator of a taxon’s sequence recurrence and distribution within public databases, independent of nomenclatural inconsistencies, offering a novel framework for assessing bacterial representation across global datasets.

Main conclusion: The genetic architecture underlying co-infection with Meloidogyne incognita and Fusarium oxysporum f. sp. phaseoli differs from that observed when these pathogens are evaluated individually. Common bean ( Phaseolus vulgaris L.) yield is threatened by simultaneous infection with the root-knot nematode Meloidogyne incognita (Mi) and Fusarium oxysporum f. sp. phaseoli (Fop). Root-knot nematodes are believed to intensify the severity of Fusarium wilt in common bean, and specific genomic regions are understood to be associated with the host response, whether conferring susceptibility or resistance. To elucidate the genetic mechanisms of this interaction, phenotypic traits were evaluated in greenhouse trials, followed by associative mapping using a genome-wide association study (GWAS) approach. The plant material consisted of a core collection of 180 common bean genotypes from the Agronomic Institute (IAC, Campinas, Brazil) diversity panel. The effects of Fop and Mi were evaluated individually and in co-infection. Associative mapping was performed using the Bayesian information and linkage disequilibrium iteratively nested keyway (BLINK) model. When plants were infected with Mi, significant SNPs were detected on chromosomes Pv07, Pv08, and Pv10 based on gall counts. SNPs were detected on Pv05, Pv06, P10, and Pv11 in association with co-infection. Regions associated with egg mass count were detected on Pv02, Pv04, and Pv05. However, co-infection revealed SNPs on chromosomes P10 and Pv11. Three SNPs were associated with Fusarium wilt–two on Pv07 and one on Pv08. The genomic regions and markers associated with resistance to Mi and Fop provide new resources for advancing understanding of host–pathogen relationships in these important pathosystems.

The common bean (Phaseolus vulgaris L.) is a crucial legume crop and an ideal evolutionary model to study adaptive diversity in wild and domesticated populations. Here, we present a common bean pan-genome based on five high-quality genomes and whole-genome reads representing 339 genotypes. It reveals ~234 Mb of additional sequences containing 6,905 protein-coding genes missing from the reference, constituting 49% of all presence/absence variants (PAVs). More non-synonymous mutations are found in PAVs than core genes, probably reflecting the lower effective population size of PAVs and fitness advantages due to the purging effect of gene loss. Our results suggest pan-genome shrinkage occurred during wild range expansion. Selection signatures provide evidence that partial or complete gene loss was a key adaptive genetic change in common bean populations with major implications for plant adaptation. The pan-genome is a valuable resource for food legume research and breeding for climate change mitigation and sustainable agriculture.

Domesticated crops have been disseminated by humans over vast geographic areas. Common bean (Phaseolus vulgaris L.) was introduced in Europe after 1492. Here, by combining whole-genome profiling, metabolic fingerprinting and phenotypic characterisation, we show that the first common bean cultigens successfully introduced into Europe were of Andean origin, after Francisco Pizarro’s expedition to northern Peru in 1529. We reveal that hybridisation, selection and recombination have shaped the genomic diversity of the European common bean in parallel with political constraints. There is clear evidence of adaptive introgression into the Mesoamerican-derived European genotypes, with 44 Andean introgressed genomic segments shared by more than 90% of European accessions and distributed across all chromosomes except PvChr11. Genomic scans for signatures of selection highlight the role of genes relevant to flowering and environmental adaptation, suggesting that introgression has been crucial for the dissemination of this tropical crop to the temperate regions of Europe.

The knowledge of the organization of the domesticated gene pool of crop species is an essential requirement to understand crop evolution, to rationalize conservation programs, and to support practical decisions in plant breeding. Here, we integrate simple sequence repeat (SSR) analysis and phenotypic characterization to investigate a globe artichoke collection that comprises most of the varieties cultivated worldwide. We show that the cultivated gene pool of globe artichoke includes five distinct genetic groups associated with the major phenotypic typologies: Catanesi (which based on our analysis corresponds to Violetti di Provenza), Spinosi, Violetti di Toscana, Romaneschi, and Macau. We observed that 17 and 11% of the molecular and phenotypic variance, respectively, is between these groups, while within groups, strong linkage disequilibrium and heterozygote excess are evident. The divergence between groups for quantitative traits correlates with the average broad-sense heritability within the groups. The phenotypic divergence between groups for both qualitative and quantitative traits is strongly and positively correlated with SSR divergence (F_ST) between groups. All this implies a low population size and strong bottleneck effects, and indicates a long history of clonal propagation and selection during the evolution of the domesticated gene pool of globe artichoke. Moreover, the comparison between molecular and phenotypic population structures suggests that harvest time, plant architecture (i.e., plant height, stem length), leaf spininess, head morphology (i.e., head shape, bract shape, spininess) together with the number of heads per plant were the main targets of selection during the evolution of the cultivated germplasm. We emphasize our findings in light of the potential exploitation of this collection for association mapping studies.

The valorization of plant genetic resources and their direct use in local markets can make a significant contribution to the preservation of agrobiodiversity, while also contributing to the sustainability of rural communities. Indeed, plant genetic resources are a precious source of genes, and they represent an important crop heritage for the quality and sensory characteristics that are required by both farmers and consumers. However, an efficient strategy of agrobiodiversity conservation is strictly connected to product marketability and to consumer preferences. In the present study, choice experiments that involved 920 consumers were carried out to determine their willingness to pay for ancient local tomato varieties (landraces) rather than commercial varieties based on their preferences, and to determine how much they valued these products. The results obtained indicate that consumers are willing to pay premium prices for ancient local tomato varieties (an additional €0.90 kg−1), thus demonstrating their increasing attention to sustainable food and the willingness to contribute to agrobiodiversity conservation and enhancement. These results provide the basis for planning strategies and programs to support the cultivation of these landraces and the development of regional and national markets to acknowledge their characteristics, which will considerably increase the effectiveness and efficiency of conservation strategies.

In legumes, pod shattering occurs when mature pods dehisce along the sutures, and detachment of the valves promotes seed dispersal. In Phaseolus vulgaris (L)., the major locus qPD5.1-Pv for pod indehiscence was identified recently. We developed a BC₄/F₄ introgression line population and narrowed the major locus down to a 22.5 kb region. Here, gene expression and a parallel histological analysis of dehiscent and indehiscent pods identified an AtMYB26 orthologue as the best candidate for loss of pod shattering, on a genomic region ~11 kb downstream of the highest associated peak. Based on mapping and expression data, we propose early and fine up-regulation of PvMYB26 in dehiscent pods. Detailed histological analysis establishes that pod indehiscence is associated with the lack of a functional abscission layer in the ventral sheath, and that the key anatomical modifications associated with pod shattering in common bean occur early during pod development. We finally propose that loss of pod shattering in legumes resulted from histological convergent evolution and that it is the result of selection at orthologous loci.

Tomato (Solanum lycopersicum L.) is a widely used model plant species for dissecting out the genomic bases of complex traits to thus provide an optimal platform for modern “-omics” studies and genome-guided breeding. Genome-wide association studies (GWAS) have become a preferred approach for screening large diverse populations and many traits. Here, we present GWAS analysis of a collection of 115 landraces and 11 vintage and modern cultivars. A total of 26 conventional descriptors, 40 traits obtained by digital phenotyping, the fruit content of six carotenoids recorded at the early ripening (breaker) and red-ripe stages and 21 climate-related variables were analyzed in the context of genetic diversity monitored in the 126 accessions. The data obtained from thorough phenotyping and the SNP diversity revealed by sequencing of ripe fruit transcripts of 120 of the tomato accessions were jointly analyzed to determine which genomic regions are implicated in the expressed phenotypic variation. This study reveals that the use of fruit RNA-Seq SNP diversity is effective not only for identification of genomic regions that underlie variation in fruit traits, but also of variation related to additional plant traits and adaptive responses to climate variation. These results allowed validation of our approach because different marker-trait associations mapped on chromosomal regions where other candidate genes for the same traits were previously reported. In addition, previously uncharacterized chromosomal regions were targeted as potentially involved in the expression of variable phenotypes, thus demonstrating that our tomato collection is a precious reservoir of diversity and an excellent tool for gene discovery.

Brazil is the largest consumer and third highest producer of common beans (Phaseolus vulgaris L.) worldwide. Since the 1980s, the commercial Carioca variety has been the most consumed in Brazil, followed by Black and Special beans. The present study evaluates genetic diversity and population structure of 185 Brazilian common bean cultivars using 2827 high-quality single-nucleotide polymorphisms (SNPs). The Andean allelic introgression in the Mesoamerican accessions was investigated, and a Carioca panel was tested using an association mapping approach. The results distinguish the Mesoamerican from the Andean accessions, with a prevalence of Mesoamerican accessions (94.6%). When considering the commercial classes, low levels of genetic differentiation were seen, and the Carioca group showed the lowest genetic diversity. However, gain in gene diversity and allelic richness was seen for the modern Carioca cultivars. A set of 1060 ‘diagnostic SNPs’ that show alternative alleles between the pure Mesoamerican and Andean accessions were identified, which allowed the identification of Andean allelic introgression events and shows that there are putative introgression segments in regions enriched with resistance genes. Finally, genome-wide association studies revealed SNPs significantly associated with flowering time, pod maturation, and growth habit, showing that the Carioca Association Panel represents a powerful tool for crop improvements.